SEARCH · Engineering Papers

Results for “hyperbolic dynamics”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

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 NTRS1991

Fractal dimension in nonhyperbolic chaotic scattering

In chaotic scattering there is a Cantor set of input-variable values of zero Lebesgue measure (i.e., zero total length) on which the scattering function is singular. For cases where the dynamics leading to chaotic scattering is nonhyperbolic (e.g., there are Kolmogorov-Arnol'd-Moser tori), the nature of this singular set is fundamentally different from that in the hyperbolic case. In particular, for the nonhyperbolic case, although the singular set has zero total length, strong evidence is presented to show that its fractal dimension is 1.

Lau, Yun-Tung↗

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 NTRS1992

A numerical method for systems of conservation laws of mixed type admitting hyperbolic flux splitting

The present treatment of elliptic regions via hyperbolic flux-splitting and high order methods proposes a flux splitting in which the corresponding Jacobians have real and positive/negative eigenvalues. While resembling the flux splitting used in hyperbolic systems, the present generalization of such splitting to elliptic regions allows the handling of mixed-type systems in a unified and heuristically stable fashion. The van der Waals fluid-dynamics equation is used. Convergence with good resolution to weak solutions for various Riemann problems are observed.

Shu, Chi-Wang↗

NASA NTRS1993

On the global dynamics of adaptive systems - A study of an elementary example

The inherent nonlinear character of adaptive systems poses serious theoretical problems for the analysis of their dynamics. On the other hand, the importance of their dynamic behavior is directly related to the practical interest in predicting such undesirable phenomena as nonlinear oscillations, abrupt transients, intermittence or a high sensitivity with respect to initial conditions. A geometrical/qualitative description of the phase portrait of a discrete-time adaptive system with unmodeled disturbances is given. For this, the motions in the phase space are referred to normally hyperbolic (structurally stable) locally invariant sets. The study is complemented with a local stability analysis of the equilibrium point and periodic solutions. The critical character of adaptive systems under rather usual working conditions is discussed. Special emphasis is put on the causes leading to intermittence. A geometric interpretation of the effects of some commonly used palliatives to this problem is given. The 'dead-zone' approach is studied in more detail. The predicted dynamics are compared with simulation results.

Espana, Martin D.↗

NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

NASA NTRS1989

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗

NASA NTRS1991

Time domain convergence properties of Lyapunov stable penalty methods

Linear hyperbolic partial differential equations are analyzed using standard techniques to show that a sequence of solutions generated by the Liapunov stable penalty equations approaches the solution of the differential-algebraic equations governing the dynamics of multibody problems arising in linear vibrations. The analysis does not require that the system be conservative and does not impose any specific integration scheme. Variational statements are derived which bound the error in approximation by the norm of the constraint violation obtained in the approximate solutions.

Kurdila, A. J.↗

NASA NTRS1988

Compact wide-aperture hyperbolic analyzers

Compact electrostatic analyzers that employ a new geometry are described. They accept collimated beams of ions through relatively wide (non-slit-like) inlet apertures and focus them at various locations, depending on the ion energy per charge, along rectangular position-sensing detectors. Each analyzer has a dynamic range of 10 or more for a fixed deflection voltage. The focal line width at a given energy is a small fraction of the inlet aperture width, and the focal lines are approximately coplanar.

Curtis, C. C.↗

NASA NTRS1982

Boundary-fitted coordinate systems for numerical solution of partial differential equations - A review

A comprehensive review of methods of numerically generating curvilinear coordinate systems with coordinate lines coincident with all boundary segments is given. Some general mathematical framework and error analysis common to such coordinate systems is also included. The general categories of generating systems are those based on conformal mapping, orthogonal systems, nearly orthogonal systems, systems produced as the solution of elliptic and hyperbolic partial differential equations, and systems generated algebraically by interpolation among the boundaries. Also covered are the control of coordinate line spacing by functions embedded in the partial differential operators of the generating system and by subsequent stretching transformation. Dynamically adaptive coordinate systems, coupled with the physical solution, and time-dependent systems that follow moving boundaries are treated. References reporting experience using such coordinate systems are reviewed as well as those covering the system development.

Thompson, J. F.↗

NASA NTRS1975

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, it is assumed that the total temperature is constant. Therefore, the energy equation can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed. It is shown that the system is strictly hyperbolic and well posed for the initial value problems. Boundary conditions are described such that the linearized system is well posed. The Hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗

NASA NTRS2020

Optimal Low Thrust Orbit Transfers for Space Telescope Refueling at SEL2

The James Webb Space Telescope (JWST), a ten billion-dollar infrared telescope with a 6.5m primary mirror to be launched in 2021, is designed to operate in a Halo orbit around the second Sun-Earth Lagrange point (SEL2) for five to ten years. At that point fuel for station keeping and attitude maneuvers will run out. Refueling missions to JWST, as well as to similar space telescope missions proposed for SEL2, could greatly enhance the “science-per-dollar” value and promote a more sustainable use of space assets. In this paper, we present a novel approach to designing fuel optimal trajectories that will allow the refueling spacecraft to arrive at the SEL2 Halo orbit with maximum final mass (i.e. fuel payload). The low thrust optimal control problem is formulated using an indirect optimization method, leading to a two-point boundary value problem with a bang-bang control structure. We make use of a hyperbolic tangent smoothing technique for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and, thus, promote convergence. The problem is posed and solved in the circular restricted three-body problem. In this dynamical system, invariant manifolds exist that can be utilized to reduce fuel consumption. The here presented methodology to this challenging and important problem in astrodynamics demonstrates a significant potential for low-cost refueling mission design.

Woollands, Robyn↗

NASA NTRS2017

B-plane Evolution Under Highly Non-Keplerian Dynamics

Juno is the second of a series of New Frontiers missions and was launched in 2011. The spacecraft was set en route to Jupiter, and its Jupiter Orbit Insertion occurred on 2016-07-04 after a five-year cruise in deep space. The mission phase just prior to this large maneuver is called approach. As navigation team, we were responsible for precise orbit determination of the Juno spacecraft to ensure a successful orbit insertion. To evaluate the navigation performance we employed the B-plane mapping, where the incoming hyperbolic velocity defines the plane perpendicular to it, and the time to go to hit the plane will directly map into the uncertainties of the target point. The B-plane is a convenient tool when assessing the navigation performance as it is not affected by non-linearities due to the gravitational pull of the targeted body. However, as the Jupiter Orbit Insertion approached we started to realize that the B-plane mapping uncertainties were significantly larger than expected. We found that the B-plane mapping time had tremendous influence on the covariance inflation in the B-plane. Because of the close approach distance during the Jupiter Orbit Insertion and the strong interaction with J2 spherical harmonic coefficient, the Juno dynamics were too far from a Keplerian one assumed for the B-plane mapping. We first discuss the analytical approach of the B-plane mapping uncertainty and perform a numerical analysis to isolate the components of the Juno dynamics that cause the B-plane covariance inflation.

Bordi, John↗

NASA NTRS1990

A study of numerical methods for hyperbolic conservation laws with stiff source terms

In the present study of the behavior of typical numerical methods in the case of a model advection equation having a parameter-dependent source term, two approaches to the incorporation of the source terms are used: MacCormack-type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. The latter are found to perform slightly better. The model scalar equation is used to show that the incorrectness of the propagation speeds of discontinuities observed in the stiff case is due to the introduction of nonequilibrium values through numerical dissipation in the advection step.

Leveque, R. J.↗

NASA NTRS1976

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗

NASA NTRS1978

Extrapolation methods for dynamic partial differential equations

Several extrapolation procedures are presented for increasing the order of accuracy in time for evolutionary partial differential equations. These formulas are based on finite difference schemes in both the spatial and temporal directions. On practical grounds the methods are restricted to schemes that are fourth order in time and either second, fourth or sixth order in space. For hyperbolic problems the second order in space methods are not useful while the fourth order methods offer no advantage over the Kreiss-Oliger method unless very fine meshes are used. Advantages are first achieved using sixth order methods in space coupled with fourth order accuracy in time. Computational results are presented confirming the analytic discussions.

Turkel, E.↗

NASA NTRS2024

Capsule_Grid User Guide

Software utility CAPSULE_GRID and several related procedures that automate grid generation for computational fluid dynamics calculations about axisymmetric atmospheric entry vehicles are described as a supplement to the documentation within the source codes. This application program and its ancillary software utilities are all in the public domain at the web site https://software.nasa.gov/software/category/all/arc/1/cfdtools with the exception of the 2D hyperbolic volume gridding step for which scripts driving either GRIDGEN1 or HYPGEN2 are available. CAPSULE_GRID was developed at NASA Ames Research Center in support of the DPLR3 real gas flow solver. Both forebodies (only) and full bodies are treated, for axisymmetric 2D flow calculations or for 3D calculations at angle of attack. The expected application is to the smooth outer mold line of a hypersonic atmospheric entry vehicle, for which a single layer of point-matched structured grid blocks normally suffices. Control files for a variety of representative cases are included with illustrations of results, and some insights into the subtleties of curvature-based grid point distribution are provided. Gridding of a capsule on a sting is supported. Additional procedures for treating arbitrarily long 2D wakes, prompted by asteroid studies and relevant to spacecraft such as Stardust, are also documented. Most recently, treating an asymmetric forebody (including an off-center nose, with or without an aft body) has been incorporated as an extra step.

spacecraft↗

NASA NTRS1977

Fast methods incorporating direct elliptic solvers for nonlinear applications in fluid dynamics

Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Martin, E. D.↗