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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Solution of steady-state, two-dimensional conservation laws by mathematical programming

A truly two-dimensional algorithm is created for solving the steady-state two-dimensional conservation-law problem. An overdetermined system of algebraic equations is obtained through discretization by finite-volume formulas. These equations are perturbed nonsingularly and are solved by an efficient geometrically oriented l(1) procedure. The basic algorithm and the theory for the linear case f(u) = u are presented, and computational results for the nonlinear case f(u) = sq u are also analyzed. It is noted that the l(1) procedure captures boundary shocks as well as oblige and zigzag interior shocks in bands that are one cell wide, and the solution values are accurate up to the edge of the shock.

Lavery, John E.↗
NASA NTRS1990

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗
NASA NTRS1992

A test suite for magnetohydrodynamical simulations

A collection is presented of MHD problems that will provide researchers who have interests in modeling MHD phenomena with a battery of tests (a 'test suite') for calibrating their numerical algorithms. Use of these tests will provide a common reference for comparison of different numerical MHD algorithms. The test suite includes both 1D and 2D problems. Taken together, these problems test the abilities of a numerical method to propagate accurately all of the MHD wave families in moving and stationary media, to capture MHD shocks and contact discontinuities, and to model accurately Lorentz force terms in multiple dimensions. An example solution of each test problem is presented. Diagnostic convergence-testing procedures that provide a quantitative evaluation of MHD algorithms are demonstrated.

Stone, James M.↗
NASA NTRS1993

A method of predicting quasi-steady aerodynamics for flutter analysis of high speed vehicles using steady CFD calculations

High speed linear aerodynamic theories like piston theory and Newtonian impact theory are relatively inexpensive to use for flutter analysis. These theories have limited areas of applicability depending on the configuration and the flow conditions. In addition, these theories lack the ability to capture viscous, shock and real gas effects. CFD methods can model all of these effects accurately, but the unsteady calculations required for flutter are expensive and often impractical. This paper describes a method for using steady CFD calculations to approximate the generalized aerodynamic forces for a flutter analysis. Example two-and three-dimensional aerodynamic force calculations are provided. In addition, a flutter analysis of a NASP-type wing will be discussed.

Scott, Robert C.↗
NASA NTRS1993

A method of predicting quasi-steady aerodynamics for flutter analysis of high speed vehicles using steady CFD calculations

High speed linear aerodynamic theories like piston theory and Newtonian impact theory are relatively inexpensive to use for flutter analysis. These theories have limited areas of applicability depending on the configuration and the flow conditions. In addition, these theories lack the ability to capture viscous, shock, and real gas effects. CFD methods can model all of these effects accurately, but the unsteady calculations required for flutter are expensive and often impractical. This paper describes a method for using steady CFD calculations to approximate the generalized aerodynamic forces for a flutter analysis. Example two-and three-dimensional aerodynamic force calculations are provided. In addition, a flutter analysis of a NASP-type wing will be discussed.

Scott, Robert C.↗
NASA NTRS2000

Transonic Flow Computations Using Nonlinear Potential Methods

This presentation describes the state of transonic flow simulation using nonlinear potential methods for external aerodynamic applications. The presentation begins with a review of the various potential equation forms (with emphasis on the full potential equation) and includes a discussion of pertinent mathematical characteristics and all derivation assumptions. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, is discussed. Key characteristics of all numerical algorithm types used for solving nonlinear potential equations, including steady, unsteady, space marching, and design methods, are described. Both spatial discretization and iteration scheme characteristics are examined. Numerical results for various aerodynamic applications are included throughout the presentation to highlight key discussion points. The presentation ends with concluding remarks and recommendations for future work. Overall. nonlinear potential solvers are efficient, highly developed and routinely used in the aerodynamic design environment for cruise conditions. Published by Elsevier Science Ltd. All rights reserved.

Holst, Terry L.↗
NASA NTRS1997

Computation of Transonic Flows Using Potential Methods

