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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Vortex methods for two- and three-dimensional flow simulations

The point vortex and vortex blob methods for two dimensional flows are presented. Several results are discussed concerning the numerical analysis of the latter scheme, e.g., the preservation of globally conserved quantities and the analysis of the spatial discretization error resulting from the convection of fixed blobs of vorticity. An application to the two dimensional mixing layer is briefly described. The contour dynamics method is also discussed. The simulation of three dimensional flows with vortex methods is discussed. A natural way to represent the vorticity is in the form of closed tubes of filaments of vorticity, although other schemes are examined. Applications to aircraft trailing vortices and to a turbulent spot in a laminar boundary layer are presented. Hybrid schemes that use an Eulerian mesh to solve the Poisson equation for the velocity field are discussed. The goal of these schemes is to avoid the high cost of the Biot-Savart integration if many vortex elements are used while enjoying most of the advantages of pure Lagrangian schemes.

Leonard, A.↗

NASA NTRS1992

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulent models. The probability density function (PDF) method offers an attractive alternative: in a PDF model, the chemical source terms are closed and do not require additional models. Because the number of computational operations grows only linearly in the Monte Carlo scheme, it is chosen over finite differencing schemes. A grid dependent Monte Carlo scheme following J.Y. Chen and W. Kollmann has been studied in the present work. It was found that in order to conserve the mass fractions absolutely, one needs to add further restrictions to the scheme, namely alpha(sub j) + gamma(sub j) = alpha(sub j - 1) + gamma(sub j + 1). A new algorithm was devised that satisfied this restriction in the case of pure diffusion or uniform flow problems. Using examples, it is shown that absolute conservation can be achieved. Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗

NASA NTRS1992

Two-dimensional hybrid continuum/particle approach for rarefied flows

A hybrid numerical technique previously developed for one-dimensional rarefied gas flows is generalized to two dimensions. The method is based on the fact that the flowfield that develops near a body in a rarefied gas typically contains local regions of continuum, transitional and free-molecular flow. By utilizing the solution technique most appropriate for each region and coupling the techniques in an interface region where both are applicable, more computationally efficient solutions can be obtained, or equivalently, more complex flowfields can be analyzed. The present method combines finite difference solution of the Navier-Stokes equations in the continuum regions, with Direct Simulation Monte Carlo in the more rarefied regions. The two schemes are coupled interactively via a general conservative flux boundary condition. The method is tested by application to the model problem of pressure-driven rarefied flow through a slit. Results show the hybrid scheme offers a speedup of a factor of nearly two for the nominal conditions considered, due to the decrease in the size of the Direct Simulation domain.

Wadsworth, Dean C.↗

NASA NTRS1983

Convergence of generalized MUSCL schemes

Semi-discrete generalizations of the second order extension of Godunov's scheme, known as the MUSCL scheme, are constructed, starting with any three point ""E" scheme. They are used to approximate scalar conservation laws in one space dimension. For convex conservation laws, each member of a wide class is proven to be a convergent approximation to the correct physical solution. Comparison with another class of high resolution convergent schemes is made.

Osher, S.↗

NASA NTRS1984

Convergence of generalized MUSCL schemes

Semi-discrete generalizations of the second order extension of Godunov's scheme, known as the MUSCL scheme, are constructed, starting with any three point E scheme. They are used to approximate scalar conservation laws in one space dimension. For convex conservation laws, each member of a wide class is proven to be a convergent approximation to the correct physical solution. Comparison with another class of high resolution convergent schemes is made.

Osher, S.↗

NASA NTRS1985

A 3-dimensional mass conserving element for compressible flows

A variety of finite element schemes has been used in the numerical approximation of compressible flows particularly in underwater acoustics. In many instances instabilities have been generated due to the lack of mass conservation. Two- and three-dimensional elements are developed which avoid these problems.

Fix, G.↗

NASA NTRS1995

Efficient implementation of weighted ENO schemes

In this paper, we further analyze, test, modify and improve the high order WENO (weighted essentially non-oscillatory) finite difference schemes of Liu, Osher and Chan. It was shown by Liu et al. that WENO schemes constructed from the r-th order (in L1 norm) ENO schemes are (r+1)-th order accurate. We propose a new way of measuring the smoothness of a numerical solution, emulating the idea of minimizing the total variation of the approximation, which results in a 5-th order WENO scheme for the case r = 3, instead of the 4-th order with the original smoothness measurement by Liu et al. This 5-th order WENO scheme is as fast as the 4-th order WENO scheme of Liu et al., and both schemes are about twice as fast as the 4-th order ENO schemes on vector supercomputers and as fast on serial and parallel computers. For Euler systems of gas dynamics, we suggest computing the weights from pressure and entropy instead of the characteristic values to simplify the costly characteristic procedure. The resulting WENO schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights, and work well for problems which do not contain strong shocks or strong reflected waves. We also prove that, for conservation laws with smooth solutions, all WENO schemes are convergent. Many numerical tests, including the 1D steady state nozzle flow problem and 2D shock entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the WENO scheme using the new smoothness measurement, in resolving complicated shock and flow structures. We have also applied Yang's artificial compression method to the WENO schemes to sharpen contact discontinuities.

