Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock capturing”

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 577 records · Page 32

A point implicit unstructured grid solver for the Euler and Navier-Stokes equations

An upwind finite element technique that uses cell centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured triangular grids. A structured grid of quadrilaterals is laid out near the body surface. For inviscid flows the method is stable at Courant numbers of over 100,000. A first order basic scheme and a higher order flux corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view of simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively generated grids. The adaptive mesh generator and the upwind Navier-Stokes solver are combined in a set of programs called LARCNESS, an acronym for Langley Adaptive Remeshing Code and Navier-Stokes Solver.

Thareja, Rajiv R.↗

Euler and Navier-Stokes solutions for hypersonic flows

An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.

Thareja, Rajiv R.↗

Compressible 2D Flow Field Interaction of Two Contra-Rotating Blades

The blades of coaxial, contra-rotating rotor systems cross each other in close proximity and at high relative speeds. This crossing event is a potential source of noise and severe blade loads. Effects of compressibility can aggravate the interaction and significantly alter the pressure field signature and phase relationships. A 2-D analysis of this phenomenon is performed by simulating two airfoils passing each other at specified speeds and vertical separation distances. Several test cases spanning a relevant range of Reynolds numbers, angles of attack, and relative Mach number are considered. The Mach number is varied to simulate the radial variation of velocity from the root to tip of a rotor blade to capture the pressure signature, lift, and drag of the airfoils. The velocity and pressure distributions on the airfoils, and in the space between the airfoils are computed before, at, and after airfoil crossing. The variations of lift and drag coefficients through the interaction are captured. The upper airfoil experiences an increase in lift followed by a very sharp drop in lift during the interaction. When relative Mach numbers are transonic, the region of interaction is greatly extended, with shock interactions occurring. The results show the complex nature of the aerodynamic and fluid dynamic impulses generated by blade-blade interactions, with implications to aeroelastic loads and aeroacoustic sources.

Barbely, Natasha L.↗

Coherence Analysis of the Space Launch System using Unsteady Pressure Sensitive Paint

Transonic buffet during atmospheric ascent is a major source of unsteady loading on launch vehicles that, in the past, has led to structural failures. Thus, determining accurate buffet forcing functions (BFFs) to properly predict the vehicle response to buffet is of vital importance in launch vehicle design. The state of the art for obtaining the BFFs relies on wind-tunnel tests where the fluctuating pressures are measured by pressure transducers (PTs) at sparse locations on a rigid, geometrically-scaled buffet model. To compute the BFFs, the outer mold line (OML) of the vehicle is mapped onto contiguous panels centered at the location of the PTs and the measured fluctuating pressures are integrated over the panels’ areas. To mitigate conservatism due to the assumption that the measured buffet pressures act in phase across each panel, coherence factors are applied, effectively reducing the integration areas and, therefore, the buffet forces. For coherence factors to provide the proper level of BFF attenuation, accurate knowledge of the spatial and temporal correlation of the buffet pressures is paramount. Unfortunately, even with hundreds of unsteady pressure transducers instrumenting the models, compromises must be made on the spatial resolution of the PTs. Typically, the PT layout aims at resolving the pressure correlation along the longitudinal axis of the vehicle, assuming that coherent structures propagate mostly in the longitudinal direction. As a result, the distribution of azimuthal/cross-stream correlation is not well known and its impact on the estimated BFFs is often neglected. To fill this and other knowledge gaps, extremely high spatial-resolution uPSP data were collected for three different production-design configurations of the Space Launch System in the NASA Ames Research Center Unitary Plan Wind Tunnel 11-Foot Transonic Wind Tunnel. Specifically, two cargo configurations, the Block 1 and Block 1B, and one crew configuration, the Block 1B crew, were surveyed at resolutions ranging from 600,000 to over 1 million uPSP measurement locations. To shed light on the temporal and spectral behavior of are presented. Several OML regions and flow features of interest are investigated, from theexpansion/shock on the Orion Multi-Purpose Crew Vehicle, to the terminal shock environment on the core stage, and the Strouhal shedding behind the boosters forward attach. The sensitivity of these environments to the vehicle attitude is examined. Furthermore, for selected panels,coherence factors that accurately capture the azimuthal coherence distribution of the buffetpressures are derived and their impact on the estimated BFFs is discussed. Finally, distributions of the local convection velocity and cross spectrum phase are presented.

buffet↗

Numerical simulation of helicopter multi-bladed rotor flow

