Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

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 361 records · Page 20

Validation of High-Fidelity CFD/CAA Framework for Launch Vehicle Acoustic Environment Simulation against Scale Model Test Data

A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed for launch vehicle liftoff acoustic environment predictions. The framework couples the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin solver developed in the same production framework, Loci/THRUST, to accurately resolve and propagate acoustic physics across the entire launch environment. Time-accurate, Hybrid RANS/LES CFD modeling is applied for predicting the acoustic generation physics at the plume source, and a high-order accurate unstructured discontinuous Galerkin (DG) method is employed to propagate acoustic waves away from the source across large distances using high-order accurate schemes. The DG solver is capable of solving 2nd, 3rd, and 4th order Euler solutions for non-linear, conservative acoustic field propagation. Initial application testing and validation has been carried out against high resolution acoustic data from the Ares Scale Model Acoustic Test (ASMAT) series to evaluate the capabilities and production readiness of the CFD/CAA system to resolve the observed spectrum of acoustic frequency content. This paper presents results from this validation and outlines efforts to mature and improve the computational simulation framework.

Liever, Peter A.↗

Validation of High-Fidelity CFD/CAA Framework for Launch Vehicle Acoustic Environment Simulation against Scale Model Test Data

A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed for launch vehicle liftoff acoustic environment predictions. The framework couples the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate Discontinuous Galerkin solver developed in the same production framework, Loci/THRUST, to accurately resolve and propagate acoustic physics across the entire launch environment. Time-accurate, Hybrid RANS/LES CFD modeling is applied for predicting the acoustic generation physics at the plume source, and a high-order accurate unstructured mesh Discontinuous Galerkin (DG) method is employed to propagate acoustic waves away from the source across large distances using high-order accurate schemes. The DG solver is capable of solving 2nd, 3rd, and 4th order Euler solutions for non-linear, conservative acoustic field propagation. Initial application testing and validation has been carried out against high resolution acoustic data from the Ares Scale Model Acoustic Test (ASMAT) series to evaluate the capabilities and production readiness of the CFD/CAA system to resolve the observed spectrum of acoustic frequency content. This paper presents results from this validation and outlines efforts to mature and improve the computational simulation framework.

Liever, Peter A.↗

Numerical simulation of instability and transition physics

The study deals with the algorithm technology used in instability and transition simulations. Discretization methods are outlined, and attention is focused on high-order finite-difference methods and high-order centered-difference formulas. One advantage of finite-difference methods over spectral methods is thought to be in implementation of nonrigorous boundary conditions. It is suggested that the next significant advances in the understanding of transition physics and the ability to predict transition will come with more physically-realistic simulations. Compressible-flow algorithms are discussed, and it is noted that with further development, exploration of bypass mechanism on simple bodies at high speed would be possible.

Streett, C. L.↗

Towards a Third-Order Accurate, Second-Derivative-Free, Shock-Capturing Finite-Volume Method for Hypersonic Flows on Tetrahedral Grids

In this paper, we report progress in the development of a third-order accurate, second-derivative-free, shock-capturing finite-volume solver for three-dimensional unstructured grids. The method is economical in the sense that the computation and storage of second derivatives are not required for third-order accuracy. It is based on point-valued numerical solutions stored at cells, gradients computed and stored at nodes, and an efficient projected-derivative formula that eliminates the need for second derivatives in a quadratic solution interpolation. The projected-derivative formula is also used to eliminate second derivatives from a high-order flux quadrature formula, so that it can be implemented conveniently in the form of a numerical flux at a face center plus a correction term. Similarly, a high-order source quadrature formula can also be implemented in the form of a cell-center point evaluation plus a similar correction term. These features make it relatively straightforward to extend an existing second-order finite-volume code to third-order. This paper reports progress of implementing the method in the NASA VULCAN-CFD code and discusses the implementation of a high-order accurate limiter for shock capturing.

Weighted Least-Squares↗

A Technique of Treating Negative Weights in WENO Schemes

