Engineering topics

Jameson, A.

Publications and source records attributed to Jameson, A..

At least 37 records · Page 2

NASA NTRS1981

Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes

A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.

Jameson, A.↗

NASA NTRS1981

An alternating-direction implicit algorithm for unsteady potential flow

An implicit finite-difference scheme is presented for the efficient computation of unsteady potential flow about airfoils. The formulation uses density and the velocity potential as dependent variables and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time-varying sheared-rectilinear coordinate transformation is employed. Calculated time-history solutions on a pulsating airfoil are compared with the results of other unsteady transonic codes, including a previous method of the authors. The present method is demonstrated to be unconditionally stable and to give accurate solutions with sharply resolved shocks.

Chipman, R.↗

NASA NTRS1981

Basic advances in the finite-volume method for transonic potential flow calculations

The finite-volume method of Jameson and Caughey provides a framework within which it is possible to calculate transonic potential flows about essentially arbitrary geometrical configurations. Improvements designed to increase the accuracy of the basic scheme and its consistency in the far field will be described. These include the incorporation of an artificial viscosity which maintains the formal second-order accuracy of the scheme in supersonic zones, and a modification of the flux balances to allow the free-stream conditions to satisfy the difference equations identically. Results of calculations illustrating the importance of these effects will be presented.

Caughey, D. A.↗

NASA NTRS1981

Accelerated finite-volume calculation of transonic potential flows

A fully conservative finite-volume algorithm is used to calculate transonic potential flows past isolated airfoils and through two-dimensional channels. The difference equations are solved by a multigrid technique which uses an alternating direction implicit method as a smoothing algorithm. The finite-volume formulation provides a framework within which to treat flows past complicated geometries, while the multigrid/alternating direction scheme provides rapid convergence of the solution to very small residuals.

Jameson, A.↗

NASA NTRS1980

Transonic flow past swept wings

FLO-22 computer program aids in numerical analysis of transonic potential flow past lifting, swept wing. FLO-22 uses relaxation method to solve finite-difference approximation of full-potential equation for transonic flow past configuration consisting of wing or arbitrary planform and dihedral extending from symmetry plane or wall. Comparisons of FLO-22 calculated results with experimental data for both conventional and super-critical transport wings show good agreement.

Caughey, D. A.↗

NASA NTRS1979

Recent progress in finite-volume calculations for wing-fuselage combinations

Progress in the application of finite-volume methods to the calculation of transonic potential flows past general wing-body combinations is reviewed. Two different methods of generating boundary-conforming grids are investigated, and the results compared to provide an estimate of solution sensitivity to grid geometry. Both conservative and quasi-conservative difference schemes are used in one of the coordinate systems. Results show that the error introduced by the quasi-conservative formulation seems to be small, although a one-dimensional analysis suggests that schemes of this type do not necessarily produce mass-conserving shocks. Comparison of calculated results with experimental data for realistic fuselage geometries clearly shows the importance of modelling the effect of fuselage geometry upon the wing pressure distribution.

Caughey, D. A.↗

NASA NTRS1979

Acceleration of transonic potential flow calculations on arbitrary meshes by the multiple grid method

A multiple grid method for transonic flow calculations is developed. The proposed scheme incorporates a generalized alternating direction method as the smoothing algorithm. Numerical experiments indicate that this multigrid alternating direction method converges rapidly and reliably for a range of cases typical of the cruising regime up to the onset of drag rise. It also appears that the method can be readily generalized to treat three-dimensional flows.

Jameson, A.↗

NASA NTRS1978

Recent experiences with three-dimensional transonic potential flow calculations

Some recent experiences with computer programs capable of solving finitie-difference approximations to the full potential equation for the transonic flow past three dimensional swept wings and simple wing-fuselage combinations are discussed. The programs used are a nonconservative program for swept wings, a quasi-conservative finite-volume program capable of treating swept wings mounted on fuselages of slowly varying circular cross section, and a fully conservative finite volume scheme capable of treating swept wings and wing-cylinder combinations. The present capabilities of these codes are reviewed. The relative merits of the conservative and nonconservative formulations are discussed, and the results of calculations including corrections for the boundary-layer displacement effect are presented.

Caughey, D. A.↗

NASA NTRS1978

Transonic flow about airfoils

Program analyzes airfoils that permit transonic flow for subsonic free-stream mach numbers. Transonic refers to aircraft speeds less than speed of sound, but close enough so that top of wing, where airflow is fastest, mach number becomes greater than 1. Program should aid design phase of new airfoil and in analysis of existing airfoils.

Bauer, F.↗

NASA NTRS1978

Preliminary study of the use of the STAR-100 computer for transonic flow calculations

A description is presented for a new explicit algorithm for solving the transonic small-disturbance potential equation on the STAR-100 computer. The STAR computer has a 'pipeline' type of processor which is very efficient in doing arithmetic operations on long vectors. Unfortunately, the successive line over-relaxation method, which represents the best method for solving the considered equation, is not suitable for vector-arithmetic procedures. The described algorithm, however, can be vectorized to make an employment of the STAR-100 possible. It is found that a use of the new algorithm for solving the equation on the STAR-100 computer can almost halve the computer time required for this type of computation when compared to successive line over-relaxation on the CYBER 175 computer.

Keller, J. D.↗

NASA NTRS1977

Numerical calculation of transonic flow past a swept wing by a finite volume method

