SEARCH · Engineering Papers

Results for “Poisson equation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1991

A fractional step solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

The time-dependent, three-dimensional incompressible Navier-Stokes equations are presently solved in generalized coordinate systems by means of a fractional-step method whose primitive variable formulation uses as dependent variables, in place of the Cartesian components of the velocity: (1) pressure (defined at the center of the computational cell), and (2) volume fluxes across the faces of the cells. The momentum equations are solved by means of an approximate factorization method. A novel 'ZEBRA' scheme incorporating four-color ordering efficiently solves the Poisson equation. Illustrative two- and three-dimensional laminar flow test cases are computed and evaluated relative to extant numerical and experimental results, and good agreement is obtained.

Rosenfeld, Moshe↗

NASA NTRS1981

Comparison of theory and in situ observations for electron and ion distributions in the near wake of the Explorer 31 and AE-C satellites

Measurements of electron density, plasma potential, and mean ion mass from the Explorer 31 satellite, and measurements of ion current, plasma potential, and ion composition from the Atmosphere Explorer C satellite were used in a comparative study with Parker's theory regarding the charged particle distribution in the near wake of an ionospheric satellite (1976). It is shown that theory and experiment agree fairly well in the angle-of-attack range between 90 and 135 deg. In the maximum rarefaction zone (between 145 and 180 deg), however, the theoretical model overestimates the measured ion depletion by several orders of magnitude. A comparison between theory and the Explorer 31 electron measurements shows that the theory again overestimates the electron depletion. These discrepancies are mainly due to the use of a steady-state theory and a single ion equation (using a mean ion mass). Improved agreement between theory and experiment can be obtained by the use of the time-dependent Vlasov-Poisson equations with separate equations for the various ion species.

Samir, U.↗

DOE OSTI2025

Physics-based stabilized finite element approximations of the Poisson–Nernst–Planck equations

We present and analyze two stabilized finite element methods for solving numerically the Poisson–Nernst–Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for the ion equations, whereas the discrete equation for the electric potential need not be stabilized. Discrete solutions stemmed from the first algorithm preserve both maximum and minimum discrete principles. For the second algorithm, its discrete solutions are conceived so that they hold discrete principles and obey an entropy law provided that an acuteness condition is imposed for meshes. Remarkably the latter is found to be unconditionally stable. We validate our methodology through transient numerical experiments that show convergence toward steady-state solutions.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Prediction of cascade performance using an incompressible Navier-Stokes technique

A fully elliptic, control volume solution of the two-dimensional incompressible Navier-Stokes equations for the prediction of cascade performance over a wide incidence range is presented. The numerical technique is based on a new pressure substitution method. A Poisson equation is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. The analysis of a double circular arc compressor cascade is presented, and the results are compared with the available experimental data at various incidence angles. Good agreement is obtained for the blade pressure distribution, boundary layer and wake profiles, skin friction coefficient, losses and outlet angles. Turbulence effects are simulated by the Low-Reynolds-Number version of the k-epsilon turbulence model.

Hobson, G. V.↗

NASA NTRS1990

Simulation Of Unsteady, Viscous, Incompressible Flow

Method for numerical solution of Navier-Stokes equations of viscous, incompressible flow developed based on use of fractional-step procedure. Accurate to second order in both space and time. Attempt made to minimize Poisson-equation difficulties by choosing pressures at centers and volume fluxes across faces of computational cells as dependent variables instead of familiar Cartesian components of velocity. Choice ensures satisfaction of discrete equation of conservation of mass to within round-off errors in any coordinate system and has favorable effects on convergence properties.

Rosenfeld, Moshe↗

NASA NTRS1997

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

NASA NTRS2008

Flow Solver for Incompressible 2-D Drive Cavity

This software solves the Navier-Stokes equations for the incompressible driven cavity flow problem. The code uses second-order finite differencing on a staggered grid using the Chorin projection method. The resulting intermediate Poisson equation is efficiently solved using the fast Fourier transform. Time stepping is done using fourth-order Runge-Kutta for stability at high Reynolds numbers. Features include check-pointing, periodic field snapshots, ongoing reporting of kinetic energy and changes between time steps, time histories at selected points, and optional streakline generation.

Kalb, Virginia L.↗

NASA NTRS1986

A coupled marching procedure for the partially parabolized Navier-Stokes equations

A coupled finite-difference formulation is described for solving a reduced form of the compressible Navier-Stokes equations by a multiple space marching procedure. The properties of the equations are discussed and the solution algorithm presented. The scheme is used to compute incompressible flows by taking the incompressible limit (M about 0.1) of the compressible formulation. A separate Poisson equation is not required for the pressure. Results are compared with experimental data and other numerical predictions for low Reynolds number channel inlet flow, flow over a rearward-facing step in a channel, and flow near the trailing edge of a flat plate.

