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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Transonic Navier-Stokes computations of strake-generated vortex interactions for a fighter-like configuration

Transonic Euler/Navier-Stokes computations are accomplished for wing-body flow fields using a computer program called Transonic Navier-Stokes (TNS). The wing-body grids are generated using a program called ZONER, which subdivides a coarse grid about a fighter-like aircraft configuration into smaller zones, which are tailored to local grid requirements. These zones can be either finely clustered for capture of viscous effects, or coarsely clustered for inviscid portions of the flow field. Different equation sets may be solved in the different zone types. This modular approach also affords the opportunity to modify a local region of the grid without recomputing the global grid. This capability speeds up the design optimization process when quick modifications to the geometry definition are desired. The solution algorithm embodied in TNS is implicit, and is capable of capturing pressure gradients associated with shocks. The algebraic turbulence model employed has proven adequate for viscous interactions with moderate separation. Results confirm that the TNS program can successfully be used to simulate transonic viscous flows about complicated 3-D geometries.

Reznick, Steve↗

NASA NTRS1988

Adaptive computational methods for aerothermal heating analysis

The development of adaptive gridding techniques for finite-element analysis of fluid dynamics equations is described. The developmental work was done with the Euler equations with concentration on shock and inviscid flow field capturing. Ultimately this methodology is to be applied to a viscous analysis for the purpose of predicting accurate aerothermal loads on complex shapes subjected to high speed flow environments. The development of local error estimate strategies as a basis for refinement strategies is discussed, as well as the refinement strategies themselves. The application of the strategies to triangular elements and a finite-element flux-corrected-transport numerical scheme are presented. The implementation of these strategies in the GIM/PAGE code for 2-D and 3-D applications is documented and demonstrated.

Price, John M.↗

NASA NTRS1989

Upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes equations. This method does not require the addition of user-specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming (1978) scheme in terms of accuracy, stability, computer time and storage requirements, and programming effort. The new algorithm has been validated by applying it to three laminar test cases, including flat-plate boundary-layer flow, hypersonic flow past a 15-deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with results obtained using the conventional Beam-Warming algorithm.

Lawrence, Scott L.↗

NASA NTRS1991

Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation

Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in a high gradient region or the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational costs. A detailed description is given of the enrichment and coarsening procedures and comparisons with alternative results and experimental data are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.

Rausch, Russ D.↗

NASA NTRS1991

Computation of compressible quasi-axisymmetric slender vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux difference splitting finite volume scheme. The developed three dimensional solver was verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic, quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown were captured. The problem was also calculated using the Euler solver of the same code; the results were compared with those of the Navier-Stokes solver. The effect of the initial swirl was investigated.

Kandil, Osama A.↗

NASA NTRS1991

Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation

Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with alternative results and experimental data to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.

Rausch, Russ D.↗

NASA NTRS1992

Parabolized reduced Navier-Stokes computational techniques

A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.

Rubin, Stanley G.↗

NASA NTRS1994

A Continuing Search for a Near-Perfect Numerical Flux Scheme

While enjoying demonstrated improvement in accuracy, efficiency, and robustness over existing schemes, the Advection Upstream Splitting Scheme (AUSM) was found to have some deficiencies in extreme cases. This recent progress towards improving the AUSM while retaining its advantageous features is described. The new scheme, termed AUSM+, features: unification of velocity and Mach number splitting; exact capture of a single stationary shock; and improvement in accuracy. A general construction of the AUSM+ scheme is layed out and then focus is on the analysis of the a scheme and its mathematical properties, heretofore unreported. Monotonicity and positivity are proved, and a CFL-like condition is given for first and second order schemes and for generalized curvilinear co-ordinates. Finally, results of numerical tests on many problems are given to confirm the capability and improvements on a variety of problems including those failed by prominent schemes.

Liou, Meng-Sing↗

NASA NTRS2002

Aeroacoustic Flow Phenomena Accurately Captured by New Computational Fluid Dynamics Method

