Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete equation of state”

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 235 records · Page 13

The Influence of Shear on Deep Convection Initiation. Part I: Theory

Abstract This article introduces a novel hypothesis for the role of vertical wind shear (“shear”) in deep convection initiation (DCI). In this hypothesis, initial moist updrafts that exceed a width and shear threshold will “root” within a progressively deeper steering current with time, increase their low-level cloud-relative flow and inflow, widen, and subsequently reduce their susceptibility to entrainment-driven dilution, evolving toward a quasi-steady self-sustaining state. In contrast, initial updrafts that do not exceed the aforementioned thresholds experience suppressed growth by shear-induced downward pressure gradient accelerations, will not root in a deep-enough steering current to increase their inflow, will narrow with time, and will succumb to entrainment-driven dilution. In the latter case, an externally driven lifting mechanism is required to sustain deep convection, and deep convection will not persist in the absence of such lifting mechanism. A theoretical model is developed from the equations of motion to further explore this hypothesis. The model indicates that shear generally suppresses DCI, raising the initial subcloud updraft width that is necessary for it to occur. However, there is a pronounced bifurcation in updraft growth in the model after the onset of convection. Sufficiently wide initial updrafts grow and eventually achieve a steady state. In contrast, insufficiently wide initial updrafts shrink with time and eventually decay completely without external support. A sharp initial updraft radius threshold discriminates between these two outcomes. Thus, consistent with our hypothesis and observations, shear inhibits DCI in some situations, but facilitates it in others.

42 ENGINEERING↗

The interior structure of Jupiter (Consequences of Pioneer 10 data)

Models of the Jovian interiors based on theoretical equations of state of hydrogen and helium supported by a few experimental points and on observed parameters such as oblateness, gravitational coefficients, heat emission, magnetic fields, are discussed. The models fall into three categories: (1) those that assume a uniform and rather low H2/He ratio throughout the planet; (2) those in which this ratio is solar and thus higher; and (3) those that take into account the lack of complete miscibility of the two elements in the condensed state. It appears now also that within the limits of error the planet is in a hydrostatic equilibrium. The large heat emission and the need for an efficient source of internal heat are confirmed, but the results do not indicate which one of the various possible mechanisms is favored, although new evolutionary models suggest that the primordial heat may be insufficient. A new red spot has been discovered. Finally, the presence of a highly eccentric and inclined magnetic field poses new problems related to the pattern of internal convection and to the possibility of a north-south asymmetry of the interior. Further analysis of the available data may throw additional light on these questions.

Smoluchowski, R.↗

Study of the agreement between binary neutron star ejecta models derived from numerical relativity simulations

Neutron star mergers have recently become a tool to study extreme gravity, nucleosynthesis, and the chemical composition of the Universe. To date, there has been one joint gravitational and electromagnetic observation of a binary neutron star merger, GW170817, as well as a solely gravitational observation, GW190425. In order to accurately identify and interpret electromagnetic signals of neutron star mergers, better models of the matter outflows generated by these mergers are required. We compare a series of ejecta models to see where they provide strong constraints on the amount of ejected mass expected from a system, and where systematic uncertainties in current models prevent us from reliably extracting information from observed events. We also examine 2396 neutron star equations of state compatible with GW170817 to see whether a given ejecta mass could be reasonably produced with a neutron star of said equation of state, and whether different ejecta models provide consistent predictions. We find that the difference between models is often comparable to or larger than the error generally assumed for these models, implying better constraints on the models are needed. In conclusion, we also note that the extrapolation of outflow models outside of their calibration window, while commonly needed to analyze gravitational wave events, is extremely unreliable and occasionally leads to completely unphysical results.

79 ASTRONOMY AND ASTROPHYSICS↗

A steady-state stratospheric aerosol model

This paper deals with the development of a one-dimensional steady-state stratospheric aerosol model and the subsequent perturbations caused by including the expected Space Shuttle particulate effluents in the model. Two approaches to the basic modeling effort have been made: in one, enough simplifying assumptions were introduced so that a more or less exact solution to the descriptive equations could be obtained; in the other, very few simplifications were made and a computer technique was used to solve the equations. The most complete form of the model contains the effects of sedimentation, diffusion, particle growth and coagulation. The results indicate that the model is capable of describing many aspects of the stratospheric aerosol layer, such as size distribution and the vertical profile of particles greater than 0.3 micron diameter.

