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.
A finite difference scheme is developed to find an approximate solution of two similar hyperbolic equations, namely a first-order plane wave and spherical wave problem. Finite difference approximations are made for both the space and time derivatives. The result is a conditionally stable equation yielding an exact solution when the Courant number is set to one.
We consider the problem of obtaining integral representation of feedback operators for damped hyperbolic control systems. We show that for the wave equation with Kelvin-Voigt damping and non-compact input operator, the feedback gain operator is Hilbert-Schmidt. This result is then used to provide an explicit integral representation for the feedback operator in terms of functional gains. Numerical results are given to illustrate the role that damping plays in the smoothness of these gains.
We present a systematic method for constructing boundary conditions (numerical and physical) of the required accuracy, for compact (Pade-like) high-order finite-difference schemes for hyperbolic systems. First a proper summation-by-parts formula is found for the approximate derivative. A 'simultaneous approximation term' is then introduced to treat the boundary conditions. This procedure leads to time-stable schemes even in the system case. An explicit construction of the fourth-order compact case is given. Numerical studies are presented to verify the efficacy of the approach.
This paper presents an hp-adaptive discontinuous Galerkin method for linear hyperbolic conservation laws. A priori and a posteriori error estimates are derived in mesh-dependent norms which reflect the dependence of the approximate solution on the element size (h) and the degree (p) of the local polynomial approximation. The a posteriori error estimate, based on the element residual method, provides bounds on the actual global error in the approximate solution. The adaptive strategy is designed to deliver an approximate solution with the specified level of error in three steps. The a posteriori estimate is used to assess the accuracy of a given approximate solution and the a priori estimate is used to predict the mesh refinements and polynomial enrichment needed to deliver the desired solution. Numerical examples demonstrate the reliability of the a posteriori error estimates and the effectiveness of the hp-adaptive strategy.
In these lecture notes we describe the construction, analysis, and application of ENO (Essentially Non-Oscillatory) and WENO (Weighted Essentially Non-Oscillatory) schemes for hyperbolic conservation laws and related Hamilton- Jacobi equations. ENO and WENO schemes are high order accurate finite difference schemes designed for problems with piecewise smooth solutions containing discontinuities. The key idea lies at the approximation level, where a nonlinear adaptive procedure is used to automatically choose the locally smoothest stencil, hence avoiding crossing discontinuities in the interpolation procedure as much as possible. ENO and WENO schemes have been quite successful in applications, especially for problems containing both shocks and complicated smooth solution structures, such as compressible turbulence simulations and aeroacoustics. These lecture notes are basically self-contained. It is our hope that with these notes and with the help of the quoted references, the reader can understand the algorithms and code them up for applications.
This paper presents a new discretization scheme for hyperbolic systems of conservations laws. It satisfies the TVD property and relies on the new high-resolution mechanism which is compatible with the genuinely multidimensional approach proposed recently. This work can be regarded as a first step towards extending the genuinely multidimensional approach to unsteady problems. Discontinuity capturing capabilities and accuracy of the scheme are verified by a set of numerical tests.
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.
The lecture material will discuss the application of one-dimensional approximate Riemann solutions and high order accurate data reconstruction as building blocks for solving multi-dimensional hyperbolic equations. This building block procedure is well-documented in the nationally available literature. The relevant stability and convergence theory using positive operator analysis will also be presented. All participants in the minisymposium will be asked to solve one or more generic test problems so that a critical comparison of accuracy can be made among differing approaches.
A non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented. The technique is based on the first principle of non-reflecting, plane wave propagation and the hyperbolicity of the Euler equation system. The NRBC is simple and effective, provided the numerical scheme maintains locally a C(sup 1) continuous solution at the boundary. Several numerical examples in 1D, 2D, and 3D space are illustrated to demonstrate its robustness in practical computations.
This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.
Interstellar meteoroids, solid particles arriving from outside our Solar System, are not easily distinguished from local meteoroids. A velocity above the escape velocity of the Sun is often used as an indicator of a possible interstellar origin. We demonstrate that the gravitational slingshot effect, resulting from the passage of local meteoroid near a planet, can produce hyperbolic meteoroids at the Earth s orbit with excess velocities comparable to those expected of interstellar meteoroids.
Ablative materials are required to protect a space vehicle from the extreme temperatures encountered during the most demanding (hyperbolic) atmospheric entry velocities, either for probes launched toward other celestial bodies, or coming back to Earth from deep space missions. To that effect, the resin-impregnated carbon ablator (RICA) is a high-temperature carbon/phenolic ablative thermal protection system (TPS) material designed to use modern and commercially viable components in its manufacture. Heritage carbon/phenolic ablators intended for this use rely on materials that are no longer in production (i.e., Galileo, Pioneer Venus); hence the development of alternatives such as RICA is necessary for future NASA planetary entry and Earth re-entry missions. RICA s capabilities were initially measured in air for Earth re-entry applications, where it was exposed to a heat flux of 14 MW/sq m for 22 seconds. Methane tests were also carried out for potential application in Saturn s moon Titan, with a nominal heat flux of 1.4 MW/sq m for up to 478 seconds. Three slightly different material formulations were manufactured and subsequently tested at the Plasma Wind Tunnel of the University of Stuttgart in Germany (PWK1) in the summer and fall of 2010. The TPS integrity was well preserved in most cases, and results show great promise.
Spacecraft modeling, a critically important portion in validating planned spacecraft activities, is currently carried out using a time consuming method of mission to mission model implementations and integration. A current project in early development, Integrated Spacecraft Analysis (ISCA), aims to remedy this hindrance by providing reusable architectures and reducing time spent integrating models with planning and sequencing tools. The principle objective of this internship was to develop a user interface for an experimental ontology-based structure visualization of navigation and attitude control system modeling software. To satisfy this, a number of tree and graph visualization tools were researched and a Java based hyperbolic graph viewer was selected for experimental adaptation. Early results show promise in the ability to organize and display large amounts of spacecraft model documentation efficiently and effectively through a web browser. This viewer serves as a conceptual implementation for future development but trials with both ISCA developers and end users should be performed to truly evaluate the effectiveness of continued development of such visualizations.
A non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented. The technique is based on the hyperbolicity of the Euler equation system and the first principle of plane (simple) wave propagation. The NRBC is simple and effective, provided the numerical scheme maintains locally a C(sup 1) continuous solution at the boundary. Several numerical examples in ID, 2D and 3D space are illustrated to demonstrate its robustness in practical computations.
Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.
Not Available
A first-order perturbation theory exists for the flyby problem with the Born approximation as the reference orbit.