ORIGINAL PAPERS

Exponential Time Differencing Method for a Reaction- Diffusion System with Free Boundary

Expand
  • 1. Department of Mathematics, University of California, San Diego, 9500 Gilman Drive, La Jolla, California, 92093-0112, USA;
    2. Department of Mathematics, University of South Carolina, Columbia, South Carolina, 29208, USA

Received date: 2022-07-31

  Revised date: 2023-02-09

  Online published: 2024-04-16

Supported by

This work is partially supported by the NSF DMS1853365.

Abstract

For reaction-diffusion equations in irregular domains with moving boundaries, the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes, while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary. It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties. Applying an implicit scheme may be able to remove the stability constraints on the time step, however, it usually requires solving a large global system of nonlinear equations for each time step, and the computational cost could be significant. Integration factor (IF) or exponential time differencing (ETD) methods are one of the popular methods for temporal partial differential equations (PDEs) among many other methods. In our paper, we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries. In particular, we rewrite all ETD schemes into a linear combination of specific ϕ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications, which offers significant computational advantages with adaptive Krylov subspaces. In addition, we extend this method by incorporating the level set method to solve the free boundary problem. The accuracy, stability, and efficiency of the developed method are demonstrated by numerical examples.

Cite this article

Shuang Liu, Xinfeng Liu . Exponential Time Differencing Method for a Reaction- Diffusion System with Free Boundary[J]. Communications on Applied Mathematics and Computation, 2024 , 6(1) : 354 -371 . DOI: 10.1007/s42967-023-00261-1

References

