Computational magnetohydrodynamics

Last updated

Computational magnetohydrodynamics (CMHD) is a rapidly developing branch of magnetohydrodynamics that uses numerical methods and algorithms to solve and analyze problems that involve electrically conducting fluids. Most of the methods used in CMHD are borrowed from the well established techniques employed in Computational fluid dynamics. The complexity mainly arises due to the presence of a magnetic field and its coupling with the fluid. One of the important issues is to numerically maintain the (conservation of magnetic flux) condition, from Maxwell's equations, to avoid the presence of unrealistic effects, namely magnetic monopoles, in the solutions.

Contents

Open-source MHD software

Closed-source MHD software

See also

Related Research Articles

<span class="mw-page-title-main">Magnetohydrodynamics</span> Model of electrically conducting fluids

Magnetohydrodynamics is a model of electrically conducting fluids that treats all interpenetrating particle species together as a single continuous medium. It is primarily concerned with the low-frequency, large-scale, magnetic behavior in plasmas and liquid metals and has applications in numerous fields including geophysics, astrophysics, and engineering.

<span class="mw-page-title-main">Computational physics</span> Numerical simulations of physical problems via computers

Computational physics is the study and implementation of numerical analysis to solve problems in physics. Historically, computational physics was the first application of modern computers in science, and is now a subset of computational science. It is sometimes regarded as a subdiscipline of theoretical physics, but others consider it an intermediate branch between theoretical and experimental physics — an area of study which supplements both theory and experiment.

<span class="mw-page-title-main">Magnetohydrodynamic drive</span> Vehicle propulsion using electromagnetic fields

A magnetohydrodynamic drive or MHD accelerator is a method for propelling vehicles using only electric and magnetic fields with no moving parts, accelerating an electrically conductive propellant with magnetohydrodynamics. The fluid is directed to the rear and as a reaction, the vehicle accelerates forward.

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

In plasma physics, the particle-in-cell (PIC) method refers to a technique used to solve a certain class of partial differential equations. In this method, individual particles in a Lagrangian frame are tracked in continuous phase space, whereas moments of the distribution such as densities and currents are computed simultaneously on Eulerian (stationary) mesh points.

A magnetohydrodynamic generator is a magnetohydrodynamic converter that transforms thermal energy and kinetic energy directly into electricity. An MHD generator, like a conventional generator, relies on moving a conductor through a magnetic field to generate electric current. The MHD generator uses hot conductive ionized gas as the moving conductor. The mechanical dynamo, in contrast, uses the motion of mechanical devices to accomplish this.

Numerical resistivity is a problem in computer simulations of ideal magnetohydrodynamics (MHD). It is a form of numerical diffusion. In near-ideal MHD systems, the magnetic field can diffuse only very slowly through the plasma or fluid of the system; it is rate-limited by the inverse of the resistivity of the fluid. In Eulerian simulations where the field is arbitrarily aligned compared to the simulation grid, the numerical diffusion rate takes the form similar to an additional resistivity, causing non-physical and sometimes bursty magnetic reconnection in the simulation. Numerical resistivity is a function of resolution, alignment of the magnetic field with the grid, and numerical method. In general, numerical resistivity will not behave isotropically, and there can be different effective numerical resistivities in different parts of the computational domain. For current (2005) simulations of the solar corona and inner heliosphere, this numerical effect can be several orders of magnitude larger than the physical resistivity of the plasma.

The Advection Upstream Splitting Method (AUSM) is a numerical method used to solve the advection equation in computational fluid dynamics. It is particularly useful for simulating compressible flows with shocks and discontinuities.

<span class="mw-page-title-main">Bram van Leer</span> Dutch mathematician

Bram van Leer is Arthur B. Modine Emeritus Professor of aerospace engineering at the University of Michigan, in Ann Arbor. He specializes in Computational fluid dynamics (CFD), fluid dynamics, and numerical analysis. His most influential work lies in CFD, a field he helped modernize from 1970 onwards. An appraisal of his early work has been given by C. Hirsch (1979)

<span class="mw-page-title-main">Plasma modeling</span>

Plasma modeling refers to solving equations of motion that describe the state of a plasma. It is generally coupled with Maxwell's equations for electromagnetic fields or Poisson's equation for electrostatic fields. There are several main types of plasma models: single particle, kinetic, fluid, hybrid kinetic/fluid, gyrokinetic and as system of many particles.

<span class="mw-page-title-main">Gauss's law for magnetism</span> Foundational law of classical magnetism

In physics, Gauss's law for magnetism is one of the four Maxwell's equations that underlie classical electrodynamics. It states that the magnetic field B has divergence equal to zero, in other words, that it is a solenoidal vector field. It is equivalent to the statement that magnetic monopoles do not exist. Rather than "magnetic charges", the basic entity for magnetism is the magnetic dipole.

<span class="mw-page-title-main">Coolfluid</span> Scientific computing environment

COOLFluiD is a component based scientific computing environment that handles high-performance computing problems with focus on complex computational fluid dynamics (CFD) involving multiphysics phenomena.

