Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multigrid methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Multigrid methods for numerical simulation of laminar diffusion flames

This paper documents the result of a computational study of multigrid methods for numerical simulation of 2D diffusion flames. The focus is on a simplified combustion model, which is assumed to be a single step, infinitely fast and irreversible chemical reaction with five species (C3H8, O2, N2, CO2 and H2O). A fully-implicit second-order hybrid scheme is developed on a staggered grid, which is stretched in the streamwise coordinate direction. A full approximation multigrid scheme (FAS) based on line distributive relaxation is developed as a fast solver for the algebraic equations arising at each time step. Convergence of the process for the simplified model problem is more than two-orders of magnitude faster than other iterative methods, and the computational results show good grid convergence, with second-order accuracy, as well as qualitatively agreement with the results of other researchers.

Liu, C.↗
NASA NTRS1996

Some Aspects of Multigrid Methods on Non-Structured Meshes

To solve a given fine mesh problem, the design of a multigrid method requires the definition of coarse levels, associated coarse grid operators and inter-grid transfer operators. For non-structured simplified meshes, these definitions can rely on the use of non-nested triangulations. These definitions can also be founded on agglomeration/aggregation techniques in a purely algebraic manner. This paper analyzes these two options, shows the connections of the volume-agglomeration method with algebraic methods and proposes a new definition of prolongation operator suitable for the application of the volume-agglomeration method to elliptic problems.

Guillard, H.↗
NASA NTRS1981

A multigrid method for the transonic full potential equation discretized with finite elements on an arbitrary body fitted mesh

A multigrid method for the acceleration of transonic potential flow calculations based on a Galerkin finite element approach is described. In order to allow the use of arbitrary body fitted meshes it is necessary to introduce nonuniform interpolation and residual weighting. Emphasis is put on the construction of these operators consistent with the finite element approximation, while standard successive line overrelaxation is used as a smoothing step. Substantial convergence acceleration is obtained and results are presented for different transonic flow configurations including shocks.

Deconinck, H.↗
NASA NTRS1991

A multi-block multigrid method for the solution of the Euler and Navier-Stokes equations for three-dimensional flows

A multi-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell-vertex method which employs central-difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques, and variable-coefficient numerical smoothing results in an efficient and robust scheme. The multi-block strategy places the block loop within the Runge-Kutta loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one- and two-block calculations in which the two-block grid is generated by splitting the one-block grid. Results are presented for both Euler and Navier-Stokes computations of wings and wing-fuselage combinations.

Atkins, H. L.↗
NASA NTRS1996

Numerical Study of Multigrid Methods with Various Smoothers for the Elliptical Grid Generation Equations

A robust solver for the elliptic grid generation equations is sought via a numerical study. The system of PDEs is discretized with finite differences, and multigrid methods are applied to the resulting nonlinear algebraic equations. Multigrid iterations are compared with respect to the robustness and efficiency. Different smoothers are tried to improve the convergence of iterations. The methods are applied to four 2D grid generation problems over a wide range of grid distortions. The results of the study help to select smoothing schemes and the overall multigrid procedures for elliptic grid generation.

Golik, W. L.↗
NASA NTRS2004

On Efficient Multigrid Methods for Materials Processing Flows with Small Particles

Multiscale modeling of materials requires simulations of multiple levels of structural hierarchy. The computational efficiency of numerical methods becomes a critical factor for simulating large physical systems with highly desperate length scales. Multigrid methods are known for their superior efficiency in representing/resolving different levels of physical details. The efficiency is achieved by employing interactively different discretizations on different scales (grids). To assist optimization of manufacturing conditions for materials processing with numerous particles (e.g., dispersion of particles, controlling flow viscosity and clusters), a new multigrid algorithm has been developed for a case of multiscale modeling of flows with small particles that have various length scales. The optimal efficiency of the algorithm is crucial for accurate predictions of the effect of processing conditions (e.g., pressure and velocity gradients) on the local flow fields that control the formation of various microstructures or clusters.

