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,675 records · Page 93

Computational Modeling of Cephalad Fluid Shift for Application to Microgravity-Induced Visual Impairment

An improved understanding of spaceflight-induced ocular pathology, including the loss of visual acuity, globe flattening, optic disk edema and distension of the optic nerve and optic nerve sheath, is of keen interest to space medicine. Cephalad fluid shift causes a profoundly altered distribution of fluid within the compartments of the head and body, and may indirectly generate phenomena that are biomechanically relevant to visual function, such as choroidal engorgement, compromised drainage of blood and cerebrospinal fluid (CSF), and altered translaminar pressure gradient posterior to the eye. The experimental body of evidence with respect to the consequences of fluid shift has not yet been able to provide a definitive picture of the sequence of events. On earth, elevated intracranial pressure (ICP) is associated with idiopathic intracranial hypertension (IIH), which can produce ocular pathologies that look similar to those seen in some astronauts returning from long-duration flight. However, the clinically observable features of the Visual Impairment and Intracranial Pressure (VIIP) syndrome in space and IIH on earth are not entirely consistent. Moreover, there are at present no experimental measurements of ICP in microgravity. By its very nature, physiological measurements in spaceflight are sparse, and the space environment does not lend itself to well-controlled experiments. In the absence of such data, numerical modeling can play a role in the investigation of biomechanical causal pathways that are suspected of involvement in VIIP. In this work, we describe the conceptual framework for modeling the altered compartmental fluid distribution that represents an equilibrium fluid distribution resulting from the loss of hydrostatic pressure gradient.

Nelson, Emily S.↗

MIRO Observation of Comet C/2002 T7 (LINEAR) Water Line Spectrum

Comet C/2002 T7 (LINEAR) was observed with the Microwave Instrument for Rosetta Orbiter (MIRO) on April 30, 2004, between 5 hr and 16 hr UT. The comet was 0.63AU distance from the Sun and 0.68AU distance from the MIRO telescope at the time of the observations. The water line involving the two lowest rotational levels at 556.936 GHz is observed at 557.070 GHz due to a large Doppler frequency shift. The detected water line spectrum is interpreted using a non local thermal equilibrium (Non-LTE) molecular excitation and radiative transfer model. Several synthetic spectra are calculated with various coma profiles that are plausible for the comet at the time of observations. The coma profile is modeled with three characteristic parameters: outgassing rate, a constant expansion velocity, and a constant gas temperature. The model calculation result shows that for the distant line observation where contributions from a large coma space is averaged, the combination of the outgassing rate and the gas expansion velocity determines the line shape while the gas temperature has a negligible effect. The comparison between the calculated spectra and the MIRO measured spectrum suggests that the outgassing rate of the comet is about 2.0x1029 molecules/second and its gas expansion velocity about 1.2 km/s at the time of the observations.

spectral line↗

Predictive Simulation of Gas Adsorption in Fixed-Beds and Limitations due to the Ill-Posed Danckwerts Boundary Condition

The 1-D axially dispersed plug flow model is a mathematical model widely used for the simulation of adsorption processes. Lumped mass transfer coefficients such as the Glueckauf linear driving force (LDF) term and the axial dispersion coefficient are generally obtained by fitting simulation results to the experimental breakthrough test data. An approach is introduced where these parameters, along with the only free parameter in the energy balance equations, are individually fit to specific test data that isolates the appropriate physics. It is shown that with this approach this model provides excellent simulation results for the C02 on zeolite SA sorbent/sorbate system; however, for the H20 on zeolite SA system, non-physical deviations from constant pattern behavior occur when fitting dispersive experimental results with a large axial dispersion coefficient. A method has also been developed that determines a priori what values of the LDF and axial dispersion terms will result in non-physical simulation results for a specific sorbent/sorbate system when using the one-dimensional axially dispersed plug flow model. A relationship between the steepness of the adsorption equilibrium isotherm as indicated by the distribution factor, the magnitude of the axial dispersion and mass transfer coefficient, and the resulting non-physical behavior is derived. This relationship is intended to provide a guide for avoiding non-physical behavior by limiting the magnitude of the axial dispersion term on the basis of the mass transfer coefficient and distribution factor.

Knox, James Clinton↗

Phase coexistence implications of violating Newton’s third law

Newton’s third law, action = reaction, is a foundational statement of classical mechanics. However, in natural and living systems, this law appears to be routinely violated for constituents interacting in a nonequilibrium environment. Here, we use computer simulations to explore the macroscopic phase behavior implications of breaking microscopic interaction reciprocity for a simple model system. We consider a binary mixture of attractive particles and introduce a parameter that is a continuous measure of the degree to which interaction reciprocity is broken. In the reciprocal limit, the species are indistinguishable, and the system phase separates into domains with distinct densities and identical compositions. Furthermore, increasing nonreciprocity is found to drive the system to explore a rich assortment of phases, including phases with strong composition asymmetries and three-phase coexistence. Many of the states induced by these forces, including traveling crystals and liquids, have no equilibrium analogs. By mapping the complete phase diagram for this model system and characterizing these unique phases, our findings offer a concrete path forward toward understanding how nonreciprocity shapes the structures found in living systems and how this might be leveraged in the design of synthetic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A high‐order discontinuous Galerkin approach for physics‐based thermospheric modeling

Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.