Coronal seismology is a technique of studying the plasma of the Sun's corona with the use of magnetohydrodynamic (MHD) waves and oscillations. Magnetohydrodynamics studies the dynamics of electrically conducting fluids - in this case the fluid is the coronal plasma. Observed properties of the waves (e.g. period, wavelength, amplitude, temporal and spatial signatures, characteristic scenarios of the wave evolution, combined with a theoretical modelling of the wave phenomena, may reflect physical parameters of the corona which are not accessible in situ, such as the coronal magnetic field strength and Alfvén velocity and coronal dissipative coefficients. Originally, the method of MHD coronal seismology was suggested by Y. Uchida in 1970 for propagating waves, and B. Roberts et al. in 1984 for standing waves, but was not practically applied until the late 90s due to a lack of necessary observational resolution. Philosophically, coronal seismology is similar to the Earth's seismology, helioseismology, and MHD spectroscopy of laboratory plasma devices. In all these approaches, waves of various kind are used to probe a medium.

Magnetohydrodynamic turbulence concerns the chaotic regimes of magnetofluid flow at high Reynolds number. Magnetohydrodynamics (MHD) deals with what is a quasi-neutral fluid with very high conductivity. The fluid approximation implies that the focus is on macro length-and-time scales which are much larger than the collision length and collision time respectively.

<span class="mw-page-title-main">Computational astrophysics</span> Methods and computing tools developed and used in astrophysics research

Computational astrophysics refers to the methods and computing tools developed and used in astrophysics research. Like computational chemistry or computational physics, it is both a specific branch of theoretical astrophysics and an interdisciplinary field relying on computer science, mathematics, and wider physics. Computational astrophysics is most often studied through an applied mathematics or astrophysics programme at PhD level.

The Pencil Code is a high-order finite-difference code for solving partial differential equations, written in Fortran 95. The code is designed for efficient computation with massive parallelization. Due to its modular structure, it can be used for a large variety of physical setups like hydro- and magnetohydrodynamics relevant for, e.g., astrophysics, geophysics, cosmology, turbulence, and combustion. Many such setups are available as ready-to-run samples. Pencil Code is free software released under the GNU GPL v2.


Astrophysical fluid dynamics is a branch of modern astronomy which deals with the motion of fluids in outer space using fluid mechanics, such as those that make up the Sun and other stars. The subject covers the fundamentals of fluid mechanics using various equations, such as continuity equations, the Navier–Stokes equations, and Euler's equations of collisional fluids. Some of the applications of astrophysical fluid dynamics include dynamics of stellar systems, accretion disks, astrophysical jets, Newtonian fluids, and the fluid dynamics of galaxies.

<span class="mw-page-title-main">MFEM</span> Open-source C++ library

MFEM is an open-source C++ library for solving partial differential equations using the finite element method, developed and maintained by researchers at the Lawrence Livermore National Laboratory and the MFEM open-source community on GitHub. MFEM is free software released under a BSD license.

James McLellan Stone is an American astrophysicist who specialises in the study of fluid dynamics. He is currently a faculty member at the school of natural sciences at the Institute for Advanced Study. Stone is also the Lyman Spitzer Jr. Professor of Theoretical Astrophysics, emeritus, and professor of astrophysical sciences and applied and computational mathematics, emeritus, at Princeton University.

William Henry Matthaeus is an American astrophysicist and plasma physicist. He is known for his research on turbulence in magnetohydrodynamics (MHD) and astrophysical plasmas, for which he was awarded the 2019 James Clerk Maxwell Prize for Plasma Physics.

References

  1. Teyssier, R (2002). "Cosmological hydrodynamics with adaptive mesh refinement. A new high resolution program called RAMSES". Astronomy and Astrophysics. 385: 337–364. arXiv: astro-ph/0111367 . Bibcode:2002A&A...385..337T. doi:10.1051/0004-6361:20011817. S2CID   5504247.
  2. Gheller, C; Wang, P; Vazza, F; Teyssier, R (28 September 2015). "Numerical cosmology on the GPU with Enzo and Ramses". Journal of Physics: Conference Series. 640 (1): 012058. arXiv: 1412.0934 . Bibcode:2015JPhCS.640a2058G. doi:10.1088/1742-6596/640/1/012058. S2CID   118194615 . Retrieved 1 July 2016.
  3. Stone, James M.; Gardiner, Thomas A.; Teuben, Peter; Hawley, John F.; Simon, Jacob B. (September 2008). "Athena: A New Code for Astrophysical MHD". The Astrophysical Journal Supplement Series. 178 (1): 137–177. arXiv: 0804.0402 . Bibcode:2008ApJS..178..137S. doi:10.1086/588755. S2CID   10934839.
  4. Vencels, Juris; Råback, Peter; Geža, Vadims (2019-01-01). "EOF-Library: Open-source Elmer FEM and OpenFOAM coupler for electromagnetics and fluid dynamics". SoftwareX . 9: 68–72. Bibcode:2019SoftX...9...68V. doi: 10.1016/j.softx.2019.01.007 . ISSN   2352-7110.

Further reading