Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 613 records · Page 34

A simplified theory of oscillating aerofoils in transonic flow - Review and extension

A simplified theory of the dynamic motion of aerofoils of finite thickness in transonic flow is presented which excludes the effect of shock waves on the aerofoil itself and, thus, is restricted to free stream Mach numbers equal to unity or above. Numerical examples are analyzed for two-dimensional steady and unsteady (including transient) aerofoil motion and three-dimensional steady and unsteady flow over delta wings. The effects of flow separation and improvements in Bernoulli's equation and the surface boundary condition are briefly discussed.

Dowell, E. H.↗

Simultaneous detection of separation and transition in surface shear layers

Flush-mounted hot-film gages have proved effective in detecting boundary-layer transition and in measuring skin friction but with limited success in detecting laminar separation and reattachment. The development of multielement micro hot-film sensors, and the recent discovery of the phase reversal phenomena associated with low-frequency dynamic shear stress signals across regions of laminar separation and turbulent reattachment, have made it possible to simultaneously and unambiguously detect these surface shear layer characteristics. Experiments were conducted on different airfoils at speeds ranging from low subsonic to transonic speeds to establish the technique for incompressible and compressible flow applications. The multielement dynamic shear stress sensor technique was successfully used to detect laminar separation, turbulent reattachment, as well as, shock induced laminar and turbulent separation.

Mangalam, Siva M.↗

Magnetized supernova remnants with cosmic rays

The effects of interstellar magnetic fields and cosmic rays on the dynamics of an SNR expanding into a warm H I gas are examined. As long as the shock wave driven by the SN explosion propagates faster than 110 km/s, the vicinity of the shock front is fully ionized, and cosmic rays are well coupled to the thermal fluid. They are first accelerated at the adiabatic front, and further compressed in the postshock cooling zone. When the shock velocity drops below 110 km/s, ion-neutral collisions in the vicinity of the shock dissipate the waves which couple cosmic rays to the thermal gas, and impede cosmic-ray acceleration. It is found that magnetic and cosmic-ray pressures together dominate over thermal pressure away from the magnetic poles. As a result, most of the shell becomes considerably thicker, and the shock wave propagates somewhat faster than in the nonmagnetic case. At late times, the transverse mass motions which take place from the poles to the equator create H I holes at the polar caps. This theory leads to a simple interpretation of the 'barrel-shaped' distribution of radio emission observed in some SNRs.

Ferriere, Katia M.↗

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗

The Postshock Chemical Lifetimes of Outflow Tracers and a Possible New Mechanism to Produce Water Ice Mantles

We have used a coupled time-dependent chemical and dynamical model to investigate the lifetime of the chemical legacy in the wake of C-type shocks. We concentrate this study on the chemistry of H2O and O2, two molecules which are predicted to have abundances that are significantly affected in shock-heated gas. Two models are presented: (1) a three-stage model of preshock, shocked, and postshock gas; and (2) a Monte Carlo cloud simulation where we explore the effects of stochastic shock activity on molecular gas over a cloud lifetime. For both models we separately examine the pure gas-phase chemistry as well as the chemistry including the interactions of molecules with grain surfaces. In agreement with previous studies, we find that shock velocities in excess of 10 km/s are required to convert all of the oxygen not locked in CO into H2O before the gas has an opportunity to cool. For pure gas phase models the lifetime of the high water abundances, or "H2O legacy," in the postshock gas is approximately (4-7) x 10(exp 5) yr, independent of the gas density. A density dependence for the lifetime of H2O is found in gas-grain models as the water molecules deplete onto grains at the depletion timescale. Through the Monte Carlo cloud simulation we demonstrate that the time-average abundance of H2O, the weighted average of the amount of time gas spends in preshock, shock, and postshock stages, is a sensitive function of the frequency of shocks. Thus we predict that the abundance of H2O, and to a lesser extent O2, can be used to trace the history of shock activity in molecular gas. We use previous large-scale surveys of molecular outflows to constrain the frequency of 10 km/s shocks in regions with varying star formation properties and discuss the observations required to test these results. We discuss the postshock lifetimes for other possible outflow tracers (e.g., SiO and CH3OH) and show that the differences between the lifetimes for various tracers can produce potentially observable chemical variations between younger and older outflows. For gas-grain models we find that the abundance of water-ice on grain surfaces can be quite large and is comparable to that observed in molecular clouds. This offers a possible alternative method to create water mantles without resorting to grain surface chemistry: gas heating and chemical modification due to a C-type shock and subsequent depletion of the gas-phase species onto grain mantles.

