Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELLIPTIC 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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Numerical solution of supersonic three-dimensional free-mixing flows using the parabolic-elliptic Navier-Stokes equations

A numerical method is presented for solving the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to three-dimensional supersonic laminar jet flow issuing parallel with a supersonic free stream. A coordinate transformation is introduced which maps the boundaries at infinity into a finite computational domain in order to eliminate difficulties associated with the imposition of free-stream boundary conditions. Results are presented for an approximate circular jet, a square jet, varying aspect ratio rectangular jets, and interacting square jets. The solution behavior varies from axisymmetric to nearly two-dimensional in character. For cases where comparisons of the present results with those obtained from shear layer calculations could be made, agreement was good.

Hirsh, R. S.↗
NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗
NASA NTRS1975

Calculation of supersonic three-dimensional free-mixing flows using the parabolic-elliptic Navier-Stokes equations

A numerical method is presented which is valid for integration of the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to the three-dimensional supersonic flow of a jet issuing into a supersonic free stream. Difficulties associated with the imposition of free-stream boundary conditions are noted, and a coordinate transformation, which maps the point at infinity onto a finite value, is introduced to alleviate these difficulties. Results are presented for calculations of a square jet and varying-aspect-ratio rectangular jets. The solution behavior varies from axisymmetry for the square jet to nearly two-dimensional for the high-aspect-ratio rectangle, although the computation always calculates the flow as though it were truly three-dimensional.

Hirsh, R. S.↗
NASA NTRS1995

A Study of Two-Equation Turbulence Models on the Elliptic Streamline Flow

Several two-equation turbulence models are compared to data from direct numerical simulations (DNS) of the homogeneous elliptic streamline flow, which combines rotation and strain. The models considered include standard two-equation models and models with corrections for rotational effects. Most of the rotational corrections modify the dissipation rate equation to account for the reduced dissipation rate in rotating turbulent flows, however, the DNS data shows that the production term in the turbulent kinetic energy equation is not modeled correctly by these models. Nonlinear relations for the Reynolds stresses are considered as a means of modifying the production term. Implications for the modeling of turbulent vortices will be discussed.

Blaisdell, Gregory A.↗
NASA NTRS1993

A Reynolds stress model for near-wall turbulence

The paper formulates a tensorially consistent near-wall second-order closure model. Redistributive terms in the Reynolds stress equations are modeled by an elliptic relaxation equation in order to represent strongly nonhomogeneous effects produced by the presence of walls; this replaces the quasi-homogeneous algebraic models that are usually employed, and avoids the need for ad hoc damping functions. The model is solved for channel flow and boundary layers with zero and adverse pressure gradients. Good predictions of Reynolds stress components, mean flow, skin friction, and displacement thickness are obtained in various comparisons to experimental and direct numerical simulation data. The model is also applied to a boundary layer flowing along a wall with a 90-deg, constant-radius, convex bend.

Durbin, P. A.↗
NASA NTRS1984

Three-dimensional viscous flows with large secondary velocity

A new system of approximation equations is derived for three-dimensional steady viscous compressible flows in which a primary-flow direction is present, but in which both transverse velocity components can be large. Previous approaches which address simplification of the steady Navier-Stokes equations are discussed, and a new approach is proposed. The transverse velocity vector which corrects a given potential flow has been decomposed into potential and rotational components. It is found that the potential-velocity vector may be assumed small, whereas the rotational-velocity vector may be assumed small, whereas the rotational velocity vector and hence the composite secondary flow can be of order unity. This assumption leads to a system of governing equations whose characteristic polynomial has a non-elliptic form for arbitrary Mach numbers. The resulting non-elliptic approximation equations can be solved as an initial/boundary-value problem. Computed results confirm the small scalar-potential approximation.

Briley, W. R.↗
NASA NTRS1977

Variants and extensions of a fast direct numerical cauchy-riemann solver, with illustrative applications

Revised and extended versions of a fast, direct (noniterative) numerical Cauchy-Riemann solver are presented for solving finite difference approximations of first order systems of partial differential equations. Although the difference operators treated are linear and elliptic, one significant application of these extended direct Cauchy-Riemann solvers is in the fast, semidirect (iterative) solution of fluid dynamic problems governed by the nonlinear mixed elliptic-hyperbolic equations of transonic flow. Different versions of the algorithms are derived and the corresponding FORTRAN computer programs for a simple example problem are described and listed. The algorithms are demonstrated to be efficient and accurate.

Martin, E. D.↗
NASA NTRS1984

Application of the Bounded Derivative Initialization Method

The utility of the bounded derivative method (BDM) as an initialization tool was examined using real data and the GLAS fourth order global barotropic model. Two constraints on initial data were derived. The first constraint based on the BDM is an equation for calculating initial divergence. The second constraint yields an additional equation of elliptic type for computing the initial height field. This equation is more general than the balance equation first proposed by Charney (1955). It also has more terms which include the contributions of initial divergence and orography. The development of a data insertion scheme based on the BDM is discussed.

Semazzi, F. H. M.↗
NASA NTRS1990

High accuracy solutions of incompressible Navier-Stokes equations

