A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number (original) (raw)
Related papers
An incompressible magnetohydrodynamics solver
A two-dimensional nonlinear incompressible magneto-hydrodynamics (MHD) code is presented to solve steady state or transient charged or neutral convection problems. The flows considered are incompressible and the divergence conditions on the velocity and magnetic fields are handled by a relaxation scheme as pseudo-iterations between the real time levels. The numerical method takes the advantage of a matrix distribution scheme that runs on structured or unstructured triangular meshes. The time-dependent algorithm utilizes a semi-implicit dual time stepping technique with multistage Runge-Kutta algorithm. It is shown that the code developed here can be accurately used for the solutions of incompressible Navier Stokes or MHD equations.
Numerical analysis of Backward-Euler discretization for simplified magnetohydrodynamic flows
Applied Mathematical Modelling, 2015
Magnetohydrodynamics (MHD) studies consider the dynamics of electrically conducting fluids. MHD are described by a set of equations, which are a combination of the Navier-Stokes equations of fluid dynamics and Maxwell's equations of electromagnetism. In most terrestrial applications, MHD flows occur at low magnetic Reynolds numbers. In this study, we apply the finite element method to time-dependent MHD flows with Backward-Euler discretization at low magnetic Reynolds number. We present a comprehensive error analysis for fully discrete approximation. Finally, the effectiveness of the method is illustrated by several numerical examples.
Computational Aspects of Magneto-fluid-dynamics Formulations
40th AIAA Plasmadynamics and Lasers Conference, 2009
There are two formulations available for modeling magnetofluiddynamics (MFD) flow fields. (1) The classical full MFD (FMFD) equations and (2) low magnetic Reynolds number approximation. In addition, a third formulation-decomposed MFD (DMFD) equations, based on decomposition of total magnetic field for solving full MFD equations has also been proposed and validated in the available literature. In this work, computational performance of the existing MFD formulations has been investigated. The governing equations are transformed to a generalized computational domain and discretized using a finite difference technique. Time-explicit multistage Runge-Kutta scheme augmented with total variation diminishing (TVD) limiters for time integration is implemented. Time-dependent flow field over a flat plate with imposed magnetic field has been considered for investigating the performance of each MFD formulation. It has been found that the performance of DMFD equations is increased compared to FMFD equations when strong magnetic field is applied. Furthermore, it is shown that the low magnetic Reynolds number approach requires minimum amount of time to provide the solution and remains valid only for small values of magnetic Reynolds numbers. Nomenclature ⎪ ⎭ ⎪ ⎬ ⎫ ⎪ ⎩ ⎪ ⎨ ⎧ = I = identity tensor J = Jacobian of transformation J r = current density vector ∞ M = free stream Mach number p = pressure q = dynamic pressure Q = field vector
A finite element method for magnetohydrodynamics
Computer Methods in Applied Mechanics and Engineering, 2001
This paper presents a ®nite element method for the solution of 3D incompressible magnetohydrodynamic (MHD)¯ows. Two important issues are thoroughly addressed. First, appropriate formulations for the magnetic governing equations and the corresponding weak variational forms are discussed. The selected B; q formulation is conservative in the sense that the local divergencefree condition of the magnetic ®eld is accounted for in the variational sense. A Galerkin-least-squares variational formulation is used allowing equal-order approximations for all unknowns. In the second issue, a solution algorithm is developed for the solution of the coupled problem which is valid for both high and low magnetic Reynolds numbers. Several numerical benchmark tests are carried out to assess the stability and accuracy of the ®nite element method and to test the behavior of the solution algorithm. Ó
Journal of Applied Mathematics, 2012
The flow of an incompressible electrically conducting viscous fluid in convergent or divergent channels under the influence of an externally applied homogeneous magnetic field is studied both analytically and numerically. Navier-Stokes equations of fluid mechanics and Maxwell’s electromagnetism equations are reduced into highly non-linear ordinary differential equation. The resulting non-linear equation has been solved analytically using a very efficient technique, namely, differential transform method (DTM). The DTM solution is compared with the results obtained by a numerical method (shooting method, coupled with fourth-order Runge-Kutta scheme). The plots have revealed the physical characteristics of flow by changing angles of the channel, Hartmann and Reynolds numbers.
Numerical Simulation of the Electrically Conducting Flows in the External Magnetic Field
2004
The paper contains the original mathematical model that describes the incompressible electrically conducting fluid flows under the influence of the electromagnetic field quasi-magneto-hydrodynamic (QMHD) equation system. Simplified inductionless approximation of QMHD system for the quasi-neutral liquid is constructed and implemented for numerical simulation of Marangoni convection in melted semiconductor suppressed by the static external magnetic field.
Theoretical and Computational Fluid Dynamics, 2015
We consider a plane channel flow of an electrically conducting fluid which is driven by a mean pressure gradient in the presence of an applied magnetic field that is streamwise periodic with zero mean. Magnetic flux expulsion and the associated bifurcation in such a configuration is explored using direct numerical simulations (DNS). The structure of the flow and magnetic fields in the Hartmann regime (where the dominant balance is through Lorentz forces) and the Poiseuille regime (where viscous effects play a significant role) are studied and detailed comparisons to the existing one-dimensional model of Kamkar and Moffatt (J. Fluid. Mech., Vol.90, pp 107-122, 1982) are drawn to evaluate the validity of the model. Comparisons show good agreement of the model with DNS in the Hartmann regime, but significant diferences arising in the Poiseuille regime when non-linear effects become important. The effects of various parameters like the magnetic Reynolds number, imposed field wavenumber etc. on the bifurcation of the flow are studied. Magnetic field line reconnections occuring during the dynamic runaway reveal a specific two-step pattern that leads to the gradual expulsion of flux in the core region.
Finite element method in applications of magnetohydrodynamics
The magnetically induced flow was examined numerically using a computational code based on the finite element method with the streamline-upwind/pressure-stabilized Petrov-Galerkin approach. The mathematical model considers an incompressible unsteady flow with a low frequency and low induction magnetic field. The validation of the magnetic force calculation was carried out on a cylindrical cavity, where the time-dependent electric potential and current density distribution can be derived analytically. The flow under the rotating magnetic field was simulated for the axisymmetric cylindrical and non-axisymmetric square cavity. The effect of the different geometries on the distribution of the time-averaged magnetic force and magnetically driven rotating flow was discussed.
Computer Methods in Applied Mechanics and Engineering, 2010
We introduce and analyze a mixed finite element method for the numerical discretization of a stationary incompressible magnetohydrodynamics problem, in two and three dimensions. The velocity field is discretized using divergence-conforming Brezzi-Douglas-Marini (BDM) elements and the magnetic field is approximated by curl-conforming Nédélec elements. The H 1-continuity of the velocity field is enforced by a DG approach. A central feature of the method is that it produces exactly divergence-free velocity approximations, and captures the strongest magnetic singularities. We prove that the energy norm error is convergent in the mesh size in general Lipschitz polyhedra under minimal regularity assumptions, and derive nearly optimal a priori error estimates for the two-dimensional case. We present a comprehensive set of numerical experiments, which indicate optimal convergence of the proposed method for twodimensional as well as three-dimensional problems.
A fully implicit numerical method for single-fluid resistive magnetohydrodynamics
Journal of Computational Physics, 2006
We present a nonlinearly implicit, conservative numerical method for integration of the single-fluid resistive MHD equations. The method uses a high-order spatial discretization that preserves the solenoidal property of the magnetic field. The fully coupled PDE system is solved implicitly in time, providing for increased interaction between physical processes as well as additional stability over explicit-time methods. A high-order adaptive time integration is employed, which in many cases enables time steps ranging from one to two orders of magnitude larger than those constrained by the explicit CFL condition. We apply the solution method to illustrative examples relevant to stiff magnetic fusion processes which challenge the efficiency of explicit methods. We provide computational evidence showing that for such problems the method is comparably accurate with explicit-time simulations, while providing a significant runtime improvement due to its increased temporal stability.