Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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,027 records · Page 57

Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools

The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.

pressurization↗

Laminar-to-turbulent flame transition and cycle-to-cycle variations in large eddy simulation of spark-ignition engines

This paper investigates the effect of laminar-to-turbulent flame transition modeling on the prediction of cycle-to-cycle variations (CCVs) in large eddy simulation (LES) of spark-ignition (SI) engines. A laminar-to-turbulent flame transition model that describes the non-equilibrium sub-filter flame speed evolution during an early stage of flame kernel growth is developed. In the present model, the flame transition is characterized by the flame kernel size at which the flame transition ends, defined here as the flame transition scale. The proposed model captures the effects that variations in a turbulent flow field have on the evolution of early-stage burning rates, through variations in the flame transition scale. The proposed flame transition model is combined with the front propagation formulation (FPF) method and a spark-ignition model to predict CCVs in a gasoline direct injection SI engine. It is found that multi-cycle LES with the proposed flame transition model reproduces experimentally-observed CCVs satisfactorily. When the transition model is not considered or when variations in the transition process are neglected, CCVs are significantly under-predicted for the case considered here. These results indicate the importance of modeling the laminar-to-turbulent flame transition and the effect of turbulence on the transition process, when predicting CCVs, under certain engine conditions. The LES results are also used to analyze sources for variations in the flame transition. It is found, for the present engine case, that the most important source is the cycle-to-cycle variation in the turbulence dissipation rate, which is used to measure the strength of turbulence in the proposed model, near a spark plug. The large-scale velocity field and the variations of the laminar flame speed due to the mixture composition and thermal stratification are also found to be important factors to contribute to the variations in the flame transition.

42 ENGINEERING↗

Non-LTE line-blanketed model atmospheres of hot stars. 2: Hot, metal-rich white dwarfs

We present several model atmospheres for a typical hot metal-rich DA white dwarf, T(sub eff) = 60,000 K, log g = 7.5. We consider pure hydrogen models, as well as models with various abundances of two typical 'trace' elements-carbon and iron. We calculte a number of Local Thermodynamic Equilibrium (LTE) and non-LTE models, taking into account the effect of numerous lines of these elements on the atmospheric structure. We demostrate that while the non-LTE effects are notvery significant for pure hydrogen models, except for describing correctly the central emission in H-alpha they are essential for predicting correctly the ionization balance of metals, such as carbon and iron. Previously reported discrepancies in LTE abundances determinations using C III and C IV lines are easily explained by non-LTE effects. We show that if the iron abundance is larger than 10(exp -5), the iron line opacity has to be considered not only for the spectrum synthesis, but also in the model construction itself. For such metal abundances, non-LTE metal line-blanketed models are needed for detailed abundance studies of hot, metal-rich white dwarfs. We also discuss the predicted Extreme Ultraviolet (EUV) spectrum and show that it is very sensitive to metal abundances, as well as to non-LTE effects.

Lanz, T.↗

Reynolds stress calculations of homogeneous turbulent shear flow with bounded energy states

Reynolds stress calculations of homogeneous turbulent shear flow are conducted with a second-order closure model modified to account for non-equilibrium vortex stretching in the dissipation rate transport equation, as recently proposed by Bernard and Speziale. As with the earlier reported k-epsilon model calculations incorporating this vortex stretching effect, a production-equals-dissipation equilibrium is obtained with bounded turbulent kinetic energy and dissipation. However, this equilibrium is not achieved until the dimensionless time greater than 60, an elapsed time that is at least twice as large as any of those considered in previous numerical and physical experiments on homogeneous shear flow. Direct quantitative comparisons between the model predictions and the results of experiments are quite favorable. In particular, it is shown that the inclusion of this non-equilibrium vortex stretching effect has the capability of explaining the significant range of production to dissipation ratios observed in experiments.

Speziale, Charles G.↗

Dynamic Oscillation and Motion of Oil‐in‐Water Emulsion Droplets

Abstract The interplay of interfacial tensions on droplets results in a range of self‐powered motions that mimic those of living systems and serve as a tunable model to understand their complex non‐equilibrium behavior. Spontaneous shape deformations and oscillations are crucial features observed in nature but difficult to incorporate in synthetic artificial systems. Here, we report sessile oil‐in‐water emulsions that exhibit rapid oscillating behavior. The oscillations depend on the nature and concentration of the surfactant, the chemical composition of the oil, and the wettability of the solid substrate. The rapid changes in the contact angle per oscillation are observed using side‐view optical microscopy. We propose that the changes in the interfacial tension of the oil droplets is due to the partitioning of the surfactant into the oil phase and the movement of self‐emulsified oil out of the parent droplets giving rise to the rhythmic variation in droplet contact‐line. The ability to control and understand droplet oscillation can help model similar oscillations in out‐of‐equilibrium systems in nature and reproduce biomimetic behavior in artificial systems for various applications, such as microfluidic lab‐on‐a‐chip and adaptive materials.

