Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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

Reactive Capture and Conversion of Carbon Dioxide to Methanol with ZnZrO 2 and Alkali-Promoted Mg 3 AlO x Mixed Oxide Catalytic Sorbents

Reactive capture and conversion (RCC) explores the use of a single-unit process to capture CO 2 and produce a product, in this case, methanol (MeOH). In this study, different configurations of a catalytic sorbent (CS) composed of ZnZrO 2 catalyst and Mg 3 AlO x sorbent with and without alkali modification are evaluated for CO 2 adsorption, steady-state catalysis with cofed CO 2 and H 2 , and transient RCC performance. A catalyst composed of a physical mixture of Mg 3 AlO x with ZnZrO 2 resulted in a slight increase in CO 2 uptake, with a low impact on the catalytic activity and RCC of the materials compared to ZnZrO 2 alone. In contrast, Na impregnation significantly increased the level of CO 2 uptake from 0.28 mmol/g (ZnZrO 2 alone) to 0.6 and 1.1 mmol/g for the CS with Na on the catalyst or Mg 3 AlO x , respectively. However, Na impregnation reduced the CO 2 conversion rate and MeOH selectivity during steady-state cofeed experiments at 300 °C and 6 bar. In contrast to steady-state catalysis conditions, RCC, which is a cyclic capture and conversion process, creates dynamic CO 2 and H 2 surface coverages, favoring CH 4 in the early stages of the conversion step and then CO and MeOH as the catalyst CO 2 coverage reduces. The highest MeOH productivity during RCC was achieved with CS that balanced the CO 2 uptake with only moderate catalyst rate reductions caused by Na addition. The optimal material, ZnZrO 2 +10%Na/Mg 3 AlO x , achieved a CO 2 uptake of 0.8 mmol/g and a MeOH productivity of 0.5 mmol/g with 100% selectivity at 260 °C and 6 bar during RCC. This marks the highest RCC MeOH productivity reported to date, although the process needs further optimization and even with optimization, may remain impractical. The results further demonstrate that optimization of catalytic sorbents under steady-state flow conditions does not easily correlate to transient capture and conversion cycles for methanol synthesis from CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
NASA NTRS1981

Advantages of multigrid methods for certifying the accuracy of PDE modeling

Numerical techniques for assessing and certifying the accuracy of the modeling of partial differential equations (PDE) to the user's specifications are analyzed. Examples of the certification process with conventional techniques are summarized for the three dimensional steady state full potential and the two dimensional steady Navier-Stokes equations using fixed grid methods (FG). The advantages of the Full Approximation Storage (FAS) scheme of the multigrid technique of A. Brandt compared with the conventional certification process of modeling PDE are illustrated in one dimension with the transformed potential equation. Inferences are drawn for how MG will improve the certification process of the numerical modeling of two and three dimensional PDE systems. Elements of the error assessment process that are common to FG and MG are analyzed.

Forester, C. K.↗
NASA NTRS2002

Magnetic Fluctuations in the Jovian Magnetosphere

The engine that drives the jovian magnetosphere is the mass added to the Io ion torus, accelerated to corotational velocities by field-aligned currents that couple the Io torus to the jovian ionosphere. The mass of the torus builds up to an amount that the magnetic forces cannot contain and the plasma, first slowly and then more rapidly, drifts outward. Numerous authors have treated this problem based first on the observations of the Pioneer 10 and 11 flybys; then on Voyager 1 and 2, and Ulysses; and finally most recently the Galileo orbiter. The initial observations revealed the now familiar magnetodisk, in which the field above and below the magnetic equator became quite radial in orientation and much less dipolar. The Galileo observations show this transformation to occur on average at 24 R(sub J) and to often be quite abrupt. These observations are consistent with outward transport of magnetized plasma that moves ever faster radially until about 50 R(sub J) on the nightside where the field lines stretch to the breaking point, reconnection occurs, and plasma and field islands are transported down the tail ultimately removing the mass from the magnetosphere that Io had deposited deep in the inner torus. The reconnection process creates empty flux tubes connected to Jupiter that are buoyant and thought to float inward and replace the flux carried out with the torus plasma. As described above, the jovian magnetosphere could very well be in a state of steady laminar circulation, but indeed it is not. The process is very unsteady and the wave levels can be very intense. The existence of these waves in turn can lead to processes that compete with the radial circulation pattern in removing plasma from the system. These waves can scatter particles so that they precipitate into the ionosphere. This process should be important in the Io torus where the atmospheric loss cone is relatively large and becomes less important as the loss cone decreases in size with radial distance. However, the Io torus is relatively quiet compared to the region outside the torus and it is not obvious without studying this scattering carefully whether the loss in the torus or out of the torus is greater and whether it can act rapidly enough to compete with the radial transport of ions to the tail in the life cycle of the mass added at Io. Closer to Io the ion cyclotron waves are most intense and possibly are associated with the losses in the Io flux tube. The waves are also diagnostic of both the Io atmospheric composition and the size and strength of the massloading process.

