Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Coupled neutronic-thermal-hydraulic simulations of the European SFR core

Within the European SFR (Sodium-cooled Fast Reactor) - Safety Measures Assessment and Research Tools (ESFR-SMART) project, steady-state coupled simulation of the ESFR core has been performed using several core analysis packages, with the objective of quantifying the coupling effect. Focus is on the fuel Doppler effect and coolant expansion effect. Standalone neutronics calculations in TRACE/PARCS (PSI), DYN3D (HZDR) and WIMS (Jacobs) showed superb agreement with the reference Serpent power distribution (root mean square - rms discrepancies of 1.3%, 1.5% and 0.7% respectively). Results for COUNTHER (CIEMAT) were also in reasonable agreement, but with a somewhat higher discrepancy of 3.7%. Temperature distributions from thermal-hydraulic calculations were also compared and are found to be in good agreement. The effect of Doppler and coolant density feedback on core power distribution was predicted by TRACE/PARCS, DYN3D and WIMS to be between 0.4% and 0.8% rms difference in assembly powers. Reactivity coefficients for perturbations in the inlet temperature, flow rate and core power were shown to be negative for these three codes, with values of roughly -0.5 pcm/C. degree, -0.3 pcm/C. degree and -3.5 pcm/C. degree respectively. A preliminary investigation of differential thermal expansion effects indicates that this may have a significant effect on core power distribution of a few %, greater than anticipated a priori and may warrant inclusion in coupled core analysis to ensure the accurate calculation of power distributions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Coupled Multiphysics Earthquake Simulation of MSRs

Earthquakes affect numerous nuclear power plants around the world and safe earthquake design requires that the physics of reactors subjected to earthquakes is well understood. This project attempts to understand the physics of molten salt reactors (MSRs) subjected to earthquakes using coupled multiphysics tools in MOOSE. A beyond-design-basis scenario of an unprotected pump trip (ULOFA) during an earthquake is analyzed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2022

COUPLED MULTIPHYSICS PRIMARY LOOP SIMULATIONS OF THE Mk1-FHR IN THE VIRTUAL TEST BED

The Virtual Test Bed (VTB) repository hosts a wide range of challenge problems to showcase modeling and simulation capabilities to support advanced reactor demonstration. A coupled multiphysics model of the Mk1 pebble-bed fluoride-salt-cooled high-temperature reactor (PB-FHR) is presented here. The analysis leverages NEAMS tools (Griffin), SAM, Pronghorn, and the MOOSE heat conduction module) for core neutronics, thermal hydraulics of the core and the primary loop, and multiscale fuel performance simulations. The analysis was entirely created by coupling standalone simulations of the reactor previously available on the VTB. All input files and documentation developed for this example are available on the VTB website: \url{mooseframework.inl.gov/virtual_test_bed/}}. This model was featured in the NRIC Tech Talk presented in December 2021.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Simple Hamiltonian for quantum simulation of strongly coupled $(2+1)D$ SU(2) lattice gauge theory on a honeycomb lattice

Here, we find a simple spin Hamiltonian to describe physical states of $(2+1)$-dimensional SU(2) lattice gauge theory on a honeycomb lattice with a truncation of the electric field representation at $j_{max}=\frac{1}{2}$. The simple spin Hamiltonian contains only local products of Pauli matrices, even though Gauss’s law has been completely integrated out.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Accuracy of the explicit energy-conserving particle-in-cell method for under-resolved simulations of capacitively coupled plasma discharges

The traditional explicit electrostatic momentum-conserving particle-in-cell algorithm requires strict resolution of the electron Debye length to deliver numerical stability and accuracy. The explicit electrostatic energy-conserving particle-in-cell algorithm alleviates this constraint with minimal modification to the traditional algorithm, retaining its simplicity, ease of parallelization, and acceleration on modern supercomputing architectures. In this article, we apply the algorithm to model a one-dimensional radio frequency capacitively coupled plasma discharge relevant to industrial applications. The energy-conserving approach closely matches the results from the momentum-conserving algorithm and retains accuracy even for cell sizes up to 8 times the electron Debye length. For even larger cells, the algorithm loses accuracy due to poor resolution of steep gradients within the radio frequency sheath. Accuracy can be recovered by adopting a non-uniform grid, which resolves the sheath and allows for cell sizes up to 32 times the electron Debye length in the quasi-neutral bulk of the discharge. The effect is an up to 8 times reduction in the number of required simulation cells, an improvement that can compound in higher-dimensional simulations. We therefore consider the explicit energy-conserving algorithm as a promising approach to significantly reduce the computational cost of full-scale device simulations and a pathway to delivering kinetic simulation capabilities of use to industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An Assessment of the Error Due to Computing Waste Isolation Pilot Plant Porosity Using the Porosity Response Surface Approach

