Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent 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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Generation of a composite grid for turbine flows and consideration of a numerical scheme

A composite grid was generated for flows in turbines. It consisted of the C-grid (or O-grid) in the immediate vicinity of the blade and the H-grid in the middle of the blade passage between the C-grids and in the upstream region. This new composite grid provides better smoothness, resolution, and orthogonality than any single grid for a typical turbine blade with a large camber and rounded leading and trailing edges. The C-H (or O-H) composite grid has an unusual grid point that is connected to more than four neighboring nodes in two dimensions (more than six neighboring nodes in three dimensions). A finite-volume lower-upper (LU) implicit scheme to be used on this grid poses no problem and requires no special treatment because each interior cell of this composite grid has only four neighboring cells in two dimensions (six cells in three dimensions). The LU implicit scheme was demonstrated to be efficient and robust for external flows in a broad flow regime and can be easily applied to internal flows and extended from two to three dimensions.

Choo, Y.↗
NASA NTRS1975

A numerical scheme to solve unstable boundary value problems

A new iterative scheme for solving boundary value problems is presented. It consists of the introduction of an artificial time dependence into a modified version of the system of equations. Then explicit forward integrations in time are followed by explicit integrations backwards in time. The method converges under much more general conditions than schemes based in forward time integrations (false transient schemes). In particular it can attain a steady state solution of an elliptical system of equations even if the solution is unstable, in which case other iterative schemes fail to converge. The simplicity of its use makes it attractive for solving large systems of nonlinear equations.

Kalnay Derivas, E.↗
NASA NTRS2006

A Scheme for finding the Front Boundary of an Interplanetary Magnetic Cloud

We developed a scheme for finding the front boundary of an interplanetary magnetic cloud (MC) based on criteria that depend on the possible existence of any one or all of six specific solar wind features. The features that the program looks for, within +/- 2 hours of a preliminarily determined time for the front boundary, estimated either by visual inspection or by an automatic MC identification scheme, are: (1) a sufficiently large directional discontinuity in the interplanetary magnetic field (IMF), (2) existence of a magnetic hole, (3) a significant proton plasma beta drop, (4) a significant proton temperature drop, (5) a marked increase in the IMF's intensity, and (6) a significant decrease in a normalized root-mean-square deviation (RMS)of the magnetic field - where the scheme was tested using 5, 10, 15, and 20 minute averages of the relevant physical quantities, in order to find the optimum average (and RMS) to use. Other criteria, besides these six, were examined and dismissed as not reliable, e.g., plasma speed. The scheme was developed specifically for aiding in forecasting the strength and timing of a geomagnetic storm due to the passage of an interplanetary MC in real-time, but can be used in post ground-data collection for imposition of consistency in choosing a MC's front boundary. The scheme has been extensively tested, first using 80 bona fide MCs over about 9 years of WIND data, and also for 121 MC-like structures as defined by a program that automatically identifies such structures over the same period. Optimum limits for various parameters in the scheme were found by statistical studies of the WIND MCs. The resulting limits can be user-adjusted for other data sets, if desired. Final testing of the 80 MCs showed that for 50 percent of the events the boundary estimates occurred within +/-10 minutes of visually determined times, 80 percent occurred within +/-30 minutes, and 91 percent occur within +/-60 minutes, and three or more individual boundary tests were passed for 88 percent of the total MCs. The scheme and its testing will be described.

Lepping, Ronald P.↗
DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

A Unified Representation Scheme for Solid Geometric Objects Using B-splines (extended Abstract)

A geometric representation scheme called the B-spline cylinder, which consists of interpolation between pairs of uniform periodic cubic B-spline curves is discussed. This approach carries a number of interesting implications. For one, a single relatively simple database schema can be used to represent a reasonably large class of objects, since the spline representation is flexible enough to allow a large domain of representable objects at very little cost in data complexity. The model is thus very storage-efficient. A second feature of such a system is that it reduces to one the number of routines which the system must support to perform a given operation on objects. Third, the scheme enables easy conversion to and from other representations. The formal definition of the cylinder entity is given. In the geometric properties of the entity are explored and several operations on such objects are defined. Some general purpose criteria for evaluating any geometric representation scheme are introduced and the B-spline cylinder scheme according to these criteria is evaluated.

Bahler, D.↗
NASA NTRS2000

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Convection-Diffusion Problems

