Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient structure”

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 1,135 records · Page 63

Stochastic inflation lattice simulations - Ultra-large scale structure of the universe

Non-Gaussian fluctuations for structure formation may arise in inflation from the nonlinear interaction of long wavelength gravitational and scalar fields. Long wavelength fields have spatial gradients, a (exp -1), small compared to the Hubble radius, and they are described in terms of classical random fields that are fed by short wavelength quantum noise. Lattice Langevin calculations are given for a toy model with a scalar field interacting with an exponential potential where one can obtain exact analytic solutions of the Fokker-Planck equation. For single scalar field models that are consistent with current microwave background fluctuations, the fluctuations are Gaussian. However, for scales much larger than our observable Universe, one expects large metric fluctuations that are non-Gaussian. This example illuminates non-Gaussian models involving multiple scalar fields which are consistent with current microwave background limits.

Salopek, D. S.↗

Oregon transect: Comparison of leaf-level reflectance with canopy-level and modelled reflectance

The Oregon Transect Ecosystem Research (OTTER) project involves the collection of a variety of remotely-sensed and in situ measurements for characterization of forest biophysical and biochemical parameters. The project includes nine study plots located along an environmental gradient in west-central Oregon, extending from the Pacific coast inland approximately 300km. These plots represent a broad range in ecosystem structure and function. Within the OTTER project, the sensitivity of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) signal to absorption by foliar biochemicals is being examined. AVIRIS data were acquired over all plots in conjunction with the four OTTER Multi-sensor Aircraft Campaigns spanning the growing season. Foilage samples were gathered during each campaign for biochemical determination (at Ames Research Center), to estimate stand-level constituency at each plot. Directional-hemispheric leaf reflectance throughout the 400-2400nm region was measured in the laboratory as an aid to interpreting concurrent AVIRIS data. Obtaining leaf spectra in this manner reduces or eliminates the confounding influences of atmosphere, canopy architecture, and reflectance by woody components, understory, and exposed soils which are present in airborne observations. These laboratory spectra were compared to simulated spectra derived by inverting the PROSPECT leaf-level canopy reflectance derived from AVIRIS data by use of the LOWTRAN-7 atmospheric radiative-transfer model.

Johnson, Lee F.↗

Solution of the three-dimensional Helmholtz equation with nonlocal boundary conditions

The Helmholtz equation is solved within a three-dimensional rectangular duct with a nonlocal radiation boundary condition at the duct exit plane. This condition accurately models the acoustic admittance at an arbitrarily-located computational boundary plane. A linear system of equations is constructed with second-order central differences for the Helmholtz operator and second-order backward differences for both local admittance conditions and the gradient term in the nonlocal radiation boundary condition. The resulting matrix equation is large, sparse, and non-Hermitian. The size and structure of the matrix makes direct solution techniques impractical; as a result, a nonstationary iterative technique is used for its solution. The theory behind the nonstationary technique is reviewed, and numerical results are presented for radiation from both a point source and a planar acoustic source. The solutions with the nonlocal boundary conditions are invariant to the location of the computational boundary, and the same nonlocal conditions are valid for all solutions. The nonlocal conditions thus provide a means of minimizing the size of three-dimensional computational domains.

Hodge, Steve L.↗

Ionospheric Instabilities and Shock Capturing Algorithms

We describe the basic physics of the gradient-drift instability and the collisional Rayleigh-Taylor instability, the two plasma instabilities which dominate the production of small-scale structure in the earth's ionosphere. We then focus our attention on the nonlinear evolution of these instabilities for two particular cases: 1) the recursive bifurcation and "freezing" of artificial ion clouds released in the ionosphere; and 2) equatorial spread F. Finally, we discuss the numerical challenges faced in attempting to numerically simulate these instabilities, and the role played by shock-capturing algorithms therein.

Zalesak, Steven T.↗

Analysis of the Radar Reflectivity of Aircraft Vortex Wakes

