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,009 records · Page 56

Estimates of total quantity of meteorites in the East Antarctic ice cap

Based on a steady-state model of the East Antarctic ice cap, and current estimates of meteorite influx, a model is developed which predicts that the steady-state number of meteorites being carried in or on the ice is at least 760,000. This large meteorite population does not require unusual influx conditions, since the cold, dry climate preserves virtually all meteorites that fall with the exception of the fragile, porous carbonaceous chondrites. Application of the model to the Greenland ice cap yields a steady-state population of about 61,000 meteorites.

Olsen, E. J.↗

Characterizing the Non-equilibrium Dynamics of Field-Driven Correlated Quantum Systems

Recent experimental advances in ultrafast phenomena have triggered renewed interest in the dynamics of correlated quantum systems away from equilibrium. We review non-equilibrium dynamical mean-field theory studies of both the transient and the steady states of a DC field-driven correlated quantum system. In particular, we focus on the non-equilibrium behavior and how it relates to the fluctuation-dissipation theorem. The fluctuation-dissipation theorem emerges as an indicator for how the system thermalizes and for how it reaches a steady state. When the system thermalizes in an infinite temperature steady state it can pass through a succession of quasi-thermal states that approximately obey the fluctuation-dissipation theorem. We also discuss the Wigner distribution and what its evolution tells us about the non-equilibrium many-body problem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Design, analysis, and tests of a shuttle-type heat-pipe-cooled leading edge

A one-half scale heat-pipe-cooled leading edge model was designed and fabricated to verify feasibility and performance of a full-scale Phase B shuttle orbiter design. Model performance was investigated experimentally by radiant heating and aerothermal tests and analytically by using a simple technique which predicts heat pipe start-up from the frozen state and also predicts transient and steady-state thermal behavior. Analytical results agree well with experimental results for start-up and steady-state heat pipe operation. The results indicate that variations in angle of attack and roll orientation had a negligible effect on heat pipe performance. The heat pipes effectively isothermalized the leading edge, and reduced peak temperatures to levels compatible with the use temperatures of superalloys. Results of these tests demonstrated the durability of the heat-pipe-cooled leading edge in withstanding earth-entry thermal and mechanical loads and indicate that the use of heat pipes in leading edge structures is a reliable concept for fully reusable hypersonic cruise and space transportation systems.

Camarda, C. J.↗

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of a turbojet engine considered as a quasi-static system

A determination of the dynamic characteristics of a typical turbojet engine with a centrifugal compressor, a sonic-flow turbine-nozzle diaphragm, and fixed area exhaust nozzle is presented. A generalized equation for transient behavior of the engine was developed; this equation was then verified by calculations using compressor and turbine performance charts extrapolated from equilibrium operating data and by experimental data obtained from an engine operated under transients in fuel flow. The results indicate that a linear differential equation for engine acceleration as a function of fuel flow and engine speed for operation near a steady-state operating condition can be written. The coefficients of this equation can be obtained either from actual transient data or, with a fair degree of accuracy, from the steady-state performance maps of the compressor and turbine, and can be corrected for altitude in the same manner that steady-state performance data are corrected.

Otto, Edward W↗

Numerical Modeling of Solidification in Space With MEPHISTO-4

