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.

Newtonian CAFE: a new ideal MHD code to study the solar atmosphere

CAFE: A New Relativistic MHD Code

CAFE-AMR: A computational MHD Solar Physics simulation tool that uses AMR

Mancha3D Code: Multipurpose Advanced Nonideal MHD Code for High-Resolution Simulations in Astrophysics

A New Implementation of the Magnetohydrodynamics-Relaxation Method for Nonlinear Force-Free Field Extrapolation in the Solar Corona

Mancha3D code: Multi-purpose Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics

Three-dimensional numerical simulation of magnetohydrodynamic-gravity waves and vortices in the solar atmosphere

A new MHD code with adaptive mesh refinement and parallelization for astrophysics

Implicit Solar Coronal Magnetohydrodynamic (MHD) Modeling with a Low-dissipation Hybridized AUSM-HLL Riemann Solver

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

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

One-dimensional, geometrically stratified semi-empirical models of the quiet-Sun photosphere and lower chromosphere

Simulating the Photospheric to Coronal Plasma Using Magnetohydrodynamic Characteristics. I. Data-driven Boundary Conditions

Investigating the effects of numerical algorithms on global magnetohydrodynamics (MHD) simulations of solar wind in the inner heliosphere

A Unified Radiative Magnetohydrodynamics Code for Lightning-Like Discharge Simulations

On the extrapolation of magneto-hydro-static equilibria on the sun

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

MPI-AMRVAC 2.0 for Solar and Astrophysical Applications

Global MHD Simulations of the Time-Dependent Corona

M3DIS - A grid of 3D radiation-hydrodynamics stellar atmosphere models for stellar surveys. I. Procedure, validation, and the Sun

Modified Path-conservative HLLEM Scheme for Magnetohydrodynamic Solar Wind Simulations