SEARCH · Engineering Papers

Results for “conservative scheme”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1976

An implicit finite-difference algorithm for hyperbolic systems in conservation-law form

An implicit finite-difference scheme is developed for the efficient numerical solution of nonlinear hyperbolic systems in conservation-law form. The algorithm is second-order time-accurate, noniterative, and in a spatially factored form. Second- or fourth-order central and second-order one-sided spatial differencing are accommodated within the solution of a block tridiagonal system of equations. Significant conceptual and computational simplifications are made for systems whose flux vectors are homogeneous functions (of degree one), e.g., the Eulerian gasdynamic equations. Conservative hybrid schemes, which switch from central to one-sided spatial differencing whenever the local characteristic speeds are of the same sign, are constructed to improve the resolution of weak solutions. Numerical solutions are presented for a nonlinear scalar model equation and the two-dimensional Eulerian gasdynamic equations.

Beam, R. M.↗

NASA NTRS1978

Recent experiences with three-dimensional transonic potential flow calculations

Some recent experiences with computer programs capable of solving finitie-difference approximations to the full potential equation for the transonic flow past three dimensional swept wings and simple wing-fuselage combinations are discussed. The programs used are a nonconservative program for swept wings, a quasi-conservative finite-volume program capable of treating swept wings mounted on fuselages of slowly varying circular cross section, and a fully conservative finite volume scheme capable of treating swept wings and wing-cylinder combinations. The present capabilities of these codes are reviewed. The relative merits of the conservative and nonconservative formulations are discussed, and the results of calculations including corrections for the boundary-layer displacement effect are presented.

Caughey, D. A.↗

NASA NTRS1988

Total Variation Diminishing (TVD) schemes of uniform accuracy

Explicit second-order accurate finite-difference schemes for the approximation of hyperbolic conservation laws are presented. These schemes are nonlinear even for the constant coefficient case. They are based on first-order upwind schemes. Their accuracy is enhanced by locally replacing the first-order one-sided differences with either second-order one-sided differences or central differences or a blend thereof. The appropriate local difference stencils are selected such that they give TVD schemes of uniform second-order accuracy in the scalar, or linear systems, case. Like conventional TVD schemes, the new schemes avoid a Gibbs phenomenon at discontinuities of the solution, but they do not switch back to first-order accuracy, in the sense of truncation error, at extrema of the solution. The performance of the new schemes is demonstrated in several numerical tests.

Hartwich, PETER-M.↗

NASA NTRS1995

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Yee, H. C.↗

NASA NTRS1986

Linearized form of implicit TVD schemes for the multidimensional Euler and Navier-Stokes equations

Linearized alternating direction implicit (ADI) forms of a class of total variation diminishing (TVD) schemes for the Euler and Navier-Stokes equations have been developed. These schemes are based on the second-order-accurate TVD schemes for hyperbolic conservation laws developed by Harten (1983, 1984). They have the property of not generating spurious oscillations across shocks and contact discontinuities. In general, shocks can be captured within 1-2 grid points. These schemes are relatively simple to understand and easy to implement into a new or existing computer code. One can modify a standard three-point central-difference code by simply changing the conventional numerical dissipation term into the one designed for the TVD scheme. For steady-state applications, the only difference in computation is that the current schemes require a more elaborate dissipation term for the explicit operator; no extra computation is required for the implicit operator. Numerical experiments with the proposed algorithms on a variety of steady-state airfoil problems illustrate the versatility of the schemes.

Yee, H. C.↗

NASA NTRS1985

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The novel implicit and unconditionally stable, high resolution Total Variation Diminishing (TVD) scheme whose application to steady state calculations is presently examined is a member of a one-parameter family of implicit, second-order accurate systems developed by Harten (1983) for the computation of weak solutions for one-dimensional hyperbolic conservation laws. The scheme will not generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments for a quasi-one-dimensional nozzle problem show that the experimentally determined stability limit correlates exactly with the theoretical stability limit for the nonlinear scalar hyberbolic conservation laws.

Yee, H. C.↗

NASA NTRS1975

Transonic potential flow calculations using conservation form

A method is presented for the solution of the full potential equation in conservation form. This assures that a weak solution satisfying proper isentropic jump conditions is obtained, giving an improved representation of shock waves in comparison with earlier nonconservative schemes. The method uses the concept of artificial viscosity to produce a stable difference scheme in the supersonic zone, and artificial time to generate a convergent iterative scheme for the solution of the difference equations. Results of calculations using the nonconservative and conservative schemes are compared.

Jameson, A.↗

NASA NTRS2018

High-Order Stabilized Finite Elements on Dynamic Meshes

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.

Anderson, W. Kyle↗

NASA NTRS2013

Studies of Inviscid Flux Schemes for Acoustics and Turbulence Problems

Five different central difference schemes, based on a conservative differencing form of the Kennedy and Gruber skew-symmetric scheme, were compared with six different upwind schemes based on primitive variable reconstruction and the Roe flux. These eleven schemes were tested on a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem and a turbulent channel flow problem. The central schemes were generally very accurate and stable, provided the grid stretching rate was kept below 10%. As near-DNS grid resolutions, the results were comparable to reference DNS calculations. At coarser grid resolutions, the need for an LES SGS model became apparent. There was a noticeable improvement moving from CD-2 to CD-4, and higher-order schemes appear to yield clear benefits on coarser grids. The UB-7 and CU-5 upwind schemes also performed very well at near-DNS grid resolutions. The UB-5 upwind scheme does not do as well, but does appear to be suitable for well-resolved DNS. The UF-2 and UB-3 upwind schemes, which have significant dissipation over a wide spectral range, appear to be poorly suited for DNS or LES.

