0
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

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

1

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

## Abstract

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.

<>

## Figures

Figure 4

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

Figure 5

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

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

Figure 7

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

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

Figure 1

Physical model

Figure 2

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

Figure 3

Strategy for calculation

## Errata

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 Proceedings Articles
Related eBook Content
Topic Collections