Abstract
In this paper we generalise to non-uniform grids of quad-tree type the Compact WENO reconstruction of Levy et al. (SIAM J Sci Comput 22(2):656–672, 2000), thus obtaining a truly two-dimensional non-oscillatory third order reconstruction with a very compact stencil and that does not involve mesh-dependent coefficients. This latter characteristic is quite valuable for its use in h-adaptive numerical schemes, since in such schemes the coefficients that depend on the disposition and sizes of the neighbouring cells (and that are present in many existing WENO-like reconstructions) would need to be recomputed after every mesh adaption. In the second part of the paper we propose a third order h-adaptive scheme with the above-mentioned reconstruction, an explicit third order TVD Runge–Kutta scheme and the entropy production error indicator proposed by Puppo and Semplice (Commun Comput Phys 10(5):1132–1160, 2011). After devising some heuristics on the choice of the parameters controlling the mesh adaption, we demonstrate with many numerical tests that the scheme can compute numerical solution whose error decays as \(\langle N\rangle ^{-3}\), where \(\langle N\rangle \) is the average number of cells used during the computation, even in the presence of shock waves, by making a very effective use of h-adaptivity and the proposed third order reconstruction.
Similar content being viewed by others
Notes
As an exception, during the first time step, if the initial condition is known analytically, it is more accurate to use the analytic expression to set the cell averages in the newly created cells.
For simplicity we assume \(\nu \) is independent on the nature of the singularity, which of course is not true in general.
References
Abgrall, R.: On essentially non-oscillatory schemes on unstructured meshes: analysis and implementation. J. Comput. Phys. 114, 45–58 (1994)
Aràndiga, F., Baeza, A., Belda, A.M., Mulet, P.: Analysis of WENO schemes for full and global accuracy. SIAM J. Numer. Anal. 49(2), 893–915 (2011)
Arvanitis, C., Delis, A.I.: Behavior of finite volume schemes for hyperbolic conservation laws on adaptive redistributed spatial grids. SIAM J. Sci. Comput. 28(5), 1927–1956 (2006). doi:10.1137/050632853
Bastian, P., Blatt, M., Dedner, A., Engwer, C., Fahlke, J., Gräser, C., Klöfkorn, R., Nolte, M., Ohlberger, M., Sander, O.: DUNE web page (2011). http://www.dune-project.org
Bastian, P., Blatt, M., Dedner, A., Engwer, C., Klöfkorn, R., Kornhuber, R., Ohlberger, M., Sander, O.: A generic grid interface for parallel and adaptive scientific computing. Part II: Implementation and tests in DUNE. Computing 82(2–3), 121–138 (2008)
Bastian, P., Blatt, M., Dedner, A., Engwer, C., Klöfkorn, R., Ohlberger, M., Sander, O.: A generic grid interface for parallel and adaptive scientific computing. Part I: Abstract framework. Computing 82(2–3), 103–119 (2008)
Berger, M.J., LeVeque, R.J.: Adaptive mesh refinement using wave-propagation algorithms for hyperbolic systems. SIAM J. Numer. Anal. 35, 2298–2316 (1998)
Burri, A., Dedner, A., Klöfkorn, R., Ohlberger, M.: An efficient implementation of an adaptive and parallel grid in dune. In: Notes on Numerical Fluid Mechanics, vol. 91, pp. 67–82 (2006). See also http://aam.mathematik.uni-freiburg.de/IAM/Research/alugrid
Čada, M., Torrilhon, M.: Compact third-order limiter functions for finite volume methods. J. Comput. Phys. 228(11), 4118–4145 (2009). doi:10.1016/j.jcp.2009.02.020
Constantinescu, E.M., Sandu, A.: Multirate timestepping methods for hyperbolic conservation laws. J. Sci. Comput. 33(3), 239–278 (2007). doi:10.1007/s10915-007-9151-y
Cravero, I., Semplice, M.: On the accuracy of WENO and CWENO reconstructions of third order on nonuniform meshes. Preprint arXiv:1503.00736
Davies-Jones, R.: Comments on “Kinematic analysis of frontogenesis associated with a nondivergent vortex”. J. Atmos. Sci. 42(19), 2073–2075 (1985)
Dumbser, M., Käser, M.: Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems. J. Comput. Phys. 221(2), 693–723 (2007)
Feng, H., Hu, F., Wang, R.: A new mapped weighted essentially non-oscillatory scheme. J. Sci. Comput. 51, 449–473 (2012)
Fürst, J.: A weighted least square scheme for compressible flows. Flow Turbul. Combust. 76(4), 331–342 (2006)
Gottlieb, S., Ketcheson, D., Shu, C.W.: Strong Stability Preserving Runge–Kutta and Multistep Time Discretizations. World Scientific Publishing, Hackensack (2011)
Henrick, A.K., Aslam, T.D., Powers, J.M.: Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points. J. Comput. Phys. 207, 542–567 (2005)
Hu, C., Shu, C.W.: Weighted essentially non-oscillatory schemes on triangular meshes. J. Comput. Phys. 150(1), 97–127 (1999)
Karni, S., Kurganov, A.: Local error analysis for approximate solutions of hyperbolic conservation laws. Adv. Comput. Math. 22(1), 79–99 (2005). doi:10.1007/s10444-005-7099-8
Ketcheson, D.I., Parsani, M., LeVeque, R.J.: High-order wave propagation algorithms for hyperbolic systems. SIAM J. Sci. Comput. 35(1), A351–A377 (2013)
Kirby, R.: On the convergence of high resolution methods with multiple time scales for hyperbolic conservation laws. Math. Comput. 72(243), 1239–1250 (2003). doi:10.1090/S0025-5718-02-01469-2
Kolb, O.: On the full and global accuracy of a compact third order WENO scheme. SIAM J. Numer. Anal. 52(5), 2335–2355 (2014)
Levy, D., Puppo, G., Russo, G.: Compact central WENO schemes for multidimensional conservation laws. SIAM J. Sci. Comput. 22(2), 656–672 (2000)
Li, W., Ren, Y.X.: High-order k-exact weno finite volume schemes for solving gas dynamic euler equations on unstructured grids. Int. J. Numer. Methods Fluids 70(6), 742–763 (2012)
Lörcher, F., Gassner, G., Munz, C.D.: A discontinuous Galerkin scheme based on a space-time expansion. I. Inviscid compressible flow in one space dimension. J. Sci. Comput. 32(2), 175–199 (2007). doi:10.1007/s10915-007-9128-x
Mandli, K.T., Ketcheson, D.I., et al.: Pyclaw software (2011). http://numerics.kaust.edu.sa/pyclaw
Ohlberger, M.: A review of a posteriori error control and adaptivity for approximations of non-linear conservation laws. Int. J. Numer. Methods Fluids 59, 333–354 (2009). doi:10.1002/fld.1686
Osher, S., Sanders, R.: Numerical approximations to nonlinear conservation laws with locally varying time and space grids. Math. Comput. 41(164), 321–336 (1983). doi:10.2307/2007679
Puppo, G.: Numerical entropy production for central schemes. SIAM J. Sci. Comput. 25(4), 1382–1415 (2003/04)
Puppo, G., Semplice, M.: Numerical entropy and adaptivity for finite volume schemes. Commun. Comput. Phys. 10(5), 1132–1160 (2011)
Puppo, G., Semplice, M.: Well-balanced high order 1d schemes on non-uniform grids and entropy residuals. Preprint arXiv:1403.4112
Rogerson, A., Meiburg, E.: A numerical study of the convergence properties of ENO schemes. J. Sci. Comput. 5(2), 151–167 (1990)
Semplice, M., Coco, A.: dune-fv software. http://www.personalweb.unito.it/matteo.semplice/codes.htm
Shi, J., Hu, C., Shu, C.W.: A technique of treating negative weights in WENO schemes. J. Comput. Phys. 175(1), 108–127 (2002)
Shu, C.: Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. In: Advanced Numerical Approximation of Nonlinear Hyperbolic Equations (Cetraro, 1997), Lecture Notes in Mathematics, vol. 1697, pp. 325–432. Springer, Berlin (1998)
Shu, C.W., Osher, S.: Efficient implementation of essentially nonoscillatory shock-capturing schemes. J. Comput. Phys. 77(2), 439–471 (1988)
Acknowledgments
This work was supported by “National Group for Scientific Computation (GNCS-INDAM)”.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Semplice, M., Coco, A. & Russo, G. Adaptive Mesh Refinement for Hyperbolic Systems Based on Third-Order Compact WENO Reconstruction. J Sci Comput 66, 692–724 (2016). https://doi.org/10.1007/s10915-015-0038-z
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-015-0038-z