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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Aqueous Alteration and Hydrogen Generation on Parent Bodies of Unequilibrated Ordinary Chondrites: Thermodynamic Modeling for the Semarkona Composition

Ordinary chondrites are the most abundant class of meteorites that could represent rocky parts of solar system bodies. However, even the most primitive unequilibrated ordinary chondrites (UOC) reveal signs of mild alteration that affected the matrix and peripheral zones of chondrules. Major chemical changes include oxidation of kamacite, alteration of glass, removal of alkalis, Al, and Si from chondrules, and formation of phases enriched in halogens, alkalis, and hydrogen. Secondary mineralogical changes include formation of magnetite, ferrous olivine, fayalite, pentlandite, awaruite, smectites, phosphates, carbonates, and carbides. Aqueous alteration is consistent with the oxygen isotope data for magnetite. The presence of secondary magnetite, Ni-rich metal alloys, and ferrous silicates in UOC implies that H2O was the oxidizing agent. However, oxidation by H2O means that H2 is produced in each oxidative pathway. In turn, production of H2, and its redistribution and possible escape should have affected total pressure, as well as the oxidation state of gas, aqueous and mineral phases in the parent body. Here we use equilibrium thermodynamic modeling to explore water-rock reactions in UOC. The chemical composition of gas, aqueous, and mineral phases is considered.

Zolotov, M. Y.↗
NASA NTRS1985

Three-dimensional viscous-shock-layer applications for the Space Shuttle orbiter

A nonorthogonal three-dimensional viscous-shock-layer code was modified and used to compute the windward surface flowfield over a Shuttle-like configuration. Code modifications were made in the methods used to evaluate shock boundary conditions and in the implementation of the equilibrium air models. Additional viscous terms were also retained in the governing equations. Effects of these modifications are demonstrated. Three-dimensional flowfield solutions obtained with the modified code showed good comparison of surface heating and pressure with actual flight data. An apparently three-dimensional effect, characterized by a drop in surface heating on the Shuttle windward centerline in the wing root area, was predicted by the present method. A preliminary review of Shuttle ground-based experimental and flight data confirmed this trend. The present three-dimensional solutions were also compared with the predictions of several approximate methods. Those comparisons showed the approximate methods to be valid for predicting heating along the Shuttle windward centerline uner both perfect gas and equilibrium air conditions.

Thompson, R. A.↗
NASA NTRS2005

Fast Whole-Engine Stirling Analysis

An experimentally validated approach is described for fast axisymmetric Stirling engine simulations. These simulations include the entire displacer interior and demonstrate it is possible to model a complete engine cycle in less than an hour. The focus of this effort was to demonstrate it is possible to produce useful Stirling engine performance results in a time-frame short enough to impact design decisions. The combination of utilizing the latest 64-bit Opteron computer processors, fiber-optical Myrinet communications, dynamic meshing, and across zone partitioning has enabled solution times at least 240 times faster than previous attempts at simulating the axisymmetric Stirling engine. A comparison of the multidimensional results, calibrated one-dimensional results, and known experimental results is shown. This preliminary comparison demonstrates that axisymmetric simulations can be very accurate, but more work remains to improve the simulations through such means as modifying the thermal equilibrium regenerator models, adding fluid-structure interactions, including radiation effects, and incorporating mechanodynamics.

Dyson, Rodger W.↗
NASA NTRS2007

Fast Whole-Engine Stirling Analysis

An experimentally validated approach is described for fast axisymmetric Stirling engine simulations. These simulations include the entire displacer interior and demonstrate it is possible to model a complete engine cycle in less than an hour. The focus of this effort was to demonstrate it is possible to produce useful Stirling engine performance results in a time-frame short enough to impact design decisions. The combination of utilizing the latest 64-bit Opteron computer processors, fiber-optical Myrinet communications, dynamic meshing, and across zone partitioning has enabled solution times at least 240 times faster than previous attempts at simulating the axisymmetric Stirling engine. A comparison of the multidimensional results, calibrated one-dimensional results, and known experimental results is shown. This preliminary comparison demonstrates that axisymmetric simulations can be very accurate, but more work remains to improve the simulations through such means as modifying the thermal equilibrium regenerator models, adding fluid-structure interactions, including radiation effects, and incorporating mechanodynamics.

