Geometry of adipocyte packing in subcutaneous tissue contributes to nonlinear tissue properties captured through a Gaussian process surrogate model

Subcutaneous tissue mechanical response is governed by the geometry and mechanical properties at the microscale and drives physiological and clinical processes such as drug delivery. Even though adipocyte packing is known to change with age, disease, and from one individual to another, the link between the geometry of the packing and the overall mechanical response of adipose tissue remains poorly understood. Here we create 1200 periodic representative volume elements (RVEs) that sample the possible space of Laguerre packings describing adipose tissue. RVE mechanics are modeled under tri-axial loading. Equilibrium configuration of RVEs is solved by minimizing an energetic potential that includes volume change contributions from adipocyte expansion, and area change contributions from collagen foam stretching. The resulting mechanical response across all RVE samples is interpolated with the aid of a Gaussian process (GP), revealing how the microscale geometry dictates the overall RVE mechanics. For example, increase in adipocyte size and increase in sphericity lead to adipose tissue softening. We showcase the use of the homogenized model in finite element simulations of drug injection by implementing a Blatz-Ko model, informed by the GP, as a custom material in the popular open-source package FEBio. These simulations show how microscale geometry can lead to vastly different injection dynamics even if the constituent parameters are held constant, highlighting the importance of characterizing individual's adipose tissue structure in the development of personalized therapies. We simulate RVEs, use a GP to interpolate and get insights about the mechanical response as a function of geometry, and incorporate the model into finite element simulations of drug delivery.