Engineering Papers⌕ Search

Engineering topics

Caughey, D. A.

Publications and source records attributed to Caughey, D. A..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Distributed parallel processing applied to an implicit multigrid Euler/Navier-Stokes algorithm

An implicit multigrid algorithm for the solution of the Euler and Navier-Stokes equations has been implemented within the framework of multiple block-structured grids in which the physical domain is spatially decomposed into several blocks and the solution is advanced in parallel on each block. Utilities have been developed to implement such a scheme in a distributed computing environment. The multi-block algorithm is designed so that the explicit residual calculation is identical to that of single-block scheme, and therefore converged solutions for both schemes must be the same. To accelerate convergence, synchronous and asynchronous multigrid strategies are implemented. Significant speedups have been achieved in a multiple processor environment, while convergence rates similar to those of the single-block scheme are observed.

Tysinger, T. L.↗
NASA NTRS1991

Diagonal implicit multigrid solution of compressible turbulent flows

The Caughey (1988) multigrid diagonal implicit algorithm for Euler equations is extended to solve the two-dimensional thin layer Navier-Stokes equations for turbulent transonic flows. To further accelerate convergence to a steady state, the implicit scheme is used within the multigrid method. Results are presented for transonic flows past airfoils, and the flow-field results are compared with other computational data and experiments demonstrating the accuracy of the method.

Varma, R. R.↗
NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the three-dimensional unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalizing procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the three-dimensional Euler equations.

Yokota, J. W.↗
NASA NTRS1987

An L-U implicit multigrid algorithm for the three-dimensional Euler equations

An LU implicit multigrid scheme is developed for the calculation of three-dimensional transonic flow through rotating cascades. This numerical method solves the unsteady Euler equations of gas dynamics in a finite-volume form. The implicit scheme makes it possible to take a much larger time step than is normally permitted in most explicit schemes, while the multigrid method is incorporated to accelerate the convergence rate for steady state calculations. Using this method, computational storage requirements are comparable to those of explicit schemes, while operation counts are considerably less than those found in the more widely-used ADI schemes.

Yokota, Jeffrey W.↗
NASA NTRS1986

Aerodynamics of Engine-Airframe Interaction

The report describes progress in research directed towards the efficient solution of the inviscid Euler and Reynolds-averaged Navier-Stokes equations for transonic flows through engine inlets, and past complete aircraft configurations, with emphasis on the flowfields in the vicinity of engine inlets. The research focusses upon the development of solution-adaptive grid procedures for these problems, and the development of multi-grid algorithms in conjunction with both, implicit and explicit time-stepping schemes for the solution of three-dimensional problems. The work includes further development of mesh systems suitable for inlet and wing-fuselage-inlet geometries using a variational approach. Work during this reporting period concentrated upon two-dimensional problems, and has been in two general areas: (1) the development of solution-adaptive procedures to cluster the grid cells in regions of high (truncation) error;and (2) the development of a multigrid scheme for solution of the two-dimensional Euler equations using a diagonalized alternating direction implicit (ADI) smoothing algorithm.

Caughey, D. A.↗
NASA NTRS1986

Transonic potential flow in hyperbolic nozzles

The full potential equation for the classical problem of transonic flow through a hyperbolic nozzle (with or without a shock wave) is solved in conservation form using the finite volume method of Jameson and Caughey (1977). Either a firstor a second-order numerical viscosity is added in the direction of the flow, explicitly, in conservation form. A multigrid alternating direction implicit method is used to solve the difference equations, and the results obtained are compared with analytical and numerical results from previous researches.

Park, M.↗
NASA NTRS1986

Computation of transonic flow about helicopter rotor blades

An inviscid, nonconservative, three-dimensional full-potential flow code, ROT22, has been developed for computing the quasi-steady flow about a lifting rotor blade. The code is valid throughout the subsonic and transonic regime. Calculations from the code are compared with detailed laser velocimeter measurements made in the tip region of a nonlifting rotor at a tip Mach number of 0.95 and zero advance ratio. In addition, comparisons are made with chordwise surface pressure measurements obtained in a wind tunnel for a nonlifting rotor blade at transonic tip speeds at advance ratios from 0.40 to 0.50. The overall agreement between theoretical calculations and experiment is very good. A typical run on a CRAY X-MP computer requires about 30 CPU seconds for one rotor position at transonic tip speed.