A pre-flight analysis of the directional solidification of BiSn with MEPHISTO-4 is presented. Simplified Bridgman growth under microgravity conditions is simulated using a two dimensional finite element model. This numerical model is a single domain, pseudo-steady state model, and includes the effects of both thermal and solutal convection. The results show that for all orientations of the applied steady state gravity vector, of magnitude 1 micro-g, the directional solidification process remains diffusion controlled. The maximum convective velocity was found to be 4.424 x 10(exp -5) cm/s for the horizontal Bridgman growth configuration. This value is an order of magnitude lower than the growth velocity. The maximum and minimum values or solute concentration in the liquid at the crystal-melt interface were 13.867 at.% and 13.722 at.%, respectively. This gives a radial segregation value of xi = 1.046% at the interface. A secondary objective of this work was to compare the results obtained to those that consider thermal convection only (no solutal convection). It was found that the convective flow patterns in simulations which included solutal convection were significantly different from those which ignored solutal convection. The level of radial segregation predicted by the current simulations is an order of magnitude lower than that found in simulations which ignore solutal convection. The final aim was to investigate the effect of g-jitter on the crystal growth process. A simulation was performed to calculate the system response to a 1 second, 100 micro-g gravity impulse acting normal to the direction of growth. This pulse is consistent with that induced by Orbiter thruster firings. The results obtained indicate that such a gravity pulse causes an increase in the level of radial solute segregation at the interface from the steady state values. The maximum value of solute concentration in the liquid was found to be 13.888 at.%, the minimum value calculated was 13.706 at.%, yielding a radial segregation value of xi = 1.31% at the interface. These values occurred 126 seconds after the pulse terminated. Thus it is anticipated that the process will remain diffusion controlled even when subjected to such g-jitter.

Simpson, James E.↗

The photosynthetic response of C 3 and C 4 bioenergy grass species to fluctuating light

Abstract Bioenergy grass species are a renewable energy source, but their productivity has not been fully realized. Improving photosynthetic efficiency has been proposed as a mechanism to increase the productivity of bioenergy grass species. Fluctuating light, experienced by all field grown crops, is known to reduce photosynthetic efficiency. This experiment aimed to evaluate the photosynthetic performance of both C 3 and C 4 bioenergy grass species under steady state and fluctuating light conditions by examining leaf gas exchange. The fluctuating light regime used here decreased carbon assimilation across all species when compared to expected steady state values. Overall, C 4 species assimilated more carbon than C 3 species during the fluctuating light regime, with both photosynthetic types assimilating about 16% less carbon than expected based on steady state measurements. Little diversity was observed in response to fluctuating light among C 3 species, and photorespiration partially contributed to the rapid decreases in net photosynthetic rates during high to low light transitions. In C 4 species, differences among the four NADP‐ME species were apparent. Diversity observed among C 4 species in this experiment provides evidence that photosynthetic efficiency in response to fluctuating light may be targeted to increase C 4 bioenergy grass productivity.

09 BIOMASS FUELS↗

DC-Link Current Minimization Control for Current Source Converter-Based Solid-State Transformer

This article proposes a fast predictive control method and a small DC-link inductor to minimize the DC-link current in current-source converter (CSC)-based solid-state transformer. The DC-link current minimization can significantly reduce power loss and improve efficiency. The challenge of this problem is on improving both steady-state and dynamic performance. PI control methods and large DC-link inductors are conventionally used in the CSC but have limited dynamic performance. A model predictive control (MPC) method is proposed to achieve switching-cycle-level settling time, and the DC-link inductor is sized for 40% ripple to enable fast current change. Importantly, this article also proposes to minimize the DC-link current by varying the current even within a line cycle under single-phase load to improve the steady-state performance, in contrast with the reduction to a constant value in the literature. The proposed MPC features a constant switching frequency without weighting factors. The MPC does not have a high computational burden and is implemented in a regular digital controller for a prototype of soft-switching solid-state transformer (S4T) with reduced conduction loss. The effectiveness of the proposed method has been experimentally verified on the SiC S4T prototype during steady-state and dynamics under different multiport power flow conditions up to 2 kV peak. Here, the DC-link current in the experiments is close to the minimum current with a short zero-vector duration, which further verifies the performance of the proposed method.

14 SOLAR ENERGY↗

Combustion of a premixed system in stagnation flow. I - Theoretical.

The theoretical steady-state and transient characteristics of a premixed combustible system in laminar axisymmetric stagnation flow are studied by way of numerical solution of the governing boundary-layer equations. Computations, which are based on the oxidation reaction of carbon monoxide in humid air, invoke a minimum of simplifying assumptions. Steady-state results include prediction of (1) profiles of the state variables, (2) extinction of blow-off conditions, (3) ignition conditions, and (4) the heat flux to the stagnation point. Computer simulations of the unsteady state provide an indication of the stability of a given steady state to perturbations. Comparisons between some of the predictions of the mathematical model and the results of laboratory experiments will be reported in Part II.

