Engineering topics

Melson, N. Duane

Publications and source records attributed to Melson, N. Duane.

At least 19 records

NASA NTRS1999

Plane Smoothers for Multiblock Grids: Computational Aspects

Standard multigrid methods are not well suited for problems with anisotropic discrete operators, which can occur, for example, on grids that are stretched in order to resolve a boundary layer. One of the most efficient approaches to yield robust methods is the combination of standard coarsening with alternating-direction plane relaxation in the three dimensions. However, this approach may be difficult to implement in codes with multiblock structured grids because there may be no natural definition of global lines or planes. This inherent obstacle limits the range of an implicit smoother to only the portion of the computational domain in the current block. This report studies in detail, both numerically and analytically, the behavior of blockwise plane smoothers in order to provide guidance to engineers who use block-structured grids. The results obtained so far show alternating-direction plane smoothers to be very robust, even on multiblock grids. In common computational fluid dynamics multiblock simulations, where the number of subdomains crossed by the line of a strong anisotropy is low (up to four), textbook multigrid convergence rates can be obtained with a small overlap of cells between neighboring blocks.

Llorente, Ignacio M.↗

NASA NTRS1998

Robust Multigrid Smoothers for Three Dimensional Elliptic Equations with Strong Anisotropies

We discuss the behavior of several plane relaxation methods as multigrid smoothers for the solution of a discrete anisotropic elliptic model problem on cell-centered grids. The methods compared are plane Jacobi with damping, plane Jacobi with partial damping, plane Gauss-Seidel, plane zebra Gauss-Seidel, and line Gauss-Seidel. Based on numerical experiments and local mode analysis, we compare the smoothing factor of the different methods in the presence of strong anisotropies. A four-color Gauss-Seidel method is found to have the best numerical and architectural properties of the methods considered in the present work. Although alternating direction plane relaxation schemes are simpler and more robust than other approaches, they are not currently used in industrial and production codes because they require the solution of a two-dimensional problem for each plane in each direction. We verify the theoretical predictions of Thole and Trottenberg that an exact solution of each plane is not necessary and that a single two-dimensional multigrid cycle gives the same result as an exact solution, in much less execution time. Parallelization of the two-dimensional multigrid cycles, the kernel of the three-dimensional implicit solver, is also discussed. Alternating-plane smoothers are found to be highly efficient multigrid smoothers for anisotropic elliptic problems.

Llorente, Ignacio M.↗

NASA NTRS1996

Seventh Copper Mountain Conference on Multigrid Methods

The Seventh Copper Mountain Conference on Multigrid Methods was held on April 2-7, 1995 at Copper Mountain, Colorado. This book is a collection of many of the papers presented at the conference and so represents the conference proceedings. NASA Langley graciously provided printing of this document so that all of the papers could be presented in a single forum. Each paper was reviewed by a member of the conference organizing committee under the coordination of the editors. The vibrancy and diversity in this field are amply expressed in these important papers, and the collection clearly shows the continuing rapid growth of the use of multigrid acceleration techniques.

Melson, N. Duane↗

NASA NTRS1996

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗

NASA NTRS1995

Efficiency and Accuracy of Time-Accurate Turbulent Navier-Stokes Computations

The accuracy and efficiency of two types of subiterations in both explicit and implicit Navier-Stokes codes are explored for unsteady laminar circular-cylinder flow and unsteady turbulent flow over an 18-percent-thick circular-arc (biconvex) airfoil. Grid and time-step studies are used to assess the numerical accuracy of the methods. Nonsubiterative time-stepping schemes and schemes with physical time subiterations are subject to time-step limitations in practice that are removed by pseudo time sub-iterations. Computations for the circular-arc airfoil indicate that a one-equation turbulence model predicts the unsteady separated flow better than an algebraic turbulence model; also, the hysteresis with Mach number of the self-excited unsteadiness due to shock and boundary-layer separation is well predicted.

