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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

A system debugging model.

Consideration of the nature of the 'debugging' process applied to a new complex system during the initial period of its life. During this period failures and errors are corrected as they occur, with resulting improvement in the subsequent performance of the system. One mathematical idealization of this process leads to the assumption that system failure rate is decreasing with time. In practice, the debugging phase is considered completed when the failure rate reaches an equilibrium or constant value. Models are formulated for this phenomenon. Maximum likelihood estimates are obtained for relevant failure rate functions and for the end of the debugging period. A conservative upper confidence bound on the stable failure rate is obtained.

Barlow, R. E.↗
DOE OSTI2023

Muted extratropical low cloud seasonal cycle is closely linked to underestimated climate sensitivity in models

Abstract A large spread in model estimates of the equilibrium climate sensitivity (ECS), defined as the global mean near-surface air-temperature increase following a doubling of atmospheric CO 2 concentration, leaves us greatly disadvantaged in guiding policy-making for climate change adaptation and mitigation. In this study, we show that the projected ECS in the latest generation of climate models is highly related to seasonal variations of extratropical low-cloud fraction (LCF) in historical simulations. Marked reduction of extratropical LCF from winter to summer is found in models with ECS > 4.75 K, in accordance with the significant reduction of extratropical LCF under a warming climate in these models. In contrast, a pronounced seasonal cycle of extratropical LCF, as supported by satellite observations, is largely absent in models with ECS < 3.3 K. The distinct seasonality in extratropical LCF in climate models is ascribed to their different prevailing cloud regimes governing the extratropical LCF variability.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Modeling and interpretation of the optical and HST UV spectrum of SN 1993J

We present spherically symmetric, non-local thermodynamic equilibrium (non-LTE), expanding model atmosphere synthetic spectra of SN 1993J and compare them to the ultraviolet and optical spectra obtained simultaneously on 1993 April 15 by the Hubble Space Telescope (HST) and at Lick Observatory. We are able to fit the optical and HST near-UV spectra with a variety of compositions, but find that our best fits are obtained with an enhanced helium abundance (Y = 0.8) and significant nonthermal ionization due to gamma rays produced by the decay of Ni-56. We find the effects of X-rays produced by thermal bremsstrahlung in the circumstellar material to have little effect on the spectrum. Our current models are unable to fit the HST far-UV region, predicting too little flux and strong, unseen Fe II features.

Baron, E.↗
NASA NTRS1974

The GISS model of the global atmosphere

A description and numerical results are presented for a global atmospheric circulation model developed at the Goddard Institute for Space Studies (GISS). The model version described is a 9-level primitive-equation model in sigma coordinates. It includes a realistic distribution of continents, oceans and topography. Detailed calculations of energy transfer by solar and terrestrial radiation make use of cloud and water vapor fields calculated by the model. The model hydrologic cycle includes two precipitation mechanisms: large-scale supersaturation and a parameterization of subgrid-scale cumulus convection. Results are presented both from a comparison of the 13th to the 43rd days (January) of one integration with climatological statistics, and from five short-range forecasting experiments. In the extended integration, the near-equilibrium January-mean model atmosphere exhibits an energy cycle in good agreement with observational estimates, together with generally realistic zonal mean fields of winds, temperature, humidity, transports, diabatic heating, evaporation, precipitation, and cloud cover.

Somerville, R. C. J.↗
NASA NTRS2014

The Virtual Observatory Service TheoSSA: Establishing a Database of Synthetic Stellar Flux Standards II. NLTE Spectral Analysis of the OB-Type Subdwarf Feige 110

Context. In the framework of the Virtual Observatory (VO), the German Astrophysical VO (GAVO) developed the registered service TheoSSA (Theoretical Stellar Spectra Access). It provides easy access to stellar spectral energy distributions (SEDs) and is intended to ingest SEDs calculated by any model-atmosphere code, generally for all effective temperatures, surface gravities, and elemental compositions. We will establish a database of SEDs of flux standards that are easily accessible via TheoSSA's web interface. Aims. The OB-type subdwarf Feige 110 is a standard star for flux calibration. State-of-the-art non-local thermodynamic equilibrium stellar-atmosphere models that consider opacities of species up to trans-iron elements will be used to provide a reliable synthetic spectrum to compare with observations. Methods. In case of Feige 110, we demonstrate that the model reproduces not only its overall continuum shape from the far-ultraviolet (FUV) to the optical wavelength range but also the numerous metal lines exhibited in its FUV spectrum. Results. We present a state-of-the-art spectral analysis of Feige 110. We determined Teff =47 250 +/- 2000 K, log g=6.00 +/- 0.20, and the abundances of He, N, P, S, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, and Ge. Ti, V, Mn, Co, Zn, and Ge were identified for the first time in this star. Upper abundance limits were derived for C, O, Si, Ca, and Sc. Conclusions. The TheoSSA database of theoretical SEDs of stellar flux standards guarantees that the flux calibration of astronomical data and cross-calibration between different instruments can be based on models and SEDs calculated with state-of-the-art model atmosphere codes.

Observatory↗
NASA NTRS1988

Numerical simulation of hypersonic inlet flows with equilibrium or finite rate chemistry