Rosen, J. M.↗

A Martian global groundwater model

A global groundwater flow model was constructed for Mars to study hydrologic response under a variety of scenarios, improving and extending earlier simple cross sectional models. The model is capable of treating both steady state and transient flow as well as permeability that is anisotropic in the horizontal dimensions. A single near surface confining layer may be included (representing in these simulations a coherent permafrost layer). Furthermore, in unconfined flow, locations of complete saturation and seepage are determined. The flow model assumes that groundwater gradients are sufficiently low that DuPuit conditions are satisfied and the flow component perpendicular to the ground surface is negligible. The flow equations were solved using a finite difference method employing 10 deg spacing of latitude and longitude.

Howard, Alan D.↗

New model atmospheres for very cool white dwarfs with mixed H/He and pure He compositions

New model atmosphere calculations for very cool white dwarfs with mixed H/He and pure He compositions are presented. The hydrogen-rich models incorporate improved cross section calculations of the collision-induced absorption by molecular hydrogen due to collisions with H2, H, and He. The effects associated with variations in the effective temperature (T(sub eff) greater than or equal to 4000 and less than or equal to 10,000 K), the surface gravity (log g greater than or equal to 7.5 and less than or equal to 9.5), and the chemical composition (N(He)/N(H) greater than or equal to O and less than or equal to 100) are investigated. Results from earlier calculations are confirmed qualitatively, but a more detailed comparison reveals large quantitative deviations. Cool white dwarfs with mixed H/He chemical compositions are shown to be easily recognizable from their predicted strong infrared flux deficiency. Pure helium model calculations are described as well. These include a modified version of the recently developed equation of state of D. Saumon and G. Chabrier. Nonideal effects brought about by various equations of state are explored in detail. For the purpose of this analysis, a model of pressure ionization based on an accurate description of the interactions in a mostly atomic helium fluid is developed. The effects of pressure ionization are shown to be the most important issue in the model calculations. A critical discussion of previous generations of pure helium model calculations is presented. Finally, broadband color indices are provided for the complete model grid.

Bergeron, P.↗

Investigation A: Numerical solutions of driven vortices of binary fluid in a cylindrical geometry

The steady state laminar motion of a viscous, incompressible and binary fluid is studied for a rotating flow in a cylindrical geometry. The mathematical model employed is a cold flow simulation of the fluid mechanics of the light-bulb concept of the gaseous core nuclear engine. A numerical treatment is developed for the rotating flow which includes a description of the nuclear fuel addition. The problem is formulated with the complete Navier-Stokes equations in order to show the interaction between the fuel addition, the main flow, and the boundary layer flow in an accurate manner. The results presented show holdup of the nuclear fuel for the case of steady fuel addition.

Shin, Y. W.↗

A vectorized, finite-volume, adaptive-grid algorithm for Navier-Stokes calculations

An adaptive grid, finite-volume method has been used to solve the Navier-Stokes equations for complete (forebody and afterbody) flowfields around blunt bodies. The code, which is applicable for axisymmetric or two-dimensional flows, allows the mesh to adjust during the computation to provide a closer spacing of mesh points in regions of high gradients, thus minimizing the number of required computational points. The solution technique is explicit, utilizing a maximum time-step advancement at each grid point to accelerate convergence to the steady state. The code has been fully vectorized for efficient solution on the CYBER 203 computer. A very flexible rezoning routine is used to concentrate mesh points anywhere in the field, either by a user-defined weighting function or by allowing high gradient regions to adjust the grid. The grid adjustment routine is implicit in nature and represents a very small portion of the total computational cost. Currently, the code runs in approximately 0.000016 seconds per grid point per iteration.

Gnoffo, P. A.↗

Interplay between charge and spin noise in the near-surface theory of decoherence and relaxation of C 3 v symmetry qutrit spin-1 centers

