New Renormalization Schemes for Conductivity Upscaling in Heterogeneous Media

Two new renormalization schemes for conductivity upscaling in heterogeneous media are presented. The schemes follow previous ones by performing the renormalization over square cells of size 2 ^d with d being the dimensionality. Contrasting with previous schemes, the two-dimensional scheme makes use of the exact 2 × 2 block-conductivity. On the basis of the structure of the exact two-dimensional block-conductivity, an analogous three-dimensional scheme is proposed. The new schemes are tested on a number of benchmark problems and are shown to be significantly more accurate than existing schemes.     

Fifth‐order Hermitian schemes for computational linear aeroacoustics

We develop a class of fifth‐order methods to solve linear acoustics and/or aeroacoustics. Based on local Hermite polynomials, we investigate three competing strategies for solving hyperbolic linear problems with a fifth‐order accuracy. A one‐dimensional (1D) analysis in the Fourier series makes it possible to classify these possibilities. Then, numerical computations based on the 1D scalar advection equation support two possibilities in order to update the discrete variable and its first and second derivatives: the first one uses a procedure similar to that of Cauchy–Kovaleskaya (the ‘Δ‐P5 scheme’); the second one relies on a semi‐discrete form and evolves in time the discrete unknowns by using a five‐stage Runge–Kutta method (the ‘RGK‐P5 scheme’). Although the RGK‐P5 scheme shares the same local spatial interpolator with the Δ‐P5 scheme, it is algebraically simpler. However, it is shown numerically that its loss of compactness reduces its domain of stability. Both schemes are then extended to bi‐dimensional acoustics and aeroacoustics. Following the methodology validated in (J. Comput. Phys. 2005; 210 :133–170; J. Comput. Phys. 2006; 217 :530–562), we build an algorithm in three stages in order to optimize the procedure of discretization. In the ‘reconstruction stage’, we define a fifth‐order local spatial interpolator based on an upwind stencil. In the ‘decomposition stage’, we decompose the time derivatives into simple wave contributions. In the ‘evolution stage’, we use these fluctuations to update either by a Cauchy–Kovaleskaya procedure or by a five‐stage Runge–Kutta algorithm, the discrete variable and its derivatives. In this way, depending on the configuration of the ‘evolution stage’, two fifth‐order upwind Hermitian schemes are constructed. The effectiveness and the exactitude of both schemes are checked by their applications to several 2D problems in acoustics and aeroacoustics. In this aim, we compare the computational cost and the computation memory requirement for each solution. The RGK‐P5 appears as the best compromise between simplicity and accuracy, while the Δ‐P5 scheme is more accurate and less CPU time consuming, despite a greater algebraic complexity. Copyright © 2007 John Wiley & Sons, Ltd.     

Differential-equation-based representation of truncation errors for accurate numerical simulation

High-order compact finite difference schemes for two-dimensional convection-diffusion-type differential equations with constant and variable convection coefficients are derived. The governing equations are employed to represent leading truncation terms, including cross-derivatives, making the overall O(h⁴) schemes conform to a 3 × 3 stencil. We show that the two-dimensional constant coefficient scheme collapses to the optimal scheme for the one-dimensional case wherein the finite difference equation yields nodally exact results. The two-dimensional schemes are tested against standard model problems, including a Navier-Stokes application. Results show that the two schemes are generally more accurate, on comparable grids, than O(h²) centred differencing and commonly used O(h) and O(h³) upwinding schemes.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

An efficient high order direct ALE ADER finite volume scheme with a posteriori limiting for hydrodynamics and magnetohydrodynamics

