Abstract:Liquid metals are employed as a coolant in liquid metal fast reactors (LMFR) and are considered breeders and coolants for future power-producing fusion reactors. The interaction of liquid metals with other coolants, such as air and water, is one of the possible occurrences in liquid metal-cooled reactors. Although the SIMMER-III computer code is currently in the validation and verification phase, it is a prospective candidate code for correctly investigating possible accidents. This study aims to determine the stability and well-posedness of the SIMMER-III code Eulerian-Eulerian two-fluid model (TFM) under any such accident scenario. This work also considers the influence of the virtual mass force and diffusion forces in momentum on TFM stability in accelerated liquid metal, steam, and non-condensable gaseous flows throughout all flow regimes. The characteristics method is used to assess the ill-posed nature of TFM for all types of accidents and accident scenarios in a hypothetical simplified system model with several components (Lead–Lithium, non-condensable gases, and water vapor). It has been discovered that the analysis findings vary from the air and water two-component two-phase flows (characteristics roots spectrum and error growth rate patterns) and are particularly sensitive to the diffusion and virtual mass coefficients. This is because liquid metals have a higher density than liquid water and steam, resulting in strong virtual mass forces and weak diffusion forces in liquid-metal and gas two-phase flows. Because of the extensive range of possible interactions between fluxes of different components, producing an accurate representation of diffusion in multi-component mixtures is difficult. Because of this, it is strongly suggested that the virtual mass coefficient and diffusion coefficients be handled more accurately for these kinds of flows. The values of these coefficients significantly affect how accurate proposed TFM predictions are, but there hasn't been much research on how to estimate them. The study also sheds light on the model's accuracy and highlights the areas where the model's predictions will be mathematically trustworthy.

FESTIM: An open-source code for hydrogen transport simulations

Kinetic surface model in FESTIM: Verification and Validation

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

An OpenFOAM multi-region solver for tritium transport modeling in fusion systems

Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations

FERMI: Fusion Energy Reactor Models Integrator

FISH: A 3D parallel MHD code for astrophysical applications

FLIPEC, an Ideal MHD free-boundary axisymmetric equilibrium solver in the presence of macroscopic flows.

Hermes-3: Multi-component plasma simulations with BOUT++

The CXSFIT spectral fitting code: Past, present and future

Self-consistent simulation of plasma scenarios for ITER using a combination of 1.5D transport codes and free-boundary equilibrium codes

SIMMER code two-fluid model well-posedness and stability for simulating liquid metal (lead–lithium eutectic), water vapor, and non-condensable gases multi-component multi-phase flows

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

A package to bridge experimental tokamak data to modelling workflows for heating and transport

TDS Simulator: A MATLAB App to model temperature-programmed hydrogen desorption

FreeMHD: validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

Development and verification of tritium transport code based on RELAP5/MOD3.3 with generic model towards COOL blanket

FUSE (Fusion Synthesis Engine): A Next Generation Framework for Integrated Design of Fusion Pilot Plants

IDEFIX: a versatile performance-portable Godunov code for astrophysical flows

FEQIS: A free-boundary equilibrium solver for integrated modeling of tokamak plasmas

TRANSP integrated modeling code for interpretive and predictive analysis of tokamak plasmas