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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

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

Insights on the continuous representations of piecewise-smooth nonlinear systems: limits of applicability and effectiveness

Dynamical systems subject to intermittent contact are often modeled with piecewise-smooth contact forces. However, the discontinuous nature of the contact can cause inaccuracies in numerical results or failure in numerical solvers. Representing the piecewise contact force with a continuous and smooth function can mitigate these problems, but not all continuous representations may be appropriate for this use. Here, five representations used by previous researchers (polynomial, rational polynomial, hyperbolic tangent, arctangent, and logarithm-arctangent functions) are studied to determine which ones most accurately capture nonlinear behaviors including super- and subharmonic resonances, multiple solutions, and chaos. The test case is a single-DOF forced Duffing oscillator with freeplay nonlinearity, solved using direct time integration. This work intends to expand on past studies by determining the limits of applicability for each representation and what numerical problems may occur.

42 ENGINEERING↗
NASA NTRS1976

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.↗
NASA NTRSDate not supplied

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↗
DOE OSTI2026

Quantum Ising model on (2+1)-dimensional anti–de Sitter space using tensor networks

We study the quantum Ising model on (2+1)-dimensional anti-de Sitter space using matrix product states (MPS) and matrix product operators (MPOs). We explore the bulk phase diagram of the theory on regular tessellations of hyperbolic space with coordination number seven and find disordered and ordered phases separated by a phase transition. We find that the boundary-boundary spin correlation function exhibits power law scaling deep in the disordered phase of the Ising model consistent with holography. At the critical point, we find the boundary entanglement entropy scales logarithmically with subsystem size but away from this, we see a linear scaling. In comparison, the full system exhibits a volume law scaling, which is expected in chaotic and/or highly connected systems. We also measure out of time ordered correlators (OTOCs) to explore the scrambling behavior of the theory.

Quantum spin models↗
NASA NTRS1975

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

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.↗
NASA NTRSDate not supplied

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↗
NASA NTRS1980

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

Self-Assembled BaTiO 3 –Au x Ag 1– x Low-Loss Hybrid Plasmonic Metamaterials with an Ordered “Nano-Domino-like” Microstructure

Metallic plasmonic hybrid nanostructures have attracted enormous research interest due to the combined physical properties coming from different material components and the broad range of applications in nanophotonic and electronic devices. However, the high loss and narrow range of property tunability of the metallic hybrid materials have limited their practical applications. In this study, a metallic alloy-based self-assembled plasmonic hybrid nanostructure, i.e., a BaTiO 3 –Au x Ag 1– x (BTO) vertically aligned nanocomposite, has been integrated by a templated growth method for low-loss plasmonic systems. Comprehensive microstructural characterizations including high-resolution scanning transmission electron microscopy (HRSTEM), energy-dispersive X-ray spectroscopy (EDS), and three-dimensional (3D) electron tomography demonstrate the formation of an ordered “nano-domino-like” morphology with Au 0.4 Ag 0.6 nanopillars as cylindrical cores and BTO as square shells. By comparing with the BTO–Au hybrid thin film, the BTO–Au 0.4 Ag 0.6 alloyed film exhibits much broader plasmon resonance, hyperbolic dispersion, low-loss, and thermally robust features in the UV–vis–NIR wavelength region. This study provides a feasible platform for a complex alloyed plasmonic hybrid material design with low-loss and highly tunable optical properties toward all-optical integrated devices.

36 MATERIALS SCIENCE↗
NASA NTRS1998

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

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

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

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

SNS: A Solution-Based Nonlinear Subspace Method for Time-Dependent Model Order Reduction

Several reduced order models have been successfully developed for nonlinear dynamical systems. To achieve a considerable speed-up, a hyper-reduction step is needed to reduce the computational complexity due to nonlinear terms. Many hyper-reduction techniques require the construction of nonlinear term basis, which introduces a computationally expensive offline phase. A novel way of constructing nonlinear term basis within the hyper-reduction process is introduced. In contrast to the traditional hyper-reduction techniques where the collection of nonlinear term snapshots is required, the SNS method avoids collecting the nonlinear term snapshots. Instead, it uses the solution snapshots that are used for building a solution basis, which enables avoiding an extra data compression of nonlinear term snapshots. As a result, the SNS method provides a more efficient offline strategy than the traditional model order reduction techniques, such as the DEIM, GNAT, and ST-GNAT methods. The SNS method is theoretically justified by the conforming subspace condition and the subspace inclusion relation. It is useful for model order reduction of large-scale nonlinear dynamical problems to reduce the offline cost. It is especially useful for ST-GNAT that has shown promising results, such as a good accuracy with a considerable online speed-up for hyperbolic problems in a recent paper by Choi and Carlberg [SIAM J. Sci. Comput., 41 (2019), pp. A26--A58], because ST-GNAT involves an expensive offline cost related to collecting nonlinear term snapshots. Error analysis for the SNS method is presented. Numerical results support that the accuracy of the solution from the SNS method is comparable to the traditional methods and a considerable speed-up (i.e., a factor of two to a hundred) is achieved in the offline phase.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A new effective theory for stochastic relativistic hydrodynamics

Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.

Mullins, Nicki [University of Illinois at Urbana-C↗
NASA NTRS2010

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

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↗