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

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ice Supersaturation Variability in Cirrus Clouds: Role of Vertical Wind Speeds and Deposition Coefficients

Abstract Aircraft measurements reveal ice supersaturation statistics in cirrus (ISSs) with broad maxima around ice saturation and pronounced variance. In this study, processes shaping ISSs in midlatitude and tropical upper tropospheric conditions are systematically investigated. Water vapor deposition and sublimation of size‐resolved ice crystal populations are simulated in an air parcel framework. Mesoscale temperature fluctuations (MTFs) due to gravity waves force the temporal evolution of supersaturation. Various levels of background wave forcing and cirrus thickness are distinguished in stochastic ensemble simulations. Kinetic limitations to ice mass growth are brought about by supersaturation‐dependent deposition coefficients that represent efficient and inefficient growth modes as a function of ice crystal size. The simulations identify a wide range of deposition coefficients in cirrus, but most values stay above 0.01 such that kinetic limitations to water uptake remain moderate. Supersaturation quenching times are long, typically 0.5–2 hr. The wave forcing thus causes a remarkably large variability in ISSs and cirrus microphysical properties except in the thickest cirrus, producing ensemble‐mean ISSs in line with in‐situ measurements. ISS variance is controlled by MTFs and increases with decreasing cirrus integral radii. In comparison, the impact of ice crystal growth rates on ISSs is small. These results contribute to efforts directed at identifying and solving issues associated with ice‐supersaturated areas and non‐equilibrium cirrus physics in global models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Tax Credits for Clean Electricity: The Distributional Impacts of Supply-Push Policies in the Power Sector

We evaluate distributional and efficiency consequences of the bulk power clean electricity tax credits authorized by the 2022 Inflation Reduction Act. To do so, we link detailed electricity capacity expansion, computable general equilibrium, microsimulation, and air pollution models to estimate economic welfare and health incidence across demographic groups. We evaluate trade-offs between policy efficiency and income progressivity by comparing the tax credits to cap-and-trade policies. The tax credits encourage increased clean electricity investment, resulting in a reallocation of capital from elsewhere in the economy, higher prices for capital and other goods, lower power prices, and lower emissions. The tax credits yield progressive outcomes for economic welfare at the expense of efficiency while all modeled policies demonstrate progressivity in health impacts. The health benefits, absent climate benefits, exceed total policy costs and provide greater benefits for low-income and historically marginalized households given coincidence of household locations and emissions exposure intensity.

distributional impacts↗
DOE OSTI2022

Radiative pulsed L-mode operation in ARC-class reactors

Abstract A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to ∼ 90 % core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with P fus > 1000 MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and radiofrequency current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solar Thermochemical Carbon Dioxide Splitting Using Ceria and Iron Aluminate Foam Devices and Simulation of a Plant System for Demonstration

An international research project has been undertaken to integrate a unique solar thermal processing reactor system with ceria and iron aluminate as active redox materials for CO2 splitting. Experimental investigations for CO2 splitting were conducted using a solar simulator and tube furnace at Niigata University, followed by demonstrations using a high-flux solar furnace (HFSF) at the National Renewable Energy Laboratory (NREL) in Golden, CO. Each experimental setup consisted of foam devices composed of reticulated porous ceramic (RPC). The RPC has a full ceria or iron aluminate body. It fabricated using the replica method and subjected to a two-step redox reaction, which iteratively separated a stream of CO2 into O2 and CO. Reactivity was evaluated using CO production per mass of the reactive material. The tubular furnace yielded a CO production of 6.41 mL/g at a reduction temperature of 1600degrees C, showcasing a higher CO production rate and total amount than those obtained from experiments conducted with solar simulators and solar furnace setups. For iron aluminate RPC, the productivity was measured as 3.57 mL/g using HFSF at a reduction temperature of 1450degrees C. These results are somewhat higher than those of the previous experiment at lower reduction temperatures of 1400degrees C-1500degrees C. Additionally, the production of CO in the case of ceria RPC was compared with the steady flow model simulation, which assumed chemical equilibrium at various levels of oxygen partial pressure during the reduction process. On the basis of these results, this study proposes a solar fuel system with an open receiver that uses a high-temperature heat transfer fluid.

carbon dioxide thermochemical splitting↗
DOE OSTI2024

