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 541 records · Page 30

Local gravity theories in conformal superspace

The conformal method for the initial value formulation of GR is one of the most powerful and widely usedtools in both analytical studies and numerical simulations of solutions. As is well-known, it exploits a hidden spatial conformal symmetry that can be extracted from the equations. But why conformal symmetry? Motivated by these questions and by well-known obstacles to quantum gravity, I look for the most general geometrodynamical symmetries compatible with a reduced physical configuration space for metric gravity. I argue that they are indeed spatial conformal diffeomorphisms. The next question is: what sort of theories most naturally incorporate this symmetry? Demanding locality for an action for metric gravity compatible with these principles determines the allowed operators, both for the purely gravitational part as well as matter couplings. The symmetries guarantee that there are two gravitational propagating physical degrees of freedom, but no explicit refoliation-invariance. The simplest such system has a geometric interpretation as a geodesic model in infinite-dimensional conformal superspace. One example of solution to the equations of motion corresponds to a static Bianchi IX spatial ansatz. The unique coupling to electromagnetism forces the electromagnetic equations to be hyperbolic, enabling us to “build” a standard space–time causal structure. There are, however, deviations from the standard Maxwell equations when space–time anisotropies become too large. Regarding quantization, with the geometric interpretation and the lack of refoliation invariance, the path integral treatment of the symmetries becomes much less involved than the similar approaches to GR. The symmetries form an (infinite-dimensional) Lie algebra, and no BFV treatment is necessary. Moreover, one can use gauge-invariant variational principles for selecting the boundary conditions of the path integral. We find that the propagator around the homogeneous solution has up to 6th order spatial derivatives, giving it plausible regularization properties as in Horava–Lifschitz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analytic representations of bath correlation functions for ohmic and superohmic spectral densities using simple poles

We present a scheme to express a bath correlation function (BCF) corresponding to a given spectral density (SD) as a sum of damped harmonic oscillations. Such a representation is needed, for example, in many open quantum system approaches. To this end we introduce a class of fit functions that enables us to model ohmic as well as superohmic behavior. We show that these functions allow for an analytic calculation of the BCF using pole expansions of the temperature dependent hyperbolic cotangent. We demonstrate how to use these functions to fit spectral densities exemplarily for cases encountered in the description of photosynthetic light harvesting complexes. Finally, we compare absorption spectra obtained for different fits with exact spectra and show that it is crucial to take properly into account the behavior at small frequencies when fitting a given SD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Variants and extensions of a fast direct numerical cauchy-riemann solver, with illustrative applications

Revised and extended versions of a fast, direct (noniterative) numerical Cauchy-Riemann solver are presented for solving finite difference approximations of first order systems of partial differential equations. Although the difference operators treated are linear and elliptic, one significant application of these extended direct Cauchy-Riemann solvers is in the fast, semidirect (iterative) solution of fluid dynamic problems governed by the nonlinear mixed elliptic-hyperbolic equations of transonic flow. Different versions of the algorithms are derived and the corresponding FORTRAN computer programs for a simple example problem are described and listed. The algorithms are demonstrated to be efficient and accurate.

Martin, E. D.↗

On the search for trajectories for a cometary nucleus sample return mission

The steps in the procedure for the generation of single trajectories for a comet nucleus sample return mission described by Stuhlinger et al. (1985) are discussed. The following assumptions are made for such a mission: (1) the spacecraft is equipped with an ion propulsion system supplied by solar cells, and (2) the comet has a low inclination and a period of about six years. The calculation of a specific trajectory incorporates the following main parameters: (1) the spacecraft mass at the beginning of the mission and the hyperbolic excess velocity, (2) the mass as a function of the injection energy, and (3) the available thrust level. The limitation of a supply of solar electric power to the region within 3 AU of the sun and the placement of permanent in-situ experiments at the comet are discussed. It is concluded that single trajectories for a nucleus sample return mission can be generated. A further step would be an integration of single solutions into classes of trajectories.

Schobert, D.↗

Optimal Low Thrust Orbit Transfers for Space Telescope Refueling at SEL2

