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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft

This report compiles the various research activities conducted under the auspices of the NASA Grant NAG3-1026, "Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft" during the period of April 1989 to April 1994. The effort involved the development of multigrid based algorithms and computer programs for the calculation of the flow and temperature fields generated by Short Take-off and Vertical Landing (STOVL) aircraft, while hovering in ground proximity. Of particular importance has been the interaction of the exhaust jets with the head wind which gives rise to the hot gas ingestion process. The objective of new STOVL designs to reduce the temperature of the gases ingested into the engine. The present work describes a solution algorithm for the multi-dimensional elliptic partial-differential equations governing fluid flow and heat transfer in general curvilinear coordinates. The solution algorithm is based on the multigrid technique which obtains rapid convergence of the iterative numerical procedure for the discrete equations. Initial efforts were concerned with the solution of the Cartesian form of the equations. This algorithm was applied to a simulated STOVL configuration in rectangular coordinates. In the next phase of the work, a computer code for general curvilinear coordinates was constructed. This was applied to model STOVL geometries on curvilinear grids. The code was also validated in model problems. In all these efforts, the standard k-Epsilon model was used.

Vanka, S. P.↗

NASA NTRS1976

Steady supersonic flowfields with embedded subsonic regions

Supersonic flow past a blunt body is considered, where the flow contains an embedded subsonic region which lies between the shock wave and the body surface and is bounded by sonic lines from the body to the shock. A numerical approach is taken, which uses a basic finite difference scheme that solves the unsteady fluid dynamic equations in integral form. The unsteady equations are everywhere hyperbolic in time so no distinction need be made between subsonic and supersonic regions. Solutions to the mixed elliptic and hyperbolic steady flow equations are approached asymptotically in time. The method is illustrated for two-dimensional flows.

Maccormack, R. W.↗

NASA NTRS1961

An Investigation of the Natural Frequencies and Mode Shapes of Liquids in Oblate Spheroidal Tanks

An experimental investigation was conducted to gain some understanding of the character of the free vibration modes of liquids In oblate spheroidal tanks applicable in missile and space vehicle systems, Measured natural frequencies were obtained for the lowest three or four, antisymmetric modes of oscillation as a function of the liquid depth for three orientations of each of several such tanks of different size and oblateness. The orientations considered were such that: (a) the equator of the spheroid was horizontal and oscillations were along a diameter of the circular liquid surface; (b) the equator of the spheroid was vertical and the oscillations were along the minor axis of the elliptical liquid surface; and (c) the equator of the spheroid was vertical and the oscillations were along the major axis of the elliptical liquid surface; The frequency data are presented as dimensionless parameters developed for each orientation to permit the application of the experimental results to the prediction of the natural frequencies of tanks of different size and oblateness. Photographs we re made of representative surface wave or mode, shapes for each orientation.

Leonard, H. Wayne↗

NASA NTRS1978

Minimum film thickness in elliptical contacts for different regimes of fluid-film lubrication

The film-parameter equations are provided for four fluid-film lubrication regimes found in elliptical contacts. These regimes are isoviscous-rigid; viscous-rigid; elastohydrodynamic of low-elastic-modulus materials, or isoviscous-elastic; and elastohydrodynamic, or viscous-elastic. The influence or lack of influence of elastic and viscous effects is the factor that distinguishes these regimes. The film-parameter equations for the respective regimes come from earlier theoretical studies by the authors on elastohydrodynamic and hydrodynamic lubrication of elliptical conjunctions. These equations are restated and the results are presented as a map of the lubrication regimes, with film-thickness contours on a log-log grid of the viscosity and elasticity parameters for five values of the ellipticity parameter. The results present a complete theoretical film-parameter solution for elliptical contacts in the four lubrication regimes.

Hamrock, B. J.↗

NASA NTRS1999

A Pseubo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

Morrison, J. H.↗

NASA NTRS1999

A Pseudo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

White, J. A.↗