One of the challenges in the computational fluid dynamics area is the accurate calculation of aeroacoustic phenomena, especially in the presence of shock waves. One such phenomenon is "transonic resonance," where an unsteady shock wave at the throat of a convergent-divergent nozzle results in the emission of acoustic tones. The space-time Conservation-Element and Solution-Element (CE/SE) method developed at the NASA Glenn Research Center can faithfully capture the shock waves, their unsteady motion, and the generated acoustic tones. The CE/SE method is a revolutionary new approach to the numerical modeling of physical phenomena where features with steep gradients (e.g., shock waves, phase transition, etc.) must coexist with those having weaker variations. The CE/SE method does not require the complex interpolation procedures (that allow for the possibility of a shock between grid cells) used by many other methods to transfer information between grid cells. These interpolation procedures can add too much numerical dissipation to the solution process. Thus, while shocks are resolved, weaker waves, such as acoustic waves, are washed out.

Blech, Richard A.↗

NASA NTRS1990

Comparison of several finite-difference methods

The present numerical study establishes that, with the MacCormack and 'two-four' methods, the solution accuracy obtainable in problems involving wave propagation, shock-wave and contact discontinuities, and viscous effects, will be strongly dependent on Courant number. The application of flux correction to the MacCormack and two-four methods is noted to significantly attenuate dispersion errors; the ensuing solutions capture the discontinuities in the shock-tube problem with improved accuracy and resolution, and are free of dispersion errors for the viscous Burgers' equation.

Biringen, S.↗

NASA NTRS1975

Numerical computation of two-dimensional viscous blunt body flows with an impinging shock

Two-dimensional, viscous, blunt body flows with an impinging shock wave are computed using a time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock wave is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which shock waves of different strengths are allowed to impinge on the flow field surrounding a circular cylinder resulting in different shock interference patterns. The two-dimensional results are compared qualitatively with existing three-dimensional experiments.

Tannehil, J. C.↗

NASA NTRS2017

Quantitative Image Analysis Techniques with High-Speed Schlieren Photography

Optical flow visualization techniques such as schlieren and shadowgraph photography are essential to understanding fluid flow when interpreting acquired wind tunnel test data. Output of the standard implementations of these visualization techniques in test facilities are often limited only to qualitative interpretation of the resulting images. Although various quantitative optical techniques have been developed, these techniques often require special equipment or are focused on obtaining very precise and accurate data about the visualized flow. These systems are not practical in small, production wind tunnel test facilities. However, high-speed photography capability has become a common upgrade to many test facilities in order to better capture images of unsteady flow phenomena such as oscillating shocks and flow separation. This paper describes novel techniques utilized by the authors to analyze captured high-speed schlieren and shadowgraph imagery from wind tunnel testing for quantification of observed unsteady flow frequency content. Such techniques have applications in parametric geometry studies and in small facilities where more specialized equipment may not be available.

Pollard, Victoria J.↗

NASA NTRS2013

Hybrid Large-Eddy/Reynolds-Averaged Simulation of a Supersonic Cavity Using VULCAN

Simulations of a supersonic recessed-cavity flow are performed using a hybrid large-eddy/Reynolds-averaged simulation approach utilizing an inflow turbulence recycling procedure and hybridized inviscid flux scheme. Calorically perfect air enters a three-dimensional domain at a free stream Mach number of 2.92. Simulations are performed to assess grid sensitivity of the solution, efficacy of the turbulence recycling, and the effect of the shock sensor used with the hybridized inviscid flux scheme. Analysis of the turbulent boundary layer upstream of the rearward-facing step for each case indicates excellent agreement with theoretical predictions. Mean velocity and pressure results are compared to Reynolds-averaged simulations and experimental data for each case and indicate good agreement on the finest grid. Simulations are repeated on a coarsened grid, and results indicate strong grid density sensitivity. Simulations are performed with and without inflow turbulence recycling on the coarse grid to isolate the effect of the recycling procedure, which is demonstrably critical to capturing the relevant shear layer dynamics. Shock sensor formulations of Ducros and Larsson are found to predict mean flow statistics equally well.

Quinlan, Jesse↗

NASA NTRS1974

Numerical computation of viscous blunt body flows with a planar impinging shock

Two- and three-dimensional, viscous blunt body flows with planar impinging shocks are computed using an explicit, time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which planar shock waves of different strengths and orientations are allowed to impinge on the flow field surroundings an infinite cylinder resulting in two- and three-dimensional shock interference patterns. The numerical results are compared with experiment.

Holst, T. L.↗

NASA NTRS2022

