A 3D-1D virtual element method for modeling root water uptake

An optimization-based strategy is proposed for coupling three-dimensional and one-dimensional problems (3D-1D coupling) in the context of soil-root interaction simulations. This strategy, originally designed to tackle generic 3D-1D coupled problems with discontinuous solutions, is here extended to the case of non-linear problems and applied, for the first time, along with a virtual element discretization of the 3D soil sample. This further enhances the capability of the method to handle geometrical complexities, allowing to easily mesh domains characterized, for instance, by the presence of stones and other impervious obstacles of arbitrary shape. A discrete-hybrid tip-tracking strategy is adopted to model both the root growth and the evolution in time of the water flux, the pressure head and the water content, both in the roots and in the surrounding soil sample. By choosing proper rules for the generation of branches, realistic root-network configurations are obtained. Several numerical examples are proposed, proving both the accuracy of the adopted method and its applicability in realistic and large scale simulations.