Rumsey, Christopher L.↗

NASA NTRS1994

Multiblock, Multigrid Solution Of Euler Equations

Method of numerical solution of Euler equations of three-dimensional flow of compressible fluid involves combination of multiblock and multigrid strategies. In multiblock strategy, flow field divided, into multiple smaller, more computationally-convenient zones and computational grid fitted to applicable flow boundaries generated in each block. In multigrid strategy used here, different quantities computed, variously, on finer or coarser grids. Minimizing cost of computation by using fewest grid points yielding acceptably accurate values of affected variable. Multigrid strategy found effective in accelerating convergence to steady state, while multiblock strategy provides geometric flexibility.

Melson, N. Duane↗

NASA NTRS1993

The Sixth Copper Mountain Conference on Multigrid Methods, part 2

The Sixth Copper Mountain Conference on Multigrid Methods was held on April 4-9, 1993, at Copper Mountain, Colorado. This book is a collection of many of the papers presented at the conference and so represents the conference proceedings. NASA Langley graciously provided printing of this document so that all of the papers could be presented in a single forum. Each paper was reviewed by a member of the conference organizing committee under the coordination of the editors. The multigrid discipline continues to expand and mature, as is evident from these proceedings. The vibrancy in this field is amply expressed in these important papers, and the collection clearly shows its rapid trend to further diversity and depth.

Melson, N. Duane↗

NASA NTRS1993

Time-accurate Navier-Stokes calculations with multigrid acceleration

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.

Melson, N. Duane↗

NASA NTRS1993

The Sixth Copper Mountain Conference on Multigrid Methods, part 1

The Sixth Copper Mountain Conference on Multigrid Methods was held on 4-9 Apr. 1993, at Copper Mountain, CO. This book is a collection of many of the papers presented at the conference and as such represents the conference proceedings. NASA LaRC graciously provided printing of this document so that all of the papers could be presented in a single forum. Each paper was reviewed by a member of the conference organizing committee under the coordination of the editors. The multigrid discipline continues to expand and mature, as is evident from these proceedings. The vibrancy in this field is amply expressed in these important papers, and the collection clearly shows its rapid trend to further diversity and depth.

Melson, N. Duane↗

NASA NTRS1993

Analyzing Transonic Flow Over A Wing And Fuselage

TAWFIVE computer program calculates transonic flow over transport-type wing and fuselage. Combines multigrid acceleration technique in iterative solution of potential equation with use of integral-form boundary-layer equations to provide software tool computationally efficient, yet accurate enough for use in design. Simplifies solution process by breaking problem into loosely coupled set of modified equations. Inviscid-flow method, using standard inviscid-flow equations, valid in "outer" region away from wing, whereas boundary-layer equations valid in thin region near solid surface of wing. Written in FORTRAN 77.

Melson, N. Duane↗

NASA NTRS1992

Design and analysis of low boom concepts at Langley Research Center

The objective of the sonic boom research in the current High Speed Research Program is to ultimately make possible overland supersonic flight by a high speed civil transport. To accomplish this objective, it is felt that results in four areas must demonstrate that such a vehicle would be acceptable by the general public, by the airframers, and by the airlines. It should be demonstrated: (1) that some waveform shape has the possibility of being acceptable to the general public; (2) that the atmosphere would not totally destroy such a waveform during propagation; (3) that a viable airplane could be built which produces such a waveform; and (4) that any performance penalty suffered by a low boom aircraft would be counteracted by the economic benefit of overland supersonic flight. The work being done at LaRC is in support of the third element listed above--the area of configuration design. The initial part of the paper will give a review of the theory being used for configuration designs and discuss two theory validation models which were built and tested within the past two years. Discussion of the wind tunnel and theoretical results (linear theory and higher order methods) and their implications for future designs will be included.

Darden, Christine M.↗

NASA NTRS1992

Comparison of jet plume shape predictions and plume influence on sonic boom signature

