Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium reconstruction”

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 19 records

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗

Equilibrium condensation in a solar nebula

In attempts to reconstruct the environment of condensation of solar system materials, particularly exemplified by certain meteorite components, the relative temperatures of the gas and the solid are of critical importance. The relationships that determine the heat balance in a circumsolar grain-gas system are examined. Fundamental considerations show that regardless of opacity or gas density, the gas will always be at a higher temperature than the solid in such regions of the system where condensation is possible. Implications of the characteristic temperature differential between the gas and the condensing solid are discussed.

Arrhenius, G.↗

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗

Using Multiplicative Algebraic Reconstruction Techniques (MART) to Derive Instrument Requirements for Computed Tomography Imaging Spectrographs

Computed Tomography Imaging Spectrographs (CTISs), which are generally slitless or large-aperture spectrographs that observe a wide, dispersed field of view in multiple dispersion angles or diffraction orders, provide a unique opportunity to capture spectral information over an extended source, such as the Sun, but require software reconstruction techniques to be fully utilized. A useful data product that can be recovered from CTIS observations is intensity maps of solar features in single spectral lines. These intensity maps can then be used to determine temperature, density, abundance, and equilibrium properties of the emitting plasma. We apply the multiplicative algebraic reconstruction technique (MART) to example data from a variety of CTIS configurations to determine the capability and limitations of the method to return spectrally-pure maps. By completing this study, we aim to establish a path to derive requirements for CTIS instruments that rely on reconstruction techniques to meet their science objectives.

Amy Winebarger↗

Spinel cataclasites in 15445 and 72435 - Petrology and criteria for equilibrium

The problem of establishing the existence of equilibrium among the coexisting phases in the rock is addressed by presenting petrographic and mineral chemistry data on a new spinel cataclasite from 15445 (clast H) and data more extensive than those previously available on two clasts in 72435. Criteria useful in reconstructing the original petrology of these and other spinel cataclasites are analyzed by considering equilibrium among the different phases, that is, the mono- or polymict nature of these cataclasized samples. Finally, the role of impact processes in disturbing the equilibria is discussed.

Baker, M. B.↗

3D Material Response Analysis of PICA Pyrolysis Experiments

Primarily interested in improving ablation modeling for use in inverse reconstruction of flight environments on ablative heat shields. Ablation model is essentially a component of the heat flux sensor, so model uncertainties lead to measurement uncertainties. Non-equilibrium processes have been known to be significant in low density ablators for a long time, but increased accuracy requirements of the reconstruction process necessitates incorporating this physical effect. Attempting to develop a pyrolysis model for implementation in material response based on the PICA data produced by Bessire and Minton. Pyrolysis gas species molar yields as a function of temperature and heating rate. Several problems encountered while trying to fit Arrhenius models to the data led to further investigation of the experimental setup.

Oliver, Brandon A.↗

The Dos and Don'ts of how to Build a Planet, Using the Moon as an Example

The bulk chemical compositions of planets may yield important clues concerning planetary origins. Failing that, bulk compositions are still important, in that they constrain calculation of planetary mineralogies and also constrain the petrogenesis of basaltic magmas. In the case of the Earth, there is little or no debate about the composition of the Earth's upper mantle. This is because our sample collections contain peridotitic xenoliths of that mantle. The most fertile of these are believed to have been little modified from their primary compositions. Using these samples and chondritic meteorites as a starting point, small perturbations on the compositions of existing samples allow useful reconstruction of the bulk silicate Earth (BSE). Elsewhere, I have argued that the next simplest case is the Eucrite Parent Body (EPB). Reconstructions based on Sc partitioning indicate that the EPB can be well approximated by a mixture of 20% eucrite and 80% equilibrium olivine. This leads to a parent body that is similar to CO (or devolatilized CM) chondrites. Partial melting experiments on CM chondrites confirm this model, because the residual solids in these experiments are dominated by olivine with minor pigonite [3]. The most difficult bodies to reconstruct are those that have undergone the most differentiation. Both the Moon and Mars may have passed through a magma ocean stage. In any event, lunar and martian basalts, unlike eucrites, were not derived from undifferentiated source regions. Reconstructions are primarily based on compositional trends within the basalts themselves with some critical assumptions: (i) Refractory lithophile elements (Ca, Al, REE, actinides) are presumed to be in chondritic relative abundances; and (ii) some major element ratio is believed to exist in a chondritic ratio (e.g., Mg/Si, Mg/Al). The most commonly used parameter is Mg/Si.