Krishna Mani, Sanjana↗

Dynamic Oscillation and Motion of Oil‐in‐Water Emulsion Droplets

The interplay of interfacial tensions on droplets results in a range of self‐powered motions that mimic those of living systems and serve as a tunable model to understand their complex non‐equilibrium behavior. Spontaneous shape deformations and oscillations are crucial features observed in nature but difficult to incorporate in synthetic artificial systems. Here, we report sessile oil‐in‐water emulsions that exhibit rapid oscillating behavior. The oscillations depend on the nature and concentration of the surfactant, the chemical composition of the oil, and the wettability of the solid substrate. The rapid changes in the contact angle per oscillation are observed using side‐view optical microscopy. We propose that the changes in the interfacial tension of the oil droplets is due to the partitioning of the surfactant into the oil phase and the movement of self‐emulsified oil out of the parent droplets giving rise to the rhythmic variation in droplet contact‐line. The ability to control and understand droplet oscillation can help model similar oscillations in out‐of‐equilibrium systems in nature and reproduce biomimetic behavior in artificial systems for various applications, such as microfluidic lab‐on‐a‐chip and adaptive materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formation of Minor Phases in a Nickel-Based Disk Superalloy

The minor phases of powder metallurgy disk superalloy LSHR were studied. Samples were consistently heat treated at three different temperatures for long times to approximate equilibrium. Additional heat treatments were also performed for shorter times, to then assess non-equilibrium conditions. Minor phases including MC carbides, M23C6 carbides, M3B2 borides, and sigma were identified. Their transformation temperatures, lattice parameters, compositions, average sizes and total area fractions were determined, and compared to estimates of an existing phase prediction software package. Parameters measured at equilibrium sometimes agreed reasonably well with software model estimates, with potential for further improvements. Results for shorter times representing non-equilibrium indicated significant potential for further extension of the software to such conditions, which are more commonly observed during heat treatments and service at high temperatures for disk applications.

Gabb, T. P.↗

A comparative study of the time-dependent standard 8-, 13- and 16-moment transport formulations of the polar wind

The ionosphere, composed of O(+), H(+), and electrons is modeled with four different transport formulations. The equations corresponding to the standard set, the 8-, 13-, and 16-moment approximations are presented, and the collision terms are expressed. Using a time-dependent technique, the ionosphere is studied between altitudes of 200 and 8600 km. The production of electrons and O(+) ions is described by a neutral atmosphere simplified photoionization scheme, while the energy deposition is supported by a downward electron heat flow of -0.005 erg/sq cm per s imposed at the topside boundary. When the models reach a steady state equilibrium, the electron solutions show differences due to the introduction of temperature anisotropies and heat flows between the components parallel and perpendicular to the magnetic field. As a corollary, the ions show structures depending on the level of approximation. A depletion of a factor of 10 is then applied to the ion densities above a certain altitude, and the development of the perturbation is followed for 1000 s for all the models.

Blelly, P. L.↗

Simple Approximation for the Ideal Reference State of Gases Adsorbed on Solid-State Surfaces

Reference states are useful as models for facilitating calculations of equilibrium constants, and they may also serve as standard states that are convenient for organizing and tabulating thermodynamic data; however, standard state conventions and appropriate reference states for adsorbed species have received less attention than those for pure substances and solutes. Here, we compare seven choices of reference states for calculations of equilibrium constants and transition state theory rate constants for flat surfaces, in particular (1) an ideal 2D harmonic oscillator, (2) an ideal rigid-molecule harmonic oscillator, (3) an ideal 2D harmonic oscillator with separable surface modes, (4) a 2D ideal gas, (5) an ideal 2D hindered translator, (6) an ideal 2D hindered translator with lowest-order barriers, and (7) a simple ideal 2D hindered translator proposed in this work. The advantage of models 5–7 is that they can treat both mobile and localized adsorbates in a consistent way, whereas models 1–3 are only appropriate for localized adsorbates, and model 4 is only appropriate for a freely translating adsorbate. Furthermore, models 6 and 7 reduce the computational cost without the user having to calculate barrier heights for diffusion. An advantage of the simple ideal 2D hindered translator is that it has a physical high-temperature limit. We also propose a reference state for nonflat surfaces. Here, the user is encouraged to choose a reference state based on the appropriateness of the model and the practicality of the calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Models of resistive wall tearing mode disruptions

