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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fast Multigrid Reduction-in-Time for Advection via Modified Semi-Lagrangian Coarse-Grid Operators

Many iterative parallel-in-time algorithms have been shown to be highly efficient for diffusion-dominated partial differential equations (PDEs) but are inefficient or even divergent when applied to advection-dominated PDEs. We consider the application of the multigrid reduction-in-time (MGRIT) algorithm to linear advection PDEs. Here, the key to efficient time integration with this method is using a coarse-grid operator that provides a sufficiently accurate approximation to the so-called ideal coarse-grid operator. For certain classes of semi-Lagrangian discretizations, we present a novel semi-Lagrangian-based coarse-grid operator that leads to fast and scalable multilevel time integration of linear advection PDEs. The coarse-grid operator is composed of a semi-Lagrangian discretization followed by a correction term, with the correction designed so that the leading-order truncation error of the composite operator is approximately equal to that of the ideal coarse-grid operator. Parallel results show substantial speed-ups over sequential time integration for variable-wave-speed advection problems in one and two spatial dimensions, and using high-order discretizations up to order five. The proposed approach establishes the first practical method that provides small and scalable MGRIT iteration counts for advection problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A compressible two-phase model for dispersed particle flows with application from dense to dilute regimes

Multiphase flows are present in many important fields ranging from multiphase explosions to chemical processing. An important subset of multiphase flow applications involves dispersed materials, such as particles, droplets, and bubbles. This work presents an Eulerian–Eulerian model for multiphase flows containing dispersed particles surrounded by a continuous media such as air or water. Following a large body of multiphase literature, the driving force for particle acceleration is modeled as a direct function of both the continuous-phase pressure gradient and the gradient of intergranular stress existing within the particle phase. While the application of these two components of driving force is well accepted in much of the literature, other models exist in which the particle-phase pressure gradient itself drives particle motion. The multiphase model treats all phases as compressible and is derived to ensure adherence to the 2nd Law of Thermodynamics. The governing equations are presented and discussed, and a characteristic analysis shows the model to be hyperbolic, with a degeneracy in the case that the intergranular stress, which is modeled as a configuration pressure, is zero. Finally, results from a two sample problems involving shock-induced particle dispersion are presented. The results agree well with experimental measurements, providing initial confidence in the proposed model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Analysis of a second-order-accurate finite-volume method for temporally-growing compressible shear layers

A finite-volume method for solving the compressible laminar Navier-Stokes equation in two dimensions is assessed for use in the simulation of transitional flows. The method is second-order-accurate, uses alternating-direction-implicit time integration, and a total-variation-diminishing smoothing operator in the inviscid flux terms. The test problem was a free shear layer with a forced periodicity in x at a specified wavelength and with the parallel mean flow specified as a hyperbolic-tangent profile. A grid refinement investigation verified this method to be second-order-accurate with less than 4 percent error in the growth rate of isolated modes on a 32 by 64 grid. The interaction between modes of similar amplitude was small, less than 1 percent. However, in cases where there was a dominant mode, numerical phase error pumped energy into all other modes. This error results in nonphysical growth rate for modes whose energy content is several orders of magnitude below the dominant mode. In the case when the dominant mode saturates, the pumping action stops and the growth rate of the next largest growing mode regains physical significance.

Atkins, H. L.↗
DOE OSTI2024

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1971

Optimum transfer between hyperbolic asymptotes with turn angle less than the maximum natural turn angle.

The impulsive transfer of minimum characteristic velocity between two given hyperbolic asymptotes associated with a real planet is considered. The constraint of remaining above the surface of the planet introduces a maximum natural turn angle. This divides the problem naively into two categories. The first is the case when the turn angle required by the given asymptotes is less than or equal to the maximum natural turn angle. This is the more reasonable situation for a swingby maneuver, and is the case considered in this paper. The resulting optimal transfers are of five basic types. In boundary conditions space, by far the largest area of transfers are either single impulse transfers or transfers through the parabolic level requiring two finite impulses. In no case has an optimal transfer been found which uses more than two finite impulses and three infinitesimal impulses. Complete results are presented for the case in which the given turn angle is less than or equal to the maximum natural turn angle.

Walton, J. M.↗
NASA NTRS1978

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗
NASA NTRS1999

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗
DOE OSTI2024

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Wilson loops and minimal area surfaces in hyperbolic space

