Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density model”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Assimilative modeling of low latitude ionosphere

In this paper we present an observation system simulation experiment for modeling low-latitude ionosphere using a 3-dimensional (3-D) global assimilative ionospheric model (GAIM). The experiment is conducted to test the effectiveness of GAIM with a 4-D variational approach (4DVAR) in estimation of the ExB drift and thermospheric wind in the magnetic meridional planes simultaneously for all longitude or local time sectors. The operational Global Positioning System (GPS) satellites and the ground-based global GPS receiver network of the International GPS Service are used in the experiment as the data assimilation source. 'The optimization of the ionospheric state (electron density) modeling is performed through a nonlinear least-squares minimization process that adjusts the dynamical forces to reduce the difference between the modeled and observed slant total electron content in the entire modeled region. The present experiment for multiple force estimations reinforces our previous assessment made through single driver estimations conducted for the ExB drift only.

data assimilation↗
NASA NTRS1974

Molecular column densities in selected model atmospheres

From an examination of predicted column densities, the following conclusions were drawn: (1) The SiO ought to be visible in carbon stars which were generated from triple alpha burning, but absent from carbon stars generated from the CNO bi-cycle. (2) Variation in the observed relative strengths of TiO and ZrO is indicative of real differences in the ratio Ti/Zr. (3) The TiO/ZrO ratio shows a small variation as C/O and effective temperature is changed. (4) Column density of silicon dicarbide (SiC2) is sensitive to abundance, temperature, and gravity; hence all relationships between the strength of SiC2 and other stellar parameters will show appreciable scatter. There is however, a substantial luminosity effect present in the SiC2 column densities. (5) Unexpectedly, SiC2 is anti-correlated with C2. (6) The presence of SiC2 in a carbon star eliminates the possibility of these stars having temperatures greater than or equal to 3000 K, or being produced through the CNO bi-cycle.

Johnson, H. R.↗
DOE OSTI2023

Limits on Leptonic TeV Emission from the Cygnus Cocoon with Swift-XRT

$\gamma$-ray observations of the Cygnus Cocoon, an extended source surrounding the Cygnus X star-forming region, suggest the presence of a cosmic ray accelerator reaching energies up to a few PeV. The very-high-energy (VHE; 0.1-100~TeV) $\gamma$-ray emission may be explained by the interaction of cosmic-ray hadrons with matter inside the Cocoon, but an origin of inverse Compton radiation by relativistic electrons cannot be ruled out. Inverse Compton $\gamma$-rays at VHE are accompanied by synchrotron radiation peaked in X-rays. Hence, X-ray observations may probe the electron population and magnetic field of the source. We observed eleven fields in or near the Cygnus Cocoon with the Neil Gehrels Swift Observatory's X-Ray Telescope (Swift-XRT) totaling 110 ksec. We fit the fields to a Galactic and extra-galactic background model and performed a log-likelihood ratio test for an additional diffuse component. We found no significant additional emission and established upper limits in each field. By assuming that the X-ray intensity traces the TeV intensity and follows an dN/dE ∝ E -2.5 spectrum, we obtained a 90% upper limit of F X < 8.7 x 10 -11 erg cm -2 s -1 or 5.2 x10 -11 erg cm -2 s -1 on the X-ray flux of the entire Cygnus Cocoon between 2 and 10 keV depending on the choice of hydrogen column density model. This suggests that no more than one quarter of the $\gamma$-ray flux at 1 TeV is produced by inverse Compton scattering, when assuming an equipartition magnetic field of ~20 μG.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Simulation and Comparison of Martian Surface Ionization Radiation

The spectrum of energetic particle radiation and corresponding doses at the surface of Mars is being characterized by the Radiation Assessment Detector (RAD), one of ten science instruments on the Mars Science Laboratory (MSL) Curiosity Rover. The time series of dose rate for the first 300 Sols after landing on Mars on August 6, 2012 is presented here. For the comparison to RAD measurements of dose rate, Martian surface ionization radiation is simulated by utilizing observed space quantities. The GCR primary radiation spectrum is calculated by using the Badhwar-O'Neill 2011 (BO11) galactic cosmic ray (GCR) model, which has been developed by utilizing all balloon and satellite GCR measurements since 1955 and the newer 1997-2012 Advanced Composition Explorer (ACE) measurements. In the BO11 model, solar modulation of the GCR primary radiation spectrum is described in terms of the international smoothed sunspot number and a time delay function. For the transport of the impingent GCR primary radiation through Mars atmosphere, a vertical distribution of atmospheric thickness at each elevation is calculated using the vertical profiles of atmospheric temperature and pressure made by Mars Global Surveyor measurements. At Gale Crater in the southern hemisphere, the seasonal variation of atmospheric thickness is accounted for the daily atmospheric pressure measurements of the MSL Rover Environmental Monitoring Station (REMS) by using low- and high-density models for cool- and warm-season, respectively. The spherically distributed atmospheric distance is traced along the slant path, and the resultant directional shielding by Martian atmosphere is coupled with Curiosity vehicle for dose estimates. We present predictions of dose rate and comparison to the RAD measurements. The simulation agrees to within +/- 20% with the RAD measurements showing clearly the variation of dose rate by heliospheric conditions, and presenting the sensitivity of dose rate by atmospheric pressure, which has been found from the RAD experiments and driven by thermal tides on Martian surface.

