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.

The Advanced Computing Hub at BSC: Improving fusion codes following modern software engineering standards

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

Best-Estimate for System Codes (BeSYC): A New Software to Perform Best-Estimate Plus Uncertainty Analyses with Thermal-Hydraulic and Safety System Codes for Both Fusion and Fission Scenarios

The Plasma Simulation Code: A modern particle-in-cell code with load-balancing and GPU support

Experience and Analysis of Scalable High-Fidelity Computational Fluid Dynamics on Modular Supercomputing Architectures

Understanding Large-Scale Plasma Simulation Challenges for Fusion Energy on Supercomputers

Scalable Flow Simulations with the Lattice Boltzmann Method

Ktirio Urban Building: A Computational Framework for City Energy Simulations Enhanced by CI/CD Innovations on EuroHPC Systems

SPH-EXA: Enhancing the Scalability of SPH codes Via an Exascale-Ready SPH Mini-App

US ATLAS and US CMS HPC and Cloud Blueprint

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

A high performance scientific cloud computing environment for materials simulations

High Performance Computing in Medical Image Analysis HuSSaR

Fusion research using Azure A100 HPC instances

Exascale Simulations of Fusion and Fission Systems

Accelerating cluster dynamics simulation of fission gas behavior in nuclear fuel on deep computing unit–based heterogeneous architecture supercomputer

HPC Alongside User-space Kubernetes

UNEDF:Advanced Scientific Computing Collaboration Transforms the Low-Energy Nuclear Many-Body Problem

UNEDF: Advanced Scientific Computing Collaboration Transforms the Low-Energy Nuclear Many-Body Problem

High-Performance I/O: HDF5 for Lattice QCD

HACC: Simulating Sky Surveys on State-of-the-Art Supercomputing Architectures