Dyson, Rodger W.↗
DOE OSTI2021

Climate reconstructions for the Last Glacial Maximum from a simple cirque glacier in Fiordland, New Zealand

Mountain glacier records offer important constraints on the timing and magnitude of climate variability during the last glacial cycle. Existing moraine chronologies from the central Southern Alps indicate maximum ice extent was achieved during marine isotope stages 3–4, followed by repeated advances of similar, but gradually declining extent during marine isotope stage 2, until onset of the glacial termination. Questions remain over the precise role of climate in driving these changes, as most existing moraine chronologies come from large, complex former valley glacier systems, where non-climatic influences such as changing bed topography and proglacial lake formation may have influenced glacier length changes. In this work, we address this problem via a new cosmogenic 10Be chronology and equilibrium line altitude reconstruction from a cirque glacier situated in Fiordland, New Zealand. Our chronology shows moraine deposition at 32 ± 11 ka, 18.7 ± 0.2 ka, 18.1 ± 0.1 ka, and c. 17.2 ± 0.3 ka. The simple geometry of the former glacier supports the role of climate in driving a net decline of regional ice volume during marine isotope stages 3–2. Close spacing and good preservation of the 19–17 ka moraines permits 2D glacier reconstruction which suggests the equilibrium line altitude remained depressed by c.1130 m (equivalent to 5.8 ± 0.6 °C colder than present) during this interval. Onset of warming after 17.2 ± 0.2 ka is consistent with climate proxy evidence for a sustained southward shift in the southern westerly winds, which may have promoted deglaciation via shifting ocean currents and promoting increases in atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

The Scientific Case for Sample Return from a Primitive Near-Earth Asteroid

This paper is predicated on the recent experimental findings that the smallest structural ferromagnesiosilica entities in collected chondritic aggregate interplanetary dust particles (IDPs) have a predictable metastable eutectic composition. Kinetically controlled gas to solid condensation of MgFe-SiO-O2-H2 vapors does NOT produce stoichiometric crystalline solids (i.e. minerals) such as predicted by equilibrium condensation models but instead yields amorphous solids that are chemically ordered at metastable eutectics in the binary phase diagrams. Therefore these condensed dust grains will have a considerable amount of 'internal free energy' that will make them highly responsive to changes in their environments. Their inherent high energy-content 'buys' time and energy for mineralogical modification and chemical readjustment in response to changing environmental conditions. The activation energy barrier for reactions in metastable eutectic solids will be lower than for crystalline solids. As a result the reactions can take place at much lower temperatures and on much shorter time scales. That is, they will also occur much earlier in the evolution of a parent body wherein heat-producing sources may either be immature or inefficient. The ensuing reaction chains from metastable equilibrium to full thermodynamic equilibration with the local environment will be a chaotic 'Ostwald cascade' but from a well defined starting point to a predictable end result. Here we will discuss the implications of metastable eutectic dust as we see them for dust properties in icy and ice-free parent bodies wherein circumstellar dust is still recognizable because of the primitive nature of these parent bodies. We point to potential engineering constraints on sample acquisition and storage during Earth transit. Among the scientific goals of a primitive Near Earth Asteroid (NEA) sample return mission will be verification of the nature of dust forming and modification processes during hierarchical accretion in the earliest protoplanets. We propose sampling an infrared P- or D-class NEA or an object showing cometary activity such as 2201 Oljato that could be an asteroid or evolved comet. Additional information is contained in the original extended abstract.

Rietmeijer, Frans J. M.↗
NASA NTRS1994

Rapid adaptation to Coriolis force perturbations of arm trajectory

