Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-order numerics”

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 37 records · Page 2

High-order state space simulation models of helicopter flight mechanics

This paper describes the formulation and validation of a high-order linearized mathematical model of helicopter flight mechanics, which includes rotor flap and lag degrees of freedom as well as inflow dynamics. The model is extracted numerically from an existing nonlinear, blade element, time simulation model. Extensive modifications in the formulation and solution process of the nonlinear model, required for a theoretically rigorous linearization, are described in detail. The validation results show that the linearized model successfully captures the coupled rotor-fuselage dynamics in the frequency band most critical for the design of advanced flight control systems. Additional results quantify the extent to which the order of the model can be reduced without loss of fidelity.

Kim, Frederick D.↗

Discussion of the potential and limitations of direct and large-eddy simulations

The full text of the discussion paper presented at the Whither Turbulence Workshop on the potential and limitations of direct and large-eddy simulations is provided. Particular emphasis is placed on discussing the role of numerics and mathematical theory in direct simulations of both compressible and incompressible flows. A variety of unresolved issues with large-eddy simulations such as their implementation in high-order finite difference codes, problems with defiltering, and modifications to accommodate integrations to solid boundaries are elaborated on. These as well as other points are discussed in detail along with the authors' views concerning the prospects for future research.

Hussaini, M. Y.↗

Spectral element methods for the incompressible Navier-Stokes equations

Spectral element methods are high-order weighted-residual techniques for partial differential equations that combine the geometric flexibility of finite element techniques with the rapid convergence rate of spectral schemes. The theoretical foundations and numerical implementation of spectral element methods for the incompressible Navier-Stokes equations are presented, considering the construction and analysis of optimal-order spectral element discretizations for elliptic and saddle (Stokes) problems, as well as the efficient solution of the resulting discrete equations by rapidly convergent tensor-product-based iterative procedures. Several examples of spectral element simulation of moderate Reynolds number unsteady flow in complex geometry are presented.

Maday, Yvon↗

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster↗

Comparison of numerical integration techniques for orbital applications

The present work gives a brief comparison of the performance of programs for integrating differential equations for orbital applications. The evaluation criteria and the method of testing are described, and the results of the test problem set are included. Integration methods that were chosen for comparison include high-order Runge-Kutta methods; a rational extrapolation method (Bulirsch and Stoer, 1966); a variable step, variable order, multistep method (Krogh, 1969); classical multistep methods of Adams and Cowell; and modified multistep methods. The high-order Runge-Kutta methods used in the comparison include RKF 7(8) and RKF 8(9) (Fehlberg, 1968), and RKS 8-10 (Shanks, 1966).

Moore, H.↗

A multi-domain spectral computation of three-dimensional laminar horseshoe vortex flow using incompressible Navier-Stokes equations

Multidomain spectral methods are presently used to numerically simulate a strut-wall intersection's laminar horseshoe vortex flow through direct solution of the three-dimensional, incompressible, time-dependent Navier-Stokes equations. Direct expansion in Chebyshev polynomials and spectral element method spatial discretization of flow dependence are used to achieve high-order accuracy, and minimal dispersion errors. Low and moderate Reynolds number results are presented to illustrate the method application.

Tan, C. S.↗

Numerical investigation of chemically reacting flows in ramjet dump combustors

The time-dependent Navier-Stokes equations, including second-order turbulence model, are numerically integrated by using four-stage Runge-Kutta scheme to predict the steady-state supersonic flow structures in ramjet dump combustors. The formulation is derived for reacting flows with finite-rate chemistry. In the present study, it is firstly attempted to assess the accuracy of existing high-order turbulence model in supersonic flows. The comparison shows reasonable agreement between calculated and measured data in terms of velocity distributions. It is indicated that a modified constant C-mu for calculating turbulent eddy viscosity is needed in the supersonic flow regime and the adaptive meshing is preferred to capture the recirculation zone. In the reacting flow calculation, the results from a test case of hydrogen and air combustion at premixed conditon show that the rearward facing step is able to increase flow residence time and stabilize the flame in supersonic flows.

Hsieh, Kwang-Chung↗

Utility of a finite element solution algorithm for initial-value problems

The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.

Baker, A. J.↗

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗

Dynamics of Pluto

The present study of the Pluto orbit yielded by the LONGSTOP 1B, 100-Myr numerical integration of the outer planets has given attention to the 3:2 resonance in mean motion with Neptune up to degree 2 in eccentricities and inclinations. Confirmations are obtained for both the 19,900-year period longitudinal libration and the 3.78-Myr period libration of the Pluto pericenter argument. Because chaos originates with small divisors, possible resonances have been searched for among the frequencies associated with the critical arguments. The macroscopic stability of the orbit of Pluto is explainable in terms of the association of high-order resonances with small chaotic regions.

Milani, A.↗

Minimum-fuel control of high-order systems

The minimum-fuel control problem is of special interest in various space systems. To date, solutions of minimum-fuel control problems have been carried out for relatively low-order systems. Space structures, however, are generally characterized by a large number of degrees of freedom, so that minimum-fuel control of such systems requires a new approach. In the independent modal-space control (IMSC) method, the control laws are designed in the modal space for each mode independently. The minimum-fuel problem reduces to that of a set of independent second-order systems, so that minimum-fuel control is possible. This paper shows how the IMSC method can be used to control a space structure with a minimum amount of fuel. A numerical example is presented.

Shenhar, J.↗

Finite-dimensional approximation for optimal fixed-order compensation of distributed parameter systems