The AdS/CFT correspondence relates Wilson loops in N = 4 SYM theory to minimal area surfaces in AdS space. If the loop is a plane curve the minimal surface lives in hyperbolic space H 3 (or equivalently Euclidean AdS 3 space). We argue that finding the area of such extremal surface can be easily done if we solve the following problem: given two real periodic functions V 0,1 (s), V 0,1 (s +2π) = V 0,1 (s) a third periodic function V 2 (s) is to be found such that all solutions to the equation &#x2212;<!-- - --> &#x2202;<!-- ∂ --> S 2 &#x03D5;<!-- φ --> + [ V 0 + 1 2 ( &#x03BB;<!-- λ --> + 1 &#x03BB;<!-- λ --> ) V 1 + i 2 ( &#x03BB;<!-- λ --> &#x2212;<!-- - --> 1 &#x03BB;<!-- λ --> ) V 2 ] &#x03D5;<!-- φ --> = 0 are anti-periodic in s ϵ [0, 2π] for any value of λ. This problem is equivalent to the statement that the monodromy matrix is trivial. It can be restated as that of finding a one complex parameter family of curves X(λ, s) where X(λ = 1, s) is the given shape of the Wilson loop and such that the Schwarzian derivative {X(λ, s), s} is meromorphic in λ with only two simple poles. We present a formula for the area in terms of the functions V 0,1,2 and discuss solutions to these equivalent problems in terms of theta functions. Finally, we also consider the near circular Wilson loop clarifying its integrability properties and rederiving its area using the methods described in this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

Computation of space shuttle flow fields using noncentered finite-difference schemes.

Second- and third-order, noncentered finite-difference schemes are described for the numerical solution of the hyperbolic equations of fluid dynamics. The advantages of noncentered methods over the more conventional centered schemes are: simpler programming logic, nonhomogeneous terms are easily included, and generalization to multidimensional problems is direct. Second- and third-order methods are compared with regard to dissipative and dispersive errors and shock-capturing ability. These schemes are then used in a shock-capturing technique to determine the inviscid, supersonic flow field surrounding space shuttle vehicles (SSV). Resulting flow fields about typical pointed and blunted, delta-winged SSVs at angle of attack are presented and compared with experiment.

Kutler, P.↗
NASA NTRS1980

Technology requirements for a generic aerocapture system

The technology requirements for the design of a generic aerocapture vehicle system are summarized. These spacecraft have the capability of completely eliminating fuel-costly retropropulsion for planetary orbit capture through a single aerodynamically controlled atmospheric braking pass from a hyperbolic trajectory into a near circular orbit. This generic system has application at both the inner and outer planets. Spacecraft design integration, navigation, communications, and aerothermal protection system design problems were assessed in the technology requirements study and are discussed in this paper.

Cruz, M. I.↗
NASA NTRS1985

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The novel implicit and unconditionally stable, high resolution Total Variation Diminishing (TVD) scheme whose application to steady state calculations is presently examined is a member of a one-parameter family of implicit, second-order accurate systems developed by Harten (1983) for the computation of weak solutions for one-dimensional hyperbolic conservation laws. The scheme will not generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments for a quasi-one-dimensional nozzle problem show that the experimentally determined stability limit correlates exactly with the theoretical stability limit for the nonlinear scalar hyberbolic conservation laws.

Yee, H. C.↗
NASA NTRS1988

Far field boundary conditions for computation over long time

A new class of computational far field boundary conditions for hyperbolic partial differential equations is developed. These boundary conditions combine properties of absorbing boundary conditions for transient solutions and properties of far field boundary conditions for steady-state problems. The conditions can be used to limit the computational domain when both traveling waves and evanescent waves are present. Boundary conditions for scalar wave equations are derived and analyzed. Extensions to systems of equations are discussed and results from numerical experiments are presented.

Engquist, Bjorn↗
DOE OSTI2019

Strong Stability Preserving Integrating Factor Two-Step Runge-Kutta Methods

