Some applications of the concept of minimized integrated exponential error for low dispersion and low dissipation

Several techniques to optimize parameters that regulate dispersion and dissipation effects in finite difference schemes have been devised in our previous works. They all use the notion that dissipation neutralizes dispersion. These techniques are the minimized integrated square difference error (MISDE) and the minimized integrated exponential error for low dispersion and low dissipation (MIEELDLD). It is shown in this work based on several numerical schemes tested that the technique of MIEELDLD is more accurate than MISDE to optimize the parameters that regulate dispersion and dissipation effects with the aim of improving the shock‐capturing properties of numerical methods. First, we consider the family of third‐order schemes proposed by Takacs. We use the techniques MISDE and MIEELDLD to optimize two parameters, namely, the cfl number and another variable which also controls dispersion and dissipation. Second, these two techniques are used to optimize a numerical scheme proposed by Gadd. Moreover, we compute the optimal cfl for some multi‐level schemes in 1D. Numerical tests for some of these numerical schemes mentioned above are performed at different cfl numbers and it is shown that the results obtained are dependent on the cfl number chosen. The errors from the numerical results have been quantified into dispersion and dissipation using a technique devised by Takacs. Finally, we make use of a composite scheme made of corrected Lax–Friedrichs and the two‐step Lax–Friedrichs schemes like the CFLF4 scheme at its optimal cfl number, to solve some problems in 2D, namely: solid body rotation test, acoustics and the circular Riemann problem. Copyright © 2011 John Wiley & Sons, Ltd.     

Optimised composite numerical schemes in 2‐D for hyperbolic conservation laws

One of the techniques available for optimising parameters that regulate dispersion and dissipation effects in finite difference schemes is the concept of minimised integrated exponential error for low dispersion and low dissipation. In this paper, we work essentially with the two‐dimensional (2D) Corrected Lax–Friedrichs and Lax–Friedrichs schemes applied to the 2D scalar advection equation. We examine the shock‐capturing properties of these two numerical schemes, and observe that these methods are quite effective from the point of being able to control computational noise and having a large range of stability. To improve the shock‐capturing efficiency of these two methods, we derive composite methods using the idea of predictor/corrector or a linear combination of the two schemes. The optimal cfl number for some of these composite schemes are computed. Some numerical experiments are carried out in two dimensions such as cylindrical explosion, shock‐focusing, dam‐break and Riemann gas dynamics tests. The modified equations of some of the composite schemes when applied to the 2D scalar advection equation are obtained. We also perform some convergence tests to obtain the order of accuracy and show that better results in terms of shock‐capturing property are obtained when the optimal cfl obtained using minimised integrated exponential error for low dispersion and low dissipation is used. Copyright © 2011 John Wiley & Sons, Ltd.     

Entropy dissipation scheme and minimums‐increase‐and‐maximums‐decrease slope limiter

In this paper, we continue to study the entropy dissipation scheme developed in former. We start with a numerical study of the scheme without the entropy dissipation term on the linear advection equation, which shows that the scheme is stable and numerical dissipation and numerical dispersion free for smooth solutions. However, the numerical results for discontinuous solutions show nonlinear instabilities near jump discontinuities. This is because the scheme enforces two related conservation properties in the computation. With this study, we design a so‐called ‘minimums‐increase‐and‐maximums‐decrease’ slope limiter in the reconstruction step of the scheme and delete the entropy dissipation in the linear fields and reduce the entropy dissipation terms in the nonlinear fields. Numerical experiments show improvements of the designed scheme compared with the results presented in former. However, the minimums‐increase‐and‐maximums‐decrease limiter is still not perfect yet, and better slope limiters are still sought. Copyright © 2011 John Wiley & Sons, Ltd.     

A new weighting method for improving the WENO‐Z scheme

The calculation of the weight of each substencil is very important for a weighted essentially nonoscillatory (WENO) scheme to obtain high‐order accuracy in smooth regions and keep the essentially nonoscillatory property near discontinuities. The weighting function introduced in the WENO‐Z scheme provides a straightforward method to analyze the accuracy order in smooth regions. In this paper, we construct a new sixth‐order global smoothness indicator (GSI‐6) and a function about GSI‐6 and the local smoothness indicators (IS_k) to calculate the weights. The analysis and numerical results show that, with the new weights, the scheme satisfies the sufficient condition for the fifth‐order convergence in smooth regions even at critical points. Meanwhile, it can also maintain low dissipation for discontinuous solutions due to relative large weights assigned to discontinuous substencils.     

