Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK LAYER”

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 721 records · Page 40

Simulations of Hayabusa2 Atmospheric Entry and Comparisons with Data from the Imaging Campaign

Results from flow and radiation simulations for the Hayabusa2 Sample Return Capsule (SRC) along its best-estimated trajectory (BET) are presented. Flow fields are first computed at several points along the best-estimated trajectory for a freestream consisting of N2, O2, Ar, CO2, and H2O; the last three are trace species in Earth’s atmosphere. The thermal response of an assumed ablative material is estimated and flow fields are recomputed at each time point with modified surface temperature distributions. The spectral radiance, from the combination of an aeroheated capsule and the high-temperature shock layer around it, is computed on an imaging plane for a given view angle and slant range at each time point. The radiance is transported through the intervening atmosphere between the Hayabusa2 SRC and the observing aircraft and upon integration over solid angle representing the field of view of the instruments converted to a spectral irradiance. The spectral irradiance is convolved with the instrument function to obtain results that can be directly compared with measurement. Fair to good comparisons with measurements required shifting the time coordinate of measurements by 7.5 s. The disagreement between prediction and measurement is substantial after about 72.5 s along the BET. Additional computations seem to support the idea of flow transition over the capsule.

Atmospheric entry↗

Assessment of Particle Drag Models Using GBK Experimental Data

Developing a simulation capability to model dust-laden flows requires the ability to accurately predict the velocity of a particle as it travels through a shock layer, nozzle plume, or the flow of an experimental facility. The primary mechanism that determines the particle velocity is the drag force acting on the particle. The drag force is typically expressed in terms of the frontal area of the particle, its velocity relative to the surrounding fluid, the density of the surrounding fluid, and a non-dimensional drag coefficient. As part of the ESM project, an integrated CFD-particle trajectory code named US3D-DUST is being developed]. This code uses a Lagrangian-based framework that computes the trajectories of individual particles in the flow. Several of the most commonly-used spherical drag models have been implemented in US3D-DUST. In this study, the US3D-DUST code will be used to simulate particle flows in the GBK facility. Differences in particle velocities predicted using the different drag models will be compared against recent GBK data

Dusty flow↗

Simulations of Hayabusa2 Atmospheric Entry and Comparisons with Data from the Imaging Campaign

Results from flow and radiation simulations for the Hayabusa2 Sample Return Capsule (SRC) along its best-estimated trajectory (BET) are presented. Flow fields are first computed at several points along the best-estimated trajectory for a freestream consisting of N2, O2, Ar, CO2, and H2O; the last three are trace species in Earth’s atmosphere. The thermal response of an assumed ablative material is estimated and flow fields are recomputed at each time point with modified surface temperature distributions. The spectral radiance, from the combination of an aeroheated capsule and the high-temperature shock layer around it, is computed on an imaging plane for a given view angle and slant range at each time point. The radiance is transported through the intervening atmosphere between the Hayabusa2 SRC and the observing aircraft and upon integration over solid angle representing the field of view of the instruments converted to a spectral irradiance. The spectral irradiance is convolved with the instrument function to obtain results that can be directly compared with measurement. Fair to good comparisons with measurements required shifting the time coordinate of measurements by 7.5 s. The disagreement between prediction and measurement is substantial after about 72.5 s along the BET. Additional computations seem to support the idea of flow transition over the capsule.

Atmospheric entry↗

New Features of the NEQAIR Radiation Code

