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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Turbulence Model Predictions of Strongly Curved Flow in a U-Duct

The ability of three types of turbulence models to accurately predict the effects of curvature on the flow in a U-duct is studied. An explicit algebraic stress model performs slightly better than one- or two-equation linear eddy viscosity models, although it is necessary to fully account for the variation of the production-to-dissipation-rate ratio in the algebraic stress model formulation. In their original formulations, none of these turbulence models fully captures the suppressed turbulence near the convex wall, whereas a full Reynolds stress model does. Some of the underlying assumptions used in the development of algebraic stress models are investigated and compared with the computed flowfield from the full Reynolds stress model. Through this analysis, the assumption of Reynolds stress anisotropy equilibrium used in the algebraic stress model formulation is found to be incorrect in regions of strong curvature. By the accounting for the local variation of the principal axes of the strain rate tensor, the explicit algebraic stress model correctly predicts the suppressed turbulence in the outer part of the boundary layer near the convex wall.

Rumsey, Christopher L.↗
NASA NTRS2011

Molybdenum Valence in Basaltic Silicate Melts: Effects of Temperature and Pressure

The metal-silicate partitioning behavior of molybdenum has been used as a test for equilibrium core formation hypotheses [for example, 1-6]. However, current models that apply experimental data to equilibrium core-mantle differentiation infer the oxidation state of molybdenum from solubility data or from multivariable coefficients from metal-silicate partitioning data [1,3,7]. Molybdenum, a multi-valent element with a valence transition near the fO2 of interest for core formation (approx.IW-2) will be sensitive to changes in fO2 of the system and silicate melt structure. In a silicate melt, Mo can occur in either 4+ or 6+ valence state, and Mo(6+) can be either octahedrally or tetrahedrally coordinated. Here we present X-ray absorption near edge structure (XANES) measurements of Mo valence in basaltic run products at a range of P, T, and fO2 and further quantify the valence transition of Mo.

Danielson, L. R.↗
NASA NTRS2018

Recent Updates in the SMAP Level-4 Soil Moisture Algorithm

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level-4 Soil Moisture (L4_SM) product. The L4_SM data are published with a mean latency of ~2.5 days from the time of observation and provide global, three-hourly, 9-km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. The L4_SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF). In 2018, the L4_SM algorithm was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high-quality, space-borne remote sensing observations. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of surface soil moisture from below under non-equilibrium conditions, which brings the model's surface soil moisture more in line with the SMAP Level-2 and in situ soil moisture. Moreover, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4, and the "catchment deficit" model variable was removed from the EnKF state vector to avoid degrading the model's groundwater estimates.Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation-minus-forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4_SM product and the lessons learned from the effort.

Reichle, Rolf H.↗
NASA NTRS2013

Chandra Grating Spectroscopy of Three Hot White Dwarfs

High-resolution soft X-ray spectroscopic observations of single hot white dwarfs are scarce. With the Chandra Low-Energy Transmission Grating, we have observed two white dwarfs, one is of spectral type DA (LB1919) and the other is a non-DA of spectral type PG1159 (PG1520+525). The spectra of both stars are analyzed, together with an archival Chandra spectrum of another DA white dwarf (GD246). Aims. The soft X-ray spectra of the two DA white dwarfs are investigated in order to study the effect of gravitational settling and radiative levitation of metals in their photospheres. LB1919 is of interest because it has a significantly lower metallicity than DAs with otherwise similar atmospheric parameters. GD246 is the only white dwarf known that shows identifiable individual iron lines in the soft X-ray range. For the PG1159 star, a precise effective temperature determination is performed in order to confine the position of the blue edge of the GW Vir instability region in the HRD. Methods. The Chandra spectra are analyzed with chemically homogeneous as well as stratified NLTE model atmospheres that assume equilibrium between gravitational settling and radiative acceleration of chemical elements. Archival EUV and UV spectra obtained with EUVE, FUSE, and HST are utilized to support the analysis. Results. No metals could be identified in LB1919. All observations are compatible with a pure hydrogen atmosphere. This is in stark contrast to the vast majority of hot DA white dwarfs that exhibit light and heavy metals and to the stratified models that predict significant metal abundances in the atmosphere. For GD246 we find that neither stratified nor homogeneous models can fit the Chandra spectrum. The Chandra spectrum of PG1520+525 constrains the effective temperature to T(sub eff) = 150 000 +/- 10 000 K. Therefore, this nonpulsating star together with the pulsating prototype of the GWVir class (PG1159-035) defines the location of the blue edge of the GWVir instability region. The result is in accordance with predictions from nonadiabatic stellar pulsation models. Such models are therefore reliable tools to investigate the interior structure of GW Vir variables. Conclusions. Our soft X-ray study reveals that the understanding of metal abundances in hot DA white dwarf atmospheres is still incomplete. On the other hand, model atmospheres of hydrogen-deficient PG1159-type stars are reliable and reproduce well the observed spectra from soft X-ray to optical wavelengths.

