Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hypersonic Cone”

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.

484 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Techniques for Transition and Surface Temperature Measurements on Projectiles at Hypersonic Velocities- A Status Report

A research effort to advance techniques for determining transition location and measuring surface temperatures on graphite-tipped projectiles in hypersonic flight in a ballistic range is described. Projectiles were launched at muzzle velocities of approx. 4.7 km/sec into air at pressures of 190-570 Torr. Most launches had maximum pitch and yaw angles of 2.5-5 degrees at pressures of 380 Torr and above and 3-6 degrees at pressures of 190-380 Torr. Arcjet-ablated and machined, bead-blasted projectiles were launched; special cleaning techniques had to be developed for the latter class of projectiles. Improved methods of using helium to remove the radiating gas cap around the projectiles at the locations where ICCD (intensified charge coupled device) camera images were taken are described. Two ICCD cameras with a wavelength sensitivity range of 480-870 nm have been used in this program for several years to obtain images. In the last year, a third camera, with a wavelength sensitivity range of 1.5-5 microns [in the infrared (IR)], has been added. ICCD and IR camera images of hemisphere nose and 70 degree sphere-cone nose projectiles at velocities of 4.0-4.7 km/sec are presented. The ICCD images clearly show a region of steep temperature rise indicative of transition from laminar to turbulent flow. Preliminary temperature data for the graphite projectile noses are presented.

Wilder, M. C.↗
NASA NTRS2016

Nano-ADEPT Aeroloads Wind Tunnel Test

A wind tunnel test of the Adaptable Deployable Entry and Placement Technology (ADEPT) was conducted in April 2015 at the US Army's 7 by10 Foot Wind Tunnel located at NASA Ames Research Center. Key geometric features of the fabric test article were a 0.7 meter deployed base diameter, a 70 degree half-angle forebody cone angle, eight ribs, and a nose-to-base radius ratio of 0.7. The primary objective of this wind tunnel test was to obtain static deflected shape and pressure distributions while varying pretension at dynamic pressures and angles of attack relevant to entry conditions at Earth, Mars, and Venus. Other objectives included obtaining aerodynamic force and moment data and determining the presence and magnitude of any dynamic aeroelastic behavior (buzz/flutter) in the fabric trailing edge. All instrumentation systems worked as planned and a rich data set was obtained. This paper describes the test articles, instrumentation systems, data products, and test results. Four notable conclusions are drawn. First, test data support adopting a pre-tension lower bound of 10 foot pounds per inch for Nano-ADEPT mission applications in order to minimize the impact of static deflection. Second, test results indicate that the fabric conditioning process needs to be reevaluated. Third, no flutter/buzz of the fabric was observed for any test condition and should also not occur at hypersonic speeds. Fourth, translating one of the gores caused ADEPT to generate lift without the need for a center of gravity offset. At hypersonic speeds, the lift generated by actuating ADEPT gores could be used for vehicle control.

system engineering↗
DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗
NASA NTRS2019

The Multi-Mission Earth Entry Vehicle for Sample Return Missions – Past, Present, and Future

The Multi-Mission Earth Entry Vehicle (MMEEV) is an enabling technology developed at NASA’s Langley Research Center (LaRC) over the last two decades for returning samples to Earth across a wide array of space science missions. Currently, the MMEEV is being considered for NASA’s Mars Sample Return (MSR) mission. The original vehicle concept, the Earth Entry Vehicle (EEV), was innovated at LaRC in 1998 as a robust solution to return Mars soil samples to Earth under stringent backward contamination requirements. These backward contamination requirements drove the EEV to have higher reliability than any capsule previously designed for a return-to-Earth sample return mission. The EEV achieved this high reliability by employing a passive (no active systems) vehicle architecture optimized for fault tolerance in a compact, low-mass configuration that is extensible to virtually any sample return mission. The original EEV concept utilized a carbon-carbon primary structure with high-density carbon phenolic thermal protection system. The capsule had a 60-degree sphere-cone forebody and a backshell geometry uniquely tailored to produce aerodynamics that would passively re-orient the vehicle if it entered the atmosphere with an off-nominal attitude. Contrary to every other sample return capsule conceived at the time, the EEV was designed to land without a parachute. The vehicle incorporated an integral energy-absorbing crushable structure that protected the Mars sample for landings on surfaces ranging from soft soil to solid concrete. This paper describes 20 years of technological advancements LaRC has incorporated into the EEV architecture to evolve it from the original, MSR-enabling vehicle, to a true multi-mission capability relevant to any sample return mission. The vehicle’s unique Integrated Composite Stiffener Structure (ICoSS) has been optimized for specific strength - supporting high-G atmospheric entries with steep entry angles that produce precise landing footprints on the ground. The vehicle geometry has been refined through wind tunnel testing and computational fluid dynamics simulations to improve the vehicle’s aerodynamic stability and robustness to off-nominal conditions from hypersonic to subsonic flight. The resulting configuration of the current MMEEV architecture is described, with details provided on its sample carrying capacity and atmospheric entry trajectory capabilities. The upgraded vehicle performance is mapped into current space science objectives, showing how the MMEEV supports future sample return missions and continues to be an enabling technology for NASA’s vision to return samples from Mars.

