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

Observation of Universal Quench Dynamics and Townes Soliton Formation from Modulational Instability in Two-Dimensional Bose Gases

We experimentally study universal nonequilibrium dynamics of two-dimensional atomic Bose gases quenched from repulsive to attractive interactions. Here, we observe the manifestation of modulational instability that, instead of causing collapse, fragments a large two-dimensional superfluid into multiple wave packets universally around a threshold atom number necessary for the formation of Townes solitons. We confirm that the density distributions of quench-induced solitary waves are in excellent agreement with the stationary Townes profiles. Furthermore, our density measurements in the space and time domain reveal detailed information about this dynamical process, from the hyperbolic growth of density waves, the formation of solitons, to the subsequent collision and collapse dynamics, demonstrating multiple universal behaviors in an attractive many-body system in association with the formation of a quasistationary state.

74 ATOMIC AND MOLECULAR PHYSICS↗

Non-intrusive nonlinear model reduction via machine learning approximations to low-dimensional operators

Abstract Although projection-based reduced-order models (ROMs) for parameterized nonlinear dynamical systems have demonstrated exciting results across a range of applications, their broad adoption has been limited by their intrusivity: implementing such a reduced-order model typically requires significant modifications to the underlying simulation code. To address this, we propose a method that enables traditionally intrusive reduced-order models to be accurately approximated in a non-intrusive manner. Specifically, the approach approximates the low-dimensional operators associated with projection-based reduced-order models (ROMs) using modern machine-learning regression techniques. The only requirement of the simulation code is the ability to export the velocity given the state and parameters; this functionality is used to train the approximated low-dimensional operators. In addition to enabling nonintrusivity, we demonstrate that the approach also leads to very low computational complexity, achieving up to $$10^3{\times }$$ 10 3 × in run time. We demonstrate the effectiveness of the proposed technique on two types of PDEs. The domain of applications include both parabolic and hyperbolic PDEs, regardless of the dimension of full-order models (FOMs).

42 ENGINEERING↗

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimum solar electric interplanetary mission opportunities from 1975 to 1990

A collection of optimum trajectory and spacecraft data is presented for unmanned interplanetary missions from 1975 to 1990 using solar electric propulsion. Data are presented for one-way flyby and orbiter missions from Earth to Venus, Mars, Jupiter, Saturn, Uranus, Neptune, and Pluto. The solar system model assumes planetary ephemerides which very closely approximate the true motion of the planets. Direct and indirect flight profiles are investigated. Data are presented for two representative flight times for each mission. The launch vehicle is the Titan 3 B (core)/Centaur, and a constant jet exhaust speed solar electric propulsion system having a specific mass of 30 kg/kw is completely optimized in terms of power level and jet exhaust speed to yield maximum net spacecraft mass. The hyperbolic excess speeds at departure and arrival and the launch date are optimized for each mission. For orbiter missions, a chemical retro stage is used to brake the spacecraft into a highly eccentric capture orbit about the target planet.

Mann, F. I.↗

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗

A simple hyperbolic model for communication in parallel processing environments

We introduce a model for communication costs in parallel processing environments called the 'hyperbolic model,' which generalizes two-parameter dedicated-link models in an analytically simple way. Dedicated interprocessor links parameterized by a latency and a transfer rate that are independent of load are assumed by many existing communication models; such models are unrealistic for workstation networks. The communication system is modeled as a directed communication graph in which terminal nodes represent the application processes that initiate the sending and receiving of the information and in which internal nodes, called communication blocks (CBs), reflect the layered structure of the underlying communication architecture. The direction of graph edges specifies the flow of the information carried through messages. Each CB is characterized by a two-parameter hyperbolic function of the message size that represents the service time needed for processing the message. The parameters are evaluated in the limits of very large and very small messages. Rules are given for reducing a communication graph consisting of many to an equivalent two-parameter form, while maintaining an approximation for the service time that is exact in both large and small limits. The model is validated on a dedicated Ethernet network of workstations by experiments with communication subprograms arising in scientific applications, for which a tight fit of the model predictions with actual measurements of the communication and synchronization time between end processes is demonstrated. The model is then used to evaluate the performance of two simple parallel scientific applications from partial differential equations: domain decomposition and time-parallel multigrid. In an appropriate limit, we also show the compatibility of the hyperbolic model with the recently proposed LogP model.

