SEARCH · Search NASA

Results for “hyperbolic problems”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗

NASA NTRS1983

A study of dynamic energy equations for Stirling cycle analysis

An analytical and computer study of the dynamic energy equations that describe the physical phenomena that occurs in a Stirling cycle engine. The basic problem is set up in terms of a set o hyperbolic partial differential equations. The characteristic lines are determined. The equations are then transformed to ordinary differential equations that are valid along characteristic lines. Computer programs to solve the differential equations and to plot pertinent factors are described.

Larson, V. H.↗

NASA NTRS1984

The time evolution of spectral discretizations of hyperbolic systems

A Chebyshev collocation spectral method, applied to hyperbolic systems is considered, particularly for those initial boundary value problems which possess only solutions tending to zero at large times. It is shown that the numerical solutions of the system will also vanish at infinity, if and only if, the numerical solution of a scalar equation of the same type does. This result is then generalized for other spectral approximations.

Lustman, L.↗

NASA NTRS1989

An Eigenvalue Analysis of finite-difference approximations for hyperbolic IBVPs

The eigenvalue spectrum associated with a linear finite-difference approximation plays a crucial role in the stability analysis and in the actual computational performance of the discrete approximation. The eigenvalue spectrum associated with the Lax-Wendroff scheme applied to a model hyperbolic equation was investigated. For an initial-boundary-value problem (IBVP) on a finite domain, the eigenvalue or normal mode analysis is analytically intractable. A study of auxiliary problems (Dirichlet and quarter-plane) leads to asymptotic estimates of the eigenvalue spectrum and to an identification of individual modes as either benign or unstable. The asymptotic analysis establishes an intuitive as well as quantitative connection between the algebraic tests in the theory of Gustafsson, Kreiss, and Sundstrom and Lax-Richtmyer L(sub 2) stability on a finite domain.

Warming, Robert F.↗

NASA NTRS1990

An eigenvalue analysis of finite-difference approximations for hyperbolic IBVPs

Warming, Robert F.↗

NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗

NASA NTRS1990

Consistency and convergence for numerical radiation conditions

The problem of imposing radiation conditions at artificial boundaries for the numerical simulation of wave propagation is considered. Emphasis is on the behavior and analysis of the error which results from the restriction of the domain. The theory of error estimation is briefly outlined for boundary conditions. Use is made of the asymptotic analysis of propagating wave groups to derive and analyze boundary operators. For dissipative problems this leads to local, accurate conditions, but falls short in the hyperbolic case. A numerical experiment on the solution of the wave equation with cylindrical symmetry is described. A unified presentation of a number of conditions which have been proposed in the literature is given and the time dependence of the error which results from their use is displayed. The results are in qualitative agreement with theoretical considerations. It was found, however, that for this model problem it is particularly difficult to force the error to decay rapidly in time.

Hagstrom, Thomas↗

NASA NTRS1973

Resonant oscillations of intermediate frequency in a stratified atmosphere.

A class of solutions to a model of forced oscillations in a rotating stratified atmospheric layer is derived and analyzed. The basic model is found to reduce to a boundary value problem with a second-order linear partial differential equation of the hyperbolic type for this range of forcing frequencies. The forced solutions are shown to exhibit resonances with the normal modes of oscillation of the layer. The characteristics of the resonant modes are analyzed in terms of mean tropospheric values of temperature, temperature lapse, wind speed, horizontal and vertical wind shears, latitude, and the frequency and horizontal wavelength of the forcing mechanism. These solutions are compared with solutions to the model for a different (subinertial) range of forcing frequencies. This comparison leads to an elliptic boundary value problem. The solutions in that case do not exhibit the same type of resonance and generally decay away from the region of forcing.

Thompson, O. E.↗

NASA NTRS1991

Finite-volume application of high-order ENO schemes to two-dimensional boundary-value problems

Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.

Casper, Jay↗

NASA NTRS1975

