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 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1978

New formulation of de Sitter's theory of motion for Jupiter I-IV. I - Equations of motion and the disturbing function

Elliptic orbits are substituted for circular orbits in the first approximation, in an analysis of the common retrograde motion of Jupiter's satellites. A modification of the de Sitter theory, made possible by extended observations of the satellites, is presented with attention to that aspect of the theory which eliminates small divisors at all stages of the solution. The convergence problem is circumvented by use of Poincare's canonical relative coordinates. In addition, modified Delaunay variables and their associated Poincare variables are applied to the disturbing function, which is expanded by means of generalized Newcomb operators.

Aksnes, K.↗

NASA NTRS1978

Additional aspects of elastohydrodynamic lubrication

Elastohydrodynamic lubrication (EHL) for materials of low elastic modulus is considered. Engineering applications in which EHL is important for low-elastic-modulus materials include seals, human joints, tires, and elastomeric-material machine elements. Theoretical solutions of the problem of fully flooded line contacts in the case of low-elastic-materials are discussed. The equation regarding dimensionless minimum film thickness for fully flooded elliptical contacts for low elastic modulus is compared with the corresponding equation for materials of high elastic modulus. The powers of the dimensionless speed parameter U are quite similar in both equations, but the power of the dimensionless load parameter W is much more significant for low-elastic-modulus materials. Attention is given to aspects of dimensionless grouping, the isoviscous rigid regime, the viscous rigid regime, the isoviscous elastic regime, and the viscous elastic regime.

Hamrock, B. J.↗

NASA NTRS1982

The gravito-electrodynamics of charged dust in planetary magnetospheres

The orbital dynamics of small electrically charged dust grains within the corotating regions of planetary magnetospheres is considered. Equations are derived for the elliptical epicyclic motion of positively and negatively charged particles about the guiding center in an equilibrium circular orbit under the influence of small perturbations. Those orbits that are stable to the perturbations are found to have a ratio of the semiaxes of the epicycle between 1/2 and 1, depending on the specific charge, and a gyration frequency about the guiding center between the Kepler frequency and the grain gyrofrequency in a nonrotating frame. In the magnetospheres of Jupiter and Saturn, where the grains are expected to be negatively charged and move in a prograde sense, the guiding centers are predicted to have speeds intermediate to the Kepler speed and the corotation speed and thus may undergo a 1:1 magneto-gravitational resonance with a neighboring satellite. Results may be used in the interpretation of the waves in the F ring of Saturn in terms of the dust size distribution.

Mendis, D. A.↗

NASA NTRS1982

Aerothermal modeling

The objectives, approach, and status of a program to develop the computational fluid dynamics tools needed to improve combustor design and analysis are outlined. The calculation procedure selected consists of a finite difference solution of the time averaged, steady state, primitive variable, elliptic form of the Reynolds equations. Standard TEACH type numerics are used to solve the resulting equations. These include hybrid differencing, SIMPLE algorithm for the pressure field, line by line iterative solution using the ADI method and the tridiagonal matrix algorithm (TDMA). Convergence is facilitated by using under relaxation. The physical processes are modeled by a two equation eddy viscosity model for turbulence; combustion is represented by a simple, irreversible, one step chemical reaction whose rate is influenced only by the time scale of the turbulence. The model evaluation procedure is also described.

Sturgess, G. J.↗

NASA NTRS1987

Algebraic multigrid

The state of the art in algebraic multgrid (AMG) methods is discussed. The interaction between the relaxation process and the coarse grid correction necessary for proper behavior of the solution probes is discussed in detail. Sufficient conditions on relaxation and interpolation for the convergence of the V-cycle are given. The relaxation used in AMG, what smoothing means in an algebraic setting, and how it relates to the existing theory are considered. Some properties of the coarse grid operator are discussed, and results on the convergence of two-level and multilevel convergence are given. Details of an algorithm particularly studied for problems obtained by discretizing a single elliptic, second order partial differential equation are given. Results of experiments with such problems using both finite difference and finite element discretizations are presented.

Ruge, J. W.↗

NASA NTRS1989

A numerical analysis of the effects of conjugate heat transfer, vapor compressibility, and viscous dissipation in heat pipes

The effects of conjugate heat transfer, vapor compressibility, and viscous dissipation in heat pipes are discussed. The accuracy of the partially parabolic versus the elliptic presentation of the governing equations is also examined. The results show that the axial wall conduction has a tendency to make the temperature distribution more uniform for heat pipes with large ratios of pipe wall to effective liquid-wick thermal conductivity. The compressible and incompressible models show very close agreement for the total pressure drop, while the local pressure variations along the heat pipe are quite different for these two models when the radial Reynolds number at the interface is high.