Stoica, Ion↗

Aerocapture: An Enabling Technology for Flagship-Class Uranus Orbiter and Probe Mission

Introduction: The current decadal survey published by the National Academies of Sciences has informed National Aeronautics and Space Administration (NASA) to prioritize the study of the Ice Giants, especially Uranus. To gather the required data that addresses the science questions raised in this survey, a mission to Uranus with an orbiter and atmospheric probe must be designed. The Uranus Orbiter and Probe (UOP) study, which the survey identified as the flagship mission of this decade, proposes a 2031 launch to take advantage of a Jupiter fly-by and utilizes a fully propulsive orbit insertion design with an Earth-to-Uranus transit times ranging from 13 to 15 years. This fully propulsive orbit insertion at Uranus will be very fuel expensive (wet mass percentages of around 60-70\%) thereby leaving less mass for the scientific payload and additional planetary probes. In addition, scientists are more interested in visiting Uranus before 2049, when the Spring Equinox will occur, as it allows studying Uranus seasons not seen during Voyager 2's flyby in 1987. A NASA Flagship-Class mission would require at least 10 years of lead time prior to launch thereby making the 2031 launch to take advantage of the Jupiter fly-by extremely challenging. The consequence of missing the Jupiter fly-by and launching in the late 2030s is the challenge of a fully-propulsive mission like UOP to have a feasible alternative interplanetary trajectory that reaches Uranus before 2050. As an alternative, to address the shortcomings of the fully propulsive mission, one can design a mission to Uranus using aerocapture. What is aerocapture: Aerocapture is an atmospheric maneuver that uses aerodynamic forces generated by flight within a planetary atmosphere to decelerate and achieve orbit insertion. Using aerocapture allows one to provide the change in velocity (Delta V) needed to slow down from the approach hyperbolic trajectory to achieve the desired captured orbit around the target planet using aerodynamic forces of the vehicle (lift and drag). Using aerodynamic forces instead of fully-propulsive maneuvers results in significant savings in the fuel. Furthermore, aerocapture can also allow one to consider interplanetary trajectories with higher approach hyperbolic velocities, thus reducing the mission transit times. Aerocapture as an enabling technology: To use aerocapture as an enabling technology for the Uranus exploration mission, one would require an integrated system-level design, including a Thermal Protection System (TPS), hardware needed for aerodynamic modulation, and autonomous Guidance, Navigation, and Control (GNC) systems. Aerocapture has yet to be demonstrated, despite considering it for several past missions. Recent advancements in TPS and GNC capabilities show the potential for using rigid, heritage entry vehicle configurations already flown at other planetary bodies for Uranus aerocapture. Aerocapture can be a robust technology that can deliver spacecraft to Uranus science orbits while substantially increasing on-orbit payload mass (more than 40\%) that can enable robust atmospheric entry probes. Additionally, the aerocapture maneuver would reduce the interplanetary transit time by 2–5 years (15-30\%) relative to fully-propulsive orbit insertion. Recent work has shown that one can conduct a flagship-class mission in a shorter period than fully-propulsive missions if using aerocapture. What does aerocapture bring in for a Uranus Mission: Using aerocapture for a Uranus orbiter and atmospheric probe mission can provide one with considerable propellant savings. Spacecraft in interplanetary trajectories to Uranus typically need an Delta V in orders of kilometers per second to insert into science orbit. One would require thousands of kilograms of fuel to achieve such a Delta V using a traditional fully propulsive maneuver, leaving less mass for payload during the mission launch. Aerocapture can reduce the propulsion needs by dissipating energy in the sizable atmosphere of Uranus without a significant mass increase due to the need for an aeroshell. One can use the mass savings achieved using the aerocapture to reduce the launch vehicle requirements. In addition, one can also have additional science instruments on the orbiter or create a robust instrumentation suite on an atmospheric probe that can significantly increase the science outcome of the Uranus exploration mission. Furthermore, since aerocapture performance is relatively insensitive to increases in hyperbolic excess velocity, one can design the interplanetary trajectory to arrive at Uranus faster, reducing the interplanetary transit time and operations cost. All these savings, achieved using aerocapture, could help fit a larger class mission, such as the Uranus mission within a smaller capital, e.g., a Flagship-class orbiter mission in a New Frontiers class capital. Summary: This talk will provide an overview of how aerocapture can enable the Uranus exploration mission. Specifically, this talk will discuss the latest advancements made in the Uranus aerocapture study, such as investigating interplanetary trajectories with higher hyperbolic approach velocities and their implications on the aero heating and the TPS design, incorporating FNPAG (an advanced numerical-predictor guidance) and comparison of multiple navigation approaches. In addition, this talk will focus on mechanical design that can house more than one atmospheric probe and the six degrees of freedom simulation of aerocapture at Uranus. Findings from a recent NASA Space Technology Mission Directorate (STMD)-funded activity studying the aerocapture as an enabling technology for a Uranus orbiter will be presented. Using the science payload recommended by the Decadal Survey for Uranus exploration, this work shows many improvements over the baseline fully-propulsive mission. These improvements include a shorter cruise phase, flexibility in launch opportunities late into the 2030s while reaching Uranus before the 2050 equinox for the desired science opportunities, and lower propellant mass needs. This talk will highlight how aerocapture can be utilized for Uranus science orbit insertion using a lower-risk, heritage entry vehicle configuration used extensively as a Mars entry, descent, and landing vehicle. Furthermore, this talk will explore how the demonstration of aerocapture at Earth can benefit the aerocapture-enabled Uranus mission.

