Growth models for tree stems and vines

When obstacles are present, the model takes the form of a differential inclusion with state constraints. At each time t, a cone of admissible reactions is determined by the minimization of an elastic deformation energy. The main theorem shows that local solutions exist and can be prolonged globally in time, except when a specific “breakdown configuration” is reached. Approximate solutions are constructed by an operator-splitting technique. Some numerical simulations are provided at the end of the paper.