J M Corliss↗
NASA NTRS2012

Secondary Instability of Second Modes in Hypersonic Boundary Layers

Second mode disturbances dominate the primary instability stage of transition in a number of hypersonic flow configurations. The highest amplification rates of second mode disturbances are usually associated with 2D (or axisymmetric) perturbations and, therefore, a likely scenario for the onset of the three-dimensionality required for laminar-turbulent transition corresponds to the parametric amplification of 3D secondary instabilities in the presence of 2D, finite amplitude second mode disturbances. The secondary instability of second mode disturbances is studied for selected canonical flow configurations. The basic state for the secondary instability analysis is obtained by tracking the linear and nonlinear evolution of 2D, second mode disturbances using nonlinear parabolized stability equations. Unlike in previous studies, the selection of primary disturbances used for the secondary instability analysis was based on their potential relevance to transition in a low disturbance environment and the effects of nonlinearity on the evolution of primary disturbances was accounted for. Strongly nonlinear effects related to the self-interaction of second mode disturbances lead to an upstream shift in the upper branch neutral location. Secondary instability computations confirm the previously known dominance of subharmonic modes at relatively small primary amplitudes. However, for the Purdue Mach 6 compression cone configuration, it was shown that a strong fundamental secondary instability can exist for a range of initial amplitudes of the most amplified second mode disturbance, indicating that the exclusive focus on subharmonic modes in the previous applications of secondary instability theory to second mode primary instability may not have been fully justified.

Li, Fei↗
NASA NTRSDate not supplied

Isolator Shock Train Response to Dual-Mode Scramjet Throttling-Like Forcing

Typical dual-mode scramjet designs employ an isolator component situated between the inlet and combustor. The isolator fluidically buffers the inlet from the combustor via the formation of a shock train, an interconnected series of shock-boundary layer interactions. The isolator is required during ramjet-engine operation and at the lower edges of the scramjet operational envelope. In an operational engine, transient behavior imposes fluidic dynamics due to attitude adjustments, combustor throttling, aberrant resonances within the combustor, etc. To accommodate these transients, the shock train adjusts on both temporal and spatial scales. Sufficiently large downstream pressure increases in the combustor play a critical role during unstart, whereby strong compression waves propagate upstream in the isolator to ultimately force the shock train out the inlet, if it cannot stabilize to the pressure increase. This critically reduces the engine mass flow rate, potentially causing a catastrophic reduction in thrust and loss of the vehicle. A nearly instantaneous unstart can, in principle, be predicted by a ratio of the normal shock pressure jump in an inviscid isolator. However, lower-magnitude perturbations, below this bound, also induce unstarts. Understanding the transient behavior inherent in unsteady isolator operation is critical to reducing safety margin in isolator length and refining engine control algorithms to improve engine performance. The Isolator Dynamics Research Lab (IDRL) is a cold-flow, direct-connect isolator facility ideal for assessing the fundamental fluid dynamics of shock trains. The internal flow of a scramjet isolator is modeled using a 2D, symmetric Mach 2.5 nozzle followed by a 695 mm long, 25.4x50.8 mm rectangular isolator test section. A blockage cone on a high-speed ball screw actuator downstream of the isolator is pulsed to simulate the throttling behavior in an operational combustor. Wall flush high-speed static pressure transducers detect the streamwise shock train within the isolator. Time delays between the onset of the throttling and the response motion of the shock train are analyzed to assess the controlling scales and parameters. From this analysis, the response time of the isolator shock train system is found to depend primarily on the isolator back pressure time-rate-of-change but quickly saturates to a maximum propagation speed along the isolator. These results inform dual-mode scramjet engine designers of the critical parameters and limits of an unstart resistant engine.

scramjet↗
NASA NTRS1989

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗
NASA NTRS2020

Shape Optimization of Vortex Generators to Control Mack Mode Amplification