L2Roe: a low dissipation version of Roe's approximate Riemann solver for low Mach numbers

A modification of the Roe scheme called L²Roe for low dissipation low Mach Roe is presented. It reduces the dissipation of kinetic energy at the highest resolved wave numbers in a low Mach number test case of decaying isotropic turbulence. This is achieved by scaling the jumps in all discrete velocity components within the numerical flux function. An asymptotic analysis is used to show the correct pressure scaling at low Mach numbers and to identify the reduced numerical dissipation in that regime. Furthermore, the analysis allows a comparison with two other schemes that employ different scaling of discrete velocity jumps, namely, LMRoe and a method of Thornber et al. To this end, we present for the first time an asymptotic analysis of the last method. Numerical tests on cases ranging from low Mach number (M_∞=0.001) to hypersonic (M_∞=5) viscous flows are used to illustrate the differences between the methods and to show the correct behavior of L²Roe. No conflict is observed between the reduced numerical dissipation and the accuracy or stability of the scheme in any of the investigated test cases. Copyright © 2015 John Wiley & Sons, Ltd.     

Optimized sixth‐order monotonicity‐preserving scheme by nonlinear spectral analysis

In this paper, sixth‐order monotonicity‐preserving optimized scheme (OMP6) for the numerical solution of conservation laws is developed on the basis of the dispersion and dissipation optimization and monotonicity‐preserving technique. The nonlinear spectral analysis method is developed and is used for the purpose of minimizing the dispersion errors and controlling the dissipation errors. The new scheme (OMP6) is simple in expression and is easy for use in CFD codes. The suitability and accuracy of this new scheme have been tested through a set of one‐dimensional, two‐dimensional, and three‐dimensional tests, including the one‐dimensional Shu–Osher problem, the two‐dimensional double Mach reflection, and the Rayleigh–Taylor instability problem, and the three‐dimensional direct numerical simulation of decaying compressible isotropic turbulence. All numerical tests show that the new scheme has robust shock capturing capability and high resolution for the small‐scale waves due to fewer numerical dispersion and dissipation errors. Moreover, the new scheme has higher computational efficiency than the well‐used WENO schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

A physically based flux limiter for QUICK calculations of advective scalar transport

Transient, advective transport of a contaminant into a clean domain will exhibit a moving sharp front that separates contaminated and clean regions. Due to ‘numerical diffusion’—the combined effects of ‘cross‐wind diffusion’ and ‘artificial dispersion’—a numerical solution based on a first‐order (upwind) treatment will smear out the sharp front. The use of higher‐order schemes, e.g. QUICK (quadratic upwinding) reduces the smearing but can introduce non‐physical oscillations in the solution. A common approach to reduce numerical diffusion without oscillations is to use a scheme that blends low‐order and high‐order approximations of the advective transport. Typically, the blending is based on a parameter that measures the local monotonicity in the predicted scalar field. In this paper, an alternative approach is proposed for use in scalar transport problems where physical bounds C_Low?C?C_High on the scalar are known a priori. For this class of problems, the proposed scheme switches from a QUICK approximation to an upwind approximation whenever the predicted upwind nodal value falls outside of the physical range [C_Low, C_High]. On two‐dimensional steady‐state and one‐dimensional transient test problems predictions obtained with the proposed scheme are essentially indistinguishable from those obtained with monotonic flux‐limiter schemes. An analysis of the modified equation explains the observed performance of first‐ and second‐order time‐stepping schemes in predicting the advective transport of a step. In application to the transient two‐dimensional problem of contaminate transport into a streambed, predictions obtained with the proposed flux‐limiter scheme agree with those obtained with a scheme from the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical analysis and construction of limiter of high resolution difference scheme

