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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Numerical calculations of flow fields

The solutions to the equations of motion for inviscid fluid flow around a pointed elliptic cone at incidence are presented. The numerical method used, MacCormack's second order preferential predictor-corrector finite difference approximation, is applied to the fluid flow equations derived in conservation-law form. The entropy boundary condition, hitherto unused for elliptic cone problems, is investigated and compared to reflection boundary condition solutions. The stagnation streamline movement of the inclined elliptic cone is noted and surface pressure coefficients are plotted. Also presented are solutions for an elliptic cone and a circular cone at zero incidence and a circular cone at a small angle of attack. Comparisons are made between these present solutions and previously published theory.

Anderson, D. M.↗
NASA NTRS1974

Autonomous orbital navigation using Kepler's equation

A simple method of determining the six elements of elliptic satellite orbits has been developed for use aboard manned and unmanned spacecraft orbiting the earth, moon, or any planet. The system requires the use of a horizon sensor or other device for determining the local vertical, a precision clock or timing device, and Apollo-type navigation equipment including an inertial measurement unit (IMU), a digital computer, and a coupling data unit. The three elements defining the in-plane motion are obtained from simultaneous measurements of central angle traversed around the planet and elapsed flight time using a linearization of Kepler's equation about a reference orbit. It is shown how Kalman filter theory may also be used to determine the in-plane orbital elements. The three elements defining the orbit orientation are obtained from position angles in celestial coordinates derived from the IMU with the spacecraft vertically oriented after alignment of the IMU to a known inertial coordinate frame.

Boltz, F. W.↗
NASA NTRS1996

Multigrid Techniques for Highly Indefinite Equations

A multigrid method for the solution of finite difference approximations of elliptic PDE's is introduced. A parallelizable version of it, suitable for two and multi level analysis, is also defined, and serves as a theoretical tool for deriving a suitable implementation for the main version. For indefinite Helmholtz equations, this analysis provides a suitable mesh size for the coarsest grid used. Numerical experiments show that the method is applicable to diffusion equations with discontinuous coefficients and highly indefinite Helmholtz equations.

Shapira, Yair↗
NASA NTRS1993

Multigrid methods for differential equations with highly oscillatory coefficients

New coarse grid multigrid operators for problems with highly oscillatory coefficients are developed. These types of operators are necessary when the characters of the differential equations on coarser grids or longer wavelengths are different from that on the fine grid. Elliptic problems for composite materials and different classes of hyperbolic problems are practical examples. The new coarse grid operators can be constructed directly based on the homogenized differential operators or hierarchically computed from the finest grid. Convergence analysis based on the homogenization theory is given for elliptic problems with periodic coefficients and some hyperbolic problems. These are classes of equations for which there exists a fairly complete theory for the interaction between shorter and longer wavelengths in the problems. Numerical examples are presented.

Engquist, Bjorn↗
NASA NTRS1973

Wind-driven currents in a shallow lake or sea.

A technique is given for obtaining exact analytical solutions for wind-driven currents in a shallow lake or sea with a large class of bottom topographies which may or may not contain an island. Welander's second-order partial differential equation for the lake's surface displacement is used for determining the lake circulation. The technique, together with some results based on conformal mapping, is used to reduce the problem of finding solutions to Welander's equation for a closed body of water to the problem of solving an ordinary differential equation. The method is applied to an elliptically shaped lake and a circular lake containing an eccentrically located circular island.

Goldstein, M.↗
NASA NTRS1980

Effects of supply conditions on film thickness in lubricated Hertzian contacts

A generalization of the hydrodynamic expression for Hertzian contacts is described and various methods for calculating the thickness of the oil film winter steady-state, isothermal conditions are given. This is important for engineering applications such as gears and bearings because these results are closer to real operating conditions. Theories of lubrication are discussed, and the mathematics involved are presented using approximately 30 equations and 13 figures. For lubricated, linear, elliptical or point Hertzian contacts it is demonstrated how to calculate the thickness of the oil film at the center of the contact for steady-state isothermal conditions.

Dalmaz, G.↗
NASA NTRS1981

A general time element for orbit integration in Cartesian coordinates

Two techniques are discussed for increasing the accuracy of the numerical integration of eccentric orbits in Cartesian coordinates. One involves the use of an independent variable different from time; this increases the efficiency of the numerical integration. The other uses a time element, which reduces the in-track error. A general expression is given of a time element valid for an arbitrary independent variable. It is pointed out that this time element makes it possible to switch the independent variable merely by applying a scaling factor; there is no need to change the differential equations of the motion. Eccentric, true, and elliptic anomalies are used as independent variables in the case of a transfer orbit for a geosynchronous orbit. The elliptic anomaly is shown to perform much better than the other classical anomalies.

Janin, G.↗
NASA NTRS1982

Recent applications of the transonic wing analysis computer code, TWING

An evaluation of the transonic-wing-analysis computer code TWING is given. TWING utilizes a fully implicit approximate factorization iteration scheme to solve the full potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations were analyzed, and the limits of applicability of this code was evaluated. Comparisons of computed results were made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative full potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.