Bergin, Edwin A.↗

NASA Ultra-Sensitive Miniature Accelerometer

Using micro-machined silicon technology, an ultra-sensitive miniature acce.,rometer can be constructed which meets the requirements for microgravity experiments in the space environment.Such an accelerometer will have a full scale sensitivity of 1C2 g a resolution of lC8 g, low cross axis sensitivity, and low temperature sensitivity. Mass of the device is approximately five grams and its footprint is 2 cm x 2 cm. Innovative features of the accelerometer, which are patented, are: electrostatic caging to withstand handling shock up to 150 g, in-situ calibration, in situ performance characterization, and both static and dynamic compensation. The transducer operates on a force balance principle wherein the displacement of the proof mass is monitored by measuring tunneling electron current flow between a conductive tip, and a fixed platen. The four major parts of the accelerometer are tip die, incorporating the tunneling tip and four field plates for controlling pitch and roll of the proof mass; two proof mass dies, attached to the surrounding frame by sets of four leg" springs; and a force plate die. The four parts are fuse-bonded into a complete assembly. External electrical connections are made at bond pads on the front surface of the force plate die. Materials and processes used in the construction of the transducer are compatible with volume production.

Zavracky, Paul M.↗

Scramjet Engine Flowpath Development for the Hyper-X Mach 10 Flight Test

The third flight of the Hyper-X Research Vehicle at nearly Mach 10 in November 2004 demonstrated the potential for airbreathing propulsion at hypersonic speeds. The flowpath of the scramjet engine used in the Hyper-X Research Vehicle was developed using computational fluid dynamic based design method in combination with a systematic series of ground tests in the NASA HyPulse Shock Tunnel at conditions duplicating Mach 10 flight. The test model was a true scale height and length of the internal flight engine flowpath that was partial width and truncated fore and aft. Tests were conducted using a mixture of hydrogen augmented with two-percent silane. Silane gas served as an ignition aid during the short duration of the pulse facility tests. The engine inflow values of pressure, temperature, and Mach number were parametrically varied during the ground test entries to broaden the database over the expected uncertainty in the flight conditions. A review of the ground test technique and comparisons of the ground test pressures along with selected data from the third flight are presented.

Rogers, R. C.↗

Optical Diagnostic Imaging of Multi-Rocket Plume-Induced Base Flow Environments

Non-intrusive optical diagnostic imaging was used for the first time to visualize multi-rocket plume-induced reacting base flows to simulate launch vehicle ascent from sea-level to 250,000 ft. In particular, planar laser induced florescence (PLIF) and infrared (IR) imaging were implemented for the first time to visualize and quantify base flow and rocket plume environments from sub-scale, short-duration propulsion models within a shock tunnel facility. This report discusses the successful imaging diagnostic methods for capturing base flow features and dynamics as a function of altitude. Important base flow and plume features were captured with PLIF and IR diagnostics to develop a conceptual base flow physics model. This imaging data specifically provides insight into the Space Launch System vehicle core-stage and Exploration Upper Stage base environments and further validates short-duration ground test techniques and computational modeling.

Mehta, Manish↗

Generation of a Reduced-Order Model of an Unmanned Combat Air Vehicle Using Indicial Response Functions

Within the NATO STO AVT-251 Task Group, a generic Unmanned Combat Air Vehicle (UCAV) planform was redesigned based on requirements derived from parts of the flight envelope of the defined mission. Because of the lambda-shape planform and associated flow phenomena including shocks and vortices, the aerodynamic design process relies heavily on high-fidelity predictions from computational fluid dynamics (CFD). These simulations enable accurate flight dynamics predictions, allowing for the identification of any potential performance issues early in the design process. Potentially, a complete flight dynamics model can be derived from CFD calculations provided that the underlying reduced order model (ROM) captures the underlying physics contained within the CFD results. This work investigates the creation of a nonlinear ROM using indicial response functions. The response functions are obtained using a grid motion approach that separates the effects of angle of attack and pitch rate. This approach is demonstrated using three different CFD codes from various organizations: ENSOLV at Royal Netherlands Aerospace Centre, USM3D at NASA Langley Research Center, and Kestrel at the United States Air Force Academy. ROM predictions were generated for a manoeuver resembling a low-speed pull-up. The predictions are found to be sensitive to the quality of the steady-state solutions from which step-response calculations were started, as well as the convergence of the step responses themselves. Nevertheless, the indicial response method is shown to provide aerodynamic predictions with acceptable accuracy including transient effects, and as such is a promising method for computationally efficient flight dynamics predictions.