Russell, Christopher T.↗
NASA NTRS1973

Initial results from the NASA Lewis Bumpy Torus experiment

Initial results were obtained from low power operation of the NASA Lewis Bumpy Torus experiment, in which a steady-state ion heating method based on the modified Penning discharge is applied in a bumpy torus confinement geometry. The magnet facility consists of 12 superconducting coils, each 19 cm i.d. and capable of 3.0 T, equally spaced in a toroidal array 1.52 m in major diameter. A 18 cm i.d. anode ring is located at each of the 12 midplanes and is maintained at high positive potentials by a dc power supply. Initial observations indicate electron temperatures from 10 to 150 eV, and ion kinetic temperatures from 200 eV to 1200 eV. Two modes of operation were observed, which depend on background pressure, and have different radial density profiles. Steady state neutron production was observed. The ion heating process in the bumpy torus appears to parallel closely the mechanism observed when the modified Penning discharge was operated in a simple magnetic mirror field.

Roth, J. R.↗
NASA NTRS1989

Nonlinear combustion instability model in two- to three-dimensions

A model is developed on using quasi-steady state equations for the injection, atomization, and vaporization processes that consider the instantaneous local pressures, temperatures, density, and gas velocity vector associated with a 3-D wave. The coupled sets of equations are solved to determine frequency and growth rate of the oscillations as a function of wave amplitude. Calculated results are used to explain the observed instabilities in a 40 Klb LOX/Methane combustor.

Priem, Richard J.↗
NASA NTRS1993

Low power arcjet system spacecraft impacts

Application of electrothermal arcjets on communications satellites requires assessment of integration concerns identified by the user community. Perceived risks include plume contamination of spacecraft materials, induced arcing or electrostatic discharges between differentially charged spacecraft surfaces, and conducted and radiated electromagnetic interference (EMI) for both steady state and transient conditions. A Space Act agreement between Martin Marietta Astro Space, the Rocket Research Company, and NASA's Lewis Research Center was established to experimentally examine these issues. Spacecraft materials were exposed to an arcjet plume for 40 hours, representing 40 weeks of actual spacecraft life, and contamination was characterized by changes in surface properties. With the exception of the change in emittance of one sample, all measurable changes in surface properties resulted in acceptable end of life characteristics. Charged spacecraft samples were benignly and consistently reduced to ground potential during exposure to the powered arcjet plume, suggesting that the arcjet could act as a charge control device on spacecraft. Steady state EMI signatures obtained using two different power processing units were similar to emissions measured in a previous test. Emissions measured in UHF, S, C, Ku and Ka bands obtained a null result which verified previous work in the UHF, S, and C bands. Characteristics of conducted and radiated transient emissions appear within standard spacecraft susceptibility criteria.

Pencil, Eric J.↗
NASA NTRS1986

Structural dynamic measurement practices for turbomachinery at the NASA Lewis Research Center

Methods developed for measuring blade and rotor-shaft system response include optical systems, transient instruments, and special digital data processing equipment. Optical methods offer some distinct benefits for blade vibration measurement. Transient and steady state measurements of the response of rotor-shaft systems strongly affect analytical methods development. Digital computing systems allow processing of large volumes of high speed data from rotating blade sets. Also, digital systems develop useful vibration response signatures from randomly excited systems. Research facilities include the spin rig facility and the transient rotor response lab.

Kiraly, L. J.↗
DOE OSTI2021

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

Friction stir welding (FSW) has found increased applications in automotive and aerospace industries due to its advantages of solid-state bonding, no fusion and melting, and versatility in various working conditions and material combinations. However, the relationship among processing parameters, material properties, and bonding extent and fidelity remains largely empirical, primarily because of the lack of the mechanistic understanding of the tool-workpiece frictional behavior that affects our subsequent understanding of microstructural evolution and interface bonding formation. While the tool-workpiece stick-slip condition is believed to dictate the resulting torque and heat generation rate during the welding process, it remains rare and elusive to conduct a quantitative experimental measurement of such interfacial field. On the other hand, numerical simulations based on Computational Fluid Dynamics (CFD) rely on ad hoc assumptions of interfacial pressure and shear-stress conditions, but predictions can only be validated via the medium- and far-range temperature field which is known to be insensitive to the interfacial frictional behavior. This work first presents a comparison among two CFD-based simulation methodologies and the Coupled Eulerian Lagrangian (CEL) model in finite element method, the last of which uses the Coulomb friction so that the stick-slip is naturally developed. Based on the Hill-Bower similarity relationship in the contact analysis, an analytical model is developed here to prove why a constant stick-slip fraction will be developed in the steady state, to correlate the stick-slip fraction to processing parameters such as the tool spin rate, and further to derive dimensionless functions for torque and heat-generation-rate predictions. Pros and cons of various numerical approaches in predicting stick-slip are discussed, and our analytical model has been found to agree well with our numerical simulation and literature experimental results. These analyses provide the critical strain-rate and temperature fields that are needed for the bonding analysis in our future work.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations

Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Parameterization of planetary wave breaking in the middle atmosphere

