Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservative schemes”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗
NASA NTRS1984

Improved finite difference schemes for transonic potential calculations

Engquist and Osher (1980) have introduced a finite difference scheme for solving the transonic small disturbance equation, taking into account cases in which only compression shocks are admitted. Osher et al. (1983) studied a class of schemes for the full potential equation. It is proved that these schemes satisfy a new discrete 'entropy inequality' which rules out expansion shocks. However, the conducted analysis is restricted to steady two-dimensional flows. The present investigation is concerned with the adoption of a heuristic approach. The full potential equation in conservation form is solved with the aid of a modified artificial density method, based on flux biasing. It is shown that, with the current scheme, expansion shocks are not possible.

Hafez, M.↗
NASA NTRS1992

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗
DOE OSTI2021

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. Here, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗
DOE OSTI2021

Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa

In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

A comparison of several computational techniques for solving some common aeronomic problems

Several numerical integration techniques for solving common aeronomic problems involving species rate equations are compared for speed and accuracy. A newer technique that defines families of species that are nearly conserved is found to be superior to an iterative technique when both methods are applied to simple test problems. The 'conservation' technique is also found to be more economical than the more complex Gear (1969) integration scheme for comparable accuracy.

Turco, R. P.↗
NASA NTRS1983

Application of second-order-accurate Total Variation Diminishing (TVD) schemes to the Euler equations in general geometries

A one-parameter family of explicit and implicit second-order-accurate, entropy satisfying, total variation diminishing (TVD) schemes was developed by Harten. These TVD schemes were the property of not generating spurious oscillations for one-dimensional nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. Application of these methods to one- and two-dimensional fluid flows containing shocks (in Cartesian coordinates) yields highly accurate nonoscillatory numerical solutions. The goal of this work is to expand these methods to the multidimensional Euler equations in generalized coordinate systems. Some numerical results of shock waves impinging on cylindrical bodies are compared with MacCormack's method.

Yee, H. C.↗
NASA NTRS1986

Analysis of Scramjet Inlets

NASCRIN analyzes two-dimensional flow fields in supersonic combustion ramjet (scramjet) inlets. Solves two-dimensional Euler or Navier-Strokes equations in conservative form by unsplit, explicit, two-step finitedifference method. More recent explicit/implicit, two-step scheme also incorporated for viscous flow analysis. Algebraic, two-layer eddy-viscosity model used for turbulent flow calculations.

Kumar, A.↗
NASA NTRS1984

Finite difference simulation of unsteady interactive flows

The development of time-dependent numerical simulations of unsteady interactive flows of an aerodynamic nature is reviewed with emphasis on compressible flows at flight Reynolds numbers and noniterative schemes based on Navier-Stokes equations. The importance of writing the equations in strong conservation-law form for a generalized body-oriented coordinate system is pointed out. The discussion covers time and length scales and numerical methods currently in use. Some computed results are presented and compared with experimental data.

Deiwert, G. S.↗
NASA NTRS1988

A longitudinally-patched grid approach with applications to high speed flows

The combined use of an implicit, upwind finite-volume scheme with an 'equivalent' gamma formulation (for real gas effects) and with a three-dimensional conservative patched grid formulation is discussed. Results are presented over a wide Mach number range on both single and patched grids.

Walters, R. W.↗
NASA NTRS1991

The effect of bleed configuration on shock/boundary layer interactions

The flow field characteristics are simulated numerically in an oblique shock wave/laminar boundary layer interaction for three different bleed configurations. The numerical solution for the flow field is obtained for the strong conservation-law form of the two-dimensional compressible Navier-Stokes equations using an implicit scheme. The computations model the flow in the interaction region and inside the bleed slot for an impinging oblique shock of sufficient strength to cause boundary layer separation on a flat plate boundary layer. The computed results are presented for a normal bleed slot at three different locations, upstream, across and downstream of the shock impingement point. The detailed flow characteristics in the interaction zone and inside the bleed slot are compared for the three cases. The resulting surface pressure and shear stress distributions are also compared for the three bleed cases and for the shock wave boundary layer interaction without bleed.