Decoherence and relaxation of solid-state defect qutrits near a crystal surface, where they are commonly used as quantum sensors, originate from charge and magnetic field noise. A complete theory requires a formalism for decoherence and relaxation that includes all Hamiltonian terms allowed by the defect's point-group symmetry. This formalism, presented here for the C 3 v symmetry of a spin-1 defect in a diamond, silicon carbide, or similar host, relies on a Lindblad dynamical equation and clarifies the relative contributions of charge and spin noise to relaxation and decoherence, along with their dependence on the defect spin's depth and resonant frequencies. The calculations agree with the experimental measurements of Sangtawesin [], and corroborate the importance of charge noise. Published by the American Physical Society 2024

Materials Science↗

Boundary conditions for implicit solutions to the compressible Navier-Stokes equations in finite computational domains

The admissible boundary conditions for the unsteady Navier-Stokes equations are insufficient to determine completely the flow variables at grid points on the boundaries of a finite computational domain. A general method is presented in which the boundary conditions are supplemented by conservatively differenced forms of the flow equations. For implicit algorithms, the method provides a fully implicit coupling between boundary and interior points. Implicit boundary-point computation schemes are given for most types of flow boundary encountered in internal and external flow problems. Convergence to steady state is generally faster than with locally explicit boundary-point schemes.

Thomas, P. D.↗

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie, close to the Hugoniot. Calcite 3 to 2 transition, upon release, was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90% devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

Combustion of liquid fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both sub- and super-critical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influences on the fluid transport, gas/liquid interface thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibits a significant variation near the critical burning pressure, mainly as a result of reduced mass-diffusion rate and latent heat of vaporization with increased pressure. The influence of droplet size on the burning characteristics is also noted.

Shuen, J. S.↗

Vaporization and combustion of fuel droplets at supercritical conditions

Vaporization and combustion liquid-fuel droplets in both sub- and super-critical environments have been examined. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. In addition, the dynamic responses of droplet vaporization and combustion to ambient-pressure oscillations are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗

Gauge-fixing quantum density operators at scale

We provide a theory, algorithms, and simulations of nonequilibrium quantum systems using a one-dimensional (1D) completely positive (CP), matrix-product (MP) density-operator (𝜌) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MP⁢𝜌 factorization naturally emerges. In this setting, we analytically and numerically examine the virtual gauge freedoms associated with the representation of quantum density operators. Based on this perspective, we simplify algorithms in certain limits to speed up the integration of the canonical-form master-equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance in the near term. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, because of loss from environmental interactions, nonlocal complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional nonequilibrium settings, beyond qubit engineering.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie close to the Hugoniot. Calcite 3 to 2 transition, upon release was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90 percent devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

Combustion of liquid-fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both subcritical and supercritical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates variable thermophysical properties, finite-rate chemical kinetics, and a full treatment of liquid-vapor phase equilibrium at the drop surface. The governing equations and associated interfacial boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures in the range of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influence on fluid transport, gas-liquid interfacial thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibit a considerable change in the combustion mechanism at the critical pressure, mainly as a result of reduced mass diffusivity and latent heat of vaporization with increased pressure.

Shuen, J. S.↗

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The 2-7 May 1998 Storm

A complete description of a self-consistent model of magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves and back on waves are considered self-consistently by solving both equations on a global magnetospheric scale under nonsteady state conditions. The developed model is employed to simulate the entire 2-7 May 1998 storm period. First, the trapped number fluxes of the ring current protons are calculated and presented along with comparison with the data measured by the three- dimensional hot plasma instrument Polar/HYDRA. Incorporating in the model the wave-particle interaction leads to much better agreement between the experimental data and the model results. Second, examining of the wave (MLT, L shell) distributions produced by the model during the storm progress reveals an essential intensification of the wave emission about 2 days after the main phase of the storm. This result is well consistent with the earlier ground-based observations. Finally, the theoretical shapes and the occurrence rates of the wave power spectral densities are studied. It is found that about 2 days after the storm s main phase on 4 May, mainly non-Gaussian shapes of power spectral densities are produced.

Khazanov, G. V.↗