The James Webb Space Telescope (JWST), a ten billion-dollar infrared telescope with a 6.5m primary mirror to be launched in 2021, is designed to operate in a Halo orbit around the second Sun-Earth Lagrange point (SEL2) for five to ten years. At that point fuel for station keeping and attitude maneuvers will run out. Refueling missions to JWST, as well as to similar space telescope missions proposed for SEL2, could greatly enhance the “science-per-dollar” value and promote a more sustainable use of space assets. In this paper, we present a novel approach to designing fuel optimal trajectories that will allow the refueling spacecraft to arrive at the SEL2 Halo orbit with maximum final mass (i.e. fuel payload). The low thrust optimal control problem is formulated using an indirect optimization method, leading to a two-point boundary value problem with a bang-bang control structure. We make use of a hyperbolic tangent smoothing technique for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and, thus, promote convergence. The problem is posed and solved in the circular restricted three-body problem. In this dynamical system, invariant manifolds exist that can be utilized to reduce fuel consumption. The here presented methodology to this challenging and important problem in astrodynamics demonstrates a significant potential for low-cost refueling mission design.

Woollands, Robyn↗

Impact of Magnetohydrodynamic (MHD) Induced Drag on High Ballistic Coefficient Aerocapture at Neptune

This paper presents analyses illustrating the benefits of magnetohydrodynamic (MHD) drag force control for a capsule performing aerocapture at Neptune relative to the same capsule utilizing conventional aerodynamic-only control across a range of design and flight conditions. MHD technology utilizes an onboard mag-net to create a magnetic field that interacts with the ionized plasma surrounding the capsule at hypersonic velocities, thus creating an electron current and inducing a Lorentz force used to increase capsule drag. The study assumes a 4.57-meterMars Science Laboratory (MSL) 70-degree sphere-cone capsule aerodynamic shape arrives hyperbolically at Neptune and aerocaptures into a synchronous Tri-ton orbit (3986 km x 430,000 km orbit inclined 157 degrees) that offers repeated flyby opportunities. Hyperbolic arrival velocity and the capsule entry mass (i.e., ballistic coefficient) are varied over a wide range of values to study the impact of these variables on key flight performance and sizing parameters. Results show reduced heat rates, heat loads, and thermal protection system mass for MHD cases at all conditions, with the benefit of MHD relative to conventional aerodynamic-only options increasing as ballistic coefficient increases, suggesting the payoff for MHD is significantly better for high-mass entry system design concepts with high ballistic coefficients.

Entry Descent and Landing↗

Unraveling the Mechanics: Computational Modeling of Residual Stresses in U10Mo Fuel with Aluminum Cladding

LEU (Low Enriched Uranium) plate-type fuel elements consist of a high-density, low-enrichment U–Mo alloy-based fuel foil encapsulated in an aluminum cladding and is fabricated through the Hot Isostatic Pressing (HIP) technique. Understanding the mechanism of possible failure modes of nuclear fuel is critical to mitigate potential consequences. One of the major contributing factors in various failure modes is stress and it greatly affects the integrity under temperature, pressure, and irradiation. These stresses could originate from various sources, including manufacturing, and operation, and can lead to deformation, cracking, or even complete failure of fuel elements. This study focuses on the development of a computational model that accurately predicts the residual stresses generated in the U10Mo fuel during fabrication. It has been observed that cladding creep has a substantial impact on the residual stresses in the U10Mo fuel post the HIP fabrication process. Furthermore, during the HIP bonding process the fuel plate system is heated to a temperature of 560 oC and as a result, the aluminum cladding transitions from (Al 6061-T6 to -O). This presents a challenge in capturing the change in material properties accurately. Therefore, a temperature dependent creep model such as hyperbolic sine creep model is considered to estimate the creep properties of the Aluminum cladding. The proposed calibrated creep model accurately predicts the residual stresses in the U10Mo fuel foil and agrees well with the experiments.

42 ENGINEERING↗

Powered Explicit Guidance Modifications and Enhancements for Space Launch System Block-1 and Block-1B Vehicles