Thomas, James↗
NASA NTRS1993

Numerical solution of flame sheet problems with and without multigrid methods

Flame sheet problems are on the natural route to the numerical solution of multidimensional flames, which, in turn, are important in many engineering applications. In order to model the structure of flames more accurately, we use the vorticity-velocity formulation of the fluid flow equations, as opposed to the streamfunction-vorticity approach. The numerical solution of the resulting nonlinear coupled elliptic partial differential equations involves a pseudo transient process and a steady state Newton iteration. Rather than working with dimensionless variables, we introduce scale factors that can yield significant savings in the execution time. In this context, we also investigate the applicability and performance of several multigrid methods, focusing on nonlinear damped Newton multigrid, using either one way or correction schemes.

Douglas, Craig C.↗
NASA NTRS1982

Numerical boundary condition procedures and multigrid methods; Proceedings of the Symposium, NASA Ames Research Center, Moffett Field, CA, October 19-22, 1981

Papers presented in this volume provide an overview of recent work on numerical boundary condition procedures and multigrid methods. The topics discussed include implicit boundary conditions for the solution of the parabolized Navier-Stokes equations for supersonic flows; far field boundary conditions for compressible flows; and influence of boundary approximations and conditions on finite-difference solutions. Papers are also presented on fully implicit shock tracking and on the stability of two-dimensional hyperbolic initial boundary value problems for explicit and implicit schemes.

Source record↗
NASA NTRS1986

Multigrid method for a vortex breakdown simulation

An inviscid model for a steady axisymmetric flow with swirl was studied. The governing equation is a nonlinear elliptic equation which has more than one solution for a certain range of the swirl parameter. The physically interesting solutions have closed streamlines that look like vortex breakdown (bubble-like solutions). A multigrid method is used to find these solutions. Using an FMG algorithm (nested iteration), the problem is solved in just a few multigrid cycles.

Taasan, S.↗
NASA NTRS1987

An unstructured multigrid method for elliptic problems

A multigrid algorithm for implementation on unstructured meshes is proposed. The algorithm uses a sequence of unnested grids and requires the development of efficient inter-grid interpolation procedures. It is demonstrated how elliptic problems can be solved in this fashion by using Jacobi smoothers.

Lohner, R.↗
NASA NTRS1987

Unstructured multigrid methods

The use of the multigrid procedure with a sequence of unnested coarser grids is discussed. Validity of the procedure is assessed by considering the solution of a single linear elliptic equation. It is demonstrated how a scheme with the optimum order of operations can be constructed. Application to the solution of the Euler equations is considered.

Loehner, R.↗
NASA NTRS1992

Efficient computation of inviscid flow fields around complex configurations using a multiblock multigrid method

The applicability of a multigrid technique to block-structured, body-fitted meshes is examined focusing on three different strategies. In the first strategy data are exchanged between blocks in each stage of a five-stage Runge-Kutta time-stepping scheme which keeps a possible time lag between blocks to a minimum, but requires a large amount of I/O operations and storage. The second strategy is based on performing a complete Runge-Kutta cycle within a block before switching to the next. In the third strategy both a complete Runge-Kutta cycle and the residual evaluation for the restriction operator are done within a block, allowing a minimum of I/O and storage. The inviscid flow around a wing-body/engine-pylon configuration was computed on a mesh consisting of 11 computational blocks. It was found that both the first and the second strategies delivered converged results, but the third failed due to larger time lag between blocks.

Rossow, C.-C.↗
NASA NTRS1989

A cell-vertex multigrid method for the Navier-Stokes equations

A cell-vertex scheme for the Navier-Stokes equations, which is based on central difference approximations and Runge-Kutta time stepping, is described. Using local time stepping, implicit residual smoothing, a multigrid method, and carefully controlled artificial dissipative terms, very good convergence rates are obtained for a wide range of two- and three-dimensional flows over airfoils and wings. The accuracy of the code is examined by grid refinement studies and comparison with experimental data. For an accurate prediction of turbulent flows with strong separations, a modified version of the nonequilibrium turbulence model of Johnson and King is introduced, which is well suited for an implementation into three-dimensional Navier-Stokes codes. It is shown that the solutions for three-dimensional flows with strong separations can be dramatically improved, when a nonequilibrium model of turbulence is used.