NASA NTRS1975

A fast semidirect method for computing transonic aerodynamic flows

A fast, semidirect, iterative computational method, previously introduced for finite-difference solution of subsonic and slightly supercritical flow over airfoils, is extended both to apply to strongly supercritical conditions and to include full second-order accuracy in computing inviscid flows over airfoils. The nonlinear small-disturbance equations are solved iteratively by a direct, linear, elliptic solver. General, fully conservative, type-dependent difference equations are formulated, including parabolic- and shock-point transition operators that provide consistency with the integral conservation laws. These equations specialize to either first-order or to fully second-order-accurate equations. Various free parameters are evaluated for rapid convergence of the first-order scheme. Resulting pressure distributions and computing times are compared with the improved Murman-Cole line-relaxation method.

Martin, E. D.↗

NASA NTRS1977

On two-dimensional packets of capillary-gravity waves

The motion of a two-dimensional packet of capillary-gravity waves on water of finite depth is studied. The evolution of a packet is described by two partial differential equations: the nonlinear Schroedinger equation with a forcing term and a linear equation, which is of either elliptic or hyperbolic type depending on whether the group velocity of the capillary-gravity wave is less than or greater than the velocity of long gravity waves. These equations are used to examine the stability of the Stokes capillary-gravity wave train. The analysis reveals the existence of a resonant interaction between a capillary-gravity wave and a long gravity wave. The interaction requires that the liquid depth be small in comparison with the wavelength of the (long) gravity waves and the evolution equations describing the dynamics of this interaction are derived.

Djordjevic, V. D.↗

NASA NTRS1976

Simplified solution for point contact deformation between two elastic solids

A linear-regression by the method of least squares is made on the geometric variables that occur in the equation for point contact deformation. The ellipticity and the complete eliptic integrals of the first and second kind are expressed as a function of the x, y-plane principal radii. The ellipticity was varied from 1 (circular contact) to 10 (a configuration approaching line contact). These simplified equations enable one to calculate easily the point-contact deformation to within 3 percent without resorting to charts or numerical methods.

Brewe, D. E.↗

NASA NTRS1982

Boundary-fitted coordinate systems for numerical solution of partial differential equations - A review

A comprehensive review of methods of numerically generating curvilinear coordinate systems with coordinate lines coincident with all boundary segments is given. Some general mathematical framework and error analysis common to such coordinate systems is also included. The general categories of generating systems are those based on conformal mapping, orthogonal systems, nearly orthogonal systems, systems produced as the solution of elliptic and hyperbolic partial differential equations, and systems generated algebraically by interpolation among the boundaries. Also covered are the control of coordinate line spacing by functions embedded in the partial differential operators of the generating system and by subsequent stretching transformation. Dynamically adaptive coordinate systems, coupled with the physical solution, and time-dependent systems that follow moving boundaries are treated. References reporting experience using such coordinate systems are reviewed as well as those covering the system development.

Thompson, J. F.↗

NASA NTRS1983

Adaptive grids generated by elliptic systems

It is pointed out that a finite difference grid which moves with the solution of the partial differential equation being solved can improve the accuracy and efficiency of a numerical algorithm. This technique is particularly advantageous in the solution of problems involving boundary layers or shocks where a poorly chosen grid may give a numerical solution which is useless because of poor resolution or extreme oscillations. The present investigation is concerned with the development of a scheme which does not excessively distort the grid. The grid generation algorithm is based on the numerical solution of a system of elliptic differential equations. Holst and Brown (1981) have used a preliminary solution to move points on the boundary of the physical region and then resolved the problem on a new grid generated by an elliptic system. In the current investigation, the grid movement and the solution will develop simultaneously. The solution is used to modify the generating equations, thereby controlling the grid point distributions.

Mastin, C. W.↗

NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗

NASA NTRS1991

Faster Algorithm For Computation Of Incompressible Flow