Kim, Myung-Hee Y.↗
NASA NTRS2005

Combined Radar and Radiometer Analysis of Precipitation Profiles for a Parametric Retrieval Algorithm

A methodology to analyze precipitation profiles using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) and precipitation radar (PR) is proposed. Rainfall profiles are retrieved from PR measurements, defined as the best-fit solution selected from precalculated profiles by cloud-resolving models (CRMs), under explicitly defined assumptions of drop size distribution (DSD) and ice hydrometeor models. The PR path-integrated attenuation (PIA), where available, is further used to adjust DSD in a manner that is similar to the PR operational algorithm. Combined with the TMI-retrieved nonraining geophysical parameters, the three-dimensional structure of the geophysical parameters is obtained across the satellite-observed domains. Microwave brightness temperatures are then computed for a comparison with TMI observations to examine if the radar-retrieved rainfall is consistent in the radiometric measurement space. The inconsistency in microwave brightness temperatures is reduced by iterating the retrieval procedure with updated assumptions of the DSD and ice-density models. The proposed methodology is expected to refine the a priori rain profile database and error models for use by parametric passive microwave algorithms, aimed at the Global Precipitation Measurement (GPM) mission, as well as a future TRMM algorithms.

Masunaga, Hirohiko↗
NASA NTRS2001

Analytical Models for Variable Density Multilayer Insulation Used in Cryogenic Storage

A unique multilayer insulation concept for orbital cryogenic storage was experimentally evaluated at NASA Marshall Space Flight Center (MSFC) using the Multipurpose Hydrogen Test Bed (MHTB). A combination of foam/Multi layer Insulation (MLI) was used. The MLI (45 layers of Double Aluminized Mylar (DAM) with Dacron net spacers) was designed for an on-orbit storage period of 45 days and included several unique features such as: a variable layer density and larger but fewer DAM perforations for venting during ascent to orbit. The focus of this paper is on analytical modeling of the variable density MLI performance during orbital coast periods. The foam/MLI combination model is considered to have five segments. The first segment represents the foam layer. The second, third, and fourth segments represent the three layers of MLI with different layer densities and number of shields. Finally, the last segment is considered to be a shroud that surrounds the last MLI layer. The hot boundary temperature is allowed to vary from 164 K to 305 K. To simulate MLI performance, two approaches are considered. In the first approach, the variable density MLI is modeled layer by layer while in the second approach, a semi-empirical model is applied. Both models account for thermal radiation between shields, gas conduction, and solid conduction through the separator materials. The heat flux values predicted by each approach are compared for different boundary temperatures and MLI systems with 30, 45, 60, and 75 layers.

Hedayat, A.↗
NASA NTRS2023

Structured Type III Radio Bursts Observed in Interplanetary Space

Context . The last few decades have seen numerous studies dedicated to fine structures of type III radio bursts observed in the meter–decameter wavelengths. Most of the explanations of the structured radio emission involve the propagation of electron beams through the strongly inhomogeneous plasma in the low corona. To date, only a few type III bursts with fine structures, observed at hecto-kilometric wavelengths, have been reported. Aims . We report here the existence of numerous structured type III radio bursts observed during the STEREO era by all three WAVES instruments on board STEREO A, B, and Wind. The aim of the study is to report and classify structured type III bursts, and to present the characteristics of their fine structures. The final goal is to try to understand the physical mechanism responsible for the generation of structured radio emission. Methods . In this study we used data from all available spacecraft, specifically STEREO and Wind. We employed 1D density models to obtain the speed of the source of type III radio emission, the electron beam. We also performed a spectral analysis of the fine structures in order to compare their characteristics with the metric–decametric fine structures. Results . The presented similarities of the type III fine structures in the metric to decametric and interplanetary wavelengths indicate that the physical processes responsible for the generation of structured type III radio bursts could be the same, at heights from the low corona to the interplanetary range. We show that the observed structuring and intermittent nature of the type III bursts can be explained by the variation in the level of density fluctuations, at different distances from the Sun.

