Engineering Papers⌕ Search

Engineering topics

Thompson, J. F.

Publications and source records attributed to Thompson, J. F..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Adaptive EAGLE dynamic solution adaptation and grid quality enhancement

In the effort described here, the elliptic grid generation procedure in the EAGLE grid code was separated from the main code into a subroutine, and a new subroutine which evaluates several grid quality measures at each grid point was added. The elliptic grid routine can now be called, either by a computational fluid dynamics (CFD) code to generate a new adaptive grid based on flow variables and quality measures through multiple adaptation, or by the EAGLE main code to generate a grid based on quality measure variables through static adaptation. Arrays of flow variables can be read into the EAGLE grid code for use in static adaptation as well. These major changes in the EAGLE adaptive grid system make it easier to convert any CFD code that operates on a block-structured grid (or single-block grid) into a multiple adaptive code.

Luong, Phu Vinh↗
NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

Efforts in transferring computational work from the LRC computer to the IRIS Graphics Workstation at MSU are reported and the computation of a conservative solution of a simple hyperbolic equation on an overlapping grid is discussed. Several conclusions concerning computations on overlapping grids are apparent. Problems only occur when there is a major difference in grid spacing on the individual component grids. In the case of hyperbolic equations, it is necessary that both interpolation and extrapolation be applied at the grid boundaries. When interpolated values are used at outflow boundary points, excessive oscillations in the numerical solution may be the result. The same conclusions would be valid for more complicated systems of hyperbolic equations such as the Euler equations for inviscid flow. Some of the solution values would be extrapolated at the overlap boundary, the exact number depending on the number of characteristics pointing out of the overlap region. It is also possible that similar boundary conditions may be needed for some parabolic equations such as high Reynolds number viscous flow equations. Efforts were also expended on the development of three-dimensional conservative interpolation procedures. Finally, the investigation of grid smoothing procedures were initiated during this reporting period. It was decided that the first grid smoothing algorithms will be based on the concepts of variational grid generation.

Mastin, C. W.↗
NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

A natural grid is defined on any parameterized curve or surface by selecting an equi-spaced set of parameter values. Redistributing the grid points can be accomplished by defining a new parameterization. Reparameterization techniques are introduced and applied in the construction of computational grids.

Mastin, C. W.↗
NASA NTRS1985

A Survey of Composite Grid Generation for General Three-dimensional Sections

The generation and use of composite grids for general three-dimensional physical boundary configurations is discussed, and the availability of several codes or procedures is noted. With the composite framework, the physical region is segmented into sub-regions, each bounded by six curved sides, and a grid is generated in each sub-region. These grids may be joined at the interfaces between the sub-regions with various degrees of continuity. This structure allows codes to be constructed to operate on rectangular blocks in computational space, so that existing solution procedures can be readily incorporated in the construction of codes for general configurations. Numerical grid generation is an integral part of the numerical solution of partial differential equations and is one of the pacing items in the development of codes for general configurations. The numerically generated grid frees the computational simulation from restriction to certain boundary shapes and allows general codes to be written in which the boundary shape is specified simply by input. The numerically generated grid allows all computation to be done on a fixed square grid in the computational space, which is always rectangular by construction.

Thompson, J. F.↗
NASA NTRS1984

Transformation of two and three-dimensional regions by elliptic systems

The research during this period continued to expand the class of numerical algorithms that can be accurately and efficiently implemented on overlapping grids. Whereas previous calculations have been used to solve elliptic equations and to find the steady-state solution of parabolic equations, the present work is aimed towards developing time-accurate solution techniques for parabolic and hyperbolic equations. The primary difficulty here is in the correct treatment of the interior boundary nodes that must be updated at each iteration. The implementation of explicit methods is straightforward. However, the common practice of lagging these values when using an implicit methods leads to inconsistencies in the difference equation. One way to avoid this problem is to alternately calculate with an implicit and an explicit method on each subgrid. With this procedure, the explicit method generates boundary values at the next time level which are then used by the implicit step. It can be shown that when a backward implicit method is combined with a forward explicit method, the composite method is second order accurate and unconditionally stable for linear problems. A second area in which progress can be reported is in the distribution of grid points on curves and surfaces.