Smith, H. W.↗

Orbit-modulated transport and sources in time-dependent plasmas

The orbit-modulation framework for steady-state ion-orbit boundary fluxes is reformulated to apply to time-dependent plasmas. The new formulation can quantify the extent to which orbit fluxes are affected by the time dependence of the distribution function, and may be applied both to time-dependent physical scenarios and to computationally expensive simulations that are unable to reach a true steady state. As the system tends to a time-independent state, the orbit-loss fluxes approach their steady-state values, which are determined entirely by upstream transport (turbulent and collisional) and upstream sources and sinks. This result holds for the fluxes of any moments that are expressible in terms of the invariants that determine the orbits, including parallel toroidal angular momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Studies of pressure-velocity coupling schemes for analysis of incompressible and compressible flows

Two pressure-velocity coupling schemes, both of which solve the fully implicit discretized equations governing the flow of fluids were examined, and the capability of performing large Reynolds number, low Mach number compressible flow calculations were assessed. The semi-implicit iterative SIMPLE algorithm is extended to handle transient compressible flow calculations. This extension takes into account a strong coupling between the pressure and temperature through a correction procedure, based on the equations of state. Results obtained from the extended SIMPLE algorithm are then compared to similar results obtained from the non-iterative PISO algorithm. Both time-dependent and steady state calculations were performed using an axisymmetric 2:1 pipe expansion geometry and laminar flow conditions corresponding to Reynolds number of 1000 and Mach number of 2.0. For calculations simulating a time-dependent compression/expansion process, both schemes exhibit transient features in excellent agreement with each other, and moreover, the PISO method shows a significant computational time reduction of 60 percent over the SIMPLE scheme, regardless of the time step size or grid size employed. The effects of numerical diffusion are shown to be significant in these calculations. For steady state compressible flows, however, the SIMPLE algorithm displays increasing computational efficiency over the PISO method as the time step sizes employed to reach steady state conditions are decreased.

Burgreen, Gregory Wayne↗

Discovering SPP Thermomechanical Pathways

Final report for SPPSi project in Thrust One Friction extrusion is a thermomechanical process that combines conventional extrusion with the action of a rotating die. The plastic deformation of the material being sheared and extruded is the primary source of process heat and it produces strain distributions unlike those resulting from conventional extrusion. This paper proposes an improved strain analysis that evaluates three main strain components in a series of rate-controlled friction extrusions in which the steady state was achieved. Cylindrical AA1100 extrusion billets with two embedded markers were extruded to wire with a 10:1 diametral reduction. The shape change of the embedded markers was determined via serial, transverse sectioning and quantitative metallography of the extruded wires. Three mutually orthogonal strain components (longitudinal, circumferential, and radial) were calculated at different positions along each extruded wire from the marker's shape change. The development of strain from the initial transient to the steady state is discussed. The variation of the steady-state strain with different process parameters is correlated with the die advance per revolution.

42 ENGINEERING↗

Wind Loads on Flat Plate Photovoltaic Array Fields

The aerodynamic forces resulting from winds acting on flat plate photovoltaic arrays were investigated. Local pressure distributions and total aerodynamic forces on the arrays are shown. Design loads are presented to cover the conditions of array angles relative to the ground from 20 deg to 60 deg, variable array spacings, a ground clearance gap up to 1.2 m (4 ft) and array slant heights of 2.4 m (8 ft) and 4.8 m (16 ft). Several means of alleviating the wind loads on the arrays are detailed. The expected reduction of the steady state wind velocity with the use of fences as a load alleviation device are indicated to be in excess of a factor of three for some conditions. This yields steady state wind load reductions as much as a factor of ten compared to the load incurred if no fence is used to protect the arrays. This steady state wind load reduction is offset by the increase in turbulence due to the fence but still an overall load reduction of 2.5 can be realized. Other load alleviation devices suggested are the installation of air gaps in the arrays, blocking the flow under the arrays and rounding the edges of the array. A wind tunnel test plan to supplement the theoretical study and to evaluate the load alleviation devices is outlined.

Miller, R.↗