A parameterization of planetary wave breaking in the middle atmosphere has been developed and tested in a numerical model which includes governing equations for a single wave and the zonal-mean state. The parameterization is based on the assumption that wave breaking represents a steady-state equilibrium between the flux of wave activity and its dissipation by nonlinear processes, and that the latter can be represented as linear damping of the primary wave. With this and the additional assumption that the effect of breaking is to prevent further amplitude growth, the required dissipation rate is readily obtained from the steady-state equation for wave activity; diffusivity coefficients then follow from the dissipation rate. The assumptions made in the derivation are equivalent to those commonly used in parameterizations for gravity wave breaking, but the formulation in terms of wave activity helps highlight the central role of the wave group velocity in determining the dissipation rate. Comparison of model results with nonlinear calculations of wave breaking and with diagnostic determinations of stratospheric diffusion coefficients reveals remarkably good agreement, and suggests that the parameterization could be useful for simulating inexpensively, but realistically, the effects of planetary wave transport.

Garcia, Rolando R.↗
DOE OSTI2023

Cross-flow separation characteristics and piloting of graphene oxide nanofiltration membrane sheets and tubes for kraft black liquor concentration

Dewatering of weak black liquor (WBL) in the kraft cycle by evaporation is highly energy intensive. Membranes are an attractive alternative for energy-efficient dewatering, but existing commercial polymeric or ceramic membranes are either degraded in BL or have high capital costs. Our recent works have demonstrated the engineering of graphene oxide (GO) nanofiltration membranes, their stability and promising performance in BL conditions, and preliminary scale-up into sheets and tubes. Here, we describe in detail the separation characteristics of GO membrane sheets and tubes under real BL conditions and crossflow operation. Recycle-mode piloting of a GO tubular membrane showed average “production flux” of 16 L/m2/h (LMH) and high rejections of lignin (98.3%), total solids (66%), and total organic carbon (83%), with no signs of irreversible fouling identified. A corresponding GO sheet membrane produced an average flux of ~25 LMH and maintained high lignin rejection of ~97% during a slipstream pilot at a kraft mill site using WBL with ~16 wt% total solids (TS). Finally, we piloted a Dow/DuPont XUS1808 polyamide composite reverse osmosis (RO) membrane for last-mile processing of the GO nanofiltration membrane permeate. The RO membrane showed a steady state flux of 19 LMH at 65 bar and produced ~0.02 wt% TS water product, which is highly suitable for reuse in pulp washing operations in the kraft process. The results have strong positive implications for the industrial application of GO membranes in BL concentration and other related applications.

Materials Science↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
NASA NTRS1995

Adaptive Optimization of Aircraft Engine Performance Using Neural Networks

Preliminary results are presented on the development of an adaptive neural network based control algorithm to enhance aircraft engine performance. This work builds upon a previous National Aeronautics and Space Administration (NASA) effort known as Performance Seeking Control (PSC). PSC is an adaptive control algorithm which contains a model of the aircraft's propulsion system which is updated on-line to match the operation of the aircraft's actual propulsion system. Information from the on-line model is used to adapt the control system during flight to allow optimal operation of the aircraft's propulsion system (inlet, engine, and nozzle) to improve aircraft engine performance without compromising reliability or operability. Performance Seeking Control has been shown to yield reductions in fuel flow, increases in thrust, and reductions in engine fan turbine inlet temperature. The neural network based adaptive control, like PSC, will contain a model of the propulsion system which will be used to calculate optimal control commands on-line. Hopes are that it will be able to provide some additional benefits above and beyond those of PSC. The PSC algorithm is computationally intensive, it is valid only at near steady-state flight conditions, and it has no way to adapt or learn on-line. These issues are being addressed in the development of the optimal neural controller. Specialized neural network processing hardware is being developed to run the software, the algorithm will be valid at steady-state and transient conditions, and will take advantage of the on-line learning capability of neural networks. Future plans include testing the neural network software and hardware prototype against an aircraft engine simulation. In this paper, the proposed neural network software and hardware is described and preliminary neural network training results are presented.