1. Forward reaching movements made during body rotation generate tangential Coriolis forces that are proportional to the cross product of the angular velocity of rotation and the linear velocity of the arm. Coriolis forces are inertial forces that do not involve mechanical contact. Virtually no constant centrifugal forces will be present in the background when motion of the arm generates transient Coriolis forces if the radius of body rotation is small. 2. We measured the trajectories of arm movements made in darkness to a visual target that was extinguished as movement began. The reaching movements were made prerotation, during rotation at 10 rpm in a fully enclosed rotating room, and postrotation. During testing the subject was seated at the center of the room and pointed radially. Neither visual nor tactile feedback about movement accuracy was present. 3. In experiment 1, subjects reached at a fast or slow rate and their hands made contact with a horizontal surface at the end of the reach. Their initial perrotary movements were highly significantly deviated relative to prerotation in both trajectories and end-points in the direction of the transient Coriolis forces that had been generated during the reaches. Despite the absence of visual and tactile feedback about reaching accuracy, all subjects rapidly regained straight movement trajectories and accurate endpoints. Postrotation, transient errors of opposite sign were present for both trajectories and endpoints. 4. In a second experiment the conditions were identical except that subjects pointed just above the location of the extinguished target so that no surface contact was involved. All subjects showed significant initial perrotation deviations of trajectories and endpoints in the direction of the transient Coriolis forces. With repeated reaches the trajectories, as viewed from above, again became straight, but there was only partial restoration of endpoint accuracy, so that subjects reached in a straight line to the wrong place. Aftereffects of opposite sign were transiently present in the postrotary movements. 5. These observations fail to support current equilibrium point models, both alpha and lambda, of movement control. Such theories would not predict endpoint errors under our experimental conditions, in which the Coriolis force is absent at the beginning and end of a movement. Our results indicate that detailed aspects of movement trajectory are being continuously monitored on the basis of proprioceptive feedback in relation to motor commands. Adaptive compensations can be initiated after one perturbation despite the absence of either visual or tactile feedback about movement trajectory and endpoint error. Moreover, movement trajectory and end-point can be remapped independently.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗
DOE OSTI2023

Variable renewable energy deployment in low-emission scenarios: The role of technology cost and value

While rapid deployment of variable renewable energy (VRE) technologies, namely wind and solar PV, is often projected in 2C pathways generated by integrated assessment models, there is a wide range in projected VRE deployment by mid-century. Such differences could be the result of differences in assumptions about future technology costs and/or differences in model approaches for capturing other aspects of technology competitiveness. Here we introduce a consistent competitiveness metric, profitability-adjusted levelized cost of electricity (or PLCOE), to an integrated assessment model (EPPA) to evaluate the representation of technology competition, including VRE, in low-emission scenarios. We show that representing the value of technology (alongside cost) may significantly impact VRE deployment relative to scenarios without such an adjustment. In addition, we show that varying VRE costs by about 35% in 2050 results in differences in VRE deployment that span much of the range in outcomes (over the same period) observed in likely 2C scenarios assessed by the IPCC, suggesting that both cost and value are key drivers of VRE deployment in such scenarios. Given the central role that VRE technologies play in the electricity mix across most scenarios, we also find that alternative cost assumptions for VRE technologies can lead to changes in electricity prices, the associated demand for electricity, and total final and primary energy consumption. However, the demand for fuels other than electricity is relatively insensitive to VRE assumptions in the 2C scenarios considered here.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

The Role of Deep Convection and Large-scale Circulation in Driving Model Spread in Low Cloud Feedback and Equilibrium Climate Sensitivity

This project aims to advance the understanding of the processes that drive the large uncertainties in climate change projections, use observations to constrain model physics and reduce the inter-model spread in equilibrium climate sensitivity (ECS). There are three major goals: 1) Characterize the representation of the physical pathways that link deep convection, large-scale circulation and low cloud feedback in CMIP6 model simulations and determine the relative contribution of each pathway to the CMIP6 model spread in low cloud feedback and ECS; 2) Use process-oriented diagnostics and multiple observations to evaluate CMIP6 model performance in capturing the observed cloud-circulation relation and deep convection characteristics including convective transition statistics and the bulk properties of mesoscale convective systems (MCSs). Error decomposition in CMIP6 models will be performed. 3) Conduct E3SM short-range hindcasts following the DOE Cloud-Associated Parameterizations Testbed (CAPT) protocol to pinpoint specific model parameters/processes that are crucial to the representation of deep convection, circulation, clouds and the pathways that connect them. We will modify convective parameters in E3SM and analyze the perturbed physics experiments (PPEs) to isolate model parameters that are critical to the uncertainty of ECS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Hydrazine Energy Storage: Displacing N 2 H 4 from the Metal Coordination Sphere

