Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hyperbolic Systems”

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 433 records · Page 24

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

The proper modeling of nonequilibrium gas dynamics is required in certain regimes of hypersonic flow. For inviscid flow this gives a system of conservation laws coupled with source terms representing the chemistry. Often a wide range of time scales is present in the problem, leading to numerical difficulties as in stiff systems of ordinary differential equations. Stability can be achieved by using implicit methods, but other numerical difficulties are observed. The behavior of typical numerical methods on a simple advection equation with a parameter-dependent source term was studied. Two approaches to incorporate the source term were utilized: MacCormack type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. Various comparisons over a wide range of parameter values were made. In the stiff case where the solution contains discontinuities, incorrect numerical propagation speeds are observed with all of the methods considered. This phenomenon is studied and explained.

Leveque, R. J.↗

Colossal Optical Anisotropy from Atomic-Scale Modulations

Materials with large birefringence (Δn, where n is the refractive index) are sought after for polarization control (e.g., in wave plates, polarizing beam splitters, etc.), nonlinear optics, micromanipulation, and as a platform for unconventional light–matter coupling, such as hyperbolic phonon polaritons. Layered 2D materials can feature some of the largest optical anisotropy; however, their use in most optical systems is limited because their optical axis is out of the plane of the layers and the layers are weakly attached. This work demonstrates that a bulk crystal with subtle periodic modulations in its structure—Sr 9/8 TiS 3 —is transparent and positive-uniaxial, with extraordinary index n e = 4.5 and ordinary index n o = 2.4 in the mid- to far-infrared. The excess Sr, compared to stoichiometric SrTiS 3 , results in the formation of TiS 6 trigonal-prismatic units that break the chains of face-sharing TiS 6 octahedra in SrTiS 3 into periodic blocks of five TiS 6 octahedral units. The additional electrons introduced by the excess Sr form highly oriented electron clouds, which selectively boost the extraordinary index n e and result in record birefringence (Δn > 2.1 with low loss). The connection between subtle structural modulations and large changes in refractive index suggests new categories of anisotropic materials and also tunable optical materials with large refractive-index modulation.

36 MATERIALS SCIENCE↗

On the modeling and nonlinear dynamics of autonomous Silva-Young type chaotic oscillators with flat power spectrum

The nonlinear dynamics of fourth-order Silva-Young type chaotic oscillators with flat power spectrum recently introduced by Tamaseviciute and collaborators is considered. In this type of oscillators, a pair of semiconductor diodes in an anti-parallel connection acts as the nonlinear component necessary for generating chaotic oscillations. Based on the Shockley diode equation and an appropriate selection of the state variables, a smooth mathematical model (involving hyperbolic sine and cosine functions) is derived for a better description of both the regular and chaotic dynamics of the system. The complex behavior of the oscillator is characterized in terms of its parameters by using time series, bifurcation diagrams, Lyapunov exponents' plots, Poincaré sections, and frequency spectra. It is shown that the onset of chaos is achieved via the classical period-doubling and symmetry restoring crisis scenarios. Some PSPICE simulations of the nonlinear dynamics of the oscillator are presented in order to confirm the ability of the proposed mathematical model to accurately describe/predict both the regular and chaotic behaviors of the oscillator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of aerothermodynamic environment of a Titan aerocapture vehicle

The feasibility of an aerocapture vehicle mission has been emphasized recently for inner and outer planetary missions. Aerocapture involves a system concept which utilizes aerodynamic drag to acquire the velocity reduction necessary to obtain a closed planetary orbit from a hyperbolic flyby trajectory. It has been proposed to use the atmosphere of Titan for braking into a Saturn orbit. This approach for a Saturn orbital mission is expected to cut the interplanetary cruise travel time to Saturn from 8 to 3.5 years. In connection with the preparation of such a mission, it will be necessary to provide a complete analysis of the aerodynamic environment of the Titan aerocapture vehicle. The main objective of the present investigation is, therefore, to determine the extent of convective and radiative heating for the aerocapture vehicle under different entry conditions. This can be essentially accomplished by assessing the heating rates in the stagnation and windward regions of an equivalent body.

Tiwari, S. N.↗