Characterization of the space shuttle reaction control system engine

A computer program was developed and written in FORTRAN 5 which predicts the transient and steady state performance and heat transfer characteristics of a pulsing GO2/GH2 rocket engine. This program predicts the dynamic flow and ignition characteristics which, when combined in a quasi-steady state manner with the combustion and mixing analysis program, will provide the thrust and specific impulse of the engine as a function of time. The program also predicts the transient and steady state heat transfer characteristics of the engine using various cooling concepts. The computer program, test case, and documentation are presented. The program is applicable to any system capable of utilizing the FORTRAN 4 or FORTRAN 5 language.

Wilson, M. S.↗

Flux Balancing Control of Ungapped Nanocrystalline Core-Based Transformer in Dual Active Bridge Converters

Transformer flux dc bias is a critical issue, impacting the reliable operation of dual active bridge (DAB) converters especially when ungapped high permeability nanocrystalline core is used. Steady state current dc bias can easily saturate ungapped nanocrystalline transformers, and it is even more dangerous in transient conditions. A dc bias model is proposed to analyze steady state dc bias in different load conditions. The magnetizing current detection is necessary for closed-loop control of dc bias; conductors' position is shown to impact sensor noise, and it is analyzed in detail. To deal with both steady state and transient dc bias, a unified flux balancing control (UFBC) is proposed introducing a predictive bias suppression (PBS) method with closed-loop flux balancing control (CFBC). With the proposed PBS, both primary/secondary current balance and flux balance can be achieved within one switching cycle using UFBC. Power characteristics and interaction between power control and flux balancing control of DAB converters are analyzed, and the CFBC needs to work in a low bandwidth due to sensor bandwidth limitation and interaction between power control and flux balancing control. Lastly, the UFBC is verified on a 300-kW cascaded DAB converter prototype.

42 ENGINEERING↗

Origin, transport, and losses of energetic He(+) and He(2+) ions in the magnetosphere of the Earth - AMPTE/CCE observations

Data from the ion charge-energy-mass spectrometer CHEM flown on AMPTE/CCE spacecraft are used to investigate the origin, transport, and losses of energetic He(+) and He(2+) ions in the earth's magnetosphere. The L profiles of the average ion phase space density f were determined as a function of the magnetic momentum. It is shown that the L profiles have an inner part, where f increases with L for both He(+) adn He(2+) and where steady-state conditions are fulfilled. The outer boundary L(lim) of this region is located at a distance that depends on the ion species and the geomagnetic activity level. Steady-state conditions continue outside L(lim) for He(+) ions, while the He(2+) ion distribution outside L(lim) is strongly influenced by ion convection causing a lack of steady-state conditions. It is concluded that solar wind is the origin of the He(2+), while a mixed origin is suggested for the He(+) ions, in which the major contribution is from the solar wind via charge exchange production from the He(2+) ions.

Kremser, G.↗

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

Comparing Constant and Transient Membrane Transport Parameters for Use in Wave Desalination Models

Directly pressurizing seawater for desalination with reverse osmosis membranes via wave motion is a promising and sustainable method for producing freshwater in coastal regions. However, such a system could result in significant pressure fluctuations and a departure from conventional steady-state desalination operations. This study sought to assess if membrane transport parameters (apparent water and salt permeability) should be modeled as transient or constant in solution–diffusion-based modeling efforts of dynamically operated desalination systems, such as those coupled to wave power. Two approaches were used to model membrane transport parameters: one considered each parameter to be a function of the net driving pressure of the system, and the other assumed they were constant across all conditions. A pilot-scale system was used to conduct steady-state and controlled ramping experiments. Data from steady-state experiments were used to calculate transient and constant transport parameters. Parameter combinations were used in a simulation model to predict water flux and effective permeate salinity, and simulation outcomes were compared against experimental ramping results. The transient relationships for both water and salt permeability produced the most accurate results for water flux and comparable results for effective permeate salinity. Development of such relationships would be unique to a specific system but could be valuable in modeling wave-driven desalination systems across the wide range of operating conditions they experience.

13 HYDRO ENERGY↗