Open Access
ARTICLE
The Chebyshev Tau Spectral Method for the Solution of the Linear Stability Equations for Rayleigh-Bénard Convection with Melting
Departamento de Termofluidos, Facultad de Ingeniería, Universidad Nacional Autónoma de México, Mexico D.F. C.P. 04510, ravila@servidor.unam.mx
Centro de Investigación en Energía, Universidad Nacional Autónoma de México, AP. P. 34, 62580, Temixco Mor. México, erm@mazatl.cie.unam.mx
Center for Aerospace Research & Education, University of California, Irvine
Computer Modeling in Engineering & Sciences 2009, 51(1), 73-92. https://doi.org/10.3970/cmes.2009.051.073
Abstract
A Chebyshev Tau numerical algorithm is presented to solve the perturbation equations that result from the linear stability analysis of the convective motion of a fluid layer that appears when an unconfined solid melts in the presence of gravity. The system of equations that describe the phenomenon constitute an eigenvalue problem whose accurate solution requires a robust method. We solve the equations with our method and briefly describe examples of the results. In the limit where the liquid-solid interface recedes at zero velocity the Rayleigh-Bénard solution is recovered. We show that the critical Rayleigh number Rac and the critical wave number ac are monotonically decreasing functions of the rate of melting of the solid. We conclude that the parameters Rac and ac are independent functions of the Prandtl number in the range 1 ≤ Pr ≤ 10,000. We also show that as the Pr number is reduced, Pr < 1, the critical parameters are nonmonotonic functions of the rate of melting.Keywords
Cite This Article
This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.