Short Communication: A numerical method for diffusive transport with moving boundaries and discontinuous material properties
✍ Scribed by Li, Zhilin; McTigue, David F.; Heine, Jan T.
- Publisher
- John Wiley and Sons
- Year
- 1997
- Tongue
- English
- Weight
- 116 KB
- Volume
- 21
- Category
- Article
- ISSN
- 0363-9061
No coin nor oath required. For personal study only.
✦ Synopsis
A method for the numerical simulation of diffusive transport with moving boundaries is developed and tested. The variable domain is mapped onto a fixed region, which introduces a term of convective form to the transformed governing equation. The resulting convection/diffusion equation is solved by a finitedifference method. An 'Immersed Interface' Method (IIM) is introduced in order to retain second-order accuracy near discontinuities in material properties, where the solution is not smooth. The method performs well in benchmark calculations against an analytical solution. The IIM scheme is capable of treating a strong discontinuity in the gradient, and it is readily extended to two or three dimensions. The methods are illustrated through a calculation for the temperature profile in a growing continental ice sheet, in which the thermal properties are discontinuous at the rock/ice interface.