Simon, Donald L.↗
DOE OSTI2024

Effects of boundary hydraulics, dissolved oxygen, and dissolved organic carbon on growth and death dynamics of aerobic microbes in riverbed dune-induced hyporheic zones

Surface and groundwater interact in the hyporheic zone beneath and adjacent to rivers in the presence of a diverse microbial community. Heterotrophic bacteria mediate a range of environmentally important reactions, yet few studies have quantified bacterial growth and death dynamics in the hyporheic zone, and none have systematically analyzed their response to variations in hydraulic or chemical conditions. Here we used MODFLOW and SEAM3D to simulate hydraulics; dissolved oxygen (DO) and dissolved organic carbon (DOC) transport; and aerobic microbial metabolism, growth, and death in hyporheic zones induced by riverbed dunes. We ran simulations both with and without growth/death processes, and varied hydraulic parameters and DO/DOC boundary concentrations. Microbial biomass reached steady state (t = 3 days) in every simulation, at which time there was greater biomass and DOC biodegradation rates in the hyporheic flowcell (300% and 85% higher for the base case, respectively) when accounting for microbial growth dynamics. This occurred as microbial biomass tailored its spatial distribution to the availability of DO and DOC, demonstrating the importance of simulating growth/death processes. Biomass generally increased with hyporheic flow cell area as upwelling groundwater decreased. When varying surface water DO and DOC source concentrations relative to the base case, the greatest effect on biomass occurred when increasing DOC and decreasing DO. We determined minimum DO and DOC steady-state concentrations required for microbial growth, but the minimums were not absolute or related by stoichiometry. Increasing DOC created a smaller area of microbes with higher concentrations relative to the base case. Increasing DO slightly increased the area occupied by microbes while keeping the total biomass nearly constant. Overall, microbial growth and death dynamics depend on DO and DOC availability in the hyporheic zone, which is dependent on DOC/DO boundary concentrations and hyporheic flow paths, and in turn the hydraulic interaction between surface water and groundwater.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Predictive modeling of a subcritical pulverized-coal power plant for optimization: Parameter estimation, validation, and application

As renewable power generation deployment increases, fossil fuel plants are increasingly required to operate more flexibly. Many coal-fired power plants were originally designed to operate at base load and do not operate optimally at partial load. Predictive first-principles plant-wide models can be employed to identify opportunities for flexibility improvements and diagnose low-load operating issues. This paper describes the application of the Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES) to model and optimize flexible power plant operations. The key benefits of using IDAES are that it provides an open-source, fully equation-oriented modeling framework for efficient modular model construction, reuse, and customization, together with a mathematical optimization framework leveraging powerful, state-of-the-art solvers. The process systems engineering workflow from predictive process simulation to parameter estimation, model validation, and plant optimization is applicable to a variety of existing and next-generation energy systems as well as other chemical and environmental processes. Here, to demonstrate this capability, a physics-based, steady-state model was developed to improve full- and part-load performance of the Escalante Generating Station, a 245 MWe (net) subcritical pulverized coal-fired power plant owned and operated by Tri-State Generation and Transmission Association. Specifically, sixty-nine model parameters were simultaneously estimated from several months of operating data enabling prediction of flow rates, temperatures, pressures, and steam quality throughout the plant. The validated model was leveraged by Escalante to reduce the minimum operating load from 90 MW to 50 MW by diagnosing a low-load water-hammer issue, enabling coal usage and emissions reductions during periods of low power demand. Additionally, opportunities for heat rate reduction (i.e., efficiency improvement) through a steeper sliding-pressure approach to load-following and optimization of other boiler operating variables were also identified and quantified. For example, a potential efficiency improvement of 0.7 percentage points was observed at half-load operation.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1992

Subsonic flight test evaluation of a propulsion system parameter estimation process for the F100 engine

An adaptive-performance-seeking control system which optimizes the quasi-steady-state performance of the F-15 propulsion system is discussed. This paper presents flight- and ground-test evaluations of the propulsion system parameter-estimation process used by the performance seeking control system. The estimator consists of a compact propulsion system model and an extended Kalman filter. The extended Kalman filter estimates five engine component deviation parameters from measured inputs. The compact model uses measurements and Kalman-filter estimates as inputs to predict unmeasured propulsion parameters such as net propulsive force and fan stall margin. The ability to track trends and estimate absolute values of propulsion system parameters was demonstrated. For example, thrust stand results show a good correlation especially in trends between the performance seeking control estimated and measured thrust.

Orme, John S.↗