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 379 records · Page 21

The Prospective Search for Highly Ionized Technetium in Hot (pre-) White Dwarfs

The discovery of technetium (Tc) in the atmospheres of red giants by Merrill (1952, ApJ, 116, 21) constituted convincing proof that s-process nucleosynthesis is indeed occurring in evolved stars. In principle, Tc should still be present in the atmospheres of hot post-AGB stars and (pre-) white dwarfs although, due to radioactive decay, it should be present in decreasing quantities along post-AGB evolution. The recent discovery of a large number of trans-iron group elements in hot white dwarfs with atomic numbers in the range A = 30-56 (Zn to Ba) raises the prospect that Tc (A = 43) may also be detected. However, this is currently not feasible because no atomic data exist for ionization stages beyond Tc ii. As an initial step, we calculated atomic energy levels and oscillator strengths of Tc iv-vi and used these data to compute non-local thermodynamic equilibrium (NLTE) model atmospheres to estimate at which minimum abundance level Tc could be detected. We show that Tc lines can be found in ultraviolet spectra of hot white dwarfs provided Tc is as abundant as other detected trans-Fe elements. We find that radiative levitation can keep Tc in large, easily detectable quantities in the atmosphere. A direct identification of Tc lines is still not feasible because wavelength positions cannot be computed with necessary precision. Laboratory measurements are necessary to overcome this problem. Our results suggest that such efforts are beneficial to the astrophysical community.

prospective↗

Stellar Laboratories

For the spectral analysis of high-resolution and high signal-to-noise (SN) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. Aims: To identify molybdenum lines in the ultraviolet (UV) spectra of the DA-type white dwarf G191B2B and the DO-type white dwarf RE 0503289 and, to determine their photospheric Mo abundances, reliable Mo iv-vii oscillator strengths are used. Methods: We newly calculated Mo iv-vii oscillator strengths to consider their radiative and collisional bound-bound transitions indetail in our NLTE stellar-atmosphere models for the analysis of Mo lines exhibited in high-resolution and high SN UV observations of RE 0503289.Results. We identified 12 Mo v and nine Mo vi lines in the UV spectrum of RE 0503289 and measured a photospheric Mo abundance of 1.2 3.0 104(mass fraction, 22 500 56 400 times the solar abundance). In addition, from the As v and Sn iv resonance lines,we measured mass fractions of arsenic (0.51.3 105, about 300 1200 times solar) and tin (1.33.2 104, about 14 300 35 200 times solar). For G191B2B, upper limits were determined for the abundances of Mo (5.3 107, 100 times solar) and, in addition, for Kr (1.1106, 10 times solar) and Xe (1.7107, 10 times solar). The arsenic abundance was determined (2.35.9 107, about 21 53 times solar). A new, registered German Astrophysical Virtual Observatory (GAVO) service, TOSS, has been constructed to provide weighted oscillator strengths and transition probabilities.Conclusions. Reliable measurements and calculations of atomic data are a prerequisite for stellar-atmosphere modeling. Observed Mo v-vi line profiles in the UV spectrum of the white dwarf RE 0503289 were well reproduced with our newly calculated oscillator strengths. For the first time, this allowed the photospheric Mo abundance in a white dwarf to be determined.

Rauch, T.↗

Vii. New Kr IV - VII Oscillator Strengths and an Improved Spectral Analysis of the Hot, Hydrogen-deficient Do-type White Dwarf RE 0503-289

For the spectral analysis of high-resolution and high signal-to-noise (SN) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. Aims. New Krivvii oscillator strengths for a large number of lines enable us to construct more detailed model atoms for our NLTEmodel-atmosphere calculations. This enables us to search for additional Kr lines in observed spectra and to improve Kr abundance determinations. Methods. We calculated Krivvii oscillator strengths to consider radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of Kr lines that are exhibited in high-resolution and high SN ultraviolet (UV)observations of the hot white dwarf RE 0503.

Rauch, T.↗

Regime-Based Evaluation of Cloudiness in CMIP5 Models