Aerocapture↗

Quantum dynamics of cold atomic gas with SU(1,1) symmetry

Motivated by recent advances in quantum dynamics, we investigate the dynamics of the system with SU(1,1) symmetry. Instead of performing the time-ordered integral for the evolution operator of the time-dependent Hamiltonian, we show that the time evolution operator can be expressed as an SU(1,1) group element. Since the SU(1,1) group describes the “rotation” on a hyperbolic surface, the dynamics can be visualized on a Poincaré disk, a stereographic projection of the upper hyperboloid. As an example, we present the trajectory of the revival of Bose-Einstein condensation and that of the scale-invariant Fermi gas on the Poincaré disk. Further considering quantum gases in an oscillating lattice, we also study the dynamics of the system with time-dependent single-particle dispersion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Highly-scalable GPU-accelerated compressible reacting flow solver for modeling high-speed flows

Emerging supercomputing systems utilize a combination of central processing units (CPUs) and graphics processing units (GPUs) in an effort to reach exascale capabilities while minimizing the energy footprint of operating such systems. Such heterogeneous machines introduce new challenges for fluids solvers because the hardware architecture and operation of a GPU are fundamentally different from conventional CPUs. In this work, a general approach for efficient implementation of finite-volume based reacting flow solvers on such heterogeneous systems is presented. Three main challenges, namely, data access pattern, thread divergence, and thread safety, are addressed. Since compressible reacting flows require special methods to deal with chemical reactions, hyperbolic and nonlinear convection terms, and the presence of turbulence, specific algorithms that ensure GPU-based efficiency are developed. The approach is demonstrated on the widely available OpenFOAM open source software by modifying core algorithms for GPU accessibility. The scalability of the resulting solver, is demonstrated using practical test cases, including flow through a scramjet engine and the dynamics of a rotating detonation engine. Here, the solver provides near-ideal scaleup on a large number of GPUs (>3000), and extremely efficient use of the GPUs, with throughput nearly a constant even when processing a large number of control volumes.