This paper demonstrates the potential to use shape optimization for the design of vortex generators in an axisymmetric boundary layer. This shape optimization increases the amplitude of stationary streaks created by streamwise vorticity, with the goal of reducing the amplification of Mack mode instabilities. The test case under consideration matches a trajectory point during the ascent phase of the HIFiRE-1 flight experiment. Wall-mounted vortex generators are added to the cone and their shape is optimized to control the amplification of Mack mode instabilities that are known to initiate laminar-turbulent transition in this flow. An empirical objective function is developed based on the previous studies of optimal streaks. A constraint on the maximum streak amplitude is added to avoid the initiation of bypass transition. The shape optimization is conducted using SU2, an open-source suite for multiphysics simulation and design. A significant improvement is observed in an integral metric of the streak amplitude, while the maximum streak amplitude is maintained close to the baseline level. The qualitative features of the optimized geometry are discussed.

Vortex generators↗
NASA NTRSDate not supplied

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗
NASA NTRS2011

Axisymmetric Implementation for 3D-Based DSMC Codes

The primary objective in developing NASA s DSMC Analysis Code (DAC) was to provide a high fidelity modeling tool for 3D rarefied flows such as vacuum plume impingement and hypersonic re-entry flows [1]. The initial implementation has been expanded over time to offer other capabilities including a novel axisymmetric implementation. Because of the inherently 3D nature of DAC, this axisymmetric implementation uses a 3D Cartesian domain and 3D surfaces. Molecules are moved in all three dimensions but their movements are limited by physical walls to a small wedge centered on the plane of symmetry (Figure 1). Unfortunately, far from the axis of symmetry, the cell size in the direction perpendicular to the plane of symmetry (the Z-direction) may become large compared to the flow mean free path. This frequently results in inaccuracies in these regions of the domain. A new axisymmetric implementation is presented which aims to solve this issue by using Bird s approach for the molecular movement while preserving the 3D nature of the DAC software [2]. First, the computational domain is similar to that previously used such that a wedge must still be used to define the inflow surface and solid walls within the domain. As before molecules are created inside the inflow wedge triangles but they are now rotated back to the symmetry plane. During the move step, molecules are moved in 3D but instead of interacting with the wedge walls, the molecules are rotated back to the plane of symmetry at the end of the move step. This new implementation was tested for multiple flows over axisymmetric shapes, including a sphere, a cone, a double cone and a hollow cylinder. Comparisons to previous DSMC solutions and experiments, when available, are made.

Stewart, Benedicte↗
NASA NTRS2018

Görtler Instability and Its Control via Surface Suction over an Axisymmetric Cone at Mach 6

The characteristics of Görtler instability over an axisymmetric cone with an aft concave section are studied via linear and nonlinear instability analysis and direct numerical simulations. Several options for the cone geometry have been investigated numerically, subject to a fixed forecone section and constraints on the maximum cone diameter, overall cone length, and minimum N-factor for the most amplified Görtler modes. Computations show that it is possible to design a cone with a peak N-factor of Nmax > 8 at the target Reynolds number of 12.1×10(exp 6) per meter, corresponding to the maximum quiet Reynolds number in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. Direct numerical simulations show that an array of roughness elements corresponding to a peak roughness height of 0.1006 mm at the center can excite Görtler vortices that evolve into sufficiently strong streamwise streaks that may break down via high-frequency secondary instability. Thus, the selected axisymmetric configuration of interest should provide an acceptable baseline to investigate the feasibility of several aspects of laminar-flow control via boundary-layer suction. The apparatus that is being used for measurements in the Boeing/AFOSR Quiet Tunnel is described, along with some preliminary experimental results.

Gortler vortices↗
NASA NTRS1997

Experiments in Aircraft Roll-Yaw Control using Forebody Tangential Blowing

Advantages of flight at high angles of attack include increased maneuverability and lift capabilities. These are beneficial not only for fighter aircraft, but also for future supersonic and hypersonic transport aircraft during take-off and landing. At high angles of attack the aerodynamics of the vehicle are dominated by separation, vortex shedding and possibly vortex breakdown. These phenomena severely compromise the effectiveness of conventional control surfaces. As a result, controlled flight at high angles of attack is not feasible for current aircraft configurations. Alternate means to augment the control of the vehicle at these flight regimes are therefore necessary. The present work investigates the augmentation of an aircraft flight control system by the injection of a thin sheet of air tangentially to the forebody of the vehicle. This method, known as Forebody Tangential Blowing (FTB), has been proposed as an effective means of increasing the controllability of aircraft at high angles of attack. The idea is based on the fact that a small amount of air is sufficient to change the separation lines on the forebody. As a consequence, the strength and position of the vortices are altered causing a change on the aerodynamic loads. Although a very effective actuator, forebody tangential blowing is also highly non-linear which makes its use for aircraft control very difficult. In this work, the feasibility of using FTB to control the roll-yaw motion of a wind tunnel model was demonstrated both through simulations and experimentally. The wind tunnel model used in the experiments consists of a wing-body configuration incorporating a delta wing with 70-degree sweep angle and a cone-cylinder fuselage. The model is equipped with forebody slots through which blowing is applied. There are no movable control surfaces, therefore blowing is the only form of actuation. Experiments were conducted at a nominal angle of attack of 45 degrees. A unique apparatus that constrains the model to two degrees-of-freedom, roll and yaw, was designed and built. The apparatus was used to conduct dynamic experiments which showed that the system was unstable, its natural motion divergent. A model for the unsteady aerodynamic loads was developed based on the basic physics of the flow and results from flow visualization experiments. Parameters of the aerodynamic model were identified from experimental data. The model was validated using data from dynamic experiments. The aerodynamic model completes the equations of motion of the system which were used in the design of control laws using blowing as the only actuator. The unsteady aerodynamic model was implemented as part of the real-time vehicle control system. A control strategy using asymmetric blowing was demonstrated experimentally. A discrete vortex method was developed to help understand the main physics of the flow. The method correctly captures the interactions between forebody and wing vortices. Moreover, the trends in static loads and flow structure are correctly represented. Flow visualization results revealed the vortical structure of the flow to be asymmetric even for symmetric flight conditions. The effects of blowing, and roll and yaw angles on the flow structure were determined. It is shown that superimposing symmetric and asymmetric blowing has a linearizing effect on the actuator characteristics. Transient responses of roll and yaw moments to step input blowing were characterized, and their differences were explained based on the physical mechanisms through which these loads are generated.