Radespiel, R.↗
NASA NTRS1993

High order finite difference and multigrid methods for spatially evolving instability in a planar channel

The fourth-order finite-difference scheme with fully implicit time-marching presently used to computationally study the spatial instability of planar Poiseuille flow incorporates a novel treatment for outflow boundary conditions that renders the buffer area as short as one wavelength. A semicoarsening multigrid method accelerates convergence for the implicit scheme at each time step; a line-distributive relaxation is developed as a robust fast solver that is efficient for anisotropic grids. Computational cost is no greater than that of explicit schemes, and excellent agreement with linear theory is obtained.

Liu, C.↗
NASA NTRS1993

Viscous analysis of three-dimensional rotor flows using a multigrid method

A three-dimensional code for rotating blade-row flow analysis was developed. The space discretization uses a cell-centered scheme with eigenvalues scaling for the artificial dissipation. The computational efficiency of a four-stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full-multigrid method. An application is presented for the NASA rotor 67 transonic fan. Due to the blade stagger and twist, a zonal, non-periodic H-type grid is used to minimize the mesh skewness. The calculation is validated by comparing it with experiments in the range from the maximum flow rate to a near-stall condition. A detailed study of the flow structure near peak efficiency and near stall is presented by means of pressure distribution and particle traces inside boundary layers.

Arnone, A.↗
NASA NTRS1993

Multigrid methods for flow transition in three-dimensional boundary layers with surface roughness

The efficient multilevel adaptive method has been successfully applied to perform direct numerical simulations (DNS) of flow transition in 3-D channels and 3-D boundary layers with 2-D and 3-D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semi-coarsening multigrid method associated with line distributive relaxation scheme, and an improved outflow boundary-condition treatment, which needs only a very short buffer domain to damp all order-one wave reflections, are developed. These approaches make the multigrid DNS code very accurate and efficient. This allows us not only to be able to do spatial DNS for the 3-D channel and flat plate at low computational costs, but also to do spatial DNS for transition in the 3-D boundary layer with 3-D single and multiple roughness elements, which would have extremely high computational costs with conventional methods. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments. The contribution of isolated and distributed roughness to transition is analyzed.

Liu, Chaoqun↗
NASA NTRS1983

A multigrid method for the Euler equations

A multigrid algorithm has been developed for the numerical solution of the steady two-dimensional Euler equations. Flux vector splitting and one-sided differencing are employed to define the spatial discretization. Newton's method is used to solve the nonlinear equations, and a multigrid solver is used on each linear problem. The relaxation scheme for the linear problems is symmetric Gauss-Seidel. Standard restriction and interpolation operators are employed. Local mode analysis is used to predict the convergence rate of the multigrid process on the linear problems. Computed results for transonic flows over airfoils are presented.

Jespersen, D. C.↗
NASA NTRS1983

Preconditioners for the spectral multigrid method

The systems of algebraic equations which arise from spectral discretizations of elliptic equations are full and direct solutions of them are rarely feasible. Iterative methods are an attractive alternative because Fourier transform techniques enable the discrete matrix-vector products to be computed with nearly the same efficiency as is possible for corresponding but sparse finite difference discretizations. For realistic Dirichlet problems preconditioning is essential for acceptable convergence rates. A brief description of Chebyshev spectral approximations and spectral multigrid methods for elliptic problems is given. A survey of preconditioners for Dirichlet problems based on second-order finite difference methods is made. New preconditioning techniques based on higher order finite differences and on the spectral matrix itself are presented. The preconditioners are analyzed in terms of their spectra and numerical examples are presented.

Phillips, T. N.↗