NASA is currently building the Space Launch System (SLS) Block-1 launch vehicle for the Exploration Mission 1 (EM-1) test flight. NASA is also currently designing the next evolution of SLS, the Block-1B. The Block-1 and Block-1B vehicles will use the Powered Explicit Guidance (PEG) algorithm (of Space Shuttle heritage) for closed loop guidance. To accommodate vehicle capabilities and design for future evolutions of SLS, modifications were made to PEG for Block-1 to handle multi-phase burns, provide PEG updated propulsion information, and react to a core stage engine out. In addition, due to the relatively low thrust-to-weight ratio of the Exploration Upper Stage (EUS) and EUS carrying out Lunar Vicinity and Earth Escape missions, certain enhancements to the Block-1 PEG algorithm are needed to perform Block-1B missions to account for long burn arcs and target translunar and hyperbolic orbits. This paper describes the design and implementation of modifications to the Block-1 PEG algorithm as compared to Space Shuttle. Furthermore, this paper illustrates challenges posed by the Block-1B vehicle and the required PEG enhancements. These improvements make PEG capable for use on the SLS Block-1B vehicle as part of the Guidance, Navigation, and Control (GN&C) System.

Von der Porten, Paul↗

Framework for X-ray mirror surface shape fitting

For accurate characterization of grazing-incidence X-ray mirrors, we present a comprehensive framework to fit measured surface shapes (either slope or height) of X-ray mirrors used in synchrotron radiation and free-electron laser facilities. We summarize the closed-form expressions of some typical surface shapes of X-ray mirrors including elliptic cylinders, hyperbolic cylinders, ellipsoids, hyperboloids, and diaboloids. This framework is composed of four layers: definition of standard shapes with closed-form expressions, generation of theoretical surface with pose parameters (six degrees of freedom defining an object's position and orientation relative to a coordinate system), parameter optimization with the ability to select which parameters are fit and which are held constant, and the development of user-friendly fitting function wrappers for particular fitting tasks. A few practical fitting examples are demonstrated to verify the effectiveness of the proposed fitting framework. We discuss the physical meanings of the fitting parameters, and provide several examples using the elliptic cylinder and ellipsoid shapes to highlight some features of the framework. Moreover, we provide the presented framework as open-source codes (MATLAB and Python codes available at https://github.com/nsls2omf/xmf) to the community to encourage academic collaboration and further improvements.

36 MATERIALS SCIENCE↗

Dynamical evolution of the Oort cometary cloud

The dynamical evolution of comets in the Oort cloud under the influence of stellar perturbations has been modeled using Monte Carlo techniques. It is shown that the cloud has been depleted over the history of the solar system. Comets are lost from the cloud by direct ejection due to close stellar encounters, diffusion of aphelia to distances beyond the sun's sphere of influence, or diffusion of perihelia into the planetary region where Jupiter and Saturn perturbations either eject them on hyperbolic trajectories or capture them to short-period orbits. The population of the cloud is estimated to be 1.0 - 1.5 x 10 to the 12th comets and the total mass is on the order of 1.9 earth masses. In addition to random passing stars, less frequent encounters with giant molecular clouds may play a significant role in randomizing the orbits of comets in the cloud and reducing the effective radius of the sun's sphere of influence.

Weissman, P. R.↗

TIGGERC: Turbomachinery interactive grid generator energy distributor and restart code

A two dimensional multi-block grid generator was developed for a new design and analysis system for studying multi-blade row turbomachinery problems with an axisymmetric viscous/inviscid 'average passage' through flow code. TIGGERC is a mouse driven, fully interactive grid generation program which can be used to modify boundary coordinates and grid packing. TIGGERC generates grids using a hyperbolic tangent or algebraic distribution of grid points on the block boundaries and the interior points of each block grid are distributed using a transfinite interpolation approach. TIGGERC generates a blocked axisymmetric H grid, C grid, I grid, or O grid for studying turbomachinery flow problems. TIGGERC was developed for operation on small high speed graphic workstations.

Miller, David P.↗

Robust control of systems with real parameter uncertainty and unmodelled dynamics

Two significant contributions have been made during this research period in the research 'Robust Control of Systems with Real Parameter Uncertainty and Unmodelled Dynamics' under NASA Research Grant NAG-1-1102. They are: (1) a fast algorithm for computing the optimal H(sub infinity) norm for the four-block, the two block, or the one-block optimal H(sub infinity) optimization problem; and (2) a construction of an optimal H infinity controller without numerical difficulty. In using GD (Glover and Doyle) or DGKF (Doyle, Glover, Khargonekar, and Francis) approach to solve the standard H infinity norm which required bisection search. In this research period, we developed a very fast iterative algorithm for this computation. Our algorithm was developed based on hyperbolic interpolations which is much faster than any existing algorithm. The lower bound of the parameter, gamma, in the H infinity Riccati equation for solution existence is shown to be the square root of the supremum over all frequencies of the maximum eigenvalue of a given transfer matrix which can be computed easily. The lower band of gamma such that the H infinity Riccati equation has positive semidefinite solution can be also obtained by hyperbolic interpolation search. Another significant result in this research period is the elimination of the numerical difficulties arising in the construction of an optimal H infinity controller by directly applying the Glover and Doyle's state-space formulas. With the fast iterative algorithm for the computation of the optimal H infinity norm and the reliable construction of an optimal H infinity controller, we are ready to apply these tools in the design of robust controllers for the systems with unmodelled uncertainties. These tools will be also very useful when we consider systems with structured uncertainties.

Chang, Bor-Chin↗

Three dimensional flows around airfoils with shocks

The present work describes a mathematical model and numerical scheme for the computer-aided calculation of two- and three-dimensional transonic flow over an isolated yawed wing with oblique shock waves and a trailing vortex sheet. The flow is modeled by the potential equation for irrotational flow which is hyperbolic at supersonic points and elliptic at subsonic points. A coordinate-invariant difference scheme is used in which retarded difference formulas are constructed to conform with the local flow direction. The resulting 'rotated' difference scheme allows complete flexibility in the choice of a coordinate system. Shock waves are located automatically in the form of compression bands spread over a few mesh widths. The scheme has proven to be stable and convergent throughout the transonic range. Calculations have been performed for Mach number up to 1.2 and yaw angles up to 60 deg, the likely operating range of a yawed-wing transport designed to fly at supersonic speeds. Calculations becomes less accurate towards upper end of range, because the difference scheme is first-order accurate in the supersonic range.

Jameson, A.↗

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

Aerocapture Technology Development for Planetary Science - Update

Within NASA's Science Mission Directorate is a technological program dedicated to improving the cost, mass, and trip time of future scientific missions throughout the Solar System. The In-Space Propulsion Technology (ISPT) Program, established in 2001, is charged with advancing propulsion systems used in space from Technology Readiness Level (TRL) 3 to TRL6, and with planning activities leading to flight readiness. The program's content has changed considerably since inception, as the program has refocused its priorities. One of the technologies that has remained in the ISPT portfolio through these changes is Aerocapture. Aerocapture is the use of a planetary body's atmosphere to slow a vehicle from hyperbolic velocity to a low-energy orbit suitable for science. Prospective use of this technology has repeatedly shown huge mass savings for missions of interest in planetary exploration, at Titan, Neptune, Venus, and Mars. With launch vehicle costs rising, these savings could be the key to mission viability. This paper provides an update on the current state of the Aerocapture technology development effort, summarizes some recent key findings, and highlights hardware developments that are ready for application to Aerocapture vehicles and entry probes alike. Description of Investments: The Aerocapture technology area within the ISPT program has utilized the expertise around NASA to perform Phase A-level studies of future missions, to identify technology gaps that need to be filled to achieve flight readiness. A 2002 study of the Titan Explorer mission concept showed that the combination of Aerocapture and a Solar Electric Propulsion system could deliver a lander and orbiter to Titan in half the time and on a smaller, less expensive launch vehicle, compared to a mission using chemical propulsion for the interplanetary injection and orbit insertion. The study also identified no component technology breakthroughs necessary to implement Aerocapture on such a mission. Similar studies of Aerocapture applications at Neptune, Venus, and Mars were studied in 2003 through 2005. All showed significant performance improvements for the missions studied. Findings from these studies were used to guide the technology development tasks originally solicited in a 2002 NASA ROSS Research Announcement. The tasks are now in their final year and have provided numerous improvements in modeling and hardware, for use in proposals or new mission starts. Major Accomplishments: Since validation of the Aerocapture maneuver requires a space flight, ground developments have focused on modeling and environment prediction, materials, and sensors. Lockheed Martin has designed and built a 2-meter Carbon-Carbon aeroshell "hot structure." The article utilizes co-cured stiffening ribs and advanced insulation to achieve large scale, and up to a 40% reduction in areal density over the Genesis probe construction. This concept would be an efficient solution for probes that experience heat rates near 800-1000 W/cm(exp 2), such as at Venus and Earth. Applied Research Associates has extensively tested a family of efficient ablative TPS materials that provide solutions for a range of heating conditions. These materials are being applied to high-temperature structures built by ATK Space Systems, led by Langley Research Center. One-meter aeroshells will be thermally tested to validate construction and demonstrate higher bondline temperatures, which can lead to mass savings of up to 30% over traditional heatshields. Ames Research Center has developed aeroshell instrumentation that could measure environmental conditions and material performance during atmospheric entry. Instruments to measure TPS recession, heat flux, and catalycity could be combined with traditional sensors to provide a "plug-and-play" system for minimal mass and power, that would acquire flight data for model improvement and risk reduction on future missions. Improved atmospheric and aerothermodynamic models ha also been a major focus of the program. Next Steps: Aerocapture is one of five technologies in competition for a flight validation opportunity through the New Millennium Program. If selected, a fully autonomous vehicle will perform an Aerocapture at Earth in 2010, and flight data will be used to validate the guidance system and the TPS material for science mission infusion.

Munk, Michelle M.↗

Comment on 'Continuum modes in rotating plasmas: General equations and continuous spectra for large aspect ratio tokamaks'[Phys. Plasmas 18, 092103 (2011)]

It is shown that some of the main results of the recent paper by Lakhin and Ilgisonis [Phys. Plasmas 18, 092103 (2011)], viz. the derivation of the equations for the continuous spectra of poloidally and toroidally rotating plasmas and their special solution for large aspect ratio tokamaks with large parallel flows were obtained before by Goedbloed, Belieen, van der Holst, and Keppens [Phys. Plasmas 11, 28 (2004)]. A further rearrangement of the system of equations for the coupled Alfven and slow continuous spectra clearly exhibits: (a) coupling through a single tangential derivative, which is a generalization of the geodesic curvature; (b) the 'transonic' transitions of the equilibrium, which need to be carefully examined in order to avoid entering hyperbolic flow regimes where the stability formalism breaks down. A critical discussion is devoted to the implications of this failure, which is generally missed in the tokamak literature, possibly as a result of the wide-spread use of the sonic Mach number of gas dynamics, which is an irrelevant and misleading parameter in 'transonic' magnetohydrodynamics. Once this obstacle in understanding is removed, further application of the theory of trans-slow Alfven continuum instabilities to both tokamaks, with possible implications for the L-H transition, and astrophysical objects like 'fat' accretion disks, with a possible new route to magnetohydrodynamic turbulence, becomes feasible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multigrid Reduction in Time for Chaotic and Hyperbolic Problems (Final Report)

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallelin-time methods, with the most common method, parareal, diverging in many cases. The current state-of-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, this work shall continue researching the directions from the current collaboration with Dr. Falgout, which are laid out in the work Toward Parallel in Time for Chaotic Dynamical Systems and showed the first known results of a parallel-in-time speedup for a chaotic problem. This work outlines two key improvements to XBraid for chaotic problems, the so-called “theta” and “delta-correction” methods. Here, these two improvements will be implemented in a high-performance but general way in XBraid and explored for more complicated problems. We will additionally research, as time allows, improvements to these techniques, as well as multigrid relaxation techniques based on Least Squares Shadowing (LSS by Wang) and a nonintrusive block tridiagonal solver based on MGRIT, called TriMGRIT.

97 MATHEMATICS AND COMPUTING↗

Multigrid Reduction in Time for Chaotic and Hyperbolic Problems (Final Report)

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallel-in-time methods, with the most common method, parareal, diverging in many cases. The current state of-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, this work shall continue researching the directions from the current collaboration with Dr. Falgout, which are laid out in the work Toward Parallel in Time for Chaotic Dynamical Systems and showed the first known results of a parallel-in-time speedup for a chaotic problem. This work outlines two key improvements to XBraid for chaotic problems, the so-called “theta” and “delta-correction” methods. Here, these two improvements will be further researched and improved (including with a new relaxation method inspired by on Least Squares Shadowing (LSS)) and explored for more complicated problems.

97 MATHEMATICS AND COMPUTING↗