Adamczak, J.↗
NASA NTRS1988

Parametric study of grid size, time step and turbulence modeling on Navier-Stokes computations over airfoils

An upwind-biased implicit approximate factorization algorithm is applied to several steady and unsteady turbulent flows. The thin layer form of the compressible Navier-Stokes equation is used. Both the flux vector splitting and flux difference splitting methods are used to determine fluxes, and the results are compared. Flux difference splitting predicts results more accurately than flux vector splitting on a given mesh size, but, in its present implementation, is more severely limited by the maximum CFL number for unsteady time accurate flows. Physical aspects of the computations are also examined. An equilibrium turbulent boundary layer model computes generally better steady and unsteady results than a nonequilibrium model when there is little to no boundary layer separation. Conversely, when a significant region of separation exists, the nonequilibrium model performs in better agreement with experiment.

Rumsey, Christopher L.↗
NASA NTRS1982

An investigation of equilibrium concepts

A different approach to modeling of the thermochemistry of rocket engine combustion phenomena is presented. The methodology described is based on the hypothesis of a new variational principle applicable to compressible fluid mechanics. This hypothesis is extended to treat the thermochemical behavior of a reacting (equilibrium) gas in an open system.

Prozan, R. J.↗
NASA NTRS1977

Minimum-time acceleration of aircraft turbofan engines

Minimum-time accelerations of the F100 turbofan engine are presented. A piecewise-linear engine model, having three state variables and four control variables, is used to obtain the minimum-time solutions. The linear model which applies at a given time in the trajectory is determined by calculating a normalized distance from the current state to the equilibrium state associated with each linear model. The linear model associated with the closest equilibrium point is then used. The control histories for the minimum-time solutions are used as input to a nonlinear simulation of the F100 engine to verify the accuracy of the piecewise linear solutions.

Teren, F.↗
NASA NTRS1977

Minimum-time acceleration of aircraft turbofan engines

Minimum-time accelerations of the F100 turbofan engine are presented. A piecewise-linear engine model, having three state variables and four control variables, is used to obtain the minimum-time solutions. The linear model which applies at a given time in the trajectory is determined by calculating a normalized 'distance' from the current state to the equilibrium state associated with each linear model. The linear model associated with the closest equilibrium point is then used. The control histories for the minimum-time solutions are used as input to a nonlinear simulation of the F100 engine to verify the accuracy of the piecewise-linear solutions.

Teren, F.↗
NASA NTRS1982

Preliminary modeling of high pressure partial melting - Implications for early lunar differentiation

A quantitative model is presented for equilibrium crystallization of lunar magmas at high pressure, which, viewed in reverse, is equivalent to equilibrium partial melting. Data are extrapolated from simple systems, mare basalt systems, and from experiments on terrestrial compositions. Liquidus phase equilibria in a portion of the CaO-MgO-Al2O3-SiO2 system serve as a guide for predicting phase equilibria in the natural system. High pressure partial melting sequences are calculated for several proposed whole moon compositions at nominal pressure of 10 and 20 kb. Results show that compositions with chondritic Ca/Al have quite different melting paths from those with subchondritic Ca/Al, and as a result the types of magmas produced are quite different. If the moon possesses a chondritic Ca/Al ratio, models which rely on small degrees of partial melting to produce the early lunar magnetism are suspect.

Longhi, J.↗
NASA NTRS1992

Lumley's energy cascade dissipation rate model for boundary-free turbulent shear flows

True dissipation occurs mainly at the highest wavenumbers where the eddy sizes are comparatively small. These high wavenumbers receive their energy through the spectral cascade of energy starting with the largest eddies spilling energy into the smaller eddies, passing through each wavenumber until it is dissipated at the microscopic scale. However, a small percentage of the energy does not spill continuously through the cascade but is instantly passed to the higher wavenumbers. Consequently, the smallest eddies receive a certain amount of energy almost immediately. As the spectral energy cascade continues, the highest wavenumber needs a certain time to receive all the energy which has been transferred from the largest eddies. As such, there is a time delay, of the order of tau, between the generation of energy by the largest eddies and the eventual dissipation of this energy. For equilibrium turbulence at high Reynolds numbers, there is a wide range where energy is neither produced by the large eddies nor dissipated by viscosity, but is conserved and passed from wavenumber to higher wavenumbers. The rate at which energy cascades from one wavenumber to another is proportional to the energy contained within that wavenumber. This rate is constant and has been used in the past as a dissipation rate of turbulent kinetic energy. However, this is true only in steady, equilibrium turbulence. Most dissipation models contend that the production of dissipation is proportional to the production of energy and that the destruction of dissipation is proportional to the destruction of energy. In essence, these models state that the change in the dissipation rate is proportional to the change in the kinetic energy. This assumption is obviously incorrect for the case where there is no production of turbulent energy, yet energy continues to cascade from large to small eddies. If the time lag between the onset on the energy cascade to the destruction of energy at the microscale can be modeled, then there will be a better representation of the dissipation process. Development of an energy cascade time scale equation is discussed.

