Three-dimensional flow simulations for supersonic mixed-compression inlets at incidence
Previously cited in issue 07, p. 965, Accession no. A82-19778
Engineering topics
Publications and source records attributed to Hoffman, J. D..
Previously cited in issue 07, p. 965, Accession no. A82-19778
A streamwise marching procedure, approximately 200 times faster than a full Navier-Stokes procedure with comparable accuracy, is presented for solving problems of compressible viscous subsonic flow. Results are presented and compared with experimental data for the cases of developing turbulent flow in a circular pipe; turbulent flow in a two-dimensional S-duct; and turbulent flow in a typical subsonic diffuser. Prior to each main marching step, a preliminary marching step is taken in which the integral continuity equation and an uncoupled form of the streamwise momentum equation are solved simultaneously to obtain the viscous pressure correction. During the main marching step the equations for continuity, streamwise momentum, cross-flow momentum, and energy are solved simultaneously as a coupled system using an implicit finite-difference method, with the viscous pressure correction treated as a source term. The analysis may be used for flows with both favorable and adverse pressure gradients and to predict the location of flow separation.
A computer program was developed which is capable of calculating the flow field in the supersonic portion of a mixed compression aircraft inlet operating at angle of attack. The supersonic core flow is computed using a second-order three dimensional method-of-characteristics algorithm. The bow shock and the internal shock train are treated discretely using a three dimensional shock fitting procedure. The boundary layer flows are computed using a second-order implicit finite difference method. The shock wave-boundary layer interaction is computed using an integral formulation. The general structure of the computer program is discussed, and a brief description of each subroutine is given. All program input parameters are defined, and a brief discussion on interpretation of the output is provided. A number of sample cases, complete with data listings, are provided.
A new streamwise marching procedure was developed and coded for compressible viscous subsonic flow in planar or axisymmetric ducts with or without centerbodies. The continuity, streamwise momentum, cross-flow momentum, and energy equations are written in generalized orthogonal curvilinear coordinates. To allow the use of a marching procedure, second derivatives in the streamwise momentum equation are written as the sum of a known two dimensional imposed pressure field and an unknown one dimensional viscous correction. For turbulent flow, the Reynolds stress and heat flux terms are modeled using two-layer eddy viscosity turbulence models.
The flow field in supersonic mixed compression aircraft inlets at angle of attack is calculated. A zonal modeling technique is employed to obtain the solution which divides the flow field into different computational regions. The computational regions consist of a supersonic core flow, boundary layer flows adjacent to both the forebody/centerbody and cowl contours, and flow in the shock wave boundary layer interaction regions. The zonal modeling analysis is described and some computational results are presented. The governing equations for the supersonic core flow form a hyperbolic system of partial differential equations. The equations for the characteristic surfaces and the compatibility equations applicable along these surfaces are derived. The characteristic surfaces are the stream surfaces, which are surfaces composed of streamlines, and the wave surfaces, which are surfaces tangent to a Mach conoid. The compatibility equations are expressed as directional derivatives along streamlines and bicharacteristics, which are the lines of tangency between a wave surface and a Mach conoid.
An analysis is presented for calculating the steady three-dimensional flow field in supersonic mixed-compression inlets at incidence. A zonal modeling approach is employed to obtain the solution. The supersonic core flow is computed using a second-order pentahedral bicharacteristic algorithm. The bow shock wave and the reflected internal shock train are determined using a three-dimensional discrete shock fitting procedure. The boundary layer flow adjacent to both the centerbody and the cowl is computed using a second-order implicit finite difference method. The flow in a shock wave-boundary layer interaction region is computed using an integral formulation. The culmination of the present research effort is the development of a production-type computer program capable of analyzing flow in a variety of mixed-compression aircraft inlets. Numerical results and experimental correlations are presented to illustrate application of the analysis.
Program uses method of characteristics for steady three-dimensional flow to calculate flow field in supersonic portion of mixed-compression aircraft inlet at non-zero angle of attack. Results agree well with experimental data except in regions of high viscous interaction. Flow field for variety of mixed-compression inlets can be calculated. Input includes geometry and attack of inlet. Output consists of list of parameters, solution planes, and description of shock waves. Program is written in FORTRAN IV for batch execution on CDC 6000-series.
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.
The results of weakly viscous flow analysis are presented. The flow field, including molecular transport, is computed with the aid of a bicharacteristic method. The bow shock wave and the internal shock wave are computed with the aid of a three-dimensional shock wave fitting procedure. Characteristic equations are presented, and numerical integration procedure is discussed. Here, an inverse marching scheme is employed in which the solution is obtained on space-like planes of constant x and on space curves defined by the intersections of the internal shock wave with the solid boundaries. The distance between solution planes is arrived at by the Courant-Friedrichs-Lewy stability criterion.
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.
The calculation procedure is based on the method of characteristics for steady three-dimensional flow. The bow shock wave and the internal shock wave system were computed using a discrete shock wave fitting procedure. The general structure of the computer program is discussed, and a brief description of each subroutine is given. All program input parameters are defined, and a brief discussion on interpretation of the output is provided. A number of sample cases, complete with data deck listings, are presented.
A second-order numerical method of characteristics based on a bicharacteristics scheme was developed for the evaluation of steady, supersonic, nonequilibrium, chemically reacting flows. The absolute accuracy and order of accuracy of the method was shown by comparisons with spherical source flow and axisymmetric, nonequilibrium flows. Numerical results were obtained for elliptical and super-elliptical nozzles. These results illustrate the complex nature of three-dimensional flows and the inadequacy of quasi-three-dimensional methods which neglect cross flows.
The accuracy and order of three bicharacteristic schemes based on the numerical method of characteristics for three-dimensional, steady, isentropic, supersonic flow were compared based on numerical calculations. Comparisons were made for source flows and for flows in contoured axisymmetric nozzles. The first two methods were tetrahedral line networks, differing only in the manner in which the cross derivatives were evaluated. The third method was the pentahedral line network proposed by Butler. All three methods were numerically stable and produced efficient and accurate solutions. A discussion of the relative merits and disadvantages of each method is presented.
Swirling flow equations in converging nozzles, comparing analytical and experimental data
Efficiency of conical thrust nozzles with various lengths and cone angles, determining optimization coefficient
Optimization analysis for axisymmetric rocket motor nozzle design based on assumptions for gas particle flow
Damping effects of ring-spoke baffles on acoustic wave instabilities in rocket combustion chamber
Optimum thrust nozzle contours for chemically reacting gas flows, obtaining set of partial differential equations for gas dynamic properties