SEARCH · Engineering Papers

Results for “hyperbolic”

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 847 records · Page 47

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.

Harten, A.↗

NASA NTRS1986

Regularity estimates up to the boundary for elliptic systems of difference equations

Regularity estimates up to the boundary for solutions of elliptic systems of finite difference equations were proved. The regularity estimates, obtained for boundary fitted coordinate systems on domains with smooth boundary, involve discrete Sobolev norms and are proved using pseudo-difference operators to treat systems with variable coefficients. The elliptic systems of difference equations and the boundary conditions which are considered are very general in form. The regularity of a regular elliptic system of difference equations was proved equivalent to the nonexistence of eigensolutions. The regularity estimates obtained are analogous to those in the theory of elliptic systems of partial differential equations, and to the results of Gustafsson, Kreiss, and Sundstrom (1972) and others for hyperbolic difference equations.

Strikwerda, J. C.↗

NASA NTRS1986

Navier-Stokes computations for exotic airfoils

A hyperbolic grid generator with improvements for handling sharp corners and concave surfaces combined with a Navier-Stokes flow solver is applied to complex two dimensional airfoil configurations. Steady separated flow about an iced leading edge of an airfoil is presented. Unsteady viscous separated flows past an airfoil at two angles of attack with a spoiler deployed at 60 deg are compared with experiment. The spoiler computations are performed with two different topological maps of the physical domain to the computational domain. Graphical techniques for the static and unsteady display of the flow fields are discussed.

Barth, T. J.↗

NASA NTRS1985

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching-step sizes with the implicit shock fitting procedure. A correlation is found between the marching-step size, Reynolds number and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗

NASA NTRS1986

A new approximate LU factorization scheme for the Reynolds-averaged Navier-Stokes equations

A new approximate LU factorization scheme is developed to solve the steady state Reynolds-averaged Navier-Stokes (NS) equations. Central differencing is used for both implicit and explicit operator and special care is taken to obtain well-conditioned factors on the implicit side. The scheme is then analyzed and optimized according to a simple linear analysis. It is unconditionally stable for the model hyperbolic equation in both two- and three-dimensions. However, the requirement for well-conditioned factors has essentially limited the effective time step the scheme can achieve. Supersonic and transonic three-dimensinal flows past a hemisphere cylinder are computed to demonstrate the convergence characteristics of the scheme. A good convergence rate is achieved for the inviscid case. Finally, an explicit eigenvector annihilation procedure is adopted successfully to remove the stiffness caused by the fine grid spacing for viscous flows.

Pan, D.↗

NASA NTRS1985

Stability of finite-difference models containing two boundaries or interfaces

The stability of finite difference models of hyperbolic initial boundary value problems is connected with the propagation and reflection of parasitic waves. Wave propagation ideas are applied to models containing two boundaries or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Various known instability phenomena are accounted for in a unified way. Results show: (1) dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface; (2) algebraic GKS instabilities can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present; and (3) GKS-stability and P-stability can be established in certain problems by showing that all numerical reflection coefficients have modulus less than 1.

Trefethen, L. N.↗

NASA NTRS1985

Topics in spectral methods

After detailing the construction of spectral approximations to time-dependent mixed initial boundary value problems, a study is conducted of differential equations of the form 'partial derivative of u/partial derivative of t = Lu + f', where for each t, u(t) belongs to a Hilbert space such that u satisfies homogeneous boundary conditions. For the sake of simplicity, it is assumed that L is an unbounded, time-independent linear operator. Attention is given to Fourier methods of both Galerkin and pseudospectral method types, the Galerkin method, the pseudospectral Chebyshev and Legendre methods, the error equation, hyperbolic partial differentiation equations, and time discretization and iterative methods.

Gottlieb, D.↗

NASA NTRS1986

On a large time-step high resolution scheme

This paper presents a class of new second-order accurate (2K + 3)-point explicit schemes for the computation of weak solutions of hyperbolic conservation laws, that are total-variation-diminishing under a Courant-Friedrichs-Lewy restriction of K. These highly nonlinear schemes are obtained by applying a nonoscillatory first-order accurate (2K + 1)-point scheme to a modified flux. The derived second-order accurate schemes achieve high resolutions, while retaining the robustness of the original first-order accurate scheme.

Harten, A.↗

NASA NTRS1986

Technology and the Hubble Space Telescope

The components of the Hubble Space Telescope which is a 13.1 m long, 4.27 m in diameter, and weighs 11,000 kg are described. The telescope will be placed in a 600-km circular orbit of 28.5 deg inclination and it has an expected service life of 15 years. The optical telescope assembly contains the hyperbolic Ritchey-Chretin type primary and secondary mirrors and the optical control system. The telescope is to have a field of view of 28 arcmin in diameter. The use of the fine guidance sensors for pointing stability is examined. The optical and scientific system of the telescope are contained in the support system module. A wide-field/planetary camera, ESA-furnished faint object camera, the high-speed photometer, faint-object spectrograph, high resolution spectrograph, and ESA-furnished solar arrays are located on the telescope. The maintenance and refurbishment of the telescope are discussed. The precise lambda/20 optics, stringent cleanliness, and stable pointing of the telescope produce good image quality for astronomical and solar system science studies.