The longest-lived code for predicting shock layer radiation, NEQAIR, is now in its 5th decade of service. Substantial changes to the code have been made over the previous decade, the most recent report of which was at the 5th Workshop on Radiation in High Temperature Gases in 2014, for the version referred to as NEQAIR14. This paper will review some of the improvements made to the NEQAIR code since then, which is now at v15.2. Some of these features are discussed briefly below. NEQAIR15 and subsequent versions have enabled parallel evaluation of multiple lines of sight. This is accomplished by utilizing the HDF5 file format and placing multiple lines into a single file, LOS.h5, which is used for both input and output. This approach enables straightforward parallel execution both over the number of lines of sight and the number of points per line. For large problems, runtime reduces linearly with the number of nodes deployed since each line is processed independently by a subset of MPI ranks. Three applications of the multi-line solver are discussed. The first has to do with performing loosely coupled radiation-flowfield solutions. In this case the computed absorption and emission coefficients are used to evaluate the total energy absorbed or emitted at each point, allowing evaluation of the volumetric source term in the flowfield. The second computation is for obtaining heat flux from nonuniform flows, which require integration over spherical co-ordinates. These are of particular interest for evaluating radiation on the vehicle backshell. This 3D option improves the angular integration scheme and allows adaptive line selection that together reduce the number of lines required by about an order of magnitude. The final application is for remote observation, which is essentially the 3D integration problem over a small solid angle. For all three of these computations, data can be stored in the HDF5 file which allows a NEQAIR run to be restarted when it times out, or to add atmospheric absorption or instrument scan functions. An additional level of parallelism is enabled in NEQAIR15.2 using GPU routines. The GPU parallelism has realized up to 8x speed-up when running on a single core but diminishes as CPU parallelism is increased. For running multi-line simulations, it may be easier to reserve a large number of CPU nodes than to obtain the number of GPU nodes required for similar performance. A GUI, known as NEQTPY, allows for reading and creating input files, running NEQAIR, and displaying results. A significant feature of NEQTPY is the ability to perform spectral fits to data. The fits can operate on a single line spectrum (radiance vs. wavelength) or a 3D input file with multiple columns of data. Other new features include improved constants, additional species, more detailed non-Boltzmann modelling, advanced user controls, the ability to read and calculate spectra from HITRAN datafiles, photodissociation and photoionization cross-sections. A “fast” automatic grid option may reduce the size and time of spectral calculations while still maintaining good accuracy for total heat flux.

Brett A Cruden↗

The NASA Ames Aerothermodynamics Branch

The Ames Aerothermodynamics branch provides NASA and the country an integrated modeling, simulation, and test capability in the areas of aerothermodynamics, shock layer radiation, and entry system design. With roots dating back to the founding of Ames in the 1940s, the branch has contributed to every major NASA mission with an atmospheric entry component and actively partners with commercial space companies to assist with entry system development and technology transfer. This presentation provides a brief overview of the branch, its core competencies, and recent accomplishments.

Jeffrey Hill↗

New Developments In NASA’s Entry Systems Modeling Project

This paper describes recent developments for modeling entry, descent, and landing (EDL)of spacecraft in support of NASA’s exploration missions. Mission-specific research and model development for entry systems occurs across the Agency (e.g., within flight programs like Artemis/Orion and Mars Sample Return), however the aim of this paper is to discuss the research conducted by NASA’s Entry Systems Modeling (ESM) Project, which serves as the Agency’s only effort dedicated to advancing modeling capabilities that cross-cut multiple technical disciplines, missions and destinations. The ESM portfolio is developed to address the specific needs expressed by a cross-section of NASA stakeholders, including flight and research projects, technical leadership, and subject matter experts. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Thermal Protection System (TPS) material modeling, (2) Shock layer kinetics and radiation, (3) Aerosciences, and (4) Guidance, navigation, and control. In addition to the core technical areas, special topics are rolled into the project portfolio as specific demands arise. Current special topics include TPS Certification by Analysis, improving understanding of woven TPS material performance; Hypersonic Wake Flows, assessing and improving predictions of base flows; and the MEDLI2 Deep Dive, furthering analysis of data obtained during Mars2020 mission’s entry and descent at Mars. Key results from each of these areas are presented in this paper, along with associated references to serve as a roadmap for other EDL researchers to access NASA’s publications.

Aaron M Brandis↗

The Effect of Drag Model on Heatshield Recession due to Particle Impacts for Martian Spacecraft