Arieli, R.↗
NASA NTRS1985

Multi-grid calculation of transonic potential flows

The finite-volume method discussed by Jameson and Caughey (1977), and Caughey and Jameson (1979, 1980) has made it possible to calculate the transonic potential flow past any configuration for which a suitable boundary-conforming coordinate grid can be constructed. However, computations for practical three-dimensional problems have remained quite expensive in terms of the required computer time. The reason for this is primarily related to the large number of grid cells necessary for adequate resolution in these complex three-dimensional problems, taking into account the large number of iterations required to achieve even modest convergence on these fine grids. The present chapter provides a description of work directed at removing this latter difficulty by making use of the multigrid method. Attention is given to finite-volume formulation, multigrid iteration, geometrical aspects, and computed results.

Caughey, D. A.↗
NASA NTRS1984

An implict LU scheme for the Euler equations applied to arbitrary cascades

An implicit scheme for solving the Euler equations is derived and demonstrated. The alternating-direction implicit (ADI) technique is modified, using two implicit-operator factors corresponding to lower-block-diagonal (L) or upper-block-diagonal (U) algebraic systems which can be easily inverted. The resulting LU scheme is implemented in finite-volume mode and applied to 2D subsonic and transonic cascade flows with differing degrees of geometric complexity. The results are presented graphically and found to be in good agreement with those of other numerical and analytical approaches. The LU method is also 2.0-3.4 times faster than ADI, suggesting its value in calculating 3D problems.

Buratynski, E. K.↗
NASA NTRS1983

Comparison of calculated and measured pressures on straight and swept-tip model rotor blades

Using the quasi-steady, full potential code, ROT22, pressures were calculated on straight and swept tip model helicopter rotor blades at advance ratios of 0.40 and 0.45, and into the transonic tip speed range. The calculated pressures were compared with values measured in the tip regions of the model blades. Good agreement was found over a wide range of azimuth angles when the shocks on the blade were not too strong. However, strong shocks persisted longer than predicted by ROT22 when the blade was in the second quadrant. Since the unsteady flow effects present at high advance ratios primarily affect shock waves, the underprediction of shock strengths is attributed to the simplifying, quasi-steady, assumption made in ROT22.

Tauber, M. E.↗
NASA NTRS1983

Multi-grid calculation of three-dimensional transonic potential flows

A multi-grid algorithm has been developed to speed the iterative convergence of calculations for the transonic potential flow past swept wings and wing-fuselage combinations. The method is based upon a fully-conservative, finite-volume approximation to the steady potential equation which is second-order accurate everywhere in the flow field except near shock waves. The multi-grid scheme is incorporated within the framework of an alternating successive-line-overrelaxation (SLOR) solver of the difference equations. Computed results confirm the second-order accuracy of the scheme, and demonstrate the effectiveness of the multi-grid procedure.

Caughey, D. A.↗
NASA NTRS1981

On various treatments of potential equations at shocks

The potential equation that governs transonic inviscid flowfields is elliptic in subsonic regions and hyperbolic in supersonic regions. The transition from supersonic to subsonic flow may occur either continuously or across a surface of discontinuity known as a shock. Because of the changing type of governing equation and the presence of the surfaces of discontinuity, special treatments of the governing equations are needed at shocks to preserve mass conservation and to ensure satisfaction of the governing equation. In the so-called shock-capturing techniques, shocks are captured by smearing the discontinuity over several mesh cells in the computational domain. Most existing conservative schemes conserve mass flux throughout the flowfield, including shock boundaries; however, they also introduce zero-order errors in the approximation to the governing equation. It is shown that the zero-order errors, which do not diminish as the mesh spacing approaches zero, can cause discrepancies in the prediction of shock strength and location. These parameters are extremely sensitive features of transonic inviscid flowfield calculations. It is also shown that the zero-order errors can be avoided by developing higher-order schemes that properly model the governing equation at shocks.

Chen, L. T.↗
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.↗