Research Papers: Porous Media

Moving Boundary-Moving Mesh Analysis of Freezing Process in Water-Saturated Porous Media Using a Combined Transfinite Interpolation and PDE Mapping Methods

[+] Author and Article Information
P. Rattanadecho1

Faculty of Engineering,  Thammasat University (Rangsit Campus), Pathumthani 12121, Thailandratphadu@engr.tu.ac.th

S. Wongwises

Department of Mechanical Engineering,  King Mongkut's University of Technology Thonburi, 91 Suksawas 48, Rasburana, Bangkok 10140, Thailand


Corresponding author.

J. Heat Transfer 130(1), 012601 (Jan 25, 2008) (10 pages) doi:10.1115/1.2780177 History: Received June 14, 2006; Revised April 20, 2007; Published January 25, 2008

This paper couples the grid generation algorithm with the heat transport equations and applies them to simulate the thermal behavior of freezing process in water-saturated porous media. Focus is placed on establishing a computationally efficient approach for solving moving boundary heat transfer problem, in two-dimensional structured grids, with specific application to an undirectional solidification problem. Preliminary grids are first generated by an algebraic method, based on a transfinite interpolation method, with subsequent refinement using a partial differential equation (PDE) mapping (parabolic grid generation) method. A preliminary case study indicates successful implementation of the numerical procedure. A two-dimensional solidification model is then validated against available analytical solution and experimental results and subsequently used as a tool for efficient computational prototyping. The results of the problem are in good agreement with available analytical solution and experimental results.

Copyright © 2008 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 2

The parametric domain with f(u,w) specified on planes of constant u,w

Grahic Jump Location
Figure 3

Strategy for calculation

Grahic Jump Location
Figure 4

Validation test for a planar freezing font in a phase-change slab

Grahic Jump Location
Figure 5

The interface deformation in computational domain with different numerical grids: (a) 100×50 grids (b) 100×100 grids

Grahic Jump Location
Figure 6

Grid simulating the deformation of an interface: (a) freezing time of 30s, (b) freezing time of 60s, (c) freezing time of 90s, (d) freezing time of 120s, (e) freezing time of 150s, and (f) freezing time of 180s

Grahic Jump Location
Figure 7

Comparison of experimental data and simulated freezing front from present numerical study

Grahic Jump Location
Figure 8

The simulations of temperature distribution (Unit: °C): (a) freezing time of 30s, (b) freezing time of 60s, (c) freezing time of 90s, (d) freezing time of 120s, (e) freezing time of 150s, and (f) freezing time of 180s




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In