Search NASA⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Application of two-point implicit central-difference methods to hyperbolic systems

This paper presents a general solution algorithm for the set of difference equations that arise when two-point central differences are used to approximate the flux difference terms in systems of hyperbolic differential equations. The general algorithm eliminates the weak points associated with the nonstandard algorithm reported by Wornom and Hafez (1986). The disadvantages of their algorithm relate to its implementation. It consists of separate algorithms for subsonic, supersonic, sonic and shock cells, applied individually, which presents a major bookkeeping problem when multiple sonic and shock cells are present. The general algorithm eliminates this problem and introduces an improved shock treatment which produces shocks with at most one interior shock point.

Wornom, Stephen F.↗
DOE OSTI2015

A hyperbolic-equation system approach for magnetized electron fluids in quasi-neutral plasmas

A new approach using a hyperbolic-equation system (HES) is proposed to solve for the electron fluids in quasi-neutral plasmas. The HES approach avoids treatments of cross-diffusion terms which cause numerical instabilities in conventional approaches using an elliptic equation (EE). A test calculation reveals that the HES approach can robustly solve problems of strong magnetic confinement by using an upwind method. The computation time of the HES approach is compared with that of the EE approach in terms of the size of the problem and the strength of magnetic confinement. The results indicate that the HES approach can be used to solve problems in a simple structured mesh without increasing computational time compared to the EE approach and that it features fast convergence in conditions of strong magnetic confinement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Multidomain approach for calculating compressible flows

A multidomain approach for calculating compressible flows by using unsteady or pseudo-unsteady methods is presented. This approach is based on a general technique of connecting together two domains in which hyperbolic systems (that may differ) are solved with the aid of compatibility relations associated with these systems. Some examples of this approach's application to calculating transonic flows in ideal fluids are shown, particularly the adjustment of shock waves. The approach is then applied to treating a shock/boundary layer interaction problem in a transonic channel.

Cambier, L.↗
DOE OSTI2018

Spectral Analysis of Linear Models of Viscoelasticity