Duncan, B. S.↗
NASA NTRS1987

Navier-Stokes and viscous shock-layer solutions for radiating hypersonic flows

Results are presented from the Navier-Stokes and viscous shock-layer (VSL) calculations with nonequilibrium and equilibrium chemistry, respectively. These calculations contain coupling to the Aerotherm radiation code RAD. A simplified form of the electron energy equation is used to obtain an electron temperature in the Navier-Stokes calculations. The radiation in the flowfield is calculated using this temperature. The Navier-Stokes code is used at high altitude only, whereas the VSL code is employed for the entire entry period to make estimates of the radiative and convective heating to the Fire II vehicle. Results from the Navier-Stokes code have also been compared with the predictions of Lee and Kawamura, who used gray-gas radiation model and thin-layer Navier-Stokes equations. Quite good agreement is obtained between the measured and computed values of radiative and convective heating from the VSL code in th medium-to-low altitude flight regime of the Fire II vehicle. At high altitudes, the Navier-Stokes calculations considerably overpredict the Fire II flight data for radiative intensity. This is attributed to the deficiencies in the Aerotherm radiation model when used for low-density flight conditions. This model contains the thermal equilibrium assumption and precludes accounting for the collision-limiting phenomenon at high altitudes. Present Navier-Stokes calculations highlight the effect of these assumptions on radiative heating calculations for such conditions.

Gupta, Roop N.↗
NASA NTRS1986

The ortho/para ratio of water vapor in Comet Halley

The ortho/para ratio of H2O is shown to be an invariant in the cometary coma. The dependence of ortho-para ratio on temperature in thermal equilibrium is given, and the nuclear-spin-temperature is defined. Its relation to the physical temperature of the cometary ices is discussed, and the prospects for using the observed ortho/para ratio to infer properties of the cometary nucleus are explored. The ortho/para ratio in Halley's comet is derived from high resolution infrared spectra of near 2.7 microns wavelength. On UT December 24.1, 1985 it was 2.73 + or - 0.17, and on UT March 22.7, 1986 it was 3.23 + or - 0.37. The nuclear-spin-temperature was 35 K (+9 K, -5 K) pre-perihelion, and less than 40 K post-perihelion, at the 67% confidence limit. Both numbers are consistent with modeled values of the equilibrium temperature of the cometary nucleus at aphelion (47 K). However, at the 95% confidence limit they are also fully consistent with temperatures less than 50 K, corresponding to an ortho/para ratio of about 3.0.

Mumma, Michael J.↗
NASA NTRS2003

Various Numerical Applications on Tropical Convective Systems Using a Cloud Resolving Model

In recent years, increasing attention has been given to cloud resolving models (CRMs or cloud ensemble models-CEMs) for their ability to simulate the radiative-convective system, which plays a significant role in determining the regional heat and moisture budgets in the Tropics. The growing popularity of CRM usage can be credited to its inclusion of crucial and physically relatively realistic features such as explicit cloud-scale dynamics, sophisticated microphysical processes, and explicit cloud-radiation interaction. On the other hand, impacts of the environmental conditions (for example, the large-scale wind fields, heat and moisture advections as well as sea surface temperature) on the convective system can also be plausibly investigated using the CRMs with imposed explicit forcing. In this paper, by basically using a Goddard Cumulus Ensemble (GCE) model, three different studies on tropical convective systems are briefly presented. Each of these studies serves a different goal as well as uses a different approach. In the first study, which uses more of an idealized approach, the respective impacts of the large-scale horizontal wind shear and surface fluxes on the modeled tropical quasi-equilibrium states of temperature and water vapor are examined. In this 2-D study, the imposed large-scale horizontal wind shear is ideally either nudged (wind shear maintained strong) or mixed (wind shear weakened), while the minimum surface wind speed used for computing surface fluxes varies among various numerical experiments. For the second study, a handful of real tropical episodes (TRMM Kwajalein Experiment - KWAJEX, 1999; TRMM South China Sea Monsoon Experiment - SCSMEX, 1998) have been simulated such that several major atmospheric characteristics such as the rainfall amount and its associated stratiform contribution, the Qlheat and Q2/moisture budgets are investigated. In this study, the observed large-scale heat and moisture advections are continuously applied to the 2-D model. The modeled cloud generated from such an approach is termed continuously forced convection or continuous large-scale forced convection. A third study, which focuses on the respective impact of atmospheric components on upper Ocean heat and salt budgets, will be presented in the end. Unlike the two previous 2-D studies, this study employs the 3-D GCE-simulated diabatic source terms (using TOGA COARE observations) - radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the Ocean mixed-layer (OML) model.