I.IntroductionAdjustingandgoverningproperlythenumericaldissipationanddispersionarethekeytotheconstructionofhighresolutionandnon-oscillationscheme.SinceVanLcerll]introducednuxlimiterandobtainedhighresolutionandnon-oscillationscheme,choisillgandconstructingsuitablelimiterbecomesanimportantwayofdesigninghighresolutiollalld11on-oscillationscheme,andadjustsnumericaldissipationeffects,dispersioneffectsaswellasgroupvelocityeffects.Roe121'l'l,Chakravarthy'andOSherl4]andetal.providedvariouslimitersan…     

A robust low‐dissipation AUSM‐family scheme for numerical shock stability on unstructured grids

The simple low‐dissipation advection upwind splitting method (SLAU) scheme is a parameter‐free, low‐dissipation upwind scheme that has been applied in a wide range of aerodynamic numerical simulations. In spite of its successful applications, the SLAU scheme could be showing shock instabilities on unstructured grids, as many other contact resolved upwind schemes. Therefore, a hybrid upwind flux scheme is devised for improving the shock stability of SLAU scheme, without compromising on accuracy and low Mach number performance. Numerical flux function of the hybrid scheme is written in a general form, in which only the scalar dissipation term is different from that of the SLAU scheme. The hybrid dissipation term is defined by using a differentiable multidimensional‐shock‐detection pressure weight function, and the dissipation term of SLAU scheme is combined with that of the Van Leer scheme. Furthermore, the hybrid dissipation term is only applied for the solution of momentum fluxes in numerical flux function. Based on the numerical test results, the hybrid scheme is deemed to be a successful improvement on the shock stability of SLAU scheme, without compromising on the efficiency and accuracy. Copyright © 2016 John Wiley & Sons, Ltd.     

A finite volume‐based high‐order,Cartesian cut‐cell method for wave propagation

Computational aeroacoustics requires numerical techniques capable of yielding low artificial dispersion and dissipation to preserve the amplitude and the frequency characteristics of the physical processes. Furthermore, for engineering applications, the techniques need to handle irregular geometries associated with realistic configurations. We address these issues by developing an optimized prefactored compact finite volume (OPC‐fv) scheme along with a Cartesian cut‐cell technique. The OPC‐fv scheme seeks to minimize numerical dispersion and dissipation while satisfying the conservation laws. The cut‐cell approach treats irregularly shaped boundaries using divide‐and‐merge procedures for the Cartesian cells while maintaining a desirable level of accuracy. We assess these techniques using several canonical test problems, involving different levels of physical and geometric complexities. Richardson extrapolation is an effective tool for evaluating solutions of no high gradients or discontinuities, and is used to evaluate the performance of the solution technique. It is demonstrated that while the cut‐cell method has a modest effect on the order of accuracy, it is a robust method. The combined OPC‐fv scheme and the Cartesian cut‐cell technique offer good accuracy as well as geometric flexibility. Copyright © 2008 John Wiley & Sons, Ltd.     

Improved quadratic isogeometric element simulation of one-dimensional elastic wave propagation with central difference method

Two improved isogeometric quadratic elements and the central difference scheme are used to formulate the solution procedures of transient wave propagation problems. In the proposed procedures, the lumped matrices corresponding to the isogeometric elements are obtained. The stability conditions of the solution procedures are also acquired. The dispersion analysis is conducted to obtain the optimal Courant-Friedrichs-Lewy (CFL) number or time-step sizes corresponding to the spatial isogeometric elements. The dispersion analysis shows that the isogeometric quadratic element of the fourth-order dispersion error (called the isogeometric analysis (IGA)-f quadratic element) provides far more desirable numerical dissipation/dispersion than the element of the second-order dispersion error (called the IGA-s quadratic element) when appropriate time-step sizes are selected. The numerical simulations of one-dimensional (1D) transient wave propagation problems demonstrate the effectiveness of the proposed solution procedures.     

关于无振荡、无自由参数有限元格式的研究

利用双曲守恒律方程的Ｔａｙｌｏｒ弱解表达式，建立了有限元法修正方程，选择合适的展开式系数能得到一系列数值格式．通过稳定性分析研究了格式的稳定性、色散误差与有限元修正方程导数项系数之间的关系，该关系与差分法的ＮＮＤ格式一致．在选定格式下，通过ＣＦＬ数可控制有限元离散解的振荡而使格式不含自由参数．最后，用数值算例验证了这一关系，并在二、三维欧拉方程作了推广应用．     