A spacecraft entering the Martian atmosphere during a dust storm may experience recession to the heatshield due to dust particle impacts. Aerodynamic drag is the primary force that determines the trajectory of the dust particles through the shock layer. This paper examines the effect of particle drag model on the heatshield recession. Three particle drag models are assessed including two that are intended to be applicable over a wide range of particle flow conditions. Particle trajectories are computed in conditions measured during the 2007 major global dust storm. It was found that accounting for Knudsen number and compressibility effects made a large difference in the estimated particle impact velocity. The two drag models that were valid for transitional, compressible particle flow environments predicted only slightly-different amounts of heatshield recession due to dust particle impacts. A brief description of the effects of non-spherical particles on drag coefficients is provided.

Heatshield↗

The Effect of Drag Model on Heatshield Recession due to Particle Impacts for Martian Spacecraft

A spacecraft entering the Martian atmosphere during a dust storm may experience recession to the heatshield due to dust particle impacts. Aerodynamic drag is the primary force that determines the trajectory of the dust particles through the shock layer. This paper examines the effect of particle drag model on the heatshield recession. Three particle drag models are assessed including two that are intended to be applicable over a wide range of particle flow conditions. Particle trajectories are computed in conditions measured during the 2007 major global dust storm. It was found that accounting for Knudsen number and compressibility effects made a large difference in the estimated particle impact velocity. The two drag models that were valid for transitional, compressible particle flow environments predicted only slightly-different amounts of heatshield recession due to dust particle impacts. A brief description of the effects of non-spherical particles on drag coefficients is provided.

Heatshield↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Improved Chemistry and Attenuation Models for Communication Blackout Simulation During Mars 2020 Entry

As a blunt body enters a planetary atmosphere, a plasma forms in the hypersonic shock layer and attenuates radio communication causing signal blackout for some duration of the entry sequence. In our previous work,1 computational fluid dynamics (CFD) was applied to model the entry flow around the Mars 2020 spacecraft, including ionization and electron density throughout the flow field, and predict ultra-high frequency (UHF) radio wave attenuation due to electrons. In total, 17 chemical species and their spatial profiles are modelled around the Mars 2020 spacecraft at 11 different points in time during entry. Although the simulation predicted the onset of attenuation well, the timing of the end of the predicted blackout window significantly preceded the end time observed during the 2021 landing. The present work seeks to improve the attenuation model by accounting for the fact that electrons undergo collisions with heavier species in the flow, which is an effect that was neglected in previous analyses. It is determined that including electron collisions increases the overall magnitude of attenuation predicted especially towards the end of the measured attenuation period, improving qualitative agreement between predicted and measured attenuation to both spacecraft receiving the signal from Mars 2020. To explore the remaining uncertainty in signal attenuation predictions further, a sensitivity study is performed to investigate the impact of associative ionization and electron-impact ionization rate coefficients on the electron density predicted by CFD and on the resulting attenuation predictions. These coefficients are believed to contain up to order-of-magnitude uncertainty, and therefore may significantly affect the number density of electrons throughout the flow field. Variations in associative ionization coefficients demonstrate significant impact on the magnitude of attenuation due to variation in the electron density coming from associative ionization. However, the start and end times of the predicted signal attenuation period are only slightly impacted by said variation.

Eve Papajak↗

Simulations of OSIRIS-REx Atmospheric Entry for Comparisons with Data from the Imaging Campaign

Results from pre-flight flow and radiation simulations for the OSIRIS-REx Sample Return Capsule along a candidate entry trajectory are presented. These results are from pre-flight simulations in support of the observation campaign. Details of the methodology for computing temporal variation of spectral irradiance on an instrument onboard an aerial observation platform are presented along with a sampling of results. An attempt is made to use the methodology presented to assess the level of elemental sodium in the shock layer. The final paper will be based on the actual flight trajectory and will include comparisons of predictions with calibrated spectrometric data.

Spectrometric Imaging↗

Improved Chemistry and Attenuation Models for Communication Blackout Simulation During Mars 2020 Entry