In recent years, high accuracy finite difference approximations were developed for partial differential equations of elliptic type, with particular emphasis on the convection-diffusion equation. These approximations are of compact type, have a local truncation error of fourth order, and allow the use of standard iterative schemes to solve the resulting systems of algebraic equations. These high accuracy approximations are extended to the solution of Navier-Stokes equations. Solutions are obtained for the model problem of driven cavity and are compared with solutions obtained using other approximations and those obtained by other authors. It is discovered that the high order approximations do indeed produce high accuracy solutions and have a potential for use in solving important problems of viscous fluid flows.

Gupta, Murli M.↗
NASA NTRS1986

Solution of separated flows using a simultaneous iteration technique

Generalization and improvement of an earlier work developed for studying separated flows using boundary layer type equations are presented. The improvements include extensions to a general coordinate system and use of a more general zonal technique for solving the coupled equations. In order to be able to consider arbitrary geometries, second order accurate (in space) conservative differences are generated by considering the integral formulation of the governing equations in a general coordinate system. The general coordinate system is handled in as general a manner as possible to allow for the use of either analytically or numerically generated coordinate systems. A marching procedure was used for solving the Partially Parabolized Navier-Stokes (PPNS) equations in the viscous region coupled in a fully implicit manner with the elliptic inviscid equation. To test the algorithm and compare to other solutions, solutions for flow over a flat plate and flow past the symmetrical 12 percent thick Joukowski airfoil (J012) at zero angle of attack were presented.

Source record↗
NASA NTRS1986

Efficient inverse solution of Kepler's equation

A bicubic polynomial approximation to Kepler's equation for elliptic orbits is shown to provide accurate starting values for efficient numerical solution of this equation for eccentric anomaly. In the approximate equation, a cubic in mean anomaly is set equal to a cubic in eccentric anomaly. The coefficients in the two cubics are obtained as functions of eccentricity by specifying values of function and slope at the midpoint and the endpoint of the complete interval (0 to pi). The initial estimate of eccentric anomaly to use in an iteration formula is obtained by evaluating the cubic in mean anomaly and finding the single real root of the cubic in eccentric anomaly. Numerical results are presented which indicate that the estimate accuracy of this method is roughly an order of magnitude better than that of other recently-reported formulas.

Boltz, Frederick W.↗
NASA NTRS2012

Detached-Eddy Simulation Based on the V2-F Model

Detached-eddy simulation (DES) based on the v(sup 2)-f Reynolds-averaged Navier-Stokes (RANS) model is developed and tested. The v(sup 2)-f model incorporates the anisotropy of near-wall turbulence which is absent in other RANS models commonly used in the DES community. The v(sup 2)-f RANS model is modified in order the proposed v(sup 2)-f-based DES formulation reduces to a transport equation for the subgrid-scale kinetic energy isotropic turbulence. First, three coefficients in the elliptic relaxation equation are modified, which is tested in channel flows with friction Reynolds number up to 2000. Then, the proposed v(sup 2)-f DES model formulation is derived. The constant, C(sub DES), required in the DES formulation was calibrated by simulating both decaying and statistically-steady isotropic turbulence. After C(sub DES) was calibrated, the v(sub 2)-f DES formulation is tested for flow around a circular cylinder at a Reynolds number of 3900, in which case turbulence develops after separation. Simulations indicate that this model represents the turbulent wake nearly as accurately as the dynamic Smagorinsky model. Spalart-Allmaras-based DES is also included in the cylinder flow simulation for comparison.

Jee, Sol Keun↗
NASA NTRS1979

Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods

The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.

Periaux, J.↗
NASA NTRS1991

One shot methods for optimal control of distributed parameter systems 1: Finite dimensional control

The efficient numerical treatment of optimal control problems governed by elliptic partial differential equations (PDEs) and systems of elliptic PDEs, where the control is finite dimensional is discussed. Distributed control as well as boundary control cases are discussed. The main characteristic of the new methods is that they are designed to solve the full optimization problem directly, rather than accelerating a descent method by an efficient multigrid solver for the equations involved. The methods use the adjoint state in order to achieve efficient smoother and a robust coarsening strategy. The main idea is the treatment of the control variables on appropriate scales, i.e., control variables that correspond to smooth functions are solved for on coarse grids depending on the smoothness of these functions. Solution of the control problems is achieved with the cost of solving the constraint equations about two to three times (by a multigrid solver). Numerical examples demonstrate the effectiveness of the method proposed in distributed control case, pointwise control and boundary control problems.

Taasan, Shlomo↗
NASA NTRS1978

Viscous flow analysis in mixed flow rotors

A method for analyzing the viscous flow through turbomachine rotors is presented. The field analysis is based upon the solution of the full Navier-Stokes equations over the rotor blade-to-blade stream channels. An Alternating-Direction-Implicit method is employed to carry out the necessary numerical integration of the elliptic governing equations. The flow analysis may be applied to various types of turbomachine rotors. Preliminarily, only the case of laminar flows are considered in this paper. The flow characteristics within the rotors of a mixed flow turbine and a radial bladed compressor are investigated over a wide range of operating conditions. Excellent results are obtained when compared with existing experimental data. The method of this analysis is quite general and can deal with a wide range of applications. Possible modification of the present study to deal with turbulent flow cases are also identified.

Khalil, I.↗