Improved algorithm yields faster numerical solutions of Navier-Stokes equations of steady or unsteady three-dimensional flow of incompressible fluid. In artificial-compressibility method, unsteady flow treated as incompressible in advancing from one time step to next, but at each time step (or in steady state), fluid treated as having variable compressibility enabling propagation of flow field, and subiterations performed in increments of pseudotime until effects of compressibility subside. Directly couples pressure and velocity fields at same time step and converts elliptic incompressible Navier-Stokes equations to hyperbolic form more amenable to numerical integration.

Rogers, S. E.↗

NASA NTRS2002

Nonlinear Viscoelastic Mechanics of Cross-linked Rubbers

The paper develops a general theory for finite rubber viscoelasticity, and specifies it in the form, convenient for solving problems important for rubber, tire and space industries. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory has been developed for arbitrary nonisothermal deformations of viscoelastic solids. In this theory, the constitutive equations are presented as the sum of known equilibrium (rubber elastic) and non-equilibrium (liquid polymer viscoelastic) terms. These equations are then simplified using several modeling arguments. Stability constraints for the proposed constitutive equations are also discussed. It is shown that only strong ellipticity criteria are applicable for assessing stability of the equations governing viscoelastic solids.

Freed, Alan D.↗

NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

NASA NTRS1989

Efficient grid generation

Because the governing equations in fluid dynamics contain partial differentials and are too difficult in most cases to solve analytically, these differentials are generally replaced by finite difference terms. These terms contain terms in the solution at nearby states. This procedure discretizes the field into a finite number of states. These states, when plotted, form a grid, or mesh, of points. It is at these states, or field points, that the solution is found. The optimum choice of states, the x, y, z coordinate values, minimizes error and computational time. But the process of finding these states is made more difficult by complex boundaries, and by the need to control step size differences between the states, that is, the need to control the spacing of field points. One solution technique uses a different set of state variables, which define a different coordinate system, to generate the grid more easily. A new method, developed by Dr. Joseph Steger, combines elliptic and hyperbolic partial differential equations into a mapping function between the physical and computational coordinate systems. This system of equations offers more control than either equation provides alone. The Steger algorithm was modified in order to allow bodies with stronger concavities to be used, offering the possibility of generating a single grid about multiple bodies. Work was also done on identifying areas where grid breakdown occurs.

Seki, Rycichi↗

NASA NTRS1986

Numerical simulation of vortical flow over an elliptical-body missile at high angles of attack

Numerical solutions to the Reynolds-averaged Navier-Stokes equations are given for the flow about an elliptical body missile (3:1 ellipse) at a Mach number of 2.5 and a unit Reynolds number of 6.56 x 10 to the 6th/m. At high angles of attack, the flow is dominated by large-scale free vortices which occur in the lee-side flow field due to crossflow boundary-layer separation. Emphasis is focused on the accurate prediction of the lee-side vortical flow. Solutions are presented for both symmetric and asymmetric (body rolled 45 deg) configurations at 10 deg and 20 deg angle of attack. The computed results are compared with experimental surface pressure coefficients and vapor-screen photographs. Excellent agreement is obtained in all cases.

Newsome, R. N.↗

NASA NTRS1975

Analysis of atmospheric flow over a surface protrusion using the turbulence kinetic energy equation

Atmospheric flow fields resulting from a semi-elliptical surface obstruction in an otherwise horizontally homogeneous statistically stationary flow are modelled with the boundary-layer/Boussinesq-approximation of the governing equation of fluid mechanics. The turbulence kinetic energy equation is used to determine the dissipative effects of turbulent shear on the mean flow. Mean-flow results are compared with those given in a previous paper where the same problem was attacked using a Prandtl mixing-length hypothesis. Iso-lines of turbulence kinetic energy and turbulence intensity are plotted in the plane of the flow. They highlight regions of high turbulence intensity in the stagnation zone and sharp gradients in intensity along the transition from adverse to favourable pressure gradient.

Frost, W.↗