As a blunt body enters a planetary atmosphere, a plasma forms in the hypersonic shock layer and attenuates radio communication causing signal blackout for some duration of the entry sequence. In our previous work,1 computational fluid dynamics (CFD) was applied to model the entry flow around the Mars 2020 spacecraft, including ionization and electron density throughout the flow field, and predict ultra-high frequency (UHF) radio wave attenuation due to electrons. In total, 17 chemical species and their spatial profiles are modelled around the Mars 2020 spacecraft at 11 different points in time during entry. Although the simulation predicted the onset of attenuation well, the timing of the end of the predicted blackout window significantly preceded the end time observed during the 2021 landing. The present work seeks to improve the attenuation model by accounting for the fact that electrons undergo collisions with heavier species in the flow, which is an effect that was neglected in previous analyses. It is determined that including electron collisions increases the overall magnitude of attenuation predicted especially towards the end of the measured attenuation period, improving qualitative agreement between predicted and measured attenuation to both spacecraft receiving the signal from Mars 2020. To explore the remaining uncertainty in signal attenuation predictions further, a sensitivity study is performed to investigate the impact of associative ionization and electron-impact ionization rate coefficients on the electron density predicted by CFD and on the resulting attenuation predictions. These coefficients are believed to contain up to order-of-magnitude uncertainty, and therefore may significantly affect the number density of electrons throughout the flow field. Variations in associative ionization coefficients demonstrate significant impact on the magnitude of attenuation due to variation in the electron density coming from associative ionization. However, the start and end times of the predicted signal attenuation period are only slightly impacted by said variation.

Eve Papajak↗

Testing of Two Mars Powered Descent Vehicle Concepts in the Langley Unitary Plan Wind Tunnel

Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). Each model was designed to accommodate up to eight nozzles, with the HIAD model having five different nozzle configurations to investigate the effects of nozzle location, cant angle, and area ratio. The models were tested with high pressure air as the nozzle plume gas, and included the following instrumentation: high-speed video, discrete steady state and high-frequency pressure, pressure sensitive paint, and a new flow-through force and moment balance for the HIAD model. The high-speed imagery showed the overall expected growth of the shock layer to increasing thrust levels. The discrete stagnation pressure data on the HIAD model with four and eight blowing nozzles was sensitive to thrust coefficient with the nozzles canted at 0 degrees, whereas having nozzles canted outward 20 degrees and/or being located closer to the heatshield shoulder largely removed that sensitivity. The CobraMRV model stagnation pressure was more sensitive to tunnel Mach number and sideslip angle, due to the nozzle arrangement and plume interference. Pressure sensitive paint data quality was compromised by paint damage from frequent model changes, especially for the HIAD model. However, pressure coefficient data on the CobraMRV model showed the same trends as the discrete pressure measurements, as well as a sensitivity to non-zero sideslip angles. The aerodynamic force coefficients were derived from the pressure sensitive paint data. On the HIAD models, the primary force coefficient decreased with increasing thrust due to the nozzle plumes blocking flow to the heatshield area surrounding the nozzle exist. On the CobraMRV model, the force coefficient was relatively insensitive to thrust coefficient at the lower Mach number when sideslip angle was 0 degrees The force coefficient decreases with increasing thrust coefficient when the sideslip angle is 10 degrees. Balance data quality was negatively impacted by thermal drift issues that were not apparent in pre-test calibration measurements, thus preventing usable test data. Results and lessons learned will be used to take further technology development steps, including more advanced ground test techniques and flight testing.

Karl T Edquist↗

Aerothermal Analysis and Environment Predictions for the Mars Sample Retrieval Lander (SRL)

Mars Sample Retrieval Lander (SRL) is the next mission to Mars, and an integral part of the proposed Mars Sample Return (MSR) Program. Aerothermal Analysis of the SRL capsule takes advantage of the design and analysis of the Mars Science Laboratory and the Mars 2020 missions, findings from the MEDLI and MEDLI2 heatshield instrumentation campaigns, and developments in predictive capabilities over the last 20 years. In particular, SRL is being designed to enter the Mars atmosphere at velocities as high as 8 km/s, which would be the highest for a Mars entry, and is expected to encounter additional shock layer radiation physics compared to previous missions to Mars. This paper presents the status of analysis including the overall methodology, models and assumptions of the aerothermal environment predictions, with a focus on differences from the approaches and modeling used for Mars 2020 and MSL.