Hamed, A.↗
NASA NTRS1992

A non-linearly stable implicit finite element algorithm for hypersonic aerodynamics

A generalized curvilinear coordinate Taylor weak statement implicit finite element algorithm is developed for the two-dimensional and axisymmetric compressible Navier-Stokes equations for ideal and reacting gases. For accurate hypersonic simulation, air is modeled as a mixture of five perfect gases, i.e., molecular and atomic oxygen and nitrogen as well as nitric oxide. The associated pressure is then determined via Newton solution of the classical chemical equilibrium equation system. The directional semidiscretization is achieved using an optimal metric data Galerkin finite element weak statement, on a developed 'companion conservation law system', permitting classical test and trial space definitions. Utilizing an implicit Runge-Kutta scheme, the terminal algorithm is then nonlinearly stable, and second-order accurate in space and time on arbitrary curvilinear coordinates. Subsequently, a matrix tensor product factorization procedure permits an efficient numerical linear algebra handling for large Courant numbers. For ideal- and real-gas hypersonic flows, the algorithm generates essentially nonoscillatory numerical solutions in the presence of strong detached shocks and boundary layer-inviscid flow interactions.

Iannelli, G. S.↗
NASA NTRS1992

Unsteady transonic Euler solutions using finite elements

A finite element solution of the unsteady Euler equations is presented and demonstrated for 2D airfoil configurations oscillating in transonic flows. Computations are performed by spatially discretizing the conservation equations using the Galerkin weighted residual method and then employing a multistage Runge-Kutta scheme to march forward in time. A mesh deformation scheme has been developed to efficiently move interior points in a smooth fashion as the airfoil undergoes rigid body pitch and plunge motion. Both steady and unsteady results are presented, and a comparison is made with solutions obtained using finite-volume techniques. The effects of using either a lumped or consistent mass matrix are presented; the finite element method provides an accurate solution for unsteady transonic flows about isolated airfoils.

Davis, Gary A.↗
NASA NTRS1994

Lox droplet vaporization in a supercritical forced convective environment

A systematic investigation has been conducted to study the effects of ambient flow conditions (i.e. pressure and velocity) on supercritical droplet gasification in a forced-convective environment. The model is based on the time-dependent conservation equations in axisymmetric coordinates, and accommodates thermodynamic nonidealities and transport anomalies. In addition, an efficient scheme for evaluating thermophysical properties over the entire range of fluid thermodynamic states is established. The analysis allows a thorough examination of droplet behavior during its entire lifetime, including transient gasification, dynamic deformation, and shattering. A parametric study of droplet vaporization rate in terms of ambient pressure and Reynolds number is also conducted.

Hsiao, Chia-Chun↗
NASA NTRS1994

The Flux-integral Method for Multidimensional Convection and Diffusion

The flux-integral method is a procedure for constructing an explicit, single-step, forward-in-time, conservative, control volume update of the unsteady, multidimensional convection-diffusion equation. The convective plus diffusive flux at each face of a control-volume cell is estimated by integrating the transported variable and its face-normal derivative over the volume swept out by the convecting velocity field. This yields a unique description of the fluxes, whereas other conservative methods rely on nonunique, arbitrary pseudoflux-difference splitting procedures. The accuracy of the resulting scheme depends on the form of the subcell interpolation assumed, given cell-average data. Cellwise constant behavior results in a (very artificially diffusive) first-order convection scheme. Second-order convection-diffusion schemes correspond to cellwise linear (or bilinear) subcell interpolation. Cellwise quadratic subcell interpolants generate a highly accurate convection-diffusion scheme with excellent phase accuracy. Under constant-coefficient conditions, this is a uniformly third-order polynomial interpolation algorithm (UTOPIA).

Leonard, B. P.↗
NASA NTRS1995

Results of two new methods for aeroacoustics benchmark problems

Two new methods for the numerical solution of conservation laws (the Euler equations in particular) are presented: a uniformly second-order accurate upwind scheme and a third-order accurate centered scheme. Results of these schemes are shown for problems 1, 2, and 5 of this workshop's benchmark problems.

Huynh, H. T.↗