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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗
NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗
NASA NTRS2004

Effects of Surfactant Contamination on the Next Generation Gas Trap for the ISS Internal Thermal Control System

The current dual-membrane gas trap is designed to remove non-condensed gas bubbles from the Internal Thermal Control System (ITCS) coolant on board the International Space Station (ISS). To date it has successfully served its purpose of preventing gas bubbles from causing depriming, overspeed, and shutdown of the ITCS pump. However, contamination in the ITCS coolant has adversely affected the gas venting rate and lifetime of the gas trap, warranting a development effort for a next-generation gas trap. Previous testing has shown that a hydrophobic-only design is capable of performing even better than the current dual-membrane design for both steady-state gas removal and gas slug removal in clean deionized water. This paper presents results of testing to evaluate the effects of surfactant contamination on the steady-state performance of the hydrophobic-only design.

Leimkuehler, Thomas O.↗
DOE OSTI2025

Validation of Pronghorn for Natural-Circulation Molten Salt Loops

This paper presents the development and validation of a high-fidelity thermal-hydraulic model of a molten salt natural circulation flow loop, designed for integration within a digital twin framework. The study evaluates the performance of Idaho National Laboratory’s Pronghorn against experimental data from Texas A&M University Molten Salt Flow Loop (MSFL) four Hitec-salt test benchmark data. Natural circulation of high-Prandtl-number fluids exhibits complex, counter-intuitive flow patterns that make pointwise thermocouple readings unreliable. Experimental work at TAMU’s MSFL provides benchmark data, including flow visualization at a test-section and centerline steady-state temperature measurements along the loop. Validation includes four single-phase natural circulation test cases with Hitec salt. Key metrics include flow profile agreement and steady-state temperature accuracy. Pronghorn results for two-dimensional single-phase agree qualitatively with the experimental flow profile. This paper illustrates the importance of Computational Fluid Dynamics (CFD) in elucidating the behavior of high-Prandtl-number thermal-hydraulics, along with how misleading centerline temperature measurements can be. Pronghorn reproduces the axial and radial stratification that makes single thermocouple readings unreliable. Future research will focus on reduced-order modeling techniques to enable rapid simulation suitable for real-time digital twin applications. The validated cases provide a basis for developing reduced-order surrogates aimed at real-time digital-twin applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Geometric decoherence time in Lindbladian dynamics

The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometric decoherence time, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-$\frac{1}{2}$ entropy—exactly equal across any bipartition for pure states—breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: Its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail—residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local 𝑍 dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Furthermore, our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Time-Resolved Structural Kinetics of an Organic Mixed Ionic–Electronic Conductor

The structure and packing of organic mixed ionic–electronic conductors have an especially significant effect on transport properties. In operating devices, this structure is not fixed but is responsive to changes in electrochemical potential, ion intercalation, and solvent swelling. Toward this end, the steady-state and transient structure of the model organic mixed conductor, poly(3,4-ethylenedioxythiophene) polystyrene sulfonate (PEDOT:PSS), is characterized using multimodal time-resolved operando techniques. Steady-state operando X-ray scattering reveals a doping-induced lamellar expansion of 1.6 Å followed by 0.4 Å relaxation at high doping levels. Time-resolved operando X-ray scattering reveals asymmetric rates of lamellar structural change during doping and dedoping that do not directly depend on potential or charging transients. Further, time-resolved spectroscopy establishes a link between structural transients and the complex kinetics of electronic charge carrier subpopulations, in particular the polaron–bipolaron equilibrium. These findings provide insight into the factors limiting the response time of organic mixed-conductor-based devices, and present the first real-time observation of the structural changes during doping and dedoping of a conjugated polymer system via X-ray scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Secure Inverters Against Malicious Setpoints

The next generation of grid-interactive inverters brings a communication feature that allows data sharing from utility supervisory controllers and smart devices that are connected to the same network. This feature enhances the control capabilities of grid-interactive inverters to provide services beyond active power injection. However, communication networks entail more vulnerable surfaces to malicious attacks that may result in modifying active and reactive power setpoints and causing weak-grid conditions or abnormal inverter operation. In this paper, steady-state and the dynamic behavior of the inverter for the incoming setpoints are analyzed to detect false data injection attacks and provide device-level security. The steady-state behavior of the inverter in the operating region is determined from the grid parameters such as the grid voltage and the grid impedance. These estimations are accomplished by the proposed self-security technique through a low-frequency signal injection-based approach combined with the recursive least square method. Moreover, a reduced fourth-order inverter model is used as the dynamic reference model, and grid parameters as well as the incoming setpoints are implemented to the reference model to verify whether the dynamic behavior of the inverter is inside the permissible region of operation. The validity and performance of the proposed method are verified experimentally through Allen-Bradley Powerflex 755 three-phase inverter and a 12 kW NHR 9410 regenerative power grid emulator. The results show that the self-secure smart-inverter is able to accept or reject the incoming commands and thus is protected from malicious cyber-physical attacks.

