Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STEADY STATE”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Instability of steady-state mixed-state symmetry-protected topological order to strong-to-weak spontaneous symmetry breaking

Recent experimental progress in controlling open quantum systems enables the pursuit of mixed-state nonequilibrium quantum phases. We investigate whether open quantum systems hosting mixed-state symmetry-protected topological states as steady states retain this property under symmetric perturbations. Focusing on the decohered cluster state – a mixed-state symmetry-protected topological state protected by a combined strong and weak symmetry – we construct a parent Lindbladian that hosts it as a steady state. This Lindbladian can be mapped onto exactly solvable reaction-diffusion dynamics, even in the presence of certain perturbations, allowing us to solve the parent Lindbladian in detail and reveal previously-unknown steady states. Using both analytical and numerical methods, we find that typical symmetric perturbations cause strong-to-weak spontaneous symmetry breaking at arbitrarily small perturbations, destabilize the steady-state mixed-state symmetry-protected topological order. However, when perturbations introduce only weak symmetry defects, the steady-state mixed-state symmetry-protected topological order remains stable. Additionally, we construct a quantum channel which replicates the essential physics of the Lindbladian and can be efficiently simulated using only Clifford gates, Pauli measurements, and feedback.

Shah, Jeet [University of Maryland, College Park, ↗
DOE OSTI2025

HFIR LEU High Density Silicide Dispersion Optimized Design Steady-State Heat Transfer Analyses

Steady-state heat transfer simulations of the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) with the low-enriched uranium (LEU) high-density silicide dispersion Optimized fuel design were performed to support comprehensive performance and safety metric studies concerning this design. The LEU Optimized design operates at 95 MW to maintain HFIR’s current highly enriched uranium (HEU) core performance level at 85 MW. Full cycle Mode 1 full flow Case 1 (inlet temperature), Case 2 (flux-to-flow), and Case 3 (inlet pressure) safety limit analyses were performed to assess the margins to critical heat flux. Under the prescribed conditions, this LEU design meets the safety limit and limiting control setting requirements outlined in HFIR’s documented safety analysis; however, the safety margins are less than those for the 85 MW HEU core, and several assumptions were made where fuel fabrication and qualification data are currently lacking for the silicide fuel design. Effects of changes to pertinent fuel fabrication assumptions and uncertainty factors on thermal safety margins were also evaluated, showing that the margins are sensitive to many of these parameters. Power and pressure perturbations were also performed, indicating that significant steady-state thermal margins could be gained by increasing the coolant inlet pressure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Generating entangled steady states in multistable open quantum systems via initial state control

Entanglement underpins the power of quantum technologies, yet it is fragile and typically destroyed by dissipation. Paradoxically, the same dissipation, when carefully engineered, can drive a system toward robust entangled steady states. However, this engineering task is nontrivial, as dissipative many-body systems are complex, particularly when they support multiple steady states. Here, we derive analytic expressions that predict how the steady state of a system evolving under a Lindblad equation depends on the initial state, without requiring integration of the dynamics. These results extend Refs. [V. V. Albert and L. Jiang, Phys. Rev. A 89, 022118 (2014); V. V. Albert et al., Phys. Rev. X 6, 041031 (2016)], showing that while the steady-state manifold is determined by the Liouvillian kernel, the weights within it depend on both the Liouvillian and the initial state. We identify a special class of Liouvillians for which the steady state depends only on the initial overlap with the kernel. Our framework provides analytical insight and a computationally efficient tool for predicting steady states in open quantum systems. As an application, we propose schemes to generate metrologically useful entangled steady states in spin ensembles via balanced collective decay.

Dissipative dynamics↗
DOE OSTI2025

Efficient Simulation of Cascading Outages Using an Energy Function-Embedded Quasi-Steady-State Model

Here, this paper proposed an energy function-embedded quasi-steady-state model for efficient simulation of cascading outages on a power grid while addressing transient stability concerns. Compared to quasi-steady-state models, the proposed model incorporates short-term dynamic simulation and an energy function method to efficiently evaluate the transient stability of a power grid together with outage propagation without transient stability simulation. Cascading outage simulation using the proposed model conducts three steps for each disturbance such as a line outage. First, it performs time-domain simulation for a short term to obtain a post-disturbance trajectory. Second, along the trajectory, the system state with the local maximum potential energy is found and used as the initial point to search for a relevant unstable equilibrium by Newton's method. Third, the transient energy margin is estimated based on this unstable equilibrium to predict an out-of-step condition with generators. The proposed energy function-embedded quasi-steady-state model is tested in terms of its accuracy and time performance on an NPCC 140-bus power system and compared to a quasi-steady-state model embedding transient stability simulation.