Hydrogen carriers, such as hydrazine (N 2 H 4 ), may facilitate long duration energy storage, a vital component for resilient grids by enabling more renewable energy generation. Lanthanide coordination chemistry with N 2 H 4 as well as efforts to displace N 2 H 4 from the metal coordination sphere to develop an efficient catalytic production cycle were detailed. Modeling the equilibrium of different ligand coordination, it was predicted that strong sigma donor molecules would be required to displace N 2 H 4 . Monitoring competition experiments with nuclear magnetic resonance confirmed that trimethyl phosphine oxide, dimethylformamide, and dimethyl sulfoxide displaced N 2 H 4 in large or small lanthanide complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗
DOE OSTI2026

Investigation of the Surface Chemistry of Magnesium and Magnesium Alloy Powders

Magnesium (Mg) and Mg alloy powders are generally difficult to sinter at atmospheric pressure in inert atmospheres purportedly due to the presence of powder particle surface oxide layers. To better understand sintering of Mg and Mg alloy powders, x-ray photoelectron spectroscopy (XPS) was performed to characterize the surface chemistry of gas atomized pure Mg and Mg alloys AZ91D and WE43 powders. These results are used to quantitatively determine the surface chemistry of these powders. Additionally, equilibrium thermodynamic modeling of these Mg/Mg alloys materials exposed to oxygen down to part per billion concentrations are presented for relative comparison to experimental XPS results. Analysis of these results are intended to improve sintering of Mg and Mg alloy powder materials through more quantitative knowledge of powder surface chemistry.

25 ENERGY STORAGE↗
DOE OSTI2024

Melting and solidification dynamics during laser melting of reaction-based metal matrix composites uncovered by in-situ synchrotron X-ray diffraction

Laser additive manufacturing (AM) of reaction-based metal matrix composites (MMCs) involves highly complex and non-equilibrium material transformation behavior, including melting, dissolution, precipitation, and solidification. Yet, the dynamics and interplay of these phase transformation processes remain poorly understood, posing substantial challenges in identifying the microstructure formation mechanism, and predicting and controlling the microstructure in the printed parts. Here, in this study, we performed the in-situ X-ray diffraction experiment to characterize the phase evolution dynamics of the 316L + 10 vol.%TiC system during laser melting, which provides direct and quantitative insights of the complex phase reaction and evolution dynamics under rapid heating and cooling conditions relevant to additive manufacturing of reaction-based MMCs. Further in-depth thermodynamic and kinetic calculations revealed that most of the phase evolution behavior observed in the in-situ X-ray diffraction experiment cannot be solely explained by widely used equilibrium thermodynamic models, and diffusion-controlled nonequilibrium dissolution and precipitation kinetics must be considered to elucidate the complex phase evolution behavior, including incomplete TiC dissolution, and three-step TiC precipitation. The three distinct types of precipitates generate unique hierarchical TiC micro- and nanostructures, which enhances the yield strength from 513 MPa to 877 MPa by 71 %, tensile strength from 628 MPa to 1054 MPa by 68 %, and Young's modulus from 193 GPa to 221 GPa by 14 %. The findings of our research provide the knowledge foundation for the design of unique microstructures and advanced MMC materials through laser AM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the hot surface ignition of a wall-stagnating spray flame