The utility of numerical methods for predicting transonic flows over wings and bodies is well established. The computer program FLO22, based on a method presented earlier, has actually been widely used to calculate the aerodynamic performance of wings of transport aircraft. Provided that a correction is made for the displacement effect of the viscous boundary layer, this code has been found to give predictions which are accurate enough to serve as a useful design guide. The main disadvantages of the scheme used in FLO22 are the use of nonconservative difference formulas, which result in a failure to satisfy conservation of mass across shock waves, and the difficulty of finding suitable transformations of coordinates to permit the treatment of more complex geometric configurations. The method described here is an attempt to overcome these shortcomings, while retaining the successful features of the previous method. The basic idea is to use a discrete approximation which directly represents a balance of the mass flow through small volume elements. This leads to a relatively simple treatment of the potential flow equation in conservation form.

Jameson, A.↗

NASA NTRS1977

Preliminary study of the use of the STAR-100 computer for transonic flow calculations

An explicit method for solving the transonic small-disturbance potential equation is presented. This algorithm, which is suitable for the new vector-processor computers such as the CDC STAR-100, is compared to successive line over-relaxation (SLOR) on a simple test problem. The convergence rate of the explicit scheme is slower than that of SLOR, however, the efficiency of the explicit scheme on the STAR-100 computer is sufficient to overcome the slower convergence rate and allow an overall speedup compared to SLOR on the CYBER 175 computer.

Keller, J. D.↗

NASA NTRS1977

Numerical calculation of the transonic flow past a swept wing

A numerical method is presented for analyzing the transonic potential flow past a lifting, swept wing. A finite difference approximation to the full potential equation is solved in a coordinate system which is nearly conformally mapped from the physical space in planes parallel to the symmetry plane, and reduces the wing surface to a portion of one boundary of the computational grid. A coordinate invariant, rotated difference scheme is used, and the difference equations are solved by relaxation. The method is capable of treating wings of arbitrary planform and dihedral, although approximations in treating the tips and vortex sheet make its accuracy suspect for wings of small aspect ratio. Comparisons of calculated results with experimental data are shown for examples of both conventional and supercritical transport wings. Agreement is good for both types, but it was found necessary to account for the displacement effect of the boundary layer for the supercritical wing, presumably because of its greater sensitivity to changes in effective geometry.

Jameson, A.↗

NASA NTRS1977

Implicit approximate-factorization schemes for the efficient solution of steady transonic flow problems

Implicit approximate-factorization algorithms (AF) are developed for the solution of steady-state transonic flow problems. The performance of the AF solution method is evaluated relative to that of the standard solution method for transonic flow problems, successive line over-relaxation (SLOR). Both methods are applied to the solution of the nonlinear, two-dimensional transonic small-disturbance equation. Results indicate that the AF method requires substantially less computer time than SLOR to solve the nonlinear finite-difference matrix equation for a transonic flow field. This increase in computational efficiency is achieved with no appreciable increase in computer storage or coding complexity.

Ballhaus, W. F.↗

NASA NTRS1977

A finite volume method for transonic potential flow calculations

It is proposed to solve the exact transonic potential flow equation on a mesh constructed from small volume elements, which can be conveniently packed around any reasonably smooth configuration. The calculation is performed on two sets of interlocking cells. The velocity and density are calculated in the primary cells, and a flux balance is then established in the secondary cells. The scheme is desymmetrized by the addition of artificial viscosity in the supersonic zone. Some results are included for a swept wing and a wing-cylinder combination.

Jameson, A.↗

NASA NTRS1977

Calculation of transonic potential flowfields about complex, three-dimensional configurations

Methods for extending iterative, finite-difference calculations of transonic potential flowfields to complex three-dimensional configurations are discussed. One particularly attractive approach is to use relatively simple conformal mappings in combination with shearing transformations to generate computational domains that are nearly-conformally mapped from the physical space in one family of coordinate surfaces, and which map the complex boundaries to grid surfaces. The application of such a method to a general wing-body combination or to a multi-bladed fan is discussed. A transformation to map the wing-fuselage or fan-hub combination to a convenient computational domain is proposed. The transformation is useful in its own right for treating the two-dimensional problems of flow past a profile in a wind tunnel or through a cascade. Some results of preliminary calculations are presented.

Caughey, D. A.↗

NASA NTRS1976

A brief description of the Jameson-Caughey NYU transonic swept-wing computer program: FLO 22

A computer program for analyzing inviscid, isentropic, transonic flow past 3-D swept configurations is presented. Some basic aspects of the program are: (1) the free-stream Mach number is restricted only by the isentropic assumption; (2) weak shock waves are automatically located wherever they occur in the flow; (3) the finite-difference form of the full equation for the velocity potential is solved by the method of relaxation, after the flow exterior to the airfoil is mapped to the upper half plane; (4) the mapping procedure allows exact satisfaction of the boundary conditions and use of supersonic free stream velocities; (5) the finite difference operator is locally rotated in supersonic flow regions so as to properly account for the domain of dependence; and (6) the relaxation algorithm was stabilized using criteria from a time-like analogy.

Jameson, A.↗

NASA NTRS1976

Accelerated iterative calculation of transonic nacelle flowfields

A method is presented for the calculation of inviscid, supercritical flowfields about axisymmetric inlet cowls. A finite-difference calculation is performed in a simple, rectangular domain obtained from the nacelle geometry by a nearly-conformal mapping procedure. Type-dependent finite-differences are constructed using a coordinate-independent, 'rotated' differencing scheme. Methods of accelerating convergence of the iterative solution are demonstrated including a hybrid fast-Poisson-solver/relaxation scheme and an extrapolated relaxation procedure. Calculated pressure distributions are compared with experimental data for a variety of Mach numbers and mass-flow ratios, and show generally good agreement.

Caughey, D. A.↗