High order accurate weighted essentially non-oscillatory (WENO) schemes have recently been developed for finite difference and finite volume methods both in structural and in unstructured meshes. A key idea in WENO scheme is a linear combination of lower order fluxes or reconstructions to obtain a high order approximation. The combination coefficients, also called linear weights, are determined by local geometry of the mesh and order of accuracy and may become negative. WENO procedures cannot be applied directly to obtain a stable scheme if negative linear weights are present. Previous strategy for handling this difficulty is by either regrouping of stencils or reducing the order of accuracy to get rid of the negative linear weights. In this paper we present a simple and effective technique for handling negative linear weights without a need to get rid of them.

Shi, Jing↗

Microgravity Experiments Being Developed for Microscopic Study of Colloidal Spheres

Microscopic spheres suspended in liquid become highly ordered under the proper conditions. Such collections of particles, called colloidal suspensions or colloids, are the subject of a series of ongoing microgravity experiments at the NASA Lewis Research Center. By studying the way these colloidal suspensions order themselves, scientists can better understand how atoms of a liquid become ordered to form a solid. In addition, highly ordered colloids have special properties that may make them useful in future hightech applications. Work is underway at Lewis to develop an optical microscope to view these colloidal suspensions sphere by sphere in microgravity.

Rogers, Richard B.↗

A spectral multi-domain technique with application to generalized curvilinear coordinates

Spectral collocation methods have proven to be efficient discretization schemes for many aerodynamic and fluid mechanic problems. The high order accuracy and resolution shown by these methods allows one to obtain engineering accuracy solutions on coarse meshes, or alternatively, to obtain solutions with very small error. One drawback to these techniques was the requirement that a complicated physical domain must map into a simple computational domain for discretization. This mapping must be smooth if the high order accuracy and expontential convergence rates associated with spectral methods are to be preserved. Additionally even smooth stretching transformations can decrease the accuracy of a spectral method, if the stretching is severe. A further difficulty with spectral methods was in their implementation on parallel processing computers, where efficient spectral algorithms were lacking. The above restrictions are overcome by splitting the domain into regions, each of which preserve the advantages of spectral collocation, and allow the ratio of the mesh spacing between regions to be several orders of magnitude higher than allowable in a single domain. Such stretchings would be required to resolve the thin viscous region in an external aerodynamic problem. Adjoining regions are interfaced by enforcing a global flux balance which preserves high-order continuity of the solution, regardless of the type of the equations being solved.

Macaraeg, M. G.↗

Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction

High order accurate finite-volume schemes for solving the Euler equations of gasdynamics are developed. Central to the development of these methods are the construction of a k-exact reconstruction operator given cell-averaged quantities and the use of high order flux quadrature formulas. General polygonal control volumes (with curved boundary edges) are considered. The formulations presented make no explicit assumption as to complexity or convexity of control volumes. Numerical examples are presented for Ringleb flow to validate the methodology.

Barth, Timothy J.↗

Modeling with Liapunov functions.

Behavior of high-order linear control systems analyzed using liapunov second method, by finding low-order model with closely approximate response behavior of high order linear control systems analyzed, using Liapunov second method, by finding low order model with closely approximate response

LIAPUNOV FUNCTION↗

Liquid Crystal on Silicon Wavefront Corrector

A low cost, high resolution, liquid crystal on silicon, spatial light modulator has been developed for the correction of huge aberrations in an optical system where the polarization dependence and the chromatic nature are tolerated. However, the overall system performance suggests that this device is also suitable for real time correction of aberration in human eyes. This device has a resolution of 1024 x 768, and is driven by an XGA display driver. The effective stroke length of the device is 700 nm and 2000 nm for the visible and IR regions of the device, respectively. The response speeds are 50 Hz and 5 Hz, respectively, which are fast enough for real time adaptive optics for aberrations in human eyes. By modulating a wavefront of 2 pi, this device can correct for arbitrary high order wavefront aberrations since the 2-D pixel array is independently controlled by the driver. The high resolution and high accuracy of the device allow for diffraction limited correction of the tip and tilt or defocus without an additional correction loop. We have shown that for every wave of aberration, an 8 step blazed grating is required to achieve high diffraction efficiency around 80%. In light of this, up to 125 waves peak to valley of tip and tilt can be corrected if we choose the simplest aberration. Corrections of 34 waves of aberration, including high order Zernicke terms in a high magnification telescope, to diffraction limited performance (residual wavefront aberration less than 1/30 lambda at 632.8 nm) have been observed at high efficiency.