Ulysses dust measurements near Jupiter

Submicrometer- to micrometer-sized particles were recorded by the Ulysses dust detector within 40 days of the Jupiter flyby. Nine impacts were recorded within 50 Jupiter radii with most of them recorded after closest approach. Three of these impacts are consistent with particles on prograde orbits around Jupiter and the rest are believed to have resulted from gravitationally focused interplanetary dust. From the ratio of the impact rate before the Jupiter flyby to the impact rate after the Jupiter flyby it is concluded that interplanetary dust particles at the distance of Jupiter move on mostly retrograde orbits. On 10 March 1992, Ulysses passed through an intense dust stream. The dust detector recorded 126 impacts within 26 hours. The stream particles were moving on highly inclined and apparently hyperbolic orbits with perihelion distances of greater than 5 astronomical units. Interplanetary dust is lost rather quickly from the solar system through collisions and other mechanisms and must be almost continuously replenished to maintain observed abundances. Dust flux measurements, therefore, give evidence of the recent rates of production from sources such as comets, asteroids, and moons, as well as the possible presence of interstellar grains.

Gruen, Eberhard↗

A One-year, Short-Stay Crewed Mars Mission Using Bimodal Nuclear Thermal Electric Propulsion (BNTEP) - A Preliminary Assessment

A crewed mission to Mars poses a significant challenge in dealing with the physiological issues that arise with the crew being exposed to a near zero-gravity environment as well as significant solar and galactic radiation for such a long duration. While long surface stay missions exceeding 500 days are the ultimate goal for human Mars exploration, short round trip, short surface stay missions could be an important intermediate step that would allow NASA to demonstrate technology as well as study the physiological effects on the crew. However, for a 1-year round trip mission, the outbound and inbound hyperbolic velocity at Earth and Mars can be very large resulting in a significant propellant requirement for a high thrust system like Nuclear Thermal Propulsion (NTP). Similarly, a low thrust Nuclear Electric Propulsion (NEP) system requires high electrical power levels (10 megawatts electric (MWe) or more), plus advanced power conversion technology to achieve the lower specific mass values needed for such a mission. A Bimodal Nuclear Thermal Electric Propulsion (BNTEP) system is examined here that uses three high thrust Bimodal Nuclear Thermal Rocket (BNTR) engines allowing short departure and capture maneuvers. The engines also generate electrical power that drives a low thrust Electric Propulsion (EP) system used for efficient interplanetary transit. This combined system can help reduce the total launch mass, system and operational requirements that would otherwise be required for equivalent NEP or Solar Electric Propulsion (SEP) mission. The BNTEP system is a hybrid propulsion concept where the BNTR reactors operate in two separate modes. During high-thrust mode operation, each BNTR provides 10's of kilo-Newtons of thrust at reasonably high specific impulse (Isp) of 900 seconds for impulsive transplanetary injection and orbital insertion maneuvers. When in power generation/EP mode, the BNTR reactors are coupled to a Brayton power conversion system allowing each reactor to generate 100's of kWe of electrical power to a very high Isp (3000 s) EP thruster system for sustained vehicle acceleration and deceleration in heliocentric space.

Burke, Laura M.↗

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

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme.

Yee, H. C.↗

Implicit total variation diminishing (TVD) schemes for steady-state calculations

The application of a new implicit unconditionally stable high resolution total variation diminishing (TVD) scheme to steady state calculations. It is a member of a one parameter family of explicit and implicit second order accurate schemes developed by Harten for the computation of weak solutions of hyperbolic conservation laws. This scheme is guaranteed not to generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments show that this scheme not only has a rapid convergence rate, but also generates a highly resolved approximation to the steady state solution. A detailed implementation of the implicit scheme for the one and two dimensional compressible inviscid equations of gas dynamics is presented. Some numerical computations of one and two dimensional fluid flows containing shocks demonstrate the efficiency and accuracy of this new scheme. Previously announced in STAR as N83-23085

Yee, H. C.↗

Interface procedures for overlapping grids

