Engineering topics
Publications and source records attributed to Jameson, A..
Some recent developments in numerical methods for calculating solutions to the transonic potential flow equation are reviewed, including (1) the construction of the stable coordinate independent difference schemes; (2) the use of conservation form to insure proper shock jump conditions; (3) analysis of the relaxation method by the time dependent analogy; (4) accelerated iterative schemes.
A review is presented of some recently developed numerical methods for the solution of nonlinear equations of mixed type. The methods considered use finite difference approximations to the differential equation. Central difference formulas are employed in the subsonic zone and upwind difference formulas are used in the supersonic zone. The relaxation method for the small disturbance equation is discussed and a description is given of difference schemes for the potential flow equation in quasi-linear form. Attention is also given to difference schemes for the potential flow equation in conservation form, the analysis of relaxation schemes by the time dependent analogy, the accelerated iterative method, and three-dimensional calculations.
The development of relaxation methods to calculate transonic flows is discussed. Rather accurate predictions can be made for a number of flows of interest using the transonic potential flow equation, which may be derived from the Euler equations for inviscid compressible flow by introducing the assumption that the flow is irrotational. Once the choice of a mathematical model has been settled, the numerical procedure for actually computing a solution contains two main elements: the construction of a discrete approximation which converges to the solution of the continuous problem in the limit as the mesh width is reduced to zero, and the solution of the resulting set of nonlinear difference equations by a convergent iterative scheme. The choice of an appropriate coordinate system and its influence on the accuracy of the discrete approximation is discussed. Several applications of the general method are described.
The use of a fast elliptic solver in combination with relaxation is presented as an effective way to accelerate the convergence of transonic flow calculations, particularly when a marching scheme can be used to treat the supersonic zone in the relaxation process.
A mathematical theory for the design and analysis of supercritical wing sections was previously presented. Examples and computer programs showing how this method works were included. The work on transonics is presented in a more definitive form. For design, a better model of the trailing edge is introduced which should eliminate a loss of fifteen or twenty percent in lift experienced with previous heavily aft loaded models, which is attributed to boundary layer separation. How drag creep can be reduced at off-design conditions is indicated. A rotated finite difference scheme is presented that enables the application of Murman's method of analysis in more or less arbitrary curvilinear coordinate systems. This allows the use of supersonic as well as subsonic free stream Mach numbers and to capture shock waves as far back on an airfoil as desired. Moreover, it leads to an effective three dimensional program for the computation of transonic flow past an oblique wing. In the case of two dimensional flow, the method is extended to take into account the displacement thickness computed by a semi-empirical turbulent boundary layer correction.
A method is presented for the solution of the full potential equation in conservation form. This assures that a weak solution satisfying proper isentropic jump conditions is obtained, giving an improved representation of shock waves in comparison with earlier nonconservative schemes. The method uses the concept of artificial viscosity to produce a stable difference scheme in the supersonic zone, and artificial time to generate a convergent iterative scheme for the solution of the difference equations. Results of calculations using the nonconservative and conservative schemes are compared.
The numerical aspect of theoretical work on transonics and supercritical wing sections are compiled. A model of the trailing edge is introduced which eliminates the loss of 15 to 20 percent experienced with heavily aft-loaded models, and it is indicated how drag creep can be reduced at off-design conditions. A rotated finite difference scheme is presented which can handle supersonic as well as subsonic free stream Mach numbers and leads to an effective three-dimensional program for the computation of transonic flow past an oblique wing. In the case of two-dimensional flow, the method is extended to take into account the displacement thickness computed by a semiempirical turbulent boundary layer correction. A series of supercritical wing sections is discussed together with comparisons between experimental and theoretical data. Computer programs and a brief manual for their operation are listed. It is shown that the programs furnish a physically adequate computer simulation of the compressible flows that arise in problems of transonic aerodynamics.
A new method of calculating transonic flows based on a 'rotated' difference scheme is described. It is suitable for the calculation of both two- and three-dimensional flows without restriction on the speed at infinity and is well adapted to computer use. The Murman procedure is modified to eliminate any assumptions about the direction of flow when constructing the difference scheme. The proper directional property is obtained by rotating the difference scheme to conform with the local stream direction. In the hyperbolic region retarded difference formulas are used for all contributions to the streamwise second derivative, producing a correctly oriented positive artificial viscosity. In the absence of a simple implicit scheme in the hyperbolic and elliptic regions, the concept of iterations as steps in artificial time is introduced. Computer testing of this procedure provides numerical confirmation of the existence and uniqueness of weak solutions of the potential equation when a suitable entropy inequality is enforced.
The present work describes a mathematical model and numerical scheme for the computer-aided calculation of two- and three-dimensional transonic flow over an isolated yawed wing with oblique shock waves and a trailing vortex sheet. The flow is modeled by the potential equation for irrotational flow which is hyperbolic at supersonic points and elliptic at subsonic points. A coordinate-invariant difference scheme is used in which retarded difference formulas are constructed to conform with the local flow direction. The resulting 'rotated' difference scheme allows complete flexibility in the choice of a coordinate system. Shock waves are located automatically in the form of compression bands spread over a few mesh widths. The scheme has proven to be stable and convergent throughout the transonic range. Calculations have been performed for Mach number up to 1.2 and yaw angles up to 60 deg, the likely operating range of a yawed-wing transport designed to fly at supersonic speeds. Calculations becomes less accurate towards upper end of range, because the difference scheme is first-order accurate in the supersonic range.
A finite-difference relaxation method is presented for numerical solution of the full potential equation and exact boundary conditions for general axisymmetric bodies is inviscid, steady transonic flow. Body-normal coordinates are used in the nose region and sheared cylindrical coordinates are used on the afterbody to accommodate corners such as boattails and flares. An improved difference scheme is used which does not require that the flow be nearly alined with a coordinate direction in supersonic regions, and which treats either subsonic or supersonic free streams. Numerical results are illustrated for some simple classical shapes such as spheres and ellipsoids, and for more practical shapes like tangent-ogives with boattails. Special attention is given to bodies which have been studied for area-rule applications. Agreement with available experimental results is good in cases where viscous effects and wind-tunnel wall interference are not important.
Results are presented of calculations of the three dimensional steady transonic flow over a finite yawed wing. The full potential flow equation is solved in a transformed coordinate system which permits the boundary conditions to be satisfied exactly. The correct differential properties are enforced by rotating the difference scheme to conform with the flow direction, and fast convergence is assured by simulating a time dependent equation designed to settle quickly to a steady state. Computed lift drag ratios are consistent with the results of wind tunnel tests of a yawed wing conducted by R. T. Jones (1972).-
Computer programs predict effect of slipstream-wing flow interaction on aerodynamic characteristics of deflected slipstream and tilt aircraft. One program calculates lift distribution, lift, and drag of wing in wide slipstream. Results permit development of simplified lifting surface theory for circular jet.
Theoretical analysis of wing lift in wide and circular stripstreams, and engineering method for predicting practical V/STOL configuration characteristics
Numerical analysis of propeller-wing flow interaction with slipstreams