Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-order meshes”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗
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 NTRS1994

Development of an adaptive hp-version finite element method for computational optimal control

In this research effort, the usefulness of hp-version finite elements and adaptive solution-refinement techniques in generating numerical solutions to optimal control problems has been investigated. Under NAG-939, a general FORTRAN code was developed which approximated solutions to optimal control problems with control constraints and state constraints. Within that methodology, to get high-order accuracy in solutions, the finite element mesh would have to be refined repeatedly through bisection of the entire mesh in a given phase. In the current research effort, the order of the shape functions in each element has been made a variable, giving more flexibility in error reduction and smoothing. Similarly, individual elements can each be subdivided into many pieces, depending on the local error indicator, while other parts of the mesh remain coarsely discretized. The problem remains to reduce and smooth the error while still keeping computational effort reasonable enough to calculate time histories in a short enough time for on-board applications.

Hodges, Dewey H.↗
NASA NTRS1998

Electromagnetics and Fluid Dynamics

Previous efforts focused on developing tools for design of low observables were sustained. The final product was the maturation of a high-order accurate finite-volume based code to solve Maxwell's equations. One of the primary achievements was the development and implementation of efficient filtering techniques which enhance the robustness of high-order and optimized schemes without significant adverse impact on accuracy. This has eliminated the stability barrier which restrains the common use of high-order schemes for conservative wave propagation phenomena on curvilinear meshes. A study was performed of crossing shock interactions under conditions of increasing interaction strength and asymmetry. In the first category, the observed computed topological bifurcations were correlated with the formation of various lines of coalescence and divergence evident in experimental and computed surf-ace oil maps. ne flow structure arising from asymmetric interactions was investigated with particular emphasis on: 1) vorticity dynamics, 2) shock-structure and 3) sidewall vortex loading. Several efforts of the prior year were successfully published in archival journals. The high-order algorithms developed for CEM have been implemented into the FDL3DI CFD code are presently undergoing extensive testing. Preliminary results are highly encouraging.

Gaitonde, Datta↗
NASA NTRS1999

Parallelization of an Object-Oriented Unstructured Aeroacoustics Solver

A computational aeroacoustics code based on the discontinuous Galerkin method is ported to several parallel platforms using MPI. The discontinuous Galerkin method is a compact high-order method that retains its accuracy and robustness on non-smooth unstructured meshes. In its semi-discrete form, the discontinuous Galerkin method can be combined with explicit time marching methods making it well suited to time accurate computations. The compact nature of the discontinuous Galerkin method also makes it well suited for distributed memory parallel platforms. The original serial code was written using an object-oriented approach and was previously optimized for cache-based machines. The port to parallel platforms was achieved simply by treating partition boundaries as a type of boundary condition. Code modifications were minimal because boundary conditions were abstractions in the original program. Scalability results are presented for the SCI Origin, IBM SP2, and clusters of SGI and Sun workstations. Slightly superlinear speedup is achieved on a fixed-size problem on the Origin, due to cache effects.

Baggag, Abdelkader↗
NASA NTRS1992

Theory and implementation of high-order adaptive hp methods for analysis of incompressible viscous flows

An account is given of 'smart' algorithms for CFD which change in structure and performance with during flow calculations to accommodate changing properties of the solution. Such algorithms prominently include adaptive FEM methods, which are designed to adjust mesh parameters for the control of numerical error; attention is presently given to those which change the mesh size h and the local spectral order p in order to achieve high accuracies with minimal numbers of degrees of freedom. These 'hp methods' produce exponentially convergent approximations through which flow features are resolved by automatically distributing element sizes and spectral orders. This leads to calculation of the local (elementwise) error in the approximation.

Oden, J. T.↗
NASA NTRS1993

A finite-volume high-order ENO scheme for two-dimensional hyperbolic systems