In the space-time conservation element and solution element (CE/SE) method, the independent marching variables used comprise not only the mesh value of the physical dependent variables but also, in contrast to it typical numerical method, the Mesh values of the spatial derivatives of the physical variables The use of the extra marching variables results from the need to construct the two-level explicit and nondissipative schemes which are at the core of the CE/SE development. It also results from the need to minimize the stencil while maintaining accuracy. In this paper using the 1D(sub (alpha)-mu) scheme as an example, the effect of this added complication on consistency, accuracy and operation count is assessed. As part of this effort, an equivalent yet more efficient form of the alpha-mu scheme in which the independent marching variables are the local fluxes tied to each mesh point is introduced. Also, the intriguing relations that exist among the alpha-mu. Leapfrog, and DuFort-Frankel schemes are further explored. In addition, the redundance of the Leapfrog, DUFort-Frankel, and Lax scheme and the remedy for this redundance are discussed. This paper is concluded with the construction and evaluation of a CE/SE solver for the inviscid Burger equation.

Chang, Sin-Chung↗
NASA NTRS1989

The design and application of upwind schemes on unstructured meshes

Solution and mesh generation algorithms for solving the Euler equations on unstructured meshes consisting of triangle and quadrilateral control volumes are presented. Cell-centered and mesh-vertex upwind finite-volume schemes are developed which utilize multi-dimensional monotone linear reconstruction procedures. These algorithms differ from existing algorithms (even on structured meshes). Numerical results in two dimensions are presented.

Barth, Timothy J.↗
NASA NTRS2009

Construction of Low Dissipative High Order Well-Balanced Filter Schemes for Non-Equilibrium Flows

The goal of this paper is to generalize the well-balanced approach for non-equilibrium flow studied by Wang et al. [26] to a class of low dissipative high order shock-capturing filter schemes and to explore more advantages of well-balanced schemes in reacting flows. The class of filter schemes developed by Yee et al. [30], Sjoegreen & Yee [24] and Yee & Sjoegreen [35] consist of two steps, a full time step of spatially high order non-dissipative base scheme and an adaptive nonlinear filter containing shock-capturing dissipation. A good property of the filter scheme is that the base scheme and the filter are stand alone modules in designing. Therefore, the idea of designing a well-balanced filter scheme is straightforward, i.e., choosing a well-balanced base scheme with a well-balanced filter (both with high order). A typical class of these schemes shown in this paper is the high order central difference schemes/predictor-corrector (PC) schemes with a high order well-balanced WENO filter. The new filter scheme with the well-balanced property will gather the features of both filter methods and well-balanced properties: it can preserve certain steady state solutions exactly; it is able to capture small perturbations, e.g., turbulence fluctuations; it adaptively controls numerical dissipation. Thus it shows high accuracy, efficiency and stability in shock/turbulence interactions. Numerical examples containing 1D and 2D smooth problems, 1D stationary contact discontinuity problem and 1D turbulence/shock interactions are included to verify the improved accuracy, in addition to the well-balanced behavior.

Wang, Wei↗
NASA NTRS1994

Compositional Variation in Apollo 16 Impact-Melt Breccias and Inferences for the Geology and Bombardment History of the Central Highlands of the Moon

High-precision data for the concentrations of a number of lithophile and siderophile elements were obtained on multiple subsamples from 109 impact-melt rocks and breccias (mostly crystalline) from the Apollo 16 site. Compositions of nearly all Apollo 16 melt rocks fall on one of two trends of increasing Sm concentration with increasing Sc concentration. The Eastern trend (lower Sm/Sc, Mg/Fe, and Sm/Yb ratios) consists of compositional groups 3 and 4 of previous classification schemes. These melt rocks are feldspathic, poor in incompatible and siderophile elements, and appear to have provenance in the Descartes formation to the east of the site. The Western trend (higher Sm/Sc. Mg/Fe, and Sm/ Yb ratios) consists of compositional groups 1 and 2. These relatively mafic, KREEP-bearing breccias are a major component (approx.35%) of the Cayley plains west of the site and are unusual, compared to otherwise similar melt breccias from other sites, in having high concentrations of Fe-Ni metal ( 1-2 %). The metal is the carrier of the low-Ir/Au (approx. 0.3 x chondritic) siderophile-element signature that is characteristic of the Apollo 16 site. Four compositionally distinct groups (1M, 1F, 2DB, and 2NR) of Western-trend melt breccias occur that are each represented by at least six samples. Compositional group 1 or previous classification schemes (the 'poikilitic' or 'LKFM' melt breccias) can be subdivided into two groups. Group 1M (represented by six samples, including 60315) is characterized by lower Al2O3 concentrations, higher MgO and alkali concentrations, and higher Mg/Fe and Cr/Sc ratios than group 1F (represented by fifteen samples, including 65015). Group 1M also has siderophile-element concentrations averaging about twice those of group lF and Ir/Au and Ir/Ni ratios that are even lower than those of other Western-trend melt rocks (Ir/Au = 0.24 +/- 0.03. CI-normalized). At the mafic extreme of group 2 ('VHA' melt breccias), the melt lithology occurring as clasts in feldspathic fragmental breccias from North Ray crater (group 2NR) is compositionally distinct from the melt lithology ofdimict breccias from the Cayley plains (group 2DB) in having higher concentrations of Sc, Cr, and heavy rare earth elements and lower concentrations of siderophile elements. The distinct siderophile-element signature (high absolute abundances, low Ir/Au ratio) suggest that the four groups ofmafic melt breccia are all somehow related. Ratios ofsome lithophile elements also suggest that they are more closely related to each other than then, are to melt breccias from other Apoll sites. However, none of the breccia compositions can be related to any of the others by any simple process of igneous fractionation or mixing involving common lunar materials. Thus, the origin of the four groups of mafic melt breccia is enigmatic. If they were produced in only one or two impacts, then a mechanism exists for generating regimes of impact-melt breccia in a single impact that are substantially different from each other in composition. For various reasons, including the problem of delivering large volumes of four different types of melt to the Apollo 16 site, it is unlikely that any of these breccias were produced in basin-forming impacts. If they were produced in as many as four crater-forming impacts, then the unusual siderophile-element signature is difficult to explain. Possible explanations are (1) the four groups of melt breccia all contain metal from a single, earlier impact, (2) they were each formed by related metal-rich meteoroids, or (3) some common postimpact process has resulted in metal of similar composition in each of four melt pools. Within a compositional group, most intrasample and intersample variation in lithophile element concentrations is caused by differences among samples in the proportion of a component of normative anorthosite or noritic anorthosite. In most cases, this compositional variation probably reflects variation in clast abundance. For group 2DB (and probably 2NR), differences in abundance of a component of ferroan anorthosite (estimated Al2O3 approx. 32%) accounts for the compositional variation. For groups 1M and 1F, the anorthositic component is more mafic (estimated Al203 approx. 26%). Some group-2 samples may be related by a troctolitic component of varying abundance.

