Perturbations of a Spacecraft Orbit During an Hyperbolic Flyby
A first-order perturbation theory exists for the flyby problem with the Born approximation as the reference orbit.
perturbations Hyperbolic flyby Celestrial Mechanic↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Geometric control of hyperbolic exciton-polariton condensate dimers
Coupled many-body quantum systems exhibit rich emergent physics with diverse stationary and dynamical behaviours. By engineering platforms with tunable and distinct coupling mechanisms, new insights emerge into the collective behaviour of coupled many body systems. Particles can be exchanged via evanescent or ballistic coupling: the former, based on proximity, yields large spectral splitting, while the latter requires strict phase-matching, analogous to phase-coupled harmonic oscillators and has a smaller impact on the energy landscape. We demonstrate an all-optically tunable quantum fluid dimer based on exciton-polariton condensates in a photonic crystal waveguide with hyperbolic (saddle-like) dispersion. Varying the dimer’s angle relative to the grating tunes the coupling from evanescent to ballistic. We directly observe spectral features and mass flow shaped by the saddle dispersion. This work highlights photonic crystals as powerful platforms to explore condensed matter phenomena lying at the interface between delay-coupled nonlinear oscillators and tight binding physics.
Georgakilas, Ioannis [IBM Research-Zurich, Rüschli↗
An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case
Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.
Goodrich, John W.↗
High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs
This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.
Nielsen, Tanner B.↗
High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes
The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.
Barth, Timothy↗
Mixed gauge-global symmetries, elliptic modes, and black hole thermodynamics in Hořava-Lifshitz gravity
In Hořava-Lifshitz gravity, a putative consistent theory of quantum gravity for which there is evidence for both black hole thermodynamics and a holographic construction, spacetime is endowed with a preferred dynamical spacelike foliation. The theory has a leaf reparameterization symmetry that is neither global nor local gauge, hyperbolic and elliptic equations of motion, a lack of splittability, and universal horizon black hole solutions. The reparameterization symmetry is “mixed”: it is a local symmetry in one coordinate yet global on each leaf. More broadly it is an example of both unfree and projectable gauge symmetries. The mixed symmetry and associated charge has not yet been accounted for in calculations of universal horizon thermodynamics in Hořava-Lifshitz gravity. This has led to problems, in particular the failure of the first law in a class of asymptotically AdS solutions where the normal to the leaves of the foliation is not aligned with the time translation Killing vector at infinity. We show how the dynamics of the charge corresponding to this symmetry coupled with the other features above resolves this issue. We then briefly comment how this mixed symmetry, the corresponding charge, and the elliptic equations of motion also conspire to evade recent holographic arguments for only local gauge fields in consistent theories of quantum gravity due to the lack of splittability of the elliptic equation and associated mode.
Global Symmetries↗
Development and application of unified algorithms for problems in computational science
A framework is presented for developing computationally unified numerical algorithms for solving nonlinear equations that arise in modeling various problems in mathematical physics. The concept of computational unification is an attempt to encompass efficient solution procedures for computing various nonlinear phenomena that may occur in a given problem. For example, in Computational Fluid Dynamics (CFD), a unified algorithm will be one that allows for solutions to subsonic (elliptic), transonic (mixed elliptic-hyperbolic), and supersonic (hyperbolic) flows for both steady and unsteady problems. The objectives are: development of superior unified algorithms emphasizing accuracy and efficiency aspects; development of codes based on selected algorithms leading to validation; application of mature codes to realistic problems; and extension/application of CFD-based algorithms to problems in other areas of mathematical physics. The ultimate objective is to achieve integration of multidisciplinary technologies to enhance synergism in the design process through computational simulation. Specific unified algorithms for a hierarchy of gas dynamics equations and their applications to two other areas: electromagnetic scattering, and laser-materials interaction accounting for melting.
Shankar, Vijaya↗
A Trajectory-Driven Algorithm for Differentiating SRB Measures on Unstable Manifolds
Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.
97 MATHEMATICS AND COMPUTING↗
Mathematical modelling of the concave front in the adjacent high explosive detonation problem
This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.
acceleration↗
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.↗
A probability of encounter with interstellar comets and the likelihood of their existence
A theory of the probability of encounter of the sun with an interstellar comet at a distance comparable to the earth-sun distance is formulated, and a general expression is derived establishing the relationship among the influx rate of interstellar comets, the perihelion distance, the space density of the comets, the Maxwellian distribution of comet velocities in the interstellar cloud, and the cloud's systematic velocity relative to the sun. An upper limit of .0006 solar mass per cubic parsec is established for the space density of interstellar comets. The theoretical distribution of semimajor axes of interstellar comets shows that a strong hyperbolic excess must be present in the orbits of a majority of interstellar comets regardless of the dynamical characteristics of the comet cloud, except when the cloud is moving along with the sun and the distribution of individual velocities has a very low dispersion. This case, however, implies a possibility of capture by the sun and thus becomes a problem of an Oort-type cloud.
Sekanina, Z.↗
Pressure Distribution Measurements at High Speed and Oblique Incidence of Flow
The present report contains the results of a series of observations obtained for a wing of symmetrical profile for different angles of yaw. The shock tunnel with 0.4m x 0.4m cross section, of the Institute for Gas Dynamics at L.F.A. - Braunchweig was available for this work. The profile used was a 9 percent, thick hyperbolic profile with maximum thickness 40 percent aft, which was calculated by the method of F. Ringleb (2). Since the conformal transformation of this profile is known, the theoretical pressure distribution could be determined exactly for the case of incompressible flow. Then by the use of the Prandtl rule one may extend the comparison of theory and experiment to the case of higher velocity of incident flow.
Lippisch, A.↗
Looking at extremal black holes from very far away
Near-extremal black holes are subject to large quantum effects, which modify their low-temperature thermodynamic behavior. Hitherto, these quantum effects were analyzed by separating the geometry into the near-horizon region and its exterior. It is desirable to understand and reproduce such corrections from the full higher-dimensional asymptotically flat or AdS geometry’s perspective. We address this question in this article and fill this gap. Specifically, we find off-shell eigenmodes of the quadratic fluctuation operator of the Euclidean gravitational dynamics, with eigenvalues that vanish linearly with temperature. We illustrate this for BTZ and neutral black holes with hyperbolic horizons in AdS in Einstein-Hilbert theory, and for the charged black holes in Einstein-Maxwell theory. The linear scaling with Matsubara frequency, which is a distinctive feature of the modes, together with the fact that their wavefunctions localize close to the horizon as we approach extremality, identifies them as responsible for the aforementioned quantum effects. We provide a contour prescription to deal with the sign indefiniteness of the Euclidean Einstein-Maxwell action, which we derive to aid our analysis. We also resolve a technical puzzle regarding modes associated with rotational isometries in stationary black hole spacetimes.
AdS-CFT Correspondence↗
Where Are the Interstellar Comets?
The existing theories of comet dynamics held that large numbers of comets are ejected to interstellar space and that other stellar systems would do the same; however, no comet on a hyperbolic orbit has ever been observed. Calculations developed by Duncan et al. are used to show that most comets thought to be ejected would actually be captured by the Oort Cloud.
comets interstellar comets Oort Cloud orbital dyna↗
An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD
The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Dynamical capture of field stars by interstellar clouds
Evidence suggesting the existence of old stars in young galactic clusters, has motivated an investigation of a model in which such stars are captured from the field during the gas cloud collapse phase of star formation. A hyperbolic field star cannot be gravitationally captured by a static cloud since the star's total energy is conserved. However, if the cloud's potential is an explicit function of time, the energy of a penetrating star can be decreased. As a first step in modeling the dynamical capture of field stars, capture by freely falling, homogeneous spherical cloud is considered. This simplified model is solved, and a prediction is obtained for number of captured field stars in terms of the cloud parameters and dispersion velocities of the field stars.
Whitman, P. G.↗
Numerical methods for systems of conservation laws of mixed type using flux splitting
The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.
Shu, Chi-Wang↗
A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods
The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts 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 waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic 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 flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.
Yee, H. C.↗