42 ENGINEERING↗

The Space Transportation System summer environment on launch pad

This paper describes a 2D flow and thermalanalysis to determine the solar effect on the Space Shuttle launch components subsequent to the external tank (ET) loading operation in extremely hot conditions. An existing CFD code Parabolic Hyperbolic or Elliptical Numerical Integration Code Series was used in the study. The analysis was done for a 2D slice between planes perpendicular to the longitudinal axis of the STS and passing through the lower portions of the Redesigned Solid Rocket Motors (RSRMs), the ET, and the wing of the Orbiter. The results are presented as local and average values of the heat transfer coefficient, and the Nusselt number, and the surface temperature around the RSRMs and the ET. Solar heating effects increased the surface temperatures of the RSRMs by 9-11 F. Higher prelaunch surface temperatures measured on the east and west RSRMs (in the inboard region between the RSRMs and the ET) during 19 most recent launches of the STS are correlated as a function of the ambient temperature.

Ahmad, R. A.↗

A neutral-atmosphere composition experiment for the Atmosphere Explorer-C, -D, and -E.

The neutral-atmosphere composition experiment instrumentation is designed to obtain in-situ measurements of neutral thermosphere composition from Atmosphere Explorer-C, -D, and -E. The system is based on previously flown OGO-6 and San Marco-3 composition instruments. The mass-spectrometer sensor includes a gold-plated thermalizing chamber and ion source, a hyperbolic rod quadrupole analyzer, and an off-axis electron multiplier. Automatic ion-source sensitivity control and pulse-counting techniques provide density measurement capability from approximately 125 to 1000 km altitude. The normal operating mode includes measurement at all masses in the range of 1 to 44 amu, with emphasis on hydrogen, helium, oxygen, nitrogen, and argon.

Pelz, D. T.↗

Nozzle fabrication technique

This invention relates to techniques for fabricating hour glass throat or convergent divergent nozzle shapes, and more particularly to new and improved techniques for forming rocket nozzles from electrically conductive material and forming cooling channels in the wall thereof. The concept of positioning a block of electrically conductive material so that its axis is set at a predetermined skew angle with relation to a travelling electron discharge machine electrode and thereafter revolving the body about its own axis to generate a hyperbolic surface of revolution, either internal or external is novel. The method will generate a rocket nozzle which may be provided with cooling channels using the same control and positioning system. The configuration of the cooling channels so produced are unique and novel. Also the method is adaptable to nonmetallic material using analogous cutting tools, such as, water jet, laser, abrasive wire and hot wire.

Wells, Dennis L.↗

On least squares approximations to indefinite problems of the mixed type

A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.

Fix, G. J.↗

Science aspects of a mission beyond the planets

A mission out of the planetary system, launched about the year 2000, could provide valuable data concerning characteristics of the heliopause, the interstellar medium, stellar distances (by parallax measurements), low-energy cosmic rays, interplanetary gas distribution, and mass of the solar system. Secondary objectives include investigation of Pluto. Candidate science measurements, instruments, and instrument development needs are discussed. The mission should extend from 400 to 1000 AU from the sun. A heliocentric hyperbolic escape velocity of 50-100 km/sec or more is needed to attain this distance within a reasonable mission duration (20-50 years). The trajectory should be toward the incoming interstellar gas. For a year 2000 launch, a Pluto encounter and orbiter can be included. A second mission targeted parallel to the solar axis would also be worthwhile.

Jaffe, L. D.↗

Design of 3D Oxide–Metal Hybrid Metamaterial for Tailorable Light–Matter Interactions in Visible and Near-Infrared Region