A new matrix dissipation model for central scheme

In this paper, a simple and efficient improvement to the famous Swanson–Turkel matrix dissipation model for the central scheme is proposed. In the new matrix dissipation model, the accuracy is improved by eliminating the second‐difference dissipation added to the characteristic fields representing the vorticity waves. This strategy is proposed based on analyzing the flow‐physics about shock‐vortex interaction using the Rankine–Hugoniot jump condition. In this paper, the behavior of central scheme for rotational flow is also theoretically and numerically analyzed. Results show a newfound problem of the original scalar and matrix dissipation models, in which for rotational flow excessive second‐difference dissipation is added due to the pressure‐based shock sensor. With current new matrix dissipation model improved accuracy is obtained at minimal cost overhead, especially, in the highly vortical region where the second‐difference dissipation is reduced. At the same time, it preserves the excellent shock capturing capability and convergence speed of original method. Numerical properties of this new matrix dissipation model are validated with a series of numerical experiments and results comparison with original model verifies improved performance of current method. Copyright © 2013 John Wiley & Sons, Ltd.     

Implementing a high‐order accurate implicit operator scheme for solving steady incompressible viscous flows using artificial compressibility method

This paper uses a fourth‐order compact finite‐difference scheme for solving steady incompressible flows. The high‐order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two‐dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier–Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth‐order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block‐tridiagonal matrix inversions. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. In this study, the high‐order compact implicit operator scheme is also extended for computing three‐dimensional incompressible flows. The accuracy and efficiency of this high‐order compact method are demonstrated for different incompressible flow problems. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated. Test cases considered herein for validating the results are incompressible flows in a 2‐D backward facing step, a 2‐D cavity and a 3‐D cavity at different flow conditions. Results obtained for these cases are in good agreement with the available numerical and experimental results. The study shows that the scheme is robust, efficient and accurate for solving incompressible flow problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Depth‐integrated free‐surface flow with a two‐layer non‐hydrostatic formulation

This paper presents the derivation of a depth‐integrated wave propagation and runup model from a system of governing equations for two‐layer non‐hydrostatic flows. The governing equations are transformed into an equivalent, depth‐integrated system, which separately describes the flux‐dominated and dispersion‐dominated processes. The depth‐integrated system reproduces the linear dispersion relation within a 5 error for water depth parameter up to kd = 11, while allowing direct implementation of a momentum conservation scheme to model wave breaking and a moving‐waterline technique for runup calculation. A staggered finite‐difference scheme discretizes the governing equations in the horizontal dimension and the Keller box scheme reconstructs the non‐hydrostatic terms in the vertical direction. An semi‐implicit scheme integrates the depth‐integrated flow in time with the non‐hydrostatic pressure determined from a Poisson‐type equation. The model is verified with solitary wave propagation in a channel of uniform depth and validated with previous laboratory experiments for wave transformation over a submerged bar, a plane beach, and fringing reefs. Copyright © 2011 John Wiley & Sons, Ltd.     

Combined compact difference method for solving the incompressible Navier–Stokes equations

This paper presents a numerical method for solving the two‐dimensional unsteady incompressible Navier–Stokes equations in a vorticity–velocity formulation. The method is applicable for simulating the nonlinear wave interaction in a two‐dimensional boundary layer flow. It is based on combined compact difference schemes of up to 12th order for discretization of the spatial derivatives on equidistant grids and a fourth‐order five‐ to six‐alternating‐stage Runge–Kutta method for temporal integration. The spatial and temporal schemes are optimized together for the first derivative in a downstream direction to achieve a better spectral resolution. In this method, the dispersion and dissipation errors have been minimized to simulate physical waves accurately. At the same time, the schemes can efficiently suppress numerical grid‐mesh oscillations. The results of test calculations on coarse grids are in good agreement with the linear stability theory and comparable with other works. The accuracy and the efficiency of the current code indicate its potential to be extended to three‐dimensional cases in which full boundary layer transition happens. Copyright © 2011 John Wiley & Sons, Ltd.     