Disruptions are a serious issue in tokamaks. In a disruption, the thermal energy is lost by means of an instability which could be a resistive wall tearing mode (RWTM). During precursors to a disruption, the plasma edge region cools, causing the current to contract. Model sequences of contracted current equilibria are given, and their stability is calculated. A linear stability study shows that there is a maximum value of edge q a ≈ 3 for RWTMs to occur. Further, this also implies a minimum rational surface radius normalized to plasma radius from RWTMs to be unstable. Nonlinear simulations are performed using a similar model sequence derived from an equilibrium reconstruction. There is a striking difference in the results, depending on whether the wall is ideal or resistive. With an ideal wall, the perturbations saturate at moderate amplitude, causing a minor disruption without a thermal quench. With a resistive wall, there is a major disruption with a thermal quench, if the edge q a ≤ 3. There is a sharp transition in nonlinear behavior at q a = 3. This is consistent with the linear model and with experiments. If disruptions are caused by RWTMs, then devices with highly conducting walls, such as the International Tokamak Experimental Reactor will experience much milder, tolerable, disruptions than presently predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent phenomena in correlated materials

Understanding time-dependent processes and light-matter interaction in strongly correlated materials, and the interplay between electronic, orbital, vibrational, and spin degrees of freedom, is a cornerstone of condensed matter. These mechanisms can be proven by measuring the response of the systems to time-dependent perturbations. The corresponding time scales are dictated by the way light couples to the different excitations, and how these excitations exchange energy and momentum. Our research advances our understanding of these processes, and the interpretation of different equilibrium and time-resolved spectroscopies. Our project encompasses two main themes: (i) developing and refining computational techniques to study non-equilibrium spectroscopies including non-perturbative effects and (ii) applications to non-equilibrium phenomena. We have developed a new computational approach that works directly in the time domain: by including all the degrees of freedom involved in the scattering process (e.g. electrons, photons, neutrons), we solve the time dependent problem: a faithful numerical simulation of the experiment. By measuring the energy and momentum of the outgoing particles, we can extract information about the energy and momentum absorbed by the system. Prior to our work, people attempting to model and calculate non-equilibrium spectral functions relied on a description of the scattering cross section based on a formulation in the frequency domain, a treatment that is extremely cumbersome and complex. Our technique works in and out of equilibrium and can reproduce spectra by several spectroscopic techniques, such as time-resolved photoemission, neutron scattering, Raman, X-ray spectroscopies (RIXS, Auger, XAS, XMCD), and, by not relying on analytical approximations, yields results that reveal novel overlooked transient mechanisms. These tools provide sorely needed intuition for understanding the phenomenology of strongly correlated materials and will help experimentalists in identifying signatures of relevant excitations in pump-probe experiments, such as those conducted in DOE supported facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structure of confined laminar spray diffusion flames: Numerical investigation

The structure of confined laminar spray diffusion flames is investigated numerically by solving the gas-phase conservation equations for mass species, continuity, momentum, and energy and the liquid-phase equations for droplet position, velocity, size, and temperature. A one-step global reaction scheme along with six equilibrium reactions are employed to model the flame chemistry. Monodisperse as well as polydisperse sprays are considered. The numerical results demonstrate that liquid spray flames substantially differ from gaseous flames in their structure, i.e., temperature, concentration, and velocity fields, shape, and dimensions under the same conditions. Spray flames are predicted to be taller and narrower than their counterpart gaseous ones and their shapes are almost cylindrical. This is in agreement with experimental observations. The numerical computations also show that the use of the equilibrium reactions with the one-step reaction scheme decreases the flame temperature compared to the one-step reaction scheme without the equilibrium reactions and more importantly increases the surface area of the flame zone due to a phenomenon termed 'equilibrium broadening.' The spray flames also possess a finite thickness with minimal overlap of the fuel and oxygen species. A case for which a fuel-mixture consisting of 20 to 80 percent gas-liquid by mass is introduced into the combustor is also investigated and compared with predictions using only gaseous or liquid fuel.

Mawid, M. A.↗

Nonhomologous contraction and equilibria of self-gravitating, magnetic interstellar clouds embedded in an intercloud medium: Star formation. II - Results

