Engineering Papers⌕ Search

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.

At least 163 records · Page 9

iharm3D: Vectorized General Relativistic Magnetohydrodynamics

iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM (“High Accuracy Relativistic Magnetohydrodynamics”) algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most directly derived from Ryan et al. (2015) but with radiative transfer portions removed. HARM is a conservative finite-volume scheme for solving the equations of ideal GRMHD, a hyperbolic system of partial differential equations, on a logically Cartesian mesh in arbitrary coordinates.

79 ASTRONOMY AND ASTROPHYSICS↗

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytical investigation of the dynamics of tethered constellations in Earth orbit (phase 2)

Simulation of two short distance crawling maneuvers of the elevator both with and without environmental perturbations acting upon the system is discussed. These simulation runs were performed in order to provide results useful for the interpretation of the data from the tests, on the ground, of a scaled down engineering model of the elevator. In these simulation runs the elevator crawls along the tether in accordance to the developed mirror image motion control law (MIMCL). Results from the simulation of the 4 km long maneuver run were compared to those obtained by adopting the modified hyperbolic tangent control law (MHTCL). A preprocessor was developed for setting up the initial conditions of a tethered system with L platforms, M longitudinal dampers, and N lumped masses (platforms plus tether beads). A short test run of the 4-platform system with 3 longitudinal dampers and 10 lumper mass without any perturbation acting upon the system is illustrated. In support of the Tether Applications Working Group, SAO is preparing a catalog of tether simulations, has prepared a set of simulator test cases, obtained results from SKYHOOK, and solicited results from simulators at other institutions, and prepared a paper on a specific analytic solution.

Lorenzini, Enrico C.↗

A Robust Absorbing Boundary Condition for Compressible Flows

An absorbing non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented with theoretical proof. This paper is a continuation and improvement of a previous paper by the author. The absorbing NRBC technique is based on a first principle of non reflecting, which contains the essential physics that a plane wave solution of the Euler equations remains intact across the boundary. The technique is theoretically shown to work for a large class of finite volume approaches. When combined with the hyperbolic conservation laws, the NRBC is simple, robust and truly multi-dimensional; no additional implementation is needed except the prescribed physical boundary conditions. Several numerical examples in multi-dimensional spaces using two different finite volume schemes are illustrated to demonstrate its robustness in practical computations. Limitations and remedies of the technique are also discussed.

Loh, Ching Y.↗

Adiabatic behavior of an elliptical vortex in a time-dependent external strain flow

The adiabatic behavior of two-dimensional (2D) fluid vortices subject to slowly varying, time-dependent external strain flows is studied theoretically, and experimentally using magnetized non-neutral electron plasmas. Here the E × B drift dynamics of the electrons perpendicular to the magnetic field are analogous to the motion of fluid vorticity under the 2D Euler equations describing ideal fluids. A low-dimensional elliptical vortex patch model is used along with a WKB approximation to derive a formula for the breaking of an adiabatic invariant due to the changing external flow. The invariant is interpreted as the amplitude of a perturbative oscillation about a stable fixed point corresponding to a steady elliptical vortex. Smooth, hyperbolic tangent and piecewise linear ramp functions are considered for the external strain time dependence. Standard exponential breaking is observed in the former case, whereas the latter exhibits a power-law breaking curve with periodic modulation. It is found that a driving term in the equations of motion contributes the majority of the breaking whereas the frequency variation plays a weaker, but significant role. The experimental data agree closely with the theoretical model. The most significant deviation is due to inviscid damping behavior associated with the smooth edges of the experimental vortices, which tends to reduce the amplitude of the oscillation. Finally, the results are compared and contrasted with other related experimental, numerical, and theoretical work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamical evolution of the Oort cometary cloud

The dynamical evolution of comets in the Oort cloud under the influence of stellar perturbations has been modeled using Monte Carlo techniques. It is shown that the cloud has been depleted over the history of the solar system. Comets are lost from the cloud by direct ejection due to close stellar encounters, diffusion of aphelia to distances beyond the sun's sphere of influence, or diffusion of perihelia into the planetary region where Jupiter and Saturn perturbations either eject them on hyperbolic trajectories or capture them to short-period orbits. The population of the cloud is estimated to be 1.0 - 1.5 x 10 to the 12th comets and the total mass is on the order of 1.9 earth masses. In addition to random passing stars, less frequent encounters with giant molecular clouds may play a significant role in randomizing the orbits of comets in the cloud and reducing the effective radius of the sun's sphere of influence.