Mastin, E. W.↗
NASA NTRS1984

A vectorized solution for incompressible flow

An algorithm is developed to obtain solutions to the unsteady Reynolds-averaged incompressible Navier-Stokes equations in general curvilinear coordinates on a vector processor. The governing equations are in nonconservative form with the velocity and pressure as dependent variables. Two momentum equations and the Poisson equation for pressure form a set of three governing equations for three flow field unknowns: u, v, and p. The governing equations and boundary conditions are expressed in terms of boundary-conforming curvilinear coordinates, and a checkerboard SOR iteration is used to solve the governing equations. Several possible sequences for a checkerboard SOR iteration are investigated for finding the best overall convergence rate. The efficiency and capability of the present algorithm was assessed using the example of an 18 percent thick NACA 66(3)018 airfoil at zero degree angle of attack for chord Reynolds number range 1000-40,000.

Patel, N. R.↗
NASA NTRS1983

Transformation of two and three-dimensional regions by elliptic systems

Finite difference methods for composite grids were analyzed. It was observed that linear interpolation between grids would suffice only where low order accuracy was required. In the context of fluid flow, this would be in regions where the flow was essentially free stream. Higher order interpolation schemes were also investigated. The well known quadratic and cubic interpolating polynomials would increase the formal accuracy of the overall numerical algorithm. However, it can also be shown that the stability of the algorithm may be adversely affected. Further numerical results are needed in order to assess the nature of this instability induced by the interpolation procedure. Finally, error analysis and the order of difference expressions on general curvilinear coordinates are discussed.

Thompson, J. F.↗
NASA NTRS1983

Numerical solution of the Navier-Stokes equations for arbitrary 2-dimensional multi-element airfoils

Numerical solutions of the Navier-Stokes equations, with an algebraic turbulence model, for time-dependent two dimensional flow about multi-element airfoils were developed. Fundamental to these solutions was the use of numerically-generated boundary-conforming curvilinear coordinate systems to allow bodies of arbitrary shape to be treated. A general two dimensional grid generation code for multiple-body configuration was written as a part of this project and made available through the COSMIC code library.

Thompson, J. F.↗
NASA NTRS1983

A survey of grid generation techniques in computational fluid dynamics

The use of numerically-generated body-conforming coordinate systems has become well-established in computational fluid dynamics, allowing general codes to be written which can treat boundaries of arbitrary shape. There is now a rapidly increasing interest in the dynamic coupling of the grid with the physical solution so that the grid lines continually move to achieve concentration in regions of strong variation of the physical solution and alignment with shocks, flame fronts, etc. Such adaptive coordinate systems have been developed using a variety of techniques for sensing and reacting to the developing regions of concentration and/or alignment. This paper reviews the techniques and applications of these dynamically adaptive systems and other systems.

Thompson, J. F.↗
NASA NTRS1983

Adaptive grids generated by elliptic systems

It is pointed out that a finite difference grid which moves with the solution of the partial differential equation being solved can improve the accuracy and efficiency of a numerical algorithm. This technique is particularly advantageous in the solution of problems involving boundary layers or shocks where a poorly chosen grid may give a numerical solution which is useless because of poor resolution or extreme oscillations. The present investigation is concerned with the development of a scheme which does not excessively distort the grid. The grid generation algorithm is based on the numerical solution of a system of elliptic differential equations. Holst and Brown (1981) have used a preliminary solution to move points on the boundary of the physical region and then resolved the problem on a new grid generated by an elliptic system. In the current investigation, the grid movement and the solution will develop simultaneously. The solution is used to modify the generating equations, thereby controlling the grid point distributions.

Mastin, C. W.↗
NASA NTRS1983

Quasiconformal mappings and grid generation

