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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Interpretation of waverider performance data using computational fluid dynamics

A computational study was conducted to better understand experimental results obtained from wind-tunnel tests of a Mach 4 conical-flow-derived waverider and a comparative reference configuration, which showed that the aerodynamic performance of the reference configuration was slightly better than that of the waverider. The computational results showed that the predicted surface pressure values and the integrated lift and drag coefficients were much lower for the reference model because the reference model bottom is an expansion surface. However, the lift-drag ratios for the reference model were higher due to a relatively low drag for a comparable amount of lift. The results also showed that the reference model exhibited the same shock attachment characteristics as the conical-flow-derived waverider, and is therefore also a waverider. The shock attachment characteristic gives the waverider a performance advantage over conventional hypersonic vehicles, and the results suggest that altering the bottom surface does not cause significant performance degradation. Flowfield solutions also show that the conical-flow waverider model has better propulsion/airframe integration characteristics than the reference configuration. The results also suggest that generating flowfields other than conical ones may be used to design waveriders with improved aerodynamic performance.

COMPUTATIONAL FLUID DYNAMICS↗
NASA NTRS1981

Supersonic Stall Flutter of High Speed Fans

An analytical model is developed for predicting the onset of supersonic stall bending flutter in axial flow compressors. The analysis is based on a modified two dimensional, compressible, unsteady actuator disk theory. It is applied to a rotor blade row by considering a cascade of airfoils whose geometry and dynamic response coincide with those of a rotor blade element at 85 percent of the span height (measured from the hub). The rotor blades are assumed to be unshrouded (i.e., free standing) and to vibrate in their first flexural mode. The effects of shock waves and flow separation are included in the model through quasi-steady, empirical, rotor total-pressure-loss and deviation-angle correlations. The actuator disk model predicts the unsteady aerodynamic force acting on the cascade blading as a function of the steady flow field entering the cascade and the geometry and dynamic response of the cascade. Calculations show that the present model predicts the existence of a bending flutter mode at supersonic inlet Mach numbers. This flutter mode is suppressed by increasing the reduced frequency of the system or by reducing the steady state aerodynamic loading on the cascade. The validity of the model for predicting flutter is demonstrated by correlating the measured flutter boundary of a high speed fan stage with its predicted boundary. This correlation uses a level of damping for the blade row (i.e., the log decrement of the rotor system) that is estimated from the experimental flutter data. The predicted flutter boundary is shown to be in good agreement with the measured boundary.

Adamczyk, J. J.↗
NASA NTRS1988

Origin of petrofabrics and magnetic anisotropy in ordinary chondrites

Three-dimensional finite strain and magnetic susceptibility anisotropy have been determined for 15 ordinary chondrites. The axes of strain and magnetic ellipsoids roughly correlate in both magnitude and orientation. The shapes of these ellipsoids are generally oblate spheroids that define a dominant foliation and a weak lineation. These characteristics suggest deformation involving uniaxial compaction. The degree of uniaxial deformation correlates with intensity of shock, as indicated by optical, TEM and chemical criteria. These data, plus the lack of a relationship between foliation and metamorphic history, indicate that dynamic processes, i.e., impacts, produced planar deformation fabrics in chondrites.

Sneyd, Deana S.↗
NASA NTRS1988

Numerical methods for turbulent flow

It has generally become accepted that the Navier-Strokes equations predict the dynamic behavior of turbulent as well as laminar flows of a fluid at a point in space away form a discontinuity such as a shock wave. Turbulence is also closely related to the phenomena of non-uniqueness of solutions of the Navier-Strokes equations. These second order, nonlinear partial differential equations can be solved analytically for only a few simple flows. Turbulent flow fields are much to complex to lend themselves to these few analytical methods. Numerical methods, therefore, offer the only possibility of achieving a solution of turbulent flow equations. In spite of recent advances in computer technology, the direct solution, by discrete methods, of the Navier-Strokes equations for turbulent flow fields is today, and in the foreseeable future, impossible. Thus the only economically feasible way to solve practical turbulent flow problems numerically is to use statistically averaged equations governing mean-flow quantities. The objective is to study some recent developments relating to the use of numerical methods to study turbulent flow.

Turner, James C., Jr.↗
NASA NTRS1989

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗
NASA NTRS1989

Calculation of nonequilibrium hydrogen-air reactions with implicit flux vector splitting method