Korotev, Randy L.↗
DOE OSTI2023

High precision tests of QCD without scale or scheme ambiguities: The 40th anniversary of the Brodsky–Lepage–Mackenzie method

A key issue in making precise predictions in QCD is the uncertainty in setting the renormalization scale μ r and thus determining the correct values of the QCD running coupling α s (μ r ) at each order in the perturbative expansion of a QCD observable. It has often been conventional to simply set the renormalization scale to the typical scale of the process Q and vary it in the range μ r $\in$ [Q/2, 2Q] in order to estimate the theoretical error. This is the practice of Conventional Scale Setting (CSS). The resulting CSS prediction will however depend on the theorist’s choice of renormalization scheme and the resulting pQCD series will diverge factorially. It will also disagree with renormalization scale setting used in QED and electroweak theory thus precluding grand unification. A solution to the renormalization scale-setting problem is offered by the Principle of Maximum Conformality (PMC), which provides a systematic way to eliminate the renormalization scale-and-scheme dependence in perturbative calculations. The PMC method has rigorous theoretical foundations, it satisfies Renormalization Group Invariance (RGI) and preserves all self-consistency conditions derived from the renormalization group. The PMC cancels the renormalon growth, reduces to the Gell-Mann–Low scheme in the N c → 0 Abelian limit and leads to scale- and scheme-invariant results. The PMC has now been successfully applied to many high-energy processes. In this article we summarize recent developments and results in solving the renormalization scale and scheme ambiguities in perturbative QCD. In particular, we present a recently developed method the PMC ∞ and its applications, comparing the results with CSS. The method preserves the property of renormalizable SU(N)/U(1) gauge theories defined as Intrinsic Conformality (iCF). This property underlies the scale invariance of physical observables and leads to a remarkably efficient method to solve the conventional renormalization scale ambiguity at every order in pQCD. This new method reflects the underlying conformal properties displayed by pQCD at NNLO, eliminates the scheme dependence of pQCD predictions and is consistent with the general properties of the PMC. A new method to identify conformal and β-terms, which can be applied either to numerical or to theoretical calculations is also shown. We present results for the thrust and C-parameter distributions in e + e - annihilation showing errors and comparison with the CSS. We also show results for a recent innovative comparison between the CSS and the PMC ∞ applied to the thrust distribution investigating both the QCD conformal window and the QED N c → 0 limit. In order to determine the thrust distribution along the entire renormalization group flow from the highest energies to zero energy, we consider the number of flavors near the upper boundary of the conformal window. In this flavor-number regime the theory develops a perturbative infrared interacting fixed point. These results show that PMC ∞ leads to higher precision and introduces new interesting features in the PMC. In fact, this method preserves with continuity the position of the peak, showing perfect agreement with the experimental data already at NNLO. We also show a detailed comparison of the PMC ∞ with the other PMC approaches: the multi-scale-setting approach (PMCm) and the single-scale-setting approach (PMCs) by comparing their predictions for three important fully integrated quantities R e+e- , R$_Τ$ and Γ (H→$b\bar{b}$) up to the four-loop accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1987