The finite-volume approach is presently used to obtain a 2D, high-order accurate and basically nonoscillatory shock-capture method whose high-order spatial accuracy is obtained by means of a piecewise polynomial approximation of the solution from cell averages. Attention is given to a high-order spatial operator that is able to both retain high-order accuracy in smooth regions and avoid the oscillations that are associated with interpolations across steep gradients. The operator is extended to hyperbolic systems of equations and curvilinear meshes.

Casper, Jay↗
NASA NTRS1999

A High-Order Method Using Unstructured Grids for the Aeroacoustic Analysis of Realistic Aircraft Configurations

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.

Atkins, Harold L.↗
NASA NTRS1986

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

RIACS

Topics considered include: high-performance computing; cognitive and perceptual prostheses (computational aids designed to leverage human abilities); autonomous systems. Also included: development of a 3D unstructured grid code based on a finite volume formulation and applied to the Navier-stokes equations; Cartesian grid methods for complex geometry; multigrid methods for solving elliptic problems on unstructured grids; algebraic non-overlapping domain decomposition methods for compressible fluid flow problems on unstructured meshes; numerical methods for the compressible navier-stokes equations with application to aerodynamic flows; research in aerodynamic shape optimization; S-HARP: a parallel dynamic spectral partitioner; numerical schemes for the Hamilton-Jacobi and level set equations on triangulated domains; application of high-order shock capturing schemes to direct simulation of turbulence; multicast technology; network testbeds; supercomputer consolidation project.

Oliger, Joseph↗
NASA NTRS1989

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↗
NASA NTRS1989

A spectral multi-domain technique applied to high-speed chemically reacting flows

The first applications of a spectral multidomain method for viscous compressible flow is presented. The method imposes a global flux balance condition at the interface so that high-order continuity of the solution is preserved. The global flux balance is imposed in terms of a spectral integral of the discrete equations across adjoining domains. Since the discretized equations interior to each domain solved are uncoupled from each other, and since the interface relation has a block structure, the solution scheme can be adapted to the particular requirements of each subdomain. The spectral multidomain technique presented is well-suited for the multiple scales associated with the chemically reacting and transition flows in hypersonic research. A nonstaggered multidomain discretization is used for the chemically reacting flow calculation, and the first implementation of a staggered multidomain mesh is presented for accurately solving the stability equation for a viscous compressible fluid.

Macaraeg, Michele G.↗
NASA NTRS1975

High-order numerical solutions using cubic splines

The cubic spline collocation procedure for the numerical solution of partial differential equations was reformulated so that the accuracy of the second-derivative approximation is improved and parallels that previously obtained for lower derivative terms. The final result is a numerical procedure having overall third-order accuracy for a nonuniform mesh and overall fourth-order accuracy for a uniform mesh. Application of the technique was made to the Burger's equation, to the flow around a linear corner, to the potential flow over a circular cylinder, and to boundary layer problems. The results confirmed the higher-order accuracy of the spline method and suggest that accurate solutions for more practical flow problems can be obtained with relatively coarse nonuniform meshes.

Rubin, S. G.↗
NASA NTRS1986

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

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially nonoscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.

Casper, J.↗
NASA NTRS1993

Some aspects of high-order numerical solutions of the linear convection equation with forced boundary conditions

A six-stage low-storage Runge-Kutta time-marching method is presented and shown to be an efficient method for use with high-accuracy spatial difference operators for wave propagation problems. The accuracy of the method for inhomogeneous ordinary differential equations is demonstrated through numerical solutions of the linear convection equation with forced boundary conditions. Numerical experiments are presented simulating a sine wave and a Gaussian pulse propagating into and through the domain. For practical levels of mesh refinement corresponding to roughly ten points per wavelength, the six-stage Runge-Kutta method is more accurate than the popular fourth-order Runge-Kutta method. Further numerical experiments are presented which show that the numerical boundary scheme at an inflow boundary can be a significant source of error when high-accuracy spatial discretizations are used.

Zingg, D. W.↗