Hossen, Tareq↗
DOE OSTI2023

Two-point Correlation Function Studies for the Milky Way: Discovery of Spatial Clustering from Disk Excitations and Substructure

We introduce a two-particle correlation function (2PCF) for the Milky Way, constructed to probe spatial correlations in the orthogonal directions of the stellar disk in the Galactic cylindrical coordinates of R, ϕ, and z. We use this new tool to probe the structure and dynamics of the Galaxy using the carefully selected set of solar neighborhood stars (d ≲ 3 kpc) from Gaia Data Release 2 that we previously employed for studies of axial symmetry breaking in stellar number counts. We make additional, extensive tests, comparing to reference numerical simulations, to ensure our control over possibly confounding systematic effects. Supposing either axial or north–south symmetry, we divide this data set into two nominally symmetric sectors and construct the 2PCF, in the manner of the Landy–Szalay estimator, from the Gaia data. In so doing, working well away from the midplane region in which the spiral arms appear, we have discovered distinct symmetry-breaking patterns in the 2PCF in its orthogonal directions, thus establishing the existence of correlations in stellar number counts alone at subkiloparsec length scales for the very first time. In particular, we observe extensive wavelike structures of amplitude greatly in excess of what we would estimate if the system were in a steady state. We study the variations in these patterns across the Galactic disk, and with increasing ∣z∣, and we show how our results complement other observations of non-steady-state effects near the Sun, such as vertical asymmetries in stellar number counts and the Gaia snail.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Study of monopropellants for electrothermal thrusters: Design and fabrication task summary report

The feasibility of operating small thrust level electrothermal thrusters with monopropellants other than MIL-grade hydrazine was studied. Analytical study, design, and fabrication of demonstration thrusters was performed, and an evaluation test program was initiated to evaluate monopropellants with freezing points lower than MIL-grade hydrazine, and to determine their applicability to electrothermal thrusters for spacecraft attitude control. Five demonstration thrusters were fabricated to determine the feasibility of operation with monomethylhydrazine, Aerozine-50, 77 percent hydrazine-23 percent hydrazine azide, and a mixture of hydrazine monopropellants consisting of 35 percent hydrazine-50 percent monomethylhydrazine-15 percent ammonia. The present thruster is designed to produce a steady-state thrust level of 0.344 N at 1.724 x 1 million N/M sq feed pressure. Vacuum specific impulse goals were set at 1961 N-s/kg steady-state and 1716 N-s/kg pulsed-mode.

Kuenzly, J. D.↗
NASA NTRS1977

A study of the unification of ground and inflight wind criteria for the space shuttle

Wind data measured on a 444-m tower located near Oklahoma City, Oklahoma, were used to investigate the structure of wind and turbulence in the height interval from 150 m to 1 km. The structure of wind and turbulence in the layer not encompassed by the existing design criteria was researched and ways for unifying the ground and inflight criteria were investigated. Steady state vertical wind profiles, directional wind component envelopes, wind shear, wind direction change, gust factor, and turbulence spectra are encompassed. A method is proposed for the specification of steady state wind profiles and shear envelopes for use in the region between 150 m and 1 km without altering in any way the existing design criteria. The data analyzed did not indicate a need for changing the existing criteria in regards to wind direction change, gust factor, or spectra of turbulence.

Scoggins, J. R.↗
NASA NTRS1987

Nonplanar interface morphologies during unidirectional solidification of a binary alloy. II - Three-dimensional computations

A finite difference method is used to obtain three-dimensional steady-state solutions for nonplanar interface morphologies in order to study the situation of equal thermal properties in the crystal and melt with negligible latent heat release. Stable steady-state solutions corresponding to two-dimensional bands and three-dimensional hexagonal nodes, as well as to rectangular interface planiforms, are found using a model of an aluminum-chromium alloy with a distribution coefficient of greater than one. Hexagonal nodes are predicted near the onset of instability, in agreement with weakly nonlinear theory.

