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 specific problem objective was to locate the sources of radiated electric field from lightning using an overdetermined set of measurements of time-of-arrival. A similar problem exists for epicentral sources in earthquake location, acoustic sources of thunder, and terrestrial navigation using LORAN and GPS.
This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques, we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.
Momentum-exchange/electrodynamic reboost (MXER) tether technology is currently being pursued to dramatically lower the launch mass and cost of interplanetary scientific spacecraft. A spacecraft boosted from LEO to a high-energy orbit by a MXER tether has most of the orbital energy it needs to escape the Earth's gravity well. However, the final targeting of the spacecraft to its eventual trajectory, and some of the unique issues brought on by the tether boost, are the subjects of this paper.
This lecture considers a-posteriori error estimates for the numerical solution of conservation laws with time invariant constraints such as those arising in magnetohydrodynamics (MHD) and gravitational physics. Using standard duality arguments, a-posteriori error estimates for the discontinuous Galerkin finite element method are then presented for MHD with solenoidal constraint. From these estimates, a procedure for adaptive discretization is outlined. A taxonomy of Green's functions for the linearized MHD operator is given which characterizes the domain of dependence for pointwise errors. The extension to other constrained systems such as the Einstein equations of gravitational physics are then considered. Finally, future directions and open problems are discussed.
Explore the source record for details and available documents.
Space probes in the Solar System have experienced unexpected changes in velocity known as the flyby anomaly [1], as well as shifts in acceleration referred to as the Pioneer anomaly [2-4]. In the case of Earth flybys, ESA s Rosetta spacecraft experienced the flyby effect and NASA s Galileo and NEAR satellites did the same, although MESSENGER did not possibly due to a latitudinal property of gravity assists. Measurements indicate that both anomalies exist, and explanations have varied from the unconventional to suggestions that new physics in the form of dark matter might be the cause of both [5]. Although dark matter has been studied for over 30 years, there is as yet no strong experimental evidence supporting it [6]. The existence of dark matter will certainly have a significant impact upon ideas regarding the origin of the Solar System. Hence, the subject is very relevant to planetary science. We will point out here that one of the fundamental problems in science, including planetary physics, is consistency. Using the well-known virial theorem in astrophysics, it will be shown that present-day concepts of orbital mechanics and cosmology are not consistent for reasons having to do with the flyby anomaly. Therefore, the basic solution regarding the anomalies should begin with addressing the inconsistencies first before introducing new physics.
This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.
The scientific desire for global coverage of moons such as Jupiter's Galilean moons or Saturn's Titan has invariably led to the design of orbiter missions. These orbiter missions require a large amount of propellant needed to insert into orbit around such small bodies, and for a given launch vehicle, the additional propellant mass takes away from mass that could otherwise be used for scientific instrumentation on a multiple flyby-only mission. This paper will present methods--expanding upon techniques developed for the design of the Cassini prime and extended missions--to obtain near global moon coverage through multiple flybys. Furthermore we will show with proper instrument suite selection, a flyby-only mission can provide science return similar (and in some cases greater) to that of an orbiter mission.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
The research work has advanced the inversion-based guidance theory for: systems with non-hyperbolic internal dynamics; systems with parameter jumps; and systems where a redesign of the output trajectory is desired. A technique to achieve output tracking for nonminimum phase linear systems with non-hyperbolic and near non-hyperbolic internal dynamics was developed. This approach integrated stable inversion techniques, that achieve exact-tracking, with approximation techniques, that modify the internal dynamics to achieve desirable performance. Such modification of the internal dynamics was used (a) to remove non-hyperbolicity which is an obstruction to applying stable inversion techniques and (b) to reduce large preactuation times needed to apply stable inversion for near non-hyperbolic cases. The method was applied to an example helicopter hover control problem with near non-hyperbolic internal dynamics for illustrating the trade-off between exact tracking and reduction of preactuation time. Future work will extend these results to guidance of nonlinear non-hyperbolic systems. The exact output tracking problem for systems with parameter jumps was considered. Necessary and sufficient conditions were derived for the elimination of switching-introduced output transient. While previous works had studied this problem by developing a regulator that maintains exact tracking through parameter jumps (switches), such techniques are, however, only applicable to minimum-phase systems. In contrast, our approach is also applicable to nonminimum-phase systems and leads to bounded but possibly non-causal solutions. In addition, for the case when the reference trajectories are generated by an exosystem, we developed an exact-tracking controller which could be written in a feedback form. As in standard regulator theory, we also obtained a linear map from the states of the exosystem to the desired system state, which was defined via a matrix differential equation.
Some results of the International Heliophysical Year (IHY) Coordinated Investigation Program (CIP) number 65 Meteors in the Earth Atmosphere and Meteoroids in the Solar System are presented. The problem of hyperbolic and near-parabolic orbits is discussed. Some possibilities for the solution of this problem can be obtained from the radar observation of faint meteors. The limiting magnitude of the Kharkov, Ukraine, radar observation program in the 1970 s was +12, resulting in a very large number of meteors being detected. 250,000 orbits down to even fainter limiting magnitude were determined in the 1972-78 period in Kharkov (out of them 7,000 are hyperbolic). The hypothesis of hyperbolic meteors was confirmed. In some radar meteor observations 1 10% of meteors are hyperbolic meteors. Though the Advanced Meteor Orbit Radar (AMOR, New Zealand) and Canadian Meteor Orbit Radar (CMOR, Canada) have accumulated millions of meteor orbits, there are difficulties in comparing the radar observational data obtained from these three sites (New Zealand, Canada, Kharkov). A new global program International Space Weather Initiative (ISWI) has begun in 2010 (http://www.iswi-secretariat.org). Today it is necessary to create the unified radar catalogue of nearparabolic and hyperbolic meteor orbits in the framework of the ISWI, or any other different way, in collaboration of Ukraine, Canada, New Zealand, the USA and, possibly, Japan. Involvement of the Virtual Meteor Observatory (Netherlands) and Meteor Data Centre (Slovakia) is desirable too. International unified radar catalogue of near-parabolic and hyperbolic meteor orbits will aid to a major advance in our understanding of the ecology of meteoroids within the Solar System and beyond.
In this paper, we develop second- and third-order non-oscillatory shock-capturing hyperbolic residual distribution schemes for irregular triangular grids, extending our second- and third-order schemes to discontinuous problems. We present extended first-order N- and Rusanov-scheme formulations for hyperbolic advection-diffusion system, and demonstrate that the hyperbolic diffusion term does not affect the solution of inviscid problems for vanishingly small viscous coefficient. We then propose second- and third-order blended hyperbolic residual-distribution schemes with the extended first-order Rusanov-scheme. We show that these proposed schemes are extremely accurate in predicting non-oscillatory solutions for discontinuous problems. We also propose a characteristics-based nonlinear wave sensor for accurately detecting shocks, compression, and expansion regions. Using this proposed sensor, we demonstrate that the developed hyperbolic blended schemes do not produce entropy-violating solutions (unphysical stocks). We then verify the design order of accuracy of these blended schemes on irregular triangular grids.
In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.