In this paper, we present a new family of direct arbitrary–Lagrangian–Eulerian (ALE) finite volume schemes for the solution of hyperbolic balance laws on unstructured meshes in multiple space dimensions. The scheme is designed to be high‐order accurate both in space and time, and the mesh motion, which provides the new mesh configuration at the next time step, is taken into account in the final finite volume scheme that is based directly on a space‐time conservation formulation of the governing PDE system. To improve the computational efficiency of the algorithm, high order of accuracy in space is achieved using the a posteriori MOOD limiting strategy that allows the reconstruction procedure to be carried out with only one reconstruction stencil for any order of accuracy. We rely on an element‐local space‐time Galerkin finite element predictor on moving curved meshes to obtain a high‐order accurate one‐step time discretization, while the mesh velocity is computed by means of a suitable nodal solver algorithm that might also be supplemented with a local rezoning procedure to improve the mesh quality. Next, the old mesh configuration at time level tⁿ is connected to the new one at t^n + 1 by straight edges, hence providing unstructured space‐time control volumes, on the boundary of which the numerical flux has to be integrated. Here, we adopt a quadrature‐free integration, in which the space‐time boundaries of the control volumes are split into simplex sub‐elements that yield constant space‐time normal vectors and Jacobian matrices. In this way, the integrals over the simplex sub‐elements can be evaluated once and for all analytically during a preprocessing step. We apply the new high‐order direct ALE algorithm to the Euler equations of compressible gas dynamics (also referred to as hydrodynamics equations) as well as to the magnetohydrodynamics equations and we solve a set of classical test problems in two and three space dimensions. Numerical convergence rates are provided up to fifth order of accuracy in 2D and 3D for both hyperbolic systems considered in this paper. Finally, the efficiency of the new method is measured and carefully compared against the original formulation of the algorithm that makes use of a WENO reconstruction technique and Gaussian quadrature formulae for the flux integration: depending on the test problem, the new class of very efficient direct ALE schemes proposed in this paper can run up to ≈12 times faster in the 3D case. Copyright © 2016 John Wiley & Sons, Ltd.     

An enthalpy‐preserving shock‐capturing term for residual distribution schemes

In this contribution, we investigate strategies to perform shock‐capturing computation of steady hypersonic flow fields by means of residual distribution schemes. The ultimate objective is the computation of flow solutions for which the correct upstream enthalpy value is recovered in the postshock region. To this end, the parallelism existing between the classical Bx scheme and the stabilized finite element techniques is exploited. The simple Lax‐Friedrichs dissipation term is leveraged to build two new residual distribution schemes. Upon testing on both inviscid and viscous steady problems, solutions obtained with one of the two schemes are shown to recover the correct upstream total enthalpy level in the postshock region. This last scheme provides also improved wall pressure and skin friction predictions; heat transfer predictions are, unfortunately, similar to those offered by the Bx scheme. A conjecture for explaining this behavior is exposed.     

Improvement of multistep WENO scheme and its extension to higher orders of accuracy

This paper presents an efficient procedure for overcoming the deficiency of weighted essentially non‐oscillatory schemes near discontinuities. Through a thorough incorporation of smoothness indicators into the weights definition, up to ninth‐order accurate multistep methods are devised, providing weighted essentially non‐oscillatory schemes with enhanced order of convergence at transition points from smooth regions to a discontinuity, while maintaining stability and the essentially non‐oscillatory behavior. We also provide a detailed analysis of the resolution power and show that the solution enhancements of the new method at smooth regions come from their ability to render smoothness indicators closer to uniformity. The new scheme exhibits similar fidelity as other multistep schemes; however, with superior characteristics in terms of robustness and efficiency, as no logical statements or mapping function is needed. Extensions to higher orders of accuracy present no extra complexity. Numerical solutions of linear advection problems and nonlinear hyperbolic conservation laws are used to demonstrate the scheme's improved behavior for shock‐capturing problems. Copyright © 2016 John Wiley & Sons, Ltd.     

求解无粘可压Euler方程组的虚拟流方法

首先将三阶Godunov型半离散中心迎风格式推广到四阶,之后再将该新的四阶半离散中心迎风格式与Level Set方法以及虚拟流方法结合起来,成功地处理了非反应激波问题和多介质流中的爆轰间断问题。由于Level Set函数能隐式地追踪到界面的位置,而虚拟流的构造能隐式地捕捉到界面的边界条件,故而本文的方法可以很自然地推广到多维情况。     

A new well‐balanced Hermite weighted essentially non‐oscillatory scheme for shallow water equations

Hermite weighted essentially non‐oscillatory (HWENO) methods were introduced in the literature, in the context of Euler equations for gas dynamics, to obtain high‐order accuracy schemes characterized by high compactness (e.g. Qiu and Shu, J. Comput. Phys. 2003; 193 :115). For example, classical fifth‐order weighted essentially non‐oscillatory (WENO) reconstructions are based on a five‐cell stencil whereas the corresponding HWENO reconstructions are based on a narrower three‐cell stencil. The compactness of the schemes allows easier treatment of the boundary conditions and of the internal interfaces. To obtain this compactness in HWENO schemes both the conservative variables and their first derivatives are evolved in time, whereas in the original WENO schemes only the conservative variables are evolved. In this work, an HWENO method is applied for the first time to the shallow water equations (SWEs), including the source term due to the bottom slope, to obtain a fourth‐order accurate well‐balanced compact scheme. Time integration is performed by a strong stability preserving the Runge–Kutta method, which is a five‐step and fourth‐order accurate method. Besides the classical SWE, the non‐homogeneous equations describing the time and space evolution of the conservative variable derivatives are considered here. An original, well‐balanced treatment of the source term involved in such equations is developed and tested. Several standard one‐dimensional test cases are used to verify the high‐order accuracy, the C‐property and the good resolution properties of the model. Copyright © 2010 John Wiley & Sons, Ltd.     