A numerical method for solving the Euler equations for multiblade rotors has been developed. The computed results, including one nonlifting forward-flight and several lifting hovering cases, are compared with experimental data. For forward-flight case, the unsteady growth and decay of the shock agree well with experimental results. For hovering case, at low and moderate transonic tip speed, the results show good agreement with experimental data in the tip region. The near-wake effects can be captured without wake modeling.

Chen, C. L.↗

Laser absorption flow diagnostics for use in hypersonic ground-based and flight experiments

Consideration is given to laboratory ground-based scale systems for measuring OH and NO and diode laser systems for measuring O2 and H2O. Several diagnostic instruments which are currently being used in a scramjet combustor test rig of the Ames 16-inch shock tunnel are described. Data on mole fraction, temperature, density, pressure, and velocity is obtained using specific data reduction techniques. These measurements provide critical data for evaluation of the propulsion system performance which include mass capture for quantification of fuel equivalence ratio (O2), flow contamination and emission levels (NO), combustion progress (OH), and injector mixing and combustion efficiency (H2O).

Cavolowsky, John A.↗

Equilibrium radiative heating tables for aerobraking in the Martian atmosphere

Studies currently underway for Mars missions often envision the use of aerobraking for orbital capture at Mars. These missions generally involve blunt-nosed vehicles to dissipate the excess energy of the interplanetary transfer. Radiative heating may be of importance in these blunt-body flows because of the highly energetic shock layer around the blunt nose. In addition, the Martian atmosphere contains CO2, whose dissociation products are known to include strong radiators. An inviscid, equilibrium, stagnation point, radiation-coupled flow-field code has been developed for investigating blunt-body atmospheric entry. The method has been compared with ground-based and flight data for air, and reasonable agreement has been found. In the present work, the method was applied to a matrix of conditions in the Martian atmosphere. These conditions encompass most trajectories of interest for Mars exploration spacecraft. The predicted equilibrium radiative heating to the stagnation point of the vehicle is presented.

Hartung, Lin C.↗

Isotopic Evidence for a Martian Regolith Component in Martian Meteorites

Noble gas measurements in gas-rich impact-melt (GRIM) glasses in EET79001 shergottite showed that their elemental and isotopic composition is similar to that of the Martian atmosphere [1-3]. The GRIM glasses contain large amounts of Martian atmospheric gases. Those measurements further suggested that the Kr isotopic composition of Martian atmosphere is approximately similar to that of solar Kr. The (80)Kr(sub n) - (80)Kr(sub M) mixing ratio in the Martian atmosphere reported here is approximately 3%. These neutron-capture reactions presumably occurred in the glass-precursor regolith materials containing Sm- and Br- bearing mineral phases near the EET79001/ Shergotty sites on Mars. The irradiated materials were mobilized into host rock voids either during shock-melting or possibly by earlier aeolian / fluvial activity.

Rao, M. N.↗

Neutrinos from type II supernovae - The first 100 milliseconds

The collapse of a 1.17 solar mass iron core is numerically followed through infall to 100 ms past core bounce, and the emergent neutrino spectra during each phase are highlighted. It is found that, even with fairly optimistic conditions for producing a strong, sustained core-bounce shock wave, the prompt shock stalls within 9 ms of core bounce at a radius of less than 250 km. It appears that a radical change in the character of the progenitor core or in our understanding of the relevant physics of stellar collapse is needed before the direct mechanism for type II supernovae can become viable. Expanding the number of neutrino types from one to six magnifies the debilitating effect of neutrino loss on shock propagation. At shock breakout, prompt bursts of all neutrino types are observed. The luminosities of the nonelectron types show a sudden turn-on in luminosity while that of the electron neutrinos steadily increases throughout infall as a result of accelerating electron capture.

Myra, Eric S.↗

A finite-volume Euler solver for computing rotary-wing aerodynamics on unstructured meshes

An unstructured-grid solver for the unsteady Euler equations has been developed for predicting the aerodynamics of helicopter rotor blades. This flow solver is a finite-volume scheme that computes flow quantities at the vertices of the mesh. Special treatments are used for the flux differencing and boundary conditions in order to compute rotary-wing flowfields, and these are detailed in the paper. The unstructured-grid solver permits adaptive grid refinement in order to improve the resolution of flow features such as shocks, rotor wakes and acoustic waves. These capabilities are demonstrated in the paper. Example calculations are presented for two hovering rotors. In both cases, adaptive-grid refinement is used to resolve high gradients near the rotor surface and also to capture the vortical regions in the rotor wake. The computed results show good agreement with experimental results for surface airloads and wake geometry.

Strawn, Roger C.↗

The interaction of interstellar gas with stellar density wave packets

