A visual incompressible magneto-hydrodynamics solver with radiation, mass, and heat transfer

摘要

A visual two-dimensional (2D) nonlinear magneto-hydrodynamics (MHD) code that is able to solve steady state or transient charged or neutral convection problems under the radiation, mass, and heat transfer effects is presented. The flows considered are incompressible and the divergence conditions on the velocity and magnetic fields are handled by similar relaxation schemes in the form of pseudo-iterations between the real time levels. The numerical method utilizes a matrix distribution scheme that runs on structured or unstructured triangular meshes. The time-dependent algorithm developed here utilizes a semi-implicit dual time stepping technique with multistage Runge-Kutta (RK) algorithm. It is possible for the user to choose different normalizations (natural, forced, Boussinesq, Prandtl, double-diffusive and radiation convection) automatically. The code is visual and runs interactively with the user. The graphics algorithms work multithreaded and allow the user to follow certain flow features (color graphs, vector graphs, one-dimensional profiles) during runs, see (Comput. Fluids 2007; 36:961-973) for details. With the code presented here nonlinear steady or time-dependent evolution of heated and stratified neutral and charged liquids, convection of mixture of neutral and charged gases, double-diffusive and salinity natural convection flows with internal heat generation/absorption and radiative heat transfer flows can be investigated. In addition, the numerical method (combining concentration, radiation, heat transfer, and MHD effects) takes the advantage of local time stepping and employs simplified residual jacobian matrix to increase pseudo-convergence rate. This code is currently being improved to simulate three-dimensional problems with parallel processing.