Engineering↗

Modeling the decay of energy containing eddies: A source of solar wind heating

To understand the solar wind heating and acceleration mechanisms one needs to understand the decay of energy containing eddies. With this goal in mind, attempts have been made to extend the fluid dynamic phenomenology of large scale quasi-equilibrium to the case of magnetohydrodynamics. Matthaeus et al. have proposed a model for the inhomogeneous transport and decay of five mean variables, namely, two mean square Elsasser variables z(exp 2) (sub +/-) their correlation lengths, and the difference between the kinetic and magnetic energies. We test the validity of this model in the simplified case of homogeneous turbulence simulated in a periodic box. We propose a class of models and show that they may fit the simulation satisfactorily. Analytic solutions of this class of model reveal their inherent properties and demonstrate the difficulties associated with finite cross helicity. It is noted that adjustments are required to make the simplest models, which are based upon isotropic turbulence, scale properly with respect to the strength of the mean magnetic field. This can be interpreted as due to anisotropic turbulence, which can be modelled by simple parameterization in the phenomenology.

Hossain, M.↗

Bounded energy states in homogeneous turbulent shear flow: An alternative view

The equilibrium structure of homogeneous turbulent shear flow is investigated from a theoretical standpoint. Existing turbulence models, in apparent agreement with physical and numerical experiments, predict an unbounded exponential time growth of the turbulent kinetic energy and dissipation rate; only the anisotropy tensor and turbulent time scale reach a structural equilibrium. It is shown that if vortex stretching is accounted for in the dissipation rate transport equation, then there can exist equilibrium solutions, with bounded energy states, where the turbulence production is balanced by its dissipation. Illustrative calculations are present for a k-epsilon model modified to account for vortex stretching. The calculations indicate an initial exponential time growth of the turbulent kinetic energy and dissipation rate for elapsed times that are as large as those considered in any of the previously conducted physical or numerical experiments on homogeneous shear flow. However, vortex stretching eventually takes over and forces a production-equals-dissipation equilibrium with bounded energy states. The validity of this result is further supported by an independent theoretical argument. It is concluded that the generally accepted structural equilibrium for homogeneous shear flow with unbounded component energies is in need of re-examination.

Bernard, Peter S.↗

Liquid rocket combustor computer code development

The Advanced Rocket Injector/Combustor Code (ARICC) that has been developed to model the complete chemical/fluid/thermal processes occurring inside rocket combustion chambers are highlighted. The code, derived from the CONCHAS-SPRAY code originally developed at Los Alamos National Laboratory incorporates powerful features such as the ability to model complex injector combustion chamber geometries, Lagrangian tracking of droplets, full chemical equilibrium and kinetic reactions for multiple species, a fractional volume of fluid (VOF) description of liquid jet injection in addition to the gaseous phase fluid dynamics, and turbulent mass, energy, and momentum transport. Atomization and droplet dynamic models from earlier generation codes are transplated into the present code. Currently, ARICC is specialized for liquid oxygen/hydrogen propellants, although other fuel/oxidizer pairs can be easily substituted.

Liang, P. Y.↗

Analytic modeling of a spray diffusion flame

A detailed model for a spray diffusion flame is described. The model is based on the boundary layer form of the equations of motion, with droplet transport accounted for using a discretized droplet size distribution function. Interphase transport of mass and energy are accounted for, with a flame-sheet model used to describe the combustion process on a droplet scale. Near dynamic equilibrium is assumed for the description of droplet transport; droplets can diffuse relative to the gas phase. Gas-phase mixing is accounted for using a two-equation turbulence model; buoyancy effects are included, with a temperature fluctuation equation used to account for buoyancy effects on turbulence structure. Thermal radiation from gas-phase CO2 and H2O is included. Gas-phase chemical kinetics are modeled using a 20-reaction, 10-species version of the advanced quasi-global chemical kinetics formulation. Results are compared with data for a vaporizing Freon spray and a pentane spray flame. It is shown that the computational approach provides a reasonably valid picture of the overall development of a spray diffusion flame, and, furthermore, provides a useful tool for the parametric examination of the spray combustion process.

Harsha, P. T.↗

2D Quantum Transport Modeling in Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density- gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Quantum simulations are focused on MIT 25, 50 and 90 nm "well- tempered" MOSFETs and compared to classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are quantitatively consistent with I D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and sub-threshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗

Survey of the Relationship Between Theory and Experiment for Convective Heat Transfer From Rocket Combustion Gases