Jiang, Guang-Shan↗

NASA NTRS1984

The large-time behavior of the scalar, genuinely nonlinear Lax-Friedrichs scheme

The Lax-Friedrichs scheme, approximating the scalar, genuinely nonlinear conservation law u sub t + f sub x (u) = O where f(u) is, say, strictly double dot f dot a sub asterisk O is studied. The divided differences of the numerical solution at time t do not exceed 2 (t dot a sub asterisk) to the -1. This one-sided Lipschitz boundedness is in complete agreement with the corresponding estimate one has in the differential case; in particular, it is independent of the initial amplitude in sharp contrast to linear problems. It guarantees the entropy compactness of the scheme in this case, as well as providing a quantitative insight into the large-time behavior of the numerical computation.

Tadmor, E.↗

NASA NTRS1983

The large-time behavior of the scalar, genuinely nonlinear Lax-Friedrichs scheme

The Lax-Friedrichs scheme, approximating the scalar, genuinely nonlinear conservation law u sub t + f sub x (u) = 0 where f(u) is, say, strictly convex double dot f dot a sub asterisk 0 is studied. The divided differences of the numerical solution at time t do not exceed 2 (t dot a sub asterisk) to the -1. This one-sided Lipschitz boundedness is in complete agreement with the corresponding estimate one has in the differential case; in particular, it is independent of the initial amplitude in sharp contrast to liner problems. It guarantees the entropy compactness of the scheme in this case, as well as providing a quantitive insight into the large-time behavior of the numerical computation.

Tadmor, E.↗

NASA NTRS1975

Implicit approximate-factorization schemes for the low-frequency transonic equation

Two- and three-level implicit finite-difference algorithms for the low-frequency transonic small disturbance-equation are constructed using approximate factorization techniques. The schemes are unconditionally stable for the model linear problem. For nonlinear mixed flows, the schemes maintain stability by the use of conservatively switched difference operators for which stability is maintained only if shock propagation is restricted to be less than one spatial grid point per time step. The shock-capturing properties of the schemes were studied for various shock motions that might be encountered in problems of engineering interest. Computed results for a model airfoil problem that produces a flow field similar to that about a helicopter rotor in forward flight show the development of a shock wave and its subsequent propagation upstream off the front of the airfoil.

Ballhaus, W. F.↗

NASA NTRS2005

Explicit Von Neumann Stability Conditions for the c-tau Scheme: A Basic Scheme in the Development of the CE-SE Courant Number Insensitive Schemes

As part of the continuous development of the space-time conservation element and solution element (CE-SE) method, recently a set of so call ed "Courant number insensitive schemes" has been proposed. The key advantage of these new schemes is that the numerical dissipation associa ted with them generally does not increase as the Courant number decre ases. As such, they can be applied to problems with large Courant number disparities (such as what commonly occurs in Navier-Stokes problem s) without incurring excessive numerical dissipation.

Chang, Sin-Chung↗

NASA NTRS2017

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗

NASA NTRS2014

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

NASA NTRS1990

Finite-element schemes for extended integrations of atmospheric models

Of the two finite-element models presently used to investigate the effect of the conservation of integral invariants, by means of finite-element discretization schemes of the shallow-water equations, in order to serve as a paradigm of long-term atmospheric-model integrations, the first employs rectangular elements and conserves total energy, while the second uses triangular elements and a high-accuracy two-stage Numerov-Galerkin method. Attention is given to critical times for numerical nonlinear instability, as well as to the determination of the critical degree of dissipation entailed by the achievement of stable long-term integrations. Relative computational efficiency and accuracy comparisons are presented for the two finite-element schemes.

Steppeler, J.↗

NASA NTRS1978

A finite element method for transmission in non-uniform ducts without flow Comparison with the method of weighted residuals

A finite element method is derived for studying the transmission of sound in nonuniform ducts without flow. The method is based on a weighted residuals technique and eight noded isoparametric elements. The method employs two computational schemes: the Helmholtz equation found through conservation equations, and the conservation equations themselves. Solutions are found for transmission and reflection coefficients for both two-dimensional and circular geometries. Numerical solutions are compared with those for a method of weighted residuals in the form of a modified Galerkin method, showing good agreement.

Astley, R. J.↗

NASA NTRS1983

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.

Phillips, T. N.↗

NASA NTRS1984

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Phillips, T. N.↗

NASA NTRS2002

A New Low Dissipative High Order Schemes for MHD Equations

The goal of this talk 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 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 curvilinear grids.

Yee, H. C.↗