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 379 records · Page 21

Far field boundary conditions for computation over long time

A new class of computational far field boundary conditions for hyperbolic partial differential equations is developed. These boundary conditions combine properties of absorbing boundary conditions for transient solutions and properties of far field boundary conditions for steady-state problems. The conditions can be used to limit the computational domain when both traveling waves and evanescent waves are present. Boundary conditions for scalar wave equations are derived and analyzed. Extensions to systems of equations are discussed and results from numerical experiments are presented.

Engquist, Bjorn↗

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗

A-Posteriori Error Estimation for Hyperbolic Conservation Laws with Constraint

This lecture considers a-posteriori error estimates for the numerical solution of conservation laws with time invariant constraints such as those arising in magnetohydrodynamics (MHD) and gravitational physics. Using standard duality arguments, a-posteriori error estimates for the discontinuous Galerkin finite element method are then presented for MHD with solenoidal constraint. From these estimates, a procedure for adaptive discretization is outlined. A taxonomy of Green's functions for the linearized MHD operator is given which characterizes the domain of dependence for pointwise errors. The extension to other constrained systems such as the Einstein equations of gravitational physics are then considered. Finally, future directions and open problems are discussed.

Barth, Timothy↗

Geometry and properties of the plasma in a quadrupole ion pump

The diameter of the plasma in a quadrupole ionization gauge and quadrupole ion pump has been measured as a function of pressure. The radio-frequency quadrupole systems are tuned in such a way that electrons remain on stable trajectories in the center of the tube. Being emitted at one end, the electrons encounter an inverted field at the other end so that they travel on spiral paths back and forth until they collide with the rest gas molecules, producing ions which are collected as an ion current at a cylinder surrounding the four hyperbolically shaped quadrupole electrodes.

Schwarz, H.↗

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic, and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching step sizes with the implicit shock fitting procedure. A correlation is found between the marching step size, Reynolds number, and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching-step sizes with the implicit shock fitting procedure. A correlation is found between the marching-step size, Reynolds number and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗

The evolution of comet orbits

The origin of comets and the evolution of their orbits are discussed. Factors considered include: the law of survival of comets against ejection on hyperbolic orbits; short-period comets are not created by single close encounters of near-parabolic comets with Jupiter; observable long-period comets do not evolve into observable short-period comets; unobservable long-period comets with perihelia near Jupiter can evolve into observable short-period comets; long-period comets cannot have been formed or created within the planetary region of the solar system (excluding the effects of stellar perturbations); it is possible that some of the short-period comets could have been formed inside the orbit of Neptune; circularly-restricted three-body problem, and its associated Jacobi integral, are not valid approximations to use in studying origin and evolution of comets.

Everhart, E.↗

Thermal Protection System to Enable Ice Giant Aerocapture Mission for Delivering Both an Orbiter and an In Situ Probe

