Application of the implicit MacCormack scheme to the parabolized Navier-Stokes equations
Explore the source record for details and available documents.
Lawrence, S. L.↗
Engineering topics
Tannehill, J. C.
Publications and source records attributed to Tannehill, J. C..
Numerical computation of space shuttle orbiter flow field
The complete inviscid viscous real gas flow around the Space Shuttle Orbiter was computed. Real gas effects are important in predicting the reentry environment around the Orbiter because the high temperatures within the shock layer cause the air to dissociate and ionize, thus invalidating the perfect gas assumption. It is shown that real gas effects has a significant influence on the Orbiter aerodynamics. The approach utilizes a time dependent Navier-Stokes code to compute the subsonic nose portion of the flow field. This nose solution provides the initial conditions for the parabolized Navier-Stokes (PNS) code.
Tannehill, J. C.↗
Numerical solution of Space Shuttle Orbiter flow field including real gas effects
The hypersonic, laminar flow around the Space Shuttle Orbiter has been computed for both an ideal gas (gamma = 1.2) and equilibrium air using a real-gas, parabolized Navier-Stokes code. This code employs a generalized coordinate transformation; hence, it places no restrictions on the orientation of the solution surfaces. The initial solution in the nose region was computed using a 3-D, real-gas, time-dependent Navier-Stokes code. The thermodynamic and transport properties of equilibrium air were obtained from either approximate curve fits or a table look-up procedure. Numerical results are presented for flight conditions corresponding to the STS-3 trajectory. The computed surface pressures and convective heating rates are compared with data from the STS-3 flight.
Prabhu, D. K.↗
Application of the implicit MacCormack scheme to the parabolized Navier-Stokes equations
MacCormack's implicit finite-difference scheme was used to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method for solving the PNS equations does not require the inversion of block tridiagonal systems of algebraic equations and permits the original explicit MacCormack scheme to be employed in those regions where implicit treatment is not needed. The advantages and disadvantages of the present adaptation are discussed in relation to those of the conventional Beam-Warming scheme for a flat plate boundary layer test case. Comparisons are made for accuracy, stability, computer time, computer storage, and ease of implementation. The present method was also applied to a second test case of hypersonic laminar flow over a 15% compression corner. The computed results compare favorably with experiment and a numerical solution of the complete Navier-Stokes equations.
Lawrence, J. L.↗
Numerical computation of viscous flow around bodies and wings moving at supersonic speeds
Research in aerodynamics is discussed. The development of equilibrium air curve fits; computation of hypersonic rarefield leading edge flows; computation of 2-D and 3-D blunt body laminar flows with an impinging shock; development of a two-dimensional or axisymmetric real gas blunt body code; a study of an over-relaxation procedure forthe MacCormack finite-difference scheme; computation of 2-D blunt body turbulent flows with an impinging shock; computation of supersonic viscous flow over delta wings at high angles of attack; and computation of the Space Shuttle Orbiter flowfield are discussed.
Tannehill, J. C.↗
Numerical solution of Space Shuttle Orbiter flowfield
Previously cited in issue 12, p. 1850, Accession no. A82-27083
Rakich, J. V.↗
Application of the implicit MacCormack scheme to the PNS equations
The two-dimensional parabolized Navier-Stokes equations are solved using MacCormack's (1981) implicit finite-difference scheme. It is shown that this method for solving the parabolized Navier-Stokes equations does not require the inversion of block tridiagonal systems of algebraic equations and allows the original explicit scheme to be employed in those regions where implicit treatment is not needed. The finite-difference algorithm is discussed and the computational results for two laminar test cases are presented. Results obtained using this method for the case of a flat plate boundary layer are compared with those obtained using the conventional Beam-Warming scheme, as well as those obtained from a boundary layer code. The computed results for a more severe test of the method, the hypersonic flow past a 15 deg compression corner, are found to compare favorably with experiment and a numerical solution of the complete Navier-Stokes equations.
Lawrence, S. L.↗
Numerical solution of Space Shuttle Orbiter flow field
The supersonic, viscous laminar flow around the Space Shuttle Orbiter forebody has been computed with a parabolized Navier-Stokes code using a generalized coordinate transformation. The initial solution for the nose part of the Orbiter geometry was obtained with a three-dimensional time-dependent Navier-Stokes solver. It was necessary to employ a wind axis oriented coordinate system to obtain the initial solution with the time-dependent code. The generalized PNS technique used in this study allows the solution to be marched from the given initial data surface to any desired surface downstream. A grid point clustering scheme was employed to accurately describe the body shape by clustering points at the wing tip and at the wing body juncture. The computed heat transfer coefficients, pressure coefficients, and shock shapes are compared with the available experimental data for 0 degrees and 30 degrees angle of attack.
Venkatapathy, E.↗
Numerical solution of supersonic viscous flow over blunt delta wings
A general parabolized Navier-Stokes code has been developed to compute the steady supersonic viscous flow around arbitrary body shapes at high angles of attack. A nonorthogonal three-dimensional coordinate frame permits the code to march with solution surfaces which are the most appropriate. The code has been used to calculate the laminar flow over a slab delta wing with 70 deg sweep at angles of attack up to 41.5 deg and Mach numbers of 6.8 and 9.6. The computed shock shapes, surface pressures and heat transfer coefficients are compared with experiment and show agreement.
Tannehill, J. C.↗
Calculation of supersonic viscous flow over delta wings with sharp subsonic leading edges
Two complementary procedures have been developed to calculate the viscous supersonic flow over conical shapes at large angles of attack, with application to cones and delta wings. In the first approach the flow is assumed to be conical and the governing equations are solved at a given Reynolds number with a time-marching explicit finite-difference algorithm. In the second method the parabolized Navier-Stokes equations are solved with a space-marching implicit noniterative finite-difference algorithm. This latter approach is not restricted to conical shapes and provides a large improvement in computational efficiency over published methods. Results from the two procedures agree very well with each other and with available experimental data.
Vigneron, Y. C.↗
Numerical solution of two-dimensional turbulent blunt body flows with an impinging shock
An implicit finite-difference method has been developed to compute two-dimensional, turbulent, blunt body flows with an impinging shock wave. The full time-averaged Navier-Stokes equations are solved with algebraic eddy viscosity and turbulent Prandtl number models employed for shear stress and heat flux. The irregular-shaped bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied. A Type III turbulent shock interference flow field has been computed and the numerical results compare favorably with existing experimental data. In addition, comparisons are made between the present implicit code and a previous explicit code.
Tannehill, J. C.↗
Calculation of supersonic viscous flow over delta wings with sharp subsonic leading edges
Two complementary procedures were developed to calculate the viscous supersonic flow over conical shapes at large angles of attack, with application to cones and delta wings. In the first approach the flow is assumed to be conical and the governing equations are solved at a given Reynolds number with a time-marching explicit finite-difference algorithm. In the second method the parabolized Navier-Stokes equations are solved with a space-marching implicit noniterative finite-difference algorithm. This latter approach is not restricted to conical shapes and provides a large improvement in computational efficiency over published methods. Results from the two procedures agree very well with each other and with available experimental data.
Vigneron, Y. C.↗
On improving the iterative convergence properties of an implicit approximate-factorization finite difference algorithm
The iterative convergence properties of an approximate-factorization implicit finite-difference algorithm are analyzed both theoretically and numerically. Modifications to the base algorithm were made to remove the inconsistency in the original implementation of artificial dissipation. In this way, the steady-state solution became independent of the time-step, and much larger time-steps can be used stably. To accelerate the iterative convergence, large time-steps and a cyclic sequence of time-steps were used. For a model transonic flow problem governed by the Euler equations, convergence was achieved with 10 times fewer time-steps using the modified differencing scheme. A particular form of instability due to variable coefficients is also analyzed.
Desideri, J. A.↗
Navier-Stokes calculations for laminar and turbulent hypersonic flow over indented nosetips
A time-accurate finite-difference Navier-Stokes code has been used to calculate the viscous flow over a severely indented blunt body in a supersonic stream. An algebraic turbulence model is used and the results are compared with experimental data from wind-tunnel tests. Qualitative agreement is obtained for the surface pressure distribution and flow-field structure, including the separated bubble in the indented region. However, uncertainties still exist in the heating calculations, which are attributed to the turbulence model. For both laminar and turbulent calculations, the flow exhibits a fundamental unsteady character at a frequency of about 50 kHz.
Rakich, J. V.↗
Over-relaxation applied to the MacCormack finite-difference scheme
An over-relaxation procedure is applied to the MacCormack finite-difference scheme in order to reduce the computation time required to obtain a steady-state solution. The implementation of this acceleration procedure to an existing computer program using the regular MacCormack method is extremely simple and does not require additional storage. The over-relaxation procedure does not alter the steady-state solution, which is second-order accurate. The method is first applied to Burgers' equation. A stability condition and an expression for the increase in the rate of convergence are derived. The method is then applied to the calculation of the hypersonic viscous flow over a flat plate, using the complete Navier-Stokes equations, and the inviscid flow over a wedge. Reductions in computing time by factors of 3 and 1.5, respectively, are obtained by over-relaxation.
Desideri, J.-A.↗
Comparison of a two-dimensional shock impingement computation with experiment
Results of computations of two-dimensional viscous blunt-body flowfields with an impinging shock wave, with a time-dependent finite-difference method employed to solve the complete set of Navier-Stokes equations, are compared with experimental results. The experimental results were obtained in a 20-inch hypersonic tunnel with a planar shock impinging on the cylindrical leading edge of a fin, hence with the shock parallel to the centerline of the leading edge, so that type III and type IV interference patterns were generated. Close agreement is found. The overall effects of smoothing and grid size on the calculations are determined. A 31 x 51 mesh is adequate for wall pressure values (except in peaked regions).
Tannehill, J. C.↗
Numerical computation of three-dimensional blunt body flow fields with an impinging shock
A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.
Holst, T. L.↗
Improved curve fits for the thermodynamic properties of equilibrium air suitable for numerical computation using time-dependent or shock-capturing methods, part 1
Simplified curve fits for the thermodynamic properties of equilibrium air were devised for use in either the time-dependent or shock-capturing computational methods. For the time-dependent method, curve fits were developed for p = p(e, rho), a = a(e, rho), and T = T(e, rho). For the shock-capturing method, curve fits were developed for h = h(p, rho) and T = T(p, rho). The ranges of validity for these curves fits were for temperatures up to 25,000 K and densities from 10 to the minus 7th power to 10 to the 3d power amagats. These approximate curve fits are considered particularly useful when employed on advanced computers such as the Burroughs ILLIAC 4 or the CDC STAR.
Tannehill, J. C.↗