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 793 records · Page 44

Low-dimensional Representation of Error Covariance

Ensemble and reduced-rank approaches to prediction and assimilation rely on low-dimensional approximations of the estimation error covariances. Here stability properties of the forecast/analysis cycle for linear, time-independent systems are used to identify factors that cause the steady-state analysis error covariance to admit a low-dimensional representation. A useful measure of forecast/analysis cycle stability is the bound matrix, a function of the dynamics, observation operator and assimilation method. Upper and lower estimates for the steady-state analysis error covariance matrix eigenvalues are derived from the bound matrix. The estimates generalize to time-dependent systems. If much of the steady-state analysis error variance is due to a few dominant modes, the leading eigenvectors of the bound matrix approximate those of the steady-state analysis error covariance matrix. The analytical results are illustrated in two numerical examples where the Kalman filter is carried to steady state. The first example uses the dynamics of a generalized advection equation exhibiting nonmodal transient growth. Failure to observe growing modes leads to increased steady-state analysis error variances. Leading eigenvectors of the steady-state analysis error covariance matrix are well approximated by leading eigenvectors of the bound matrix. The second example uses the dynamics of a damped baroclinic wave model. The leading eigenvectors of a lowest-order approximation of the bound matrix are shown to approximate well the leading eigenvectors of the steady-state analysis error covariance matrix.

Tippett, Michael K.↗

