A limiting strategy for the back and forth error compensation and correction method for solving advection equations

We further study the properties of the back and forth error compensation and correction (BFECC) method for advection equations such as those related to the level set method and for solving Hamilton-Jacobi equations on unstructured meshes. In particular, we develop a new limiting strategy which requires another backward advection in time so that overshoots/undershoots on the new time level get exposed when they are transformed back to compare with the solution on the old time level. This new technique is very simple to implement even for unstructured meshes and is able to eliminate artifacts induced by jump discontinuities in derivatives of the solution as well as by jump discontinuities in the solution itself (even if the solution has large gradients in the vicinities of a jump). Typically, a formal second order method for solving a time dependent Hamilton-Jacobi equation requires quadratic interpolation in space. A BFECC method on the other hand only requires linear interpolation in each step, thus is local and easy to implement even for unstructured meshes.

[1]  Ignacio Llamas,et al.  FlowFixer: Using BFECC for Fluid Simulation , 2005, NPH.

[2]  Chi-Tien Lin,et al.  High-resolution Non-oscillatory Central Schemes for Hamilton-jacobi Equations , 2022 .

[3]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

[4]  R. Abgrall.  Numerical discretization of the first‐order Hamilton‐Jacobi equation on triangular meshes , 1996 .

[5]  Theodore Kim,et al.  A simple boiling module , 2007, SCA '07.

[6]  T. Dupont,et al.  Back and forth error compensation and correction methods for removing errors induced by uneven gradients of the level set function , 2003 .

[7]  Yingjie Liu,et al.  Back and forth error compensation and correction methods for semi-lagrangian schemes with application to level set interface computations , 2006, Math. Comput..

[8]  Wang Hai-bing,et al.  High-order essentially non-oscillatory schemes for Hamilton-Jacobi equations , 2006 .

[9]  S. Bryson,et al.  High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations , 2013 .

[10]  Ronald Fedkiw,et al.  An Unconditionally Stable MacCormack Method , 2008, J. Sci. Comput..

[11]  Ronald Fedkiw,et al.  An unconditionally stable fully conservative semi-Lagrangian method , 2010, J. Comput. Phys..

[12]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[13]  Eitan Tadmor,et al.  New High-Resolution Semi-discrete Central Schemes for Hamilton—Jacobi Equations , 2000 .

[14]  Irakli V. Gugushvili,et al.  Semi-Lagrangian Method for Advection Equation on GPU in Unstructured R3 Mesh for Fluid Dynamics Application , 2009 .

[15]  Alexander Kurganov,et al.  Semidiscrete Central-Upwind Schemes for Hyperbolic Conservation Laws and Hamilton-Jacobi Equations , 2001, SIAM J. Sci. Comput..

[16]  R. Courant,et al.  On the solution of nonlinear hyperbolic differential equations by finite differences , 1952 .

[17]  S. Zalesak,et al.  Fully Multidimensional Flux-Corrected Transport. , 1978 .

[18]  Ignacio Llamas,et al.  Advections with Significantly Reduced Dissipation and Diffusion , 2007, IEEE Transactions on Visualization and Computer Graphics.

[19]  J. Boris,et al.  Flux-Corrected Transport , 1997 .

[20]  Peter D. Lax,et al.  On the stability of difference approximations to solutions of hyperbolic equations with variable coefficients , 1961 .

[21]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[22]  Chi-Tien Lin,et al.  High-Resolution Nonoscillatory Central Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[23]  Steve Bryson,et al.  Semi-discrete central-upwind schemes with reduced dissipation for Hamilton–Jacobi equations , 2005 .

[24]  J. Strain.  Semi-Lagrangian Methods for Level Set Equations , 1999 .

[25]  Chi-Wang Shu,et al.  High-Order WENO Schemes for Hamilton-Jacobi Equations on Triangular Meshes , 2003, SIAM J. Sci. Comput..

[26]  Ignacio Llamas,et al.  Simulation of bubbles in foam with the volume control method , 2007, SIGGRAPH 2007.

[27]  S. Osher,et al.  Uniformly high order accuracy essentially non-oscillatory schemes III , 1987 .

[28]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[29]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[30]  Todd F. Dupont,et al.  Back and Forth Error Compensation and Correction Methods for Semi-Lagranging Schemes with Application to Interface Computation Using Level Set Method , 2004 .

[31]  Alexander Kurganov,et al.  Adaptive Central-Upwind Schemes for Hamilton–Jacobi Equations with Nonconvex Hamiltonians , 2006, J. Sci. Comput..

[32]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[33]  Frank Losasso,et al.  A fast and accurate semi-Lagrangian particle level set method , 2005 .

[34]  Qinjun Kang,et al.  Lattice Boltzmann method on quadtree grids. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[36]  S. Zalesak.  Introduction to “Flux-Corrected Transport. I. SHASTA, A Fluid Transport Algorithm That Works” , 1997 .