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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Reshaping of Truncated Pd Nanocubes: Energetic and Kinetic Analysis Integrating Transmission Electron Microscopy with Atomistic-Level and Coarse-Grained Modeling

Stability against reshaping of metallic fcc nanocrystals synthesized with tailored far-from-equilibrium shapes is key to maintaining optimal properties for applications such as catalysis. Yet Arrhenius analysis of experimental reshaping kinetics, and appropriate theory and simulation, is lacking. As such, we use TEM to monitor the reshaping of Pd nanocubes of ~25 nm side length between 410 °C (over ~4.5 h) and 440 °C (over ~0.25 h), extracting a high effective energy barrier of E eff ≈ 4.6 eV. We also provide an analytic determination of the energy variation along the optimal pathway for reshaping that involves transfer of atoms across the nanocube surface from edges or corners to form new layers on side {100} facets. The effective barrier from this analysis is shown to increase strongly with the degree of truncation of edges and corners in the synthesized nanocube. Theory matches experiment for the appropriate degree of truncation. In addition, we perform simulations of a stochastic atomistic-level model incorporating a realistic description of diffusive hopping for undercoordinated surface atoms, thereby providing a visualization of the initial reshaping process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Acoustic oscillations of Jupiter

We calculate low-degree l acoustic (p-mode) oscillations of the giant planet Jupiter. We use two Jupiter models, both of which contain a rocky core, an ice mantle, and a hydrogen-rich envelope. One of the models incorporates the 'plasma phase transition' of hydrogen, a first-order phase transition that yields a density discontinuity in the hydrogen-rich envelope. Our calculations of the acoustic oscillation spectrum of Jupiter have taken account of the effects of the rapid rotation of the planet: the Coriolis force, centrifugal force, and deformation of the equilibrium state. The p-mode frequency spectra of the Jupiter models are summarized in several echelle diagrams. We show that the frequency spectra of the low-degree p-modes are strongly affected by the existence of the rocky core and by avoided crossings among high radial-order p-modes with different l's.

Lee, Umin↗
DOE OSTI2024

Super Resolving Unrolled Neural Networks for Remote Sensing

In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.

47 OTHER INSTRUMENTATION↗
NASA NTRS1993

Radiative interactions in molecular gases under local and nonlocal thermodynamic equilibrium conditions

Basic formulations, analyses, and numerical procedures are presented to investigate radiative heat interactions in diatomic and polyatomic gases under local and nonlocal thermodynamic equilibrium conditions. Essential governing equations are presented for both gray and nongray gases. Information is provided on absorption models, relaxation times, and transfer equations. Radiative flux equations are developed which are applicable under local and nonlocal thermodynamic equilibrium conditions. The problem is solved for fully developed laminar incompressible flows between two parallel plates under the boundary condition of a uniform surface heat flux. For specific applications, three diatomic and three polyatomic gases are considered. The results are obtained numerically by employing the method of variation of parameters. The results are compared under local and nonlocal thermodynamic equilibrium conditions at different temperature and pressure conditions. Both gray and nongray studies are conducted extensively for all molecular gases considered. The particular gases selected for this investigation are CO, NO, OH, CO2, H2O, and CH4. The temperature and pressure range considered are 300-2000 K and 0.1-10 atmosphere, respectively. In general, results demonstrate that the gray gas approximation overestimates the effect of radiative interaction for all conditions. The conditions of NLTE, however, result in underestimation of radiative interactions. The method developed for this study can be extended to solve complex problems of radiative heat transfer involving nonequilibrium phenomena.

Tiwari, S. N.↗
NASA NTRS1977

Departures from hydrostatic equilibrium in the global distribution of thermospheric argon

This report describes the application of a three-dimensional two-component model to the distribution of argon in the upper atmosphere of the earth. The Cira (1972) model is used to specify the background gas density and the temperature as functions of altitude, latitude, local time, and day of the year. Emphasis is placed on examining the effects of departures from hydrostatic equilibrium on the argon density distribution. Such effects are found to be significant, producing a strong summer argon bulge. Comparisons of the model results with observational data are discussed.

Straus, J. M.↗
NASA NTRS1985

Multiple-scale turbulence closure modeling of confined recirculating flows

A multiple-scale turbulence closure scheme is developed for the numerical predictions of confined recirculating flows. This model is based on the multiple-time-scale concepts of Hanjalic et al. (1980) and takes into account the non-equilibrium spectra energy transfer mechanism. Problems concerning new formulation of energy transfer rate equations and subsequent model coefficient redefinition and energy spectrum partition are discussed. Comparisons are made with several experiments of internal recirculating flows for the purpose of model validation. Numerical results using the present model show significant improvement of predictive capability over that obtained with the single-scale k-epsilon model and show promising potential for complex turbulent flow predictions.

Chen, C. P.↗
NASA NTRS1990

Numerical simulation of a catastrophe model for coronal mass ejections