Energy Partition at Collisionless Supercritical Quasiperpendicular Shocks

Collisionless shocks in astrophysical plasmas are important thermalizers, converting some of the incident flow energy into thermal energy, and non-thermalizers, partitioning that energy in unequal ways to different particle species, sub-populations thereof, and field components. This partition problem, or equivalently the shock equation of state, lies at the heart of shock physics. Here we employ systematically a framework to capture all the incident and downstream energy fluxes at two example traversals of the Earth’s bow shock by the Magnetospheric Multiscale Mission. Here and traditionally such data has to be augmented by information from other spacecraft, e.g., to provide more accurate measurements of the cold solar wind beam. With some care and fortuitous choices, the energy fluxes are constant, including instantaneous measurements through the shock layer. The dominant incident proton ram energy is converted primarily into downstream proton enthalpy flux, the majority of which is actually carried by a small fraction of suprathermal protons. Fluctuations include both real and instrumental effects. Separating these, resolving the solar wind beam, and other considerations point the way to a dedicated mission to solve this energy partition problem across a full range of plasma and shock conditions.

Steven J Schwartz↗

NASA NTRS2003

Supersonic Test of the 10-Inch Bifurcated Two-Stage Supersonic Inlet (BTSSI)

Experimental results from testing of a novel supersonic inlet model in NASA Glenn Research Center's 10- by 10-foot supersonic wind tunnel are presented. The patented inlet concept, called Two-Stage Supersonic Inlet (TSSI), incorporates a large cavity, or throat slot, in the supersonic diffuser intended to enhance the stability of the normal shock. The present embodiment of the concept is a bifurcated twin-duct) design. During the course of testing an unusual 'semi-started' mode of operation was encountered. The inlet was able to spill up to 30 percent of the captured airstream without fully expelling the normal shock. In this mode, the total pressure recovery dropped approximately 6 percent without increasing steady-state distortion. Dynamic instrumentation at the cowl lip station indicates the semi-start mode may be a series of unstart/restart cycles with frequency ranging from 0.2 to 20 Hz. Engine face total pressure measurements indicate a modest impact due to this event. However, since the current test article does not have a representative subsonic diffuser (and is in fact separated), it is unclear how this mode of operation would effect an engine. Further investigation of this phenomenon is required before it is fully understood. Prior testing of the TSSI concept allowed extension of fully started inlet operation to regions of significantly reduced supply flow without reducing recovery. The test article was a smaller scale than the present test and was a single duct design. In the present test, the expanded range of stable operation with high recovery was not realized.

Carlin, C. M.↗

NASA NTRS1991

Computation of steady and unsteady compressible quasi-axisymmetric vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux-difference splitting finite-volume scheme. The developed three-dimensional solver has been verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, isolated quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown have been captured. The problem has also been calculated using the Euler solver of the same code and the results are compared with those of the Navier-Stokes solver. The effect of the initial swirl has been tentatively studied.

Kandil, Osama A.↗

NASA NTRS2008

Forebody and Inlet Design for the HIFiRE 2 Flight Test

A forebody and inlet have been designed for the HIFiRE 2 scramjet flight test. The test will explore the operating, performance, and stability characteristics of a simple hydrocarbon-fueled scramjet combustor as it transitions from dual-mode to scramjet-mode operation and during supersonic combustion at Mach 8+ flight conditions. Requirements for the compression system were derived from inlet starting and combustor inflow requirements as well as physical size constraints. The design process is described. A planar, fixed geometry, mixed compression concept was used to produce laterally uniform flow at the inlet entrance and a conservative amount of internal contraction with respect to inlet starting. A grid sensitivity study was performed so that important flow physics caused by three-dimensional shock boundary layer interactions could be captured with confidence. Results from low Mach number operability studies, nominal trajectory cases, and high dynamic pressure heat load cases are discussed. The forebody and inlet solutions provide information for on-going combustor calculations, mass capture across the trajectory for fuel system design, and surface heating rates for thermal/structural analysis. The design has a one freestream Mach number margin for inlet starting, exceeds the high Mach number combustor entrance pressure requirement, produces high quality flow at the inlet exit for all Mach numbers and vehicle attitudes in the design space, and fits inside the booster shroud.

Ferlemann, Paul G.↗