A finite difference scheme is developed for constructing quasiconformal mappings for arbitrary simply and doubly connected regions. Computational grids are generated to reduce elliptic equations to canonical form. Examples of conformal mappings on surfaces are also included.

Mastin, C. W.↗
NASA NTRS1982

Numerical Solution for Navier-Stokes Equations

Carefully selected blend of computational techniques solves complete set of equations for viscous, unsteady, hypersonic flow in general curvilinear coordinates. New algorithm has tested computation of axially directed flow about blunt body having shape similar to that of such practical bodies as wide-body aircraft or artillery shells. Method offers significant computational advantages because of conservation-law form of equations and because it reduces amount of metric data required.

Warsi, Z. U. A.↗
NASA NTRS1982

Boundary-fitted coordinate systems for numerical solution of partial differential equations - A review

A comprehensive review of methods of numerically generating curvilinear coordinate systems with coordinate lines coincident with all boundary segments is given. Some general mathematical framework and error analysis common to such coordinate systems is also included. The general categories of generating systems are those based on conformal mapping, orthogonal systems, nearly orthogonal systems, systems produced as the solution of elliptic and hyperbolic partial differential equations, and systems generated algebraically by interpolation among the boundaries. Also covered are the control of coordinate line spacing by functions embedded in the partial differential operators of the generating system and by subsequent stretching transformation. Dynamically adaptive coordinate systems, coupled with the physical solution, and time-dependent systems that follow moving boundaries are treated. References reporting experience using such coordinate systems are reviewed as well as those covering the system development.

Thompson, J. F.↗
NASA NTRS1982

Approximate factorization with an elliptic pressure solver for incompressible flow

Two-dimensional curvilinear coordinates are used to solve the incompressible Navier-Stokes equations, in conjunction with approximate factorization for the solution of the momentum equation and the successive overrelaxation by lines method for the solution of a Poisson equation for the pressure. The combined algorithm, although not fully explicit, is marginally stable at Reynolds numbers lower than 10,000 and time increments of 0.01. Pressure distributions calculated for attack angles of zero and 6 deg are of the same shape as the experimental curves, but are shifted to one side.

Bernard, R. S.↗
NASA NTRS1982

Numerical solution of potential flow about arbitrary 2-dimensional multiple bodies

A procedure for the finite-difference numerical solution of the lifting potential flow about any number of arbitrarily shaped bodies is given. The solution is based on a technique of automatic numerical generation of a curvilinear coordinate system having coordinate lines coincident with the contours of all bodies in the field, regardless of their shapes and number. The effects of all numerical parameters involved are analyzed and appropriate values are recommended. Comparisons with analytic solutions for single Karman-Trefftz airfoils and a circular cylinder pair show excellent agreement. The technique of application of the boundary-fitted coordinate systems to the numerical solution of partial differential equations is illustrated.

Thompson, J. F.↗
NASA NTRS1982

Numerical grid generation; Proceedings of the Symposium on Numerical Generation of Curvilinear Coordinate Systems and Their Use in the Numerical Solution of Partial Differential Equations, Nashville, TN, April 13-16, 1982

General curvilinear coordinate systems are considered along with the error induced by coordinate systems, basic differential models for coordinate generation, elliptic grid generation, conformal grid generation, algebraic grid generation, orthogonal grid generation, patched coordinate systems, and solid mechanics applications of boundary fitted coordinate systems. Attention is given to coordinate system control and adaptive meshes, the application of body conforming curvilinear grids for finite difference solution of external flow, the use of solution adaptive grids in solving partial differential equations, adaptive gridding for finite difference solutions to heat and mass transfer problems, and the application of curvilinear coordinate generation techniques to the computation of internal flows. Other topics explored are related to the solution of nonlinear water wave problems using boundary-fitted coordinate systems, the numerical modeling of estuarine hydrodynamics on a boundary-fitted coordinate system, and conformal grid generation for multielement airfoils.

Thompson, J. F.↗