Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hyperbolic systems”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Highly efficient antenna system using a corrugated horn and scanning hyperbolic reflector

In a horn-reflector antenna system for producing a spherical aperture phase front, a corrugated conical horn illuminates a section of a hyperbolic reflector to produce a spherical aperture phase front. This front produces a far field beam with low sidelobes and high beam efficiency. The system is insensitive to frequency and polarization changes, and is also insensitive to orientation about the axis of the conical horn for beam scanning.

Green, K. A.↗
NASA NTRS1981

On a fourth order accurate implicit finite difference scheme for hyperbolic conservation laws. II - Five-point schemes

This paper presents a family of two-level five-point implicit schemes for the solution of one-dimensional systems of hyperbolic conservation laws, which generalized the Crank-Nicholson scheme to fourth order accuracy (4-4) in both time and space. These 4-4 schemes are nondissipative and unconditionally stable. Special attention is given to the system of linear equations associated with these 4-4 implicit schemes. The regularity of this system is analyzed and efficiency of solution-algorithms is examined. A two-datum representation of these 4-4 implicit schemes brings about a compactification of the stencil to three mesh points at each time-level. This compact two-datum representation is particularly useful in deriving boundary treatments. Numerical results are presented to illustrate some properties of the proposed scheme.

Harten, A.↗
NASA NTRS1981

On a fourth order accurate implicit finite difference scheme for hyperbolic conservation laws. I - Nonstiff strongly dynamic problems

An implicit finite difference method of fourth order accuracy in space and time is introduced for the numerical solution of one-dimensional systems of hyperbolic conservation laws. The basic form of the method is a two-level scheme which is unconditionally stable and nondissipative. The scheme uses only three mesh points at level t and three mesh points at level t + delta t. The dissipative version of the basic method given is conditionally stable under the CFL (Courant-Friedrichs-Lewy) condition. This version is particularly useful for the numerical solution of problems with strong but nonstiff dynamic features, where the CFL restriction is reasonable on accuracy grounds. Numerical results are provided to illustrate properties of the proposed method.

Harten, A.↗
NASA NTRS2016

Hyperbolic Method for Dispersive PDEs: Same High-Order of Accuracy for Solution, Gradient, and Hessian

In this paper, we introduce a new hyperbolic first-order system for general dispersive partial differential equations (PDEs). We then extend the proposed system to general advection-diffusion-dispersion PDEs. We apply the fourth-order RD scheme of Ref. 1 to the proposed hyperbolic system, and solve time-dependent dispersive equations, including the classical two-soliton KdV and a dispersive shock case. We demonstrate that the predicted results, including the gradient and Hessian (second derivative), are in a very good agreement with the exact solutions. We then show that the RD scheme applied to the proposed system accurately captures dispersive shocks without numerical oscillations. We also verify that the solution, gradient and Hessian are predicted with equal order of accuracy.

Mazaheri, Alireza↗
DOE OSTI2026

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗
NASA NTRS1980

Dynamics of micrometeoroids

The spatial and orbital distributions of meteoroids of mass less than 10 to the -8 g are discussed. Pioneer 8 and 9, HEOS-2 and Helios 1 observations of the smallest meteoroids leaving the solar system in hyperbolic orbits due to radiation pressure (beta particles), intermediate mass particles arriving from the heliocentric spacecraft apex direction (alpha meteoroids) and the largest micrometeoroids observed with no preferred directionality (sporadic meteoroids) are reviewed, and evidence for the three groups of dynamically different micrometeoroids in lunar microcrater statistics is indicated. The collisional model of the origin of dust particles from the dispersion of sporadic meteoroids of cometary origin due to the Poynting-Robertson effect and collisions is then applied to account for the dynamics of the beta and alpha micrometeoroids.

Gruen, G.↗
NASA NTRS1987

Construction of explicit and implicit symmetric TVD schemes and their applications

A one-parameter family of explicit and implicit total variation diminishing (TVD) schemes is developed which permits incorporation of an expanded group of slope and flux limiters. The numerical technique is intended for use in calculations which include a time-differencing scheme and an optional Lax-Wendroff scheme. Methods of extending the TVD models to nonlinear scalar equations and systems of hyperbolic conservation equations are described. Sample results are presented from calculations of shocked flows around NACA 0012 and NACA 0018 airfoils.

Yee, H. C.↗
NASA NTRS1987

Extension of a semi-implicit shock-capturing algorithm for 3-D fully coupled, chemically reacting flows in generalized coordinates

A semiimplicit high-resolution shock-capturing method for multidimensional systems of hyperbolic conservation laws with stiff source terms has been developed by Yee and Shinn (1987). The goal of this work is to extend this method to solve the three-dimensional fully coupled Navier-Stokes equations for a hypersonic chemically reacting flow in generalized coordinates. In this formulation, the global continuity equation was replaced by all the species continuity equations. The shock-capturing technique is a second-order-accurate, symmetric total-variation-diminishing method which accounts fully and directly for the coupling among the fluid and all the species. To verify the current approach, it was implemented into an existing computer code which contained the MacCormack method. Test results for a five-species reacting flow are shown to be oscillation-free around the shock, and the time spent per iteration only doubles when compared to the result using classical way of supplying numerical dissipation. The extra computation is more than justified by the elimination of spurious oscillation and nonlinear instability associated with the classical shock-capturing schemes in computing hypersonic reacting flows.

Shinn, J. L.↗
NASA NTRS1987

Vectorizable implicit algorithms for the flux-difference split, three-dimensional Navier-Stokes equations

The computational efficiency of four vectorizable implicit algorithms is assessed when applied to calculate steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in general coordinates. Two of these algorithms are characterized as hybrid schemes; that is, they combine some approximate factorization in two coordinate directions with relaxation in the remaining spatial direction. The other two algorithms utilize an approximate factorization approach which yields two-factor algorithms for three-dimensional systems. All four algorithms are implemented in identical high-resolution upwind schemes for the flux-difference split Navier-Stokes equations. These highly nonlinear schemes are obtained by extending an implicit Total Variation Diminishing (TVD) scheme recently developed for linear one-dimensional systems of hyperbolic conservation laws to the three-dimensional Navier-Stokes equations. The computation of vortical flow over a sharp-edged, thin delta wing has been chosen as a common numerical test case. The convergence of the algorithms is discussed and the accuracy of the computed flow-field results is assessed. The validity of the present results are demonstrated by a comparison with experimental data.

Hartwich, P. M.↗
NASA NTRS1989

Recent progress on essentially non-oscillatory shock capturing schemes

An account is given of the construction of efficient implementations of 'essentially nonoscillatory' (ENO) schemes that approximate systems of hyperbolic conservation laws. ENO schemes use a local adaptive stencil to automatically obtain information from regions of smoothness when the solution develops discontinuities. Approximations employing ENOs can thereby obtain uniformly high accuracy to the very onset of discontinuities, while retaining a sharp and essentially nonoscillatory shock transition. For ease of implementation, ENO schemes applying the adaptive stencil concept to the numerical fluxes and employing a TVD Runge-Kutta-type time discretization are constructed.

Osher, Stanley↗
NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗
NASA NTRS1992

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗
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 NTRS1997

Guidance of Nonlinear Nonminimum-Phase Dynamic Systems

The first two years research work has advanced the inversion-based guidance theory for: (1) systems with non-hyperbolic internal dynamics; (2) systems with parameter jumps; (3) systems where a redesign of the output trajectory is desired; and (4) the generation of recovery guidance maneuvers.

Devasia, Santosh↗
NASA NTRS1999

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗
NASA NTRS2004

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