Van Tend and Kuperus (1978) proposed a simple catastrophe model for magnetically driving coronal mass ejections, prominence eruptions, and two-ribbon flares. Their model, which is based on simple circuit concepts, suggests that a stable configuration containing a current filament will lose equilibrium when the filament current exceeds a critical value. Here, a two-dimensional numerical simulation is used to test how the Van Tend-Kuperus model works in an ideal MHD fluid. The simulation exhibits the expected loss of mechanical equilibrium near the predicted critical value, but the current filament moves only a short distance upward before coming to rest at a new equilibrium. However, this new equilibrium contains a current sheet which is resistively unstable to magnetic reconnection, and if magnetic reconnection occurs rapidly, the filament can continue to move upward at Alfvenic speeds.

Forbes, T. G.↗
DOE OSTI2021

Ultrafast conformational dynamics of Rydberg-excited N -methyl piperidine

We have observed the ultrafast conformational dynamics of electronically excited N-methyl piperidine (NMP) using time-resolved Rydberg fingerprint spectroscopy. Optical excitation at various wavelengths ranging from 212 nm to 229 nm leads to the 3s or 3p Rydberg states and induces coherent oscillatory motions with periods of about 700 fs. Furthermore, these coherent motions survive the internal conversion from 3p to 3s but then dephase on a time scale of a few oscillations. Intramolecular vibrational energy redistribution on a picosecond time scale leads to an equilibrium between two conformeric structures that are separated in binding energy by 0.09 eV. Model calculations using the DFT-SIC method are in excellent agreement with the experiments and identify the conformers as the chair and twist structures of NMP. The analysis of the equilibrium parameters at long time delays as a function of excitation wavelength allows for the extraction of thermodynamic parameters for the conformeric transformation. We derive an enthalpy of the chair to twist reaction in the 3s excited state of 62 meV with an entropy of 19.70 J mol –1 K –1 . An activation energy of 276 meV is also obtained with a kinetic model.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2020

Passive Roll Stabilization of the Near Earth Asteroid Scout Solar Sail Mission

The Near Earth Asteroid (NEA) Scout is a small satellite, solar sail mission set to launch on Artemis I. Analysis of the NEA Scout solar radiation pressure model determined that, for certain solar incidence angles, there exists at least one locally stable equilibrium point about the “roll” axis, or the axis normal to the solar sail plane. This analysis is extended to three other geometrically similar solar sail models, and a Lyapunov stability analysis is conducted demonstrating that the stable equilibria are locally stable for the majority of roll angles in a wide range of solar incidence angles. Under certain assumptions, it is possible that only two axis control is required for non-spinning solar sails. If the roll axis is controlled, an active rate damping controller can be designed to asymptotically stabilize the system. Simulation results are presented that follow the stability analysis, where in the homogenous case, the system enters a limit cycle and, in the actively damped case, the system asymptotically converges to the equilibrium point.

Ivan Rodrigues Bertaska↗
DOE OSTI2023

Thermal & Electrochemical Power Plant Design and Cost Estimation

Public textbook for "Thermal & Electrochemical Power Plant Design and Cost Estimation: Version#1" This public textbook is an extension of class notes from Carnegie Mellon University courses: Energy System Modeling (24-722) and Fuel Cell Systems (24-262), taught by Dr. Nicholas Siefert between 2010-2021. Textbook includes some references to class notes from Dr. Shawn Litster, Department of Mechanical Engineering, Carnegie Mellon University. Textbook covers the equilibrium and nonequilibrium thermodynamics of power systems as well as an overview of system and economic modeling of these systems. There is in-depth coverage of (a) entropy generation, (b) exergy and (c) the redox state of molecules in equilibrium with the natural environment. This textbook is integrated with other materials (such as lecture slides, solved homeworks, and solved exams) that will be posted to the PowerShare: Energy Systems Modeling group on EDX. Publication Number: DOE/NETL-2023/3913

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

On the Dependence of Simulated Convection on Domain Size in CRMs

Abstract We present a heuristic model to explain the suppression of deep convection in convection‐resolving models (CRMs) with a small number of grid columns, such as those used in super‐parameterized or multi‐scale modeling framework (MMF) general circulation models (GCM) of the atmosphere. Domains with few grid columns require greater instability to sustain convection because they force a large convective fraction, driving strong compensating subsidence warming. Updraft dilution, which is stronger for reduced horizontal grid spacing, enhances this effect. Thus, suppression of deep convection in CRMs with few grid columns can be reduced by increasing grid spacing. Radiative‐convective equilibrium simulations using standalone CRM simulations with the System for Atmospheric Modeling (SAM) and using GCM‐coupled CRM simulations with the Energy Exascale Earth System Model (E3SM)‐MMF confirm the heuristic model results.

CRM↗
NASA NTRS1991

The chromospheres of late-type stars. I - Eridani as a test case of multiline modelling