The concept of Cloud Regimes (CRs) is used to develop a framework for evaluating the cloudiness of 12 fifth Coupled Model Intercomparison Project (CMIP5) models. Reference CRs come from existing global International Satellite Cloud Climatology Project (ISCCP) weather states. The evaluation is made possible by the implementation in several CMIP5 models of the ISCCP simulator generating for each gridcell daily joint histograms of cloud optical thickness and cloud top pressure. Model performance is assessed with several metrics such as CR global cloud fraction (CF), CR relative frequency of occurrence (RFO), their product (long-term average total cloud amount [TCA]), cross-correlations of CR RFO maps, and a metric of resemblance between model and ISCCP CRs. In terms of CR global RFO, arguably the most fundamental metric, the models perform unsatisfactorily overall, except for CRs representing thick storm clouds. Because model CR CF is internally constrained by our method, RFO discrepancies yield also substantial TCA errors. Our findings support previous studies showing that CMIP5 models underestimate cloudiness. The multi-model mean performs well in matching observed RFO maps for many CRs, but is not the best for this or other metrics. When overall performance across all CRs is assessed, some models, despite their shortcomings, apparently outperform Moderate Resolution Imaging Spectroradiometer (MODIS) cloud observations evaluated against ISCCP as if they were another model output. Lastly, cloud simulation performance is contrasted with each model's equilibrium climate sensitivity (ECS) in order to gain insight on whether good cloud simulation pairs with particular values of this parameter.

cloud climatology↗

VIII. New Zr IV-VII, Xe IV-V, and Xe VII Oscillator Strengths and the Al, Zr, and Xe Abundances in the Hot White Dwarfs G191-B2B and RE0503-289

Context: For the spectral analysis of high-resolution and high-signal-to-noise spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. Aims: To search for zirconium and xenon lines in the ultraviolet (UV) spectra of G191−B2B and RE 0503−289, new Zr iv-vii, Xe iv-v, and Xe vii oscillator strengths were calculated. This allows, for the first time, determination of the Zr abundance in white dwarf (WD) stars and improvement of the Xe abundance determinations.Methods: We calculated Zr iv-vii, Xe iv-v, and Xe vii oscillator strengths to consider radiative and collisional bound-bound transitions of Zr and Xe in our NLTE stellar-atmosphere models for the analysis of their lines exhibited in UV observations of the hot WDs G191−B2B and RE 0503−289. Results: We identified one new Zr iv, 14 new Zr v, and ten new Zr vi lines in the spectrum of RE 0503−289. Zr was detected for the first time in a WD. We measured a Zr abundance of −3.5 +/- 0.2 (logarithmic mass fraction, approx. 11 500 times solar). We identified five new Xe vi lines and determined a Xe abundance of −3.9 +/- 0.2 (approx. 7500 times solar). We determined a preliminary photospheric Al abundance of −4.3 +/- 0.2 (solar) in RE 0503−289. In the spectra of G191−B2B, no Zr line was identified. The strongest Zr iv line (1598.948 Å) in our model gave an upper limit of −5.6 +/- 0.3 (approx. 100 times solar). No Xe line was identified in the UV spectrum of G191−B2B and we confirmed the previously determined upper limit of −6.8 +/- 0.3 (ten times solar). Conclusions: Precise measurements and calculations of atomic data are a prerequisite for advanced NLTE stellar-atmosphere modeling. Observed Zr iv-vi and Xe vi-vii line profiles in the UV spectrum of RE 0503−289 were simultaneously well reproduced with our newly calculated oscillator strengths.

Rauch, T.↗

Mercury-T: A New Code to Study Tidally Evolving Multi-Planet Systems. Applications to Kepler-62

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass 0:1 M(sub ⨀), Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water..

Dynamical evolution and stability↗

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC↗

A Library of ATMO Forward Model Transmission Spectra for Hot Jupiter Exoplanets

We present a grid of forward model transmission spectra, adopting an isothermal temperature-pressure profile, alongside corresponding equilibrium chemical abundances for 117 observationally significant hot exoplanets (equilibrium temperatures of 547-2710 K). This model grid has been developed using a 1D radiative-convective-chemical equilibrium model termed ATMO, with up-to-date high-temperature opacities. We present an interpretation of observations of 10 exoplanets, including best-fitting parameters and X(exp 2) maps. In agreement with previous works, we find a continuum from clear to hazy/cloudy atmospheres for this sample of hot Jupiters. The data for all the 10 planets are consistent with subsolar to solar C/O ratio, 0.005 to 10 times solar metallicity and water rather than methane-dominated infrared spectra. We then explore the range of simulated atmospheric spectra for different exoplanets, based on characteristics such as temperature, metallicity, C/O ratio, haziness and cloudiness. We find a transition value for the metallicity between 10 and 50 times solar, which leads to substantial changes in the transmission spectra. We also find a transition value of C/O ratio, from water to carbon species dominated infrared spectra, as found by previous works, revealing a temperature dependence of this transition point ranging from approximately 0.56 to approximately 1-1.3 for equilibrium temperatures from approximately 900 to approximately 2600 K. We highlight the potential of the spectral features of HCN and C2H2 to constrain the metallicities and C/O ratios of planets, using James Webb Space Telescope (JWST) observations. Finally, our entire grid (approximately 460 000 simulations) is publicly available and can be used directly with the JWST simulator PandExo for planning observations.