Liu, Xuezong↗

NASA NTRS1985

A space-marching method for incompressible Navier-Stokes equations

This paper deals with the development of a space-marching method for incompressible flows. The method solves the continuity and momentum equations as a coupled system at each streamwise station. The character of the system of equations has been changed from elliptic to hyperbolic/parabolic in order to enable the equations to be marched in space. The present method has many advantages compared to the existing parabolic or space-marching methods for incompressible flow: (1) it avoids the solution of Poisson equations, (2) it conserves the mass flow with no additional computation, (3) it does not require the specification of an assumed pressure field when used in the prediction of duct flows. The present method can capture strong secondary velocities and strong transverse pressure gradients. Predictions of the flow through straight and curved ducts are in good agreement with analytical and experimental results.

Pouagare, M.↗

NASA NTRS1979

Numerical prediction of three-dimensional juncture region flow using the parabolic Navier-Stokes equations

A numerical solution algorithm is established for prediction of subsonic turbulent three-dimensional flows in aerodynamic configuration juncture regions. A turbulence closure model is established using the complete Reynolds stress. Pressure coupling is accomplished using the concepts of complementary and particular solutions to a Poisson equation. Specifications for data input juncture geometry modification are presented.

Baker, A. J.↗

NASA NTRS1999

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

NASA NTRS1990

Computation of turbine flowfields with a Navier-Stokes code

A new technique has been developed for the solution of the incompressible Navier-Stokes equations. The numerical technique, derived from a pressure substitution method (PSM), overcomes many of the deficiencies of the pressure crrection method. This technique allows for the direct solution of the actual pressure in the form of a Poisson equation which is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. In two-dimensions a turbine flowfield, including heat transfer, has been computed with this method and the prediction of the cascade performance is presented. The extension of the pressure correction method for the solution of three-dimensional flows is also presented for laminar flow in an S-shaped duct and turbulent flow in the end-wall region of a turbine cascade.

Hobson, G. V.↗

NASA NTRS1991

A pressure based method for the solution of viscous incompressible turbomachinery flows

A new technique was developed for the solution of the incompressible Navier-Stokes equations. The numerical technique, derived from a pressure substitution method (PSM), overcomes many of the deficiencies of the pressure correction method. This technique allows for the direct solution of the actual pressure in the form of a Poisson equation which is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. Two dimensional internal flows are computed with this method. The prediction of cascade performance is presented. The extention of the pressure correction method for the solution of three dimensional flows is also presented.

Hobson, Garth Victor↗

NASA NTRS1987

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.

Erlebacher, G.↗

NASA NTRS1988

Spectral multigrid methods for the solution of homogeneous turbulence problems

Erlebacher, G.↗

DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Perturbation solution of the Navier-Stokes equations and its relation to the Lighthill-Curle solution of aerodynamic sound

The aerodynamic sound described by the Lighthill-Curle solution is reexamined using a method of matched asymptotic expansions. The governing Navier-Stokes equations written in nondimensional form are expanded for a small Mach number. First- and second-order solutions for the pressure field are obtained, and the singular nature of the expansion at large distances is indicated. The nearfield pressure is governed by the Poisson equation, whereas the farfield equations describe a linear wave system in a dissipative medium. The pseudosound is related to the incompressible Reynolds stresses associated with a solenoidal velocity field, the velocity, the pressure perturbation, and their derivatives on the boundaries. A uniformly valid first-order solution for the pressure is obtained. It is shown that viscosity, thermal conductivity, and entropy in the flow do not contribute to the first-order noise generation, while the viscous stress contributes to noise only from some boundaries. The application of the proposed perturbation method to a subsonically moving surface and a hot jet is discussed.

Pan, Y. S.↗

NASA NTRS1983

Numerical simulation of the leading-edge separation vortex for a wing and strake-wing configuration

In the present investigation, the 'thin layer' Navier-Stokes equations are used to compute the flow about a delta wing and a strake-delta wing configuration. Both configurations possess blunt noses and rounded leading edges. The computational grid about these configurations is generated using a newly developed three-dimensional grid-generation code which solves a set of Poisson equations with spherical coordinate variables using an alternating direction implicit (ADI) scheme. Computational results are obtained for an isolated wing with 60 deg sweep and a strake-wing configuration with 80-60 deg sweep. Computations for angles of attack in the range from 6 to 30 deg are presented. Attention is given to the effect of the angle of attack, the Mach number, and the Reynolds number. Theoretical results are compared with experimental data.

Fujii, K.↗