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.

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Dynamic Optimization Of Tokamak Plasmas Via Control Parameterization And The Time-Scaling Transformation

Gyrokinetic Particle Simulation of Microturbulence for General Magnetic Geometry and Experimental Profiles

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

M3D-K Simulations of Sawteeth and Energetic Particle Transport in Tokamak Plasmas

Analysis of Liquid Metal MHD Fluid Flow and Heat Transfer Using the KAT Code

Gtc Simulation of Ideal Ballooning Mode in Tokamak Plasmas

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

A new flux coordinates-based solver for fixed-boundary tokamak equilibrium with toroidal flow

An Adaptive Moving Mesh Finite Difference Scheme for Tokamak Magneto-Hydrodynamic Simulations

A mimetic finite difference based quasi-static magnetohydrodynamic solver for force-free plasmas in tokamak disruptions

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

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

GSPD: An algorithm for time-dependent tokamak equilibria design

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

Flexible, integrated modeling of tokamak stability, transport, equilibrium, and pedestal physics

Targeting realistic geometry in Tokamak code Gysela

Numerical Simulation on the Evolution of Plasma Parameters in Tokamak Fusion Reactor

Numerical Analysis of Electromagneto-Mechanical Coupling Dynamic Response of A Typical Tokomak Structure

ATEQ: Adaptive Toroidal Equilibrium code