Numerical error control for second-order explicit TVD scheme with limiters in advection simulation

This paper analyzes the causes of the numerical errors in terms of numerical diffusion and compression arising from the use of explicit second-order total variation diminishing (TVD) schemes in one-dimensional advection simulation. It demonstrates that different TVD limiters may have very different performances in different advection simulations, because of the so-called numerical diffusion and compression. The accuracy of the computed results is found to depend on not only the limiter functions themselves but also the advection features such as the concentration distribution, advection velocity and time step, etc. According to such relations, the effective ranges of the MIN_MOD, Van Leer and SUPERBEE limiters are characterized by introducing a dimensionless parameter which reflects the key features of advections, aiming to provide an approach to select a proper TVD limiter in advection simulation.

[1]  Reinhard Hinkelmann.  Efficient Numerical Methods and Information-Processing Techniques for Modeling Hydro- and Environmental Systems , 2005 .

[2]  Lian-xia Li,et al.  An improved r-factor algorithm for TVD schemes , 2008 .

[3]  A. Harten.  High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[4]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[5]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

[6]  Xiaolin Wang,et al.  Solving the depth-integrated solute transport equation with a TVD-MacCormack scheme , 2010, Environ. Model. Softw..

[7]  T. E. Short,et al.  An exact peak capturing and Oscillation‐Free Scheme to solve advection‐dispersion transport equations , 1992 .

[8]  M. Darwish,et al.  TVD schemes for unstructured grids , 2003 .

[9]  Michael Sprenger,et al.  Transport of Stratospheric Air Masses to the Nepal Climate Observatory–Pyramid (Himalaya; 5079 m MSL): A Synoptic-Scale Investigation , 2012 .

[10]  Charles Hirsch.  The Mathematical Models for Fluid Flow Simulations at Various Levels of Approximation , 2007 .

[11]  V. G. Ferreira,et al.  A bounded upwinding scheme for computing convection-dominated transport problems , 2012 .

[12]  Nicolas G. Wright,et al.  2D Process-Based Morphodynamic Model for Flooding by Noncohesive Dyke Breach , 2014 .

[13]  Luis A. Lemus,et al.  VOF/FVM prediction and experimental validation for shear-thinning fluid column collapse , 2015, Comput. Math. Appl..

[14]  P. Sweby.  High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[15]  Kamal El Kadi Abderrezzak,et al.  A physical, movable-bed model for non-uniform sediment transport, fluvial erosion and bank failure in rivers , 2014 .

[16]  Weiming Wu,et al.  Depth-Averaged Two-Dimensional Model of Unsteady Flow and Sediment Transport due to Noncohesive Embankment Break/Breaching , 2012 .

[17]  Maurizio Brocchini,et al.  A multi-purpose, intra-wave, shallow water hydro-morphodynamic solver , 2012 .

[18]  Jingming Hou,et al.  A model for overland flow and associated processes within the Hydroinformatics Modelling System , 2014 .

[19]  Chongam Kim,et al.  Multi-dimensional limiting process for hyperbolic conservation laws on unstructured grids , 2005, J. Comput. Phys..

[20]  Stéphane Clain,et al.  Monoslope and multislope MUSCL methods for unstructured meshes , 2010, J. Comput. Phys..

[21]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[22]  George Huang,et al.  An Essentially Nonoscillatory High-Order Padé-Type (ENO-Padé) Scheme , 2002 .

[23]  Sarma L. Rani,et al.  Quantification of numerical diffusivity due to TVD schemes in the advection equation , 2014, J. Comput. Phys..

[24]  Peter K. Kitanidis,et al.  PARTICLE-TRACKING EQUATIONS FOR THE SOLUTION OF THE ADVECTION-DISPERSION EQUATION WITH VARIABLE COEFFICIENTS , 1994 .

[25]  Qiuhua Liang,et al.  Multislope MUSCL method applied to solve shallow water equations , 2014, Comput. Math. Appl..

[26]  John A. Trangenstein,et al.  Numerical Solution of Hyperbolic Partial Differential Equations , 2009 .

[27]  Christian Feigenwinter,et al.  A new mass conservation approach to the study of CO2 advection in an alpine forest , 2009 .

[28]  Jingming Hou,et al.  A new TVD method for advection simulation on 2D unstructured grids , 2013 .

[29]  Yang Liu,et al.  An oscillation-free high order TVD/CBC-based upwind scheme for convection discretization , 2011, Numerical Algorithms.

[30]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[31]  Jingming Hou,et al.  Improved total variation diminishing schemes for advection simulation on arbitrary grids , 2012 .

[32]  Pramod Kumar Sharma,et al.  Contaminant transport through fractured-porous media: An experimental study , 2014 .