Michel P.C. van Rooij↗

Effect of anisotropic mesh adaptation on surface pressure predictions for atmospheric entry simulations.

Recently, a metric-based mesh adaptation capability has been developed for the advanced compressible flow solver US3D. Metric-based mesh adaptation is a solution informed mesh adaptation strategy that typically improves the resolution of the mesh along anisotropic flow features, like strong shock waves or shear layers. An anisotropic error indicator is reconstructed based on a provided Computational Fluid Dynamics (CFD) solution which prescribes the locally desired element size and orientation. The benefits of this mesh adaptation strategy is twofold: first, it simplifies the manual meshing labor that typically needs to be carried out by the CFD engineer because the mesh adaptation is informed by the provided solution. Second, the computational cost is kept to a minimum while improving the resolution of the simulation. The aim here is to demonstrate how anisotropic mesh adaptation can be employed to improve the predictions of surface pressure for atmospheric entry vehicles. In this work, we compare CFD simulations that use this newly developed anisotropic mesh adaptation capability for US3D with pressure data that was collected by Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) during the atmospheric entry phase of Mars2020.

D. Ekelschot↗

Effect of Anisotropic Mesh Adaptation on Surface Pressure Predictions for Atmospheric Entry Simulations

A metric-based mesh adaptation capability has been developed for the advanced compressible flow solver US3D. Metric-based mesh adaptation is a solution informed mesh adaptation strategy that improves the resolution of the mesh along anisotropic flow features, like strong shock waves or shear layers. An anisotropic error indicator is reconstructed based on a provided Computational Fluid Dynamics (CFD) solution which prescribes the locally desired element size and orientation. The benefits of this mesh adaptation strategy is twofold: first, it simplifies the manual meshing labor that typically needs to be carried out by the CFD engineer because the mesh adaptation is informed by the provided solution. Second, the computational cost is kept to a minimum while improving the resolution of the simulation. The aim here is to demonstrate how anisotropic mesh adaptation can be employed to improve the predictions of surface pressure for atmospheric entry vehicles. In this work, we compare CFD simulations that use this newly developed anisotropic mesh adaptation capability for US3D with pressure data that was collected by Mars Entry, Descent, and Landing Instrumentation (MEDLI2) suite.

Dirk Ekelschot↗

On Multi-Dimensional Unstructured Mesh Adaption

Anisotropic unstructured mesh adaption is developed for a truly multi-dimensional upwind fluctuation splitting scheme, as applied to scalar advection-diffusion. The adaption is performed locally using edge swapping, point insertion/deletion, and nodal displacements. Comparisons are made versus the current state of the art for aggressive anisotropic unstructured adaption, which is based on a posteriori error estimates. Demonstration of both schemes to model problems, with features representative of compressible gas dynamics, show the present method to be superior to the a posteriori adaption for linear advection. The performance of the two methods is more similar when applied to nonlinear advection, with a difference in the treatment of shocks. The a posteriori adaption can excessively cluster points to a shock, while the present multi-dimensional scheme tends to merely align with a shock, using fewer nodes. As a consequence of this alignment tendency, an implementation of eigenvalue limiting for the suppression of expansion shocks is developed for the multi-dimensional distribution scheme. The differences in the treatment of shocks by the adaption schemes, along with the inherently low levels of artificial dissipation in the fluctuation splitting solver, suggest the present method is a strong candidate for applications to compressible gas dynamics.

Wood, William A.↗

Study of shock-induced combustion using an implicit TVD scheme

The supersonic combustion flowfields associated with various hypersonic propulsion systems, such as the ram accelerator, the oblique detonation wave engine, and the scramjet, are being investigated using a new computational fluid dynamics (CFD) code. The code solves the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. It employs an iterative method and a second order differencing scheme to improve computational efficiency. The code is currently being applied to study shock wave/boundary layer interactions in premixed combustible gases, and to investigate the ram accelerator concept. Results obtained for a ram accelerator configuration indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outward and downstream. The combustion process creates a high pressure region over the back of the projectile resulting in a net positive thrust forward.

Yungster, Shayne↗

Rarefied gas dynamics; Proceedings of the Eighth International Symposium, Stanford University, Stanford, Calif., July 10-14, 1972

