×

Inflow-implicit/outflow-explicit finite volume methods for solving advection equations. (English) Zbl 1310.65103

Summary: We introduce a new class of methods for solving non-stationary advection equations. The new methods are based on finite volume space discretizations and a semi-implicit discretization in time. Its basic idea is that outflow from a cell is treated explicitly while inflow is treated implicitly. This is natural, since we know what is outflowing from a cell at the old time step but we leave the method to resolve a system of equations determined by the inflows to a cell to obtain the solution values at the new time step. The matrix of the system in our inflow-implicit/outflow-explicit (IIOE) method is determined by the inflow fluxes which results in an M-matrix yielding favorable stability properties for the scheme. Since the explicit (outflow) part is not always dominated by the implicit (inflow) part and thus some oscillations can occur, we build a stabilization based on the upstream weighted averages with coefficients determined by the flux-corrected transport approach yielding high resolution versions, S\(^1\)IIOE and S\(^2\)IIOE, of the basic scheme. We prove that our new method is exact for any choice of a discrete time step on uniform rectangular grids in the case of constant velocity transport of quadratic functions in any dimension. We also show its formal second order accuracy in space and time for 1D advection problems with variable velocity. Although designed for non-divergence free velocity fields, we show that the basic IIOE scheme is locally mass conservative in case of divergence free velocity. Finally, we show L\(^2\)-stability for divergence free velocity in 1D on periodic domains independent of the choice of the time step, and L\(^\infty\)-stability for the stabilized high resolution variant of the scheme. Numerical comparisons with the purely explicit schemes like the fully explicit up-wind and the Lax-Wendroff schemes were discussed in [K. Mikula and M. Ohlberger, Springer Proc. Math. 4, 683–691 (2011; Zbl 1246.65167); “A new inflow-implicit/outflow-explicit finite volume method for solving variable velocity advection equations”, Preprint (2010)] where the basic IIOE was originally introduced. There it has been shown that the new scheme has good properties with respect to a balance of precision and CPU time related to a possible choice of larger time steps in our scheme. In this contribution we compare the new scheme and its stabilized variants with widely used fully implicit up-wind method. In this comparison our new schemes show better behavior with respect to stability and precision of computations for time steps several times exceeding the CFL stability condition. Our new stabilized methods are L\(^\infty\) stable, second order accurate for any smooth solution and with accuracy of order 2/3 for solutions with moving discontinuities. This is opposite to implicit up-wind schemes which have accuracy order 1/2 only. All these properties hold for any choice of time step thus making our new method attractive for practical applications.

MSC:

65M08 Finite volume methods for initial value and initial-boundary value problems involving PDEs

Citations:

Zbl 1246.65167

Software:

SHASTA
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Balažovjech, M.; Mikula, K., A higher order scheme for a tangentially stabilized plane curve shortening flow with a driving force, SIAM J. Sci. Comput., 33, 5, 2277-2294 (2011) · Zbl 1276.65048
[2] Boris, J. P.; Book, D. L., Flux-corrected transport: I. SHASTA, a fluid transport algorithm that works, J. Comput. Phys., 11, 38-69 (1973) · Zbl 0251.76004
[3] Bourgine, P.; Čunderlík, R.; Drblíková, O.; Mikula, K.; Peyriéras, N.; Remešíková, M.; Rizzi, B.; Sarti, A., 4D embryogenesis image analysis using PDE methods of image processing, Kybernetika, 46, 2 (2010) · Zbl 1198.94020
[4] Caselles, V.; Kimmel, R.; Sapiro, G., Geodesic active contours, Int. J. Comput. Vis., 22, 61-79 (1997) · Zbl 0894.68131
[5] Eymard, R.; Gallouet, T.; Herbin, R., The Finite Volume Methods, (Ciarlet, Ph.; Lions, J. L., Handbook for Numerical Analysis, vol. 7 (2000), Elsevier) · Zbl 0843.65068
[6] Eymard, R.; Handlovičová, A.; Mikula, K., Study of a finite volume scheme for the regularised mean curvature flow level set equation, IMA J. Numer. Anal., 31, 813-846 (2011) · Zbl 1241.65072
[7] Frolkovič, P.; Mikula, K., Flux-based level set method: a finite volume method for evolving interfaces, Appl. Numer. Math., 57, 4, 436-454 (2007) · Zbl 1120.65093
[8] Frolkovič, P.; Mikula, K., High-resolution flux-based level set method, SIAM J. Sci. Comput., 29, 2, 579-597 (2007) · Zbl 1141.76041
[9] Kichenassamy, S.; Kumar, A.; Olver, P.; Tannenbaum, A.; Yezzi, A., Conformal curvature flows: from phase transitions to active vision, Arch. Ration. Mech. Anal., 134, 275-301 (1996) · Zbl 0937.53029
[10] Kuzmin, D., Linearity-preserving flux correction and convergence acceleration for constrained Galerkin schemes, J. Comput. Appl. Math. (2011)
[11] LeVeque, R. J., Finite Volume Methods for Hyperbolic Problems, Cambridge Texts in Applied Mathematics (2002), Cambridge University Press · Zbl 1010.65040
[12] Mikula, K.; Ohlberger, M., A new level set method for motion in normal direction based on a semi-implicit forward-backward diffusion approach, SIAM J. Sci. Comput., 32, 3, 1527-1544 (2010) · Zbl 1223.65072
[13] Mikula, K.; Ohlberger, M., Inflow-implicit/outflow-explicit scheme for solving advection equations, (Fořt, J.; etal., Finite Volumes in Complex Applications VI, Problems & Perspectives. Finite Volumes in Complex Applications VI, Problems & Perspectives, Prague, June 6-10, 2011. Finite Volumes in Complex Applications VI, Problems & Perspectives. Finite Volumes in Complex Applications VI, Problems & Perspectives, Prague, June 6-10, 2011, Springer Proceedings in Mathematics, vol. 4 (2011), Springer-Verlag), 683-692 · Zbl 1246.65167
[15] Mikula, K.; Peyrieras, N.; Remešíková, M.; Sarti, A., 3D embryogenesis image segmentation by the generalized subjective surface method using the finite volume technique, (Eymard, R.; Herard, J. M., Finite Volumes for Complex Applications V: Problems and Perspectives (2008), ISTE and Wiley: ISTE and Wiley London), 585-592 · Zbl 1374.94259
[16] Mikula, K.; Ševčovič, D.; Balažovjech, M., A simple, fast and stabilized flowing finite volume method for solving general curve evolution equations, Commun. Comput. Phys., 7, 1, 195-211 (2010) · Zbl 1364.65175
[17] Sarti, A.; Malladi, R.; Sethian, J. A., Subjective surfaces, a method for completing missing boundaries, Proc. Natl. Acad. Sci. USA, 12, 6258-6263 (2000) · Zbl 0966.68214
[18] Sethian, J. A., Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Material Science (1999), Cambridge University Press: Cambridge University Press New York · Zbl 0973.76003
[19] Zalesak, S. T., Fully multidimensional flux-corrected transport algorithms for fluids, J. Comput. Phys., 31, 335-362 (1979) · Zbl 0416.76002
[20] Zanella, C.; Campana, M.; Rizzi, B.; Melani, C.; Sanguinetti, G.; Bourgine, P.; Mikula, K.; Peyrieras, N.; Sarti, A., Cells segmentation from 3-D confocal images of early zebrafish embryogenesis, IEEE Trans. Image Process., 19, 3, 770-781 (2010) · Zbl 1371.94449
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.