On the cost of approximating and recognizing a noise perturbed straight line or a quadratic curve segment in the plane

Approximation of noisy data in the plane by straight lines or elliptic or single-branch hyperbolic curve segments arises in pattern recognition, data compaction, and other problems. The efficient search for and approximation of data by such curves were examined. Recursive least-squares linear curve-fitting was used, and ellipses and hyperbolas are parameterized as quadratic functions in x and y. The error minimized by the algorithm is interpreted, and central processing unit (CPU) times for estimating parameters for fitting straight lines and quadratic curves were determined and compared. CPU time for data search was also determined for the case of straight line fitting. Quadratic curve fitting is shown to require about six times as much CPU time as does straight line fitting, and curves relating CPU time and fitting error were determined for straight line fitting. Results are derived on early sequential determination of whether or not the underlying curve is a straight line.

Cooper, D. B.↗

NASA NTRS1979

Sound radiation from hyperboloidal inlet ducts

This paper presents rigorous solutions for the problem of sound radiation from various inlet ducts including hyperboloidal (or hyperbolic) inlet ducts and circular ducts with wide flange. The numerical results include the complex conversion (or reflection) coefficients and the radiation directivity for the various incident wave modes - spinning modes as well as axisymmetric modes. The analysis utilizes hyperboloidal wave functions which are defined here as a class of eigensolutions of the wave equation for oblate spheroidal coordinates, and is valid for the whole frequency range including frequencies above and below the cutoff frequencies of duct modes involved.

Cho, Y. C.↗

NASA NTRS1980

On the practical use of high-order methods for hyperbolic systems

The paper tests a number of high order methods on a variety of dynamic problems in one, two, and three space dimensions. The problems covered include wave propagation phenomena as well as an asymptotic approach to a steady state. Consideration is given to both smooth and shocked flows. It is shown that the methods compared require only minor modifications of many existing second-order schemes. Further, the results show that significant gains can be expected from the use of fourth-order methods. Finally, spectral methods are also considered for some of the problems presented.

Turkel, E.↗

NASA NTRS1988

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed utilizing a multigrid iteration technique. The method is compared with a noniterative approach in terms of computational time and accuracy. The computational model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobsson-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via a smeared mass or striated flow extending to both surfaces in the film gap. The mixed nature of the equations (parabolic in the full film zone and hyperbolic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗

NASA NTRS1989

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

Woods, C. M.↗

NASA NTRS1992

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗

NASA NTRS1991

A staggered mesh finite difference scheme for the computation of hypersonic Euler flows

A shock capturing finite difference method for systems of hyperbolic conservation laws is presented which avoids the need to solve Riemann problems while being competitive in performance with other current methods. A staggered spatial mesh is employed, so that complicated nonlinear waves generated at cell interfaces are averaged over cell interiors at the next time level. The full method combines to form a conservative version of the modified method of characteristics. The advantages of the method are discussed, and numerical results are presented for the two-dimensional double ellipse problem.

Sanders, Richard↗

NASA NTRS1995

Second-order numerical solution of time-dependent, first-order hyperbolic equations

A finite difference scheme is developed to find an approximate solution of two similar hyperbolic equations, namely a first-order plane wave and spherical wave problem. Finite difference approximations are made for both the space and time derivatives. The result is a conditionally stable equation yielding an exact solution when the Courant number is set to one.

Shah, Patricia L.↗

NASA NTRS1959

An Analysis of the Corridor and Guidance Requirements for Supercircular Entry into Planetary Atmospheres

The analysis of super-circular entry is developed around a new dimensionless parameter which combines certain conditions at the conic perigee altitude with certain characteristics of the vehicle; this parameter conveniently determines either deceleration-limited or heating limited corridor widths for elliptic, parabolic, or hyperbolic approach trajectories. Illustrative calculations of corridor widths and the associated guidance problems are presented for Venus, Earth, Mars, Jupiter, and Titan. Generalized curves are presented for application to various entry conditions.

Chapman, Dean R.↗