Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock fitting”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Calculation of the three-dimensional flow field in supersonic inlets at angle of attack using a bicharacteristic method with discrete shock wave fitting

An analysis is presented for calculating the flow field in supersonic mixed-compression aircraft inlets operating at angle of attack. The flow field is computed by a steady three-dimensional bicharacteristic method. The bow shock wave and the reflected internal shock wave system are computed by a three-dimensional discrete shock wave fitting procedure. Viscous and thermal diffusion may be included as source terms in the bicharacteristic method. A production type computer program capable of determining the flow field in a variety of axisymmetric mixed-compression supersonic inlets is available. The results of the present analysis agree well with those produced by the two-dimensional method of characteristics when axisymmetric flow fields are computed. For three-dimensional flow fields, the results of the present analysis agree well with experimental data except in regions of high viscous interaction and boundary layer removal. The present analysis does not compute the boundary layer, nor does it account for boundary layer bleed.

Vadyak, J.↗
NASA NTRS1987

Full potential integral solution for transonic flows with and without embedded Euler domains

Two methods are presented to solve for the transonic airfoil flow problems. The first method is based on the integral equation solution of the full-potential equation in terms of the velocity field, and a Shock Capturing-Shock Fitting (SCSF) scheme has been developed. The SCSF-scheme consists of a shock-capturing part and a shock-fitting part in which shock panels are introduced at the shock location. The sock panels are fitted and crossed by using the Rankine-Hugoniot relations. The second method is based on coupling the integral equation of the full-potential equation with the pseudo time integration of Euler equations in a small embedded region around the shock. The integral solution provides the initial and boundary conditions for the Euler domain. This scheme is named as the Integral Equation-Embedded Euler (IEEE) scheme. The two methods are applied to NACA 0012 and NACA 64A010A over a wide range of Mach numbers, and the reults are in good agreement with the experimental data and other computational results. The schemes converge within a number of iterations which is one-order of magnitude less than the finite-difference schemes.

Kandil, Osama A.↗
NASA NTRS1977

Computation of the inviscid supersonic flow about cones at large angles of attack by a floating discontinuity approach

The technique of floating shock fitting is adapted to the computation of the inviscid flowfield about circular cones in a supersonic free stream at angles of attack that exceed the cone half-angle. The resulting equations are applicable over the complete range of free-stream Mach numbers, angles of attack and cone half-angles for which the bow shock is attached. A finite difference algorithm is used to obtain the solution by an unsteady relaxation approach. The bow shock, embedded cross-flow shock, and vortical singularity in the leeward symmetry plane are treated as floating discontinuities in a fixed computational mesh. Where possible, the flowfield is partitioned into windward, shoulder, and leeward regions with each region computed separately to achieve maximum computational efficiency. An alternative shock fitting technique which treats the bow shock as a computational boundary is developed and compared with the floating-fitting approach. Several surface boundary condition schemes are also analyzed.

Daywitt, J.↗
NASA NTRS1996

Interaction of a shock with a longitudinal vortex

In this paper we study the shock/longitudinal vortex interaction problem in axisymmetric geometry. Linearized analysis for small vortex strength is performed, and compared with results from a high order axisymmetric shock-fitted Euler solution obtained for this purpose. It is confirmed that for weak vortices, predictions from linear theory agree well with results from nonlinear numerical simulations at the shock location. To handle very strong longitudinal vortices, which may ultimately break the shock, we use an axisymmetric high order essentially non-oscillatory (ENO) shock capturing scheme. Comparison of shock-captured and shock-fitted results are performed in their regions of common validity. We also study the vortex breakdown as a function of Mach number ranging from 1.3 to 10, thus extending the range of existing results. For vortex strengths above a critical value. a triple point forms on the shock and a secondary shock forms to provide the necessary deceleration so that the fluid velocity can adjust to downstream conditions at the shock.

Erlebacher, Gordon↗
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 NTRS1978

Calculation of the flow field in supersonic mixed-compression inlets at angle of attack using the three-dimensional method of characteristics with discrete shock wave fitting

The influence of molecular transport is included in the computation by treating viscous and thermal diffusion terms in the governing partial differential equations as correction terms in the method of characteristics scheme. The development of a production type computer program is reported which is capable of calculating the flow field in a variety of axisymmetric mixed-compression aircraft inlets. The results agreed well with those produced by the two-dimensional method characteristics when axisymmetric flow fields are computed. For three-dimensional flow fields, the results agree well with experimental data except in regions of high viscous interaction and boundary layer removal.

Vadyak, J.↗
NASA NTRS1987

Transonic airfoil computation using the integral equation with and without embedded Euler domains