Two methods, fully- and loosely-coupled, are developed to incorporate nonequilibrium hydrogen-air chemistry into the fluid dynamic implicit flux vector splitting code (F3D). The new code (F3D/Chem) is validated against other existing codes for two cases: nozzle expansion, and shock-induced combustion around a blunt body. The shock-induced combustion case is compared also with an experimental data. The reaction rate constants are varied in an effort to reproduce the experimental data. The fully- and loosely-coupled methods are found to yield comparable results, but the computation time is shorter using the loosely-coupled method. The present method is found to reproduce results obtained using different existing codes. The experimental data was not reproduced with any selected rate coefficients set.

Lee, Seung-Ho↗
NASA NTRS1992

Transonic aeroelastic computations on wings using Navier-Stokes equations

Aeroelastic computations are made on fighter type wings in the transonic regime. The flow is modeled using the Navier-Stokes equations and is coupled with structural equations of motion. The flow equations are solved by a time accurate finite difference scheme with moving grids. The coupled aeroelastic equations of motion are solved using the linear acceleration method. The configuration adaptive dynamic grids are time-accurately generated using the aeroelastically deformed shape of the wing. Computations are made for oscillating rigid wings with moving shock waves in the presence of leading edge vortices. The computed results compare well with the experiment. Unsteady computations are made to demonstrate the shock-vortex interaction phenomenon on wings in ramp motion. Effects of flexibility and pitch rate are demonstrated for flows with vortices.

Guruswamy, Guru P.↗
NASA NTRS1992

User's Manual for Aerofcn: a FORTRAN Program to Compute Aerodynamic Parameters

The computer program AeroFcn is discussed. AeroFcn is a utility program that computes the following aerodynamic parameters: geopotential altitude, Mach number, true velocity, dynamic pressure, calibrated airspeed, equivalent airspeed, impact pressure, total pressure, total temperature, Reynolds number, speed of sound, static density, static pressure, static temperature, coefficient of dynamic viscosity, kinematic viscosity, geometric altitude, and specific energy for a standard- or a modified standard-day atmosphere using compressible flow and normal shock relations. Any two parameters that define a unique flight condition are selected, and their values are entered interactively. The remaining parameters are computed, and the solutions are stored in an output file. Multiple cases can be run, and the multiple case solutions can be stored in another output file for plotting. Parameter units, the output format, and primary constants in the atmospheric and aerodynamic equations can also be changed.

Conley, Joseph L.↗
NASA NTRS1992

Wind-tunnel based definition of the AFE aerothermodynamic environment