A new model of the lower chromosphere of the dwarf K2 star Epsilon Eridani is derived by matching flux profiles of the Ca IR triplet lines 8498 and 8542 A H-alpha and H-beta lines and the Na D lines (all observed simultaneously at the AAT), and the Ca II K line. The coupled non-LTE equations of statistical equilibrium and radiative transfer are solved under the constraint of hydrostatic equilibrium using the Carlsson (1986) code. Within the framework of the model, the Na D lines are an important photospheric diagnostic, and the Ca IR triplet lines can be used to locate the temperature minimum. The computed H-alpha and H-beta depths are highly sensitive constraints on the transition zone gradients and base pressures allowing us to derive a pressure at the base of the transition zone of 0.9 dyn/cm.

Thatcher, John D.↗
DOE OSTI2022

Description of historical and future projection simulations by the global coupled E3SMv1.0 model as used in CMIP6

Abstract. This paper documents the experimental setup and general features of the coupled historical and future climate simulations with the first version of the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv1.0). The future projected climate characteristics of E3SMv1.0 at the highest emission scenario (SSP5-8.5) designed in the Scenario Model Intercomparison Project (ScenarioMIP) and the SSP5-8.5 greenhouse gas (GHG) only forcing experiment are analyzed with a focus on regional responses of atmosphere, ocean, sea ice, and land. Due to its high equilibrium climate sensitivity (ECS of 5.3 K), E3SMv1.0 is one of the Coupled Model Intercomparison Project phase 6 (CMIP6) models with the largest surface warming by the end of the 21st century under the high-emission SSP5-8.5 scenario. The global mean precipitation change is highly correlated with the global temperature change, while the spatial pattern of the change in runoff is consistent with the precipitation changes. The oceanic mixed layer generally shoals throughout the global ocean. The annual mean Atlantic meridional overturning circulation (AMOC) is overly weak with a slower change from ∼ 11 to ∼ 6 Sv (Sverdrup) relative to other CMIP6 models. The sea ice, especially in the Northern Hemisphere, decreases rapidly with large seasonal variability. We detect a significant polar amplification in E3SMv1.0 from the atmosphere, ocean, and sea ice. Comparing the SSP5-8.5 all-forcing experiment with the GHG-only experiment, we find that the unmasking of the aerosol effects due to the decline of the aerosol loading in the future projection period causes transient accelerated warming in the all-forcing experiment in the first half of the 21st century. While the oceanic climate response is mainly controlled by the GHG forcing, the land runoff response is impacted primarily by forcings other than GHG over certain regions, e.g., southern North America, southern Africa, central Africa, and eastern Asia. However, the importance of the GHG forcing on the land runoff changes grows in the future climate projection period compared to the historical period.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1980

Wake characteristics of buildings in disturbed boundary layers

Measurements relevant to the effect of buildings on the low level atmospheric boundary layer are presented. Field measurements of velocity and turbulence in the wake of a block building 3.2 m high and 26.8 m long are presented which show an apparent increase in momentum flow above the upwind value. Velocity-deficit and turbulence-excess decay characteristics of the disturbed or nonequilibrium layer are correlated with power law exponents and apparent roughness length at various distances downstream of the disturbance. Model wake profiles from the simulated building are compared at various stations for equilibrium and nonequilibrium upstream profiles. Empirical correlations relating building wake profiles to upstream nonequilibrium parameters are presented. The relationship of the data to the smooth-rough transition is discussed, and a flow model is presented.

Logan, E., Jr.↗
NASA NTRS1986

Further shock tunnel studies of scramjet phenomena

Scramjet phenomena were studied using the shock tunnel T3 at the Australian National University. Simple two dimensional models were used with a combination of wall and central injectors. Silane as an additive to hydrogen fuel was studied over a range of temperatures and pressures to evaluate its effect as an ignition aid. The film cooling effect of surface injected hydrogen was measured over a wide range of equivalence. Heat transfer measurements without injection were repeated to confirm previous indications of heating rates lower than simple flat plate predictions for laminar boundary layers in equilibrium flow. The previous results were reproduced and the discrepancies are discussed in terms of the model geometry and departures of the flow from equilibrium. In the thrust producing mode, attempts were made to increase specific impulse with wall injection. Some preliminary tests were also performed on shock induced ignition, to investigate the possibility in flight of injecting fuel upstream of the combustion chamber, where it could mix but not burn.

Morgan, R. G.↗
NASA NTRS1993

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 nonequilibrium vortex stretching in the dissipation rate transport equation as recently proposed by Bernard and Speziale (1992). 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 now not achieved until the dimensionless time St 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 nonequilibrium vortex stretching effect has the capability of explaining the significant range of production to dissipation ratios observed in experiments.

Speziale, Charles G.↗
NASA NTRS1992

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 a residual vortex stretching term is maintained 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 presented for a k-epsilon model modified to account for net vortex stretching.

Bernard, P. S.↗
DOE OSTI2021

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