Performance of very‐high‐order upwind schemes for DNS of compressible wall‐turbulence

The purpose of the present paper is to evaluate very‐high‐order upwind schemes for the direct numerical simulation (DNS ) of compressible wall‐turbulence. We study upwind‐biased (UW ) and weighted essentially nonoscillatory (WENO ) schemes of increasingly higher order‐of‐accuracy (J. Comp. Phys. 2000; 160 :405–452), extended up to WENO 17 (AIAA Paper 2009‐1612, 2009). Analysis of the advection–diffusion equation, both as Δx→0 (consistency), and for fixed finite cell‐Reynolds‐number Re_Δx (grid‐resolution), indicates that the very‐high‐order upwind schemes have satisfactory resolution in terms of points‐per‐wavelength (PPW ). Computational results for compressible channel flow (Re∈[180, 230]; M?_CL ∈[0.35, 1.5]) are examined to assess the influence of the spatial order of accuracy and the computational grid‐resolution on predicted turbulence statistics, by comparison with existing compressible and incompressible DNS databases. Despite the use of baseline O(Δt²) time‐integration and O(Δx²) discretization of the viscous terms, comparative studies of various orders‐of‐accuracy for the convective terms demonstrate that very‐high‐order upwind schemes can reproduce all the DNS details obtained by pseudospectral schemes, on computational grids of only slightly higher density. Copyright © 2009 John Wiley & Sons, Ltd.     

Chebyshev collocation method and multidomain decomposition for the incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations in primitive variables have been solved by a pseudospectral Chebyshev method using a semi-implicit fractional step scheme. The latter has been adapted to the particular features of spectral collocation methods to develop the monodomain algorithm. In particular, pressure and velocity collocated on the same nodes are sought in a polynomial space of the same order; the cascade of scalar elliptic problems arising after the spatial collocation is solved using finite difference preconditioning. With the present procedure spurious pressure modes do not pollute the pressure field. As a natural development of the present work a multidomain extent was devised and tested. The original domain is divided into a union of patching sub-rectangles. Each scalar problem obtained after spatial collocation is solved by iterating by subdomains. For steady problems a C¹ solution is recovered at the interfaces upon convergence, ensuring a spectrally accurate solution. A number of test cases have been solved to validate the algorithm in both its single-block and multidomain configurations. The preliminary results achieved indicate that collocation methods in multidomain configurations might become a viable alternative to the spectral element technique for accurate flow prediction.     

A volume-agglomeration multirate time advancing for high Reynolds number flow simulation

A frequent configuration in computational fluid mechanics combines an explicit time advancing scheme for accuracy purposes and a computational grid with a very small portion of much smaller elements than in the remaining mesh. Two examples of such situations are the travel of a discontinuity followed by a moving mesh, and the large eddy simulation of high Reynolds number flows around bluff bodies where together very thin boundary layers and vortices of much more important size need to be captured. For such configurations, multistage explicit time advancing schemes with global time stepping are very accurate but very CPU consuming. In order to reduce this problem, the multirate (MR) time stepping approach represents an interesting improvement. The objective of such schemes, which allow to use different time steps in the computational domain, is to avoid penalizing the computational cost of the time advancement of unsteady solutions that would become large due to the use of small global time steps imposed by the smallest elements such as those constituting the boundary layers. In the present work, a new MR scheme based on control volume agglomeration is proposed for the solution of the compressible Navier-Stokes equations equipped with turbulence models. The method relies on a prediction step where large time steps are performed with an evaluation of the fluxes on macrocells for the smaller elements for stability purpose and a correction step in which small time steps are employed. The accuracy and efficiency of the proposed method are evaluated on several benchmarks flows: the problem of a moving contact discontinuity (inviscid flow), the computation with a hybrid turbulence model of flows around bluff bodies like a flow around a space probe model at Reynolds number 10⁶, a circular cylinder at Reynolds number 8.4 × 10⁶, and two tandem cylinders at Reynolds number 1.66 × 10⁵ and 1.4 × 10⁵.     