A horizontally curvilinear non‐hydrostatic model for simulating nonlinear wave motion in curved boundaries

A horizontally curvilinear non‐hydrostatic free surface model that embeds the second‐order projection method, the so‐called θ scheme, in fractional time stepping is developed to simulate nonlinear wave motion in curved boundaries. The model solves the unsteady, Navier–Stokes equations in a three‐dimensional curvilinear domain by incorporating the kinematic free surface boundary condition with a top‐layer boundary condition, which has been developed to improve the numerical accuracy and efficiency of the non‐hydrostatic model in the standard staggered grid layout. The second‐order Adams–Bashforth scheme with the third‐order spatial upwind method is implemented in discretizing advection terms. Numerical accuracy in terms of nonlinear phase speed and amplitude is verified against the nonlinear Stokes wave theory over varying wave steepness in a two‐dimensional numerical wave tank. The model is then applied to investigate the nonlinear wave characteristics in the presence of dispersion caused by reflection and diffraction in a semicircular channel. The model results agree quantitatively with superimposed analytical solutions. Finally, the model is applied to simulate nonlinear wave run‐ups caused by wave‐body interaction around a bottom‐mounted cylinder. The numerical results exhibit good agreement with experimental data and the second‐order diffraction theory. Overall, it is shown that the developed model, with only three vertical layers, is capable of accurately simulating nonlinear waves interacting within curved boundaries. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of the generalized MacCormack scheme and its extension to low Mach number flows

The low Mach number performance of the MacCormack scheme is examined. The inherent dissipation in the scheme is found to suffer from the degradation in accuracy observed with traditional, density‐based methods for compressible flows. Two specific modifications are proposed, leading to the formation of the generalized MacCormack scheme within a dual‐time framework (called GMC‐PC). The first modification involves reformulating the flux by splitting it into particle convection and acoustic parts, with the former terms treated using the traditional MacCormack discretization and the latter terms augmented by the addition of a pressure‐based artificial dissipation. The second modification involves a reformulation of the traditional nonlinear fix introduced by MacCormack in 1971, which is found to be necessary to suppress pressure oscillations at low Mach numbers. The new scheme is demonstrated to have superior performance, independent of Mach number, compared with standard MacCormack implementations using several canonical test problems. Copyright © 2017 John Wiley & Sons, Ltd.     

A low-dissipation third-order weighted essentially nonoscillatory scheme with a new reference smoothness indicator

The classical third-order weighted essentially nonoscillatory (WENO) scheme is notoriously dissipative as it loses the optimal order of accuracy at critical points and its two-point finite difference in the smoothness indicators is unable to differentiate the critical point from the discontinuity. In recent years, modifications to the smoothness indicators and weights of the classical third-order WENO scheme have been reported to reduce numerical dissipation. This article presents a new reference smoothness indicator for constructing a low-dissipation third-order WENO scheme. The new reference smoothness indicator is a nonlinear combination of the local and global stencil smoothness indicators. The resulting WENO-Rp3 scheme with the power parameter p=1.5 achieves third-order accuracy in smooth regions including critical points and has low dissipation, but numerical results show this scheme cannot keep the ENO property near discontinuities. The recommended WENO-R3 scheme (p=1) keeps the ENO property and performs better than several recently developed third-order WENO schemes.     

Baroclinic stability for a family of two‐level,semi‐implicit numerical methods for the 3D shallow water equations

The baroclinic stability of a family of two time‐level, semi‐implicit schemes for the 3D hydrostatic, Boussinesq Navier–Stokes equations (i.e. the shallow water equations), which originate from the TRIM model of Casulli and Cheng (Int. J. Numer. Methods Fluids 1992; 15 :629–648), is examined in a simple 2D horizontal–vertical domain. It is demonstrated that existing mass‐conservative low‐dissipation semi‐implicit methods, which are unconditionally stable in the inviscid limit for barotropic flows, are unstable in the same limit for baroclinic flows. Such methods can be made baroclinically stable when the integrated continuity equation is discretized with a barotropically dissipative backwards Euler scheme. A general family of two‐step predictor‐corrector schemes is proposed that have better theoretical characteristics than existing single‐step schemes. Copyright © 2006 John Wiley & Sons, Ltd.