Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An open-source hybrid unstructured mesh - CAD fusion multiphysics analysis workflow in SALAMANDER

Plasma facing components in fusion devices will endure extreme neutron and heat fluxes. To facilitate their design using simulation tools, the open-source Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX) framework is being developed to model these components with a high-fidelity multi-physics multi-dimensional approach. It can iteratively resolve couplings between all the physics at play, from neutron radiation, to thermomechanics, to near-wall plasma dynamics. This framework is based on the Multiphysics Object Oriented Simulation Environment (MOOSE), which is developed by a collaboration of US National Laboratories since 2008, for advanced nuclear, geomechanics simulations and other applications. FENIX couples numerous simulation tools, including OpenMC, the Tritium Migration Analysis Program v8, the NekRS CFD software, and most MOOSE modules. For the coupling of radiation transport and other physics, FENIX supports a hybrid workflow between Computer Assisted Design (CAD) and unstructured mesh geometries. The CAD can be generated from skinning the unstructured mesh, to enable a coarse geometry for efficient particle transport, but still resolving the local material compositions and temperature gradients. Neutron transport is performed using DAGMC on the CAD, and Cardinal, integrated in FENIX, maps tallied quantities, such as the heat deposition or the tritium generation rates, from a tally volumetric mesh to the other physics’ unstructured mesh. This coupling was exercised on a simplified tokamak geometry, coupling neutron transport with the heat conduction equation, and on a monoblock divertor problem, coupling additionally with tritium migration. Mesh convergence studies highlight the importance of the mapping conservativeness. Coupling with thermo-mechanics is further enabled by the generalization of the approach to moving meshes. The presentation will include these coupled analysis as well as an update on status of the FENIX framework.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Rotated-Droop Control for Enhanced Stability and Power Decoupling in Microgrids With Complex Line Impedances

Classical droop control in microgrids predominantly assumes inductive line impedance, which simplifies implementation but causes power coupling and steady-state errors in systems with resistive and inductive lines. This paper proposes a rotated-droop control strategy for grid-forming inverters that reformulates the power equations by incorporating the magnitude and angle of the line impedance within a rotated reference frame. This method enhances the decoupling of active and reactive power dynamics without increasing complexity or requiring communication links. A small-signal state-space model was created to capture the dynamic behavior of the system under varying impedance parameters, preserving the original droop gains by rotating the power control structure. This allows eigenvalue-based stability analysis and enhances damping and transient performance. Simulation and experimental results validated the improved power-sharing performance, faster response, and robustness of the proposed method under different impedance conditions. This approach maintains the decentralized structure of the conventional droop control while enabling greater adaptability and scalability, making it suitable for modern inverter-based microgrids with dynamic topologies.