Weissman, P. R.↗

Advances in fluid dynamics

Various papers on advances in fluid dynamics are presented. Individual topics addressed include: nonunique solutions of the Euler equations; numerical solution of the biharmonic equation using a spectral multigrid method; finite-volume shape and placement for efficiency, accuracy, and physics; NND schemes and their application to numerical simulation of complex plume flow problems; symbolic analysis and computational algorithm construction; and velocity and mass flux distribution measurements of spherical glass beads in air flow in a 90-deg vertical-to-horizontal bend. Also discussed are: hyperbolic solution adaptive mesh generation; progress and future directions in CFD; numerical simulation of shock focusing; microbubble formation and splitting in a turbulent boundary layer for turbulence reduction; numerical study of spherical Taylor-Couette flow, trajectories of vortex lines beneath a free surface or above a plane; and a mathematical model for response of the coronary circulation to high sustained gravitational force fields.

Ballhaus, W. F., Jr.↗

On two-dimensional packets of capillary-gravity waves

The motion of a two-dimensional packet of capillary-gravity waves on water of finite depth is studied. The evolution of a packet is described by two partial differential equations: the nonlinear Schroedinger equation with a forcing term and a linear equation, which is of either elliptic or hyperbolic type depending on whether the group velocity of the capillary-gravity wave is less than or greater than the velocity of long gravity waves. These equations are used to examine the stability of the Stokes capillary-gravity wave train. The analysis reveals the existence of a resonant interaction between a capillary-gravity wave and a long gravity wave. The interaction requires that the liquid depth be small in comparison with the wavelength of the (long) gravity waves and the evolution equations describing the dynamics of this interaction are derived.

Djordjevic, V. D.↗

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↗

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗

Tidal disruption and the origin of the Moon

The dynamic problem of the tidal disruption of a rocky planetismal was solved by a direct integration of the fully three-dimensional, nonlinear equations of motion. The hypothesis that any object that passes within the Roche limit is disrupted was disproven. A time dependent solution was performed numerically, treating the planetismal as a fluid with a Murnaghan equation of state in the solid regions and zero pressure otherwise. Calculations show that a rocky body which passes by the Earth on a parabolic orbit with a perigee within the Roche limit is not tidally disrupted. Objects on hyperbolic orbits would experience even less tidal disruption. The results herein do not apply to bodies with very low viscosity. It is shown, however, that tidal disruption can be ruled out as a mechanism for reducing planetismal masses. Mechanisms for forming the Moon which rely upon tidal disruption are unlikely to be correct.

Mizuno, H.↗

Robust control of systems with real parameter uncertainty and unmodelled dynamics

Two significant contributions have been made during this research period in the research 'Robust Control of Systems with Real Parameter Uncertainty and Unmodelled Dynamics' under NASA Research Grant NAG-1-1102. They are: (1) a fast algorithm for computing the optimal H(sub infinity) norm for the four-block, the two block, or the one-block optimal H(sub infinity) optimization problem; and (2) a construction of an optimal H infinity controller without numerical difficulty. In using GD (Glover and Doyle) or DGKF (Doyle, Glover, Khargonekar, and Francis) approach to solve the standard H infinity norm which required bisection search. In this research period, we developed a very fast iterative algorithm for this computation. Our algorithm was developed based on hyperbolic interpolations which is much faster than any existing algorithm. The lower bound of the parameter, gamma, in the H infinity Riccati equation for solution existence is shown to be the square root of the supremum over all frequencies of the maximum eigenvalue of a given transfer matrix which can be computed easily. The lower band of gamma such that the H infinity Riccati equation has positive semidefinite solution can be also obtained by hyperbolic interpolation search. Another significant result in this research period is the elimination of the numerical difficulties arising in the construction of an optimal H infinity controller by directly applying the Glover and Doyle's state-space formulas. With the fast iterative algorithm for the computation of the optimal H infinity norm and the reliable construction of an optimal H infinity controller, we are ready to apply these tools in the design of robust controllers for the systems with unmodelled uncertainties. These tools will be also very useful when we consider systems with structured uncertainties.