Non-Thermal Solar Wind Electron Velocity Distribution Function

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of the non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi-equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but, unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously emitted thermal fluctuations and the background turbulence leads to the formation of a non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗
NASA NTRS1971

Calculated lamellar binding. I

Van der Walls bound lamellar solids interlayer binding energy computational model, discussing talc and pyrophyllite equilibrium stacking arrangements and force constants

Phillips, J. M.↗
NASA NTRS1975

Polynomial coefficients of thermochemical data for the C-H-O-N system

Thermodynamic data are required input for the finite kinetics and equilibrium computer programs needed for modeling the combustion of hydrocarbons in the fields of energy and pollution research. Least squares determined coefficients of the curve-fitted thermodynamic data for 193 species in the C-H-O-N system are presented in card image form and are of suitable format for use by common computer programs.

Wakelyn, N. T.↗
NASA NTRS1975

Analysis of OSO data to determine the structure and energy balance of the solar chromosphere

A detailed reexamination of the temperature-density structure of the photosphere and low chromosphere shows that the middle and upper chromosphere, which directly emits most of the OSO spectrum, is sensitive to conditions in this underlying region of the atmosphere. A model of this region is based on a unified compilation of all recently published broadband flux and central intensity observations of the solar spectrum from 500 microns in the far infrared to 1220 A in the far ultraviolet. This extensive compilation includes the OSO 4 and 6 observations in the wavelength range 1400 to 1220 A. A model is presented of the quiet solar atmosphere in the height range between the temperature minimum and the upper part of the chromosphere-corona transition region. This model is based on statistical equilibrium calculations of H, He 1, He 2, Si 1, C 1, and other ions.

Avrett, E. H.↗
NASA NTRS1976

Critical examination of two-equation turbulence closure models

Comparison of the Jones-Launder, Ng-Spalding, Saffman-Wilcox, and Wilcox-Traci two-equation turbulence models has been conducted. It was shown that the Saffman-Wilcox and Wilcox-Traci dissipation-rate formulations admit straightforward integration through the viscous sublayer, whereas integration through the viscous sublayer is a more difficult issue with the Jones-Launder dissipation-function and the Ng-Spalding length-scale formulations. Numerical computations were conducted in which the models were applied to four equilibrium boundary layer flows including adverse, zero, and favorable pressure gradients. Computations of zero pressure gradient flow over a convex wall composed the final part of the comparison.

Chambers, T. L.↗
NASA NTRS1977

Departures of the electron energy distribution from a Maxwellian in hydrogen. I - Formulation and solution of the electron kinetic equation. II - Consequences

The problem of calculating the steady-state free-electron energy distribution in a hydrogen gas is considered in order to study departures of that distribution from a Maxwellian at sufficiently low degrees of ionization. A model kinetic equation is formulated and solved analytically for the one-particle electron distribution function in a steady-state partially ionized hydrogen gas, and it is shown that the formal solution can be accurately approximated by using the WKB method. The solutions obtained indicate that the high-energy tail of the distribution is susceptible to distortion by imbalanced inelastic collisions for ionization fractions not exceeding about 0.1 and that such departures from a Maxwellian can lead to significant changes in the collisional excitation and ionization rates of ground-state hydrogen atoms. Expressions for the electron-hydrogen collision rates are derived which explicitly display their dependence on the hydrogen departure coefficients. The results are applied in order to compare self-consistent predictions with those based on the a priori assumption of a Maxwellian distribution for models of the thermal ionization equilibrium of hydrogen in the optically thin limit, spectral-line formation by a gas consisting of two-level atoms, and radiative transfer in finite slabs by a gas of four-level hydrogen atoms.

Shoub, E. C.↗
NASA NTRS1977

Pyrolysis of organic compounds in the presence of ammonia The Viking Mars lander site alteration experiment

The influence of ammonia on the pyrolysis pattern of selected organic substances sorbed on an inorganic phase was investigated. The thermal degradation products were identified by gas chromatography-mass spectrometry. The feasibility of this technique was tested on a meteoritic sample. All substances examined react with ammonia at the pyrolysis temperature of 500 C, the major products being nitriles and heterocyclic compounds in which nitrogen was incorporated. Based on these results, a model for the non-equilibrium production of organic compounds on Jupiter is discussed. The investigation was performed in connection with the Viking lander molecular analysis. The results obtained indicate that the concentrations of ammonia in the retrorocket fuel exhaust would have been probably too small to produce significant changes in the Martian soil organic compounds if any were found.