The Ice Giants have been identified as high priority science destinations in the last Decadal Survey [1] and could benefit from aerocapture as the primary method for orbit insertion [2]. A mass-efficient aerocapture system will enable the delivery of an orbiter along with an atmospheric probe (for in situ measurements to anchor global data collected by the orbiter) and possibly a lander at Triton [3]. Aerocapture could be executed either using low L/D rigid aeroshell with lift modulation (LMA) [4] or using deployable aeroshell using drag modulation (DMA) [3]. Nearly two decades ago, a NASA-funded team performed Neptune-Triton aerocapture studies with a mid-L/D lifting configuration [5] for achieving orbit using LMA. This study showed aerocapture challenges. Due to very high peak entry conditions combined with very high heat-load, a suite of TPS materials was required and this suite was deemed problematic from a qualification perspective, due test facility limitations. In the past 20 years, progress made in GN&C for lift-guided entry missions such as MSL, Orion EFT1, Mars 2020 and the upcoming Artemis missions, and the expertise in blunt body aerodynamics at large scale (~ 5m) has led the EDL community to conclude that aerocapture is a “go do” engineering activity and most technologies are in hand to propose missions with aerocapture [6] [7]. Aerocapture using DMA, currently in development, is an option for Ice Giant Missions. While DMA is simpler in some sense, due to ballistic entry and no need for lift-guided maneuvering, it has challenges and it’s maturity is lower. LMA and DMA both require one or more ablative Thermal Protection System (TPS) materials for the rigid aeroshell element. The ablative TPS needs to be robust and mass efficient due to the high heat loads and size of the rigid aeroshell. Currently, there are capable ablative thermal protection materials, e.g., Heatshield for Extreme Entry Environments Technology (HEEET), 3-D woven Mid-Density Carbon- Phenolic (3MDCP), and PICA (Phenolic-Impregnated Carbon Ablator) that are mature, i.e., at TRL 6 or higher. NASA also invested in Conformal PICA that was matured to TRL 5. Our goal is to evaluate the applicability of high TRL TPS and consider other design options. We first establish bounding aerocapture trajectories for a wide range of arrival conditions and the associated aerothermal environment. Based on the environments, we then determine the predicted TPS mass for the aeroshell [4]. In this presentation, we will outline the process by which we establish bounding aerocapture trajectories for hyperbolic excess velocities ranging from 27 km/s to 35 km/s, which are shown to be a range of velocities that can reduce the trip time from ~14 years to 8 years. The above velocity range translates to ~12 km/s to ~24 km/s arrival velocities at the planetary entry interface [2]. The velocity reduction required to achieve orbit ranges between ~2.5 km/s to 9.5 km/s for both Neptune and Triton. Propulsive insertion alone, due to the amount of fuel required to achieve the required velocity reduction, limits the science returned [2]. We establish the bounding aerocapture trajectories for a low L/D (~ 0.4) configuration for three different ballistic coefficients. The ballistic coefficient range is determined from three different aeroshell diameters of 3m, 4m and 5m and with an entry system mass of 2200 kg. With the above range of design parameters, we then determine conservative/bounding estimates of aerothermal environments by using a combination of CFD simulations and stagnation point heating estimates [7]. This engineering approach allows us to first assess the TPS need vs. TPS capability and determine the applicability of existing TPS. Once an applicable suite of TPS is determined, the TPS thickness and mass are computed. We show that the TPS mass fraction can be as low as 5% to as high as 20%, depending on total trip time reduction and other design parameters for a range of TPS. This is a large range for TPS mass fraction. We show PICA and HEEET can indeed enable aerocapture missions, but the missions incur a mass penalty. TPS mass savings, can be further reduced with the use of conformal PICA. Advancing the development of Conformal PICA to make it robust across the entire aerothermal environment (peak heat-flux, pressure and shear) range will result in TPS mass fractions of < 10% for Ice Giant aerocapture missions such as the Neptune-Triton mission. Aerocapture allows for not only shortening the trip time but enables larger mass to be placed in orbit. Furthermore, probes deployed from orbit will benefit in reduced entry environments allowing for a lower risk TPS implementation as compared to mission designs where the probe is released prior to orbit insertion. One of the challenges for the Ice Giant community is to ensure mission designs that maximize science and allow flexibility in the placement of the entry probe. The traditional approach to release the probe ahead of the orbiter may not optimize returned science. In this presentation, we will make the case for mature TPS such as HEEET and PICA. While these materials can enable aerocapture missions, completing the development of conformal PICA and extending Conformal PICA to be more robust, will have significant impact to TPS mass efficiency and significantly enhance science return for future Gas- and Ice-Giant missions.

E Venkatapathy↗

Diagonalization and simultaneous symmetrization of the gas-dynamic matrices

The hyperbolic nature of the unsteady, inviscid, gas-dynamic equations implies the existence of a similarity transformation for diagonalizing an arbitrary linear combination of coefficient matrices. It is shown that the individual matrices are simultaneously symmetrized by the similarity transformation. The transformations and their norms can be applied to the well-posedness of the Cauchy problem, linear stability theory for finite-difference approximations, and simplification of block-tridiagonal systems that arise in implicit time-split algorithms.

Warming, R. F.↗

Multiple-Beam Communications Antenna

System has both fixed-spot and scanning coverage. Advanced offset-fed spacecraft antenna systems operating in 30/20 GHz frequency bands developed to provide multiple-radiating fixed-spot and regional-coverage scanning beams for use on communications satellites. Operating at higher frequencies, provides alternate frequency bands for expansion of existing satellite services and achieving frequency reuse capability for conservation of frequency spectrum. 20-GHz transmitting antenna is offset-fed, dual-reflector configuration. 30-GHz receiving antenna has orthogonal polarized feeds, hyperbolic reflectors, and two subreflectors.

Myhre, R. W.↗

Economical Third-Order Methods for Accurate Surface Heating Predictions on Simplex Element Meshes

A node-centered edge-based finite volume discretization of the compressible Navier-Stokes equations is presented with the heat flux reformulated as a first order system. A dissipation vector is derived for the reformulated system, such that the heat flux can be upgraded to O(h^3) on simplex element meshes in the same fashion as the inviscid fluxes. The method of manufactured solutions is used to demonstrate this uniform order property in heat flux. This new system is shown to efficiently produce accurate surface heating predictions on hypersonic hemisphere flow using an anisotropic simplex element mesh, achieving O(h^3) accuracy at relatively low computational cost compared to similar methods.

Computational Fluid Dynamics↗

Hyperbolic/parabolic development for the GIM-STAR code

Flow fields in supersonic inlet configurations were computed using the eliptic GIM code on the STAR computer. Spillage flow under the lower cowl was calculated to be 33% of the incoming stream. The shock/boundary layer interaction on the upper propulsive surface was computed including separation. All shocks produced by the flow system were captured. Linearized block implicit (LBI) schemes were examined to determine their application to the GIM code. Pure explicit methods have stability limitations and fully implicit schemes are inherently inefficient; however, LBI schemes show promise as an effective compromise. A quasiparabolic version of the GIM code was developed using elastical parabolized Navier-Stokes methods combined with quasitime relaxation. This scheme is referred to as quasiparabolic although it applies equally well to hyperbolic supersonic inviscid flows. Second order windward differences are used in the marching coordinate and either explicit or linear block implicit time relaxation can be incorporated.

Spradley, L. W.↗