Non-thermal radiation mechanisms↗
DOE OSTI2021

Aqueous-phase effects on ethanol decomposition over Ru-based catalysts

The effects of an aqueous phase on ethanol decomposition for hydrogen production over a Ru(0001) catalyst surface model have been investigated from first principles. Solvent effects on the reaction mechanism and kinetic parameters have been quantified with the help of a microkinetic reactor model, density functional theory, and an implicit solvation scheme (iSMS). Our calculations indicate that in both vapor- and aqueous-phase reaction environments, the ethanol decomposition starts with acetaldehyde formation on the surface, some of which further dehydrogenates to ketenyl species (CHCO), where the C–C bond cleaves to form methylidyne (CH) and CO. In the vapor phase, adsorbed CH gets hydrogenated to methane, and CO desorbs or undergoes methanation reducing the amount of hydrogen produced. In contrast, under aqueous phase reaction conditions, the methanation is inhibited, and the water–gas shift (WGS) reaction is accelerated, leading to complete conversion of CO to CO 2 and H 2 . Furthermore, calculations indicate that the observed reaction behavior under aqueous phase reforming conditions originates primarily from the higher water chemical potential, and implicit solvent models predict only a small solvation effect.

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

Venus gravity field - Pioneer Venus Orbiter navigation results

The gravity field of Venus has been modeled by a spherical harmonic expansion of the potential to degree and order seven. The estimates of these coeficients were obtained by combining information from 43 short arcs (4 hr) of line-of-sight Doppler data centered at periapsis. The data arcs were distributed in longitude and time over more than two circulations of Venus by the Pioneer Venus Orbiter subperiapsis point which was confined to the band of latitudes from 14 deg N to 17 deg N. Convergence of the solution has been assured by iterating upon the initial estimate. All estimates were performed with zero a priori information on the gravity coefficients. Since the altitude of periapsis for most of the orbits was within the sensible Venusian atmosphere, drag effects on the estimated harmonics have been removed using an exponential atmosphere density model. Estimates of the mass parameter (GM) of Venus using this dataset are also evaluated.

Williams, B. G.↗
NASA NTRS2022

Terrestrial Exospheric Dayside H-density Profile at 3 - 15 RE from UVIS/HDAC and TWINS Lyman-alpha Data Combined

Terrestrial ecliptic dayside observations of the exospheric Lyman-α column intensity between 3–15 Earth radii (RE) by UVIS/HDAC (UVIS – ultraviolet imaging spectrograph; HDAC – hydrogen-deuterium absorption cell) Lyman-α photometer at CASSINI have been analyzed to derive the neutral exospheric H-density profile at the Earth’s ecliptic dayside in this radial range. The data were measured during CASSINI’s swing-by maneuver at the Earth on 18 August 1999 and are published by Werner et al. (2004). In this study the dayside HDAC Lyman-α observations published by Werner et al. (2004) are compared to calculated Lyman-α intensities based on the 3D H-density model derived from TWINS (Two Wide-angle Imaging Neutral-atom Spectrometers) Lyman-α observations between 2008–2010 (Zoennchen et al., 2015). It was found that both Lyman-α profiles show a very similar radial dependence in particular between 3–8RE. Between 3.0–5.5RE impact distance Lyman-α observations of both TWINS and UVIS/HDAC exist at the ecliptic dayside. In this overlapping region the cross-calibration of the HDAC profile against the calculated TWINS profile was done, assuming that the exosphere there was similar for both due to comparable space weather conditions. As a result of the cross-calibration the conversion factor between counts per second and rayleigh, fc = 3.285 counts s−1R−1, is determined for these HDAC observations. Using this factor, the radial H-density profile for the Earth’s ecliptic dayside was derived from the UVIS/HDAC observations, which constrained the neutral H density there at 10RE to a value of 35 cm-3. Furthermore, a faster radial H-density decrease was found at distances above 8RE (≈r−2.37) compared to the lower distances of 3–7RE (≈r−2.37). This increased loss of neutral H above 8RE might indicate a higher rate of H ionization in the vicinity of the magnetopause at 9–11RE (near subsolar point) and beyond, because of increasing charge exchange interactions of exospheric H atoms with solar wind ions outside the magnetosphere.

Jochen H Zoennchen↗
NASA NTRS1985

Shuttle structural dynamics characteristics: The analysis and verification