Pedreiro, Nelson↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗
NASA NTRSDate not supplied

CFD Simulation of an Inward-Turning Supersonic Inlet Unstart at Flight Mach Number 1.7

Interest in inward-turning supersonic inlets renewed during NASA’s Commercial Supersonic Technology Project, as a means of reducing disturbances caused by propulsion integration. Inward-turning designs have the potential to greatly reduce the strength of shock waves emanating outward from the nacelle, and therefore their contribution to the overall sonic boom signature. For top-mounted, or over-wing installations, disturbances to aircraft aerodynamics are reduced. The nature of the inward-turning geometry results in a boundary-layer thickening effect, opposite that of the familiar thinning of an outward-turning cone boundary-layer. The application of a small amount of boundary-layer bleed at the throat has been shown to alleviate the negative effects of the boundary-layer on subsonic diffuser performance. Boundary-layer bleed, however, comes at the expense of increased complexity, and external disturbances at bleed exits that would partially defeat the purpose of the inward-turning scheme. Vortex generators were proposed as an alternative to bleed, and an initial CFD study of their effectiveness in this application forms the basis of the present work. A source-term model in the Wind-US CFD code was used to model the effects of an array of vortex generators placed at the entrance to the subsonic diffuser. Baseline results with no vortex generators compared favorably to previously reported results with no bleed at back-pressures near the design point, but stability margin in the present case was limited by a curious non-linear phenomena similar to unstart in mixed-compression inlets. While some improvement was noted with vortex generators, their effectiveness was diminished by the limited stability margin. A time-accurate calculation was subsequently performed to better understand this process. Results of this calculation are presented and instantaneous solutions compared with steady-state, subcritical results previously reported. The stability margin limitation is attributed to the use of a constant-pressure outflow boundary condition, which differed from that used in the prior simulations.

Pre-Compressor Cooling, Water Injection, Hypersoni↗
NASA NTRSDate not supplied

Global Analysis of Nonlinear Second-Mode Development in a Mach-6 Boundary Layer From High-Speed Schlieren Data

The second-mode instability on a 7◦ half-angle sharp cone at Mach 6 is analyzed using high-speed calibrated schlieren imagery at a frame rate near the expected fundamental frequency. Experiments were conducted in the NASA Langley 20-Inch Mach 6 facility at unit Reynolds number between 6.56×106 and 9.71×106 m−1. Time-resolved pixel intensity signals throughout the boundary layer are reconstructed using spatially available data in the schlieren images to recover an effective sampling rate of over 10 MHz; these are then converted to quantitative density gradients using a thin-lens-based calibration technique. A global analysis is performed on the schlieren data to investigate the nonlinear growth of the second-mode fundamental and harmonic content. Point-wise measures of the auto-bicoherence are used to identify specific triadic interactions and the locations of their highest levels of quadratic phase coupling. Significant resonance interactions between the second-mode fundamental and harmonic instabilities were found along with interactions between these and the mean flow. Bispectral mode decomposition is employed to educe the flow structures associated with these interactions. A similar analysis is performed for the power spectrum, with power spectral densities computed for each pixel’s time-series and spectral proper orthogonal decomposition employed to derive the modal structure and energy of the flow at specific frequencies. Comparisons between the bispectral quantities and second-mode power show that nonlinear interactions, particularly resonance interactions, are closely correlated with space-time modulation of disturbances during the nonlinear stage of transition.

boundary layer transition↗