MOSQUITO: An efficient finite difference scheme for numerical simulation of 2D advection

An explicit finite difference method for the treatment of the advective terms in the 2D equation of unsteady scalar transport is presented. The scheme is a conditionally stable extension to two dimensions of the popular QUICKEST scheme. It is deduced imposing the vanishing of selected components of the truncation error for the case of steady uniform flow. The method is then extended to solve the conservative form of the depth‐averaged transport equation. Details of the accuracy and stability analysis of the numerical scheme with test case results are given, together with a comparison with other existing schemes suitable for the long‐term computations needed in environmental modelling. Although with a truncation error of formal order 0(ΔxΔt, ΔyΔt, Δt²), the present scheme is shown actually to be of an accuracy comparable with schemes of third‐order in space, while requiring a smaller computational effort and/or having better stability properties. In principle, the method can be easily extended to the 3D case. Copyright © 1999 John Wiley & Sons, Ltd.     

Evaluation of various high-order-accuracy schemes with and without flux limiters

Conventional high-order schemes with reduced levels of numerical diffusion produce results with spurious oscillations in areas where steep velocity gradients exist. To prevent the development of non-physical oscillations in the solution, several monotonic schemes have been proposed. In this work, three monotonic schemes, namely Van Leer's scheme, Roe's flux limiter and the third-order SHARP scheme, are compared and evaluated against schemes without flux limiters. The latter schemes include the standard first-order upwind scheme, the second-order upwind scheme and the QUICK scheme. All the above schemes are applied to four two-dimensional problems: (i) rotation of a scalar ‘cone’ field, (ii) transport of a scalar ‘square’ field, (iii) mixing of a cold with a hot front and (iv) deformation of a scalar ‘cone’ field. These problems test the ability of the selected schemes to produce oscillation-free and accurate results in critical convective situations. The evaluation of the schemes is based on several aspects, such as accuracy, economy and complexity. The tests performed in this work reveal the merits and demerits of each scheme. It is concluded that high-order schemes with flux limiters can significantly improve the accuracy of the results.     

Implementation of some higher-order convection schemes on non-uniform grids

A generalized formulation is applied to implement the quadratic upstream interpolation (QUICK) scheme, the second-order upwind (SOU) scheme and the second-order hybrid scheme (SHYBRID) on non-uniform grids. The implementation method is simple. The accuracy and efficiency of these higher-order schemes on non-uniform grids are assessed. Three well-known bench mark convection-diffusion problems and a fluid flow problem are revisited using non-uniform grids. These are: (1) transport of a scalar tracer by a uniform velocity field; (2) heat transport in a recirculating flow; (3) two-dimensional non-linear Burgers equations; and (4) a two-dimensional incompressible Navier-Stokes flow which is similar to the classical lid-driven cavity flow. The known exact solutions of the last three problems make it possible to thoroughly evaluate accuracies of various uniform and non-uniform grids. Higher accuracy is obtained for fewer grid points on non-uniform grids. The order of accuracy of the examined schemes is maintained for some tested problems if the distribution of non-uniform grid points is properly chosen.     

Assessment of the finite volume method applied to the v2 − f model

The objective of this paper is to assess the accuracy of low‐order finite volume (FV) methods applied to the v² ? f turbulence model of Durbin (Theoret. Comput. Fluid Dyn. 1991; 3 :1–13) in the near vicinity of solid walls. We are not (like many others) concerned with the stability of solvers ‐ the topic at hand is simply whether the mathematical properties of the v² ? f model can be captured by the given, widespread, numerical method. The v² ? f model is integrated all the way up to solid walls, where steep gradients in turbulence parameters are observed. The full resolution of wall gradients imposes quite high demands on the numerical schemes and it is not evident that common (second order) FV codes can fully cope with such demands. The v² ? f model is studied in a statistically one‐dimensional, fully developed channel flow where we compare FV schemes with a highly accurate spectral element reference implementation. For the FV method a higher‐order face interpolation scheme, using Lagrange interpolation polynomials up to arbitrary order, is described. It is concluded that a regular second‐order FV scheme cannot give an accurate representation of all model parameters, independent of mesh density. To match the spectral element solution an extended source treatment (we use three‐point Gauss–Lobatto quadrature), as well as a higher‐order discretization of diffusion is required. Furthermore, it is found that the location of the first internal node need to be well within y⁺=1. Copyright © 2009 John Wiley & Sons, Ltd.     