techniques: spectroscopic – planets and satellit↗

Modeling of the near-electrode regions of arcjets. I - Coupling of the flowfield to the non-equilibrium boundary layers

A model for arc attachment at an anode in contact with a flowing thermal plasma is described. Current limitations are identified which arise when the convective transport of ions into the electrical, thermal and hydrodynamic boundary layers is incapable of supplying the loss of ions through three-body recombination in the volume, and catalytic recombination on the electrode surface. It is shown that the anode current where this condition is achieved can be exceeded provided there is a significant Joule heating in the boundary layer. The saturation current (current density above which the anode must emit ions) is found to substantially increase with increasing free stream velocity and anode temperature.

Cappelli, Mark A.↗

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗

Hydrodynamic Models of Line-Driven Accretion Disk Winds III: Local Ionization Equilibrium

We present time-dependent numerical hydrodynamic models of line-driven accretion disk winds in cataclysmic variable systems and calculate wind mass-loss rates and terminal velocities. The models are 2.5-dimensional, include an energy balance condition with radiative heating and cooling processes, and includes local ionization equilibrium introducing time dependence and spatial dependence on the line radiation force parameters. The radiation field is assumed to originate in an optically thick accretion disk. Wind ion populations are calculated under the assumption that local ionization equilibrium is determined by photoionization and radiative recombination, similar to a photoionized nebula. We find a steady wind flowing from the accretion disk. Radiative heating tends to maintain the temperature in the higher density wind regions near the disk surface, rather than cooling adiabatically. For a disk luminosity L (sub disk) = solar luminosity, white dwarf mass M(sub wd) = 0.6 solar mass, and white dwarf radii R(sub wd) = 0.01 solar radius, we obtain a wind mass-loss rate of M(sub wind) = 4 x 10(exp -12) solar mass yr(exp -1) and a terminal velocity of approximately 3000 km per second. These results confirm the general velocity and density structures found in our earlier constant ionization equilibrium adiabatic CV wind models. Further we establish here 2.5D numerical models that can be extended to QSO/AGN winds where the local ionization equilibrium will play a crucial role in the overall dynamics.

Pereyra, Nicolas Antonio↗

A study of plasmaspheric density distributions for diffusive equilibrium conditions

The plasmaspheric density distribution has been modeled for a range of solar cycle, seasonal and diurnal conditions with a magnetic flux tube dependent diffusive equilibrium model by using experimentally determined values of ionospheric parameters at 675 km as boundary conditions. Data is presented in terms of plasmaspheric H(+) and He(+) density contours, total flux tube content, and equatorial plasma density for a range of L-values from 1.15 to 3.0. The variation of equatorial density with L-value shows good agreement with the 1/L exp 4 dependence observed experimentally. The results show that the model predicts larger solar cycle and diurnal variation in equatorial plasma density than observed using whistler techniques. However, the whistler method requires a model to deduce the equatorial density and is therefore open to interpretation. Seasonal variations are rather artificial since in this general model no attempt has been made to match equatorial densities for flux tubes emanating from the winter and summer hemispheres.

Li, W.↗

The influence of centrifugal forces on the B field structure of an axially symmetric equilibrium magnetosphere

A model is presented of an axially symmetric pole-on magnetosphere in MHD force balance, in which both plasma thermal pressure gradients and centrifugal force are taken into account. Assuming that planetary rotation leads to differentially rotating magnetotail field lines, the deformation of magnetotail field lines under the influence of both thermal plasma pressure and centrifugal forces was calculated. Analytic solutions to the Grad-Shafranov equation are presented, which include the centrifugal force term. It is shown that the nonrotational magnetosphere with hot thermal plasma leads to a field configuration without a toroidal B(phi) component and without field-aligned Birkeland currents. The other extreme, a rapidly rotating magnetosphere with cold plasma, leads to a configuration in which plasma must be confined within a thin disk in a plane where the radial magnetic field component B(r) vanishes locally.

