Search

We present a comprehensive theoretical and computational model that explores the behavior of a thin hydrated film bonded to a non-hydrated / impermeable soft substrate in the context of surface and bulk elasticity coupled with surface diffusion kinetics. This type of coupling can manifests as an integral aspect in diverse engineering processes encountered in optical interference coatings, tissue engineering, soft electronics, and can prove important in design process for the next generation of sensors and actuators, especially as the focus is shifted to systems in smaller lengthscales. The intricate interplay between solvent diffusion and deformation of the film is governed by surface poroelasticity, and the viscoelastic deformation of the substrate. While existing methodologies offer tools for studying coupled poroelasticity involving solvent diffusion and network deformation, there exists a gap in understanding how coupled poroelastic processes occurring in a film attached to the boundary of a highly deformable solid can influence its response. In this study, we introduce a non-equilibrium thermodynamics formulation encompassing the multiphysical processes of surface poroelasticity and bulk viscoelasticity, complemented by a corresponding finite element implementation. Our approach captures the complex dynamics between the finite deformation of the substrate and solvent diffusion on the surface. This work contributes valuable insights, particularly in scenarios where the coupling of surface diffusion kinetics and substrate elasticity is an important design factor., Comment: 20 pages, 8 figues. arXiv admin note: text overlap with arXiv:2305.08805

This paper develops and benchmarks an immersed peridynamics method to simulate the deformation, damage, and failure of hyperelastic materials within a fluid-structure interaction framework. The immersed peridynamics method describes an incompressible structure immersed in a viscous incompressible fluid. It expresses the momentum equation and incompressibility constraint in Eulerian form, and it describes the structural motion and resultant forces in Lagrangian form. Coupling between Eulerian and Lagrangian variables is achieved by integral transforms with Dirac delta function kernels, as in standard immersed boundary methods. The major difference between our approach and conventional immersed boundary methods is that we use peridynamics, instead of classical continuum mechanics, to determine the structural forces. We focus on non-ordinary state-based peridynamic material descriptions that allow us to use a constitutive correspondence framework that can leverage well characterized nonlinear constitutive models of soft materials. The convergence and accuracy of our approach are compared to both conventional and immersed finite element methods using widely used benchmark problems of nonlinear incompressible elasticity. We demonstrate that the immersed peridynamics method yields comparable accuracy with similar numbers of structural degrees of freedom for several choices of the size of the peridynamic horizon. We also demonstrate that the method can generate grid-converged simulations of fluid-driven material damage growth, crack formation and propagation, and rupture under large deformations.