Mars↗

Testing of Two Mars Powered Descent Vehicle Concepts in the Langley Unitary Plan Wind Tunnel

Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). Each model was designed to accommodate up to eight nozzles, with the HIAD model having five different nozzle configurations to investigate the effects of nozzle location, cant angle, and area ratio. The models were tested with high pressure air as the nozzle plume gas, and included the following instrumentation: high-speed video, discrete steady state and high-frequency pressure, pressure sensitive paint, and a new flow-through force and moment balance for the HIAD model. The high-speed imagery showed the overall expected growth of the shock layer to increasing thrust levels. The discrete stagnation pressure data on the HIAD model with four and eight blowing nozzles was sensitive to thrust coefficient with the nozzles canted at 0 degrees, whereas having nozzles canted outward 20 degrees and/or being located closer to the heatshield shoulder largely removed that sensitivity. The CobraMRV model stagnation pressure was more sensitive to tunnel Mach number and sideslip angle, due to the nozzle arrangement and plume interference. Pressure sensitive paint data quality was compromised by paint damage from frequent model changes, especially for the HIAD model. However, pressure coefficient data on the CobraMRV model showed the same trends as the discrete pressure measurements, as well as a sensitivity to non-zero sideslip angles. The aerodynamic force coefficients were derived from the pressure sensitive paint data. On the HIAD models, the primary force coefficient decreased with increasing thrust due to the nozzle plumes blocking flow to the heatshield area surrounding the nozzle exist. On the CobraMRV model, the force coefficient was relatively insensitive to thrust coefficient at the lower Mach number when sideslip angle was 0 degrees The force coefficient decreases with increasing thrust coefficient when the sideslip angle is 10 degrees. Balance data quality was negatively impacted by thermal drift issues that were not apparent in pre-test calibration measurements, thus preventing usable test data. Results and lessons learned will be used to take further technology development steps, including more advanced ground test techniques and flight testing.

Supersonic Retropropulsion↗

Aerothermal Analysis and Environment Predictions for the Mars Sample Retrieval Lander (SRL)

Mars Sample Retrieval Lander (SRL) is the next mission to Mars, and an integral part of the proposed Mars Sample Return (MSR) Program. Aerothermal Analysis of the SRL capsule takes advantage of the design and analysis of the Mars Science Laboratory and the Mars 2020 missions, findings from the MEDLI and MEDLI2 heatshield instrumentation campaigns, and developments in predictive capabilities over the last 20 years. In particular, SRL is being designed to enter the Mars atmosphere at velocities as high as 8 km/s, which would be the highest for a Mars entry, and is expected to encounter additional shock layer radiation physics compared to previous missions to Mars. This paper presents the status of analysis including the overall methodology, models and assumptions of the aerothermal environment predictions, with a focus on differences from the approaches and modeling used for Mars 2020 and MSL.

Mars↗

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

Titan’s atmosphere is composed mostly of N 2 with a small amount of CH 4 , and so, shock layers around craft entering Titan’s atmosphere will contain a variety of molecules formed from H, C, and N atoms, including the cyanogen radical CN. Sensitivity analysis has shown that the radiative heat flux predicted by computational fluid dynamics (CFD) simulations of Titan entry has up to 14% uncertainty due to the rate coefficients for collisional (de)excitation reactions that control the population of CN in its first and second excited states. The red and violet emission bands from CN’s first and second excited states, respectively, are known to be large sources of radiative heat flux on capsules entering Titan’s atmosphere.[2, 3] So, the simulated population of CN in its first and second excited states is very important, but currently has some inherent uncertainty coming from the data for the rate coefficients for reaction 1. The goal of the present project is to provide improved rate coefficient data for these reactions from first principles quantum chemistry calculations. This work reports on preliminary electronic structure calculations generated at a large number of triatomic geometries of interest, which show multiple avoided crossings at collinear arrangements. This suggests that collisional (de)excitation of CN by N atoms is likely to proceed through these geometries.

Eric Geistfeld↗