SEARCH · Engineering Papers
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.
The existing theories of comet dynamics held that large numbers of comets are ejected to interstellar space and that other stellar systems would do the same; however, no comet on a hyperbolic orbit has ever been observed. Calculations developed by Duncan et al. are used to show that most comets thought to be ejected would actually be captured by the Oort Cloud.
The inherent nonlinear character of adaptive systems poses serious theoretical problems for the analysis of their dynamics. On the other hand, the importance of their dynamic behavior is directly related to the practical interest in predicting such undesirable phenomena as nonlinear oscillations, abrupt transients, intermittence or a high sensitivity with respect to initial conditions. A geometrical/qualitative description of the phase portrait of a discrete-time adaptive system with unmodeled disturbances is given. For this, the motions in the phase space are referred to normally hyperbolic (structurally stable) locally invariant sets. The study is complemented with a local stability analysis of the equilibrium point and periodic solutions. The critical character of adaptive systems under rather usual working conditions is discussed. Special emphasis is put on the causes leading to intermittence. A geometric interpretation of the effects of some commonly used palliatives to this problem is given. The 'dead-zone' approach is studied in more detail. The predicted dynamics are compared with simulation results.
The Pioneer Venus Orbiter Neutral Mass Spectrometer (ONMS) is designed to measure the vertical and horizontal density variations of the major neutral constituents in the upper atmosphere of Venus. The mass spectrometer sensor includes a retarding potential ion source, hyperbolic quadrupole rod analyzer, and electron multiplier detector. The supporting electronic system consists of hybrid integrated circuits to reduce weight and power. The ONMS instrument was launched aboard the Pioneer Venus Orbiter on May 20, 1978, and turned on in orbit around Venus on December 4, 1978. It has operated flawlessly for over a Venus year (243 earth days) and has returned data of the composition of the major constituents in the Venus atmosphere between the altitudes of 150 and 350 km.
It has often been proposed that a vehicle returning from Mars will use aerobraking in the Earth's atmosphere to dissipate hyperbolic excess velocity to capture into Earth orbit. Here a different system for dissipating excess velocity without expenditure of reaction mass, magnetobraking, is proposed. Magnetobraking uses the force on an electrodynamic tether in the Earth's magnetic field to produce thrust. An electrodynamic tether is deployed from the spacecraft as it approaches the Earth. The Earth's magnetic field produces a force on electrical current in the tether. If the tether is oriented perpendicularly to the Earth's magnetic field and to the direction of motion of the spacecraft, force produced by the Earth's magnetic field can be used to either brake or accelerate the spacecraft without expenditure of reaction mass. The peak acceleration on the Mars return is 0.007 m/sq sec, and the amount of braking possible is dependent on the density and current-carrying capacity of the tether, but is independent of length. A superconducting tether is required. The required critical current is shown to be within the range of superconducting technology now available in the laboratory.
NASA is currently building the Space Launch System (SLS) Block-1 launch vehicle for the Exploration Mission 1 (EM-1) test flight. The design of the next evolution of SLS, Block-1B, is well underway. The Block-1B vehicle is more capable overall than Block-1; however, the relatively low thrust-to-weight ratio of the Exploration Upper Stage (EUS) presents a challenge to the Powered Explicit Guidance (PEG) algorithm used by Block-1. To handle the long burn durations (on the order of 1000 seconds) of EUS missions, two algorithms were examined. An alternative algorithm, OPGUID, was introduced, while modifications were made to PEG. A trade study was conducted to select the guidance algorithm for future SLS vehicles. The chosen algorithm needs to support a wide variety of mission operations: ascent burns to LEO, apogee raise burns, trans-lunar injection burns, hyperbolic Earth departure burns, and contingency disposal burns using the Reaction Control System (RCS). Additionally, the algorithm must be able to respond to a single engine failure scenario. Each algorithm was scored based on pre-selected criteria, including insertion accuracy, algorithmic complexity and robustness, extensibility for potential future missions, and flight heritage. Monte Carlo analysis was used to select the final algorithm. This paper covers the design criteria, approach, and results of this trade study, showing impacts and considerations when adapting launch vehicle guidance algorithms to a broader breadth of in-space operations.
This report shows that when considered in terms of optical geometry, radially parabolic and radially hyperbolic mirrors used as mirrors for Cassegrain and parabolic antennas possess values similar to common conical horn reflector antennas.
Specific problems in the theory of relativity are often simplified by an appropriate choice of the coordinate system. Restricted conformal coordinates provide an especially simple analysis of motion with uniform acceleration, known as hyperbolic motion. Conformal coordinates x', t' may be obtained from Cartesian coordinates x, t by the transformation x'+ct'=F(x+ct) and x'-ct'=G(x-ct), where c is the velocity of light. A variable motion of the x' system is determined by the choice of the functions F and G.
This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.
Here, the effectiveness of continuous vibro-impact forcing representations for the cantilevered pipe that conveys fluid is explored and analyzed. The previously accepted forcing model utilizing a smoothened trilinear spring is estimated using three continuous forcing representations, namely, polynomial, rational polynomial, and hyperbolic tangent. The accuracy of the estimated forcing functions is investigated and analyzed by calculating the root mean square error, and bifurcation diagrams are generated and compared to the nominal system. Additionally, the dynamic response of the system is further characterized using Poincare maps, power spectra, and basins of attraction. Once all continuous forcing representations are analyzed and compared to the nominal system, the computational cost of each method is examined, and further limitations of the hyperbolic tangent method are discovered. It is proved that the hyperbolic tangent forcing representation most accurately captures the dynamic response of the pipeline, and the least accurate representation is the rational polynomial representation. Additionally, considerable computational cost is saved when employing the hyperbolic tangent representation compared to the discontinuous representation.
A gas-kinetic method for the hyperbolic-elliptic equations is presented in this paper. In the mixed type system, the co-existence and the phase transition between liquid and gas are described by the van der Waals-type equation of state (EOS). Due to the unstable mechanism for a fluid in the elliptic region, interface between the liquid and gas can be kept sharp through the condensation and evaporation process to remove the "averaged" numerical fluid away from the elliptic region, and the interface thickness depends on the numerical diffusion and stiffness of the phase change. A few examples are presented in this paper for both phase transition and multifluid interface problems.
A scheme is presented for the automatic generation of structured overset meshes ongeometries that are defined by Boundary Representation (BRep) solids. The surface meshsystem consists of face, edge and node meshes corresponding to the three respective basicBRep entity types. A cut-cell method is introduced to improve robustness of the on-geometrydetermination test for a face mesh grid point. A geometric component tagging scheme is utilizedto enhance local grid point distribution on a configuration with a large range of geometricscales. A cap grid topology is automatically utilized around the trailing edge of wing and tailtips to enhance mesh quality and to enable more effective surface coverage. Robustness of thehyperbolic surface marching method is improved by replacing the point projection scheme witha walking scheme. Relaxation of surface grid spacing at concave corners enables automatedhigh quality hyperbolic volume mesh creation. Domain connectivity is automatically performedon the surface mesh system. A variety of test cases are presented including a re-entry capsule,two models of the Juncture Flow Experiment wing-body, five rotorcraft concept vehicles, andvarious components of the High-Lift Common Research Model from the High-Lift PredictionWorkshop 4.
This is the third paper in a series in which we develop machine learning (ML) moment closure models for the radiative transfer equation. In our previous work (Huang et al. in J Comput Phys 453:110941, 2022), we proposed an approach to learn the gradient of the unclosed high order moment, which performs much better than learning the moment itself and the conventional $P_N$ closure. However, while the ML moment closure has better accuracy, it is not able to guarantee hyperbolicity and has issues with long time stability. In our second paper (Huang et al., in: Machine learning moment closure models for the radiative transfer equation II: enforcing global hyperbolicity in gradient based closures, 2021. arXiv:2105.14410), we identified a symmetrizer which leads to conditions that enforce that the gradient based ML closure is symmetrizable hyperbolic and stable over long time. The limitation of this approach is that in practice the highest moment can only be related to four, or fewer, lower moments. In this paper, we propose a new method to enforce the hyperbolicity of the ML closure model. Motivated by the observation that the coefficient matrix of the closure system is a lower Hessenberg matrix, we relate its eigenvalues to the roots of an associated polynomial. Here, we design two new neural network architectures based on this relation. The ML closure model resulting from the first neural network is weakly hyperbolic and guarantees the physical characteristic speeds, i.e., the eigenvalues are bounded by the speed of light. The second model is strictly hyperbolic and does not guarantee the boundedness of the eigenvalues. Several benchmark tests including the Gaussian source problem and the two-material problem show the good accuracy, stability and generalizability of our hyperbolic ML closure model.
Abstract Teleporting energy to remote locations is new challenge for quantum information science and technology. Developing a method for transferring local energy in laboratory systems to remote locations will enable non‐trivial energy flows in quantum networks. From the perspective of quantum information engineering, we propose a method for distributing local energy to a large number of remote nodes using hyperbolic geometry. Hyperbolic networks are suitable for energy allocation in large quantum networks since the number of nodes grows exponentially. To realise long‐range quantum energy teleportation (QET), we propose a hybrid method of quantum state telepotation and QET. By transmitting local quantum information through quantum teleportation and performing conditional operations on that information, QET can theoretically be realized independent of geographical distance. The method we present will provide new insights into new applications of future large‐scale quantum networks and potential applications of quantum physics to information engineering.
LORAN-C a highly accurate radio navigation positioning system which operates at an assigned frequency of 10 kHz, and provides phase-coded pulses to develop hyperbolic time-difference lines-of-position (LOP's) was evaluated. LORAN-C provides precise time and time interval to within plus or minus 5 microseconds of UTC. The steps taken to plan, install, operate, and maintain the LORAN-C system up to the year 2000 are discussed. Topics included in the discussion were: theory of operation, timing, chain lanning, group repetition interval, coding delay versus emission delay, chain calibration, chart verification, system accuracy, signal reliability, and future developments.
In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.
Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.
Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.