Recent studies in the theory of rarefied gases (kinetic theory, shock structure, molecular forces and interactions, evaporation and condensation, expansions and jets, flow past bodies, gas-solid interactions) are presented along with descriptions of numerous experiments in these areas. Some of the topics covered include: singular limits of solutions of Boltzmann's equations, the Boltzmann and Klimontovich formalisms with reference to turbulent resolution, higher-order kinetic theory for the Lorentz gas, intermolecular potentials from shock structure experiments, cooling of particles in aerosol beams, hypersonic-subsonic transition in spherically symmetric expanding flows, absolute velocity measurements in a rarefied gas flow by an ion time-of-flight technique, the product distribution of gas-solid reactions from absolute rate theory, and relaxation of a gas with internal degrees of freedom. Individual items are announced in this issue.

Karamcheti, K.↗

Dynamics of the outgoing turbulent boundary layer in a Mach 5 unswept compression ramp interaction

Wall pressure fluctuations have been measured under the unsteady separation shock and on the ramp face in an unswept Mach 5 compression ramp interaction. The freestream Reynolds number was 51.0 x 10 to the 6th/m, and the incoming turbulent boundary layer developed on the tunnel floor under approximately adiabatic wall temperature conditions. Standard data-acquisition methods, as well as real-time and posttest conditional sampling techniques were used. The results show that the mean and rms pressure levels are strong functions of separation shock position. At all stations on the ramp, from the corner to where the pressure reaches the theoretical inviscid value, the pressure signals have two dominant components: a low frequency component characteristic of the global unsteadiness, which correlates with the separation shock motion, and a higher frequency component associated with turbulence. The former is the major contributor to the overall signal variance.

Gramann, Richard A.↗

Algorithm For Hypersonic Flow In Chemical Equilibrium

Implicit, finite-difference, shock-capturing algorithm calculates inviscid, hypersonic flows in chemical equilibrium. Implicit formulation chosen because overcomes limitation on mathematical stability encountered in explicit formulations. For dynamical portion of problem, Euler equations written in conservation-law form in Cartesian coordinate system for two-dimensional or axisymmetric flow. For chemical portion of problem, equilibrium state of gas at each point in computational grid determined by minimizing local Gibbs free energy, subject to local conservation of molecules, atoms, ions, and total enthalpy. Major advantage: resulting algorithm naturally stable and captures strong shocks without help of artificial-dissipation terms to damp out spurious numerical oscillations.

Palmer, Grant↗

The magnetosheath electron population at lunar distance - General features

A study was made of the electron population in the earth's magnetosheath at lunar distance. The data used were collected by the Apollo 14 charged particle lunar environment experiment (CPLEE) during four inbound (dusk) and three outbound (dawn) passages of the moon through the magnetosheath. The magnetotail has a diameter of 52 earth radii, while the bow-shock cross section is about 91 earth radii. The average boundary locations computed from the complete data set are consistent with the prediction of fluid dynamics. The electron characteristics for the two least-disturbed passages are presented in detail. An examination of the energy spectra shows that a high-energy (200-2000 eV) tail is superimposed on the expected low-energy (40-200 eV) magnetosheath distribution. It is argued that the high-energy magnetosheath electron population originates at the bow shock, rather than from the plasma sheet.

Reiff, P. H.↗

Compressibility effects on dynamic stall of airfoils undergoing rapid transient pitching motion

The research was carried out in the Compressible Dynamic Stall Facility, CDSF, at the Fluid Mechanics Laboratory (FML) of NASA Ames Research Center. The facility can produce realistic nondimensional pitch rates experienced by fighter aircraft, which on model scale could be as high as 3600/sec. Nonintrusive optical techniques were used for the measurements. The highlight of the effort was the development of a new real time interferometry method known as Point Diffraction Interferometry - PDI, for use in unsteady separated flows. This can yield instantaneous flow density information (and hence pressure distributions in isentropic flows) over the airfoil. A key finding is that the dynamic stall vortex forms just as the airfoil leading edge separation bubble opens-up. A major result is the observation and quantification of multiple shocks over the airfoil near the leading edge. A quantitative analysis of the PDI images shows that pitching airfoils produce larger suction peaks than steady airfoils at the same Mach number prior to stall. The peak suction level reached just before stall develops is the same at all unsteady rates and decreases with increase in Mach number. The suction is lost once the dynamic stall vortex or vortical structure begins to convect. Based on the knowledge gained from this preliminary analysis of the data, efforts to control dynamic stall were initiated. The focus of this work was to arrive at a dynamically changing leading edge shape that produces only 'acceptable' airfoil pressure distributions over a large angle of attack range.

Chandrasekhara, M. S.↗