Subramanian, N. R.↗
NASA NTRS1983

The efficient solution of transonic wing flow fields

An evaluation of the transonic-wing-analysis computer code TWING is presented. TWING utilizes a fully implicit, approximate-factorization iteration scheme to solve the full-potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations have been analyzed, and comparisons of computed results have been made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative, full-potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations, including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.

Holst, T. L.↗
NASA NTRS1990

On the three-dimensional instability of strained vortices

The three-dimensional (3-D) instability of a two-dimensional (2-D) flow with elliptical streamlines has been proposed as a generic mechanism for the breakdown of many 2-D flows. A physical interpretation for the mechanism is presented together with an analytical treatment of the problem. It is shown that the stability of an elliptical flow is governed by an Ince equation. An analytical representation for a localized solution is given and establishes a direct link with previous computations and experiments.

Waleffe, Fabian↗
NASA NTRS1990

Direct simulation for the instability and breakup of laminar liquid jets

A direct numerical simulation method is described for predicting the deformation of laminar liquid jets. In the present nonlinear direct simulation, the convective term, which has been discarded in past linear analyses by Rayleigh and others, is included in the hydrodynamic equations. It is shown that only by maintaining full complexity of the nonlinear surface tension term accurate drop formation can be predicted. The continuity and momentum equations in the transient form are integrated on an adaptive grid, conforming the jet and surface wave shape. The equations, which are parabolic in time and elliptic in space, are solved by a TVD scheme with characteristic flux splitting. The results of the present work are discussed and compared with available measurements and other analyses. The comparison shows that among the predictions, the current 1-D direct simulation results agree best with the experimental data. Furthermore, the computer time requirements are much (an order of magnitude) smaller than those of previously reported multidimensional analyses.

Chuech, S. G.↗
NASA NTRS1992

On the use of satellites in Molniya orbit for meteorological and oceanographic observations of the high latitudes

Two types of orbits are currently used for meteorological satellites: geostationary orbit and low earth orbit. Neither orbit serves the high latitudes well. The Soviet Molniya satellites fly in a highly elliptic orbit inclined 63.4 deg to the equator. For approximately 8 h of its 12-h orbit, a satellite in Molniya orbit is synchronized with the earth such that it is nearly geostationary in the high latitudes. A meteorological/oceanographic satellite in Molniya orbit would bring nearly the same frequent imaging capability to the high latitudes which the tropics and mid-latitudes now enjoy. The uses of such a satellite include Arctic forecasting, study of polar lows, high latitude precipitation estimation, Antarctic studies, sea ice monitoring, and ozone measurements.

Kidder, Stanley Q.↗
NASA NTRS1993

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗
NASA NTRS1993

A least-squares finite element method for 3D incompressible Navier-Stokes equations

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations, and results in symmetric, positive definite algebraic system. An additional compatibility equation, i.e., the divergence of vorticity vector should be zero, is included to make the first-order system elliptic. The Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. The flow in a half of 3D cubic cavity is calculated at Re = 100, 400, and 1,000 with 50 x 52 x 25 trilinear elements. The Taylor-Gortler-like vortices are observed at Re = 1,000.

Jiang, Bo-Nan↗
NASA NTRS1990

Impact of Non-geostationary Orbits on PASS

The use of satellites in non-geostationary orbits (NGO) for PASS (Personal Access Satellite System) is discussed. The following subject areas are covered: (1) orbit parameters (circular and elliptical orbit characteristics); (2) link characteristics (Doppler shift, propagation loss, CONUS coverage antenna, multibeam antenna gain, and impact on link equations); (3) number of satellites required for continuous CONUS coverage (circular and elliptical orbit results); and (4) advantages and disadvantages of NGOs.

Estabrook, Polly↗
NASA NTRS1994

Theoretical study of the incompressible Navier-Stokes equations by the least-squares method

Usually the theoretical analysis of the Navier-Stokes equations is conducted via the Galerkin method which leads to difficult saddle-point problems. This paper demonstrates that the least-squares method is a useful alternative tool for the theoretical study of partial differential equations since it leads to minimization problems which can often be treated by an elementary technique. The principal part of the Navier-Stokes equations in the first-order velocity-pressure-vorticity formulation consists of two div-curl systems, so the three-dimensional div-curl system is thoroughly studied at first. By introducing a dummy variable and by using the least-squares method, this paper shows that the div-curl system is properly determined and elliptic, and has a unique solution. The same technique then is employed to prove that the Stokes equations are properly determined and elliptic, and that four boundary conditions on a fixed boundary are required for three-dimensional problems. This paper also shows that under four combinations of non-standard boundary conditions the solution of the Stokes equations is unique. This paper emphasizes the application of the least-squares method and the div-curl method to derive a high-order version of differential equations and additional boundary conditions. In this paper, an elementary method (integration by parts) is used to prove Friedrichs' inequalities related to the div and curl operators which play an essential role in the analysis.

Jiang, Bo-Nan↗