In controlling distributed parameter systems it is often desirable to obtain low-order, finite-dimensional controllers in order to minimize real-time computational requirements. Standard approaches to this problem employ model/controller reduction techniques in conjunction with LQG theory. In this paper we consider the finite-dimensional approximation of the infinite-dimensional Bernstein/Hyland optimal projection theory. This approach yields fixed-finite-order controllers which are optimal with respect to high-order, approximating, finite-dimensional plant models. The technique is illustrated by computing a sequence of first-order controllers for one-dimensional, single-input/single-output, parabolic (heat/diffusion) and hereditary systems using spline-based, Ritz-Galerkin, finite element approximation. Numerical studies indicate convergence of the feedback gains with less than 2 percent performance degradation over full-order LQG controllers for the parabolic system and 10 percent degradation for the hereditary system.

Bernstein, Dennis S.↗

High-order low-pass implicit tangent filters for use in finite area calculations

High-order implicit tangent filters are developed. The implicit tangent filters possess a highly selective amplitude response function and they can be applied relatively close to a boundary. Comparisons are made between the implicit and the traditional or explicit filters. Numerical simulations are performed to test the response in a limited-area model. An algorithm to maximize computing efficiency is presented. All tests indicate the desirability and utility of the implicit tangent filter.

Raymond, William H.↗

An approximation technique for computing optimal fixed-order controllers for infinite-dimensional systems

The finite-dimensional approximation of the infinite-dimensional Bernstein/Hyland optimal projection theory is investigated analytically. The approach yields fixed-finite-order controllers which are optimal with respect to high-order approximating finite-dimensional plant models. The technique is illustrated by computing a sequence of first-order controllers for a one-dimensional SISO parabolic (heat/diffusion) system using a spline-based Ritz-Galerkin finite-element approximation. The numerical studies indicate convergence of the feedback gains with less than 2-percent performance degradation over full-order LQG controllers.

Bernstein, Dennis S.↗

A time accurate finite volume high resolution scheme for three dimensional Navier-Stokes equations

A time accurate, three-dimensional, finite volume, high resolution scheme for solving the compressible full Navier-Stokes equations is presented. The present derivation is based on the upwind split formulas, specifically with the application of Roe's (1981) flux difference splitting. A high-order accurate (up to the third order) upwind interpolation formula for the inviscid terms is derived to account for nonuniform meshes. For the viscous terms, discretizations consistent with the finite volume concept are described. A variant of second-order time accurate method is proposed that utilizes identical procedures in both the predictor and corrector steps. Avoiding the definition of midpoint gives a consistent and easy procedure, in the framework of finite volume discretization, for treating viscous transport terms in the curvilinear coordinates. For the boundary cells, a new treatment is introduced that not only avoids the use of 'ghost cells' and the associated problems, but also satisfies the tangency conditions exactly and allows easy definition of viscous transport terms at the first interface next to the boundary cells. Numerical tests of steady and unsteady high speed flows show that the present scheme gives accurate solutions.

Liou, Meng-Sing↗

A high-order Lagrangian-decoupling method for the incompressible Navier-Stokes equations

A high-order Lagrangian-decoupling method is presented for the unsteady convection-diffusion and incompressible Navier-Stokes equations. The method is based upon: (1) Lagrangian variational forms that reduce the convection-diffusion equation to a symmetric initial value problem; (2) implicit high-order backward-differentiation finite-difference schemes for integration along characteristics; (3) finite element or spectral element spatial discretizations; and (4) mesh-invariance procedures and high-order explicit time-stepping schemes for deducing function values at convected space-time points. The method improves upon previous finite element characteristic methods through the systematic and efficient extension to high order accuracy, and the introduction of a simple structure-preserving characteristic-foot calculation procedure which is readily implemented on modern architectures. The new method is significantly more efficient than explicit-convection schemes for the Navier-Stokes equations due to the decoupling of the convection and Stokes operators and the attendant increase in temporal stability. Numerous numerical examples are given for the convection-diffusion and Navier-Stokes equations for the particular case of a spectral element spatial discretization.

Ho, Lee-Wing↗

Understanding transition and turbulence through direct simulations

Direct simulations consist in solving the full Navier-Stokes equations, without any turbulence model, and describing all the detailed features of the flow. Usually the flows are three-dimensional and time-dependent and contain both coarse and fine structures, which makes the numerical task very challenging in terms of both the algorithm and the computational effort. Most of the work until now has involved spectral methods, which are highly accurate but not very flexible in terms of geometry or complex equations. For that reason, future work will also rely on high-order finite-difference or other methods. Direct simulations complement experimental work, and both contribute to the theory and the empirical knowledge of turbulence. Once such a simulation has been shown to be accurate, the flow field is completely known in three dimensions and time, including the pressure, the vorticity and any other quantity. On the other hand, most simulations to date solved the incompressible equations in rather simple geometries, and direct simulations will always be limited to moderate Reynolds numbers. Extensive simulations have been conducted in homogeneous turbulence, channel flows, boundary layers, and mixing layers. Much effort is devoted to addressing flows with compressibility and chemical reactions, and to new geometries such as a backward-facing step.

Spalart, P. R.↗

Numerical experiments on the accuracy of ENO and modified ENO schemes

Numerical experiments have been performed using different ENO schemes and different time discretizations, in order to assess accuracy-degeneracy phenomena of the type described by Rogerson and Meiburg (1990). A modified ENO scheme recovers the correct order of accuracy for all the test problems with smooth initial conditions and gives results comparable to the original ENO schemes for discontinuous problems. It is concluded that ENO schemes may lose the full high-order accuracy predicted by local truncation-error analysis for scalar linear conservation laws with smooth initial conditions. A modified ENO scheme can overcome this accuracy-degeneracy problem for the test problems without increasing the computational cost.

Shu, Chi-Wang↗