The solution of partial differential equations by means of either implicit or conservative explicit finite-difference methods, on the one hand, or multicomponent overlapping grid schemes on the other, is presently subjected to interpolation at grid boundaries in an attempt to eliminate the time lag. The accuracy of the transient solution of a hyperbolic or parabolic partial differential equation is dependent on the procedures used to transfer information between grids in a composite grid system. The present techniques reduce numerical solution errors.

Mastin, C. Wayne↗

Mobile Element Studies in Rocks (RAT) from Columbia Hills/West Spur at Gusev

Using elemental abundances determined by SPIRIT APX spectrometer on rocks and soils at Gusev Plains and Columbia Hills/ West Spur regions, the Athena Team discussed the aqueous geochemical implications at these sites on Mars. They suggested that these rocks were exposed to variable degrees of aqueous alteration (low to high) at Gusev crater. Earlier, we developed analytical procedures for studying aqueous geochemical behavior of fluids on rocks at Meridiani. In the present study, we apply these methods to rocks at Columbia Hills/West Spur in order to understand the significance of the Gusev rock results in reference to aqueous geochemical processes on Mars . The data analysis procedure is based on treating SO3 ("a") and Cl ("b") as two variables and tracking the relationship between "a" and "b" when the fluids undergo evaporation. This process of evaporation leads to concentration changes in these two elements finally producing salt assemblages on Martian rocks. In some cases on plotting "a"/ "b" versus "b" in salt assemblages, they yield a hyperbolic distribution. The relationship is transformed into a straight line when "a"/"b" is again plotted against 1/"b" in the system. Earlier, we used this procedure in the case of Merdiani rock abrasion tool (RAT) rocks and in this study, we discuss the application of this procedure to Gusev rocks. This study shows that the Gusev Plains rocks were exposed to low SO3/Cl solutions (sulfate-poor) for short period of time (weak interaction), whereas solutions with high SO3/Cl ratios (sulfate-rich) seem to have pervasively interacted with Columbia Hills/ West Spur rocks (strong interaction) at Gusev crater. Our conclusions seem to be consistent with the Mossbauer results given for these rocks

Rao, M. N.↗

Environmental damping and vibrational coupling of confined fluids within isolated carbon nanotubes

Abstract Because of their large surface areas, nanotubes and nanowires demonstrate exquisite mechanical coupling to their surroundings, promising advanced sensors and nanomechanical devices. However, this environmental sensitivity has resulted in several ambiguous observations of vibrational coupling across various experiments. Herein, we demonstrate a temperature-dependent Radial Breathing Mode (RBM) frequency in free-standing, electron-diffraction-assigned Double-Walled Carbon Nanotubes (DWNTs) that shows an unexpected and thermally reversible frequency downshift of 10 to 15%, for systems isolated in vacuum. An analysis based on a harmonic oscillator model assigns the distinctive frequency cusp, produced over 93 scans of 3 distinct DWNTs, along with the hyperbolic trajectory, to a reversible increase in damping from graphitic ribbons on the exterior surface. Strain-dependent coupling from self-tensioned, suspended DWNTs maintains the ratio of spring-to-damping frequencies, producing a stable saturation of RBM in the low-tension limit. In contrast, when the interior of DWNTs is subjected to a water-filling process, the RBM thermal trajectory is altered to that of a Langmuir isobar and elliptical trajectories, allowing measurement of the enthalpy of confined fluid phase change. These mechanisms and quantitative theory provide new insights into the environmental coupling of nanomechanical systems and the implications for devices and nanofluidic conduits.

36 MATERIALS SCIENCE↗

A study in three-dimensional chaotic dynamics: Granular flow and transport in a bi-axial spherical tumbler