Holzer, G.↗
NASA NTRS1978

Tropical sea surface temperature - An interactive one-dimensional atmosphere-ocean model

It is shown that the (cumulus) convective processes in the tropics may be described by a one-dimensional cloud model, while the near-surface ocean may similarly be described by a one-dimensional mixed-layer model. The coupling is achieved through a sea surface flux budget combined with the flux parameterizations implied by Monin-Obukhov similarity theory. The coupled one-dimensional atmosphere-ocean model is applied to the equilibrium situation in which all temperatures reach a steady state. For the ocean, the fluxes must vanish in equilibrium, but the atmosphere maintains a stable lapse rate by balancing cumulonimbus heating against net radiative cooling. All water precipitating from cumulonimbus clouds must have evaporated from the sea. It is shown that this equilibrium system is closed and determinable solely in terms of the solar constant.

Sarachik, E. S.↗
NASA NTRS1980

Hektor - Nature and origin of a binary asteroid

It is argued that the elongated single body model (Dunlap and Gehrels, 1969) which reproduces the rotational light curve of 624 Hektor is unstable with respect to binary fission. It is shown that a contact binary model with ellipsoidal components of equilibrium shape matches the observed amplitude. Hektor's bulk density is inferred to be about 2.5 g/cu cm.

Weidenschilling, S. J.↗
NASA NTRS1984

Amplitude distributions of the spider heartpulse in response to gravitational stimuli

The arachnid Nuctenea sclopetaria (Clerck) which possesses a neurogenic heart, measuring the heartbeat is under efferent control through a dorsal nerve arising from a brain center is discussed. It was shown that the heartrate of this spider is also modulated by an afferent input associated with small increments of gravity. A compressive force on the order of 40 micron is sufficient to elicit a threshold change in heart rate for a typical (100mg) spider. This obtains in a hyper-Gz field less than 1.001. The functional relationship between gravity and heartrate is logarithmic between the absolute threshold and at least 1.5 Gz. A model was proposed in which equilibrium and movement are maintained by changes in blood pressure. It is concluded that the arachnid equilibrium system is like a weight detector which employs a hydraulic compensatory mechanism.

Finck, A.↗
NASA NTRS1984

Protostellar formation in rotating interstellar clouds. IV Nonisothermal collapse

Radiative transfer in the Eddington approximation is included in a multidimensional, self-gravitational, hydrodynamical computer code. Details of the numerical solution and thermodynamic relations are given. Comparison calculations with previous spherically symmetrical models of protostellar collapse are used to validate the basic approach and the artifices which allow the explicit hydrodynamics code to follow the accretion of gas onto a quasi-equilibrium core. A series of axisymmetric models is used to investigate the importance of rotation in collapsing clouds, as the initial amount of angular momentum is lowered, with an emphasis on the possible formation of rings. Rings readily form even in the nonisothermal regime except for very low initial angular momenta; even these clouds may experience ring formation prior to reaching stellar densities. The models imply that other effects (such as gravitational torques or turbulent viscosity) may be necesary to avoid binary formation and thus result in a presolar nebula consistent with the assumptions of either Safronov or Cameron.

Boss, A. P.↗
NASA NTRS1986

Application of program LAURA to three-dimensional AOTV flowfields

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is a robust, finite volume, single-level storage, implicit upwind differencing algorithm which has been documented and tested on several three-dimensional blunt-body flows. The algorithm can run at unlimited Courant numbers (relaxing the steady-state equations) but requires the inversion of only a 5 x 5 matrix per computational cell. An alternating directional sweep Gauss-Seidel substitution strategy is used to relax the governing equations. At present, the Euler and thin-layer Navier-Stokes equations using Sutherland's law for viscosity have been modeled for a perfect gas, equilibrium air, and nonequilibrium air chemistry neglecting diffusion. The equilibrium and nonequilibrium air chemistry options have been described in a companion paper. Good comparisons with experimental data and another calculation method for pressure distributions, aerodynamic coefficients, and heat-transfer distributions have been demonstrated for three-dimensional blunt-body flows.

Gnoffo, P. A.↗