[1] Al-Mohy, A.H., Higham, N.J.:Computing the action of the matrix exponential with an application to exponential integrators. SIAM J. Sci. Comput. 33(2), 488-511 (2011)
[2] Apel, T., Sändig, A.M., Whiteman, J. R.:Graded mesh refinement and error estimates for finite element solutions of elliptic boundary value problems in non-smooth domains. Math. Methods Appl. Sci. 19(1), 63-85 (1996)
[3] Aslam, T.D.:A partial differential equation approach to multidimensional extrapolation. J. Comput. Phys. 193(1), 349-355 (2004)
[4] Barnett, G.A.:A robust RBF-FD formulation based on polyharmonic splines and polynomials. PhD Thesis. Citeseer (2015)
[5] Beylkin, G., Keiser, J.M., Vozovoi, L.:A new class of time discretization schemes for the solution of nonlinear PDEs. J. Comput. Phys. 147, 362-387 (1998)
[6] Braess, D.:The contraction number of a multigrid method for solving the Poisson equation. Numer. Math. 37(3), 387-404 (1981)
[7] Bunting, G., Du, Y., Krakowski, K.:Spreading speed revisited:analysis of a free boundary model. Net. Heterog. Med. 7(4), 583 (2012)
[8] Burrage, K., Butcher, J.C.:Stability criteria for implicit Runge-Kutta methods. SIAM J. Numer. Anal. 16(1), 46-57 (1979)
[9] Buzbee, B.L., Dorr, F.W., George, J.A., Golub, G.H.:The direct solution of the discrete Poisson equation on irregular regions. SIAM J. Numer. Anal. 8(4), 722-736 (1971)
[10] Caffarelli, L.A., Salsa, S.:A Geometric Approach to Free Boundary Problems, volume 68. American Mathematical Society, Providence, RI (2005)
[11] Caliari, M., Kandolf, P., Ostermann, A., Rainer, S.:The Leja method revisited:backward error analysis for the matrix exponential. SIAM J. Sci. Comput. 38(3), A1639-A1661 (2016)
[12] Cao, Y., Faghri, A., Chang, W.S.:A numerical analysis of Stefan problems for generalized multi-dimensional phase-change structures using the enthalpy transforming model. Intern. J. Heat. Mass. Trans. 32(7), 1289-1298 (1989)
[13] Chen, H., Min, C., Gibou, F.:A numerical scheme for the Stefan problem on adaptive Cartesian grids with supralinear convergence rate. J. Comput. Phys. 228(16), 5803-5818 (2009)
[14] Chen, S., Merriman, B., Osher, S., Smereka, P.:A simple level set method for solving Stefan problems. J. Comput. Phys. 135(1), 8-29 (1997)
[15] Cox, S.M., Matthews, P.C.:Exponential time differencing for stiff systems. J. Comput. Phys. 176(2), 430-455 (2002)
[16] Du, Q., Ju, L., Li, X., Qiao, Z.:Maximum bound principles for a class of semilinear parabolic equations and exponential time-differencing schemes. SIAM Rev. 63(2), 317-359 (2021)
[17] Du, Q., Zhu, W.:Stability analysis and applications of the exponential time differencing schemes. J. Comput. Math. 22, 200 (2004)
[18] Du, Q., Zhu, W.:Modified exponential time differencing schemes:analysis and applications. BIT. Numer. Math. 45, 307-328 (2005)
[19] Du, Y., Lin, Z.:Spreading-vanishing dichotomy in the diffusive logistic model with a free boundary. SIAM J. Math. Anal. 42(1), 377-405 (2010)
[20] Du, Y., Matano, H., Wang, K.:Regularity and asymptotic behavior of nonlinear Stefan problems. Archive Rat. Mech. Anal. 212(3), 957-1010 (2014)
[21] Franke, R.:Scattered data interpolation:tests of some methods. Math. Comput. 38(157), 181-200 (1982)
[22] Gaudreault, S., Pudykiewicz, J.A.:An efficient exponential time integration method for the numerical solution of the shallow water equations on the sphere. J. Comput. Phys. 322, 827-848 (2016)
[23] Gibou, F., Fedkiw, R.:A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem. J. Comput. Phys. 202(2), 577-601 (2005)
[24] Gibou, F., Fedkiw, R.P., Cheng, L.T., Kang, M.:A second-order-accurate symmetric discretization of the Poisson equation on irregular domains. J. Comput. Phys. 176(1), 205-227 (2002)
[25] Hairer, E., Wanner, G.:Stiff differential equations solved by Radau methods. J. Comput. Appl. Math. 111(1/2), 93-111 (1999)
[26] Hardy, R.L.:Multiquadric equations of topography and other irregular surfaces. J. Geophys. Res. 76(8), 1905-1915 (1971)
[27] Hochbruck, M., Lubich, C.:On Krylov subspace approximations to the matrix exponential operator. SIAM J. Numer. Anal. 34(5), 1911-1925 (1997)
[28] Hou, T.Y., Lowengrub, J.S., Shelley, M.J.:Removing the stiffness from interfacial flows with surface tension. J. Comput. Phys. 114, 312-338 (1994)
[29] Jiang, G.-S., Peng, D.:Weighted ENO schemes for Hamilton-Jacobi equations. SIAM J. Sci. Comput. 21(6), 2126-2143 (2000)
[30] Jiang, K., Ju, L., Li, J., Li, X.:Unconditionally stable exponential time differencing schemes for the mass-conserving Allen-Cahn equation with nonlocal and local effects. Numer. Methods Partial Differential Equations 38(6), 1636-1657 (2021)
[31] Johansen, H., Colella, P:A Cartesian grid embedded boundary method for Poisson's equation on irregular domains. J. Comput. Phys. 147(1), 60-85 (1998)
[32] Jomaa, Z., Macaskill, C.:The embedded finite difference method for the Poisson equation in a domain with an irregular boundary and Dirichlet boundary conditions. J. Comput. Phys. 202(2), 488-506 (2005)
[33] Jou, H.J., Leo, P.H., Lowengrub, J.S.:Microstructual evolution in inhomogeneous elastic media. J. Comput. Phys. 131, 109 (1997)
[34] Kassam, A.-K., Trefethen, L.N.:Fourth-order time stepping for stiff PDEs. SIAM J. Sci. Comput. 26, 1214-1233 (2005)
[35] Kreiss, H.-O., Petersson, N.A.:A second order accurate embedded boundary method for the wave equation with Dirichlet data. SIAM J. Sci. Comput. 27, 1141-1167 (2006)
[36] Kreiss, H.-O., Petersson, N.A., Yström, J.:Difference approximations for the second order wave equation. SIAM J. Numer. Anal. 40(5), 1940-1967 (2002)
[37] Kreiss, H.-O., Petersson, N.A., Yström, J.:Difference approximations of the Neumann problem for the second order wave equation. SIAM J. Numer. Anal. 42(3), 1292-1323 (2004)
[38] Lai, M.-C.:A note on finite difference discretizations for Poisson equation on a disk. Numer. Methods Partial Differential Equations 17(3), 199-203 (2001)
[39] Leo, P.H., Lowengrub, J.S., Nie, Q.:Microstructural evolution in orthotropic elastic media. J. Comput. Phys. 157, 44-88 (2000)
[40] Li, J., Ju, L., Cai, Y., Feng, X.:Unconditionally maximum bound principle preserving linear schemes for the conservative Allen-Cahn equation with nonlocal constraint. J. Sci. Comput. 87(3), 1-32 (2021)
[41] Liu, S.:Numerical methods for a class of reaction-diffusion equations with free boundaries. PhD Thesis. University of South Carolina (2019)
[42] Liu, X.D., Fedkiw, R.P., Kang, M.:A boundary condition capturing method for Poisson's equation on irregular domains. J. Comput. Phys. 160(1), 151-178 (2000)
[43] Louis, A.:Acceleration of convergence for finite element solutions of the Poisson equation. Numer. Math. 33(1), 43-53 (1979)
[44] Lowrie, R.B.:A comparison of implicit time integration methods for nonlinear relaxation and diffusion. J. Comput. Phys. 196(2), 566-590 (2004)
[45] Luan, V.T., Pudykiewicz, J.A., Reynolds, D.R.:Further development of efficient and accurate time integration schemes for meteorological models. J. Comput. Phys. 376, 817-837 (2019)
[46] Min, C., Gibou, F., Ceniceros, H.D.:A supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids. J. Comput. Phys. 218(1), 123-140 (2006)
[47] Ng, Y.T., Chen, H., Min, C., Gibou, F.:Guidelines for Poisson solvers on irregular domains with Dirichlet boundary conditions using the ghost fluid method. J. Sci. Comput. 41(2), 300-320 (2009)
[48] Niesen, J., Wright, W.M.:Algorithm 919:a Krylov subspace algorithm for evaluating the ϕ-functions appearing in exponential integrators. ACM Trans. Math. Soft. (TOMS) 38(3), 1-19 (2012)
[49] Oevermann, M., Klein, R.:A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces. J. Comput. Phys. 219(2), 749-769 (2006)
[50] Osher, S., Fedkiw, R.P.:Level set methods:an overview and some recent results. J. Comput. Phys. 169(2), 463-502 (2001)
[51] Osher, S., Fedkiw, R.P.:Level Set Methods and Dynamic Implicit Surfaces. Springer Verlag, New York (2002)
[52] Peng, D., Merriman, B., Osher, S., Zhao, H., Kang, M.:A PDE-based fast local level set method. J. Comput. Phys. 155(2), 410-438 (1999)
[53] Peng, Z., Appelö, D., Liu, S.:Universal AMG accelerated embedded boundary method without small cell stiffness. arXiv:2204.06083 (2022)
[54] Perrone, N., Kao, R.:A general finite difference method for arbitrary meshes. Comput. Struct. 5(1), 45-57 (1975)
[55] Piqueras, M.-A., Company, R., Jódar, L.:A front-fixing numerical method for a free boundary nonlinear diffusion logistic population model. J. Comput. Appl. Math. 309, 473-481 (2017)
[56] Rubinshteĭ, L.I.:The Stefan Problem. American Mathematical Society (1971)
[57] Sethian, J.A.:Fast marching methods. SIAM Rev. 41(2), 199-235 (1999)
[58] Sudicky, E.A.:The Laplace transform Galerkin technique:a time-continuous finite element theory and application to mass transport in groundwater. Wat. Res. Res. 25(8), 1833-1846 (1989)
[59] Sussman, M., Smereka, P., Osher, S.:A level set approach for computing solutions to incompressible two-phase flow. J. Comput. Phys. 114(1), 146-159 (1994)
[60] Tayler, A.B.:Free and moving boundary problems. By J. Crank. Clarendon, Oxford, 1984. 425 pp. £ 45.00. J. Fluid Mech. 158, 532-533 (1985)
[61] Yoon, G., Min, C.:Analyses on the finite difference method by Gibou, et al. for Poisson equation. J. Comput. Phys. 280, 184-194 (2015)
[62] Zapata, M.A.U., Hernández-López, F.J.:A GPU parallel finite volume method for a 3D Poisson equation on arbitrary geometries. Intern. J. Comb. Optimizat. Prob. Inform. 9(1), 3 (2018)
Options
Outlines

/