Mitchell, R. E.↗

NASA NTRS1985

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

NASA NTRS1985

A high-precision finite element method for shock-tube calculations

A two-pass explicit scheme is developed in order to exploit some of the capabilities of finite difference modeling (FDM) for finite element modeling (FEM), which offers the opportunity to account for any type of geometry in fluid flow modeling. Features of the first-order upwind and the Lax-Wendroff high precision explicit finite difference algorithms are reviewed. A flux limiter is developed for FEM to serve as an analog for the single limiter function which has been defined for the various FDMs. It is shown that an antidiffusive limiter must be introduced into the weighting function which normally multiplies the time-derivative term in the variational equation. The two-pass scheme which results is demonstrated to be the equivalent of FDMs with five-point support. However, the present scheme is valid only for one-dimensional calculations and linear shape functions for shock tube flow phenomena. Further work is required for its use with nonlinear hyperbolic systems.

Hughes, T. J. R.↗

NASA NTRS1985

The LDEF interplanetary dust experiment

Explorer 46 data and the Long Duration Exposure Facility (LDEF) interplanetary dust experiment are examined. Analysis of the Explorer 46 data reveals the existence of particles of 0.1 micron and a mass of 1 x 10 to the -16th gm, the injection of the submicron particles directly by a comet (injection mass of about 5000 tons), and a submicron particle lifetime of about three years. The applications of LDEF data to particles in hyperbolic orbits, particle swarms, morning-to-evening asymmetry, the effects of the earth's orbit eccentricity, and the presence of interstellar dust are discussed. The effects of space debris on data collection are considered.

Singer, S. F.↗

NASA NTRS1985

Mass input into and output from the meteoritic complex

The consequences which arise from the mutual collisions occurring between interplanetary meteoroids, the Poynting-Robertson (PR) effect and the radiation pressure ejection of small meteoroids are examined. The size distribution and flux of micrometeoroids at 1 AU are derived and the dependence of spatial density on distance from the sun is established. The following conclusions are made: (1) the lifetimes of meteoroids with masses approximately greater than 0.00001 g are dominated by catastrophic collisions; (2) after bering crushed by collisions, 70 to 85 percent of this mass will be in the form of zodiacal light particles (with masses in the range of 10 to the -10th g to 10 to the -5th g) which will in part be transported by the PR effect towards the sun where they will evaporate; (3) the 15 to 30 percent of the collisional fragments which have masses approximately less than 10 to the -10th g will, for the most part, be injected into hyperbolic orbits by radiation pressure.

Gruen, E.↗

NASA NTRS1985

Information content in spectral calculations

Analytical procedures for extracting piecewise smooth solutions of hyperbolic systems from raw oscillatory data obtained by pseudospectral methods are developed. The validity of the approach is demonstrated for the case of linear problems with constant coefficients, and plausibility arguments are presented which indicate its applicability to nonlinear operators when the steady state has been achieved. Numerical results for the development of an oblique shock when a wedge is inserted at zero angle of attack into a uniform supersonic flow of an ideal gas (the time-dependent two-dimensional Euler equations discretized in space by the pseudospectral Chebyshev method) are presented in tables and graphs.

Abarbanel, S.↗

NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗

NASA NTRS1986

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme.

Solomonoff, A.↗

NASA NTRS1986

On the Maxwellian distribution, symmetric form, and entropy conservation for the Euler equations

The Euler equations of gas dynamics have some very interesting properties in that the flux vector is a homogeneous function of the unknowns and the equations can be cast in symmetric hyperbolic form and satisfy the entropy conservation. The Euler equations are the moments of the Boltzmann equation of the kinetic theory of gases when the velocity distribution function is a Maxwellian. The present paper shows the relationship between the symmetrizability and the Maxwellian velocity distribution. The entropy conservation is in terms of the H-function, which is a slight modification of the H-function first introduced by Boltzmann in his famous H-theorem. In view of the H-theorem, it is suggested that the development of total H-diminishing (THD) numerical methods may be more profitable than the usual total variation diminishing (TVD) methods for obtaining wiggle-free solutions.

Deshpande, S. M.↗

NASA NTRS1986

A spline-based parameter estimation technique for static models of elastic structures

The problem of identifying the spatially varying coefficient of elasticity using an observed solution to the forward problem is considered. Under appropriate conditions this problem can be treated as a first order hyperbolic equation in the unknown coefficient. Some continuous dependence results are developed for this problem and a spline-based technique is proposed for approximating the unknown coefficient, based on these results. The convergence of the numerical scheme is established and error estimates obtained.

Dutt, P.↗