Pouch, John↗

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗

A practical assessment of spectral accuracy for hyperbolic problems with discontinuities

Physical and transform space filtering has been applied to the Fourier spectral collocation solution of the constant coefficient scalar wave equation with a discontinuous initial condition. High order accuracy can be extracted from the unfiltered solution. Smooth, high order Fourier space filtering gives expected polynomial order solutions away from the discontinuity. Spectral accuracy is observed with the physical space filter of Gottlieb and Tadmor.

Kopriva, David A.↗

A practical assessment of spectral accuracy for hyperbolic problems with discontinuities

Physical and transform space filtering has been applied to the Fourier spectral collocation solution of the constant coefficient scalar wave equation with a discontinuous initial condition. High order accuracy can be extracted from the unfiltered solution. Smooth, high order Fourier space filtering gives expected polynomial order solutions away from the discontinuity. Spectral accuracy is observed with the physical space filter of Gottlieb and Tadmor.

Kopriva, David A.↗

Closed-Loop System Identification Experience for Flight Control Law and Flying Qualities Evaluation of a High Performance Fighter Aircraft

This paper highlights some of the results and issues associated with estimating models to evaluate control law design methods and design criteria for advanced high performance aircraft. Experimental fighter aircraft such as the NASA-High Alpha Research Vehicle (HARV) have the capability to maneuver at very high angles of attack where nonlinear aerodynamics often predominate. HARV is an experimental F/A-18, configured with thrust vectoring and conformal actuated nose strakes. Identifying closed-loop models for this type of aircraft can be made difficult by nonlinearities and high order characteristics of the system. In this paper, only lateral-directional axes are considered since the lateral-directional control law was specifically designed to produce classical airplane responses normally expected with low-order, rigid-body systems. Evaluation of the control design methodology was made using low-order equivalent systems determined from flight and simulation. This allowed comparison of the closed-loop rigid-body dynamics achieved in flight with that designed in simulation. In flight, the On Board Excitation System was used to apply optimal inputs to lateral stick and pedals at five angles at attack : 5, 20, 30, 45, and 60 degrees. Data analysis and closed-loop model identification were done using frequency domain maximum likelihood. The structure of identified models was a linear state-space model reflecting classical 4th-order airplane dynamics. Input time delays associated with the high-order controller and aircraft system were accounted for in data preprocessing. A comparison of flight estimated models with small perturbation linear design models highlighted nonlinearities in the system and indicated that the closed-loop rigid-body dynamics were sensitive to input amplitudes at 20 and 30 degrees angle of attack.

Murphy, Patrick C.↗

Closed-Loop System Identification Experience for Flight Control Law and Flying Qualities Evaluation of a High Performance Fighter Aircraft

This paper highlights some of the results and issues associated with estimating models to evaluate control law design methods and design criteria for advanced high performance aircraft. Experimental fighter aircraft such as the NASA High Alpha Research Vehicle (HARV) have the capability to maneuver at very high angles of attack where nonlinear aerodynamics often predominate. HARV is an experimental F/A-18, configured with thrust vectoring and conformal actuated nose strakes. Identifying closed-loop models for this type of aircraft can be made difficult by nonlinearities and high-order characteristics of the system. In this paper only lateral-directional axes are considered since the lateral-directional control law was specifically designed to produce classical airplane responses normally expected with low-order, rigid-body systems. Evaluation of the control design methodology was made using low-order equivalent systems determined from flight and simulation. This allowed comparison of the closed-loop rigid-body dynamics achieved in flight with that designed in simulation. In flight, the On Board Excitation System was used to apply optimal inputs to lateral stick and pedals at five angles of attack: 5, 20, 30, 45, and 60 degrees. Data analysis and closed-loop model identification were done using frequency domain maximum likelihood. The structure of the identified models was a linear state-space model reflecting classical 4th-order airplane dynamics. Input time delays associated with the high-order controller and aircraft system were accounted for in data preprocessing. A comparison of flight estimated models with small perturbation linear design models highlighted nonlinearities in the system and indicated that the estimated closed-loop rigid-body dynamics were sensitive to input amplitudes at 20 and 30 degrees angle of attack.