The space shuttle introduced a new dimension in the complexity of the structural dynamics of a space vehicle. The four-body configuration exhibited structural frequencies as low as 2 hertz with a model density on the order of 10 modes per hertz. In the verification process, certain mode shapes and frequencies were identified by the users as more important than others and, as such, the test objectives were oriented toward experimentally extracting those modes and frequencies for analysis and test correlation purposes. To provide the necessary experimental data, a series of ground vibration tests (GVT's) was conducted using test articles ranging from the 1/4-scale structural replica of the space shuttle to the full-scale vehicle. The vibration test and analysis program revealed that the mode shapes and frequency correlations below 10 hertz were good. The quality of correlation of modes between 10 and 20 hertz ranged from good to fair and that of modes above 20 hertz ranged from poor to good. Since the most important modes, based on user preference, were below 10 hertz, it was judged that the shuttle structural dynamic models were adequate for flight certifications.

Modlin, C. T., Jr.↗
DOE OSTI2019

The role of excited-state ionization on the effective ionization of neutral neon using non-perturbative atomic data

With the aim of providing accurate atomic data for plasma edge temperature and density modeling with neutral neon, new non-perturbative R-Matrix with Pseudostates (RMPS) calculations are performed for both the electron-impact excitation and ionization of neutral neon. The new calculations agree well with past measurements and calculations for both electron-impact excitation and ionization. For ionization, the addition of pseudostates from $5\bar{d}$ to $12\bar{i}$ allows us to extract new cross sections and rate coefficients for all of the excited states through 2p 5 5p. The new RMPS ionization calculation gives a 20% increase in the cross section for the $2{p}^{5}5s$ ionization compared with previous non-perturbative calculations and shows classically-expected n 4 -scaled convergence for the np and nd valence shells. Ionization from these excited states contributes significantly to the effective ionization rate coefficient: a factor of three more than the ground-state ionization rate coefficient at an electron temperature of 25 eV and density of 1014 cm -3 . An accurate, non-perturbative set of atomic rates, both LS- and LSJ-resolved, for neutral neon across a wide range of temperatures and densities is compiled.

Optics↗
NASA NTRSDate not supplied

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗
NASA NTRSDate not supplied

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗
DOE OSTI2018

Charged-cell periodic DFT simulations via an impurity model based on density embedding: Application to the ionization potential of liquid water

Calculations of charged systems in periodic boundary conditions (PBC) are problematic because there are spurious interactions between the charges in different periodic images that can affect the physical picture. In addition, the intuitive limit of Coulomb interactions decaying to zero as the interacting charges are placed at infinite separation no longer applies, and for example total energies become undefined. Leveraging subsystem density functional theory (also known as density embedding) we define an impurity model that embeds a finite neutral or charged subsystem within an extended (infinite) surrounding subsystem. The combination of the impurity model and a consistent choice of the Coulomb reference provides us with an algorithm for evaluating the ionization potential (IP) in extended systems. Additionally, we demonstrate our approach in a pilot calculation of the IP of liquid water, based on a configuration from a prior ab initio molecular dynamics (AIMD) simulation of liquid water (Genova et al., J. Chem. Phys. 2016, 144, 234105). The calculations with the impurity model capture the broadening on the ionization energies introduced by the interactions between the water molecules. Furthermore, the calculated average IP value (10.5 eV) compare favorably to experiments (9.9-10.06 eV) and very recent simulations based on the GW approximation (10.55 eV), while at the same time outperforming density embedding calculations carried out with a naïve handling of the electrostatic interactions (about 7 eV).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Defect Thermodynamics and Transport Properties of Proton Conducting Perovskite Electrode and Electrolyte Materials Evaluated Based on Density Functional Theory Modeling

Both electron-rich and electron-poor perovskite oxides have been used in solid oxide cell applications as electrode and electrolyte materials. The rich oxygen defect chemistry and its coupling to temperature, hydrogen-steam or oxygen-steam gas pressure, or to the applied potentials creates enormous complexities for modeling performance and degradation of the materials. Herein, density functional theory-based thermodynamic modeling was carried out to describe the defect chemistry and transport properties of the proton-conducting electrolyte BaZr1-xYxO3-δ (x≤0.1) and of the triple-conducting perovskite (La,Ba)(Fe,M)O3-δ (M=Y and Zr). The defect thermodynamics of intrinsic point defects and the hydrogen-related defect reactions were solved in integrated defect models and further used to predict the Brouwer diagram and the transport properties of the functional perovskites. For the electron-poor electrolytes BaZr0.9Y0.1O3-δ, the developed model has been used to describe the experimental transport properties in the SOC operating conditions. Specifically, the roles played by the acceptor-bound holes and the intrinsic and hydrogen point defects upon the conductivities of holes, protons, and oxygen vacancies under the hydrogen-rich and oxygen-rich conditions at various humidity levels were demonstrated. A defect modeling tool was also developed for the triple-conducting perovskite (La,Ba)(Fe,M)O3-δ (M=Y and Zr) to examine magnetic effects and hydride defects in defect equilibria.

defect thermodynamics↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