The proposed paper will describe the state of the art associated with numerical solution of the full or exact velocity potential equation for solving transonic, external-aerodynamic flows. The presentation will begin with a review of the literature emphasizing research activities of the past decade. Next, the various forms of the full or exact velocity potential equation, the equation's corresponding mathematical characteristics, and the derivation assumptions will be presented and described in detail. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, will be presented and discussed relative to the more complete Euler or Navier-Stokes formulations. The technical presentation will continue with a description of recently developed full potential numerical approach characteristics. This description will include governing equation nondimensionalization, physical-to-computational-domain mapping procedures, a limited description of grid generation requirements, the spatial discretization scheme, numerical implementation of boundary conditions, and the iteration scheme. The next portion of the presentation will present and discuss numerical results for several two- and three-dimensional aerodynamic applications. Included in the results section will be a discussion and demonstration of a typical grid refinement analysis for determining spatial convergence of the numerical solution and level of solution accuracy. Computer timings for a variety of full potential applications will be compared and contrasted with similar results for the Euler equation formulation. Finally. the presentation will end with concluding remarks and recommendations for future work.

Hoist, Terry L.↗
NASA NTRS2013

Overview of the Aeroelastic Prediction Workshop

The Aeroelastic Prediction Workshop brought together an international community of computational fluid dynamicists as a step in defining the state of the art in computational aeroelasticity. This workshop's technical focus was prediction of unsteady pressure distributions resulting from forced motion, benchmarking the results first using unforced system data. The most challenging aspects of the physics were identified as capturing oscillatory shock behavior, dynamic shock-induced separated flow and tunnel wall boundary layer influences. The majority of the participants used unsteady Reynolds-averaged Navier Stokes codes. These codes were exercised at transonic Mach numbers for three configurations and comparisons were made with existing experimental data. Substantial variations were observed among the computational solutions as well as differences relative to the experimental data. Contributing issues to these differences include wall effects and wall modeling, non-standardized convergence criteria, inclusion of static aeroelastic deflection, methodology for oscillatory solutions, post-processing methods. Contributing issues pertaining principally to the experimental data sets include the position of the model relative to the tunnel wall, splitter plate size, wind tunnel expansion slot configuration, spacing and location of pressure instrumentation, and data processing methods.

Heeg, Jennifer↗
NASA NTRS1974

Numerical computation of two dimensional viscous blunt body flows with an impinging shock, part 2

Two-dimensional viscous blunt body flows with an impinging shock have been computed using a time-dependent finite-difference method which solves the complete set of Navier-Stokes equations for a compressible flow. For low Reynolds number flows, the entire flow field, including the bow shock and impinging shock, has been captured in the computation. For higher Reynolds number flows, the bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied, while the boundary layer and interaction regions are captured as before. Using this latter shock-fitting approach, a Type III shock interaction flow field has been computed with flow conditions corresponding to the space shuttle orbiter freestream conditions at 61 km (200,000 ft).

Holst, T. L.↗
NASA NTRS1994

Prediction of Undsteady Flows in Turbomachinery Using the Linearized Euler Equations on Deforming Grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. The model accounts for blade loading, blade geometry, shock motion, and wake motion. Assuming that the unsteadiness in the flow is small relative to the nonlinear mean solution, the unsteady Euler equations can be linearized about the mean flow. This yields a set of linear variable coefficient equations that describe the small amplitude harmonic motion of the fluid. These linear equations are then discretized on a computational grid and solved using standard numerical techniques. For transonic flows, however, one must use a linear discretization which is a conservative linearization of the non-linear discretized Euler equations to ensure that shock impulse loads are accurately captured. Other important features of this analysis include a continuously deforming grid which eliminates extrapolation errors and hence, increases accuracy, and a new numerically exact, nonreflecting far-field boundary condition treatment based on an eigenanalysis of the discretized equations. Computational results are presented which demonstrate the computational accuracy and efficiency of the method and demonstrate the effectiveness of the deforming grid, far-field nonreflecting boundary conditions, and shock capturing techniques. A comparison of the present unsteady flow predictions to other numerical, semi-analytical, and experimental methods shows excellent agreement. In addition, the linearized Euler method presented requires one or two orders-of-magnitude less computational time than traditional time marching techniques making the present method a viable design tool for aeroelastic analyses.

Clark, William S.↗
NASA NTRS1989

A computational study on the interaction between a vortex and a shock wave

A computational study of two-dimensional shock vortex interaction is discussed in this paper. A second order upwind finite volume method is used to solve the Euler equations in conservation form. In this method, the shock wave is captured rather than fitted so that the cases where shock vortex interaction may cause secondary shocks can also be investigated. The effects of vortex strength on the computed flow and acoustic field generated by the interaction are qualitatively evaluated.

