SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1989

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↗

NASA NTRS1989

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↗

NASA NTRS1989

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.↗

NASA NTRS1991

High-order ENO schemes applied to two- and three-dimensional compressible flow

High order essentially non-oscillatory (ENO) finite difference schemes are applied to the 2-D and 3-D compressible Euler and Navier-Stokes equations. Practical issues, such as vectorization, efficiency of coding, cost comparison with other numerical methods, and accuracy degeneracy effects, are discussed. Numerical examples are provided which are representative of computational problems of current interest in transition and turbulence physics. These require both nonoscillatory shock capturing and high resolution for detailed structures in the smooth regions and demonstrate the advantage of ENO schemes.

Shu, Chi-Wang↗

NASA NTRS1990

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↗

NASA NTRS1992

A numerical resolution study of high order essentially non-oscillatory schemes applied to incompressible flow

High order essentially non-oscillatory (ENO) schemes, originally designed for compressible flow and in general for hyperbolic conservation laws, are applied to incompressible Euler and Navier-Stokes equations with periodic boundary conditions. The projection to divergence-free velocity fields is achieved by fourth order central differences through Fast Fourier Transforms (FFT) and a mild high-order filtering. The objective of this work is to assess the resolution of ENO schemes for large scale features of the flow when a coarse grid is used and small scale features of the flow, such as shears and roll-ups, are not fully resolved. It is found that high-order ENO schemes remain stable under such situations and quantities related to large-scale features, such as the total circulation around the roll-up region, are adequately resolved.

Weinan, E.↗

NASA NTRS1991

Forward flight trim calculation and frequency response validation of a high-order helicopter simulation model

This paper describes a new trim procedure, that includes the calculation of the steady-state response of the rotor blades, and that is applicable to straight flight and steady coordinated turns. This paper also describes the results of a validation study for a high order linearized model of helicopter flight dynamics, that includes rotor, inflow, and actuator dynamics. The model is obtained by numerical perturbations of a nonlinear, blade element type mathematical model. Predicted responses are compared with flight test data for two values of flight speed. The comparison is carried out in the frequency domain. Numerical simulations show that the trim algorithm is very accurate, and preserves the periodicity of the aircraft states. The results also indicate that the predictions of the linearized model are in good agreement with flight test data, especially at medium and high frequencies.

Kim, Frederick D.↗

NASA NTRS1991

Controls for space structures

Assembly and operation of large space structures (LSS) in orbit will require robot-assisted docking and berthing of partially-assembled structures. These operations require new solutions to the problems of controls. This is true because of large transient and persistent disturbances, controller-structure interaction with unmodeled modes, poorly known structure parameters, slow actuator/sensor dynamical behavior, and excitation of nonlinear structure vibrations during control and assembly. For on-orbit assembly, controllers must start with finite element models of LSS and adapt on line to the best operating points, without compromising stability. This is not easy to do, since there are often unmodeled dynamic interactions between the controller and the structure. The indirect adaptive controllers are based on parameter estimation. Due to the large number of modes in LSS, this approach leads to very high-order control schemes with consequent poor stability and performance. In contrast, direct model reference adaptive controllers operate to force the LSS to track the desirable behavior of a chosen model. These schemes produce simple control algorithms which are easy to implement on line. One problem with their use for LSS has been that the model must be the same dimension as the LSS - i.e., quite large. A control theory based on the command generator tracker (CGT) ideas of Sobel, Mabins, Kaufman and Wen, Balas to obtain very low-order models based on adaptive algorithms was developed. Closed-loop stability for both finite element models and distributed parameter models of LSS was proved. In addition, successful numerical simulations on several LSS databases were obtained. An adaptive controller based on our theory was also implemented on a flexible robotic manipulator at Martin Marietta Astronautics. Computation schemes for controller-structure interaction with unmodeled modes, the residual mode filters or RMF, were developed. The RMF theory was modified to compensate slow actuator/sensor dynamics. These new ideas are being applied to LSS simulations to demonstrate the ease with which one can incorporate slow actuator/sensor effects into our design. It was also shown that residual mode filter compensation can be modified for small nonlinearities to produce exponentially stable closed-loop control. A theory for disturbance accommodating controllers based on reduced order models of structures was developed, and stability results for these controllers in closed-loop with large-scale finite element models of structures were obtained.

Balas, Mark↗

NASA NTRS1992

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.

Pruett, C. D.↗

NASA NTRS1993

Direct calculations of waves in fluid flows using a high-order compact difference scheme

The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expressions for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations, and the numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the CFL limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about 6 to 8 grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 to 300 time steps. Three numerical examples of waves in compressible flow are included.

Yu, Sheng-Tao↗

NASA NTRS1993

Numerical simulations of a high Mach number jet flow