An Euler shock-fitting marching code yields good agreement with semiempirically determined plume shapes, although the agreement decreases somewhat with increasing nozzle angle and the attendant increase in the nonisentropic nature of the flow. Some calculations for the low boom configuration with a simple engine indicated that, for flight at altitudes above 60,000 feet, the plume effect is dominant. This negates the advantages of a low boom design. At lower altitudes, plume effects are significant, but of the order that can be incorporated into the low boom design process.

Barger, Raymond L.↗

NASA NTRS1990

A multiblock multigrid method for the solution of the three-dimensional Euler equations

A general multiblock, multigrid method for the solution of the Euler equations has been developed. Two types of numerical methods were investigated, van Leer's flux-vector-splitting and Roe's flux-difference-splitting, with MUSCL type differencing used in both methods. An explicit two-step method and a multi-stage Runge-Kutta method have been tested. Results are presented for test cases of a channel flow, nozzle exhaust flow, and a transonic wing.

Cannizzaro, Frank E.↗

NASA NTRS1990

A multiblock multigrid three-dimensional Euler equation solver

Current aerodynamic designs are often quite complex (geometrically). Flexible computational tools are needed for the analysis of a wide range of configurations with both internal and external flows. In the past, geometrically dissimilar configurations required different analysis codes with different grid topologies in each. The duplicity of codes can be avoided with the use of a general multiblock formulation which can handle any grid topology. Rather than hard wiring the grid topology into the program, it is instead dictated by input to the program. In this work, the compressible Euler equations, written in a body-fitted finite-volume formulation, are solved using a pseudo-time-marching approach. Two upwind methods (van Leer's flux-vector-splitting and Roe's flux-differencing) were investigated. Two types of explicit solvers (a two-step predictor-corrector and a modified multistage Runge-Kutta) were used with multigrid acceleration to enhance convergence. A multiblock strategy is used to allow greater geometric flexibility. A report on simple explicit upwind schemes for solving compressible flows is included.

Cannizzaro, Frank E.↗

NASA NTRS1989

Updated users' guide for TAWFIVE with multigrid

A program for the Transonic Analysis of a Wing and Fuselage with Interacted Viscous Effects (TAWFIVE) was improved by the incorporation of multigrid and a method to specify lift coefficient rather than angle-of-attack. A finite volume full potential multigrid method is used to model the outer inviscid flow field. First order viscous effects are modeled by a 3-D integral boundary layer method. Both turbulent and laminar boundary layers are treated. Wake thickness effects are modeled using a 2-D strip method. A brief discussion of the engineering aspects of the program is given. The input, output, and use of the program are covered in detail. Sample results are given showing the effects of boundary layer corrections and the capability of the lift specification method.

Melson, N. Duane↗

NASA NTRS1988

Calculations of three-dimensional flows using the isenthalpic Euler equations with implicit flux-vector splitting

A numerical method for solving the isenthalpic form of the Euler equations is developed. The method is based on the concept of flux vector splitting in its implicit form applied to a cell centered finite volume scheme. Approximate factorization is implemented in solving the implicit part of the governing equations. Time marching to a steady state solution requires short computational times due to the relative efficiency of the basic method. Computational times are further reduced by the implementation of multigrid. Results for several basic cases are shown.

Cannizzaro, Frank E.↗

NASA NTRS1987

Multigrid acceleration of the isenthalpic form of the compressible flow equations

A numerical method for solving the isenthalpic form of the governing equations for compressible inviscid flows was developed. The method is based on the concept of flux vector splitting in its implicit form and was tested on several demanding configurations. Time marching to steady state was accelerated by the implementation of the multigrid procedure which very effectively increased the rate of convergence. High quality steady-state results were obtained for various test cases and required only short computational times due to the relative efficiency of the basic method.

Melson, N. Duane↗

NASA NTRS1987

Multigrid acceleration of the isenthalpic form of the compressible flow equations

Melson, N. Duane↗