Premium
Numerical simulations of thermal convection under the influence of an inclined magnetic field by using solenoidal bases
Author(s) -
Yarımpabuç D.,
Tarman H.I.,
Yıldırım C.
Publication year - 2014
Publication title -
mathematical methods in the applied sciences
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.719
H-Index - 65
eISSN - 1099-1476
pISSN - 0170-4214
DOI - 10.1002/mma.3034
Subject(s) - solenoidal vector field , prandtl number , vector field , mathematics , magnetic field , galerkin method , mathematical analysis , boundary value problem , mechanics , nonlinear system , convection , classical mechanics , geometry , physics , quantum mechanics
The effect of an inclined homogeneous magnetic field on thermal convection between rigid plates heated from below under the influence of gravity is numerically simulated in a computational domain with periodic horizontal extent. The numerical technique is based on solenoidal (divergence‐free) basis functions satisfying the boundary conditions for both the velocity and the induced magnetic field. Thus, the divergence‐free conditions for both velocity and magnetic field are satisfied exactly. The expansion bases for the thermal field are also constructed to satisfy the boundary conditions. The governing partial differential equations are reduced to a system of ordinary differential equations under Galerkin projection and subsequently integrated in time numerically. The projection is performed by using a dual solenoidal bases set such that the pressure term is eliminated in the process. The quasi‐steady relationship between the velocity and the induced magnetic field corresponding to the liquid metals or melts is used to generate the solenoidal bases for the magnetic field from those for the velocity field. The technique is validated in the linear case for both oblique and vertical case by reproducing the marginal stability curves for varying Chandrasekhar number. Some numerical simulations are performed for either case in the nonlinear regime for Prandtl numbers Pr = 0.05 and Pr = 0.1. Copyright © 2014 John Wiley & Sons, Ltd.