Operator splitting and discontinuous Galerkin methods for advection-reaction-diffusion problem. Application to plant root growth


  • Emilie Peynaud CIRAD, UMR AMAP, YaoundГ©, Cameroun AMAP, Univ Montpellier, CIRAD, CNRS, INRA, IRD, Montpellier, France University of YaoundГ© I, Cameroon African Cetic, University of YaoundГ© I, Cameroon http://orcid.org/0000-0002-2404-8098




Time dependent advection-reaction-diffusion, Operator splitting, Discontinuous Galerking method, Plant root growth simulation


The time dependant advection-reaction-diffusion equation is used in the C-Root model to simulate root growth. This equation can also be applied in many others applications in life sciences. In this context the unknown is related to densities and one of the important property of the problem is that the solution is non-negative for positive initial conditions. One of the difficulty at the discrete level is to preserve the positivity of the approximated solution during the simulations. In this work we solved the model using Discontinuous Galerkin elements combined with an operator splitting technique. The DG method is briefly presented then we motivated the use of the operator splitting technique by doing some numerical experiments. Those experiments showed that the same time approximation scheme may not be suitable for all the operators of the model. We validated our implementation of the splitting technique in a simple test case. Then we performed a simulation of a plagiotropic root of Eucalyptus.






Original Articles