A three-dimensional incompressible Navier-Stokes flow solver using primitive variables
Explore the source record for details and available documents.
Kwak, D.↗
SEARCH · Engineering Papers
Results for “solver”
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.
On the structure of parallelism in a highly concurrent PDE solver
A parallel multigrid algorithm for solving elliptic partial differential equations is developed and evaluated. A V-cycle multigrid method is altered to increase the degree of parallelism. A numerical analysis of the resulting concurrent-iteration multigrid algorithm is performed; its architectural implications are considered; highly parallel systems without shared memory are examined (including mesh-connected arrays, mesh-shuffle-connected systems, permutation networks, and direct VLSI embeddings); and the results of numerical experiments are presented in tables and graphs.
Gannon, D.↗
Development of an iterative boundary-layer-type solver for axisymmetric separated flows
Explore the source record for details and available documents.
Halim, A. A. M.↗
Solution of elliptic partial differential equations by fast Poisson solvers using a local relaxation factor. 2: Two-step method
A two-step semidirect procedure is developed to accelerate the one-step procedure described in NASA TP-2529. For a set of constant coefficient model problems, the acceleration factor increases from 1 to 2 as the one-step procedure convergence rate decreases from + infinity to 0. It is also shown numerically that the two-step procedure can substantially accelerate the convergence of the numerical solution of many partial differential equations (PDE's) with variable coefficients.
Chang, S. C.↗
Analysis of a parallelized nonlinear elliptic boundary value problem solver with application to reacting flows
A parallelized finite difference code based on the Newton method for systems of nonlinear elliptic boundary value problems in two dimensions is analyzed in terms of computational complexity and parallel efficiency. An approximate cost function depending on 15 dimensionless parameters is derived for algorithms based on stripwise and boxwise decompositions of the domain and a one-to-one assignment of the strip or box subdomains to processors. The sensitivity of the cost functions to the parameters is explored in regions of parameter space corresponding to model small-order systems with inexpensive function evaluations and also a coupled system of nineteen equations with very expensive function evaluations. The algorithm was implemented on the Intel Hypercube, and some experimental results for the model problems with stripwise decompositions are presented and compared with the theory. In the context of computational combustion problems, multiprocessors of either message-passing or shared-memory type may be employed with stripwise decompositions to realize speedup of O(n), where n is mesh resolution in one direction, for reasonable n.
Keyes, David E.↗
Finite element solver for 3-D compressible viscous flows
The space shuttle main engine (SSME) has extremely complex internal flow structure. The geometry of the flow domain is three-dimensional with complicated topology. The flow is compressible, viscous, and turbulent with large gradients in flow quantities and regions of recirculations. The analysis of the flow field in SSME involves several tedious steps. One is the geometrical modeling of the particular zone of the SSME being studied. Accessing the geometry definition, digitalizing it, and developing surface interpolations suitable for an interior grid generator require considerable amount of manual labor. There are several types of grid generators available with some general-purpose finite element programs. An efficient and robust computational scheme for solving 3D Navier-Stokes equations has to be implemented. Post processing software has to be adapted to visualize and analyze the computed 3D flow field. The progress made in a project to develop software for the analysis of the flow is discussed. The technical approach to the development of the finite element scheme and the relaxation procedure are discussed. The three dimensional finite element code for the compressible Navier-Stokes equations is listed.
Reddy, K. C.↗
An exterior Poisson solver using fast direct methods and boundary integral equations with applications to nonlinear potential flow
A general method is developed combining fast direct methods and boundary integral equation methods to solve Poisson's equation on irregular exterior regions. The method requires O(N log N) operations where N is the number of grid points. Error estimates are given that hold for regions with corners and other boundary irregularities. Computational results are given in the context of computational aerodynamics for a two-dimensional lifting airfoil. Solutions of boundary integral equations for lifting and nonlifting aerodynamic configurations using preconditioned conjugate gradient are examined for varying degrees of thinness.
Young, D. P.↗
Solution of elliptic PDEs by fast Poisson solvers using a local relaxation factor
A large class of two- and three-dimensional, nonseparable elliptic partial differential equations (PDEs) is presently solved by means of novel one-step (D'Yakanov-Gunn) and two-step (accelerated one-step) iterative procedures, using a local, discrete Fourier analysis. In addition to being easily implemented and applicable to a variety of boundary conditions, these procedures are found to be computationally efficient on the basis of the results of numerical comparison with other established methods, which lack the present one's: (1) insensitivity to grid cell size and aspect ratio, and (2) ease of convergence rate estimation by means of the coefficient of the PDE being solved. The two-step procedure is numerically demonstrated to outperform the one-step procedure in the case of PDEs with variable coefficients.
Chang, Sin-Chung↗
Euler solver for three-dimensional supersonic flows with subsonic pockets
A new finite-difference scheme has been developed to solve efficiently the unsteady Euler equations for three-dimensional inviscid supersonic flows with subsonic pockets. The technique utilizes planar Gauss-Seidel relaxation in the marching direction and approximate factorization n the crossflow plane. An 'infinitely large' time step is used in parts of the flowfield where the component of velocity in the marching direction is supersonic - here the Gauss-Seidel sweeps are restricted to the forward direction only, and the procedure reduces to simple space-marching; a finite time step is used in parts of the flowfield where the marching component of velocity is subsonic - here, backward and forward Gauss-Seidel sweeps are employed to allow for upstream and downstream propagation of signals, and a time-asymptotic steady state is obtained. The discretization formulas are based on finite-volume implementation of high accuracy (up to third-order) total variation diminishing formulations. Numerical solutions are obtained for an analytically defined forebody, a realistic fighter configuration, and the Space Shuttle. The results are in very good agreement with available experimental data and numerical solutions of the full-potential equation.
Chakravarthy, Sukumar R.↗
Iterative finite element solver on transputer networks
The parallelism inherent in the Conjugate Gradient method is described. The initial results of a parallel implementation on a network of twelve transputers are discussed. The high efficiencies obtained indicate that significant speedup can be obtained with larger transputer arrays if communication overhead can be kept low. To this end, a method of communication that allows large, dynamically reconfigurable transputer arrays to exchange data in log sub 4 N steps for N processors is suggested.
Danial, Albert↗
Steady and transient least square solvers for thermal problems
This paper develops a hierarchical least square solution algorithm for highly nonlinear heat transfer problems. The methodology's capability is such that both steady and transient implicit formulations can be handled. This includes problems arising from highly nonlinear heat transfer systems modeled by either finite-element or finite-difference schemes. The overall procedure developed enables localized updating, iteration, and convergence checking as well as constraint application. The localized updating can be performed at a variety of hierarchical levels, i.e., degree of freedom, substructural, material-nonlinear groups, and/or boundary groups. The choice of such partitions can be made via energy partitioning or nonlinearity levels as well as by user selection. Overall, this leads to extremely robust computational characteristics. To demonstrate the methodology, problems are drawn from nonlinear heat conduction. These are used to quantify the robust capabilities of the hierarchical least square scheme.
Padovan, Joe↗
A parallelized elliptic solver for reacting flows
A modified Newton algorithm for the solution of nonlinear elliptic boundary value problems via finite discretization methods is presented. A serial implementation of this algorithm which has recently been applied successfully to the computation of an axisymmetric over-ventilated subsonic laminar methane-air jet diffusion flame is described. Parallel implementation issues and a complexity theory are presented. Included as well are actual performance data for model systems obtained on the Intel Hypercube and a discussion of its implications for modeling realistic systems.
Keyes, David E.↗
Third-order multi-dimensional Euler/Navier-Stokes solver
A nonoscillatory multidimensional advection code is presented which is applicable to the study of very high Reynolds number steady flow, and results are obtained for the oblique advection of a step profile. The scheme is based on an exponential upwinding or linear extrapolation refinement of a previous QUICK scheme, allowing sharp monotonic resolution of the step while retaining third-order accuracy. The resulting simple flux-limiting criteria have been applied to a higher order method and a formally third-order supercompressive algorithm. Results are presented for first-order, second-order, and third-order upwinding.
Leonard, B. P.↗
A transputer based finite element solver
The feasibility of performing FEM structural-mechanics analyses on transputer systems is investigated experimentally. Transputers are programmable microprocessors equipped with local memory and point-to-point communication links; they can be joined in a large concurrent system via a programming language which supports distributed processing; this permits parallel processing at relatively low hardware cost. The computational tasks required by FEM programs are reviewed; the hardware (one PC, one master transputer, and 12 slave transputers) employed in the test calculations is described; and results demonstrating the speed and efficiency of the transputer array in assembling a global stiffness matrix and performing Gauss-Jordan matrix inversion are presented in graphs. It is predicted that larger transputer networks could approach the power of supercomputers at minicomputer costs.
Favenesi, J. A.↗
Control-volume based Navier-Stokes equation solver valid at all flow velocities
A control-volume based finite difference method to solve the Reynolds averaged Navier-Stokes equations is presented. A pressure correction equation valid at all flow velocities and a pressure staggered grid layout are used in the method. Example problems presented herein include: a developing laminar channel flow, developing laminar pipe flow, a lid-driven square cavity flow, a laminar flow through a 90-degree bent channel, a laminar polar cavity flow, and a turbulent supersonic flow over a compression ramp. A k-epsilon turbulence model supplemented with a near-wall turbulence model was used to solve the turbulent flow. It is shown that the method yields accurate computational results even when highly skewed, unequally spaced, curved grids are used. It is also shown that the method is strongly convergent for high Reynolds number flows.
Kim, S.-W.↗
Parallel linear equation solvers for finite element computations
The overall objective of this research is to develop efficient methods for the solution of linear and nonlinear systems of equations on parallel and supercomputers, and to apply these methods to the solution of problems in structural analysis. Attention has been given so far only to linear equations. The methods considered for the solution of the stiffness equation Kx=f have been Choleski factorization and the conjugate gradient iteration with SSOR and Incomplete Choleski preconditioning. More detail on these methods will be given on subsequent slides. These methods have been used to solve for the static displacements for the mast and panel focus problems in conjunction with the CSM testbed system based on NICE/SPAR.
Ortega, James M.↗
Aerodynamic applications of an efficient incompressible Navier-Stokes solver
A robust and accurate upwind scheme based on flux-difference splitting for three-dimensional, incompressible flows in general coordinates is described. Viscous effects are included through a thin-layer approximation to the Navier-Stokes equations. Turbulent flows are simulated with extensions to the Baldwin-Lomax (1978) turbulence model. The effects of accuracy and gridding on the computed flow field results are assessed. The results obtained from the code for vortical flows over a delta wing, a double delta wing, and a tangent-ogive forebody compare well with experimental data.
Hartwich, Peter M.↗
Application of direct solvers to unstructured meshes for the Euler and Navier-Stokes equations using upwind schemes
The application of Newton iteration to inviscid and viscous airfoil calculations on unstructured meshes is examined. A cell-centered finite volume scheme is employed on an unstructured mesh consisting of triangles. Roe's flux difference splitting scheme is used to compute the inviscid fluxes. Higher order accuracy is achieved by an interpolation procedure that makes use of auxiliary gradients. The efficient solution of the sparse linear system of equations which arises upon linearization in time is addressed. Results are presented for inviscid and viscous test cases. The complications which arise due to the introduction of nonlinear limiters are addressed.
Venkatakrishnan, V.↗