The present 2D simulations of plane and axisymmetric jets, which are relevant to current efforts to suppress the jet exhaust noise of prospective high-speed civil transports, were conducted by solving full Navier-Stokes equations by means of a high-order finite-difference scheme. The results obtained, which were able to generate the correct mode shape after an adjustment region of about 10 diameters, are in good agreement with linear-theory predictions of the growth of instability waves.

Hayder, M. E.↗

NASA NTRS1993

Uniform high-order spectral methods for one- and two-dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one-dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two-dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗

NASA NTRS1988

Maximum likelihood identification for large space structures

This paper examines the use of on-orbit identification based on Maximum Likelihood Estimation (MLE) to provide these high-order, high-accuracy control design models for large space structures (LSS's). First, it outlines a general MLE identification algorithm, together with a covariance-analysis procedure to assess algorithm performance in terms of systematic and stochastic errors. Next, it examines various simplifications appropriate for the LSS identification application. Simplified analytical performance results are presented, as are numerical results to support these analyses. Finally, a graphical interpretation of these results is given.

Barrett, Michael F.↗

NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗

NASA NTRS1994

An HP Adaptive Discontinuous Galerkin Method for Hyperbolic Conservation Laws

This dissertation addresses various issues for model classes of hyperbolic conservation laws. The basic approach developed in this work employs a new family of adaptive, hp-version, finite element methods based on a special discontinuous Galerkin formulation for hyperbolic problems. The discontinuous Galerkin formulation admits high-order local approximations on domains of quite general geometry, while providing a natural framework for finite element approximations and for theoretical developments. The use of hp-versions of the finite element method makes possible exponentially convergent schemes with very high accuracies in certain cases; the use of adaptive hp-schemes allows h-refinement in regions of low regularity and p-enrichment to deliver high accuracy, while keeping problem sizes manageable and dramatically smaller than many conventional approaches. The use of discontinuous Galerkin methods is uncommon in applications, but the methods rest on a reasonable mathematical basis for low-order cases and has local approximation features that can be exploited to produce very efficient schemes, especially in a parallel, multiprocessor environment. The place of this work is to first and primarily focus on a model class of linear hyperbolic conservation laws for which concrete mathematical results, methodologies, error estimates, convergence criteria, and parallel adaptive strategies can be developed, and to then briefly explore some extensions to more general cases. Next, we provide preliminaries to the study and a review of some aspects of the theory of hyperbolic conservation laws. We also provide a review of relevant literature on this subject and on the numerical analysis of these types of problems.

Bey, Kim S.↗

NASA NTRS1994

High-order centered difference methods with sharp shock resolution

In this paper we consider high-order centered finite difference approximations of hyperbolic conservation laws. We propose different ways of adding artificial viscosity to obtain sharp shock resolution. For the Riemann problem we give simple explicit formulas for obtaining stationary one and two-point shocks. This can be done for any order of accuracy. It is shown that the addition of artificial viscosity is equivalent to ensuring the Lax k-shock condition. We also show numerical experiments that verify the theoretical results.

Gustafsson, Bertil↗

NASA NTRS1994

SCISEAL: A CFD code for analysis of fluid dynamic forces in seals

A viewgraph presentation is made of the objectives, capabilities, and test results of the computer code SCISEAL. Currently, the seal code has: a finite volume, pressure-based integration scheme; colocated variables with strong conservation approach; high-order spatial differencing, up to third-order; up to second-order temporal differencing; a comprehensive set of boundary conditions; a variety of turbulence models and surface roughness treatment; moving grid formulation for arbitrary rotor whirl; rotor dynamic coefficients calculated by the circular whirl and numerical shaker methods; and small perturbation capabilities to handle centered and eccentric seals.

Athavale, Mahesh↗

NASA NTRS1994

Low-dissipation and -dispersion Runge-Kutta schemes for computational acoustics

In this paper, we investigate accurate and efficient time advancing methods for computational acoustics, where non-dissipative and non-dispersive properties are of critical importance. Our analysis pertains to the application of Runge-Kutta methods to high-order finite difference discretization. In many CFD applications multi-stage Runge-Kutta schemes have often been favored for their low storage requirements and relatively large stability limits. For computing acoustic waves, however, the stability consideration alone is not sufficient, since the Runge-Kutta schemes entail both dissipation and dispersion errors. The time step is now limited by the tolerable dissipation and dispersion errors in the computation. In the present paper, it is shown that if the traditional Runge-Kutta schemes are used for time advancing in acoustic problems, time steps greatly smaller than that allowed by the stability limit are necessary. Low-Dissipation and -Dispersion Runge-Kutta (LDDRE) schemes are proposed, based on an optimization that minimizes the dissipation and dispersion errors for wave propagation. Order optimizations of both single-step and two-step alternating schemes are considered. The proposed LDDRK schemes are remarkably more efficient than the classical Runge-Kutta schemes for acoustic computations. Moreover, low storage implementations of the optimized schemes are discussed. Special issues of implementing numerical boundary conditions in the LDDRK schemes are also addressed.

Hu, F. Q.↗