SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1978

On the accuracy and convergence of implicit numerical integration of finite element generated ordinary differential equations

A study of accuracy and convergence of linear functional finite element solution to linear parabolic and hyperbolic partial differential equations is presented. A variable-implicit integration procedure is employed for the resultant system of ordinary differential equations. Accuracy and convergence is compared for the consistent and two lumped assembly procedures for the identified initial-value matrix structure. Truncation error estimation is accomplished using Richardson extrapolation.

Baker, A. J.↗

NASA NTRS1995

The large discretization step method for time-dependent partial differential equations

A new method for the acceleration of linear and nonlinear time dependent calculations is presented. It is based on the Large Discretization Step (LDS) approximation, defined in this work, which employs an extended system of low accuracy schemes to approximate a high accuracy discrete approximation to a time dependent differential operator. Error bounds on such approximations are derived. These approximations are efficiently implemented in the LDS methods for linear and nonlinear hyperbolic equations, presented here. In these algorithms the high and low accuracy schemes are interpreted as the same discretization of a time dependent operator on fine and coarse grids, respectively. Thus, a system of correction terms and corresponding equations are derived and solved on the coarse grid to yield the fine grid accuracy. These terms are initialized by visiting the fine grid once in many coarse grid time steps. The resulting methods are very general, simple to implement and may be used to accelerate many existing time marching schemes.

Haras, Zigo↗

DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1986

Regularity estimates up to the boundary for elliptic systems of difference equations

Regularity estimates up to the boundary for solutions of elliptic systems of finite difference equations were proved. The regularity estimates, obtained for boundary fitted coordinate systems on domains with smooth boundary, involve discrete Sobolev norms and are proved using pseudo-difference operators to treat systems with variable coefficients. The elliptic systems of difference equations and the boundary conditions which are considered are very general in form. The regularity of a regular elliptic system of difference equations was proved equivalent to the nonexistence of eigensolutions. The regularity estimates obtained are analogous to those in the theory of elliptic systems of partial differential equations, and to the results of Gustafsson, Kreiss, and Sundstrom (1972) and others for hyperbolic difference equations.

Strikwerda, J. C.↗

NASA NTRS1973

The VISC code: A user's manual

The VISC code is a computer automated scheme for solving the equations describing the fully coupled viscous, radiating flow at the stagnation-point of a blunt body which may or may not be ablating. The code provides a basis for obtaining prediction of the stagnation-point heating to a body entering any planetary atmosphere at hyperbolic velocities. The code is written in FORTRAN V and is operational on both the Univac 1108 (EXEC 8) system and the CDC 7600 system. The report gives an overview of the VISC code computational logic flow, a description of the input requirements and output results and comments on the practical use of the code. As such the report forms a users manual for operation of the VISC code.

Wilson, K.↗

NASA NTRS1973

SL4 code - A user's manual

The SL-4 code is a computer automated scheme for solving the equations describing the fully-coupled viscous, radiating flow over the front face of a blunt body which may or may not be ablating. The code provides a basis for obtaining predictions of the surface beating to a body entering any planetary atmosphere at hyperbolic velocities. The code is written in FORTRAN V and is operational on both the Univac 1108 (EXEC 8) system in use at LMSC and the CDC 7600 system in use at the University of California, Berkeley. An overview of the SL-4 code computational logic flow, a description of the input requirements and output results, and comments on the practical use of the code are presented. As such this report forms a users manual for operation of the SL-4 code.

Chou, Y. S.↗

NASA NTRS1979

Aerocapture - A system design for planetary exploration

The paper presents the results of a study that developed conceptual designs of an atmospheric entry system with moderate to high L/D ratios to establish the concepts of aerocapture into low circular Mars orbits from hyperbolic flyby trajectories. The payloads considered in this study were those currently conceived for the Mars Sample Return (MSR) mission including single mission (SM) and multi-mission (MM) modes. The design concepts developed are also applicable (with small modification) for establishment of closed orbits about Venus and Saturn. The MSR mission payload requires an orbiter, lander, ascent vehicle, and rover, with a total mass of approximately 4000 kg. The emphasis during the study was placed on development of concepts for aerodynamic configuration, the aeroshell design, and a guidance, navigation, and control subsystem.

Cruz, M. I.↗

NASA NTRS1988

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.↗

NASA NTRS1982

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.↗

NASA NTRS1992

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↗

NASA NTRS2013

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.↗

NASA NTRS1983

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.↗

NASA NTRS1983

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

Yee, H. C.↗

NASA NTRS1986

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↗

NASA NTRS2007

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.↗

DOE OSTI2024

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↗

NASA NTRS1982

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.↗

NASA NTRS1999

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↗