Mcfadden, G. B.↗
NASA NTRS1988

A test of synthetic accretion disk spectra using ultraviolet flux distributions of novalike variables

Ultraviolet (UV) fluxes and other data for a sample of nine novalike cataclysmic variables are assembled from the literature. The UV fluxes and colors are compared to spectra of steady state model accretion disks, constructed using either Planck functions or model stellar atmosphere spectra as elementary radiators. Deficiencies are found with both sets of models. The restriction to steady state radial temperature profiles is relaxed for the 'stellar atmosphere' disks, but these generalized models still fail to account simultaneously for the absolute UV flux and the UV color of the observed spectra. The conclusion is that stellar atmosphere models do not reflect the physics of disks. More appropriate models probably can and definitely should be made.

Wade, Richard A.↗
NASA NTRS1994

Lifetimes and eigenstates in atmospheric chemistry

The time scales and mode of the atmosphere's response to chemical perturbations are defined by the eigenvalues and eigenvectors of the system. The eigenstates of a simplified one-box CH4-CO-OH system are analyzed. The longest time constant (smallest eigenvalue) always exceeds the lifetime defined by the steady-state loss frequency for CH4, the longest lived gas. Thus, the extent of a CH4 perturbation -- the methane response time -- is always longer than predicted by the steady-state lifetime and is independent of size of the perturbation in the linear limit. This lengthening of the atmospheric recovery time can be diagnosed by how close we are to a chemically unstable troposphere, i.e., how much OH production exceeds that minimum needed to oxidize just the global emissions of CH4, CO, and other hydrocarbons and species.

Prather, Michael J.↗
NASA NTRS2012

Heat Transfer Modeling for Rigid High-Temperature Fibrous Insulation

Combined radiation and conduction heat transfer through a high-temperature, high-porosity, rigid multiple-fiber fibrous insulation was modeled using a thermal model previously used to model heat transfer in flexible single-fiber fibrous insulation. The rigid insulation studied was alumina enhanced thermal barrier (AETB) at densities between 130 and 260 kilograms per cubic meter. The model consists of using the diffusion approximation for radiation heat transfer, a semi-empirical solid conduction model, and a standard gas conduction model. The relevant parameters needed for the heat transfer model were estimated from steady-state thermal measurements in nitrogen gas at various temperatures and environmental pressures. The heat transfer modeling methodology was evaluated by comparison with standard thermal conductivity measurements, and steady-state thermal measurements in helium and carbon dioxide gases. The heat transfer model is applicable over the temperature range of 300 to 1360 K, pressure range of 0.133 to 101.3 x 10(exp 3) Pa, and over the insulation density range of 130 to 260 kilograms per cubic meter in various gaseous environments.

Daryabeigi, Kamran↗
DOE OSTI2026

Data for KETCHUP: Parameterizing of Large-Scale Kinetic Models Using Multiple Datasets with Different Reference States

Repository for Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP can use K-FIT [2] input files. Example K-FIT input files are located in the K-FIT repository at https://github.com/maranasgroup/K-FIT.

Metabolomics↗
DOE OSTI2024

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a "multizone" technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH ($a\ast$ = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r –1 inside the Bondi radius R B , which is located at R B = 2 × 10 5 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ~ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is ≈1% of the analytical Bondi accretion rate ${\dot{M}}_{{\rm{B}}};$ and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate $\approx 0.02\dot{M}{c}^{2}$. Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ~ 10 2 –10 7 r g , and find that $\dot{M}/{\dot{M}}_{{B}}\approx {({R}_{{B}}/6\,{r}_{g})}^{-0.5}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗
NASA NTRS1972

Gas-core reactor power transient analysis.

The nuclear fuel in the gas-core reactor concept is a ball of uranium plasma radiating thermal photons. The photons are met by an inflowing hydrogen stream, which is seeded with submicron size, depleted uranium particles. A 'wall-burnout' condition exists if the thermal photons can reach the cavity liner because of insufficient absorption by the hydrogen. An analysis was conducted in order to determine the time for which the maximum steady state reactor power could be exceeded without damage to the cavity liner due to burnout. Wall-burnout time as a function of the power increase above the initial steady state condition is shown in a graph.

Kascak, A. F.↗