Morphoelastic Growth in the Presence of Nutrients at Small Strains

Abstract: A model for morphoelastic growth, that is, growth influenced by elastic stress, driven by the absorption of nutrients is considered. The model features a multiplicative decomposition of the deformation gradient into an elastic contribution and a growth tensor. While the evolution of the system is governed by an ordinary differential equation for the growth tensor on a suitable Banach space, which depends on the elastic stresses and the concentration of a nutrient field, the total deformation is given by the solution of a quasi-static equilibrium equation arising from the formal Euler-Lagrange equations of a hyperelastic variational integral. The nutrient concentration is determined by a linear elliptic reaction-diffusion equation which is formulated in Lagrangian coordinates and whose coefficients depend on the growth tensor as well as the deformation gradient accounting for the change of material properties due to elastic deformation and growth. Existence and uniqueness of solutions of this fully coupled system of differential equations is proven via a fixed-point argument.

• The paper presents a rigorous continuum mechanics framework using multiplicative decomposition to model stress-driven, nutrient-modulated growth at small strains.
• It proves local existence and uniqueness for the coupled quasi-static elastic equilibrium, nutrient diffusion, and growth evolution equations using fixed-point and Schauder theory.
• The study demonstrates numerical stability by quantifying the Lipschitz dependence of deformation and nutrient fields on the growth tensor, offering a benchmark for tissue growth simulations.

Morphoelastic Growth in the Presence of Nutrients at Small Strains

Model Formulation

The paper presents a rigorous continuum mechanical framework for morphoelastic growth at small strains, focusing on the interplay between elastic stresses and nutrient-driven growth. The core modeling assumption is the multiplicative decomposition of the deformation gradient $\nabla y = F_{el} G$, where $F_{el}$ is the elastic gradient and $G$ is the growth tensor. This decomposition, rooted in crystal plasticity, is leveraged here without the regularization generally required in large-strain regimes, capitalizing on the small-strain approximation.

Temporal scale separation is assumed: mechanical equilibration occurs rapidly, nutrient diffusion and absorption are intermediate, while growth is slowest. Accordingly, inertial effects are omitted and the system is treated as quasi-static at each step. The overall governing system comprises:

• A quasi-static equilibrium equation for elastic deformation, formulated via the first Piola-Kirchhoff stress,
• A linear elliptic reaction-diffusion equation for the nutrient field, with coefficients (diffusion, reaction rates) explicitly depending on both $G$ and $\nabla y$, capturing changes in material properties due to deformation and growth,
• An evolution ODE for $G$, incorporating local stress and nutrient concentrations.

An essential aspect is the impossibility, in general, of defining a global intermediate configuration due to the incompatibility of $G$, especially in higher dimensions. The analysis, however, is localized, focusing on neighborhoods of the identity in matrix space.

Mathematical Analysis

The main theorem proves existence and uniqueness of local-in-time solutions to this coupled system under small-strain assumptions. The equations are constructed on Banach spaces ($C^{1,\mu}$ and $C^{2,\mu}$ H\"older spaces), and rely on classical assumptions for hyperelasticity: regularity, coercivity, and frame indifference.

Existence is established through several steps:

1. Quasi-static Elastic Equilibrium: The equilibrium equation admits a unique solution for $y$, via Banach's fixed-point theorem, after linearization around the identity and application of a generalized Korn inequality. No convexity of the energy is assumed, only coercivity and suitable growth conditions.
2. Nutrient Field: Given $F_{el}$0 and $F_{el}$1, the nutrient equation is a well-posed elliptic PDE with coefficients satisfying regularity and uniform ellipticity. Classical Schauder theory provides existence, regularity, and continuous dependence on $F_{el}$2 and $F_{el}$3.
3. Growth Tensor Evolution: Local existence (in time) of solutions to the ODE for $F_{el}$4 is guaranteed by the Picard-Lindel\"of theorem, contingent on the Lipschitz continuity with respect to $F_{el}$5 of both $F_{el}$6 and $F_{el}$7. This property is rigorously established using the implicit function theorem on Banach spaces.

All mappings (from $F_{el}$8 to $F_{el}$9, and $G$0 to $G$1) are shown to possess strong regularity and are locally Lipschitz continuous, supporting the fixed-point argument and time-regularity essential for the ODE analysis.

Numerical Results and Claims

While the paper is primarily analytical, it includes specific modeling examples compliant with the assumptions. For instance, the elastic energy $G$2 is cited as a concrete functional satisfying frame indifference and coercivity.

Non-negativity of nutrients is proven robustly via the maximum principle and Hopf's lemma, under physically appropriate boundary conditions.

Lipschitz dependence of the deformation and nutrient field on the growth tensor is explicitly quantified—this is a critical property for numerical simulation stability and sensitivity analysis.

Boundary condition generality: The framework accommodates time-dependent loads and boundary conditions, as long as they are sufficiently small, demonstrating versatility for practical modeling.

Implications and Theoretical Contributions

The study establishes foundational existence and uniqueness results for morphoelastic growth models in the small-strain regime, with nutrient coupling. Unlike large-strain models—which require regularization or restrictive structural conditions—this formulation is fully nonlinear but avoids such constraints due to the small-strain setup. The results provide mathematical justification for several phenomenological models previously used in biomechanics and developmental biology.

The theoretical analysis highlights that:

• For small strains, detailed modeling of the stress-driven, nutrient-modulated growth can be mathematically justified and numerically tractable.
• The key technical obstacle—the lack of a global intermediate configuration—is mitigated locally.
• The structure of the ODE for $G$3 allows for blow-up in finite time unless additional damping mechanisms (e.g., exponential growth laws) are imposed.

Practical and Future Directions

Practically, this framework underpins simulation and analysis of early-stage growth processes in biological tissues, where strain remains modest. The explicit coupling with nutrient fields opens avenues for more realistic models of growth modulation in engineering and biomedical applications.

Possible future developments include:

• Extension to large-strain settings, though this presents considerable analytical challenges.
• Incorporation of more complex boundary conditions and external stimuli.
• Exploration of numerical algorithms for coupled PDE/ODE systems, with provable stability leveraging the Lipschitz estimates.
• Investigation of the impact of growth incompatibility in anisotropic materials and microstructure-informed modeling.

The small-strain theory may serve as a benchmark for computational schemes and for validating experimental observations in tissue engineering, plant biomechanics, and related domains.

Conclusion

This paper provides a rigorous, coupled morphoelastic growth theory—including nutrient dynamics and stress effects—anchored in small-strain continuum mechanics. By establishing existence, uniqueness, and strong regularity results for the fully nonlinear, coupled system, the work offers a mathematically sound foundation for modeling nutrient-driven growth in elastic media at early development stages. The implications extend to both theoretical biomechanics and practical simulation, with clear indications for further research in both directions.