Chang, Bor-Chin↗

Orbits Close to Asteroid 4769 Castalia

We use a radar-derived physical model of 4769 Castalia (1989 PB) to investigate close orbit dynamics around that kilometer- sized, uniformly rotating asteroid. Our methods of analysis provide a basis for systematic studies of particle dynamics close to any uniformly rotating asteroid. We establish that a Jacobi integral exists for particles orbiting this asteroid, examine the attendant zero-velocity surfaces, find families of periodic orbits, and determine their stability. All synchronous orbits and direct orbits within approx. 3 mean radii of Castalia are unstable and are subject to impact or escape from Castalia. Retrograde orbits are mostly stable and allow particles to orbit close to the asteroid surface. We derive a model which allows us to predict the escape conditions of a particle in orbit about Castalia and the (temporary) capture conditions for a hyperbolic interloper. Orbits within 1.5 km of Castalia are subject to immediate ejection from the system. Hyperbolic orbits with a V(sub infinity) less than 0.4 m/sec can potentially be captured by Castalia if their periapsis radius Is within approx. 2 km. For Castalia this capture region is small, but the results also apply to larger asteroids whose capture regions would also be larger. We determine bounds on ejecta speeds which either ensure ejecta escape or re-impact as functions of location on Castalia's surface. The speeds that ensure escape range from 0.28 to 0.84 m/sec and the speeds that ensure re-impact range from 0 to 0.18 m/sec. Speeds between these two bounds lead either to escape, re-impact, or potentially finite-time stable orbits. We develop a simple criterion which can establish whether a particle could have been ejected from the asteroid in the past and if it will impact the surface in the future.

Scheeres, D. J.↗

Nozomi Cis-Lunar Phase Orbit Determination

Japan's Institute of Space and Astronautical Science (ISAS) launched Nozomi, its first mission to the planet Mars using the newly developed M-V launch vehicle on July 3, 1998. Scientific objectives of the mission are to study the structure and dynamics of the Martian upper atmosphere and its interaction with the solar wind. Nozomi is a cooperative mission between ISAS and the National Aeronautics and Space Administration (NASA). The NASA contribution includes navigation and tracking services provided by the Jet Propulsion Laboratory (JPL). The spacecraft also serves as an engineering demonstration of basic technology for planetary exploration. One of the new technologies was a unique trajectory, developed by ISAS, which used solar gravitational perturbations at the weak stability boundary as an aid to achieve an Earth-Mars transfer orbit. This trajectory saves approximately 120 m/s of Delta V compared to direct hyperbolic insertion and is considered an enabling technology for the mission. Nozomi was the first spacecraft to employ this trajectory and provided on-orbit validation of the technique. The trajectory was achieved by initially placing the spacecraft in a highly elliptical cis-lunar phasing orbit. Six maneuvers were performed during this period to correct injection errors and target an outbound lunar swingby in September 1998. The gravity assist from the lunar swingby raised apogee to the vicinity of the weak stability boundary. After three more targeting maneuvers, Nozomi performed an inbound lunar swingby followed immediately by a powered Earth swingby in late December 1998. A 420 m/s Trans Mars Insertion (TMI) burn at the final Earth periapsis was intended to place the spacecraft on a heliocentric trajectory leading to Mars orbit insertion in October 1999. Orbit determination for Nozomi is performed in parallel by both ISAS and the Multi-Mission Navigation (MMNAV) group at JPL. This was an advantage for the mission because each group would generate solutions based on data collected from their respective tracking networks. Spacecraft events, such as sequence uplinks and maneuvers, were generally scheduled during passes at the Usuda tracking station in Japan. As a result, maneuver design and reconstruction was derived from MMNAV solutions based on JPL tracking data obtained immediately prior to or following maneuvers. Data was also exchanged between ISAS and MMNAV so orbit determination could be performed on joint data sets in support of critical targeting late in the cis-lunar phase. In this paper, information regarding the MMNAV orbit determination effort for the first six months of the mission is presented. The spacecraft trajectory is characterized first, followed by a discussion of the orbit determination estimation procedure and models. Results from selected orbit solutions are presented and compared against reconstructed trajectories. One area of emphasis in this paper is orbit determination in the vicinity of the weak stability boundary. Precise navigation was necessary to target the second lunar swingby and the powered Earth swingby. Delivery accuracy of 150 m was required for these critical encounters, but a number of factors contributed to the general degradation of orbit determination accuracy. This included the fact that the spacecraft was at apogee, at a range of 1.7 million km and moving at less than I km/sec perpendicular to the line of sight. Nozomi was also close to zero degrees declination where there are known limitations on orbit determination performance. Finally, S-band tracking data was acquired through the Nozomi backup low gain antenna. This antenna is offset from the axis of this spin stabilized spacecraft and superimposed large signatures in the Doppler and range data. These difficulties were overcome by combining long data arcs, spanning several maneuvers, with a high fidelity solar pressure model. The model included a physically accurate representation of the spacecraft structure and a high time resolution orientation model. Observation modeling included the removal of the spin induced Doppler bias, spin signature and per pass correction of range calibration errors applied for data leading up to critical events. As a result, all orbit determination goals were met. A second area of emphasis in this paper is the JPL tracking and orbit determination effort in support of the TMI maneuver. TMI occurred out of contact with ground stations and the JPL Goldstone tracking complex had the first pass following the bum. As a result, MMNAV had the responsibility to make a rapid assessment of the maneuver performance. MMNAV made the determination that a 100 m/s under bum had occurred and promptly informed ISAS via voice lines. ISAS immediately began preparations for a correction maneuver (TMIc), which had to be performed during the next Usuda pass. The near real time assessment by MMNAV provided accurate antenna frequency and pointing updates for the spacecraft acquisition at Usuda and the close coordination between the two agencies enabled the design and successful execution of the TMc maneuver. Propellant consumption during the correction burn dictated that the mission be redesigned. ISAS developed a new plan which adds 3 full solar orbits, two Earth swingbys and one lunar swingby with arrival at Mars in January 2004. The final Mars orbit will still enable the mission to achieve all of its science objectives.

