Identification of parameters of the advection-reaction-diffusion equation solved with Discontinuous Galerkin methods: application to plant root growth

It is well known that plant root systems are difficult to study in situ. Modelling is not only useful to study and simulate root growth, but also to design new experiments and help to decide what kind of data should be appropriate or not. The density modelling approach consists in following the evolution in time and space of the root density in a volume of soil. This leads to an advection-reaction-diffusion like-equation or to a system. Up to now, there is no specific field experiments to directly measure or estimate some of the model parameters. That is why, we rely on field root density data to identify those parameters by solving an inverse problem, a nonlinear least square minimization problem. The forward problem is approximated on an unstructured mesh. Since it is well known that Lagrange finite elements suffer from a lack of stability especially for a dominant advection, and motivated by the treatment of discontinuous parameters (e.g. for stratified soils), we implemented a Discontinuous Galerkin (DG) method for the space approximation. The inverse problem is solved using a Levenberg-Marquardt method. In this presentation we intend to show some simulations of root growth. Then we address the problem of identifying the parameters and we present some numerical experiments to emphasize that this problem needs a particular attention. (Texte intégral)

Saved in:
Bibliographic Details
Main Author: Peynaud, Emilie
Format: conference_item biblioteca
Language:eng
Published: Université de Perpignan
Subjects:F62 - Physiologie végétale - Croissance et développement, U10 - Informatique, mathématiques et statistiques,
Online Access:http://agritrop.cirad.fr/581132/
http://agritrop.cirad.fr/581132/1/ETAMM2016_abstract-MS9-peynaud_YD.pdf
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:It is well known that plant root systems are difficult to study in situ. Modelling is not only useful to study and simulate root growth, but also to design new experiments and help to decide what kind of data should be appropriate or not. The density modelling approach consists in following the evolution in time and space of the root density in a volume of soil. This leads to an advection-reaction-diffusion like-equation or to a system. Up to now, there is no specific field experiments to directly measure or estimate some of the model parameters. That is why, we rely on field root density data to identify those parameters by solving an inverse problem, a nonlinear least square minimization problem. The forward problem is approximated on an unstructured mesh. Since it is well known that Lagrange finite elements suffer from a lack of stability especially for a dominant advection, and motivated by the treatment of discontinuous parameters (e.g. for stratified soils), we implemented a Discontinuous Galerkin (DG) method for the space approximation. The inverse problem is solved using a Levenberg-Marquardt method. In this presentation we intend to show some simulations of root growth. Then we address the problem of identifying the parameters and we present some numerical experiments to emphasize that this problem needs a particular attention. (Texte intégral)