We study three-dimensional (3D) chaotic dynamics through an analysis of transport in a granular flow in a half-full spherical tumbler rotated sequentially about two orthogonal axes (a bi-axial “blinking” tumbler). The flow is essentially quasi-two-dimensional in any vertical slice of the sphere during rotation about a single axis, and we provide an explicit exact solution to the model in this case. Hence, the cross-sectional flow can be represented by a twist map, allowing us to express the 3D flow as a linked twist map (LTM). We prove that if the rates of rotation about each axis are equal, then (in the absence of stochasticity) particle trajectories are restricted to two-dimensional (2D) surfaces consisting of a portion of a hemispherical shell closed by a “cap''; if the rotation rates are unequal, then particles can leave the surface they start on and traverse a volume of the tumbler. The period-one structures of the governing LTM are examined in detail: analytical expressions are provided for the location of period-one curves, their extent into the bulk of the granular material, and their dependence on the protocol parameters (rates and durations of rotations). Exploiting the restriction of trajectories to 2D surfaces in the case of equal rotation rates about the axes, a method is proposed for identifying and constructing 3D Kolmogorov--Arnold--Moser (KAM) tubes around the normally elliptic period-one curves. The invariant manifold structure arising from the normally hyperbolic period-one curves is also examined. When the motion is restricted to 2D surfaces, the structure of manifolds of the hyperbolic points in the bulk differs from that corresponding to hyperbolic points in the flowing layer. Each is reminiscent of a template provided by a non-integrable perturbation to a Hamiltonian system, though the governing LTM is not. This highlights the novel 3D chaotic behaviors observed in this model dynamical system.

42 ENGINEERING↗

An efficient full potential implicit method based on characteristics for analysis of supersonic flows

A nonlinear aerodynamic prediction technique based on the full potential equation in conservation form has been developed for the treatment of supersonic flows. The method uses the theory of characteristic signal propagation to accurately simulate the flow structure, which includes shock waves and mixed elliptic-hyperbolic crossflow. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. The necessary body-fitted grid system in every marching plane is generated numerically, using an elliptic grid solver. Results are shown for conical and nonconical wing-body combinations and compared with experimental data and Euler calculations. The method demonstrates an enormous savings in execution time and memory requirements over Euler methods.

Shankar, V.↗

Human Mars Mission: Launch Window from Earth Orbit

The determination of orbital window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to the earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a delta V penalty. Usually, because of the delta V penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: 1) One impulsive maneuver from a Highly Elliptical Orbit (HEO); 2) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO); 3) One impulsive maneuver from a Low Earth Orbit (LEO); 4) Two impulsive maneuvers form LEO; and 5) Three impulsive maneuvers form LEO. The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗

Human Mars Mission: Launch Window from Earth Orbit

The determination of orbital window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a DELTA V penalty. Usually, because of the DELTA V penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: (1) One impulsive maneuver from a Highly Elliptical Orbit (HEO) (2) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO) (3) One impulsive maneuver from a Low Earth Orbit (LEO) (4) Two impulsive maneuvers from LEO (5) Three impulsive maneuvers from LEO.

Young, Archie↗

Human Exploration Missions Study Launch Window from Earth Orbit

The determination of orbital launch window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a Delta(V) penalty. Usually, because of the Delta(V) penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: (1) One impulsive maneuver from a Low Earth Orbit (LEO), (2) Two impulsive maneuvers from LEO, (3) Three impulsive maneuvers from LEO, (4) One impulsive maneuvers from a Highly Elliptical Orbit (HEO), (5) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO) The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗

Stochastic and coarse-grained two-dimensional modeling of directional particle movement

In this work, we study the evolution of interacting groups of pedestrians in two-dimensional geometries. We introduce a microscopic stochastic model that includes floor fields modeling the global flow of individual groups as well as local interaction rules. From this microscopic model we derive an analytically-tractable system of conservation laws that governs the evolution of the macroscopic pedestrian densities. Numerical simulations show good agreement between the system of conservation laws and the microscopic model, though the latter is slightly more diffusive. We conclude by deriving second-order corrections to the system of conservation laws.

97 MATHEMATICS AND COMPUTING↗

Relaxation methods for unfactored implicit upwind schemes

Relaxation methods are presented for unfactored implicit upwind schemes for hyperbolic equations. The theoretical bases are explained using linear and nonlinear scalar equations; construction of the method for the unsteady Euler equations (nonlinear system) is but a natural extension. One of the important advantages of the above methods vis a vis factored implicit schemes is the possibility of faster convergence to steady state, as illustrated by the results. Several classes of relaxation schemes such as pointwise, linewise, Gauss-Seidel, and non-Gauss-Seidel methods are discussed, along with various strategies for convergence.

Chakravarthy, S. R.↗