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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

A time-split finite-volume algorithm for three-dimensional flow-field simulation

A general finite-volume algorithm is developed for solving three-dimensional, time-dependent, compressible Navier-Stokes equations for high Reynolds number flows over an arbitrary geometry. This algorithm adapts MacCormack's (1982) explicit-implicit scheme to a time-split, three-dimensional finite-volume concept in a general coordinate system. It is shown that the thin-layer approximation in all three spatial directions significantly reduces the evaluation of viscous terms and allows the algorithm to solve more complicated geometries with all boundaries in two or all three directions. The calculated results using this method are found to be in good agreement with the experimental measurements of a blunt-fin induced shock wave and boundary-layer interaction problems. Observations of the existence of peak pressure, primary horseshoe and secondary vortices, and reversed supersonic zones show that computational fluid dynamics can effectively supplement the wind tunnel tests for aerodynamic design as well as for understanding basic fluid dynamics.

Hung, C. M.↗
NASA NTRS2003

Control of Interacting Vortex Flows at Subsonic and Transonic Speeds Using Passive Porosity

A wind tunnel experiment was conducted in the NASA Langley Research Center (LaRC) 8-foot Transonic Pressure Tunnel (TPT) to determine the effects of passive surface porosity on vortex flow interactions about a general research fighter configuration at subsonic and transonic speeds. Flow- through porosity was applied to a wind leading-edge extension (LEX) mounted to a 65 deg cropped delta wind model to promote large nose-down pitching moment increments at high angles of attack. Porosity decreased the vorticity shed from the LEX, which weakened the LEX vortex and altered the global interactions of the LEX and wing vortices at high angles of attack. Six-component forces and moments and wing upper surface static pressure distributions were obtained at free- stream Mach numbers of 0.50, 0.85, and 1.20, Reynolds number of 2.5(10(exp-6) per foot, angles of attack up to 30 deg and angles of sideslip to plus or minus 8 deg. The off-surface flow field was visualized in selected cross-planes using a laser vapor screen flow visualization technique. Test data were obtained with a centerline vertical tail and with alternate twin, wing-mounted vertical fins having 0 deg and 30 deg cant angles. In addition, the porosity of the LEX was compartmentalized to determine the sensitivity of the vortex- dominated aerodynamics to the location and level of porosity applied to the LEX.

NATO FURNISHED↗
NASA NTRS2004

Accretion and Heating of Particles by Supersonic Planetesimals

Recently, nebular shock waves have become one of the leading candidates for explaining the presence of chondrules in primitive meteorites. While shocks have been shown to be capable of explaining many of the features of chondrules, a major problem with the theory is that the source of the shocks remains unidentified. Among the suggested sources of the shocks are bow shocks created by supersonic planetesimals in the nebula. Recently, we studied the structure of the shocks that would form around such supersonic planetesimals. We found that particles that encountered the shocks at distances greater than 2 planetesimal radii from the planetesimal would cool too quickly to form the textures observed in chondrules. While the region of the shock far away from the planetesimal may not allow chondrules to form, the region closer in has not been studied in detail. In this work we consider the dynamical and thermal evolution of particles that encounter supersonic planetesimals in this region.

Ciesla, F. J.↗
NASA NTRSDate not supplied

3D Radiative MHD Modeling of the Solar Atmospheric Dynamics and Structure

Dramatic dynamical phenomena accompanied by strong thermodynamic and magnetic structuring are the primary drivers of great interest in studying the solar atmosphere with high spatial and temporal resolutions. Using current computational capabilities, it has become possible to model the magnetized solar plasma in different regimes with a high degree of realism. To study the fine structuring of the solar atmosphere and dynamics, we use 3D MHD radiative models covering all layers from the upper convection zone to the corona. Realistic 3D radiative MHD modeling of the solar magnetoconvection and atmosphere allows us to generate synthetic observables that directly link the physical properties of the solar plasma to spectroscopic observables. In the presentation, we discuss qualitative and quantitative changes of the atmospheric structure and dynamics at different layers of the solar atmosphere, properties of acoustic and surface gravity waves, sources of local heating in the chromosphere-corona transition region, formation of shocks, and high-frequency oscillations in the corona, as well as manifestation of these phenomena in the modeled observables.

SMD↗
NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗
NASA NTRS1976

Flux-corrected transport techniques for transient calculations of strongly shocked flows

New flux-corrected transport algorithms are described for solving generalized continuity equations. These techniques were developed by requiring that the finite difference formulas used ensure positivity for an initially positive convected quantity. Thus FCT is particularly valuable for fluid-like problems with strong gradients or shocks. Repeated application of the same subroutine to mass, momentum, and energy conservation equations gives a simple solution of the coupled time-dependent equations of ideal compressible fluid dynamics without introducing an artificial viscosity. FCT algorithms span Eulerian, sliding-rezone, and Lagrangian finite difference grids in several coordinate systems. The latest FCT techniques are fully vectorized for parallel/pipeline processing.

Boris, J. P.↗
NASA NTRS1979

Interplanetary magnetic fields

The paper reviews the following studies on interplanetary magnetic fields: (1) the relation of the interplanetary field to the sun and to solar-wind streams; (2) radial gradients in the solar wind and interplanetary field; (3) shocks, hydromagnetic waves and discontinuities; (4) variations of interplanetary field parameters over the solar cycle; (5) the relation between the interplanetary field and the terrestrial magnetosphere; (6) the internal dynamics of the magnetic field and solar wind plasma; and (7) the transport, acceleration, and modulation of energetic particles, including cosmic rays. It is noted that further studies will relate to the complex relations between the interplanetary field and solar structure.

Smith, E. J.↗
NASA NTRS1989

Development and validation of CNS (compressible Navier-Stokes) for hypersonic external flows

CNS, a new computational fluid dynamics procedure, has been developed to aid in hypersonic vehicle design. The code can be used to model the entire external flow around hypersonic vehicle shapes, from the captured shock at the nose to the beginning of the wake. Unlike space-marching codes, the technique allows axially separated flow regions to be modeled. Validation trials using sphere-cone data reveal good solution accuracy for the surface pressure and flowfield temperature.

Flores, Jolen↗
NASA NTRS1989

International Conference on Numerical Methods in Fluid Dynamics, 11th, Williamsburg, VA, June 27-July 1, 1988, Proceedings

Recent advances in computational fluid dynamics (CFD) are discussed in reviews and reports. Topics addressed include CFD models in plasma dynamics, parallel computation for simulation studies, CFD for hypersonic airbreathing aircraft, multigrid methods for the steady incompressible Navier-Stokes equations, upwind differencing techniques, TV stable schemes for shock-interacting flows, Euler models of hypersonic vortex flows, parallel multilevel adaptive methods, and vortex methods for slightly viscous three-dimensional flows. Consideration is given to the accuracy of node-based solutions on irregular meshes, multigrid calculations for cascades, a finite-volume-element method for planar cavity flow, parallel heterogeneous mesh refinement for advection-diffusion equations, the convergence of the spectral-viscosity method for nonlinear conservation laws, and numerical simulations of Taylor vortices in a spherical gap.

Dwoyer, D. L.↗
NASA NTRS1990

Numerical simulations of the rebound shock model for solar spicules

Using time-dependent numerical simulations, the proposed rebound shock mechanism for spicules has been examined. At temperatures above a critical value, T(c), the radiation is characteristic of the conditions in the optically thin corona and near optically thin transition region. When T less than T(c), the atmosphere has a radiative cooling time, tau(rad) characteristic of chromosphere. The spicule is initiated with a quasi-impulsive force in the low chromosphere, which drives a train of upward propagating rebound shocks along the rigid magnetic flux tube. These shocks then move the transition region upward. The material below the displaced transition region has temperatures and densities similar to those of spicules when T(c) = 20,000 K or more and tau(rad) = 500 s or more, but not when T(c) = 10,000 K, and probably not when tau(rad) = 100 s. For all the cases where the cross sectional area diverges rapidly with height, the upward velocity of the transition region is less than that of spicules. Moreover, the maximum height is less than that of average spicules. Taller, higher velocity spicules result when the magnetic field cross sectional area is constant. In all cases, the rebound shock mechanism produces substantial motions and temperature and density variations in chromospheric and transition region material. It is suggested that this may be a partial explanation for the continuous dynamic state of the lower solar atmosphere.

Sterling, Alphonse C.↗
NASA NTRS1991

Vectorization of a Monte Carlo simulation scheme for nonequilibrium gas dynamics

Significant improvement has been obtained in the numerical performance of a Monte Carlo scheme for the analysis of nonequilibrium gas dynamics through an implementation of the algorithm which takes advantage of vector hardware, as presently demonstrated through application to three different problems. These are (1) a 1D standing-shock wave; (2) the flow of an expanding gas through an axisymmetric nozzle; and (3) the hypersonic flow of Ar gas over a 3D wedge. Problem (3) is illustrative of the greatly increased number of molecules which the simulation may involve, thanks to improved algorithm performance.

Boyd, Iain D.↗
NASA NTRS1992

Recent results on Mars and Venus; Proceedings of Symposium 3 and the Topical Meetings of the Interdisciplinary Scientific Commission C (Meeting C1) and D (Meeting D1) of the COSPAR 28th Plenary Meeting, The Hague, Netherlands, June 25-July 6, 1990

In this volume, latest results on Mars and Phobos studies are discussed along with recent results on Venus and on Martian plasma environment. Papers are presented on the characteristics of aerosol phenomena in Martian atmosphere from KRFM-experiment data, Phobos map and Phobos globe, submicron-sized dust grains in the Martian environment, remote sensing of Venus atmospheric dynamics, the latitude structure of upper clouds of Venus, and Venus thermospheric response to short-term solar variations. Attention is also given to a comparison between the bow shocks and magnetotails of Venus and Mars, comparative studies of the solar wind interaction with weakly magnetized planets, the magnetic field environment of Mars as studied by Phobos-2, waves and cold plasmas near Mars, and multiple-ion effects at Martian plasma boundaries.

Shorthill, R. W.↗
NASA NTRS1990

Heat Pipe Vapor Dynamics

The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.

Issacci, Farrokh↗
NASA NTRS1995

A computational fluid dynamics simulation of the hypersonic flight of the Pegasus(TM) vehicle using an artificial viscosity model and a nonlinear filtering method

The computational fluid dynamics code, PARC3D, is tested to see if its use of non-physical artificial dissipation affects the accuracy of its results. This is accomplished by simulating a shock-laminar boundary layer interaction and several hypersonic flight conditions of the Pegasus(TM) launch vehicle using full artificial dissipation, low artificial dissipation, and the Engquist filter. Before the filter is applied to the PARC3D code, it is validated in one-dimensional and two-dimensional form in a MacCormack scheme against the Riemann and convergent duct problem. For this explicit scheme, the filter shows great improvements in accuracy and computational time as opposed to the nonfiltered solutions. However, for the implicit PARC3D code it is found that the best estimate of the Pegasus experimental heat fluxes and surface pressures is the simulation utilizing low artificial dissipation and no filter. The filter does improve accuracy over the artificially dissipative case but at a computational expense greater than that achieved by the low artificial dissipation case which has no computational time penalty and shows better results. For the shock-boundary layer simulation, the filter does well in terms of accuracy for a strong impingement shock but not as well for weaker shock strengths. Furthermore, for the latter problem the filter reduces the required computational time to convergence by 18.7 percent.

Mendoza, John Cadiz↗
NASA NTRS1998

Comprehensive Structural Dynamic Analysis of the SSME/AT Fuel Pump First-Stage Turbine Blade

A detailed structural dynamic analysis of the Pratt & Whitney high-pressure fuel pump first-stage turbine blades has been performed to identify the cause of the tip cracking found in the turbomachinery in November 1997. The analysis was also used to help evaluate potential fixes for the problem. Many of the methods available in structural dynamics were applied, including modal displacement and stress analysis, frequency and transient response to tip loading from the first-stage Blade Outer Gas Seals (BOGS), fourier analysis, and shock spectra analysis of the transient response. The primary findings were that the BOGS tip loading is impulsive in nature, thereby exciting many modes of the blade that exhibit high stress at the tip cracking location. Therefore, a proposed BOGS count change would not help the situation because a clearly identifiable resonance situation does not exist. The recommendations for the resolution of the problem are to maintain the existing BOGS count, eliminate the stress concentration in the blade due to its geometric design, and reduce the applied load on the blade by adding shiplaps in the BOGS.

Brown, A. M.↗
NASA NTRS2001

Acoustical and Flowfield Characterization of a Scaled Tabletop Rocket

An analysis of the acoustical and flowfield environment for the scaled 1-pound-force (lbf) thrust tabletop motor was performed. This tabletop motor from NASA Stennis Space Center Is composed of Plexiglas burning In gaseous oxygen with a graphite insert for the nozzle portion. The nozzle has a throat diameter of 0.2 inch and an exit diameter of 0.38 Inch. With a chamber pressure at 55 pounds per square Inch absolute (psia), a normal shock is formed immediately downstream of the nozzle exit plane as the combustion products exhaust into the ambient at atmospheric pressure. The jet characterization Is based on computational fluid dynamics (CFD) in conjunction with Kirchhoff surface integral formulation and compared with correlations developed for measured rocket noise and a pressure fluctuation scaling (PFS) method. Predictions and comparisons are made for the overall sound pressure levels (OASPL's) and spectral dependence of sound pressure level (SPL). The overall objective of this effort is to develop methods for scaling the acoustic and flowfield environment of rockets with a wide range of thrust (1 lbf to 1 million lbf).

Kandula, Max↗
NASA NTRSDate not supplied

Experimental and Computational Analysis of a Rotating Detonation Combustor

Rotating detonation combustor (RDC) research has progressed along parallel experimental and computational paths with limited opportunities for validation. Validation is a crucial step to ensure RDC simulations provide an accurate representation of the physical phenomena, and thus, can be used with confidence for design optimization in different applications. In this study, data from a quasi-two-dimensional computational fluid dynamic simulation in an annular RDC are compared with experimental measurements including particle image velocimetry at 100 kHz. Results show good agreement between measured and simulated dynamic pressures, Chapman-Jouguet speeds, and static pressure distributions. Velocity measurements at the annulus exit show good temporal and quantitative agreement with simulations including spikes from the passing shock.

Detonation↗
NASA NTRS1978

Wind tunnel evaluation of YF-12 inlet response to internal airflow disturbances with and without control

The response of terminal-shock position and static pressures in the subsonic duct of a YF-12 aircraft flight-hardware inlet to perturbations in simulated engine corrected airflow were obtained with and without inlet control. Frequency response data, obtained with inlet controls inactive, indicated the general nature of the inherent inlet dynamics, assisted in the design of controls, and provided a baseline reference for responses with active controls. All the control laws were implemented by means of a digital computer that could be programmed to behave like the flight inlet's existing analog control. The experimental controls were designed using an analytical optimization technique. The capabilities of the controls were limited primarily by the actuation hardware. The experimental controls provided somewhat better attenuation of terminal shock excursions than did the YF-13 inlet control. Controls using both the forward and aft bypass systems also provided somewhat better attenuation than those using just the forward bypass. The main advantage of using both bypasses is in the greater control flexibility that is achieved.

Cole, G. L.↗