Radar has been proposed as a way to track wake vortices to reduce aircraft spacing and tests have revealed radar echoes from aircraft wakes in clear air. The results are always interpreted qualitatively using Tatarski's theory of weak scattering by isotropic atmospheric turbulence. The goal of the present work was to predict the value of the radar cross-section (RCS) using simpler models. This is accomplished in two steps. First, the refractive index is obtained. Since the structure of the aircraft wakes is different from atmospheric turbulence, three simple mechanisms specific to vortex wakes are considered: (1) Radial density gradient in a two-dimensional vortex, (2) three-dimensional fluctuations in the vortex cores, and (3) Adiabatic transport of the atmospheric fluid in a two-dimensional oval surrounding the pair of vortices. The index of refraction is obtained more precisely for the two-dimensional mechanisms than for the three-dimensional ones. In the second step, knowing the index of refraction, a scattering analysis is performed. Tatarski's weak scattering approximation is kept but the usual assumptions of a far-field and a uniform incident wave are dropped. Neither assumption is generally valid for a wake that is coherent across the radar beam. For analytical insight, a simpler approximation that invokes, in addition to weak scattering, the far-field and wide cylindrical beam assumptions, is also developed and compared with the more general analysis. The predicted RCS values for the oval surround the vortices (mechanism C) agree with the experiments of Bilson conducted over a wide range of frequencies. However, the predictions have a cut-off away from normal incidence which is not present in the measurements. Estimates suggest that this is due to turbulence in the baroclinic vorticity generated at the boundary of the oval. The reflectivity of a vortex itself (mechanism A) is comparable to that of the oval (mechanism C) but cuts-off at frequencies lower than those considered in all the experiments to date. The RCS of a vortex happens to peak at the frequency (about 49 MHz) where atmospheric radars (known as ST radars) operate and so the present prediction could be verified in the future. Finally , we suggest that hot engine exhaust could increase RCE by 40 db and reveal vortex circulation, provided its mixing with the surroundings is prevented in the laminarising flow of the vortices.

Shariff, Karim↗

Ion Milling of Sapphire

The ion figuring system at the Marshall Space Flight Center has been successfully used for at least three previous investigations into the ion milling of metals. The research was directed toward improving the surface quality of X-ray directing optics. These studies were performed on surfaces that were already hand polished to an excellent surface quality and were intended to remove the residual unwanted figure left by those techniques. The ion milling was typically carried out on test surfaces or mandrels that were several centimeters in width and length. The good thermal conductivity of the metal samples allowed the ion beam to be directed onto the sample for an indefinite period of time. This is not true of sapphire or most electrical insulators and problems have arisen in recent attempts to ion mill thin samples of sapphire. The failure and fracture of the material was likely due to thermal stresses and the relatively low thermal conductivity of sapphire (compared to most metals), These assumed stresses actually provided the key as to how they might be monitored. A thermal gradient in the sapphire sample will induce an effective index of refraction change and because of the shape constraint and the crystal structure and simple thermal expansion, this index change will be nonuniform across the sample. In all but simple cubic crystal structures, this leads to a spatially nonuniform optical retardance induced on any polarized optical beam traversing the sample, and it is this retardance that can be monitored using standard polarimetric procedures.

Don A. Gregory↗

In-Situ Measurements of HCN and CH3CN in the Pacific Troposphere: Sources, Sinks, and Comparisons with Spectroscopic Observations

We report the first in-situ measurements of hydrogen cyanide (HCN) and acetonitrile (CH3CN) from the Pacific troposphere (0-12 km) obtained during the NASA/Trace-P mission (Feb.-April, 2001). Mean HCN and CH3CN mixing ratios of 243 (+/-118) ppt and 149 (+/-56) ppt respectively, were measured. The in-situ observations correspond to a total HCN column of 4.4-4.9 x 10(exp 15) molec. cm(exp -2) and a CH3CN column of 2.8-3.0 x 10(exp 15) molec. cm(exp -2). This HCN column is in good agreement with available spectroscopic observations. The atmospheric concentrations of HCN and CH3CN were greatly influenced by outflow of pollution from Asia. There is a linear relationship between the mixing ratios of HCN and CH3CN, and in turn these are well correlated with tracers of biomass combustion (e.g. CH3Cl, CO). Relative enhancements with respect to known tracers of biomass combustion within selected plumes in the free troposphere, and pollution episodes in the boundary layer allow an estimation of a global biomass burning source of 0.8+/-0.4 Tg (N)/y for HCN and 0.4+/-0.1 Tg (N)/y for CH3CN. In comparison, emissions from automobiles and industry are quite small (<0.05 Tg (N)/y). The vertical structure of HCN and CH3CN indicated reduced mixing ratios in the MBL (Marine Boundary Layer). Using, a simple box model, the observed gradients across the top of the MBL are used to derive an oceanic flux of 6.7 x 10(exp -15) g (N) cm(exp -2)/s for HCN and 4.8 x 10(exp -15) g (N) cm(exp -2)/s for CH3CN. An air-sea exchange model is used to conclude that this flux can be maintained if the oceans are under-saturated in HCN and CH3CN by 23% and 17%, respectively. It is inferred that oceanic loss is a dominant sink for these nitrites, and they deposit some 1.3 Tg (N) of nitrogen annually to the oceans. Assuming reaction with OH radicals and loss to the oceans as the major removal processes, a mean atmospheric residence time of 4.7 months for HCN and 5.1 months for CH3CN is calculated. A global budget analysis shows that the sources and sinks of HCN and CH3CN are roughly in balance. There are indications that biogenic sources may also be present. Mechanisms involved in nitrate formation during combustion and removal in the oceans are poorly understood.