Guo, Zhenping [Univ. of Tennessee, Knoxville, TN (↗
DOE OSTI2026

FAST-1.2.2: A Computer Code for Thermal-Mechanical Nuclear Fuel Analysis under Steady-state and Transients

Fuel Analysis under Steady-state and Transients (FAST) is the U.S. Nuclear Regulatory Commission (NRC)’s computer code that calculates the steady-state and transient response of nuclear reactor fuel rods during long-term in-reactor burnup, anticipated operational occurrences (AOOs), design basis accidents (DBAs), and dry storage conditions. The code calculates the temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include heat conduction through the fuel and other materials, heat transfer from the cladding-to-coolant, cladding elastic and plastic deformation (including creep), fuel-cladding mechanical interaction, fission gas release from the fuel, rod internal pressure, void volume, and cladding oxidation. The code contains necessary material and coolant properties, as well as clad-to-coolant heat transfer correlations, for normal operation through postulated accidents and AOOs for today’s U.S.-based light water reactor (LWR) fuel designs. FAST-1.2.2 also contains preliminary materials and models for new LWR fuel concepts, such as accident tolerant fuel (ATF), and non-LWR fuel concepts such as metallic fuels for sodium fast reactors (SFRs). FAST has been developed for use on Windows and Linux operating systems. This document describes FAST-1.2.2 and is one of a series of documents on the code; the other documents detail the material properties used by FAST as well as its integral assessment to experiments and commercial data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Initial steady-state core simulation capability or thermal and pool-type molten salt reactors, coupling reactor physics, thermal-hydraulics, and evolving chemistry

This report presents the development and validation of an initial steady-state multiphysics capability for molten salt reactors (MSRs) under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in Fiscal Year 2025. The framework integrates neutronics, thermal-hydraulics, species transport, and thermochemistry to capture the coupled dynamics of liquid-fueled systems. Implementation and testing were performed on two representative designs: the Molten Salt Reactor Experiment (MSRE), a thermal-spectrum, channeled-flow reactor, and the Lotus Molten Salt Reactor (L-MSR), a fast-spectrum, pool-type reactor. The modeling suite employs Griffin for reactor physics and depletion, Pronghorn and SAM for thermal-hydraulics, Thermochimica for chemistry, and Saline for thermophysical properties, with benchmarking and validation carried out against historical MSRE data, experimental flow-loop measurements, and reference depletion calculations from Monte Carlo codes. The framework demonstrated the ability to reproduce key reactor behaviors including temperature feedback, reactivity losses, delayed neutron precursor transport, xenon poisoning, and redox potential evolution. The results confirm the feasibility and accuracy of the coupled models in predicting steady-state and selected transient MSR behaviors. This latter ones are used in this report as a proxy indicating that the steady-state models from which the transient starts are accurate. For MSRE, validation showed good agreement with pump start-up and natural circulation tests, while for the L-MSR, benchmarking confirmed hydraulic calibration and consistency of neutronics–thermal coupling. The tools also provided new insights into species transport, noble metal deposition, and salt solidification dynamics. On the Xenon transport front, the code is validated against the steady state Xenon poisoining measurement and showed good agreement with the experimental value. Identified areas for future work include advanced void transport modeling, three-dimensional simulations, improved alloy corrosion models, and tighter integration with high-fidelity Monte Carlo codes. These developments provide a foundation for high-fidelity MSR simulations that can support reactor design optimization, safety assessments, and long-term operational strategies.

42 - ENGINEERING↗
DOE OSTI2025

Steady-state mixing state of black carbon aerosols from a particle-resolved model

Abstract. Black carbon (BC) exerts a notable warming effect due to its strong light absorption, largely influenced by its “mixing state”. However, due to computational constraints, the mixing state is challenging to accurately represent in large-scale models. In this study, we employ a particle-resolved model to simulate the evolution of BC mixing state based on field observation. Our result shows that aerosol compositions, coating thickness (CT) distribution, and optical properties of BC aerosols all exhibit a tendency toward a steady state with a characteristic timescale of less than 1 d, considerably shorter than the BC atmospheric lifetime. The rapid attainment of a steady state suggests that it is reasonable to disregard this pre-steady-state period and instead concentrate on the average properties of BC across extensive spatial and temporal scales. The distribution of CT follows an exponential linear distribution and can be characterized by a single slope parameter k. This distribution is independent of the BC core's distribution. In the model simulation, the mean CT, equivalent to the 1/k, is 62 nm, which is consistent with the statistical results indicating a mean CT of 63 nm. Utilizing the slope parameter k, which effectively characterizes the CT distribution under the steady-state simplifying assumption, the BC absorption enhancement closely corresponds to the results obtained via the particle-resolved method. This study simplifies the BC mixing state description and yields a precise evaluation of the BC optical properties, which has the potential utility for modeling efforts in the refinement of the assessment of BC's radiative effects.

Zhang, Zhouyang↗
DOE OSTI2025

Hierarchical Network Partitioning for Solution of Potential-Driven, Steady-State Nonlinear Network Flow Equations

The solution of potential-driven steady-state flow in large networks is a task which manifests in various engineering applications, such as transport of natural gas or water through pipeline networks. The resultant system of nonlinear equations depends on the network topology, and in general, there is no numerical algorithm that offers guaranteed convergence to the solution (assuming a solution exists). Some methods offer guarantees in cases where the network topology satisfies certain assumptions, but these methods fail for larger networks. On the other hand, the Newton-Raphson algorithm offers a convergence guarantee if the starting point lies close to the (unknown) solution. It would be advantageous to compute the solution of the large nonlinear system through the solution of smaller nonlinear sub-systems wherein the solution algorithms (Newton-Raphson or otherwise) are more likely to succeed. Here, this letter proposes and describes such a procedure, a hierarchical network partitioning algorithm that enables the solution of large nonlinear systems corresponding to potential-driven steady-state network flow equations.

42 ENGINEERING↗
DOE OSTI2025

Alfvén mode spectrum in the advanced steady-state tokamak scenario with negative triangularity

In the earlier studies [Zheng et al. Fundam. Plasma Phys. 10, 100051 (2024)], the negative triangularity tokamak was shown to have a unique and outstanding feature for the steady-state confinement. In addition to relieving the divertor heat load concern, it can reach a much higher beta limit in the nearly steady-state confinement, about twice the Troyon limit in the positive triangularity case. In this work, we extend the investigation to the Alfvén spectrum and toroidal Alfvén eigenmodes (TAEs). This helps to evaluate the property of energetic particle confinement in this advanced scenario of a negative triangularity tokamak. We found that, because of the large safety factor and reversed magnetic shear in the core plasma region near the magnetic axis, the TAE gap tends to be closed or have only a short radial extension between the continua. in the on-axis core region and open primarily in the off-axis region, where TAEs are shown to be present. This indicates that energetic particle confinement in the core region near the magnetic axis, where nuclear fusion reactions prevail, is not seriously affected by the excitation of TAEs in this scenario. The configuration with off-axis TAEs helps the confinement of energetic particles as compared to the configuration with near-axis TAEs. The results show that the advanced steady-state tokamak scenario with negative triangularity, because of the rapid variation of the safety factor, also has a favorable feature for energetic particle confinement.

Fourier analysis↗
DOE OSTI2025

Comparative post-irradiation examination of high burnup U-19Pu-10Zr: Assessing steady-state irradiation behavior against historical and modeled fuel performance

Here, the development of next-generation sodium-cooled fast reactors necessitates comprehensive research on metallic fuels to maximize economic performance while ensuring safe operation. In this study, we investigated the steady-state irradiation behavior of two high burnup U-19Pu-10Zr fuel pins, DP-36 and DP-40, in preparation for planned safety testing. Post-irradiation examination (PIE) was performed to quantify fuel column elongation, regions of low-density at the top of the fuel column, pin deformation, fission product distribution, fractional fission gas release, microstructural evolution, and fuel constituent redistribution. Benchmarking against existing PIE data from U-19Pu-10Zr fuel pins irradiated in EBR-II revealed consistent patterns in fuel column elongation and cladding diametral strain. However, both pins exhibited longer low-density structures, and destructive examination of DP-36 revealed more complex constituent redistribution patterns compared to previously reported data for ternary fuel pins. The steady-state irradiation of both pins was also modeled using BISON. Comparisons of PIE results with modeled predictions showed overall agreement in fractional fission gas release but consistent overestimation of axial and radial swelling due to gaseous and solid swelling models. These findings underscore the critical importance of pre-test characterization on test and sibling pins to accurately capture steady-state fuel behavior ahead of transient testing, thus establishing a baseline for post-test comparison. Additionally, these analyses identified key data gaps that warrant further investigation to improve the understanding and prediction of fuel swelling, thereby enhancing the synergy between modeling and experimental efforts in supporting accident testing.

BISON↗
DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗
DOE OSTI2024

Symbolic construction of the chemical Jacobian of quasi-steady state (QSS) chemistries for Exascale computing platforms

The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of chemical mechanisms for implementation in computational reacting flow solvers. However, for many applications, the resulting model still requires implicit methods for efficient time integration. Here, in this paper, we outline an approach to formulating the QSSA reduction that is coupled with a strategy to generate C++ source code to evaluate the net species production rates, and the chemical Jacobian. The code-generation component employs a symbolic approach enabling a simple and effective strategy to analytically compute the chemical Jacobian. For computational tractability, the symbolic approach needs to be paired with common subexpression elimination which can negatively affect memory usage. Several solutions are outlined and successfully tested on a 3D multipulse ignition problem, thus allowing portable application across chemical model sizes and GPU capabilities. The implementation of the proposed method is available at https://github.com/AMReX-Combustion/PelePhysics under an open-source license.

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

Electrochemical Modeling of PID Leakage Current of PV Modules: Steady-State Current, Transient Current, and RC-Equivalent Circuit

Potential-induced degradation (PID) remains a significant reliability concern for photovoltaic (PV) modules, arising when a voltage difference between the module frame and the solar cells drives unintended leakage current through the glass-encapsulant stack. Although PID ultimately manifests as PID-s, PID-p, or PID-c, the underlying behavior of the leakage current-its magnitude and time dependence-requires clearer electrochemical interpretation. Traditional explanations attribute the initial transient current to bulk capacitive elements of the glass, encapsulant, and antireflection coatings, and the steady-state current according to their effective ohmic resistance. More recent studies, however, indicate that electrochemical charge-transfer processes at the encapsulant-metallization interface can play a dominant role in defining the leakage-current path. This paper develops a unified electrochemical framework for modeling PID leakage current. First, an RC-equivalent circuit is formulated by combining conventional RC elements with a Randles-type interface to capture transient leakage current through double-layer capacitance and faradaic processes at ionic-electronic boundaries. Second, the steady-state current-voltage behavior is explained using a linearized Butler-Volmer relationship, showing that the measured ohmic response corresponds to the low-overpotential limit of charge-transfer kinetics. Analytical results demonstrate that, for typical module materials-3.2-mm soda-lime glass and 0.45-mm encapsulant-the dominant modulators to PID leakage current are the glass surface resistance (under dry-surface conditions), the glass bulk capacitance, and the encapsulant resistance (under wet-surface conditions), with soda lime glass surface and EVA/POE encapsulant resistances primarily governing steady-state current. The proposed electrochemical model is validated against measured leakage-current data, showing good agreement in both the magnitude and the time-dependent evolution of PID leakage current.

14 SOLAR ENERGY↗
DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2025

Steady state radiolysis–polyhalide accumulation products [Poster]

The study of steady-state radiolysis of polyhalides is crucial as it simulates the prolonged radiation exposure that materials (e.g., alloys, molten salts) experience in reactor environments. This research provides valuable insights into the long-term effects on material properties, chemical stability, and structural integrity. It is closely linked with fundamental radiation chemistry, which aids in predicting the chemical behavior of accumulation products. Understanding the initial speciation, yield, and chemical evolution of radiolysis products is essential for managing material degradation, controlling corrosion, and ensuring reactor safety.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