We have numerically studied the interaction between interstellar gas and a gravitational perturbation produced by a fluctuation in the stellar mass distribution. For a simple equation of state, shock waves develop even if the gravitational perturbation is aperiodic and transient. In a medium which can change thermal phase, shocks do not form and the gas makes a transition directly to the cold, dense phase. The cold material has an anomalous velocity imparted to it by the gravitational perturbation and it becomes an accretion front moving through the hot gas and capturing it. The accretion fronts trigger thermal phase change on a large scale and can reproduce the extent of observed cold clouds. Furthermore their high mass is favorable to gravitational instability and subsequent star formation.

Baker, P. L.↗

Inlet performance of the integrated Langley scramjet module

The inlet concept for the Langley Scramjet Module has been developed and proven in Langley wind tunnels over a Mach number range from 2.3 to 6.0 (flight simulation of Mach 2.6 to 7.6). This modular engine concept is designed to integrate with the airframe, which results in precompression of the engine airflow by the vehicle bow shock and additional expansion of the nozzle exhaust gas by the afterbody of the vehicle. With these integration advantages, the inlet can be designed with modest contraction ratios and fixed geometry. Also, the module nozzle exit area can be equal to the capture area, which permits the cowl to be alined with the local flow producing minimum external drag. The inlet leading edges and planar compression surfaces are swept at 48 deg, which provides spillage at low Mach numbers for starting and which reduces the pressure gradient on the top surface to permit ingestion of the vehicle forebody boundary layer into the inlet without separating. Three fuel injection struts provide for the use of a short combustor having low internal cooling requirements. Schedules for mass capture ratio, contraction ratio, and total pressure recovery are well within the acceptable range for a good scramjet propulsion device. The fixed geometry, minimum external drag design has proven to be a practical, high-performance inlet concept.

Trexler, C. A.↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗

A burst from a thermonuclear runaway on an ONeMg white dwarf

Studies which examine the consequences of accretion, at rates of 10(exp -9) solar mass/yr and 10(exp -10) solar mass/yr, onto an ONeMg white dwarf with a mass of 1.35 solar masses are performed. In these studies, a Lagrangian, hydrodynamic, one-dimensional computer code was used. The code now includes a network with 89 nuclei up to Ca-40, elemental diffusion, new opacities, and new equation of state. The initial abundance distribution corresponded to a mixture that was enriched to either 25, 50, or 75 percent in products of carbon burning. The remaining material in each case is assumed to have a solar composition. The evolution of the thermonuclear runaway in the 1.35 solar mass white dwarf, with M = 10(exp -9) solar mass, produced peak temperatures in the shell source exceeding 300 million degrees. The sequence produced significant amounts of Na-22 from proton captures onto Ne-20 and significant amounts of Al-26 from proton captures on Mg-24. This sequence ejected 5.2 x 10(exp -6) solar mass moving with speeds from approximately 100 km/s to 2300 km/s. When the mass accretion rate was decreased to 10(exp -10) solar mass, the resulting thermonuclear runaway produced a shock that moved through the outer envelope of the white dwarf and raised the surface luminosity to L greater than 10(exp 7) solar luminosity and the effective temperature to values exceeding 10(exp 7) K. The interaction of the material expanding from off of the white dwarf with the accretion disk should produce a burst of gamma-rays.

Starrfield, S.↗

Design of Three-Dimensional Hypersonic Inlets with Rectangular to Elliptical Shape Transition

A methodology has been devised for the design of three-dimensional hypersonic inlets which include a rectangular to elliptical shape transition. This methodology makes extensive use of inviscid streamtracing techniques to generate a smooth shape transition from a rectangular-like capture to an elliptical throat. Highly swept leading edges and a significantly notched cowl enable use of these inlets in fixed geometry configurations. The design procedure includes a three dimensional displacement thickness calculation and uses established correlations to check for boundary layer separation due to shock wave interactions. Complete details of the design procedure are presented and the characteristics of a modular inlet with rectangular to elliptical shape transition and a design point of Mach 7.1 are examined. Comparison with a classical two-dimensional inlet optimized for maximum total pressure recovery indicates that this three-dimensional inlet demonstrates good performance even well below its design point.

Smart, M. K.↗

Investigation of a mixed compression axisymmetric inlet at Mach number 5.3