The Aeroassist Flight Experiment (AFE), scheduled to be performed in 1994, will serve as a precursor for aeroassisted space transfer vehicles (ASTV's) and is representative of entry concepts being considered for missions to Mars. Rationale for the AFE is reviewed briefly as are the various experiments carried aboard the vehicle. The approach used to determine hypersonic aerodynamic and aerothermodynamic characteristics over a wide range of simulation parameters in ground-based facilities is presented. Facilities, instrumentation and test procedures employed in the establishment of the data base are discussed. Measurements illustrating the effects of hypersonic simulation parameters, particularly normal-shock density ratio (an important parameter for hypersonic blunt bodies), and attitude on aerodynamic and aerothermodynamic characteristics are presented, and predictions from computational fluid dynamic (CFD) computer codes are compared with measurement.

Miller, Charles G.↗
NASA NTRS1999

Simulation of Non-Linear Transonic Aeroelastic Behavior on the B-2

At high subsonic flight speeds, large flexible aircraft begin to encounter unsteady airloads which are not predicted by most currently available aerodynamic analysis and design methods. With increasing speed and the development of transonic flow and shocks, viscous effects quickly become very important, and flow separation can occur. The Northrop Grumman USAF B-2 Bomber encountered a nonlinear aeroelastic Residual Pitch Oscillation (RPO) under these conditions. Simulation studies were performed with the Computational Aeroelasticity Program-Transonic Small Disturbance, Viscous (CAPTSDv) computer program to evaluate its ability to predict these nonlinear aeroelastic responses. Open and closed loop simulations were performed to assess the participation of the flight control system. Control, system actuator hysteresis characteristics were modeled and found to be a significant participant in the RPO phenomenon. Simulations were also performed for varying Mach numbers and altitudes to establish the stability boundaries and compare with flight test data. These CAPTSDv simulations compared well with flight data and revealed many potential further modeling enhancements.

Residual pitch oscillation↗
NASA NTRS2011

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are candidate materials for applications in extreme environments because of their high melting point, good mechanical properties and reasonable oxidation resistance. Unlike many ceramics, these materials have high thermal conductivity which can be advantageous, for example, to reduce thermal shock. Recently, we developed Tersoff style interatomic potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current. Results at room temperature and at elevated temperatures will be reported.

Lawson, JOhn W.↗
NASA NTRS1994

Atmospheric Disturbance Environment Definition

Traditionally, the application of atmospheric disturbance data to airplane design problems has been the domain of the structures engineer. The primary concern in this case is the design of structural components sufficient to handle transient loads induced by the most severe atmospheric "gusts" that might be encountered. The concern has resulted in a considerable body of high altitude gust acceleration data obtained with VGH recorders (airplane velocity, V, vertical acceleration, G, altitude, H) on high-flying airplanes like the U-2 (Ehernberger and Love, 1975). However, the propulsion system designer is less concerned with the accelerations of the airplane than he is with the airflow entering the system's inlet. When the airplane encounters atmospheric turbulence it responds with transient fluctuations in pitch, yaw, and roll angles. These transients, together with fluctuations in the free-stream temperature and pressure will disrupt the total pressure, temperature, Mach number and angularity of the inlet flow. For the mixed compression inlet, the result is a disturbed throat Mach number and/or shock position, and in extreme cases an inlet unstart can occur (cf. Section 2.1). Interest in the effects of inlet unstart on the vehicle dynamics of large, supersonic airplanes is not new. Results published by NASA in 1962 of wind tunnel studies of the problem were used in support of the United States Supersonic Transport program (SST) (White, at aI, 1963). Such studies continued into the late 1970's. However, in spite of such interest, there never was developed an atmospheric disturbance database for inlet unstart analysis to compare with that available for the structures load analysis. Missing were data for the free-stream temperature and pressure disturbances that also contribute to the unStart problem.

Tank, William G.↗
NASA NTRS2013

Modeling Inflation Using a Non-Equilibrium Equation of Exchange

Inflation is a change in the prices of goods that takes place without changes in the actual values of those goods. The Equation of Exchange, formulated clearly in a seminal paper by Irving Fisher in 1911, establishes an equilibrium relationship between the price index P (also known as "inflation"), the economy's aggregate output Q (also known as "the real gross domestic product"), the amount of money available for spending M (also known as "the money supply"), and the rate at which money is reused V (also known as "the velocity of circulation of money"). This paper offers first a qualitative discussion of what can cause these factors to change and how those causes might be controlled, then develops a quantitative model of inflation based on a non-equilibrium version of the Equation of Exchange. Causal relationships are different from equations in that the effects of changes in the causal variables take time to play out-often significant amounts of time. In the model described here, wages track prices, but only after a distributed lag. Prices change whenever the money supply, aggregate output, or the velocity of circulation of money change, but only after a distributed lag. Similarly, the money supply depends on the supplies of domestic and foreign money, which depend on the monetary base and a variety of foreign transactions, respectively. The spreading of delays mitigates the shocks of sudden changes to important inputs, but the most important aspect of this model is that delays, which often have dramatic consequences in dynamic systems, are explicitly incorporated.macroeconomics, inflation, equation of exchange, non-equilibrium, Athena Project

Chamberlain, Robert G.↗
NASA NTRS2003

The Formation of Relativistic Jets from Kerr Black Holes

We have performed the first fully three-dimensional general relativistic magnetohydrodynamics (GRMHD) simulation for Schwarzschild and Kerr black holes with a free falling corona and thin accretion disk. The initial simulation results with a Schwarzschild metric show that a jet is created as in the previous axisymmetric simulations with mirror symmetry at the equator. However, the time to form the jet is slightly longer than in the 2-D axisymmetric simulation. We expect that the dynamics of jet formation are modified due to the additional freedom in the azimuth dimension without axisymmetry with respect to the Z axis and reflection symmetry respect to the equatorial plane. The jet which is initially formed due to the twisted magnetic fields and shocks becomes a wind at the later time. The wind flows out with a much wider angle than the initial jet. The twisted magnetic fields at the earlier time were untwisted and less pinched. The accretion disk became thicker than the initial condition. Further simulations with initial perturbations will provide insights for accretion dynamics with instabilities such as magneto-rotational instability (MRI) and accretion-eject instability (AEI). These instabilities may contribute to variabilities observed in microquasars and AGN jets.

Nishikawa, K.-I.↗
NASA NTRS1994

Significance of shock structure on supersonic jet mixing noise of axisymmetric nozzles

One of the key technical elements in NASA's high speed research program is reducing the noise level to meet the federal noise regulation. The dominant noise source is associated with the supersonic jet discharged from the engine exhaust system. Whereas the turbulence mixing is largely responsible for the generation of the jet noise, a broadband shock-associated noise is also generated when the nozzle operates at conditions other than its design. For both mixing and shock noise components, because the source of the noise is embedded in the jet plume, one can expect that jet noise can be predicted from the jet flowfield computation. Mani et al. developed a unified aerodynamic/acoustic prediction scheme by applying an extension of Reichardt's aerodynamic model to compute turbulent shear stresses which are utilized in estimating the strength of the noise source. Although this method produces a fast and practical estimate of the jet noise, a modification by Khavaran et al. has led to an improvement in aerodynamic solution. The most notable feature in this work is that Reichardt's model is replaced with the computational fluid dynamics (CFD) solution of Reynolds-averaged Navier-Stokes equations. The major advantage of this work is that the essential, noise-related flow quantities such as turbulence intensity and shock strength can be better predicted. The predictions were limited to a shock-free design condition and the effect of shock structure on the jet mixing noise was not addressed. The present work is aimed at investigating this issue. Under imperfectly expanded conditions the existence of the shock cell structure and its interaction with the convecting turbulence structure may not only generate a broadband shock-associated noise but also change the turbulence structure, and thus the strength of the mixing noise source. Failure in capturing shock structures properly could lead to incorrect aeroacoustic predictions.

Kim, Chan M.↗
NASA NTRS1995

Locations of termination shock and heliopause based on Voyager plasma and magnetic field data

The locations of the termination shock and the heliopause are studied taking into account the effects of pickup protons. The study uses available plasma and magnetic field data from Voyagers over a 14-year period (1978-1991) and Voyager observation of the 1992-93 radio emission event. Outside 30 AU, pickup protons have a significant influence on dynamical structures of the outer heliosphere. The solar wind is treated as a mixture of electrons, solar wind protons, and interstellar pickup protons. If the magnitude of the interstellar magnetic field B(sub int) is given, one can quantitatively study the motion and location of the termination shock. The location is anti-correlated with the sun spot number and the shock has an average speed of approx. 24 km/s. Because B(sub int) is poorly known, additional information is needed in studying the termination shock. Cummings, et al. have used observations of anomalous cosmic rays to estimate the location of the shock. The observations of the 1991 GMIR and GMIR shock and the 1992-93 radio emission event provide another handle for the study of the termination shock and the heliopause. After its penetration through the termination shock, the GMIR shock continued to propagate in the subsonic region of the solar wind and eventually interacted with the heliopause. This interaction produces a transmitted shock propagating outward in the interstellar medium and a reflected shock propagating inward toward the sun in the subsonic solar wind. The plasma frequencies behind the reflected and the transmitted shock can be, respectively, responsible for the 2- and 3-kHz radio emissions. Taking into account the effects of pickup protons we found that the average locations of the termination shock and the heliopause in 1991-92 are at approximately 66 AU and 150 AU, respectively.

Whang, Y. C.↗
NASA NTRS1993

Alfven wave transport effects in the time evolution of parallel cosmic-ray-modified shocks

This paper presents a numerical study of the time evolution of plane, cosmic-ray modified shocks with magnetic field parallel to the shock normal, based on the diffusive shock acceleration formalism and including the effects from the finite propagation speed and energy of Alfven waves responsible for controlling the transport of the cosmic rays. The simulations discussed are based on a three-fluid model for the dynamics, but a more complete formalism is laid out for future work. The results of the simulations confirm earlier steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and the gas sound speed are comparable, i.e., when the plasma and magnetic pressures are similar. It is also clear, however, that the impact of Alfven transport effects, which tend to slow shock evolution and reduce the time asymptotic cosmic-ray pressure in the shock, is strongly dependent upon uncertain details in the transport models. Both cosmic-ray advection tied to streaming Alfven waves and dissipation of wave energy are important to include in the models. Further, Alfven transport properties on both sides of the shock are also influential.

Jones, T. W.↗
NASA NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