Computational Study of Axisymmetric Off-Design Nozzle Flows

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles operating off-design at transonic Mach numbers have been completed. These computations span the very difficult transonic flight regime with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined, including the Explicit Algebraic Stress model. Computations have been completed at freestream Mach numbers of 0.9 and 1.2, and nozzle pressure ratios (NPR) of 4 and 6. Calculations completed with variable time-stepping (steady-state) did not converge to a true steady-state solution. Calculations obtained using constant timestepping (timeaccurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was the result of using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show reasonable agreement with experimental data. The SST turbulence model demonstrates the best overall agreement with experimental data.

DalBello, Teryn↗

Computations of Internal and External Axisymmetric Nozzle Aerodynamics at Transonic Speeds

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles have been completed in support of NASA's Next Generation Launch Technology Program to investigate the effects of high-speed nozzle geometries on the nozzle internal flow and the surrounding boattail regions. These computations span the very difficult transonic flight regime, with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined in Wind, including an Explicit Algebraic Stress model (EASM). Computations on two nozzle geometries have been completed at freestream Mach numbers ranging from 0.6 to 0.9, driven by nozzle pressure ratios (NPR) ranging from 2.9 to 5. Results obtained on converging-only geometry indicate reasonable agreement to experimental data, with the EASM and Shear Stress Transport (SST) turbulence models providing the best agreement. Calculations completed on a converging-diverging geometry involving large-scale internal flow separation did not converge to a true steady-state solution when run with variable timestepping (steady-state). Calculations obtained using constant timestepping (time-accurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was found to be the result of difficulties in using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show fairly good agreement with experimental data. The SST turbulence model demonstrates the best over-all agreement with experimental data.

Dalbello, Teryn↗

Comparisons Between Unsteady Aerodynamic Events in a Gas Turbine Generator and an Identical Compressor Rig

Extensive testing done on a T55-L-712 turboshaft engine compressor in a compressor test rig is being followed by engine tests in progress as part of the Army Non-Recoverable Stall Program. Goals include a greater understanding of the gas turbine engine start cycle and compressor/engine operation in the regions 'beyond' the normal compressor stall line (rotating stall/surge). Rig steady state instrumentation consisted of 497 steady state pressure sensors and 153 temperature sensors. Engine instrumentation was placed in similar radial/axial locations and consists of 122 steady state pressure sensors and 65 temperature sensors. High response rig instrumentation consisted of 34 wall static pressure transducers. Rig and engine high response pressure transducers were located in the same axial/radial/circumferential locations in front of the first three stages. Additional engine high response instrumentation was placed in mach probes in front of the engine and on the compressor hub. This instrumentation allows for the generation of detailed stage characteristics, overall compressor mapping, and detailed analysis of dynamic compressor events.

Owen, A. Karl↗

The Monolithic Heat Pipe Microreactor Reference Plant Model

This work introduces a reference plant model for a generic monolithic heat-pipe-cooled microreactor. The model will serve as a springboard to develop future evaluation models in the licensing process of similar microreactor designs at the U.S. Nuclear Regulatory Commission. This model has been developed with the Comprehensive Reactor Analysis Bundle and its specifications are based on open literature publications for the eVinci TM design. BlueCRAB is the U.S. Nu- clear Regulatory Commission non-light-water reactor analysis system based on the Multiphysics Object-Oriented Simulation Environment framework, which can couple the Griffin, BISON, and Sockeye applications to resolve the various physics that are essential for the safety analysis of this type of reactor system. The core specifications includes tristructural isotropic fuel, graphite monolith, graphite reflectors, and drums composed of graphite and B 4 C. No moderator or burnable poison pins are used in the design. The fuel enrichment is reduced to control excess reactivity in the core. This core design is not optimized and only serves for testing purposes, since the primary objective of this work is to exercise the multiphysics coupling for this type of reactor system. A three dimensional (3D) core heterogeneous Griffin discrete ordinates (SN) transport model allows the precise calculation of the flux distribution and pin powers. Griffin transfers the power density distribution and obtains a temperature distribution to and from BISON. The BISON model com- putes the 3D core temperature distribution and is coupled to 876 Sockeye subapplications running a heat pipe model. This 3D conduction model is coupled to the various heat pipes via heat flux boundary conditions. The model includes a small gap between the heat pipe and the monolith. Convective heat transfer boundaries with either ambient temperature or condenser temperature as heat sinks are imposed at the model boundaries. The 2D Sockeye heat pipe model uses a vapor- only methodology, which provides the needed resolution for transient calculations and allows the determination of various heat pipe limits. This approach is superior to the superconductor model traditionally used in steady-state calculations. BlueCRAB computes steady-state power and temperature distributions that serve as the initial condition for a loss-of-heat-sink transient simulation. The steady-state results show significant peaking due to the position of the control drum, but this is a characteristic of the particular design used, which is not optimized at this stage. The transient results show the reactor power slowly stabilizing towards a 3% power level after the partial loss of secondary heat removal. Several recriticalities are observed due to cooling through the secondary system but the reactor is self-stabilizing and behaves as expected.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Preliminary Thermal-hydraulics and Selected Safety Analysis for Holos-Quad Reactor Design

This report presents thermal hydraulics and safety analyses of the 10 MWe Holos-Quad micro-reactor design developed by HolosGen LLC. These analyses were executed under steady-state and station black-out (SBO) transient conditions using the System Analysis Module (SAM), with a focus on ensuring peak fuel temperatures do not exceed the safety thresholds of 1,250°C for steady state normal operating conditions and 1,600°C for transient conditions. For steady-state temperature predictions, a coupled 1D fluid-to-3D solid heat conduction methodology was used. The steady-state simulations included a model of the central core region where power peaking occurs to predict the maximum steady-state fuel temperature, and a coarse-mesh full-core model to provide the initial condition for the consequent transient SBO decay heat removal simulation. The maximum fuel temperature predicted from the central core region simulations under the normal operating condition was 1,222°C, which is below the design limit 1,250°C. To simulate a SBO transient, a full core model was necessary because of the non-axial symmetric core design and heat loss from the outer core structures to the environment. To reduce the required computational costs, a novel scaled subassembly approach was used to model the full core domain. With the full core steady-state temperature distribution provided as the initial condition, the passive decay heat removal of the core was simulated with the scaled subassembly full core model under SBO transient conditions. This analysis showed that the peak fuel temperatures of the Holos-Quad core remained below their steady-state values during the transient, therefore much lower than the safety limit 1,600°C. This is mainly due to the large thermal inertia of the graphite matrix, relatively low power density, and the large surface-to-volume-ratio of the core. The preliminary analyses presented in this report confirm the inherent safety of the Holos-Quad micro-reactor thermal-hydraulics design, as the peak fuel temperatures under both normal and off-normal operational conditions remain within design and safety limits.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The compensated Kalman filter.

This paper introduces the compensated Kalman filter, a suboptimal state estimator which can be used to eliminate steady-state bias errors when it is used in conjunction with the mismatched steady-state (asymptotic) time-invariant Kalman-Bucy filter. The uncompensated mismatched steady state Kalman-Bucy filter exhibits bias errors whenever the nominal plant parameters used in the filter design are different from the actual plant parameters. The approach used relies on the utilization of the residual (innovations) process of the mismatched filter to estimate, via a Kalman-Bucy filter, the state estimation errors and subsequent improvements of the state estimate. The compensated Kalman filter augments the mismatched steady state Kalman-Bucy filrby the introduction of additional dynamics and feedforward integral compensation channels.

Athans, M.↗

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗

Noise and Dynamical Pattern Selection

In pattern-forming systems, such as Rayleigh-Benard convection or directional solidification, a large number of linearly stable, patterned steady states exist when the basic, simple steady state is unstable. Which of these steady states will be realized in a given experiment appears to depend on unobservable details of the system's initial conditions. We show, however, that weak, Gaussian white noise drives such a system toward a preferred wave number which depends only on the system parameters and is independent of initial conditions. We give a prescription for calculating this wave number, analytically near the onset of instability and numerically otherwise.

Kurtze, Douglas A.↗

Tests of an alternating current propulsion subsystem for electric vehicles on a road load simulator

The test results of a breadboard version of an ac electric-vehicle propulsion subsystem are presented. The breadboard was installed in the NASA Lewis Research Center Road Load Simulator facility and tested under steady-state and transient conditions. Steady-state tests were run to characterize the system and component efficiencies over the complete speed-torque range within the capability of the propulsion subsystem in the motoring mode of operation. Transient tests were performed to determine the energy consumption of the breadboard over the acceleration and cruise portions of SAE J227 and driving schedules B, C, and D. Tests in the regenerative mode were limited to the low-gear-speed range of the two speed transaxle used in the subsystem. The maximum steady-state subsystem efficiency observed for the breadboard was 81.5 percent in the high-gear-speed range in the motoring mode, and 76 percent in the regenerative braking mode (low gear). The subsystem energy efficiency during the transient tests ranged from 49.2 percent for schedule B to 68.4 percent for Schedule D.

Stenger, F. J.↗

Overview of Wendelstein 7-X high-performance operation

The Wendelstein 7-X (W7-X) stellarator has completed two consecutive experimental campaigns OP 2.2 (Sep.-Dec. 2024) and OP 2.3 (Feb.-May 2025) under a new operational strategy enabling more than one year of uninterrupted device availability. This approach, supported by exceptionally high subsystem reliability, allowed sustained high-efficiency plasma operations with up to 80–100 discharges per day across a broad range of magnetic configurations. Several key technical upgrades-most notably the first operation of a 1.5 MW class steady-state gyrotron, a new steady-state pellet injector, and advanced real-time feedback control systems significantly enhanced heating, fueling, and plasma control capabilities. Together, these improvements enabled major advances in long-pulse performance, high-β operation, and confinement optimization. Long-pulse discharges achieved 1.8 GJ of injected energy under fully detached divertor conditions, while reduced-field scenarios facilitated record volume-averaged β values approaching 3%. High-performance plasmas with centrally peaked density profiles, created via neutral beam injection (NBI) or sustained pellet fueling, demonstrated strongly reduced turbulent transport and stellarator-record fusion triple products. Complementary studies of power exhaust and divertor heat loads revealed the role of scrape-off-layer drift physics in shaping strike-line patterns under attached conditions. Together, the results from OP 2.2 and OP 2.3 significantly expand the operational space of W7-X and strengthen its role as a leading platform for steady-state stellarator research and reactor-relevant plasma scenarios.

Wendelstein 7-X↗

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗

Design of Experiments for Dynamic Test Runs in Solvent-Based CO 2 Capture Pilot Plants

Test runs in the pilot plants consume significant resources, and therefore, the learning from test runs should be maximized. Test runs conducted in the pilot plants are often steady state. It takes several hours for reaching steady-state in the pilot plants, and thus, the duration of the test runs needs to be long even for collecting few steady-state data points. On the other hand, a large number of measurements can be collected through dynamic test runs in a short span of time. This paper presents a systematic design of dynamic experiments (DoDEs) for identifiability of model parameters, which is achieved by persistently exciting the inputs signals. A pseudorandom binary sequence (PRBS) is designed as the input signal for DoDE due to its efficiency in obtaining sufficient spectral content. However, due to the long sequence size of the PRBS signal, a Schroeder-phase input signal, which is a multisine signal, is also designed. Tests for both types of signals are run in the Pilot Solvent Test Unit (PSTU) at the National Carbon Capture Center in Wilsonville, Alabama. The transient data are used to solve dynamic data reconciliation and parameter estimation problem. The estimated parameters are found to be not only superior to those estimated from using data collected from hundreds of steady-state test runs in a nonreactive (air–water) system, but the parameters could be estimated by using the dynamic data collected for about 24 h from the pilot plant for the MEA-H 2 O–CO 2 system.

CO2 capture↗

Ignition and subsequent transition to flame spread in a microgravity environment

The fire safety strategy in a spacecraft is (1) to detect any fire as early as possible, (2) to keep any fire as small as possible, and (3) to extinguish any fire as quickly as possible. This suggests that a material which undergoes a momentary, localized ignition might be tolerable but a material which permits a transition to flame spread would significantly increase the fire hazard. Therefore, it is important to understand how the transition from localized ignition to flame spread occurs and what parameters significantly affect the transition. The fundamental processes involved in ignition and flame spread have been extensively studied, but they have been studied separately. Some of the steady state flame models start from ignition to reach a steady state, but since the objective of such a calculation is to obtain the steady state flame spread rate, the calculation through the transition process is made without high accuracy to save computational time. We have studied the transition from a small localized ignition at the center of a thermally thin paper in a microgravity environment. The configuration for that study was axisymmetric, but more general versions of the numerical scheme have been developed by including the effects of a slow, external flow in both two and three dimensions. By exploiting the non-buoyant nature of the flow, it is possible to achieve resolution of fractions of millimeters for 3D flow domains on the order of 10 centimeters. Because the calculations are time dependent, we can study the evolution of multiple flame fronts originating from a localized ignition source. The interaction of these fronts determines whether or not they will eventually achieve steady state spread. Most flame spread studies in microgravity consider two-dimensional flame spread initiated by ignition at one end of a sample strip with or against a slow external flow. In this configuration there is only one flame front. A more realistic scenario involves separate, oppositely directed fronts in two dimensions, or a continuous, radially directed front in three dimensions. We present here some results of both the two and three dimensional codes.

Kashiwagi, Takashi↗

Current Pulses Momentarily Enhance Thermoelectric Cooling

The rates of cooling afforded by thermoelectric (Peltier) devices can be increased for short times by applying pulses of electric current greater than the currents that yield maximum steady-state cooling. It has been proposed to utilize such momentary enhancements of cooling in applications in which diode lasers and other semiconductor devices are required to operate for times of the order of milliseconds at temperatures too low to be easily obtainable in the steady state. In a typical contemplated application, a semiconductor device would be in contact with the final (coldest) somewhat taller stage of a multistage thermoelectric cooler. Steady current would be applied to the stages to produce steady cooling. Pulsed current would then be applied, enhancing the cooling of the top stage momentarily. The principles of operation are straightforward: In a thermoelectric device, the cooling occurs only at a junction at one end of the thermoelectric legs, at a rate proportional to the applied current. However, Joule heating occurs throughout the device at a rate proportional to the current squared. Hence, in the steady state, the steady temperature difference that the device can sustain increases with current only to the point beyond which the Joule heating dominates. If a pulse of current greater than the optimum current (the current for maximum steady cooling) is applied, then the junction becomes momentarily cooled below its lowest steady temperature until thermal conduction brings the resulting pulse of Joule heat to the junction and thereby heats the junction above its lowest steady temperature. A theoretical and experimental study of such transient thermoelectric cooling followed by transient Joule heating in response to current pulses has been performed. The figure presents results from one of the experiments. The study established the essential parameters that characterize the pulse cooling effect, including the minimum temperature achieved, the maximum temperature overshoot, the time to reach minimum temperature, the time while cooled, and the time between pulses. It was found that at large pulse amplitude, the amount of pulse supercooling is about a fourth of the maximum steady-state temperature difference. For the particular thermoelectric device used in one set of the experiments, the practical optimum pulse amplitude was found to be about 3 times the optimum steady-state current. In a further experiment, a pulse cooler was integrated into a small commercial thermoelectric threestage cooler and found to provide several degrees of additional cooling for a time long enough to operate a semiconductor laser in a gas sensor.

Snyder, G. Jeffrey↗

A Dynamic Model for the Interaction Between an Insoluble Particle and an Advancing Solid/Liquid Interface

Most models that describe the interaction of an insoluble particle with an advancing solid-liquid interface are based on the assumption of steady state. However, as demonstrated by experimental work, the process does not reach steady state until the particle is pushed for a while by the interface. In this work, a dynamic mathematical model was developed. The dynamic model demonstrates that this interaction is essentially non-steady state and that steady state eventually occurs only when solidification is conducted at sub-critical velocities. The model was tested for three systems: aluminum-zirconia particles, succinonitrilepolystyrene particles, and biphenyl-glass particles. The calculated values for critical velocity of the pushing/engulfment transition were in same range with the experimental ones.

Catalina, A. V.↗

Dynamic Process Intensification via Data-Driven Dynamic Optimization: Concept and Application to Ternary Distillation

Process intensification is a design philosophy aimed at making chemical processes safer and more efficient. Its implementation often results in significant modifications to the design and structure of the process, with several conventional unit operations occurring in the same physical device. Traditionally, process intensification has focused on steady-state operation. In our previous works, we introduced dynamic process intensification (DPI) as a new intensification paradigm based on operational changes for conventional or intensified units. DPI is predicated on switching operation between two auxiliary steady states selected via a steady-state optimization calculation that ensures that the system generates, on average and over time, the same products as in nominal steady-state operation, but with favorable economics. This paper extends the DPI concept and introduces a novel dynamic optimization-based DPI strategy (Do-DPI) that involves imposing a true cyclic operation rather than switching between two discrete states. We discuss its implementation using surrogate dynamic models learned via system identification. Here, an extensive case study concerning a ternary distillation column separating a canonical hydrocarbon mixture shows that Do-DPI can reduce energy use by more than 4% relative to steady-state operation, with no significant deviations in product quality and production rate.

42 ENGINEERING↗

Device-Level Multidimensional Thermal Dynamics With Implications for Current and Future Wide Bandgap Electronics

Researchers have been extensively studying wide-bandgap (WBG) semiconductor materials such as gallium nitride (GaN) with an aim to accomplish an improvement in size, weight, and power of power electronics beyond current devices based on silicon (Si). However, the increased operating power densities and reduced areal footprints of WBG device technologies result in significant levels of self-heating that can ultimately restrict device operation through performance degradation, reliability issues, and failure. Typically, self-heating in WBG devices is studied using a single measurement technique while operating the device under steady-state direct current measurement conditions. However, for switching applications, this steady-state thermal characterization may lose significance since the high power dissipation occurs during fast transient switching events. Therefore, it can be useful to probe the WBG devices under transient measurement conditions in order to better understand the thermal dynamics of these systems in practical applications. In this work, the transient thermal dynamics of an AlGaN/GaN high electron mobility transistor (HEMT) were studied using thermoreflectance thermal imaging and Raman thermometry. Also, the proper use of iterative pulsed measurement schemes such as thermoreflectance thermal imaging to determine the steady-state operating temperature of devices is discussed. These studies are followed with subsequent transient thermal characterization to accurately probe the self-heating from steady-state down to submicrosecond pulse conditions using both thermoreflectance thermal imaging and Raman thermometry with temporal resolutions down to 15 ns.

42 ENGINEERING↗