The paper studies the equilibrium of self-gravitating isothermal models of interstellar clouds with a frozen-in magnetic field linked smoothly to the field of a hot tenuous intercloud medium. Equilibrium states are determined which can be reached by clouds contracting nonhomologously from a spherical uniform initial state. Three free parameters characterize the problem: a dimensionless initial radius related to the Jeans length of the cloud, the initial ratio of the magnetic and gas pressures in the cloud, and the initial ratio of the intercloud and cloud pressures. The dependence of the solutions on each of these parameters is investigated. It is found that: (1) the frozen-in field causes the cloud to become oblate with its major axis normal to the field lines; (2) the flattening increases as the magnetic-field strength, gravitational forces, or intercloud pressure increases; (3) increasing intercloud pressure eventually leads to gravitational collapse; and (4) the cloud can reach equilibrium only if its radius does not exceed some critical value. The observed inefficiency of the star-formation process within massive clouds is examined and explained in terms of magnetic phenomena in a collapsing cloud.

Mouschovias, T. C.↗

Deposition control of model glasses with surface-mediated orientational order

We introduce a minimal model of solid-forming anisotropic molecules that displays, in thermal equilibrium, surface orientational order without bulk orientational order. The model reproduces the nonequilibrium behavior of recent experiments in which a bulk nonequilibrium structure grown by deposition contains regions of orientational order characteristic of the surface equilibrium. This order is deposited, in general, in a nonuniform way because of the emergence of a growth-poisoning mechanism that causes equilibrated surfaces to grow slower than non-equilibrated surfaces. We use evolutionary methods to design oscillatory protocols able to grow nonequilibrium structures with uniform order, demonstrating the potential of protocol design for the fabrication of this class of materials.

Whitelam, Stephen (ORCID:0000000200866803)↗

Reversible Interfacial Charge Transfer and Delayed Emission in InP/ZnSe/ZnS Quantum Dots with Hexadecanethiol

The results in this paper show that holes are rapidly and reversibly transferred from red-emitting InP/ZnSe/ZnS quantum dots (QDs) to adsorbed hexadecanethiol (HDT) forming an equilibrium between the thiols and the QD valence band. Photoexcitation results in populations of holes in the valence band and in slightly higher-energy shell-localized traps. Trap to valence band hole tunneling results in a photoluminescence risetime having time constants varying from 300 ps to 2 ns. The presence of adsorbed HDT eliminates the slower risetime component, indicating that hole transfer from the shell-localized traps that are closest to the particle surface efficiently competes with tunneling to the QD core. This shows that the interfacial charge transfer equilibrium is established in less than 2 ns. The population of the shell-localized traps corresponds to a reservoir of hole states that eventually tunnel to the core-localized valence band, resulting in delayed emission. The amount of delayed emission increases rapidly with ZnSe shell thickness and is slightly blue-shifted from the prompt photoluminescence. We propose an energetic model in which the HDT/valence band equilibrium is affected by the extent of valence band quantum confinement and an electric field produced by core–shell interfacial dipoles. Furthermore, this model explains the core size, shell thickness, and photoluminescence (PL) wavelength dependence of this equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mapping Interfacial Solution Forces that Drive Fluorescent Nanoparticle Incorporation into Crystals

Interfacial solution structure governs processes ranging from catalytic and electrochemical reactions to particle aggregation and crystallization. Here, we design bio-inspired nanocomposites by mapping interfacial forces and controlling nanoparticle incorporation during crystal growth. Our analytical model contains measured kinetic barriers for surface approaches and equilibrium binding constants with the crystal surface. We validate this model using fluorescent silica nanoparticles and calcite. Our results show that surface chemistry dictates incorporation pathways: methoxy groups are least favorable, hydroxyls reduce kinetic barriers by interacting favorably with hydration layers, carboxylates bind strongly but must overcome a kinetic barrier, and amines outperform all kinetically and thermodynamically. This framework resolves the interplay between interfacial solution structures, particle forces and dynamics, and growth kinetics, enabling predicative control of nanocomposite formation and multiplexed mixed-particle systems.

Zhang, Mingyi [Pacific Northwest National Laborato↗

Climate warming due to increasing atmospheric CO2 - Simulations with a multilayer coupled atmosphere-ocean seasonal energy balance model

The transient response of the climate to increasing CO2 is studied using a modified version of the multilayer energy balance model of Peng et al. (1982). The main characteristics of the model are described. Latitudinal and seasonal distributions of planetary albedo, latitude-time distributions of zonal mean temperatures, and latitudinal distributions of evaporation, water vapor transport, and snow cover generated from the model and derived from actual observations are analyzed and compared. It is observed that in response to an atmospheric doubling of CO2, the model reaches within 1/e of the equilibrium response of global mean surface temperature in 9-35 years for the probable range of vertical heat diffusivity in the ocean. For CO2 increases projected by the National Research Council (1983), the model's transient response in annually and globally averaged surface temperatures is 60-75 percent of the corresponding equilibrium response, and the disequilibrium increases with increasing heat diffusivity of the ocean.

Li, Peng↗