Jones, J. H.↗

High-Fidelity Simulations of HyMETS Arc-Jet Flows for PICA-N Modeling

Arc-jet testing is a fundamental tool in the screening of thermal protection systems under atmospheric entry conditions. In this work, high-fidelity computational fluid dynamics simulations coupled with machine learning methods for the Hypersonic Materials Environmental Test System arc-jet facility have been carried out. This effort improves the MEDLI-2 data reconstruction fidelity by understanding the impact of the NuSil coating over the PICA thermocouple plugs. Thermochemical non-equilibrium models are employed to simulate the flows inside the arc-jet nozzle and chamber and machine learning techniques are used to calibrate the arc-jet inflow conditions. The material response is simulated with the recession and in-depth numerical results compared to experimental measurements.

STMD↗

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗

Assessment of Pioneer Venus Entry Heating with Coupled Radiation and Ablation

The objective of this work is to reconstruct thermocouple temperature measurements taken at the stagnation point of the Pioneer Venus large probe using uncoupled and coupled ablation calculations. The time history of the large probe forebody entry heating was obtained from a trajectory-based analysis that employed a material response solver to model equilibrium and nonequilibrium ablation chemistry at the surface. In the uncoupled approach, equilibrium ablation chemistry was assumed in the material response calculations, but the material response outputs were not returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer due to ablation (e.g., convective blockage and radiative energy absorption) did not influence the flowfield solutions. In the coupled approach, nonequilibrium ablation chemistry was performed and the material response outputs were returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer were accounted for in the flowfield solutions. Both uncoupled and coupled approaches produced reconstructed temperatures that were within 10% of the flight measurements through the entry heat pulse, with the coupled results comparing more favorably with the measurements. Overall, the reconstructed thermocouple temperatures in this work significantly improved on the results of previous reconstructions for the large probe.

Jannuel V V Cabrera↗

The Paradox of a Wet (High H2O) and Dry (Low H2O/Ce) Mantle: High Water Concentrations in Mantle Garnet Pyroxenites from Hawaii

