Publications des agents du Cirad


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

Peynaud E.. 2018. Biomath, 7 (2) : 19 p..

DOI: 10.11145/j.biomath.2018.12.037

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.

Mots-clés : eucalyptus; racine; croissance; modèle de simulation; mathématique

Documents associés

Article (b-revue à comité de lecture)

Agents Cirad, auteurs de cette publication :