Two transonic computational schemes which are based on the Integral Equation Formulation of the full potential equation were presented. The first scheme is a Shock Capturing-Shock Fitting (SCSF) scheme which uses the full potential equation throughout with the exception of the shock wave where the Rankine-Hugoniot relations are used to cross and fit the shock. The second scheme is an Integral Equation with Embedded Euler (IEEE) scheme which uses the full potential equation with an embedded region where the Euler equations are used. The two schemes are applied to several transonic airfoil flows and the results were compared with numerous computational results and experimental domains with fine grids. The SCSF-scheme is restricted to flows with weak shock, while the IEEE-scheme can handle strong shocks. Currently, the IEEE scheme is applied to other transonic flows with strong shocks as well as to unsteady pitching oscillations.

Kandil, Osama A.↗
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 NTRS1981

Implicit boundary conditions for the solution of the parabolized Navier-Stokes equations for supersonic flows

A fully implicit set of boundary conditions is developed for the solution of the parabolized Navier-Stokes equations for supersonic flow in two dimensions. Shock fitting is employed at the shock and the body has no-slip and specified temperature conditions. A specified heat transfer condition at the wall can be handled in a similar manner. In addition, the shock location is advanced in space in a fully implicit manner by utilizing the Rankine-Hugoniot conditions along with global conservation of mass.

Barnett, M.↗
NASA NTRS1989

Multidomain spectral solution of shock-turbulence interactions

The use of a fitted-shock multidomain spectral method for solving the time-dependent Euler equations of gasdynamics is described. The multidomain method allows short spatial scale features near the shock to be resolved throughout the calculation. Examples presented are of a shock-plane wave, shock-hot spot and shock-vortex street interaction.

Kopriva, David A.↗
NASA NTRS1976

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗
NASA NTRS1987

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗
NASA NTRS1989

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗
NASA NTRS1988

Calculation of unsteady Euler flows in turbomachinery using the linearized Euler equations

The model used by the present unsteady cascade flow calculation method is based on the linearized unsteady Euler equations and accounts for blade loading and geometry, shock motion, and wake motion. The steady flow that must be determined prior to the unsteady flow is ascertained by means of a Newton iteration procedure. A noteworthy feature of the procedure is the use of shock-fitting to determine steady and unsteady shock positions; the use of the Euler equations in conjunction with the Rankine-Hugoniot shock-jump conditions is found to directly model the generation of entropy and vorticity at shocks. Results are presented for both channel and cascade flows.

Hall, Kenneth C.↗
NASA NTRS1984

Spectral methods for compressible flow problems

Recent results concerning numerical simulation of shock waves using spectral methods are reviewed. Shock fitting techniques were discussed as well as shock capturing techniques with finite difference artificial viscosity. Also the notion of the information contained in the numerical results obtained by spectral methods and show how this information is recovered was discussed.

Gottlieb, D.↗
NASA NTRS1991

Computational methods for shock waves in three-dimensional supersonic flow

A central artificial-viscosity and an upwind-biased difference method are contrived to solve the Euler equations for flowfields over typical spacecrafts. The spatial discretization is based on either nodal or cell-vertex formulation in the domain extending from free stream to the end of the vehicle. The outer boundary is treated as a bow shock in the first method but is placed in the free stream in the second, which captures both bow and internal shocks using an approximate Riemann solver based on high-order extrapolation to the cell face. These methods were tested for the Shuttle and Hermes orbiters at wind-tunnel conditions and angles of attack ranging from 0 to 60 deg. The artificial-viscosity method incorporated with a shock-fitting procedure shows smeared crossflow and wing-shock positions and required 15 percent more CPU per node than the upwind method. Greater flexibility and robustness is demonstrated by the latter on a fixed grid for all cases considered.

Li, Chien-Peng↗
NASA NTRS1985

A piston-driven shock in the solar corona

A solar flare that occurred on the west limb at 20:38 UT on Mar. 25, 1981 generated a massive, rapidly expanding optical coronal transient, which moved outward with an approximately constant velocity of 800 km/s. An associated magnetohydrodynamic shock traveled out ahead of the transient with a velocity estimated to be approximately 1000 km/s. The optical and radio data on the transient and shock fit well with general theories concerning piston-driven shocks and with current MHD models for propagation of such shocks through the solar corona.

Maxwell, A.↗
NASA NTRS1993

Adjoint methods for aerodynamic wing design

A model inverse design problem is used to investigate the effect of flow discontinuities on the optimization process. The optimization involves finding the cross-sectional area distribution of a duct that produces velocities that closely match a targeted velocity distribution. Quasi-one-dimensional flow theory is used, and the target is chosen to have a shock wave in its distribution. The objective function which quantifies the difference between the targeted and calculated velocity distributions may become non-smooth due to the interaction between the shock and the discretization of the flowfield. This paper offers two techniques to resolve the resulting problems for the optimization algorithms. The first, shock-fitting, involves careful integration of the objective function through the shock wave. The second, coordinate straining with shock penalty, uses a coordinate transformation to align the calculated shock with the target and then adds a penalty proportional to the square of the distance between the shocks. The techniques are tested using several popular sensitivity and optimization methods, including finite-differences, and direct and adjoint discrete sensitivity methods. Two optimization strategies, Gauss-Newton and sequential quadratic programming (SQP), are used to drive the objective function to a minimum.

Grossman, Bernard↗