Shie, C.-L.↗
NASA NTRS2003

A Proposed Mechanism for the Thermal Denaturation of a Recombinant Bacillus Halmapalus Alpha-amylase - the Effect of Calcium Ions

The thermal stability of a recombinant alpha-amylase from Bacillus halmapalus alpha-amylase (BHA) has been investigated using circular dichroism spectroscopy (CD) and differential scanning calorimetry (DSC). This alpha-amylase is homologous to other Bacillus alpha-amylases where previous crystallographic studies have identified the existence of 3 calcium binding sites in the structure. Denaturation of BHA is irreversible with a Tm of approximately 89 C, and DSC thermograms can be described using a one-step irreversible model. A 5 C increase in T(sub m) in the presence of 10 fold excess CaCl2 was observed. However, a concomitant increase in the tendency to aggregate was also observed. The presence of 30-40 fold excess calcium chelator (EDTA or EGTA) results in a large destabilization of BHA corresponding to about 40 C lower T(sub m), as determined by both CD and DSC. Ten fold excess EGTA reveals complex DSC thermograms corresponding to both reversible and irreversible transitions, which possibly originate from different populations of BHA:calcium complexes. The observations in the present study have, in combination with structural information of homologous alpha-amylases, provided the basis for the proposal of a simple denaturation mechanism of BHA. The proposed mechanism describes the irreversible thermal denaturation of different BHA:calcium complexes and the calcium binding equilibrium involved. Furthermore, the model accounts for a temperature induced reversible structural change associated with calcium binding.

Nielsen, Anders D.↗
NASA NTRS2023

Modeling the Chromosphere and Transition Region of Planet-hosting Star GJ 436

Ahead of upcoming space missions intending to conduct observations of low-mass stars in the ultraviolet (UV) spectral region it becomes imperative to simultaneously conduct atmospheric modeling from the UV to the visible (VIS) and near-infrared (NIR). Investigations on extended spectral regions will help to improve the overall understanding of the diversity of spectral lines arising from very different atmospheric temperature regions. Here we investigate atmosphere models with a chromosphere and transition region for the M2.5V star GJ 436, which hosts a close-in Hot Neptune. The atmosphere models are guided by observed spectral features from the UV to the VIS/NIR originating in the chromosphere and transition region of GJ 436. High-resolution observations from the Hubble Space Telescope and Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs (CARMENES) are used to obtain an appropriate model spectrum for the investigated M dwarf. We use a large set of atomic species considered in nonlocal thermodynamic equilibrium conditions within our PHOENIX model computations to approximate the physics within the low-density atmospheric regions. In order to obtain an overall match for the nonsimultaneous observations, it is necessary to apply a linear combination of two model spectra, where one of them better reproduces the UV lines while the other better represents the lines from the VIS/NIR range. This is needed to adequately handle different activity states across the observations.

stars↗
NASA NTRS1973

Dissociation rate of bromine diatomics in an argon heat bath

The evolution of a collection of 300 K bromine diatomics embedded in a heat bath of argon atoms at 1800 K was studied by computer, and a dissociation-rate constant for the reaction Br2 + BR + Ar yields Br + Ar was determined. Previously published probability distributions for energy and angular momentum transfers in classical three-dimensional Br2-Ar collisions were used in conjunction with a newly developed Monte Carlo scheme for this purpose. Results are compared with experimental shock-tube data and the predictions of several other theoretical models. A departure from equilibrium is obtained which is significantly greater than that predicted by any of these other theories.

Razner, R.↗
NASA NTRS1977

Some fundamental aspects of the fatigue and fracture response of composite materials

A simple analytical model, based on the concept that the constraining plies surrounding cracks through the thickness of a given ply or plies control the spacing of those cracks by controlling the surrounding stress state, appears to agree well with experimental data. Such a model can be used to evaluate critically various constraint situations including the prediction of the ability of a given constraint layer to arrest flaw development. The analytical model, which uses an equilibrium element approach in the present case, is sensitive to lamina and laminate elastic properties, flaw size, stacking sequence, and an interlaminar stress transfer distance which is estimated from experimental observations. The analysis is easily extended to study flaw growth and plastic region formation.

Reifsnider, K. L.↗