Singh, Hanwant B.↗

African Easterly Jet: Structure and Maintenance

This article investigates the African Easterly Jet (AEJ), its structure and the forcings contributing to its maintenance, critically revisiting previous work which attributed the maintenance of the jet to soil moisture gradients over tropical Africa. A state-of-the-art global model in a high-end computer framework is used to produce a 3-member 73-year ensemble run forced by observed SST to represent the Control run. The AEJ as produced by the Control is compared with the representation of the AEJ in the European Center for Medium Range Forecast Reanalyses (ERA-40) and other observational data sets and found very realistic. Five Experiments are then performed, each represented by sets of 3-member 22 year long (1980-2001) ensemble runs. The goal of the Experiments is to investigate the role of meridional soil moisture gradients, different land surface properties and orography. Unlike previous studies, which have suppressed soil moisture gradients within a highly idealized framework (i.e., the so-called bucket model), terrestrial evaporation control is here achieved with a highly sophisticated landsurface treatment and with an extensively tested and complex methodology. The results show that the AEJ is suppressed by a combination of absence of meridional evaporation gradients over Africa and constant vegetation, even if the individual forcings taken separately do not lead to the AEJ disappearance, but only its modification. Moreover, the suppression of orography also leads to a different circulation in which there is no AEJ. This work suggests that it is not just soil moisture gradients, but a unique combination of geographical features present only in northern tropical Africa, which causes and maintains the jet.

Wu, Man-Li C.↗

HD139614: the Interferometric Case for a Group-Ib Pre-Transitional Young Disk

The Herbig Ae star HD139614 is a group-Ib object, which featureless SED indicates disk flaring and a possible pre-transitional evolutionary stage. We present mid- and near-IR interferometric results collected with MIDI, AMBER and PIONIER with the aim of constraining the spatial structure of the 0.1-10 AU disk region and assess its possible multi-component structure. A two-component disk model composed of an optically thin 2-AU wide inner disk and an outer temperature-gradient disk starting at 5.6 AU reproduces well the observations. This is an additional argument to the idea that group-I HAeBe inner disks could be already in the disk-clearing transient stage. HD139614 will become a prime target for mid-IR interferometric imaging with the second-generation instrument MATISSE of the VLTI.

Stellar Interferometry↗

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the shock cell structure and noise of supersonic jets

A linear solution modeling the shock cell structure of an axisymmetric supersonic jet operated at off-design conditions is developed by the method of multiple-scales. The model solution takes into account the gradual spatial change of the mean flow in the downstream direction. Turbulence in the mixing layer of the jet has the tendency of smoothing out the sharp velocity and density gradients induced by the shocks. To simulate this effect, eddy viscosity terms are incorporated in the model. It is known that the interaction between the quasi-periodic shock cells and the downstream propagating large turbulence structures in the mixing layer of the jet is responsible for the generation of broadband shock associated noise. Experimentally, the dominant part of this noise has been found to originate from the part of the jet near the end of the potential core. Calculated shock cell spacing at the end of the jet core according to the present model is used to estimate the peak frequencies of the shock associated noise for a range of observation angles. Very favorable agreement with experimental measurements is found.

Tam, C. K. W.↗

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Simple atmospheric perturbation models for sonic-boom-signature distortion studies

Sonic-boom propagation from flight level to ground is influenced by wind and speed-of-sound variations resulting from temperature changes in both the mean atmospheric structure and small-scale perturbations. Meteorological behavior generally produces complex combinations of atmospheric perturbations in the form of turbulence, wind shears, up- and down-drafts and various wave behaviors. Differences between the speed of sound at the ground and at flight level will influence the threshold flight Mach number for which the sonic boom first reaches the ground as well as the width of the resulting sonic-boom carpet. Mean atmospheric temperature and wind structure as a function of altitude vary with location and time of year. These average properties of the atmosphere are well-documented and have been used in many sonic-boom propagation assessments. In contrast, smaller scale atmospheric perturbations are also known to modulate the shape and amplitude of sonic-boom signatures reaching the ground, but specific perturbation models have not been established for evaluating their effects on sonic-boom propagation. The purpose of this paper is to present simple examples of atmospheric vertical temperature gradients, wind shears, and wave motions that can guide preliminary assessments of nonturbulent atmospheric perturbation effects on sonic-boom propagation to the ground. The use of simple discrete atmospheric perturbation structures can facilitate the interpretation of the resulting sonic-boom propagation anomalies as well as intercomparisons among varied flight conditions and propagation models.

Ehernberger, L. J.↗

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