A Parallel Non-Overlapping Domain-Decomposition Algorithm for Compressible Fluid Flow Problems on Triangulated Domains

This paper considers an algebraic preconditioning algorithm for hyperbolic-elliptic fluid flow problems. The algorithm is based on a parallel non-overlapping Schur complement domain-decomposition technique for triangulated domains. In the Schur complement technique, the triangulation is first partitioned into a number of non-overlapping subdomains and interfaces. This suggests a reordering of triangulation vertices which separates subdomain and interface solution unknowns. The reordering induces a natural 2 x 2 block partitioning of the discretization matrix. Exact LU factorization of this block system yields a Schur complement matrix which couples subdomains and the interface together. The remaining sections of this paper present a family of approximate techniques for both constructing and applying the Schur complement as a domain-decomposition preconditioner. The approximate Schur complement serves as an algebraic coarse space operator, thus avoiding the known difficulties associated with the direct formation of a coarse space discretization. In developing Schur complement approximations, particular attention has been given to improving sequential and parallel efficiency of implementations without significantly degrading the quality of the preconditioner. A computer code based on these developments has been tested on the IBM SP2 using MPI message passing protocol. A number of 2-D calculations are presented for both scalar advection-diffusion equations as well as the Euler equations governing compressible fluid flow to demonstrate performance of the preconditioning algorithm.

Barth, Timothy J.↗

Total Variation Diminishing (TVD) schemes of uniform accuracy

Explicit second-order accurate finite-difference schemes for the approximation of hyperbolic conservation laws are presented. These schemes are nonlinear even for the constant coefficient case. They are based on first-order upwind schemes. Their accuracy is enhanced by locally replacing the first-order one-sided differences with either second-order one-sided differences or central differences or a blend thereof. The appropriate local difference stencils are selected such that they give TVD schemes of uniform second-order accuracy in the scalar, or linear systems, case. Like conventional TVD schemes, the new schemes avoid a Gibbs phenomenon at discontinuities of the solution, but they do not switch back to first-order accuracy, in the sense of truncation error, at extrema of the solution. The performance of the new schemes is demonstrated in several numerical tests.

Hartwich, PETER-M.↗

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗

Upwind and symmetric shock-capturing schemes

The development of numerical methods for hyperbolic conservation laws has been a rapidly growing area for the last ten years. Many of the fundamental concepts and state-of-the-art developments can only be found in meeting proceedings or internal reports. This review paper attempts to give an overview and a unified formulation of a class of shock-capturing methods. Special emphasis is on the construction of the basic nonlinear scalar second-order schemes and the methods of extending these nonlinear scalar schemes to nonlinear systems via the extact Riemann solver, approximate Riemann solvers, and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows is discussed. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas dynamics problems.

Yee, H. C.↗

Figure and Dimension Metrology of Extremely Lightweight X-Ray Mirrors for Space Astronomy Applications

The International X-ray Observatory (IXO) is the next major space X-ray observatory, performing both imaging and spectroscopic studies of all kinds of objects in the Universe. It is a collaborative mission of the National Aeronautics and Space Administration of the United States, the European Space Agency, and Japan Aerospace Exploration Agency. It is to be launched into a Sun-Earth L2 orbit in 2021. One of the most challenging aspects of the mission is the construction of a flight mirror assembly capable focusing X-rays in the band of 0.1 to 40 keY with an angular resolution of better than 5 arc-seconds and with an effective collection area of more than 3 sq m. The mirror assembly will consist of approximately 15,000 parabolic and hyperbolic mirror segments, each of which is approximately 200mm by 300mm with a thickness of 0.4mm. The manufacture and qualification of these mirror segments and their integration into the giant mirror assembly have been the objectives of a vigorous technology development program at NASA's Goddard Space Flight Center. Each of these mirror segments needs to be measured and qualified for both optical figure and mechanical dimensions. In this talk, I will describe the technology program with a particular emphasis on a measurement system we are developing to meet those requirements, including the use of coordinate measuring machines, Fizeau interferometers, and custom-designed, and -built null lens. This system is capable of measuring highly off-axis aspherical or cylindrical mirrors with repeatability, accuracy, and speed.

Zhang, William W.↗

VAPRE: A tool to support the design optimization of planetary entry probe missions

The importance of the atmospheric data collected by planetary entry probes has been emphasized since the results of Galileo in 1995 and are reflected in the interest . The results add to our fundamental understanding about the formation and the evolution of our solar system. Of special interest are the Outer Planets Saturn, Uranus and Neptune, serving as missing puzzle pieces to the existing data sets. In order to maximize the data and science return from those missions, the tool VAPRE is developed to facilitate the design optimization of planetary entry probe missions by increasing the trade-space. Using a data set of interplanetary trajectories, it assesses the change in availability and accessibility of entry sites on the planets with respect to their latitudes depending on the hyperbolic arrival velocity. In this paper, we introduce the motivation behind the implementation of VAPRE, the science and engineering input gathered to bridge the knowledge between both fields, the design structure of the implementation as well as the significance and future development and potential of the tool VAPRE for future mission concept developments for planetary probe and lander missions.

Simon, Amy↗