Positive‐definite q‐families of continuous subcell Darcy‐flux CVD(MPFA) finite‐volume schemes and the mixed finite element method

A new family of locally conservative cell‐centred flux‐continuous schemes is presented for solving the porous media general‐tensor pressure equation. A general geometry‐permeability tensor approximation is introduced that is piecewise constant over the subcells of the control volumes and ensures that the local discrete general tensor is elliptic. A family of control‐volume distributed subcell flux‐continuous schemes are defined in terms of the quadrature parametrization q (Multigrid Methods. Birkhauser: Basel, 1993; Proceedings of the 4th European Conference on the Mathematics of Oil Recovery, Norway, June 1994; Comput. Geosci. 1998; 2 :259–290), where the local position of flux continuity defines the quadrature point and each particular scheme. The subcell tensor approximation ensures that a symmetric positive‐definite (SPD) discretization matrix is obtained for the base member (q=1) of the formulation. The physical‐space schemes are shown to be non‐symmetric for general quadrilateral cells. Conditions for discrete ellipticity of the non‐symmetric schemes are derived with respect to the local symmetric part of the tensor. The relationship with the mixed finite element method is given for both the physical‐space and subcell‐space q‐families of schemes. M‐matrix monotonicity conditions for these schemes are summarized. A numerical convergence study of the schemes shows that while the physical‐space schemes are the most accurate, the subcell tensor approximation reduces solution errors when compared with earlier cell‐wise constant tensor schemes and that subcell tensor approximation using the control‐volume face geometry yields the best SPD scheme results. A particular quadrature point is found to improve numerical convergence of the subcell schemes for the cases tested. Copyright © 2007 John Wiley & Sons, Ltd.     

Study of a staggered fourth‐order compact scheme for unsteady incompressible viscous flows

The objective of this paper is the development and assessment of a fourth‐order compact scheme for unsteady incompressible viscous flows. A brief review of the main developments of compact and high‐order schemes for incompressible flows is given. A numerical method is then presented for the simulation of unsteady incompressible flows based on fourth‐order compact discretization with physical boundary conditions implemented directly into the scheme. The equations are discretized on a staggered Cartesian non‐uniform grid and preserve a form of kinetic energy in the inviscid limit when a skew‐symmetric form of the convective terms is used. The accuracy and efficiency of the method are demonstrated in several inviscid and viscous flow problems. Results obtained with different combinations of second‐ and fourth‐order spatial discretizations and together with either the skew‐symmetric or divergence form of the convective term are compared. The performance of these schemes is further demonstrated by two challenging flow problems, linear instability in plane channel flow and a two‐dimensional dipole–wall interaction. Results show that the compact scheme is efficient and that the divergence and skew‐symmetric forms of the convective terms produce very similar results. In some but not all cases, a gain in accuracy and computational time is obtained with a high‐order discretization of only the convective and diffusive terms. Finally, the benefits of compact schemes with respect to second‐order schemes is discussed in the case of the fully developed turbulent channel flow. Copyright © 2008 John Wiley & Sons, Ltd.     

A critical evaluation of seven discretization schemes for convection–diffusion equations

A comparative study of seven discretization schemes for the equations describing convection-diffusion transport phenomena is presented. The (differencing) schemes considered are the conventional central- and upwind-difference schemes, together with the Leonard,¹ Leonard upwind¹ and Leonard super upwind difference¹ schemes. Also tested are the so called locally exact difference scheme² and the quadratic-upstream difference scheme.^3,4 In multidimensional problems errors arise from ‘false-diffusion’ and function approximations. It is asserted that false diffusion is essentially a multidimensional source of error. No mesh constraints are associated with errors in function approximation and discretization. Hence errors associated with discretization only may be investigated via one-dimensional problems. Thus, although the above schemes have been tested for one- and two-dimensional flows with sources, only the former are presented here. For 1D flows, the Leonard super upwind difference scheme and the locally exact scheme are shown to be far superior in accuracy to the others at all Peclet numbers and for most source distributions, for the test cases considered. Furthermore, the latter is shown to be considerably cheaper in computational terms than the former. The stability of the schemes and their CPU time requirements are also discussed.     

Some explicit triangular finite element schemes for the Euler equations

Several explicit schemes are presented for triangular P₀ and P₁ finite elements. A first-order accurate upwind P₀ scheme is compared to a FLIC type method. A second-order accurate Richtmyer scheme is constructed. Applications are given for the Euler system of conservation laws in the 2-dimensional case.