Partially Explicit Time Discretization for Nonlinear Time Fractional Diffusion Equations

Nonlinear time fractional partial differential equations are widely used in modeling and simulations. In many applications, there are high contrast changes in media properties. For solving these problems, one often uses coarse spatial grid for spatial resolution. For temporal discretization, implicit methods are often used. For implicit methods, though the time step can be relatively large, the equations are difficult to compute due to the nonlinearity and the fact that one deals with large-scale systems. On the other hand, the discrete system in explicit methods are easier to compute but it requires small time steps. In this work, we propose the partially explicit scheme following earlier works on developing partially explicit methods for nonlinear diffusion equations. In this scheme, the diffusion term is treated partially explicitly and the reaction term is treated fully explicitly. With the appropriate construction of spaces and stability analysis, we find that the required time step in our proposed scheme scales as the coarse mesh size, which creates a great saving in computing. The main novelty of this work is the extension of our earlier works for diffusion equations to time fractional diffusion equations. For the case of fractional diffusion equations, the constraints on time steps are more severe and the proposed methods alleviate this since the time step in partially explicit method scales as the coarse mesh size. We present stability results. Numerical results are presented where we compare our proposed partially explicit methods with a fully implicit approach. We show that our proposed approach provides similar results, while treating many degrees of freedom in nonlinear terms explicitly.

[1]  Yang Yang,et al.  Third-order conservative sign-preserving and steady-state-preserving time integrations and applications in stiff multispecies and multireaction detonations , 2019, J. Comput. Phys..

[2]  Yalchin Efendiev,et al.  Non-local multi-continua upscaling for flows in heterogeneous fractured media , 2017, J. Comput. Phys..

[3]  H. Owhadi,et al.  Metric‐based upscaling , 2007 .

[4]  Y. Efendiev,et al.  Generalized Multiscale Finite Element Methods. Nonlinear Elliptic Equations , 2013, 1304.5188.

[5]  R. Hilfer Applications Of Fractional Calculus In Physics , 2000 .

[6]  A. Abdulle Explicit Methods for Stiff Stochastic Differential Equations , 2012, CSE 2012.

[7]  Yalchin Efendiev,et al.  Contrast-independent partially explicit time discretizations for multiscale flow problems , 2021, J. Comput. Phys..

[8]  Willem Hundsdorfer,et al.  IMEX extensions of linear multistep methods with general monotonicity and boundedness properties , 2007, J. Comput. Phys..

[9]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[10]  Ioannis G. Kevrekidis,et al.  General Tooth Boundary Conditions for Equation Free Modeling , 2006, SIAM J. Sci. Comput..

[11]  Thomas Y. Hou,et al.  A Model Reduction Method for Multiscale Elliptic Pdes with Random Coefficients Using an Optimization Approach , 2018, Multiscale Model. Simul..

[12]  Zdzislaw Jackiewicz,et al.  Highly stable implicit-explicit Runge-Kutta methods , 2017 .

[13]  Yalchin Efendiev,et al.  Generalized Multiscale Finite Element Methods for Wave Propagation in Heterogeneous Media , 2013, Multiscale Model. Simul..

[14]  Eric T. Chung,et al.  Constraint Energy Minimizing Generalized Multiscale Finite Element Method for high-contrast linear elasticity problem , 2018 .

[15]  Qin Li,et al.  Exploring the Locally Low Dimensional Structure in Solving Random Elliptic PDEs , 2016, Multiscale Model. Simul..

[16]  Giovanni Samaey,et al.  Patch dynamics with buffers for homogenization problems , 2006, J. Comput. Phys..

[17]  Jens Eggers,et al.  The Explicit-Implicit-Null method: Removing the numerical instability of PDEs , 2014, J. Comput. Phys..

[18]  K. B. Oldham,et al.  The Fractional Calculus: Theory and Applications of Differentiation and Integration to Arbitrary Order , 1974 .

[19]  Willem Hundsdorfer,et al.  Stability of implicit-explicit linear multistep methods , 1997 .

[20]  Yalchin Efendiev,et al.  Numerical Homogenization of Monotone Elliptic Operators , 2003, Multiscale Model. Simul..

[21]  Yalchin Efendiev,et al.  Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods , 2016, J. Comput. Phys..

[22]  Yalchin Efendiev,et al.  Contrast-independent partially explicit time discretizations for nonlinear multiscale problems , 2021, Mathematics.

[23]  Thomas Y. Hou,et al.  An Adaptive Fast Solver for a General Class of Positive Definite Matrices Via Energy Decomposition , 2017, Multiscale Model. Simul..

[24]  Steven J. Ruuth Implicit-explicit methods for reaction-diffusion problems in pattern formation , 1995 .

[25]  E. Weinan,et al.  Effectiveness of implicit methods for stiff stochastic differential equations , 2008 .

[26]  Yalchin Efendiev,et al.  Partially explicit time discretization for time fractional diffusion equation , 2021, Fractional Calculus and Applied Analysis.

[27]  Yalchin Efendiev,et al.  Generalized multiscale finite element methods (GMsFEM) , 2013, J. Comput. Phys..

[28]  Björn Engquist,et al.  Heterogeneous multiscale methods for stiff ordinary differential equations , 2005, Math. Comput..

[29]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

[30]  Yalchin Efendiev,et al.  An Efficient Hierarchical Multiscale Finite Element Method for Stokes Equations in Slowly Varying Media , 2013, Multiscale Model. Simul..

[31]  Ying Li,et al.  Local discontinuous Galerkin methods with implicit-explicit multistep time-marching for solving the nonlinear Cahn-Hilliard equation , 2019, J. Comput. Phys..

[32]  H. Tchelepi,et al.  Multi-scale finite-volume method for elliptic problems in subsurface flow simulation , 2003 .

[33]  Yalchin Efendiev,et al.  Homogenization of nonlinear random parabolic operators , 2005, Advances in Differential Equations.

[34]  H. Srivastava,et al.  Theory and Applications of Fractional Differential Equations , 2006 .

[35]  P. Henning,et al.  A localized orthogonal decomposition method for semi-linear elliptic problems , 2012, 1211.3551.

[36]  Chuanju Xu,et al.  Finite difference/spectral approximations for the time-fractional diffusion equation , 2007, J. Comput. Phys..

[37]  Adrian Sandu,et al.  EPIRK-W and EPIRK-K Time Discretization Methods , 2017, J. Sci. Comput..

[38]  Yalchin Efendiev,et al.  Constraint energy minimizing generalized multiscale finite element method in the mixed formulation , 2017, Computational Geosciences.

[39]  J. M. Aldaz Strengthened Cauchy-Schwarz and Hölder inequalities , 2013, 1302.2254.

[40]  I. Podlubny Fractional differential equations : an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications , 1999 .

[41]  Yalchin Efendiev,et al.  Mixed Generalized Multiscale Finite Element Methods and Applications , 2014, Multiscale Model. Simul..

[42]  Björn Engquist,et al.  A multiscale method for highly oscillatory ordinary differential equations with resonance , 2008, Math. Comput..

[43]  Frédéric Legoll,et al.  An MsFEM Type Approach for Perforated Domains , 2013, Multiscale Model. Simul..

[44]  Yalchin Efendiev,et al.  Fast online generalized multiscale finite element method using constraint energy minimization , 2017, J. Comput. Phys..