The hypersonic diffuser portion of an uncooled high performance mixed compression, axisymmetric inlet suitable for subsonic burning engines was designed and tested. Performance of a model with a 25.4-cm capture diameter was measured in a wind tunnel and the results were compared with theoretical predictions calculated by a comprehensive computer program. All tests were conducted at a Mach number of 5.3 at a total temperature of 667 K and a total pressure of 11.57 atm. The angle of attack ranged from 0 to + or - 3 deg. Performance at angle of attack remained high. Reasonably high performance in the throat (maximum throat pitot-pressure recovery of 77 percent and an average value of 58 percent) was obtained at 0 deg angle of attack with relatively large amounts of boundary-layer bleed (11 to 22 percent of the capture mass flow). The computer program used in the design of this inlet is considered marginally adequate for predicting hypersonic inlet flow fields. Although the program as it now exists is very useful, an improved computer program that more accurately predicts the boundary layer and the shock-wave-boundary-layer interaction and accounts for boundary-layer bleed should be developed for reliability predicting hypersonic inlet flow fields.

Latham, E. A.↗

A new technique for ground simulation of hypervelocity debris

A series of hypervelocity damage experiments were preformed on spacecraft materials. These experiments employed a technique which accelerates micro flyer plates simulating space debris traveling at 3 to 8 km/sec. The apparatus used to propel the micro flyer plates was compact and fit well into a space environmental chamber equipped with instrumentation capable of analyzing the vapor ejected from the sample. Mechanical damage to the sample was also characterized using optical and scanning electron microscpopy. Data for this work was obtained from hypervelocity impacts on a polysulfone resin and a graphite polysulfone composite. Polysulfone was selected because it was flown on the Long Duration Exposure Facility (LDEF) which spent several years in low earth orbit (LEO). Chemistry of the vapor produced by the impact was analyzed with a time of flight mass spectrometer, (TOFMS). This represents the first time that ejected vapors from hypervelocity collisions were trapped and analyzed with a mass spectrometer. With this approach we are able to study changes in the vapor chemistry as a function of time after impact, obtain a velocity measurement of the vapor, and estimate a temperature of the surface at time of impact using dynamic gas equations. Samples of the vapor plume may be captured and examined by transmission electron microscopy. Studies were also conducted to determine mechanical damage to a graphite polysulfone composite and a polysulfone resin. Impact craters were examined under optical and scanning electron microscopes. The collision craters in the matrix were typical of those shown in conventional shock experiments. However, the hypervelocity collisions with the graphite polysulfone composite were remarkably different than those with the resin.

Roybal, R.↗

The intact capture of hypervelocity dust particles using underdense foams

The impact of a hypervelocity projectile (greater than 3 km/s) is a process that subjects both the impactor and the impacted material to a large transient pressure distribution. The resultant stresses cause a large degree of fragmentation, melting, vaporization, and ionization (for normal densities). The pressure regime magnitude, however, is directly related to the density relationship between the projectile and target materials. As a consequence, a high-density impactor on a low-density target will experience the lowest level of damage. Historically, there have been three different approaches toward achieving the lowest possible target density. The first employs a projectile impinging on a foil or film of moderate density, but whose thickness is much less than the particle diameter. This results in the particle experiencing a pressure transient with both a short duration and a greatly reduced destructive effect. A succession of these films, spaced to allow nondestructive energy dissipation between impacts, will reduce the impactor's kinetic energy without allowing its internal energy to rise to the point where destruction of the projectile mass will occur. An added advantage to this method is that it yields the possibility of regions within the captured particle where a minimum of thermal modification has taken place. Polymer foams have been employed as the primary method of capturing particles with minimum degradation. The manufacture of extremely low bulk density materials is usually achieved by the introduction of voids into the material base. It must be noted, however, that a foam structure only has a true bulk density of the mixture at sizes much larger than the cell size, since for impact processes this is of paramount importance. The scale at which the bulk density must still be close to that of the mixture is approximately equal to the impactor. When this density criterion is met, shock pressures during impact are minimized, which in turn maximizes the probability of survival for the impacting particle. The primary objectives of the experiment are to (1) Examine the morphology of primary and secondary hypervelocity impact craters. Primary attention will be paid to craters caused by ejecta during hypervelocity impacts of different substrates. (2) Determine the size distribution of ejecta by means of witness plates and collect fragments of ejecta from craters by means of momentum-sensitive mcropore foam. (3) Assess the directionality of the flux by means of penetration-hole alignment of thin films placed above the cells. (4) Capture intact the particles that perforated the thin film and entered the cell. Capture media consisted of both previously flight-tested micropore foams and aerogel. The foams had different latent heats of fusion and, accordingly, will capture particles over a range of momenta. Aerogel was incorporated into the cells to determine the minimum diameter than can be captured intact.

Maag, Carl R.↗