Water dissolved as trace amounts in anhydrous minerals has a large influence on the melting behavior and physical properties of the mantle. The water concentration of the oceanic mantle is inferred from the analyses of Mid-Ocean Ridge Basalt (MORB) and Oceanic Island Basalt (OIB). but there is little data from actual mantle samples. Moreover, enriched mineralogies (pyroxenites, eclogites) are thought as important sources of heterogeneity in the mantle, but their water concentrations and their effect on the water budget and cycling in the mantle are virtually unknown. Here, we analyzed by FTIR water in garnet clinopyroxenite xenoliths from Salt Lake Crater, Oahu, Hawaii. These pyroxenites are high-pressure (>20kb) crystal fractionates from alkalic melts. The clinopyroxenes (cpx) have 260 to 576 ppm wt H2O, with the least differentiated samples (Mg#>0.8) in the 400-500 ppm range. Orthopyroxene (opx) contain 117-265 ppm H2O, about half of that of cpx, consistent with other natural sample studies, but lower than cpx/opx equilibrium from experimental data. The pyroxenite cpx and opx H2O concentrations are at the high-end of on-and off-craton peridotite xenolith concentrations and those of Hawaiian spinel peridotites. In contrast, garnet has extremely low water contents (<5ppm H2O). There is no correlation between H2O in cpx and lithophile element concentrations. Phlogopite is present in some samples, and its modal abundance shows a positive correlation in Mg# with cpx, implying equilibrium. However, there is no correlation between H2O concentrations and or the presence of phlogopite. These data imply that cpx and opx may be at water saturation, far lower than experimental data suggest. Reconstructed bulk rock pyroxenite H2O ranges from 200-460 ppm (average 331 +/- 75 ppm), 2 to 8 times higher than H2O estimates for the MORB source (50-200 ppm), but in the range of E-MORB, OIB and the source of rejuvenated Hawaiian magmas. The average bulk rock pyroxenite H2O/Ce is 69 +/-35, lower than estimates of the MORB source (approx 150) or FOZO, C (200-250) mantle component, but consistent with "dry" EM sources (<100). These data suggest that a metasomatized, refertilized oceanic lithosphere that contains pyroxenitic veins (e.g. the lower part of an oceanic plate, where ascending melts can become trapped and crystallize), will have both higher water concentrations and low H2O/Ce, and may contribute to EM-type OIB sources, like that of Samoa basalts. Therefore, a low H2O/Ce mantle source may not necessarily be "dry".

Peslier, Anne H.↗

Evidence for Dry Carbonatite Metasomatism in the Oceanic Lithosphere from Peridotite Xenoliths of Samoa and Lanzarote

Water in Earth’s mantle affects processes like magmatism and plate tectonics. Experiments show that CO2-rich fluids lower the water solubility in olivine, implying that CO2-rich melts/fluids may dehydrate the lithosphere during metasomatism. To test this hypothesis, we report water concentrations (by polarized FTIR) of olivines, orthopyroxenes (OPX) and clinopyroxenes (CPX) from Savai’i (Samoa) and Lanzarote (Canary Islands) peridotite xenoliths with evidence of carbonatite metasomatism. Savai’i peridotites are highly depleted harzburgites and dunites with spinel Cr# (Cr/(Cr+Al)) ranging from 0.4 to 0.76 (estimated degree of melting: 191.5%). Strong Light Rare Earth Element (LREE) enrichments with Ti and Zr depletions in OPX and CO2-rich fluid inclusions (via Raman spectroscopy) are consistent with carbonatite metasomatism. Olivine, OPX and reconstructed bulk rock water concentrations (0.67-3.8, 17-89 and 4-26 ppm H2O, respectively) are low and show no apparent relationship with extent of carbonatite metasomatism. Calculated water concentrations of melts in equilibrium with Savai’i OPX (OPX/melt partitioning of water 0.0063 to 0.011) are, on average (0.540.32 wt% H2O), lower than host Samoan lavas (0.63 to 1.5 wt% H2O), despite the LREE enrichments in OPX. Lanzarote peridotites are also highly depleted (degree of melting from spinel Cr#: 171.8%).Water concentrations are low in olivines (1.7-5.3 ppm H2O) and variable in pyroxenes (OPX: 42-103 ppm H2O; CPX: 105-301 ppm H2O), and show no apparent correlation with indicators of carbonatite metasomatism. Both Savai’i and Lanzarote peridotites show negative correlations between water and degree of melting (i.e. Mg/(Mg+Fe), Cr#), suggesting melt depletion rather than metasomatism may have influenced their water concentrations. Calculated water concentrations of melts in equilibrium with Lanzarote CPX (average 1.90.75 wt% H2O; CPX/melt partitioning of water 0.011 to 0.012) are similar to those for Western Canaries lavas (average 1.80.31 wt%; CPX/melt partitioning of water 0.016 to 0.021) inferred from their CPX phenocrysts. However, calculated Ce concentrations in such melts (352 to 378 ppm; CPX/melt partitioning of Ce 0.07) are an order of magnitude greater than the lavas, and similar to carbonatites. This leads to H2O/Ce to be an order of magnitude lower in the inferred melts (26 to 57) than estimates for Western Canary lavas (280150). These low H2O/Ce ratios may suggest H2O loss from CPX during ascent, but the lack of strong water diffusion gradients in Lanzarote minerals does not support this. Instead we hypothesize that carbonatite metasomatism resulted in greater enrichment of Ce over H2O. Assuming carbonatite magmas are water rich, this implies a lower partitioning of water between minerals and melts during metasomatism, as suggested by experiments. Our data suggests carbonatite metasomatism does not result in significant re-hydration of the lithosphere, in contrast to silicate metasomatism as previously observed in Hawaiian peridotites.

Ashley, Aaron↗

Assessing the Impact of Laurentide Ice-sheet Topography on Glacial Climate

Simulations of past climates require altered boundary conditions to account for known shifts in the Earth system. For the Last Glacial Maximum (LGM) and subsequent deglaciation, the existence of large Northern Hemisphere ice sheets caused profound changes in surface topography and albedo. While ice-sheet extent is fairly well known, numerous conflicting reconstructions of ice-sheet topography suggest that precision in this boundary condition is lacking. Here we use a high-resolution and oxygen-isotopeenabled fully coupled global circulation model (GCM) (GISS ModelE2-R), along with two different reconstructions of the Laurentide Ice Sheet (LIS) that provide maximum and minimum estimates of LIS elevation, to assess the range of climate variability in response to uncertainty in this boundary condition.We present this comparison at two equilibrium time slices: the LGM, when differences in ice-sheet topography are maximized, and 14 ka, when differences in maximum ice-sheet height are smaller but still exist. Overall, we find significant differences in the climate response to LIS topography, with the larger LIS resulting in enhanced Atlantic Meridional Overturning Circulation and warmer surface air temperatures, particularly over northeastern Asia and the North Pacific. These up- and downstream effects are associated with differences in the development of planetary waves in the upper atmosphere, with the larger LIS resulting in a weaker trough over northeastern Asia that leads to the warmer temperatures and decreased albedo from snow and sea-ice cover. Differences between the 14 ka simulations are similar in spatial extent but smaller in magnitude, suggesting that climate is responding primarily to the larger difference in maximum LIS elevation in the LGM simulations. These results suggest that such uncertainty in ice-sheet boundary conditions alone may significantly impact the results of paleoclimate simulations and their ability to successfully simulate past climates, with implications for estimating climate sensitivity to greenhouse gas forcing utilizing past climate states.

glaciology↗

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

Does sedimentary organic delta 13C record variations in quaternary ocean [CO2(aq)]?

Ocean surface water [CO2(aq)] variations based on glacial/interglacial changes in sediment delta 13Corg are shown to compare favorably with reconstructions based on ice core [CO2]. In particular, an approximate 80 microatmospheres increase in atmospheric pCO2 during the last glacial-interglacial transition is calculated to correspond to a 3-4 micromolar increase in ocean surface water [CO2(aq)] at atmospheric equilibrium. A widespread marine delta 13Corg decrease of 1-2% accompanied this event and was not preceded by an equivalent isotopic change in surface water total dissolved inorganic carbon. These observations support the hypothesis that [CO2(aq)] influences photosynthetic isotope fractionation between marine inorganic and organic carbon pools, and therefore that plankton/sediment delta 13Corg may serve as a proxy for surface water [CO2(aq)].

NASA Program Exobiology↗

InSight Aerothermal Environment Assessment

The InSight spacecraft was proposed to be a build-to-print copy of the Phoenix vehicle due to the knowledge that the lander payload would be similar and the trajectory would be similar. However, the InSight aerothermal analysts, based on tests performed in CO2 during the Mars Science Laboratory mission (MSL) and completion of Russian databases, considered radiative heat flux to the aftbody from the wake for the first time for a US Mars mission. The combined convective and radiative heat flux was used to determine if the as-flown Phoenix thermal protection system (TPS) design would be sufficient for InSight. All analyses showed that the design would be adequate. Once the InSight lander was successfully delivered to Mars on November 26, 2018, work began to reconstruct the atmosphere and trajectory in order to evaluate the aerothermal environments that were actually encountered by the spacecraft and to compare them to the design environments.The best estimated trajectory (BET) reconstructed for the InSight atmospheric entry fell between the two trajectories considered for the design, when looking at the velocity versus altitude values. The maximum heat rate design trajectory (MHR) flew at a higher velocity and the maximum heat load design trajectory (MHL) flew at a lower velocity than the BET. For TPS sizing, the MHL trajectory drove the design. Reconstruction has shown that the BET flew for a shorter time than either of the design environments, hence total heat load on the vehicle should have been less than used in design. Utilizing the BET, both DPLR and LAURA were first run to analyze the convective heating on the vehicle with no angle of attack. Both codes were run with axisymmetric, laminar flow in radiative equilibrium and vibrational non-equilibrium with a surface emissivity of 0.8. Eight species Mitcheltree chemistry was assumed with CO2, CO, N2, O2, NO, C, N, and O. Both codes agreed within 1% on the forebody and had the expected differences on the aftbody. The NEQAIR and HARA codes were used to analyze the radiative heating on the vehicle using full spherical ray-tracing. The codes agreed within 5% on most aftbody points of interest.The LAURA code was then used to evaluate the conditions at angle of attack at the peak heating and peak pressure times. Boundary layer properties were investigated to confirm that the flow over the forebody was laminar for the flight.Comparisons of the aerothermal heating determined for the reconstructed trajectory to the design trajectories showed that the as-flown conditions were less severe than design

Beck, Robin↗

InSight's Reconstructed Aerothermal Environments

The InSight spacecraft was proposed to be a build-to-print copy of the Phoenix vehicle due to the knowledge that the lander payload would be similar and the trajectory would be similar. However, the InSight aerothermal analysts, based on tests performed in CO2 during the Mars Science Laboratory mission (MSL) and completion of Russian databases, considered radiative heat flux to the aftbody from the wake for the first time for a US Mars mission. The combined convective and radiative heat flux was used to determine if the as-flown Phoenix thermal protection system (TPS) design would be sufficient for InSight. All analyses showed that the design would be adequate. Once the InSight lander was successfully delivered to Mars on November 26, 2018, work began to reconstruct the atmosphere and trajectory in order to evaluate the aerothermal environments that were actually encountered by the spacecraft and to compare them to the design environments.The best estimated trajectory (BET) reconstructed for the InSight atmospheric entry fell between the two trajectories considered for the design, when looking at the velocity versus altitude values. The maximum heat rate design trajectory (MHR) flew at a higher velocity and the maximum heat load design trajectory (MHL) flew at a lower velocity than the BET. For TPS sizing, the MHL trajectory drove the design. Reconstruction has shown that the BET flew for a shorter time than either of the design environments, hence total heat load on the vehicle should have been less than used in design. Utilizing the BET, both DPLR and LAURA were first run to analyze the convective heating on the vehicle with no angle of attack. Both codes were run with axisymmetric, laminar flow in radiative equilibrium and vibrational non-equilibrium with a surface emissivity of 0.8. Eight species Mitcheltree chemistry was assumed with CO2, CO, N2, O2, NO, C, N, and O. Both codes agreed within 1% on the forebody and had the expected differences on the aftbody. The NEQAIR and HARA codes were used to analyze the radiative heating on the vehicle using full spherical ray-tracing. The codes agreed within 5% on most aftbody points of interest.The LAURA code was then used to evaluate the conditions at angle of attack at the peak heating and peak pressure times. Boundary layer properties were investigated to confirm that the flow over the forebody was laminar for the flight.Comparisons of the aerothermal heating determined for the reconstructed trajectory to the design trajectories showed that the as-flown conditions were less severe than design

radiation↗

Microstructures of Hibonite From an ALH A77307 (CO3.0) CAI: Evidence for Evaporative Loss of Calcium

Hibonite is a comparatively rare, primary phase found in some CAIs from different chondrite groups and is also common in Wark-Lovering rims [1]. Hibonite is predicted to be one of the earliest refractory phases to form by equilibrium condensation from a cooling gas of solar composition [2] and, therefore, can be a potential recorder of very early solar system processes. In this study, we describe the microstructures of hibonite from one CAI in ALH A77307 (CO3.0) using FIB/TEM techniques in order to reconstruct its formational history.

Han, Jangmi↗