Faghri, Amir↗

NASA NTRS1995

Introduction to multigrid methods

These notes were written for an introductory course on the application of multigrid methods to elliptic and hyperbolic partial differential equations for engineers, physicists and applied mathematicians. The use of more advanced mathematical tools, such as functional analysis, is avoided. The course is intended to be accessible to a wide audience of users of computational methods. We restrict ourselves to finite volume and finite difference discretization. The basic principles are given. Smoothing methods and Fourier smoothing analysis are reviewed. The fundamental multigrid algorithm is studied. The smoothing and coarse grid approximation properties are discussed. Multigrid schedules and structured programming of multigrid algorithms are treated. Robustness and efficiency are considered.

Wesseling, P.↗

NASA NTRS1994

Numerical Computation of Flame Spread over a Thin Solid in Forced Concurrent Flow with Gas-phase Radiation

Excerpts from a paper describing the numerical examination of concurrent-flow flame spread over a thin solid in purely forced flow with gas-phase radiation are presented. The computational model solves the two-dimensional, elliptic, steady, and laminar conservation equations for mass, momentum, energy, and chemical species. Gas-phase combustion is modeled via a one-step, second order finite rate Arrhenius reaction. Gas-phase radiation considering gray non-scattering medium is solved by a S-N discrete ordinates method. A simplified solid phase treatment assumes a zeroth order pyrolysis relation and includes radiative interaction between the surface and the gas phase.

Jiang, Ching-Biau↗

NASA NTRS1992

Computation of Turbulent Recirculating Flow in Channels, and for Impingement Cooling

Fully elliptic forms of the transport equations have been solved numerically for two flow configurations. The first is turbulent flow in a channel with transverse rectangular ribs, and the second is impingement cooling of a plane surface. Both flows are relevant to proposed designs for active cooling of hypersonic vehicles using supercritical hydrogen as the coolant. Flow downstream of an abrupt pipe expansion and of a backward-facing step were also solved with various near-wall turbulence models as benchmark problems. A simple form of periodicity boundary condition was used for the channel flow with transverse rectangular ribs. The effects of various parameters on heat transfer in channel flow with transverse ribs and in impingement cooling were investigated using the Yap modified Jones and Launder low Reynolds number k-epsilon turbulence model. For the channel flow, predictions were in adequate agreement with experiment for constant property flow, with the results for friction superior to those for heat transfer. For impingement cooling, the agreement with experiment was generally good, but the results suggest that improved modelling of the dissipation rate of turbulence kinetic energy is required in order to obtain improved heat transfer prediction, especially near the stagnation point. The k-epsilon turbulence model was used to predict the mean flow and heat transfer for constant and variable property flows. The effect of variable properties for channel flow was investigated using the same turbulence model, but comparison with experiment yielded no clear conclusions. Also, the wall function method was modified for use in the variable properties flow with a non-adiabatic surface, and an empirical model is suggested to correctly account for the behavior of the viscous sublayer with heating.

Chang, Byong Hoon↗

NASA NTRS1996

Homotopy Solutions of Kepler's Equations

Kepler's Equation is solved using an integrative algorithm developed using homotropy theory. The solution approach is applicable to both elliptic and hyperbolic forms of Kepler's Equation. The results from the proposed algorithm compare quite favorably with those from existing iterative schemes.

Fitz-Coy, Norman↗

NASA NTRS1984

Stress-intensity factor equations for cracks in three-dimensional finite bodies subjected to tension and bending loads

Stress intensity factor equations are presented for an embedded elliptical crack, a semielliptical surface crack, a quarter elliptical corner crack, a semielliptical surface crack along the bore of a circular hole, and a quarter elliptical corner crack at the edge of a circular hole in finite plates. The plates were subjected to either remote tension or bending loads. The stress intensity factors used to develop these equations were obtained from previous three dimensional finite element analyses of these crack configurations. The equations give stress intensity factors as a function of parametric angle, crack depth, crack length, plate thickness, and, where applicable, hole radius. The ratio of crack depth to plate thickness ranged from 0 to 1, the ratio of crack depth to crack length ranged from 0.2 to 2, and the ratio of hole radius to plate thickness ranged from 0.5 to 2. The effects of plate width on stress intensity variation along the crack front were also included.

Newman, J. C., Jr.↗

NASA NTRS1980

Fully flooded elastohydrodynamic lubricated elliptical contacts