Ryne, Mark↗

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↗

A spectral multidomain method for the solution of hyperbolic systems

A multidomain Chebyshev spectral collocation method for solving hyperbolic partial differential equations were developed. Though spectral methods are global methods, an attractive idea is to break a computational domain into several domains, and a way to handle the interfaces is described. The multidomain approach offers advantages over the use of a single Chebyshev grid. It allows complex geometries to be covered, and local refinement can be used to resolve important features. For steady state problems it reduces the stiffness associated with the use of explicit time integration as a relaxation scheme. Furthermore, the proposed method remains spectrally accurate. Results showing performance of the method on one dimensional linear models and one and two dimensional nonlinear gas dynamics problems are presented.

Kopriva, D.↗

Higher Order Parametric Excitation Modes for Spaceborne Quadrupole Mass Spectrometers

This paper describes a technique to significantly improve upon the mass peak shape and mass resolution of spaceborne quadrupole mass spectrometers (QMSs) through higher order auxiliary excitation of the quadrupole field. Using a novel multiresonant tank circuit, additional frequency components can be used to drive modulating voltages on the quadrupole rods in a practical manner, suitable for both improved commercial applications and spaceflight instruments. Auxiliary excitation at frequencies near twice that of the fundamental quadrupole RF frequency provides the advantages of previously studied parametric excitation techniques, but with the added benefit of increased sensed excitation amplitude dynamic range and the ability to operate voltage scan lines through the center of upper stability islands. Using a field programmable gate array, the amplitudes and frequencies of all QMS signals are digitally generated and managed, providing a robust and stable voltage control system. These techniques are experimentally verified through an interface with a commercial Pfeiffer QMG422 quadrupole rod system.When operating through the center of a stability island formed from higher order auxiliary excitation, approximately 50% and 400% improvements in 1% mass resolution and peak stability were measured, respectively, when compared with traditional QMS operation. Although tested with a circular rod system, the presented techniques have the potential to improve the performance of both circular and hyperbolic rod geometry QMS sensors.

Gershman, D. J.↗