The ignition of fuel sprays due to interaction with hot surfaces is an important phenomenon in the safety analysis of many engineering systems. We perform a parametric study of the hot surface ignition (HSI) of a fuel spray approaching a heated surface caused by the accidental leakage of a fuel line. Here, to this end, we employ a one-dimensional Eulerian-Eulerian formulation with a non-equilibrium evaporation model and a realistic chemical mechanism to describe $n$-dodecane fuel chemistry. We first describe and analyze the phenomenology of the unsteady processes leading to ignition using non-dimensionalized quantities. Through consideration of the temporal development of the most reactive mixture, we demonstrate that ignition occurs at a fuel-lean composition in a premixed region near the hot surface. Using non-dimensional parameters identified from the governing equations, we perform a parametric study of the time, location and local mixture composition at ignition and determine the ignition limits. We then identify the most important parametric sensitivities for physical analysis using a data-driven classification method. Our analysis demonstrates a contraction of the ignition limits with increased Stokes number and a regime of parametric insensitivity of igniting mixture composition. We also show that at high Damköhler numbers, the ignition location conforms to the parametric behavior of the thermal boundary layer, whereas at low Damköhler numbers approaching the ignition limit it reaches a near-unity value of the quenching Peclet number. We then compare the demonstrated parametric dependencies to the results of the quasi-steady asymptotic ignition literature, showing that our results are consistent with those obtained analytically within the limitations imposed by the simplified formulation of the latter.

42 ENGINEERING↗
DOE OSTI2020

Co-benefits of subnationally differentiated carbon pricing policies in China: Alleviation of heavy PM 2.5 pollution and improvement in environmental equity

Carbon policies differentiated according to regional disparities have been recognized and recommended. However, the air quality co-benefits of differentiated policies remain unclear. Here we combine a multi-regional dynamic computable general equilibrium (CGE) model and an extend response surface model (ERSM) to investigate the impact of nationally uniform and subnationally differentiated carbon pricing policies on fine particulate matter (PM 2.5 ) concentrations in China. All policies examined lead to the same accumulated national CO 2 emissions (2020–2050), with the 2030 emissions attaining China's national abatement target. We find that, when subnational policies differentiated according to provincial PM 2.5 concentrations are implemented, the higher-than-average carbon price over more polluted provinces results in twice as much reduction in CO 2 and air pollutant emissions in 2050 as the national policy. As a result, the subnational policies cause larger PM 2.5 concentration reductions in these provinces (9%–18%) than the national policy (7%–11%). Also, the subnational policies eliminate high PM 2.5 exposure of over 45 μg/m 3 which 12% people suffer from under the national policy. Overall, the subnational policies substantially reduce regional disparity in PM 2.5 pollution and hence improve environmental equity. The results suggest that subnationally differentiated carbon policies are a promising instrument to mitigate severe pollution and promote environmental equity.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

CFD unified approach under Eulerian–Lagrangian framework for methanol and gasoline direct injection sprays in evaporative and flash boiling conditions

Innovative synthetic fuels for advanced propulsion systems, such as methanol and ammonia, and synthetic blended fuels (E00, E10, and E30), known for their high volatility, are often injected directly into combustion chambers. It follows that Eulerian–Lagrangian spray models need to accurately capture the spray collapse as a consequence of flash boiling onset and be capable of proficiently handling the preferential evaporation of multi-component fuels in evaporative scenarios. So, we performed the assessment of an Eulerian–Lagrangian CFD code for simulating methanol and E00 gasoline blend sprays in both early and late injection conditions involving flash boiling conditions and preferential evaporation. The adoption of an effervescent breakup model and of a non-equilibrium phase transition model for the discrete phase allows the adoption of a setup that is almost completely free from specific constant tuning, especially for what concerns the breakup model. We validated the simulations using experimental PLV maps of methanol and E00 sprays issued from the ECN Spray M injector. The results highlight a significantly different morphology of the methanol spray compared to the E00 one under late injection conditions. Under stratified combustion, low-volatile fuels are likely to be ignited first, and the flame propagates toward the high-volatile fuels. In conclusion, the spray collapse was also correctly reproduced, inducing the presence of a low-pressure zone and modifying the spray morphology.

E00↗