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.
In this paper we derive fully-discrete Central WENO (CWENO) schemes for approximating solutions of one dimensional Hamilton-Jacobi (HJ) equations, which combine our previous works. We introduce third and fifth-order accurate schemes, which are the first central schemes for the HJ equations of order higher than two. The core ingredient is the derivation of our schemes is a high-order CWENO reconstructions in space.
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 roper summation-by-parts formula is found for the approximate derivative. A 'simultaneous approximation term' (SAT) 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.
The paper presents Alouette 2/ISIS 1 data which furnish evidence for the stimulation of high-order D(n) resonances (n greater than 4) by topside sounders, confirming the prediction of Osherovich (1990) based on an analogy with naturally occurring narrowband magnetospheric emissions. The results indicate that observations of stimulated ionospheric emissions and naturally occurring magnetospheric emissions can be used as complementary data sets to address such fundamental questions as the nature of the excitation mechanism for these emissions and the nature of the waves (i.e., whether the waves predominantly electrostatic or with a significant magnetic component present), and the question of the dominance among the Dn, Dn(+), and Dn(-) resonances and the conditions required for all three to be present at the same time.
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.
The summation-by-parts energy norm is used to establish a new class of high-order finite-difference techniques referred to here as 'cyclo-difference' techniques. These techniques are constructed cyclically from stable subelements, and require no numerical boundary conditions; when coupled with the simultaneous approximation term (SAT) boundary treatment, they are time asymptotically stable for an arbitrary hyperbolic system. These techniques are similar to spectral element techniques and are ideally suited for parallel implementation, but do not require special collocation points or orthogonal basis functions. The principal focus is on methods of sixth-order formal accuracy or less; however, these methods could be extended in principle to any arbitrary order of accuracy.
Runge-Kutta type formulas applied to restricted three-body problem achieving high order accuracy by m-fold differentiation and simple transformation
This research describes a new methodology for the extraction of a high-order, linear time invariant model, which allows the periodicity of the helicopter response to be accurately captured. This model provides the needed level of dynamic fidelity to permit an analysis and optimization of the AFCS and HHC algorithms. The key results of this study indicate that the closed-loop HHC system has little influence on the AFCS or on the vehicle handling qualities, which indicates that the AFCS does not need modification to work with the HHC system. However, the results show that the vibration response to maneuvers must be considered during the HHC design process, and this leads to much higher required HHC loop crossover frequencies. This research also demonstrates that the transient vibration responses during maneuvers can be reduced by optimizing the closed-loop higher harmonic control algorithm using conventional control system analyses.
The use of parallel computers for numerically solving flow fields has gained much importance in recent years. This paper introduces a new high order numerical scheme for computational fluid dynamics (CFD) specifically designed for parallel computational environments. A distributed MIMD system gives the flexibility of treating different elements of the governing equations with totally different numerical schemes in different regions of the flow field. The parallel decomposition of the governing operator to be solved is the primary parallel split. The primary parallel split was studied using a hypercube like architecture having clusters of shared memory processors at each node. The approach is demonstrated using examples of simple steady state incompressible flows. Future studies should investigate the secondary split because, depending on the numerical scheme that each of the processors applies and the nature of the flow in the specific subdomain, it may be possible for a processor to seek better, or higher order, schemes for its particular subcase.
We present a cost-effective scalable real-time wavefront control architecture based on off-the-shelf graphics processing units hosted in an ultra-low latency, high-bandwidth interconnect PC cluster environment composed of modules written in the component-oriented language of nesC. The architecture enables full-matrix reconstruction of the wavefront at up to 2 KHz with latency under 250 us for the PALM-3000 adaptive optics systems, a state-of-the-art upgrade on the 5.1 meter Hale Telescope that consists of a 64 x 64 subaperture Shack-Hartmann wavefront sensor and a 3368 active actuator high order deformable mirror in series with a 241 active actuator tweeter DM. The architecture can easily scale up to support much larger AO systems at higher rates and lower latency.
We present the first fifth order, semi-discrete central upwind method for approximating solutions of multi-dimensional Hamilton-Jacobi equations. Unlike most of the commonly used high order upwind schemes, our scheme is formulated as a Godunov-type scheme. The scheme is based on the fluxes of Kurganov-Tadmor and Kurganov-Tadmor-Petrova, and is derived for an arbitrary number of space dimensions. A theorem establishing the monotonicity of these fluxes is provided. The spacial discretization is based on a weighted essentially non-oscillatory reconstruction of the derivative. The accuracy and stability properties of our scheme are demonstrated in a variety of examples. A comparison between our method and other fifth-order schemes for Hamilton-Jacobi equations shows that our method exhibits smaller errors without any increase in the complexity of the computations.
The adaptive stencil concepts of ENO (Essentially Non-Oscillatory) methods are applied to the laminar Navier-Stokes equations to yield a high-order, time-accurate algorithm with a shock-capturing capability. The method targets problems in the areas of nonlinear acoustics, compressible transition, and turbulence which, due to the presence of shocks or complex geometries, are not easily solved by spectral methods. The present approach has been implemented and tested for the full three-dimensional Navier-Stokes equations in a transformed curvilinear coordinate system. Validation results are presented for a variety of problems which verify the method's accuracy properties and shock capturing capabilities, as well as demonstrate its use as a direct simulation tool.
Simulations of the early evolution of induced unstable disturbances in a Mach 6 hypersonic boundary layer are presented for the purposes of validating a high-order discontinuous Galerkin Navier-Stokes simulation code. The simulations are modeled after experiments performed in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. In these experiments, a forcing mechanism was employed to induce reproducible disturbances in a hypersonic boundary layer providing for the controlled study of the growth and breakdown of these disturbances into turbulent spots. Simulations revealed that the form and strength of the excitation can greatly influence the growth of the disturbance. In particular, at large forcing amplitude, the simulated forcing produces large advecting transients that appear to enhance the growth of the wave packet, relative to that of low amplitude forcing. Simulation results with large amplitude forcing agree well with experimental results while results from low amplitude forcing agree with linear stability theory.
An efficient strategy for propagating sonic boom signatures from a near-field Computational Fluid Dynamics (CFD) solution to the mid-field is presented. The method is based on a high-order accurate finite-difference discretization of the 3D Euler equations on a specially designed curvilinear grid and a single sweep space marching solution algorithm. The new approach leads to more than a factor of two reduction in overall computational resources compared to the current method used to propagate near-field sonic booms to the ground. Accuracy and efficiency of the near-field to mid-field process is demonstrated using a selection of test cases from the AIAA Sonic Boom Prediction Workshops. Azimuthal dependence of nonlinear wave propagation from the near-field to mid-field is analyzed along with its effects on the ground level noise.
The development of dynamic mesh capability for turbulent flow simulations using the Streamlined Upwind Petrov-Galerkin (SUPG) discretization is described. The current work extends previous research to include high-order spatial accuracy, including the satisfaction of the discrete geometric conservation law (GCL) on curved elements. Two closely-related schemes are described and the ability of these schemes to satisfy the GCL, while also maintaining temporal accuracy and conservation is assessed. Studies indicate that although one scheme discretizes the time derivative in conservative form, both schemes exhibit temporal conservation errors that decrease according to the expected design order of accuracy. The source of the temporal conservation errors is examined, and it is demonstrated that many finite-volume and finite-element schemes can also be expected to have difficulty strictly satisfying conservation in time. The effects on conservation are examined and, while present in the simulations, are seen to be negligible for the problems considered.
We present results of the workshop's benchmark problems of the category 2. This category of problems is designed to test the nonlinear wave propagation properties of a computational scheme. We chose three high order spatially accurate algorithms for our computations. These are the Dispersion- Relation- Preserving (DRP) scheme proposed by Tam and his colleagues, a fourth order extension of the MacCormack scheme proposed by Gottlieb and Turkel and an Essentially Non Oscillatory (ENO) scheme proposed by Shu and Osher.
A method for the prediction of acoustic scatter from complex geometries is presented. The discontinuous Galerkin method provides a framework for the development of a high-order method using unstructured grids. The method's compact form contributes to its accuracy and efficiency, and makes the method well suited for distributed memory parallel computing platforms. Mesh refinement studies are presented to validate the expected convergence properties of the method, and to establish the absolute levels of a error one can expect at a given level of resolution. For a two-dimensional shear layer instability wave and for three-dimensional wave propagation, the method is demonstrated to be insensitive to mesh smoothness. Simulations of scatter from a two-dimensional slat configuration and a three-dimensional blended-wing-body demonstrate the capability of the method to efficiently treat realistic geometries.
A method is presented for curve fitting nonparametric frequency domain data so as to identify a parametric model composed of two models in parallel, where each model has dynamics in a specified portion of the frequency band. This decomposition overcomes the problem of numerical sensitivity since lower order polynomials can be used compared to existing methods which estimate the model as a single entity. Consequently, composite curve fitting is useful for frequency domain identification of high-order systems and/or systems whose dynamics are spread over a large bandwidth. The approach can be extended to identify an arbitrary number of parallel subsystems in specified frequency regimes.
Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.