The Waste Isolation Pilot Plant Performance Assessment (WIPP PA) must predict the likelihood that radionuclides will escape into the biosphere via mechanisms that depend on geohydraulic flow. Ideally, one would predict the geohydraulic flow using coupled geohydraulic and geomechanical simulations, but such coupled simulations are not computationally tractable. Instead, Sandia has historically used a look-up table of porosities for a given fluid pressure and time, called the porosity response surface, but this approach can introduce porosity errors because it largely ignores the porosity’s dependence on the past fluid pressure history. This report discusses efforts to quantify these porosity errors for both the legacy and new porosity response surfaces. Six hundred different fluid pressure histories were fed through the legacy/new geomechanical model and the legacy/new porosity response surface to generate six hundred porosity error histories. The error associated with the legacy porosity surface was substantial, while the error associated with the new porosity surface was typically small, except when fluid pressures exceeded the lithostatic pressure at the repository. In response to the errors at high pressures, a preliminary study of the WIPP PA’s sensitivity to these porosity errors was conducted. The study found that reducing the porosity errors at high pressures negligibly affected predictions of radionuclide releases. Finally, an initial machine-learned model for porosity was developed. This ML model significantly reduced the porosity error at high pressures, but sizable errors remained, so more development is necessary before coupling an ML model to the geohydraulic model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules

Advances in machine learning provide the ability to leverage data from expensive simulations of high-energy-density experiments to significantly cut down on computational time and costs associated with the search for optimal target designs. This study presents an application of cutting-edge Bayesian optimization methods to the one-dimensional (1D) design optimization of double shell graded layer targets for inertial confinement fusion experiments. This investigation attempts to reduce hydrodynamic instabilities while retaining high yields for future NIF experiments. Machine learning methods can use predictive physics simulations to identify graded layer designs from within the vast design space that demonstrate high predicted performance, including novel designs with high uncertainty in performance that may hold unexpected promise. By applying machine learning tools to the simulation design, we map the trade-off between 1D yield and instability, specifically isolating parameter ranges, which maintain high performance while showing significantly improved Rayleigh–Taylor stability over the point design. Furthermore, the groundwork laid in this study will be a useful design tool for future NIF experiments with graded layer targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Reconstruction and interpretation of ionization asymmetry in magnetic confinement via synthetic diagnostics

Abstract Strong poloidal refueling asymmetry in the DIII-D tokamak is inferred from line radiation measurements. Synthetic diagnostics in neutral transport modeling coupled to gyrokinetic simulations illuminate implications for the plasma flow profile in the scrape-off layer of single-null beam-driven discharges. Recycling occurs primarily either on the inner or outer divertor legs, depending on the toroidal magnetic field direction. By reversing the toroidal magnetic field, the observed line radiation asymmetry is nearly eliminated or reversed. It is determined that, while relatively simple physics can describe the observed ionization asymmetry, predicting the overall brightness of the hydrogenic Lyman- α signal requires detailed simulation of the plasma and resulting turbulence. To this end, kinetic plasma simulations fully coupled to comprehensive neutral transport calculations—a novel capability—provide first-principles reproduction of Lyman- α observations on DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Coupled Experimental and Simulation Approach to Investigate the Impact of Grain Growth, Amorphization, and Grain Subdivision in Accident Tolerant U 3 Si 2 Light Water Reactor Fuel

Triuranium disilicide (U 3 Si 2 ) is being considered as an alternative accident tolerant fuel (ATF) due to its higher thermal conductivity. However, there is uncertainty in its use due to a lack of irradiation data at light water reactor (LWR) conditions. In this project, we used an integrated experimental and simulation approach to answer two significant questions pertaining to U 3 Si 2 that were not investigated by the Fuel Cycle Research and Development (FCRD) and Nuclear Energy Advanced Modeling and Simulation (NEAMS) programs: (1) Will grain growth in the hotter portions of the fuel significantly impact U 3 Si 2 LWR fuel behavior? (2) Under what conditions do grain subdivision and amorphization occur in U 3 Si 2 and will either occur at LWR conditions? Our approach to answer these questions used in situ ion irradiation and annealing to investigate grain growth, grain subdivision, and amorphization, along with mesoscale simulations using the MARMOT tool. We found that while the grain boundary mobility of U 3 Si 2 is larger than that of UO 2 , grain growth in U 3 Si 2 fuel pellets will be less than in UO 2 and may be neglected due to the lower thermal conductivity and thus lower centerline temperature. We also found that amorphization will not occur in U 3 Si 2 above 600 K, such that it is not likely to occur at typical LWR operating conditions. Grain subdivision will occur at LWR conditions and will likely occur at lower burnups than in UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