An efficient numerical program incorporated with comprehensive high temperature gas property models has been developed to simulate hypersonic inlet flows. The computer program employs an implicit lower-upper time marching scheme to solve the two-dimensional Navier-Stokes equations with variable thermodynamic and transport properties. Both finite-rate and local-equilibrium approaches are adopted in the chemical reaction model for dissociation and ionization of the inlet air. In the finite rate approach, eleven species equations coupled with fluid dynamic equations are solved simultaneously. In the local-equilibrium approach, instead of solving species equations, an efficient chemical equilibrium package has been developed and incorporated into the flow code to obtain chemical compositions directly. Gas properties for the reaction products species are calculated by methods of statistical mechanics and fit to a polynomial form for C(p). In the present study, since the chemical reaction time is comparable to the flow residence time, the local-equilibrium model underpredicts the temperature in the shock layer. Significant differences of predicted chemical compositions in shock layer between finite rate and local-equilibrium approaches have been observed.

Yu, Sheng-Tao↗
NASA NTRS1994

The chromospheric structure of the cool giant star g Herculis

Non-Local Thermodynamic Equilibrium (LTE) calculations of semiempirical chromospheric models are presented for 30 g Her (M6 III). This star is one of the coolest (T(sub eff) = 3250 K) SRb (semiregular) variable stars and has a mass perhaps as great as 4 solar mass. Chromospheric features we have observed in its spectrum include Mg II h and k; C II) UV0.01, which is sensitive to electron density; Mg I lambda 2852; Ca II H, K, and IRT; Ca I lambda 4227 and lambda 6573; Al II) UV 1; and H alpha. We pay special attention to fitting the C II intersystem lines and the Mg II resonance lines but use all the other features as constraints to some extent. The equations of radiative transfer and statistical equilibrium are solved self-consistently for H I, H(-), H2, He I, C I, C II, Na I, Mg I, Mg II, Al I, Al II, Ca I, and Ca II with the equivalent two-level technique. To simplify these calculations, a one-dimensional hydrostatic, plane-parallel atmosphere is assumed. We investigate 10 separate 'classical' chromospheric models, differing most importantly in total mass column density above the temperature minimum. Synthetic spectra from these models fit some but not all of the observations. These comparisons are discussed in detail. However, we find that no single-component classical model in hydrostatic equilibrium is able to reproduce both the Mg II line profiles and the relative strengths of the CII) lines. In all these models, chromospheric emission features are formed relatively close to the star (approximately less than 0.05 R(sub *). The circumstellar environment has a thick, cool component overlying the Mg II emission region, which is relatively static and very turbulent. Finally, we find that thermalization in the Mg II h and k lines in the coolest giant stars is controlled by continuum absorption from Ca I 4p 4p3 P0 bound-free opacity and not collisional de-excitation as is the case for warmer K giants.

Luttermoser, Donald G.↗
NASA NTRS1994

On the maximum energy release in flux-rope models of eruptive flares

We determine the photospheric boundary conditions which maximize the magnetic energy released by a loss of ideal-magnetohydrodynamic (MHD) equilibrium in two-dimensional flux-rope models. In these models a loss of equilibrium causes a transition of the flux rope to a lower magnetic energy state at a higher altitude. During the transition a vertical current sheet forms below the flux rope, and reconnection in this current sheet releases additional energy. Here we compute how much energy is released by the loss of equilibrium relative to the total energy release. When the flux-rope radius is small compared to its height, it is possible to obtain general solutions of the Grad-Shafranov equation for a wide range of boundary conditions. Variational principles can then be used to find the particular boundary condition which maximizes the magnetic energy released for a given class of conditions. We apply this procedure to a class of models known as cusp-type catastrophes, and we find that the maximum energy released by the loss of equilibrium is 20.8% of the total energy release for any model in this class. If the additional restriction is imposed that the photospheric magnetic field forms a simple arcade in the absence of coronal currents, then the maximum energy release reduces to 8.6%

Forbes, T. G.↗
NASA NTRS2018

Temperature Structure in the Perseus Cluster Core Observed with Hitomi

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8-20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region _ 7′ in diameter. The SXS was operated with an energy resolution of _5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components

Aharonian, Felix↗
DOE OSTI2020

A non-cooperative meta-modeling game for automated third-party calibrating, validating and falsifying constitutive laws with parallelized adversarial attacks

The evaluation of constitutive models, especially for high-risk and high-regret engineering applications, requires efficient and rigorous third-party calibration, validation and falsification. While there are numerous efforts to develop paradigms and standard procedures to validate models, difficulties may arise due to the sequential, manual, and often biased nature of the commonly adopted calibration and validation processes, thus slowing down data collections, hampering the progress towards discovering new physics, increasing expenses and possibly leading to misinterpretations of the credibility and application ranges of proposed models. This work attempts to introduce concepts from game theory and machine learning techniques to overcome many of these existing difficulties. Here, we introduce an automated meta-modeling game where two competing AI agents systematically generate experimental data to calibrate a given constitutive model and to explore its weakness such that the experiment design and model robustness can be improved through competitions. The two agents automatically search for the Nash equilibrium of the meta-modeling game in an adversarial reinforcement learning framework without human intervention. In particular, a protagonist agent seeks to find the more effective ways to generate data for model calibrations, while an adversary agent tries to find the most devastating test scenarios that expose the weaknesses of the constitutive model calibrated by the protagonist. By capturing all possible design options of the laboratory experiments into a single decision tree, we recast the design of experiments as a game of combinatorial moves that can be resolved through deep reinforcement learning by the two competing players. Our adversarial framework emulates idealized scientific collaborations and competitions among researchers to achieve a better understanding of the application range of the learned material laws and prevent misinterpretations caused by conventional AI-based third-party validation. Numerical examples are given to demonstrate the wide applicability of the proposed meta-modeling game with adversarial attacks on both human-crafted constitutive models and machine learning models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

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

Analytical model for the density distribution in the Io plasma torus

An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Mei, YI↗
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.↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