A Comparative Analysis of Transient Finite-Strain Coupled Diffusion-Deformation Theories for Hydrogels

Jorge-Humberto Urrea-Quintero,Michele Marino,Thomas Wick,Udo Nackenhorst
DOI: https://doi.org/10.1007/s11831-024-10101-x
IF: 9.7
2024-05-02
Archives of Computational Methods in Engineering
Abstract:This work presents a comparative review and classification between some well-known thermodynamically consistent models of hydrogel behavior in a large deformation setting, specifically focusing on solvent absorption/desorption and its impact on mechanical deformation and network swelling. The proposed discussion addresses formulation aspects, general mathematical classification of the governing equations, and numerical implementation issues based on the finite element method. The theories are presented in a unified framework demonstrating that, despite not being evident in some cases, all of them follow equivalent thermodynamic arguments. A detailed computational analysis is carried out where Taylor–Hood elements are employed in the spatial discretization to satisfy the inf-sup condition and to prevent spurious numerical oscillations. The resulting discrete problems are solved using the FEniCS platform through consistent variational formulations, employing both monolithic and staggered approaches. We conduct benchmark tests on various hydrogel structures, demonstrating that major differences arise from the chosen volumetric response of the hydrogel. The significance of this choice is frequently underestimated in the state-of-the-art literature but has been shown to have substantial implications on the resulting hydrogel behavior.
computer science, interdisciplinary applications,engineering, multidisciplinary,mathematics
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is the comparison and analysis of the diffusion - deformation coupling theory of hydrogels under large - deformation conditions. Specifically, the authors aim to explore the differences between different theoretical models and their impacts on the behavior of hydrogels, especially in terms of mechanical deformation and network swelling during the solvent absorption/desorption process. ### Main Problems 1. **Consistency of Theoretical Models**: The paper shows the equivalence between different theoretical models through a unified framework, although these models may seem not obvious in some cases. The authors hope to reveal that all models follow the same thermodynamic principles through this comparison. 2. **Choice of Volume Response**: The paper emphasizes the significant impact of choosing different volume responses on the behavior of hydrogels. This choice is often overlooked, but it actually has a significant impact on the results. 3. **Numerical Implementation**: The paper also discusses the numerical implementation issues based on the finite element method (FEM), especially how to satisfy the inf - sup condition in the discrete problem to prevent numerical oscillation. ### Specific Objectives 1. **Unified Formulas and Symbols**: Provide unified formulas and symbols for some representative diffusion - deformation theories to facilitate better understanding and comparison. 2. **Mathematical Classification and Numerical Algorithms**: Conduct a mathematical classification of single - block and staggered numerical algorithms and perform a computational comparison in the variational framework. 3. **Numerical Simulation**: Study the capabilities of these theories and the robustness of FEM implementation through one - dimensional to three - dimensional numerical simulations. 4. **Comparison of Benchmark Problems**: Compare the main features of different constitutive models in a unified benchmark problem. ### Research Background Hydrogels are a class of three - dimensional polymer network materials that can absorb and retain a large amount of liquid and are widely used in biomedical, agricultural, and engineering fields. Their properties can be finely tuned by adjusting the chemical composition, cross - linking density, and liquid content. Therefore, understanding the diffusion - deformation mechanism of hydrogels under large - deformation conditions is crucial for optimizing their applications. ### Methods The paper mainly adopts two methods to describe these phenomena: - **Hybrid Theory**: Treat the porous medium as interacting layers superimposed in space, which involves complex model parameters and constitutive choices and is difficult to calibrate in practical applications. - **Macro - scale Porous Elasticity Theory**: Treat the medium as a homogeneous material, described by a coupled deformation - diffusion response, such as the classical framework proposed by Biot et al. ### Conclusions Through detailed mathematical derivations and numerical simulations, the paper reveals the similarities and differences between different theoretical models and emphasizes the importance of choosing an appropriate volume response. This provides a solid theoretical foundation and technical support for future research. ### Formulas Involved The formulas involved in the paper include but are not limited to: - **Deformation Gradient**: \[ F=\text{Grad}(x)=I + \text{Grad}(u) \] - **Right Cauchy - Green Tensor**: \[ C = F^T F \] - **Left Cauchy - Green Tensor**: \[ b = F F^T \] - **First Invariant**: \[ I_1(F)=\text{tr}(C)=\text{tr}(F^T F) \] - **Volume Change**: \[ J_f=1+\Omega c_R \] - **Polymer Volume Fraction**: \[ \phi=(1+\Omega c_R)^{-1} \] - **Free Energy Function**: \[ \Psi_R=\Psi_{\text{mix}}(c_R)+\Psi_{\text{mech}}(F, J) \] These formulas ensure the accuracy and readability of the mathematical expressions in the paper.