The problem of predicting the distribution of heat flux in a rocket engine is complicated by a flow field and thermodynamics so perturbed by as yet indescribable combustion effects that it continues to be necessary to utilize a simplified theoretical model against which to compare experimental results in order to gain insight. The most successful model applied to this problem is that of a compressible, turbulent, accelerating boundary layer flow with an arbitrary point of origin, and chemical reaction restricted to equilibrium recombination or diffusion-limited recombination on a catalytic wall. Analysis of this model was not possible without making several assumptions about the nature of the skin friction coefficient relationship and Reynold's analogy in an accelerating boundary layer flow. The paper summarizes the analysis of this model (together with detailed treatment of the special problems of variable properties, chemical reaction, and sources of transport properties), and describes several closed-form equations that have also been developed and applied to the problem. Described are results of several experiments with hot air flow made for the purpose of checking the validity of the assumptions in the boundary layer analysis. Most significant, however, are presentations of sample distributions of local heat flux measured in rocket engines operating with a wide range of conditions, propellants, injectors and sizes. When compared with predictions, these data show the significant range of deviation in heat flux that results from the combustion effects and which precludes a generally accurate prediction theory or correlation equation until the combustion process itself is correlated.

D R Bartz↗

Non-equilibrium of high-temperature nitrogen and oxygen

A stationary colloidal-radiative model for determination of production of electrons, excited atoms, and spectral line intensities in high-temperature atomic N and O is presented. The populations of several atomic levels and partition functions, and the intensities of several lines, are calculated in a wide range of conditions. Transport of radiation is included by coupling the rate equations for production of the electrons and excited atoms with the Holstein escape factors which are not constant, but depend on plasma conditions.

Kunc, J. A.↗

The solubility of olivine in basaltic liquids - An ionic model

A model is presented which enables the temperature at which olivine is in equilibrium with any alkali-depleted basaltic compound to be calculated to within + or - 30 C. It is noted that the error increases substantially when applied to terrestrial basalts which contain several weight percent alkalis. In addition the model predicts and quantifies the reduced activity of SiO4(4-) monomers due to increasing SiO2 concentrations in the melt. It is shown that the coordination of alumina in melts which precipitate olivine only appears to be dominantly octahedral, while titanium acts as a polmerizing agent by interconnecting previously isolated SiO4(4-) monomers. It is concluded that the model is sufficiently sensitive to show that there are small repulsive forces between Mg(2+) and calcium ions which are in association with normative diopside in the melt.

Herzberg, C. T.↗

A numerical experiment on the equilibrium and stability of a rotating galactic bar

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

Miller, R. H.↗

Electrochemical-Mechanical Coupling Strongly Affects the Performance of Nanopore, Thin-Film, and Solid-State Batteries

Here, the present study focuses on the electrochemical-mechanical (ECM) coupling effects of a thin-film, solid-state battery with only stiff, ceramic materials, in contrast to prior investigations that focus on individual active material particles or aggregations of particles. We model the impacts of ECM couplings including stress-transport and stress-equilibrium potential on the full-cell performance and potential mechanical failure modes of a thin-film battery conformally deposited in a nanopore scaffold, which is an experimentally achievable device. Model results indicate electrode volume changes due to lithium insertion or removal, along with mechanical boundary conditions, result in stress gradients that alter the lithium-ion flux, reduce lithium concentration gradients, and improve cell rate capability. However, the high stress levels in the cell can also lead to mechanics-related failure such as the separation of cell layers. For the parameter set in this work, stress-transport coupling has a much greater influence on rate capability than stress-potential coupling. Optimization of thin-film batteries to harness the benefits of ECM coupling effects requires leveraging geometric design and material selection. The current work underscores the need for further theoretical and experimental investigation into ECM coupling effects in thin-film batteries to enhance their understanding and design optimization.

25 ENERGY STORAGE↗

2D Quantum Mechanical Study of Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25, 50 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. Surprisingly, the self-consistent potential profile shows lower injection barrier in the channel in quantum case. These results are qualitatively consistent with ID Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗

Cyclic uniaxial and biaxial hardening of type 304 stainless steel modeled by the viscoplasticity theory based on overstress

The isotropic theory of viscoplasticity based on overstress does not use a yield surface or a loading and unloading criterion. The inelastic strain rate depends on overstress, the difference between the stress and the equilibrium stress, and is assumed to be rate dependent. Special attention is paid to the modeling of elastic regions. For the modeling of cyclic hardening, such as observed in annealed Type 304 stainless steel, and additional growth law for a scalar quantity which represents the rate independent asymptotic value of the equilibrium stress is added. It is made to increase with inelastic deformation using a new scalar measure which differentiates between nonproportional and proportional loading. The theory is applied to correlate uniaxial data under two step amplitude loading including the effect of further hardening at the high amplitude and proportional and nonproportional cyclic loadings. Results are compared with corresponding experiments.

Yao, David↗

Quantum Non-Equilibrium Dynamics Prediction of Electronic Transport Coefficients in Non-Linear Regimes

The goal of this work was to develop unique LLNL capabilities for the computational prediction of electrical, optical and thermal transport properties from quantum simulations. Using a recently developed theory of microscopic quantum transport, we were able to model compressed metallic hydrogen and it’s non-linear and out-of-equilibrium contributions in electrical and thermal transport. This theory was implemented in the open-source, real-space time-dependent density functional theory code Octopus. This model is general and widely applicable to high energy density systems that are important to LLNL and NNSA mission areas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