Remote sensing of earth terrain

A systematic approach for the identification of terrain media such as vegetation canopy, forest, and snow covered fields is developed using the optimum polarimetric classifier. The covariance matrices for the various terrain cover are computed from theoretical models of random medium by evaluating the full polarimetric scattering matrix elements. The optimal classification scheme makes use of a quadratic distance measure and is applied to classify a vegetation canopy consisting of both trees and grass. Experimentally measured data are used to validate the classification scheme. Theoretical probability of classification error using the full polarimetric matrix are compared with classification based on single features including the phase difference between the VV and HH polarization returns. It is shown that the full polarimetric results are optimal and provide better classification performance than single feature measurements.

Kong, J. A.↗
NASA NTRS1988

Remote sensing of earth terrain

A systematic approach for the identification of terrain media such as vegetation canopy, forest, and snow covered fields is developed using the optimum polarimetric classifier. The covariance matrices for the various terrain covers are computed from the theoretical models of random medium by evaluating the full polarimetric scattering matrix elements. The optimal classification scheme makes use of a quadratic distance measure and is applied to classify a vegetation canopy consisting of both trees and grass. Experimentally measured data are used to validate the classification scheme. Theoretical probability of classification error using the full polarimetric matrix are compared with classification based on single features including the phase difference between the VV and HH polarization returns. It is shown that the full polarimetric results are optimal and provide better classification performance than single feature measurements. A systematic approach is presented for obtaining the optimal polarimetric matched filter which produces maximum contrast between two scattering classes, each represented by its respective covariance matrix.

Kong, J. A.↗
NASA NTRS1988

Identification of terrain cover using the optimum polarimetric classifier

A systematic approach for the identification of terrain media such as vegetation canopy, forest, and snow-covered fields is developed using the optimum polarimetric classifier. The covariance matrices for various terrain cover are computed from theoretical models of random medium by evaluating the scattering matrix elements. The optimal classification scheme makes use of a quadratic distance measure and is applied to classify a vegetation canopy consisting of both trees and grass. Experimentally measured data are used to validate the classification scheme. Analytical and Monte Carlo simulated classification errors using the fully polarimetric feature vector are compared with classification based on single features which include the phase difference between the VV and HH polarization returns. It is shown that the full polarimetric results are optimal and provide better classification performance than single feature measurements.

Kong, J. A.↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2004

Nonlinear filtering and limiting in high order methods for ideal and non-ideal MHD

The various filtering mechanisms and base scheme options of the newly developed adaptive numerical dissipation control in spatially high order filter schemes for the ideal and non-ideal magnetohydrodynamics (MHD) equations are investigated. These filter schemes are applicable to complex unsteady MHD high-speed shock/shear/turbulence problems. They also provide a natural and efficient way for the minimization of Div(B) numerical error. The type of spatial base scheme to be used in conjunction with our filter idea is very general. For example, spectral, compact and non-compact spatially central finite difference schemes are possible candidates. The adaptive numerical dissipation mechanism consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and to leave the rest of the region free from numerical dissipation contamination. The numerical dissipation considered consists of high order linear dissipation for the suppression of high frequency oscillation and the nonlinear dissipative portion of high-resolution shock-capturing methods for discontinuity capturing. The applicable nonlinear dissipative portion of high-resolution shock-capturing methods is also very general. The objective of this paper is to investigate the performance of using compact and non-compact central base schemes in conjunction with three commonly used types of nonlinear numerical dissipation for both the ideal and non-ideal MHD. This extended abstract shows the performance of three nonlinear filters in conjunction with a sixth-order non-compact spatial central base scheme. In the final paper, the high order compact spatial central base scheme will be illustrated and compared with the non-compact base scheme. The reason for the investigation of the high order compact spatial central base scheme over the non-compact base scheme is to evaluate if additional accuracy can be gained in regions of fine scale turbulence that are away from shocks/shears.

Yee,H. C.↗
NASA NTRS1990

Viscous shock profiles and primitive formulations

Weak solutions of hyperbolic systems in primitive (non-conservation) form for which a consistent conservation form exists are considered. It is shown that primitive formulations, shock relations are not uniquely defined by the states to either side of the shock but also depend on the viscous path connecting the two. Scheme-dependent high order correction terms are derived that enforce consistent viscous shock profiles. The resulting primitive algorithm is conservative to the order of approximation. One dimensional Euler calculations of flows containing strong shocks clearly show that conservation errors in primitive flow calculations are of comparable quality.

Karni, S.↗