Abstract:It is well-known that magnetohydrodynamics (MHD) dominates the dynamic of the liquid metal flows inside the breeding blankets (BB) of future nuclear fusion plants by magnetic confinement. MHD is a multiphysics phenomenon involving both electromagnetism and incompressible fluid mechanics. From the computational point of view, the simulation of MHD flows in fusion relevant conditions entails a significant challenge. Indeed, due to the shape of the induced electrical currents inside the bulk of the fluid, high spatial resolutions are needed to capture the large gradients found in boundary lagers and 3D effects. Besides, solving the equations accurately typically requires very small time steps for the transient algorithms. Over the past few decades, some parallel MHD codes have been developed with success to simulate complex flows in increasingly realistic geometries. Among them, the MHD tools of commercial CFD platforms have attracted attention due to their relatively soft learning curve. Most of these codes are based on the so called -formulation which, by applying the divergence free condition of the current density to the Ohms law, reduces the electromagnetic part of the problem to a single Poisson equation. As a downside, the approach segregates the fluid and electromagnetic problem. In practice, this establishes important limits to the mesh element size, to the mesh quality and to the time-step needed to obtain accurate and stable solutions that maintains charge conservation at a discrete level. In this work, these limits are explored for the commercial platform ANSYS-Fluent using a test geometry under different conditions. As an alternative, a new code based on Finite Element Methods (FEM) is introduced as well. This open-source code, called GridapMHD (https://github.com/gridapapps/GridapMHD.jl), aims at solving the full set of MHD equations using a monolithic approach. GridapMHD is still in early stages of development but it has already shown promising results.

Development of a 2D code framework for high-energy-density plasmas of X-pinch

A 2-D Axisymmetric Magneto-Hydrodynamic Model of a DC Arc Plasma Torch and Its Solution Methodology

A high-order finite-difference solver for direct numerical simulations of magnetohydrodynamic turbulence

Hybrid fluid–particle modeling of shock-driven hydrodynamic instabilities in a plasma

Approach to Nonlinear Magnetohydrodynamic Simulations in Stellarator Geometry

A globally conservative finite element MHD code and its application to the study of compact torus formation, levitation and magnetic compression

Accuracy and scalability of incompressible inductionless MHD codes applied to fusion technologies

Coupling multi-fluid dynamics equipped with Landau closures to the particle-in-cell method

HOW-MHD: A High-Order WENO-Based Magnetohydrodynamic Code with a High-Order Constrained Transport Algorithm for Astrophysical Applications

Development of a tokamak magnetohydrodynamic code with the discontinuous Galerkin and Weighted Essentially Non-Oscillatory method

3-D Hybrid Simulation of Low-current Vacuum Plasma Jets under Axial Magnetic Fields

Three Dimensional Equilibrium Solutions for A Current-Carrying Reversed-Field Pinch Plasma with A Close-Fitting Conducting Shell

Two-dimensional Numerical Simulation on Evolution of Arc Plasma

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

Development of a gyrokinetic-MHD energetic particle simulation code Part I: MHD version

Particle-in-cell Simulation Method for Macroscopic Degenerate Plasmas.

Exploring the Microscopic Physical Processes of Z-pinch by a Modified Electrostatic Direct Implicit Particle-in-Cell Algorithm

Three-Dimensional Direct-Implicit Particle-in-Cell Model Using Trilinear Anisotropic Immersed-Finite-Element for Plasma Propulsion

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

The development of an implicit full f method for electromagnetic particle simulations of Alfvén waves and energetic particle physics

Kinetic simulations of low-pressure inductively coupled plasma: an implicit electromagnetic PIC/MCC model with the ADI-FDTD method