In this paper, we examine Volterra integrodifferential equations with unbounded operator coefficients in Hilbert spaces. Equations considered are abstract hyperbolic equations perturbed by terms containing Volterra integral operators. These equations can be realized as partial integrodifferential equations that appear in the theory of viscoelasticity (see [2, 5]), as Gurtin–Pipkin integrodifferential equations (see [1, 7]) that describe finite-speed heat transfer in materials with memory. They also appear in averaging problems for multiphase media (Darcy’s law.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exact Timestep for a Pairwise Coulomb Collision

Standard numerical integrators work well for many-body Coulomb repulsion problems when the timestep is much shorter than the timescale of relative position changes. However, for ‘hard’ collisions in which two particles have a near miss and exchange a lot of momentum within one timestep, they understandably struggle. This note proposes using the exact solution of Keplerian two-body orbits (usually hyperbolic) to calculate the momentum exchange with other particles: either a selection of the ‘closest’ ones or all of them.

43 PARTICLE ACCELERATORS↗
NASA NTRS1991

The direct simulation of high-speed mixing-layers without and with chemical heat release

A direct numerical simulation of high speed reacting and non-reacting flows for H2-air systems is presented. The calculations are made for a convective Mach number of 0.38 with hyperbolic tangent initial profile and finite rate chemical reactions. A higher-order numerical method is used in time accurate mode to time advance the solution to a statistical steady state. About 600 time slices of all the variables are then stored for statistical analysis. It is shown that most of the problems of high-speed combustion with air are characterized by relatively weak heat release. The present study shows that: (1) the convective speed is reduced by heat release by about 10 percent at this convective Mach number M(sub c) = 0.38; (2) the variation of the mean and rms fluctuation of temperature can be explained on the basis of temperature fluctuation between the flame temperature and the ambient; (3) the growth rate with heat release is reduced by 7 percent; and (4) the entrainment is reduced by 25 percent with heat release. These differences are small in comparison with incompressible flow dynamics, and are argued to be due to the reduced importance of heat release in comparison with the large enthalpy gradients resulting from the large-scale vortex dynamics. It is finally suggested that the problems of reduced mixing in high-speed flows are not severely complicated by heat release.

Sekar, B.↗
DOE OSTI2014

Analysis of self-similar solutions of multidimensional conservation laws

This project focused on analysis of multidimensional conservation laws, specifically on extensions to the study of self-siminar solutions, a project initiated by the PI. In addition, progress was made on an approach to studying conservation laws of very low regularity; in this research, the context was a novel problem in chromatography. Two graduate students in mathematics were supported during the grant period, and have almost completed their thesis research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Analysis and Entropy Stability of the Line-Based Discontinuous Galerkin Method

We develop a discretely entropy-stable line-based discontinuous Galerkin method for hyperbolic conservation laws based on a flux differencing technique. By using standard entropy-stable and entropy-conservative numerical flux functions, this method guarantees that the discrete integral of the entropy is non-increasing. This nonlinear entropy stability property is important for the robustness of the method, in particular when applied to problems with discontinuous solutions or when the mesh is under-resolved. This line-based method is significantly less computationally expensive than a standard DG method. Numerical results are shown demonstrating the effectiveness of the method on a variety of test cases, including Burgers’ equation and the Euler equations, in one, two, and three spatial dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

High Order Discontinuous Galerkin Method and Applications (Final Technical Report)

This project aims at the developments and improvements of high order accurate discontinuous Galerkin finite element methods for solving partial differential equations arising from DOE applications. Algorithm development, improvement, analysis, implementation and applications aspects are addressed. The partial differential equations under study involve solutions with shocks or other discontinuities, rapidly changing regions, or multiscale solutions. These include the so-called convection dominated cases, such as the hyperbolic equations (e.g. Euler equations of gas dynamics), parabolic equations with small viscosities (e.g. high Reynolds number Navier-Stokes equations), KdV type equations with small dispersion describing water waves (small third derivative terms), time-dependent bi-harmonic equations and equations with hyper-viscosity (small fourth or higher even derivative terms), problems with highly oscillatory coefficients, or problems for which the macro-scale PDE models may not be correct in localized regions and a micro-scale model must be used instead. Furthermore, we are interested in the long time solutions of such equations with waves traveling a long distance, such as the electromagnetic waves. The requirement on numerical methods is thus high order accuracy, low dissipation/dispersion errors over long time, and a clean and sharp discontinuity resolution. It would also be highly desirable to have the numerical methods to be flexible for complicated geometry and boundary conditions, and to be easy in an adaptive and massively parallel environment. The discontinuous Galerkin finite element method has this potential, however the design and analysis of the discontinuous Galerkin method for each specific class of partial differential equations would require detailed research.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Composite methods for hyperbolic equations

A composite approximation procedure combining the properties of the Lax-Wendroff and leapfrog algorithms is proposed for solving hyperbolic equations. For a one-dimensional equation, a three-step approximation consisting of a two-step Richtmeyer method followed by a leapfrog step is considered. This is a two-level scheme, so all difficulties, including storage requirements, associated with the three-level leapfrog are eliminated. For two-dimensional problems a generalization of the preceding method is used consisting of a rotated Richtmeyer method followed by a modified leapfrog step. It is found that the composite schemes are effective in reducing oscillations and nonlinear instabilities that affect the leapfrog method. The dissipation in the composite schemes is much less than in the Richtmeyer algorithm, and hence can be used for long term integrations.

Turkel, E.↗
DOE OSTI2016

A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems

In this paper a new, simple and universal formulation of the HLLEM Riemann solver (RS) is proposed that works for general conservative and non-conservative systems of hyperbolic equations. For non-conservative PDE, a path-conservative formulation of the HLLEM RS is presented for the first time in this paper. The HLLEM Riemann solver is built on top of a novel and very robust path-conservative HLL method. It thus naturally inherits the positivity properties and the entropy enforcement of the underlying HLL scheme. However, with just the slight additional cost of evaluating eigenvectors and eigenvalues of intermediate characteristic fields, we can represent linearly degenerate intermediate waves with a minimum of smearing. For conservative systems, our paper provides the easiest and most seamless path for taking a pre-existing HLL RS and quickly and effortlessly converting it to a RS that provides improved results, comparable with those of an HLLC, HLLD, Osher or Roe-type RS. This is done with minimal additional computational complexity, making our variant of the HLLEM RS also a very fast RS that can accurately represent linearly degenerate discontinuities. Our present HLLEM RS also transparently extends these advantages to non-conservative systems. For shallow water-type systems, the resulting method is proven to be well-balanced. Several test problems are presented for shallow water-type equations and two-phase flow models, as well as for gas dynamics with real equation of state, magnetohydrodynamics (MHD & RMHD), and nonlinear elasticity. Since our new formulation accommodates multiple intermediate waves and has a broader applicability than the original HLLEM method, it could alternatively be called the HLLI Riemann solver, where the “I” stands for the intermediate characteristic fields that can be accounted for. -- Highlights: •New simple and general path-conservative formulation of the HLLEM Riemann solver. •Application to general conservative and non-conservative hyperbolic systems. •Inclusion of sub-structure and resolution of intermediate characteristic fields. •Well-balanced for single- and two-layer shallow water equations and multi-phase flows. •Euler equations with real equation of state, MHD equations, nonlinear elasticity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Transformation of two and three-dimensional regions by elliptic systems

The research during this period continued to expand the class of numerical algorithms that can be accurately and efficiently implemented on overlapping grids. Whereas previous calculations have been used to solve elliptic equations and to find the steady-state solution of parabolic equations, the present work is aimed towards developing time-accurate solution techniques for parabolic and hyperbolic equations. The primary difficulty here is in the correct treatment of the interior boundary nodes that must be updated at each iteration. The implementation of explicit methods is straightforward. However, the common practice of lagging these values when using an implicit methods leads to inconsistencies in the difference equation. One way to avoid this problem is to alternately calculate with an implicit and an explicit method on each subgrid. With this procedure, the explicit method generates boundary values at the next time level which are then used by the implicit step. It can be shown that when a backward implicit method is combined with a forward explicit method, the composite method is second order accurate and unconditionally stable for linear problems. A second area in which progress can be reported is in the distribution of grid points on curves and surfaces.

Mastin, E. W.↗
NASA NTRS1997

The effects of atmospheric [CO2] on carbon isotope fractionation and magnesium incorporation into biogenic marine calcite

The influences of atmospheric carbon dioxide on the fractionation of carbon isotopes and the magnesium incorporation into biogenic marine calcite were investigated using samples of the calcareous alga Amphiroa and benthic foraminifer Sorites grown in the Biosphere 2 Ocean system under variable atmospheric CO2 concentrations (approximately 500 to 1200 ppm). Carbon isotope fractionation was studied in both the organic matter and the skeletal carbonate. Magnesium analysis was to be performed on the carbonate removed during decalcification. These data have not been collected due to technical problems. Carbon isotope data from Amphiroa yields a linear relation between [CO2] and Delta(sup 13)C(sub Corg)values suggesting that the fractionation of carbon isotopes during photosynthesis is positively correlated with atmospheric [CO2]. [CO2] and Delta(sup 13)C(sub Corg) values for Sorites produce a relation that is best described by a hyperbolic function where Delta(sup 13)C(sub Corg) values increase between 300 and 700 ppm and decrease from 700 to 1200 ppm. Further investigation of this relation and Sorites physiology is needed.

Vieira, Veronica↗
NASA NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗
NASA NTRS1975

Diagonalization and simultaneous symmetrization of the gas-dynamic matrices

The hyperbolic nature of the unsteady, inviscid, gas-dynamic equations implies the existence of a similarity transformation for diagonalizing an arbitrary linear combination of coefficient matrices. It is shown that the individual matrices are simultaneously symmetrized by the similarity transformation. The transformations and their norms can be applied to the well-posedness of the Cauchy problem, linear stability theory for finite-difference approximations, and simplification of block-tridiagonal systems that arise in implicit time-split algorithms.

Warming, R. F.↗
NASA NTRS1978

Finite element concepts in computational aerodynamics

Finite element theory was employed to establish an implicit numerical solution algorithm for the time averaged unsteady Navier-Stokes equations. Both the multidimensional and a time-split form of the algorithm were considered, the latter of particular interest for problem specification on a regular mesh. A Newton matrix iteration procedure is outlined for solving the resultant nonlinear algebraic equation systems. Multidimensional discretization procedures are discussed with emphasis on automated generation of specific nonuniform solution grids and accounting of curved surfaces. The time-split algorithm was evaluated with regards to accuracy and convergence properties for hyperbolic equations on rectangular coordinates. An overall assessment of the viability of the finite element concept for computational aerodynamics is made.

Baker, A. J.↗
NASA NTRS1988

An adaptive pseudospectral method for discontinuous problems

The accuracy of adaptively chosen, mapped polynomial approximations is studied for functions with steep gradients or discontinuities. It is shown that, for steep gradient functions, one can obtain spectral accuracy in the original coordinate system by using polynomial approximations in a transformed coordinate system with substantially fewer collocation points than are necessary using polynomial expansion directly in the original, physical, coordinate system. It is also shown that one can avoid the usual Gibbs oscillation associated with steep gradient solutions of hyperbolic pde's by approximation in suitably chosen coordinate systems. Continuous, high gradient solutions are computed with spectral accuracy (as measured in the physical coordinate system). Discontinuous solutions associated with nonlinear hyperbolic equations can be accurately computed by using an artificial viscosity chosen to smooth out the solution in the mapped, computational domain. Thus, shocks can be effectively resolved on a scale that is subgrid to the resolution available with collocation only in the physical domain. Examples with Fourier and Chebyshev collocation are given.

Augenbaum, Jeffrey M.↗
NASA NTRS1989

An adaptive pseudospectral method for discontinuous problems

The accuracy of adaptively chosen, mapped polynomial approximations is studied for functions with steep gradients or discontinuities. It is shown that, for steep gradient functions, one can obtain spectral accuracy in the original coordinate system by using polynomial approximations in a transformed coordinate system with substantially fewer collocation points than are necessary using polynomial expansion directly in the original, physical, coordinate system. It is also shown that one can avoid the usual Gibbs oscillation associated with steep gradient solutions of hyperbolic pde's by approximation in suitably chosen coordinate systems. Continuous, high gradient solutions are computed with spectral accuracy (as measured in the physical coordinate system). Discontinuous solutions associated with nonlinear hyperbolic equations can be accurately computed by using an artificial viscosity chosen to smooth out the solution in the mapped, computational domain. Thus, shocks can be effectively resolved on a scale that is subgrid to the resolution available with collocation only in the physical domain. Examples with Fourier and Chebyshev collocation are given.

Augenbaum, J. M.↗