A first-order computational algorithm for reaction-diffusion type equations via primal-dual hybrid gradient method

We propose an easy-to-implement iterative method for resolving the implicit (or semi-implicit) schemes arising in solving reaction-diffusion (RD) type equations. We formulate the nonlinear time implicit scheme as a min-max saddle point problem and then apply the primal-dual hybrid gradient (PDHG) method. Suitable precondition matrices are applied to the PDHG method to accelerate the convergence of algorithms under different circumstances. Furthermore, our method is applicable to various discrete numerical schemes with high flexibility. From various numerical examples tested in this paper, the proposed method converges properly and can efficiently produce numerical solutions with sufficient accuracy.

[1]  J. Carrillo,et al.  Structure preserving primal dual methods for gradient flows with nonlinear mobility transport distances , 2023, ArXiv.

[2]  Long Chen,et al.  Accelerated Gradient and Skew-Symmetric Splitting Methods for a Class of Monotone Operator Equations , 2023, 2303.09009.

[3]  S. Osher,et al.  High order spatial discretization for variational time implicit schemes: Wasserstein gradient flows and reaction-diffusion systems , 2023, ArXiv.

[4]  Long Chen,et al.  Transformed primal-dual methods for nonlinear saddle point systems , 2022, Journal of Numerical Mathematics.

[5]  Chi-Wang Shu,et al.  A primal-dual approach for solving conservation laws with implicit in time approximations , 2022, J. Comput. Phys..

[6]  M. H. Duong,et al.  Non-reversible processes: GENERIC, hypocoercivity and fluctuations , 2021, Nonlinearity.

[7]  Chun Liu,et al.  A second-order accurate, operator splitting scheme for reaction-diffusion systems in an energetic variational formulation , 2021, SIAM J. Sci. Comput..

[8]  Chun Liu,et al.  A structure-preserving, operator splitting scheme for reaction-diffusion equations with detailed balance , 2021, J. Comput. Phys..

[9]  Steven M. Wise,et al.  Convergence analysis of the variational operator splitting scheme for a reaction-diffusion system with detailed balance , 2021, SIAM J. Numer. Anal..

[10]  Xiaojing Ye,et al.  Numerical solution of inverse problems by weak adversarial networks , 2020, Inverse Problems.

[11]  Gang Bao,et al.  Weak Adversarial Networks for High-dimensional Partial Differential Equations , 2019, J. Comput. Phys..

[12]  P. Jimack,et al.  High Accuracy Benchmark Problems for Allen-Cahn and Cahn-Hilliard Dynamics , 2019, Communications in Computational Physics.

[13]  J. Carrillo,et al.  Primal Dual Methods for Wasserstein Gradient Flows , 2019, Foundations of Computational Mathematics.

[14]  Wuchen Li,et al.  Solving Large-Scale Optimization Problems with a Convergence Rate Independent of Grid Size , 2018, SIAM J. Numer. Anal..

[15]  Jiang Yang,et al.  The scalar auxiliary variable (SAV) approach for gradient flows , 2018, J. Comput. Phys..

[16]  Jie Shen,et al.  A New Class of Efficient and Robust Energy Stable Schemes for Gradient Flows , 2017, SIAM Rev..

[17]  Maxime Laborde,et al.  An unbalanced optimal transport splitting scheme for general advection-reaction-diffusion problems , 2017, ESAIM: Control, Optimisation and Calculus of Variations.

[18]  J. E. Guyer,et al.  Benchmark problems for numerical implementations of phase field models , 2016, 1610.00622.

[19]  Christian Clason,et al.  Primal-Dual Extragradient Methods for Nonlinear Nonsmooth PDE-Constrained Optimization , 2016, SIAM J. Optim..

[20]  Tuomo Valkonen,et al.  A primal–dual hybrid gradient method for nonlinear operators with applications to MRI , 2013, 1309.5032.

[21]  Carlos J. García-Cervera,et al.  A new approach for the numerical solution of diffusion equations with variable and degenerate mobility , 2013, J. Comput. Phys..

[22]  Keith Promislow,et al.  Variational Models of Network Formation and Ion Transport: Applications to Perfluorosulfonate Ionomer Membranes , 2012 .

[23]  Jianqin Zhou,et al.  On discrete cosine transform , 2011, ArXiv.

[24]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[25]  Yibao Li,et al.  An unconditionally stable hybrid numerical method for solving the Allen-Cahn equation , 2010, Comput. Math. Appl..

[26]  Xiaofeng Yang,et al.  Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , 2010 .

[27]  Andrea L Bertozzi,et al.  Dissipation and displacement of hotspots in reaction-diffusion models of crime , 2010, Proceedings of the National Academy of Sciences.

[28]  Jianfeng Zhu,et al.  Application of Discontinuous Galerkin Methods for Reaction-Diffusion Systems in Developmental Biology , 2009, J. Sci. Comput..

[29]  E. Favvas,et al.  What is spinodal decomposition , 2008 .

[30]  Qing Nie,et al.  Numerical methods for stiff reaction-diffusion systems , 2007 .

[31]  Arnold Neumaier,et al.  Introduction to Numerical Analysis , 2001 .

[32]  G. Russo,et al.  Implicit-explicit Runge-Kutta schemes for stiff systems of differential equations , 2000 .

[33]  Steven J. Ruuth.  Efficient Algorithms for Diffusion-Generated Motion by Mean Curvature , 1998 .

[34]  Hans Christian Öttinger,et al.  General projection operator formalism for the dynamics and thermodynamics of complex fluids , 1998 .

[35]  Miroslav Grmela,et al.  Dynamics and thermodynamics of complex fluids. I. Development of a general formalism , 1997 .

[36]  Miroslav Grmela,et al.  Dynamics and thermodynamics of complex fluids. II. Illustrations of a general formalism , 1997 .

[37]  S. Osher,et al.  Motion of multiple junctions: a level set approach , 1994 .

[38]  J. E. Pearson.  Complex Patterns in a Simple System , 1993, Science.

[39]  J. Keller,et al.  Fast reaction, slow diffusion, and curve shortening , 1989 .

[40]  G. Strang.  A proposal for toeplitz matrix calculations , 1986 .

[41]  J. Schnakenberg,et al.  Simple chemical reaction systems with limit cycle behaviour. , 1979, Journal of theoretical biology.

[42]  J. Cahn,et al.  A microscopic theory for antiphase boundary motion and its application to antiphase domain coasening , 1979 .

[43]  A. H. Sherman.  On Newton-Iterative Methods for the Solution of Systems of Nonlinear Equations , 1978 .

[44]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[45]  Keith Promislow,et al.  High accuracy solutions to energy gradient flows from material science models , 2014, J. Comput. Phys..

[46]  Mingqiang Zhu,et al.  An Efficient Primal-Dual Hybrid Gradient Algorithm For Total Variation Image Restoration , 2008 .

[47]  D. J. Eyre,et al.  An Unconditionally Stable One-Step Scheme for Gradient Systems , 1997 .

[48]  Kazufumi Ito,et al.  Optimum damping design for an abstract wave equation , 1996, Kybernetika.