Engineering topics
Publications and source records attributed to Kerlick, G. D..
Supersonic laminar flow development in a constant-area square duct exhibits as one of its distinguishing features the formation of two contrarotating secondary flow vortices centered about the corner bisector. This phenomenon does not occur in unbounded corner flow. The secondary flow causes an outward bulging of total pressure contours in the vicinity of the corner bisector for wholly attached flow conditions.
Explore the source record for details and available documents.
Steady, developing, adiabatic supersonic flow in a square duct is investigated for an inlet Mach number of 3.91 and a unit Reynolds number of 1.8 x 10 to the 6th/m. The numerical results for laminar flow show that two secondary flow cells develop in the near vicinity of the corner which are centered about the corner bisector and distort the primary flow in this region. For turbulent flow, the experimental results indicate that two secondary flow cells also develop about the corner bisector, but are directed in an opposite sense to that observed for the laminar case. Numerical results based on the Baldwin-Lomax model show that this model is incapable of predicting turbulence-generated secondary flow cells. For a suitable choice of constants, the Gessner-Emery model is able to predict the strength of these cells, but is deficient with respect to predicting their positions in the flow and their distorting influence on the primary flow. These observations are based on comparisons made in this paper between predicted and measured total pressure contours, cross flow velocity profiles, and local wall shear stress distributions.
Solutions for the inviscid and viscous supersonic conical flow in a complete rectilinear inlet consisting of four planes intersecting at arbitrary wedge and sweep angles are obtained. To compute the flow on a specially constructed mesh, a three-dimensional flow solver ARC3D is used. It is shown that a single-pass mesh re-adaptation procedure can be used to obtain improved shock capture without the necessity of modifying the flow code, provided that the flow solver works in curvilinear coordinates and is restartable. Results computed by the method show good agreement with experimental measurements and previous calculations.
Previously cited in issue 06, p. 797, Accession no. A82-17812
Transonic Small Disturbance (TSD) theory is modified to yield steady and unsteady solutions in the case of strong shock waves (where the local Mach number ahead of the shock exceeds 1.3). The modification consists of an additional adjustable nonlinear term which allows the exact Rankine-Hugoniot shock jump relation to be satisfied at all times. The modified TSD theory is applied to steady and unsteady oscillatory transonic flows, and yields results which are in good agreement with solutions of the unsteady Euler equations, and which compare favorably with solutions of the full potential equation when shock waves are weak.
An algebraic decomposition of the Jacobian matrix which relates physical and computational variables is presented. This invertible decomposition parameterizes the mesh by the physically intuitive qualities of cell orientation, cell orthogonality, cell volume, and cell aspect ratio. The decomposition can be used to analyze numerically generated curvilinear coordinate meshes and to assess the contribution of the mesh to the truncation error for any specific differential operator and algorithm. This is worked out in detail for Laplace's equation in nonconservative and conservative forms. The analysis is applied to the solution of the full potential code TAIR, showing grid plots, carpet plots, and truncation error for a NACA 0012 airfoil.
A theory to correct the transonic small disturbance (TSD) equation to treat strong shock waves in unsteady flow is developed. The technique involves the addition of higher order terms, which are formally of negligible magnitude, to the low frequency TSD equation. These terms are then chosen such that any shock waves in the flow have strengths approximately equal to the appropriate Rankine-Hugoniot shock strength. Two correcting approaches are investigated. The first is to derive a correction for the mean steady flow and then simply use this corrected form for oscillatory flows. The second is to derive a correction for both steady and oscillatory parts of the flow. This second development is the most satisfactory and comparisons of the present results with Euler equation results are generally favorable, particularly regarding shock location, although there are some discrepancies in the pressure distribution in the leading edge region.