Murphy, Patrick C.↗

Engineering of Novel Biocolloid Suspensions

Colloidal suspensions are materials with a variety of uses from cleaners and lubricants to food, cosmetics, and coatings. In addition, they can be used as a tool for testing the fundamental tenets of statistical physics. Colloidal suspensions can be synthesized from a wide variety of materials, and in the form of monodisperse particles, which can self-assemble into highly ordered colloidal crystal structures. As such they can also be used as templates for the construction of highly ordered materials. Materials design of colloids has, to date, relied on entropic self-assembly, where crystals form as result of lower free energy due to a transition to order. Here, our goal is to develop a completely new method for materials fabrication using colloidal precursors, in which the self-assembly of the ordered colloidal structures is driven by a highly controllable, attractive interaction. This will greatly increase the range of potential structures that can be fabricated with colloidal particles. In this work, we demonstrate that colloidal suspensions can be crosslinked through highly specific biological crosslinking reactions. In particular, the molecules we use are protein-carbohydrate interactions derived from the immune system. This different driving force for self-assembly will yield different and novel suspensions structures. Because the biological interactions are heterotypic (A binding to B), this chemical system can be used to make binary alloys in which the two colloid subpopulations vary in some property - size, density, volume fraction, magnetic susceptibility, etc. An additional feature of these molecules which is unique - even within the realm of biological recognition - is that the molecules bind reversibly on reasonable time-scales, which will enable the suspension to sample different configurations, and allow us to manipulate and measure the size of the suspension dynamically. Because of the wide variety of structures that can be made from these novel colloids, and because the suspension structure can be altered dynamically, we believe this biocolloid system will yield a novel set of materials with many technological applications, including sensors (both biological and non-biological), optical filters and separation media.

Hammer, D. A.↗

The construction of high-accuracy schemes for acoustic equations

An accuracy analysis of various high order schemes is performed from an interpolation point of view. The analysis indicates that classical high order finite difference schemes, which use polynomial interpolation, hold high accuracy only at nodes and are therefore not suitable for time-dependent problems. Thus, some schemes improve their numerical accuracy within grid cells by the near-minimax approximation method, but their practical significance is degraded by maintaining the same stencil as classical schemes. One-step methods in space discretization, which use piecewise polynomial interpolation and involve data at only two points, can generate a uniform accuracy over the whole grid cell and avoid spurious roots. As a result, they are more accurate and efficient than multistep methods. In particular, the Cubic-Interpolated Psuedoparticle (CIP) scheme is recommended for computational acoustics.

Tang, Lei↗

Preconditioning for first-order spectral discretization

Efficient solution of the equations from spectral discretizations is essential if the high-order accuracy of these methods is to be realized. Direct solution of these equations is rarely feasible, thus iterative techniques are required. A preconditioning scheme for first-order Chebyshev collocation operators is proposed herein, in which the central finite difference mesh is finer than the collocation mesh. Details of the proper techniques for transferring information between the meshes are given here, and the scheme is analyzed by examination of the eigenvalue spectra of the preconditioned operators. The effect of artificial viscosity required in the inversion of the finite difference operator is examined. A second preconditioning scheme, involving a high-order upwind finite difference operator of the van Leer type is also analyzed to provide a comparison with the present scheme. Finally, the performance of the present scheme is verified by application to several test problems.

Streett, C. L.↗