SEARCH · Engineering Papers
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.
In this work, we explore time-resolved dynamics of non-equilibrium magnons in Rb2MnF4, a model two-dimensional square lattice Heisenberg antiferromagnet using a novel laser-neutron pump-probe technique. This technique uses a nanosecond pulsed laser to periodically excite magnons to out-of-equilibrium states and then probe the resulting dynamics using INS.
Non-local thermodynamics equilibrium (LTE) effects in the photosphere; recent research on the chromosphere of the M and C stars; and elementary shock-waves and pulsation theories and their applications to Mira long-period variables are discussed.
In response to the need for ground testing facilities for super orbital re-entry research, a small scale facility has been set up at the University of Queensland to demonstrate the superorbital expansion tube concept. This unique device is a free piston driven, triple diaphragm, impulse shock facility which uses the enthalpy multiplication mechanism of the unsteady expansion process and the addition of a secondary shock driver to further heat the driver gas. The pilot facility has been operated to produce quasi-steady test flows in air with shock velocities in excess of 13 km/s and with a usable test flow duration of the order of 15 micro sec. an experimental condition produced in the facility with total enthalpy of 108 MJ/kg and a total pressure of 335 MPa is reported. A simple analytical flow model which accounts for non-ideal rupture of the light tertiary diaphragm and the resulting entropy increase in the test gas is discussed. It is shown that equilibrium calculations more accurately model the unsteady expansion process than calculations assuming frozen chemistry. This is because the high enthalpy flows produced in the facility can only be achieved if the chemical energy stored in the test flow during shock heating of the test gas is partially returned to the flow during the process of unsteady expansion. Measurements of heat transfer rates to a flat plate demonstrate the usability of test flow for aerothermodynamic testing and comparison of these rates with empirical calculations confirms the usable accuracy of the flow model.
We present abundance measurements from 26 emission-line-selected galaxies at z approx. 0.6-0.7. By reaching stellar masses as low as 10(exp 8) M stellar mass, these observations provide the first measurement of the intermediate-redshift mass-metallicity (MZ) relation below 10(exp 9)M stellar mass. For the portion of our sample above M is greater than 10(exp 9)M (8/26 galaxies), we find good agreement with previous measurements of the intermediate-redshift MZ relation. Compared to the local relation, we measure an evolution that corresponds to a 0.12 dex decrease in oxygen abundances at intermediate redshifts. This result confirms the trend that metallicity evolution becomes more significant toward lower stellar masses, in keeping with a downsizing scenario where low-mass galaxies evolve onto the local MZ relation at later cosmic times. We show that these galaxies follow the local fundamental metallicity relation, where objects with higher specific (mass-normalized) star formation rates (SFRs) have lower metallicities. Furthermore, we show that the galaxies in our sample lie on an extrapolation of the SFR-M∗ relation (the star-forming main sequence). Leveraging the MZ relation and star-forming main sequence (and combining our data with higher-mass measurements from the literature), we test models that assume an equilibrium between mass inflow, outflow, and star formation.We find that outflows are required to describe the data. By comparing different outflow prescriptions, we show that momentum, driven winds can describe the MZ relation; however, this model underpredicts the amount of star formation in low-mass galaxies. This disagreement may indicate that preventive feedback from gas heating has been overestimated, or it may signify a more fundamental deviation from the equilibrium assumption.
Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.
A new procedure closely linking dissociation and exchange reactions in air to the vibrational levels of the diatomic molecules has been implemented in both one- and two-dimensional versions of Direct Simulation Monte Carlo (DSMC) programs. The previous modeling of chemical reactions with DSMC was based on the continuum reaction rates for the various possible reactions. The new method is more closely related to the actual physics of dissociation and is more appropriate to the particle nature of DSMC. Two cases are presented: the relaxation to equilibrium of undissociated air initially at 10,000 K, and the axisymmetric calculation of shuttle forebody heating during reentry at 92.35 km and 7500 m/s. Although reaction rates are not used in determining the dissociations or exchange reactions, the new method produces rates which agree astonishingly well with the published rates derived from experiment. The results for gas properties and surface properties also agree well with the results produced by earlier DSMC models, equilibrium air calculations, and experiment.
The use of o-phthalaldehyde (OPA) in combination with a thiol reagent is a common method for detecting primary amines in amino acids, peptides, and proteins. Despite its widespread use, the exact reaction mechanism has been debated since the 1980s. Here, we measure the kinetics of the reaction between OPA, alanine, and a dithiol (1,4-dithiolthreitol, DTT) as a function of pH and reagent concentration. Using these new measurements and accompanying kinetic models, we find evidence that the pH dependence of the kinetics arises from both the protonation states of alanine and DTT, the hydration state of OPA, and the unproductive equilibrium with DTT, all of which are pH-dependent. Importantly, these results support the mechanism originally proposed by Sternson [Rational design and evaluation of improved o-phthalaldehyde-like fluorogenic reagents. Anal. Biochem. 1985, 144, 233–246] and Wong [Reaction of o-phthalaldehyde with alanine and thiols: kinetics and mechanism. J. Am. Chem. Soc. 1985, 107, 6421–6422], in which the primary amine first reacts with OPA, followed by a reaction with the thiol to form the fluorescent isoindole product.
The objective of this work was to assess the impact of flow field - material response cou-pling on reproducing temperatures within the Mars 2020 forebody thermal protection system.Both equilibrium and nonequilibrium coupled material response approaches were considered.Turbulent transition estimates were included in the flow field modeling based on transitiontimes inferred from flight measurements. A coupled equilibrium ablation model was foundto overpredict the in-depth temperatures near the surface by over 100 K in the stagnationregion and by about 200 K near the leeside shoulder, which was subject to turbulent heatingaugmentation. For nonequilibrium ablation, two finite rate surface chemistry models wereused to model the chemical process on the vehicle surface. The first was a model consistingof N2catalysis and CO2recombination model, and the second was a model consisting of onlyN2and O2catalysis. Near the stagnation point, the N2and CO2model was found to be withinabout 10 K from the peak in-depth temperature near the surface; however, overpredictionof the in-depth temperatures was observed early in the trajectory, similar to the equilibriumablation model. The N2and O2catalysis model was found to significantly underpredict nearsurface in-depth temperatures.
Cell performance optimization is important for improving the overall system efficiency of a redox flow battery. To gain better insights into key controlling factors of system efficiency, this work first proposed a theoretical model for a unit flow battery cell by extending a two-dimensional analytic model to a full battery cell. Such a model is then used for cell performance optimization after validating it with experimental and numerical modeling data. With the model results, the activation, equilibrium, and pump energy losses are identified as the dominant sources of battery energy losses. Further, a guideline for reducing these sources is also proposed. Following the guideline, the mass transport coefficient is shown as a key control factor of the equilibrium energy loss and Coulombic efficiency (CE). Approaches are then proposed to improve CE and the overall system efficiency. The mass transport is also revealed as the mechanism of pump rate optimization where an optimal pump rate significantly reduces the equilibrium energy loss. With both low equilibrium and pump energy losses, an optimal electrode porosity or specific area design can further improve a battery's system efficiency based on an optimal porosity predicted by the present model. The model also demonstrates distinct behaviors and overestimation in the system efficiency when reduced to a zero-dimensional model with neglected mass transport resistance. With the new model, the guideline, and new insights, this work provides a reliable and efficient tool for the evaluation and optimization of redox flow battery design in practical applications.
A climate model resolving the seasonal cycle and the two horizontal dimensions was developed and applied to several problems of current interest. Models of this type are useful when for various reasons a general circulation model experiment is not warranted or not feasible. For example, in cases where the signal to natural variability is small it may be advantageous to first consider such a statistical dynamical model because extremely long runs may be necessary in the application. In this case the simpler statistical dynamical model serves as a pilot study device. The model developed is a thermodynamic model whose solution yields the equilibrium seasonal cycle for the surface temperature field over the globe. The model is essentially a statement of the conservation of heat energy for individual columns of the earth atmosphere system. Various terms such as the infrared radiation flux to space are parameterized with earth radiation budget data from satellites such as Nimbus 6. The primary agent modulating the seasonal cycle amplitude is the heat capacity per unit area which is a strong function of surface type -- ocean surface can store 60 times more heat per unit time than land. By adjusting its few empirical parameters the model can be brought into remarkable agreement with the observed seasonal cycle. The model is then very useful for looking at the dependence of the seasonal cycle of the temperature on such externally defined variables as the Earth's orbital elements (eccentricity, tilt, precession of equinoxes) or the configuration of land-sea geography which can be changed by continental drift.
We investigate the mapping of magnetic flux tubes between the high-latitude ionosphere and the equatorial magnetosphere, using the Toffoletto and Hill (1989) theoretical model and the Tsyganenko (1987, 1989) empirical models. In all models, strong distortion of flux tubes occurs near the magnetopause and in the distant tail. These effects are illustrated for various interplanetary magnetic field directions in the Toffoletto and Hill (1989) model and for various Kp values in the Tsyganenko (1987, 1989) model. We calculate the Birkeland current that is actually present in each model (from Ampere's law), as well as the Birkeland current that would be implied if the model were in magnetostatic equilibrium with isotropic plasma pressure (from the Vasyliunas equation). Comparison of the actual and implied Birkeland currents indicates that no physically significant Birkeland current is included in any of the models tested. Each model contains a 'noise' component of Birkeland current that is, however, much smaller in magnitude than the 'real' Birkeland currents that remain to be included.
Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.
Non-equilibrium plasma discharges in spark gaps have been an increasingly studied method for alleviating cycle to cycle variation in lean and dilute combustion environments. However, ignition models that account for streamer propagation, cathode fall, and transmission line amplification over nanosecond time scales have so far not been developed. The present study develops such a model, with emphasis on the energy delivered from circuit to cylinder. Here, key pieces of the relevant physics and chemistry are summarized, simplified, and systematically coupled to one another. The set of parameters is limited to a handful of key observables and modeled using Modelica. Results show non-trivial behavior in the energy delivery characteristics of such discharges with important implications for ignition.
This paper describes theoretical models and experimental performance of a novel Zigzag Flow Reactor (ZFR) for weekly thermochemical energy storage. The ZFR reduces redox-active metal oxide (MO x ) particles at high temperature (up to ~1100 °C) under inert gas sweep. A physical model demonstrates the approach to process equilibrium by minimizing the associated exergy destruction in a finite number of reaction steps, establishing the thermodynamic requirements for a practical reactor. The model results show several cost-relevant parameter tradeoffs, and the tradeoff analysis implies a cost-optimized set of boundary conditions. Numerical models and prototypes show that the ZFR enables significant gas phase homogenization while simultaneously enabling a customizable MO x residence time in the reactor, both key requirements for approaching an equilibrium process. A scaling model demonstrates the simplicity and affordability of sizing the ZFR to grid-scale levels, with fabrication costs at least five times lower than previously proposed scalable reactor concepts. As a result, a laboratory ZFR prototype achieved an energy storage density of ~90 Wh/kg with CaAl 0.2 Mn 0.8 O 3-δ as the MO x , at temperatures of ~850 °C in >10 h of total runtime.
Two basically different models for the filament equilibrium by Kippenhahn and Schluter (1957) and Kuperus and Raadu (1974) have appeared in the literature. A further analyses by van Tend and Kuperus (1978) added the force due to the horizontal component of the background field to the Kuperus and Raadu model. In order to obtain a better model which actually describes these phenomena, the evolution of the filament has to be considered in detail. A first attempt was recently presented by Kaastra. Kaastra did not formulate the precise energy balance equations for the problem, as is done in the present work. In the present model not only the force balance, but also the energy balance of the filament is taken into account. Thus a fully closed system of equations is obtained, that describes the evolution of the filament, first in force equilibrium during the current build-up phase, then in the non-equilibrium phase before the eruption, and the eruption itself. A neutral point appears above the photospheric surface in the non-equilibrium phase, but long before the eruption. It was found that although the filament itself may be in non-equilibrium, the evolution may still be slow up to the height where the eruption takes place. The eruption of the filament itself causes a large induced electric field at the neutral point which leads to the observed flare phenomena.
Whether ocean iron fertilization in high nutrient low chlorophyll (HNLC) regions is a practical geoengineering strategy to combat climate warming has been debated because the fate of absorbed carbon by the ocean remains unclear. Here, we use an optimized oceanic biogeochemical inverse model to explore the equilibrium responses to increased productivity in major HNLC regions. We find that an overall downward shift of nutrients and carbon from the surface and intermediate to the deep waters with maximal productivity in the HNLC Southern Ocean (SO), North Pacific (NP) and eastern equatorial Pacific (EP). The efficiency of soft tissue pump is increased by 37.8%, 2.7% and 5.4% in the SO, EP and NP simulation while atmospheric CO 2 concentration is drawn down by 45.1, 3.4 and 6.7 ppm, respectively. Meanwhile, global mean oxygen is decreased by 21.1%, 1.3% and 2.4% in the SO, EP and NP simulation, accompanied by an expansion of oxygen minimum zones by 400%, 26% and 63%, respectively. Our study indicates that iron fertilization has the potential to be an effective geoengineering approach to mitigate climate warming but also has a strong impact on the development of ocean hypoxia.
While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.
Conditions of crater saturation and equilibrium have been simulated with a Monte Carlo computer program which utilizes a generating function of the form N = b x D to the -a. The following conclusions have been deduced from the results of the model: (1) saturation and equilibrium for large craters occur at crater densities substantially greater than those observed on the moon, Mars, and Mercury, an indication that the observed populations are essentially production populations, (2) generating functions with a greater than 2 lead to final crater populations with a = 2, (3) generating functions with a between 2 and 1.5 lead to crater populations with the parameter a approaching 1.0, and (4) changes in the slope of crater populations occur gradually and over a wide interval of crater densities and diameters as saturation is approached.