Emphasis is on fully flooded, elastohydrodynamic lubricated, elliptical contacts. A fully flooded conjunction is one in which the film thickness is not significantly changed when the amount of lubricant is increased. A brief description of the relevant equations used in the elastohydrodynamic lubrication of elliptical contacts is given. The most important practical aspect of the elastohydrodynamic theory is the determination of the minimum film thickness within the contact. The maintenance of a fluid film of adequate magnitude is an essential feature of the correct operation of lubricated machine elements. The results presented show the influence of contact geometry on minimum film thickness as expressed by the ellipticity parameter and the dimensionless speed, load, and materials parameters. Film thickness equations are developed for materials of high elastic modulus, such as metal, and for materials of low elastic modulus, such as rubber. In addition to the film thickness equations that are developed, plots of pressure and film thickness are presented. These theoretical solutions for film thickness have all the essential features of previously reported experimental observations based on optical interferometry. Correlation between theory and experiments is also presented.

Hamrock, B. J.↗

NASA NTRS1995

Kepler Equation solver

Kepler's Equation is solved over the entire range of elliptic motion by a fifth-order refinement of the solution of a cubic equation. This method is not iterative, and requires only four transcendental function evaluations: a square root, a cube root, and two trigonometric functions. The maximum relative error of the algorithm is less than one part in 10(exp 18), exceeding the capability of double-precision computer arithmetic. Roundoff errors in double-precision implementation of the algorithm are addressed, and procedures to avoid them are developed.

Markley, F. Landis↗

NASA NTRS1986

Parameter estimation problems for distributed systems using a multigrid method

The problem of estimating spatially varying coefficients of partial differential equations is considered from observation of the solution and of the right hand side of the equation. It is assumed that the observations are distributed in the domain and that enough observations are given. A method of discretization and an efficient multigrid method for solving the resulting discrete systems are described. Numerical results are presented for estimation of coefficients in an elliptic and a parabolic partial differential equation.

Taasan, S.↗

NASA NTRS1990

Parameter estimation problems for distributed systems using a multigrid method

Ta'asan, Shlomo↗

NASA NTRS1976

Elliptic integral solutions to a class of space flight optimization problems

This paper is initially concerned with the minimum-time, exoatmospheric flight of a rocket with constant thrust acceleration magnitude, as in the cases of nuclear and solar electric propulsion. Gravitational acceleration is assumed to be a constant scalar multiple of the radius vector, plus a correction term which is a given function of time. The solution to the state equations is obtained in terms of elliptic integrals. A method is presented for the solution of the two-point boundary-condition problem associated with orbital transfer. At most, the latter method requires iteration upon final time, angle of injection, and two other parameters which are bounded. An example problem is provided which involves a rocket with very low thrust and a spiraling trajectory of many revolutions, but an altitude change of only several hundred miles above the earth. Finally, the original elliptic integral solution is extended to a larger class of low and intermediate thrust problems with constant thrust magnitude, mass decreasing with time, and an inverse square gravitational force.

Andrus, J. F.↗

NASA NTRS1988

Nonlinear, dispersive, elliptically polarized Alfven wavaes

The derivative nonlinear Schroedinger (DNLS) equation is derived by an efficient means that employs Lagrangian variables. An expression for the stationary wave solutions of the DNLS that contains vanishing and nonvanishing and modulated and nonmodulated boundary conditions as subcases is then obtained. The solitary wave solutions for elliptically polarized quasiparallel Alfven waves in the magnetohydrodynamic limit (nonvanishing, unmodulated boundary conditions) are obtained. These converge to the Korteweg-de Vries and the modified Korteweg-de Vries solitons obtained previously for oblique propagation, but are more general. It is shown that there are no envelope solitary waves if the point at infinity is unstable to the modulational instability. The periodic solutions of the DNLS are characterized.

Kennel, C. F.↗

NASA NTRS1996

First-Order System Least Squares for the Stokes Equations, with Application to Linear Elasticity

Following our earlier work on general second-order scalar equations, here we develop a least-squares functional for the two- and three-dimensional Stokes equations, generalized slightly by allowing a pressure term in the continuity equation. By introducing a velocity flux variable and associated curl and trace equations, we are able to establish ellipticity in an H(exp 1) product norm appropriately weighted by the Reynolds number. This immediately yields optimal discretization error estimates for finite element spaces in this norm and optimal algebraic convergence estimates for multiplicative and additive multigrid methods applied to the resulting discrete systems. Both estimates are uniform in the Reynolds number. Moreover, our pressure-perturbed form of the generalized Stokes equations allows us to develop an analogous result for the Dirichlet problem for linear elasticity with estimates that are uniform in the Lame constants.

Cai, Z.↗