Meadows, Kristine R.↗
NASA NTRS1994

Comparison of coordinate-invariant and coordinate-aligned upwinding for the Euler equations

A floating-shock fitting method for the Euler equations has been developed that uses one-sided spatial differences along and across streamlines. The coordinate-invariant formulation of the spatial differences permits automatic capture of shears. Results are presented for unsteady shocked flow in a duct with a ramp, for supercritical flow over a circular cylinder, and for subsonic, transonic, and supersonic (0.3 is less than or equal to M(sub infinity) is less than 1.5) flow over airfoils. In flows with strong shears, the coordinate-invariant differencing concept appears to yield some gains in accuracy over Euler methods that rely on coordinate-aligned differencing concepts. In applications to transonic airfoils, fitted shocks have a tendency to be predicted upstream of captured shocks, regardless of whether coordinate-invariant or coordinate-aligned differencing is used. The coordinate-invariant differencing method requires between 2 and 3.5 times as much computing time as its coordinate-aligned counterpart.

Hartwich, Peter M.↗
NASA NTRS1986

Acoustic propagation using computational fluid dynamics

The propagation characteristics of several helicopter airfoil profiles have been investigated using the transonic small disturbance equation. A test case was performed to generate a moving shock that propagated off the airfoil. Various grids were then examined to determine their ability to accurately capture these propagating shock waves. Finally, the case of airfoil-vortex interactions was thoroughly studied over a wide range of Mach numbers and airfoil shapes with particular emphasis on the transonic regime; this results in a highly conplicated fluctuation of lift, drag, and pitching moment. The calculated acoustic intensity levels, along with the details of the computational flow field, provide new insights into the understanding of transonic airfoil-vortex interactions.

Baeder, J. D.↗
NASA NTRS1993

Shock-vortex interaction over a 65-degree delta wing in transonic flow

Transonic flow over a 65-deg swept-back, sharp-edged, cropped delta wing is investigated computationally using the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with an implicit, upwind, flux-difference splitting, finite-volume scheme. Coarse and fine O-H grids are used to obtain the solution. The grid consists of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively. The results are presented for an angle of attack of 20 deg Mach number of 0.85 and Reynolds number of 3.23 x 10 exp 6. With the fine grid, the results show that a system of shocks has been captured over the upper wing surface and that the leading-edge vortex core experiences an unsteady supersonic vortex breakdown after passing through a spanwise shock near the wing trailing edge. The computed results at a certain time are in good agreement with the experimental data. Topological aspects of the vortex breakdown flowfield are also presented and discussed.

Kandil, Osama A.↗
NASA NTRS1984

Molecular processes in a high temperature shock layer

Models of the shock layer encountered by an Aeroassisted Orbital Transfer Vehicle require as input accurate cross sections and rate constants for the atomic and molecular processes that characterize the shock radiation. From the estimated atomic and molecular densities in the shock layer and the expected residence time of 1 m/s, it can be expected that electron-ion collision processes will be important in the shock model. Electron capture by molecular ions followed by dissociation, e.g., O2(+) + e(-) yields 0 + 0, can be expected to be of major importance since these processes are known to have high rates (e.g., 10 to the -7th power cu/cm/sec) at room temperature. However, there have been no experimental measurements of dissociative recombination (DR) at temperatures ( 12000K) that are expected to characterize the shock layer. Indeed, even at room temperature, it is often difficult to perform experiments that determine the dependence of the translational energy and quantum yields of the product atoms on the electronic and vibrational state of the reactant molecular ions. Presented are ab initio quantum chemical studies of DR for molecular ions that are likely to be important in the atmospheric shock layer.

Guberman, S. L.↗
NASA NTRS1984

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic, and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching step sizes with the implicit shock fitting procedure. A correlation is found between the marching step size, Reynolds number, and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗
NASA NTRS1985

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching-step sizes with the implicit shock fitting procedure. A correlation is found between the marching-step size, Reynolds number and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗
NASA NTRS1992

Development of a 3-D upwind PNS code for chemically reacting hypersonic flowfields

Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.

Tannehill, J. C.↗