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 343 records · Page 19

A note on mass loss during collisions between galaxies and the formation of giant systems.

Stars will be liberated from the outer portions of galaxies which interact tidally. Using models based on the density distribution observed in the E1 galaxy NGC 3379, the mass loss is computed for hyperbolic collisions on an approximate but conservative basis. Binary collisions are then summed to show that large numbers of stars will be evaporated from elliptical members of dense clusters of galaxies. An estimate of one trillion solar masses is obtained for the amount of material dispersed in the Coma cluster, and in rich very dense systems (Morgan's Type I) the large amount of tidally dispersed debris may constitute the centrally located galaxies designated cD systems which often dominate these clusters.

Gallagher, J. S., III↗

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

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↗

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↗

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

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

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

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

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗

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↗

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↗