Numercal Simulation Of MHD Based On Finite Volume Method

Numerical simulation is one of the most important methods to study the magnetic fluid dynamics.For nearly a decade,a considerable number of researchers have developed a series of codes such as M3 D and NIMROD,for studying magnetic fluid dynamics,which have get a lot of very meaningful simulation results.These are widely used numerical procedures,mostly based on traditional discrete operator to prepare and implement the algorithm.In the future,if the fusion reactor in order to achieve long-term steady-state operation,the nonlinear effects will be a problem for researchers.Different from the traditional operator discrete algorithm,Because of its inherent conservation properties the flux difference splitting scheme(FDS)based on the finite volume method even in the nonlinear development of magnetohydrodynamics,it is also possible to satisfy the conservation of the relevant physical quantities,which guarantees the accuracy and high fidelity of the numerical results on a long time scale.This algorithm can be used in the small volume of the yuan boundary(and the nonlinear effect of the most obvious place)to get a very good high precision Riemannian approximate solution.Because of its high accuracy and numerical stability,it has been widely used in aviation design,weather forecasting,astrophysics and other fields.In this paper,based on the finite volume method flux differentialsplitting scheme,we develop a numerical capability to solve one/two dimensional ideal/resistive MHD equations by combining a multi-state Harten-Lax-Van Leer approximate Riemann solver with the hyperbolic divergence cleaning technique,high order shock-capturing reconstruction schemes,and a third order total variance diminishing Runge-Kutta time evolving scheme.And the simulation results are compared with those of the previous studies.The correctness,high precision and good numerical stability of the algorithm are verified.It is worth mentioning that the simulation of the double-torn touch instability in the flat plate shape is a good way to restore the critical phase of the magnetic field reconnection.In addition,the effect of changing the resistivity size and the distance between the two rational surfaces on the instability of the double tear is discussed.The simulation results are in good agreement with those of the previous studies.The results provide a feasible method for the numerical simulation of magnetic fluid dynamics(MHD).In this paper,based on the finite volume method,the corresponding MHD equations are solved numerically,and the simulation results are obtained.It provides a feasible method for the study of magnetic fluid dynamics(MHD),and has a certain effect on the numerical simulation of nonlinear behavior.