In this report, problems with components that feature significantly different time scales, where the stiff time-step restriction comes from a linear component, implicit-explicit (IMEX) methods alleviate this restriction if the concern is linear stability. However, when nonlinear non-inner-product stability properties are of interest, such as in the evolution of hyperbolic partial differential equations with shocks or sharp gradients, linear inner-product stability is no longer sufficient for convergence, and so strong stability preserving (SSP) methods are often needed. Where the SSP property is needed, IMEX SSP Runge-Kutta (SSP-IMEX) methods have very restrictive time-steps. An alternative to SSP-IMEX schemes is to adopt an integrating factor approach to handle the linear component exactly and step the transformed problem forward using some time-evolution method. The strong stability properties of integrating factor Runge-Kutta methods were established in Isherwood et al. (SIAM J Numer Anal 56(6):3276–3307, 2018), where it was shown that it is possible to define explicit integrating factor Runge-Kutta methods that preserve strong stability properties satisfied by each of the two components when coupled with forward Euler time-stepping. It was proved that the solution will be SSP if the transformed problem is stepped forward with an explicit SSP Runge-Kutta method that has non-decreasing abscissas. However, explicit SSP Runge-Kutta methods have an order barrier of p=4, and sometimes higher order is desired. In this work we consider explicit SSP two-step Runge-Kutta integrating factor methods to raise the order. We show that strong stability is ensured if the two-step Runge-Kutta method used to evolve the transformed problem is SSP and has non-decreasing abscissas. We find such methods up to eighth order and present their SSP coefficients. Adding a step allows us to break the fourth order barrier on explicit SSP Runge-Kutta methods; furthermore, our explicit SSP two-step Runge-Kutta methods with non-decreasing abscissas typically have larger SSP coefficients than the corresponding one-step methods. A selection of our methods are tested for convergence and demonstrate the design order. We also show, for selected methods, that the SSP time-step predicted by the theory is a lower bound of the allowable time-step for linear and nonlinear problems that satisfy the total variation diminishing (TVD) condition. We compare some of the non-decreasing abscissa SSP two-step Runge-Kutta methods to previously found methods that do not satisfy this criterion on linear and nonlinear TVD test cases to show that this non-decreasing abscissa condition is indeed necessary in practice as well as theory. We also compare these results to our SSP integrating factor Runge-Kutta methods designed in Isherwood et al. (2018).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

An approximation theory for the identification of linear thermoelastic systems

An abstract approximation framework and convergence theory for the identification of thermoelastic systems is developed. Starting from an abstract operator formulation consisting of a coupled second order hyperbolic equation of elasticity and first order parabolic equation for heat conduction, well-posedness is established using linear semigroup theory in Hilbert space, and a class of parameter estimation problems is then defined involving mild solutions. The approximation framework is based upon generic Galerkin approximation of the mild solutions, and convergence of solutions of the resulting sequence of approximating finite dimensional parameter identification problems to a solution of the original infinite dimensional inverse problem is established using approximation results for operator semigroups. An example involving the basic equations of one dimensional linear thermoelasticity and a linear spline based scheme are discussed. Numerical results indicate how the approach might be used in a study of damping mechanisms in flexible structures.

Rosen, I. G.↗
NASA NTRS1998

A Formulation of Asymptotic and Exact Boundary Conditions Using Local Operators

In this paper we describe a systematic approach for constructing asymptotic boundary conditions for isotropic wave-like equations using local operators. The conditions take a recursive form with increasing order of accuracy. In three dimensions the recursion terminates and the resulting conditions are exact for solutions which are described by finite combinations of angular spherical harmonics. First, we develop the expansion for the two-dimensional wave equation and construct a sequence of easily implementable boundary conditions. We show that in three dimensions and analogous conditions are again easily implementable in addition to being exact. Also, we provide extensions of these ideas to hyperbolic systems. Namely, Maxwell's equations for TM waves are used to demonstrate the construction. Finally, we provide numerical examples to demonstrate the effectiveness of these conditions for a model problem governed by the wave equation.

Hagstrom, T.↗
NASA NTRS1994

Stability and error estimation for Component Adaptive Grid methods

Component adaptive grid (CAG) methods for solving hyperbolic partial differential equations (PDE's) are discussed in this paper. Applying recent stability results for a class of numerical methods on uniform grids. The convergence of these methods for linear problems on component adaptive grids is established here. Furthermore, the computational error can be estimated on CAG's using the stability results. Using these estimates, the error can be controlled on CAG's. Thus, the solution can be computed efficiently on CAG's within a given error tolerance. Computational results for time dependent linear problems in one and two space dimensions are presented.

Oliger, Joseph↗