Ye, Gang↗

Comparing the results of an analytical model of the no-vent fill process with no-vent fill test results for a 4.96 cu m (175 cu ft) tank

NASA-Lewis has been investigating a no-vent fill method for refilling cryogenic storage tanks in low gravity. Analytical modeling based on analyzing the heat transfer of a droplet has successfully represented the process in 0.034 and 0.142 sq m commercial dewars using liquid nitrogen and hydrogen. Recently a large tank (4.96 sq m) was tested with hydrogen. This lightweight tank is representative of spacecraft construction. This paper presents efforts to model the large tank test data. The droplet heat transfer model is found to overpredict the tank pressure level when compared to the large tank data. A new model based on equilibrium thermodynamics has been formulated. This new model is compared to the published large scale tank's test results as well as some additional test runs with the same equipment. The results are shown to match the test results within the measurement uncertainty of the test data except for the initial transient wall cooldown where it is conservative (i.e., overpredicts the initial pressure spike found in this time frame).

Taylor, William J.↗

Comparing the results of an analytical model of the no-vent fill process with no-vent fill test results for a 4.96 cubic meters (175 cubic feet) tank

The NASA Lewis Research Center (NASA/LeRC) have been investigating a no-vent fill method for refilling cryogenic storage tanks in low gravity. Analytical modeling based on analyzing the heat transfer of a droplet has successfully represented the process in 0.034 m and 0.142 cubic m commercial dewars using liquid nitrogen and hydrogen. Recently a large tank (4.96 cubic m) was tested with hydrogen. This lightweight tank is representative of spacecraft construction. This paper presents efforts to model the large tank test data. The droplet heat transfer model is found to over predict the tank pressure level when compared to the large tank data. A new model based on equilibrium thermodynamics has been formulated. This new model is compared to the published large scale tank's test results as well as some additional test runs with the same equipment. The results are shown to match the test results within the measurement uncertainty of the test data except for the initial transient wall cooldown where it is conservative (i.e., overpredicts the initial pressure spike found in this time frame).

Taylor, William J.↗

Regional eddy vorticity transport and the equilibrium vorticity budgets of a numerical model ocean circulation

A mean vorticity budget analysis is presented of Holland's (1978) numerical ocean general circulation experiment. The stable budgets are compared with classical circulation theory to emphasize the ways in which the mesoscale motions of the model alter (or leave unaltered) classical vorticity balances. The basinwide meridional transports of vorticity by the mean flow and by the mesoscale flow in the mean are evaluated to establish the role(s) of the mesoscale in the larger scale equilibrium vorticity transports. The vorticity equation for this model fluid system is presented and the budget analysis method is described. Vorticity budgets over the selected regions and on a larger scale are given, and a summary of budget results is provided along with remarks about the utility of this type of analysis.

Harrison, D. E.↗

Modeling Inflation Using a Non-Equilibrium Equation of Exchange

Inflation is a change in the prices of goods that takes place without changes in the actual values of those goods. The Equation of Exchange, formulated clearly in a seminal paper by Irving Fisher in 1911, establishes an equilibrium relationship between the price index P (also known as "inflation"), the economy's aggregate output Q (also known as "the real gross domestic product"), the amount of money available for spending M (also known as "the money supply"), and the rate at which money is reused V (also known as "the velocity of circulation of money"). This paper offers first a qualitative discussion of what can cause these factors to change and how those causes might be controlled, then develops a quantitative model of inflation based on a non-equilibrium version of the Equation of Exchange. Causal relationships are different from equations in that the effects of changes in the causal variables take time to play out-often significant amounts of time. In the model described here, wages track prices, but only after a distributed lag. Prices change whenever the money supply, aggregate output, or the velocity of circulation of money change, but only after a distributed lag. Similarly, the money supply depends on the supplies of domestic and foreign money, which depend on the monetary base and a variety of foreign transactions, respectively. The spreading of delays mitigates the shocks of sudden changes to important inputs, but the most important aspect of this model is that delays, which often have dramatic consequences in dynamic systems, are explicitly incorporated.macroeconomics, inflation, equation of exchange, non-equilibrium, Athena Project

Chamberlain, Robert G.↗