Campo-Ossa, Daniel Dario [Univ. of Puerto Rico, Ag↗
DOE OSTI2023

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Core Design of the Holos-Quad Microreactor

The Holos-Quad micro-reactor concept, developed by HolosGen LLC, is equipped with a 22 MWt (Mega-Watt thermal) core and an integral power conversion system converting the core thermal energy into approximately 10 MWe (Mega-Watt electric). This design can be configured to support a wide range of applications. It is a very innovative high-temperature gas-cooled reactor concept using TRI-structural ISOtropic particle fuel (TRISO) distributed in graphite hexagonal blocks, cooled with helium in a direct Brayton cycle independently executed by four Subcritical Power Modules (SPMs) fitted into a hardened 40-foot container whose dimensions are in compliance with ISO shipping containers requirements. In FY2019 HolosGen LLC was awarded by the Department of Energy Advanced Research Project Agency-Energy (DOE ARPA-E) under the MEITNER funding program. As part of the MEITNER award, the Argonne National Laboratory (ANL) contributed expertise through two specialized teams: The “Design Team” and the “Resource Team”. The Design Team was dedicated to validate feasibility of the Holos-Quad core and to optimize its core design through neutronics analyses. The Resource Team was dedicated to feasibility verification via high-fidelity codes of Holos-Quad thermal-hydraulic, heat transfer, shielding, and structural aspects. This report summarizes the activities conducted by ANL Design Team. A rigorous design approach based on multi-criteria optimization and code-to-code comparison involving stochastic and high-fidelity deterministic solutions was developed and employed at several evolutionary stages of the Holos-Quad design. Several generations of the Holos-Quad core were designed within this project before converging to the current full-scale Gen 2+ design that is detailed in this report. Figure EA-1 illustrates a cross-sectional view of Gen 2+ Holos-Quad core configuration, and Figure EA-2 provides a simplified perspective view of 1-of-4 SPMs. The Holos-Quad uses four thermal-hydraulically independent SPMs locked into stationary positions during power operation, surrounded by BeO reflector and structural component fully comprised within the dimensional constraints represented by traditional ISO containers. One of the benefits of this approach is to enable transportation of each SPM promptly after irradiation in shielded containers. The core is designed to operate for approximately 8 full-power years while the reactivity controls and power conversion system enable load-following operations. The reactivity controls are represented by independent, diversified, and redundant reactivity control systems based on control drums and redundant sets of shutdown rods. The high-fidelity simulation tools were used to assess detailed power and flux distributions of the three-dimensional full-core or quarter-core of the Gen 2+ configuration. Single-physics and multi-physics simulations of the neutronics code PROTEUS and the thermal-hydraulic code System Analysis Module (SAM) were performed to analyze the Holos design configurations with detailed high-fidelity solutions. The design work performed confirmed feasibility of the Holos-Quad concept, provided realistic design description for detailed design of the operational system, and identified several core design improvements to be further considered for future reactor development activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The MUSIC Critical Benchmark and Nuclear Data

The Measurement of Uranium Subcritical and Critical (MUSIC) experiment was a series of measurements of critical and subcritical configurations of bare highly enriched uranium. The goal was to compare measurement methods, analysis techniques, and simulation methods across regimes of criticality and to provide high-quality validation of 235 U nuclear data. A benchmark evaluation of the two critical configurations of the MUSIC experiments will soon be published in the release of the International Criticality Safety Benchmark Evaluation Project Handbook. The recent execution of the experiment aids in proper quantification of model simplifications and all uncertainties associated with the experiment. Historical benchmark evaluations are heavily relied on for uranium nuclear data validation despite the fact that the same level of documentation and comparable uncertainty analysis may not be present. The MUSIC evaluation is less likely to include “unknown unknowns” that could impede accurately modeling the system. Presented are both highly detailed and very simplified models, which represent the experimental configurations accurately, aiding the users of the benchmark for nuclear data or transport code validation. The sensitivities of k eff to nuclear data and nuclear data–related uncertainties are very similar between this experiment and previous bare uranium sphere experiments. In addition, the nuclear data uncertainties to any nuclides other than 235 U are small. For all these reasons, the recently evaluated MUSIC benchmark critical configurations could prove very useful for 235 U nuclear data validation. Currently, major libraries have good agreement with the experimental results, within 200 pcm for all nuclear data libraries, and within one standard deviation of the experimental result for most. Suggested nuclear data adjustments based on MUSIC and Lady Godiva are also presented, with posterior improvements to both the agreement in k eff and the uncertainty associated with the nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

High-pressure behavior of 3.65 Å phase: Insights from Raman spectroscopy

Abstract The 3.65 Å phase [MgSi(OH)6] is a hydrous phase that is predicted to be stable in a simplified MgO-SiO2-H2O (MSH) ternary system at pressures exceeding 9 GPa. Along cold subduction zones, it is likely to transport water, bound in its crystalline lattice, into the Earth’s interior. The 3.65 Å phase consists of Mg and Si octahedral sites attached to the hydroxyl group that forms a hydrogen bond and is predicted to undergo pressure-induced symmetrization of the hydrogen bond. Therefore, in this study, we investigate the high-pressure behavior of the 3.65 Å phase using Raman spectroscopy. We have conducted five distinct compressions up to ~60 GPa using two different pressure-transmitting media—alcohol mixture and neon. At ambient conditions, we identified vibrational modes using complementary first-principles simulations based on density functional perturbation theory. Upon compression, we note that the first derivative of the vibrational modes in the lattice region stiffens, i.e., b1lattice > 0. In contrast, the hydroxyl region softens, i.e., b1OH > 0. This is indicative of the strengthening of hydrogen bonding upon compression. We noticed a significant broadening of vibrational modes related to hydroxyl groups that are indicative of proton disorder. However, within the maximum pressures explored in this study, we did not find evidence for pressure-induced symmetrization of the hydrogen bonds. We used the pressure derivative of the vibrational modes to determine the ratio of the bulk moduli and their pressure derivative. We note that the smaller bulk moduli of hydrous phases compared to the major mantle phases are compensated by significantly larger pressure derivatives of the bulk moduli for the hydrous phases. This leads to a significant reduction in the elasticity contrast between hydrous and major mantle phases. Consequently, the detection of the degree of mantle hydration is likely to be challenging at greater depths.

Geochemistry & Geophysics↗
DOE OSTI2022

Assessment of the sea surface temperature diurnal cycle in CNRM-CM6-1 based on its 1D coupled configuration

A single-column version of the CNRM-CM6-1 global climate model has been developed to ease development and validation of the boundary layer physics and air–sea coupling in a simplified environment. This framework is then used to assess the ability of the coupled model to represent the sea surface temperature (SST) diurnal cycle. To this aim, the atmospheric–ocean single-column model (AOSCM), called CNRM-CM6-1D, is implemented in a case study derived from the CINDY2011/DYNAMO campaign over the Indian Ocean, where large diurnal SST variabilities have been well documented. Comparing the AOSCM and its uncoupled components (atmospheric SCM and oceanic SCM, called OSCM) highlights the fact that the impact of coupling in the atmosphere results from both the possibility to take into account the diurnal variability of SST, which is not usually available in forcing products, and the change in mean state SST as simulated by the OSCM, with the ocean mean state not being heavily impacted by the coupling. This suggests that coupling feedbacks in the 3D model do not arise from the coupling of ocean and atmosphere vertical column physics but are more due to the large-scale dynamics resolved by the 3D model. Additionally, a sub-daily coupling frequency is needed to represent the SST diurnal variability, but the choice of the coupling time step between 15 min and 3 h does not impact the diurnal temperature range simulated much. The main drawback of a 3 h coupling is delaying the SST diurnal cycle by 5 h in asynchronous coupled models. Overall, the diurnal SST variability is reasonably well represented in CNRM-CM6-1 with a 1 h coupling time step and the upper-ocean model resolution of 1 m. This framework is shown to be a very valuable tool to develop and validate the boundary layer physics and the coupling interface. It highlights the interest to develop other atmosphere–ocean coupling case studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Flexible cryogenic air separation unit—An application for low-carbon fossil-fuel plants

In this study, the rapid integration of intermittent renewable sources into the electricity grid is driving the need for a flexible cryogenic air separation unit (ASU) coupled with a low-carbon fossil-fuel plant. However, the state-of-the-art ASU process is highly integrated and nonlinear, which can significantly restrict its ramping rate. In this work, we study the fundamental dynamics of a state-of-the-art double-column ASU, focusing on the dynamic characteristics of the highly integrated and nonlinear heat and mass transfers, and we propose basic control methods for achieving high process ramping rates. We found that the vapor–liquid countercurrent flows in the low-pressure column are critical to the cryogenic rectification of air, which governs the ramp rate of the ASU. This countercurrent flow structure is created through a complex heat integration process in the ASU. Here, this process is simplified as a countercurrent heat transfer to reduce the complexity of studying the ASU dynamics. To preserve this flow structure, the heat integration is maintained so that its heat duty follows the ASU load, using several basic controllers on the critical stream flowrates. A flow-driven dynamic ASU model is built in Aspen Plus Dynamics to capture the basic ramping dynamics. This model revealed a fundamental mismatch in the dynamics in the ASU column that causes a significant loss of O 2 product purity when ramping down the ASU and slightly increases the purity when ramping up. Based on these findings, we propose a basic control method for rapidly ramping down/up the ASU while maintaining O 2 product purity. Simulation results show the ASU basic dynamic process successfully ramps at a rate up to 10 %/min (40–100 % load) while maintaining the O 2 product purity at 95.2–95.6 mol%.

42 ENGINEERING↗
DOE OSTI2021

Modelling Nuclear Thermal Propulsion Reactor Startup Transients

The National Aeronautics and Space Administration (NASA) has set the goal of a manned mission to Mars by the year 2030 [1] and charged the national academy of sci- ences "to identify primary technical and programmatic chal- lenges, merits, and risks for maturing space nuclear propulsion technologies of interest to a future human Mars exploration mission" [2]. One relevant technology, nuclear thermal propul- sion (NTP), has notable advantages over traditional chemical rockets; most important among them is the ability to produce larger specific impulse on the order of 900s. The reduction of mission time is crucial for a manned mission to Mars to reduce the risk for the crew. Due to its higher specific impulse, NTP satisfies this need and is pursued as one technology to get humans to Mars [3, 4]. The construction of an NTP sys- tem has to negotiate several challenges laid out in Ref. [2]; one of these challenges is the need to startup the NTP sys- tem from essentially cold conditions to full power within one minute. This paper focuses on studying the neutronics and thermal-hydraulics behavior of a simplified NTP model dur- ing prescribed reactivity insertions and mass flow rate (MFR) ramps. It is the goal of this paper to investigate startup, peak material temperatures, and average specific impulse for a low enriched Uranium (LEU), ceramic and metal material (CER- MET) NTP system when varying reactivity insertion and MFR ramps.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Machine Learning-Driven Optimization of Building Enclosures for Moisture Durability and Thermal Performance

The design of moisture-durable building enclosures with low embodied carbon often involves an iterative process of selecting the materials for the specific exposure conditions to meet the performance requirements. While hygrothermal simulations are commonly used to evaluate moisture durability, they often require advanced expertise for proper implementation. Machine learning (ML) provides a promising alternative by streamlining the design process and minimizing the reliance on complex simulations. This study presents a machine learning-based approach for predicting moisture durability in residential wall assemblies. The ML model was trained to estimate the mold index and maximum moisture content of various layers under typical exposure conditions. The model achieved a high predictive accuracy, with a coefficient of determination (R²) exceeding 0.90 when compared to traditional hygrothermal simulations on materials that were not part of training the ML model. Building on these results, the ML model was developed into a practical tool for optimizing wall assembly designs. This tool allows users to automatically optimize material selections based on energy, moisture, and carbon performance criteria. By incorporating multi-objective optimization, the tool identifies configurations that minimize embodied carbon while maintaining moisture safety and code-compliant thermal performance. Additionally, it provides insights into how material choices influence assembly durability, energy efficiency, and carbon reduction. The tool will be implemented in the Building Science Advisor (BSA) to enhance its performance and provide more granularity on the results. This research highlights the potential for ML-driven tools to simplify the design of high-performance building enclosures, offering architects and engineers a faster, more efficient way to balance critical performance factors.

Salonvaara, Mikael [ORNL] (ORCID:0000000318991554)↗
DOE OSTI2021

Assessing the Energy Resilience of Office Buildings: Development and Testing of a Simplified Metric for Real Estate Stakeholders

Increasing concern over higher frequency extreme weather events is driving a push towards a more resilient built environment. In recent years there has been growing interest in understanding how to evaluate, measure, and improve building energy resilience, i.e., the ability of a building to provide energy-related services in the event of a local or regional power outage. In addition to human health and safety, many stakeholders are keenly interested in the ability of a building to allow continuity of operations and minimize business disruption. Office buildings are subject to significant economic losses when building operations are disrupted due to a power outage. We propose “occupant hours lost” (OHL) as a means to measure the business productivity lost as the result of a power outage in office buildings. OHL is determined based on indoor conditions in each space for each hour during a power outage, and then aggregated spatially and temporally to determine the whole building OHL. We used quasi-Monte Carlo parametric energy simulations to demonstrate how the OHL metric varies due to different building characteristics across different climate zones and seasons. The simulation dataset was then used to develop simple regression models for assessing the impact of ten key building characteristics on OHL. The most impactful were window-to-wall ratio and window characteristics. The regression models show promise as a simple means to assess and screen for resilience using basic building characteristics, especially for non-critical facilities where it may not be viable to conduct detailed engineering analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Modeling the Air Quality Impacts of Future Energy Scenarios

Many greenhouse gas (GHG) emission reduction measures achieve simultaneous reductions in air pollutants. Human-Earth system models can estimate such emission changes in the energy system but using them in chemistry-transport models (CTMs) to study their air quality impacts involves resource-intensive emissions processing. This is greatly simplified by an emissions scaling approach linking state-level emissions estimated by a human-Earth system model to a CTM. A scenario continuing pre-2022 energy policy in the U.S. to 2050 shows widespread air quality improvements over the 2015 baseline from SO 2 and NO x emission reductions of 50–80% from electricity generation and light-duty vehicles. Scenarios of GHG mitigation and vehicle electrification at the state and national level add further benefits. However, PM 2.5 increases from increased use of wood heating and bioenergy suggest that additional PM 2.5 management may be needed when using biofuels. In conclusion, this approach helps assess multiple future energy scenarios efficiently without sacrificing chemical detail in the air quality simulations.

air quality↗
DOE OSTI2024

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics↗
DOE OSTI2022

A preconceptual design of an inverted Stable-Salt Reactor

This paper presents a preconceptual design of an inverted stable-salt reactor (SSR), with the unique feature that the static fuel salt is accommodated in one single bulk fuel container, while the coolant salt is circulating through the core via flow tubes penetrating the fuel container.The inverted SSR is expected to substantially decrease the capital cost compared to the state-of-the-art SSR such as the SSR-W designed by Moltex Energy, because one single gas vent at the top of the fuel container is required in the inverted SSR to prevent the leakage of the radioactive fission product gases instead of one gas vent in each of the hundreds of thousands of fuel tubes in SSR-W. In addition, the operation of the inverted SSR is largely simplified. However, numerous challenges associated with the design of the inverted SSR still remain to be overcome. Preliminary neutronics and thermal-hydraulics performance characteristics of the inverted SSR are assessed in this study. A viable core configuration of the inverted SSR is identified, and the superior characteristics of the inverted SSR, including a better neutron economy and a lower maximum fuel temperature, are demonstrated in this paper. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Reexamining supercritical gas adsorption theories in nano-porous shales under geological conditions

Adsorption is an important phenomenon in surface chemistry, especially for nano-porous shales. In the shale gas-in-place, adsorbed gas could contribute up to 85%. However, adsorption is hard to quantitatively characterize due to different adsorption mechanisms, patterns, surfaces, and pore sizes. Moreover, key thermodynamic parameters, such as the enthalpy of adsorption, are challenging to determine due to uncertainties in adsorbed gas densities used in constructing absolute isotherms. In this paper, we revisit the Brunauer, Emmett, and Teller (BET) model, analytically simplify Ono-Kondo (OK) models for subsurface shales, and compare commonly used mono- and multilayer adsorption models (e.g., Langmuir, supercritical Dubinin-Radushkevich (SDR), supercritical BET (SBET), and simplified OK (OKs) models) with recently proposed pressure-dependent adsorption densities to develop a practical and reliable methodology that can be used in the supercritical state, typical for subsurface black shale conditions.Three independent data sets were used for nitrogen and methane adsorption isotherms at different temperatures. Furthermore, we demonstrate that adsorption predicted by the SDR model is comparable or lower than that predicted by the SBET model, but higher than the amounts predicted by Langmuir and OK models. The nitrogen BET method tends to underestimate the accessible SSA for methane. Measurement of isosteric heat of adsorption is suggested instead of the experimental fitting method, due to the significant difference between the calculated results by two commonly used methods. Experimental fitting and simulation methods are also briefly reviewed to guide future research on shale gas adsorption.

58 GEOSCIENCES↗
DOE OSTI2023

Exploring the parameter space of MagLIF implosions using similarity scaling. I. Theoretical framework

Magneto-inertial fusion concepts, such as the magnetized liner inertial fusion (MagLIF) platform, constitute an alternative path for achieving ignition and significant fusion yields in the laboratory. The space of experimental input parameters defining a MagLIF load is highly multi-dimensional, and the implosion itself is a complex event involving many physical processes. In the first paper of this series, we develop a simplified analytical model that identifies the main physical processes at play during a MagLIF implosion. Using non-dimensional analysis, we determine the most important dimensionless parameters characterizing MagLIF implosions and provide estimates of such parameters using typical fielded or experimentally observed quantities for MagLIF. Here, we then show that MagLIF loads can be “incompletely” similarity scaled, meaning that the experimental input parameters of MagLIF can be varied such that many (but not all) of the dimensionless quantities are conserved. Based on similarity-scaling arguments, we can explore the parameter space of MagLIF loads and estimate the performance of the scaled loads. Then, in the follow-up papers of this series, we test the similarity-scaling theory for MagLIF loads against simulations for two different scaling “vectors,” which include current scaling and rise-time scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Advancing Thermal-Hydraulic Modeling Capabilities for Molten Salt Reactors in Pronghorn

This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.

42 ENGINEERING↗