Morris, Chris↗

NASA NTRS1984

A two-point difference scheme for computing steady-state solutions to the conservative one-dimensional Euler equations

An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central-difference scheme for the flux derivatives with dissipation added at supersonic points via the retarded density concept. Application of the method to 1-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. For 1-dimensional Euler calculations, it is shown that the scheme offers two advantages over the more widely-used three-point schemes. The first is in regard to application of boundary conditions, and the second relates to the fact that the two-point algorithm is well-conditioned for large time steps.

Wornom, S. F.↗

NASA NTRS2003

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows, I: Basic Theory

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls t o limit the amount of and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvilinear grids. The four advantages of the present approach over existing MHD schemes reported in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion MHD flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike exist- ing schemes for the conservative MHD equations which suffer from ill-conditioned eigen- decompositions, the present scheme makes use of a well-conditioned eigen-decomposition obtained from a minor modification of the eigenvectors of the non-conservative MHD equations t o solve the conservative form of the MHD equations. Third, this approach of using the non-conservative eigensystem when solving the conservative equations also works well in the context of standard shock-capturing schemes for the MHD equations. Fourth, a new approach to minimize the numerical error of the divergence-free magnetic condition for high order schemes is introduced. Numerical experiments with typical MHD model problems revealed the applicability of the newly developed schemes for the MHD equations.

Sjoegreen, Bjoern↗

NASA NTRS2002

Extension of Low Dissipative High Order Hydrodynamics Schemes for MHD Equations

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD (magnetohydrodynamic) equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls to limit the amount and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvi-linear grids. The three features of the present MHD scheme over existing schemes in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion magnetized flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike existing schemes for the conservative MHD equations which suffer from ill-conditioned eigen-decompositions, the present scheme makes use of a well-conditioned eigen-decomposition to solve the conservative form of the MHD equations. This is due to, partly. the fact that the divergence of the magnetic field condition is a different type of constraint from its incompressible Navier-Stokes cousin. Third, a new approach to minimize the numerical error of the divergence free magnetic condition for high order scheme is introduced.

Yee, H. C.↗

NASA NTRS1980

Calculation of three-dimensional unsteady transonic flows past helicopter blades

A finite difference code for predicting the high speed flow over the advancing helicopter rotor is presented. The code solves the low frequency, transonic small disturbance equation and is suitable for modeling the effects of advancing blade unsteadiness on blades of nearly arbitrary planform. The method employs a quasi-conservative mixed differencing scheme and solves the resulting difference equations by an alternating direction scheme. Computed results showed good agreement with experimental blade pressure data and illustrate some of the effects of varying the rotor planform. The flow unsteadiness is shown to be an indispensible part of a transonic solution. Close to the tip at high advance ratio, cross flow effects can significantly affect the solution.

Chattot, J. J.↗

NASA NTRS1983

Computational methods for ideal compressible flow

Conservative dissipative difference schemes, the recognition and representation of discontinuities, and multidimensional methods are discussed.

Vanleer, B.↗

NASA NTRS1987

Analysis of viscous transonic flow over airfoil sections

A full Navier-Stokes solver has been used to model transonic flow over three airfoil sections. The method uses a two-dimensional, implicit, conservative finite difference scheme for solving the compressible Navier-Stokes equations. Results are presented as prescribed for the Viscous Transonic Airfoil Workshop to be held at the AIAA 25th Aerospace Sciences Meeting. The NACA 0012, RAE 2822 and Jones airfoils have been investigated for both attached and separated transonic flows. Predictions for pressure distributions, loads, skin friction coefficients, boundary layer displacement thickness and velocity profiles are included and compared with experimental data when possible. Overall, the results are in good agreement with experimental data.

Huff, Dennis L.↗

NASA NTRS1987

Analysis of viscous transonic flow over airfoil sections

Huff, Dennis L.↗

NASA NTRS1988

Computational unsteady aerodynamics for lifting surfaces

Two dimensional problems are solved using numerical techniques. Navier-Stokes equations are studied both in the vorticity-stream function formulation which appears to be the optimal choice for two dimensional problems, using a storage approach, and in the velocity pressure formulation which minimizes the number of unknowns in three dimensional problems. Analysis shows that compact centered conservative second order schemes for the vorticity equation are the most robust for high Reynolds number flows. Serious difficulties remain in the choice of turbulent models, to keep reasonable CPU efficiency.

Edwards, John W.↗

NASA NTRS1992

Crystallization, sublimation, and gas release in the interior of a porous comet nucleus

A numerical code is developed for evolutionary calculations of the thermal structure and composition of a porous comet nucleus made of water ice, in amorphous or crystalline form, other volatiles, dust, and gases trapped in amorphous ice. Bulk evaporation, crystallization, gas release, and free (Knudsen) flow of gases through the pores are taken into account. The numerical scheme yields exact conservation laws for mass and energy. The code is used to study the effect of bulk evaporation of ice in the interior of a comet nucleus during crystallization. It is found that evaporation controls the temperature distribution; the vapor prevents cooling of the crystallized layer of ice, by recondensation and release of latent heat. Thus high temperatures are maintained below the surface of the nucleus and down to depths of tens or hundreds of meters, even at large heliocentric distances, as long as crystallization goes on. Gas trapped in the ice and released during the phase transition flows both toward the interior and toward the surface and out of the nucleus. The progress of crystallization is largely determined by the contribution of gas fluxes to heat transfer.

Prialnik, Dina↗