Dielectric–metallic hybrid metamaterials exhibit extraordinary optical properties due to the light–matter interactions at the dielectric–metallic interfaces. The ability in precision control of the light–matter interactions in nanoscale is key to tailor the optical properties of hybrid metamaterials. In this work, a complex 3D framework of multilayered self-assembled BaTiO 3 (BTO)-Au hybrid thin films is demonstrated with such precision control of the light–matter interaction in nanoscale. Here, unique “bamboo-like” Au nanostructures are formed via the bilayer and trilayer stacking of BTO-Au hybrid layers with interlayers of SrTiO 3 , CeO 2 , or MgO. Different film strain states introduced by the three interlayers result in variable diameter and density of Au nanopillars. Both simulated and experimental optical data demonstrate the localized surface plasmon resonance change and hyperbolic dispersion wavelength shift in visible to near-infrared because of the effective tuning of the Au nanopillar aspect ratio and free electron density. The highly tunable optical properties along with the ferroelectric behavior and thermal robustness of the 3D hybrid film enable it to be a great candidate for multifunctional applications. This study demonstrates a unique 3D approach for precision optical property tuning and combined functionalities in oxide–metal metamaterial systems toward future integrated photonic and electronic devices.

36 MATERIALS SCIENCE↗

COMPASS Final Report: Enceladus Solar Electric Propulsion Stage

The results of the NASA Glenn Research Center (GRC) COllaborative Modeling and Parametric Assessment of Space Systems (COMPASS) internal Solar Electric Propulsion (SEP) stage design are documented in this report (Figure 1.1). The SEP Stage was designed to deliver a science probe to Saturn (the probe design was performed separately by the NASA Goddard Space Flight Center s (GSFC) Integrated Mission Design Center (IMDC)). The SEP Stage delivers the 2444 kg probe on a Saturn trajectory with a hyperbolic arrival velocity of 5.4 km/s. The design carried 30 percent mass, 10 percent power, and 6 percent propellant margins. The SEP Stage relies on the probe for substantial guidance, navigation and control (GN&C), command and data handling (C&DH), and Communications functions. The stage is configured to carry the probe and to minimize the packaging interference between the probe and the stage. The propulsion system consisted of a 1+1 (one active, one spare) configuration of gimbaled 7 kW NASA Evolutionary Xenon Thruster (NEXT) ion propulsion thrusters with a throughput of 309 kg Xe propellant. Two 9350 W GaAs triple junction (at 1 Astronomical Unit (AU), includes 10 percent margin) ultra-flex solar arrays provided power to the stage, with Li-ion batteries for launch and contingency operations power. The base structure was an Al-Li hexagonal skin-stringer frame built to withstand launch loads. A passive thermal control system consisted of heat pipes to north and south radiator panels, multilayer insulation (MLI) and heaters for the Xe tank. All systems except tanks and solar arrays were designed to be single fault tolerant.

Oleson, Steven R.↗

Loran-C expansion: Impact on precise time/time interval

On 16 May 1974, it was announced that Loran-C was chosen as the navigation system to serve the U. S. Coastal Confluence Zone. At the present time, reliable CONUS Loran-C groundwave timing coverage extends westward only about as far as Boulder, CO. The groundwave hyperbolic and timing coverage which will result from the planned CONUS expansion are illustrated. Time frames are provided. A status report on the planned reduction in Loran-C PTTI tolerances is presented.

Roeber, J. F., Jr.↗

An incompressible Navier-Stokes flow solver in three-dimensional curvilinear coordinate systems using primitive variables

An implicit, finite-difference computer code has been developed to solve the incompressible Navier-Stokes equations in a three-dimensional, curvilinear coordinate system. The pressure-field solution is based on the pseudo compressibility approach in which the time derivative pressure term is introduced into the mass conservation equation to form a set of hyperbolic equations. The solution procedure employs an implicit, approximate factorization scheme. The Reynolds stresses, that are uncoupled from the implicit scheme, are lagged by one time-step to facilitate implementing various levels of the turbulence model. Test problems for external and internal flows are computed, and the results are compared with existing experimental data. The application of this technique for general three-dimensional problems is then demonstrated.

Kwak, D.↗