A novel strategy to identify the critical conditions for growth-induced instabilities

On skin microrelief and the emergence of expression micro-wrinkles

Over the course of a life time, as a result of adaptive mechanobiological processes (e.g. ageing), or the action of external physical factors such as mechanical loading, the human skin is subjected to, and hosts complex biophysical processes. These phenomena typically operate through a complex interplay, that, ultimately, is responsible for the evolutive geometrical characteristics of the skin surface. Wrinkles are a manifestation of these effects. Although numerous theoretical models of wrinkles arising in multi-layered structures have been proposed, they typically apply to idealised geometries. In the case of skin, which can be viewed as a geometrically complex multi-layer assembly, it is pertinent to question whether the natural skin microrelief could play a significant role in conditioning the characteristics of compression-induced micro-wrinkles by acting as an array of geometrical imperfections. Here, we explore this question through the development of an anatomically-based finite strain parametric finite element model of the skin, represented as a stratum corneum layer on top of a thicker and softer substrate. Our study suggests that skin microrelief could be the dominant factor conditioning micro-wrinkle characteristics for moderate elastic modulus ratios between the two layers. Beyond stiffness ratios of 100, other factors tend to overwrite the effects of skin microrelief. Such stiffness ratio fluctuations can be induced by changes in relative humidity or particular skin conditions and can therefore have important implications for skin tribology.

Mathematical and computational modelling of skin biophysics: a review

The objective of this paper is to provide a review on some aspects of the mathematical and computational modelling of skin biophysics, with special focus on constitutive theories based on nonlinear continuum mechanics from elasticity, through anelasticity, including growth, to thermoelasticity. Microstructural and phenomenological approaches combining imaging techniques are also discussed. Finally, recent research applications on skin wrinkles will be presented to highlight the potential of physics-based modelling of skin in tackling global challenges such as ageing of the population and the associated skin degradation, diseases and traumas.

Wrinkling instabilities in soft bilayered systems

Wrinkling phenomena control the surface morphology of many technical and biological systems. While primary wrinkling has been extensively studied, experimentally, analytically and computationally, higher-order instabilities remain insufficiently understood, especially in systems with stiffness contrasts well below 100. Here, we use the model system of an elastomeric bilayer to experimentally characterize primary and secondary wrinkling at moderate stiffness contrasts. We systematically vary the film thickness and substrate prestretch to explore which parameters modulate the emergence of secondary instabilities, including period-doubling, period-tripling and wrinkle-to-fold transitions. Our experiments suggest that period-doubling is the favourable secondary instability mode and that period-tripling can emerge under disturbed boundary conditions. High substrate prestretch can suppress period-doubling and primary wrinkles immediately transform into folds. We combine analytical models with computational simulations to predict the onset of primary wrinkling, the post-buckling behaviour, secondary bifurcations and the wrinkle-to-fold transition. Understanding the mechanisms of pattern selection and identifying the critical control parameters of wrinkling will allow us to fabricate smart surfaces with tunable properties and to control undesired surface patterns like in the asthmatic airway.This article is part of the themed issue 'Patterning through instabilities in complex media: theory and applications.'

Systems biology and mechanics of growth

In contrast to inert systems, living biological systems have the advantage to adapt to their environment through growth and evolution. This transfiguration is evident during embryonic development, when the predisposed need to grow allows form to follow function. Alterations in the equilibrium state of biological systems breed disease and mutation in response to environmental triggers. The need to characterize the growth of biological systems to better understand these phenomena has motivated the continuum theory of growth and stimulated the development of computational tools in systems biology. Biological growth in development and disease is increasingly studied using the framework of morphoelasticity. Here, we demonstrate the potential for morphoelastic simulations through examples of volume, area, and length growth, inspired by tumor expansion, chronic bronchitis, brain development, intestine formation, plant shape, and myopia. We review the systems biology of living systems in light of biochemical and optical stimuli and classify different types of growth to facilitate the design of growth models for various biological systems within this generic framework. Exploring the systems biology of growth introduces a new venue to control and manipulate embryonic development, disease progression, and clinical intervention.

Patient-Specific Airway Wall Remodeling in Chronic Lung Disease

Chronic lung disease affects more than a quarter of the adult population; yet, the mechanics of the airways are poorly understood. The pathophysiology of chronic lung disease is commonly characterized by mucosal growth and smooth muscle contraction of the airways, which initiate an inward folding of the mucosal layer and progressive airflow obstruction. Since the degree of obstruction is closely correlated with the number of folds, mucosal folding has been extensively studied in idealized circular cross sections. However, airflow obstruction has never been studied in real airway geometries; the behavior of imperfect, non-cylindrical, continuously branching airways remains unknown. Here we model the effects of chronic lung disease using the nonlinear field theories of mechanics supplemented by the theory of finite growth. We perform finite element analysis of patient-specific Y-branch segments created from magnetic resonance images. We demonstrate that the mucosal folding pattern is insensitive to the specific airway geometry, but that it critically depends on the mucosal and submucosal stiffness, thickness, and loading mechanism. Our results suggests that patient-specific airway models with inherent geometric imperfections are more sensitive to obstruction than idealized circular models. Our models help to explain the pathophysiology of airway obstruction in chronic lung disease and hold promise to improve the diagnostics and treatment of asthma, bronchitis, chronic obstructive pulmonary disease, and respiratory failure.

Period-doubling and period-tripling in growing bilayered systems

Growing layers on elastic substrates are capable of creating a wide variety of surface morphologies. Moderate growth generates a regular pattern of sinusoidal wrinkles with a homogeneous energy distribution. While the critical conditions for periodic wrinkling have been extensively studied, the rich pattern formation beyond this first instability point remains poorly understood. Here we show that upon continuing growth, the energy progressively localizes and new complex morphologies emerge. Previous studies have often overlooked these secondary bifurcations; they have focused on large stiffness ratios between layer and substrate, where primary instabilities occur early, long before secondary instabilities emerge. We demonstrate that secondary bifurcations are particularly critical in the low stiffness ratio regime, where the critical conditions for primary and secondary instabilities move closer together. Amongst all possible secondary bifurcations, the mode of period-doubling plays a central role - it is energetically favorable over all other modes. Yet, we can numerically suppress period-doubling, by choosing boundary conditions, which favor alternative higher order modes. Our results suggest that in the low stiffness regime, pattern formation is highly sensitive to small imperfections: surface morphologies emerge rapidly, change spontaneously, and quickly become immensely complex. This is a common paradigm in developmental biology. Our results have significantly applications in the morphogenesis of living systems where growth is progressive and stiffness ratios are low.

On the mechanics of growing thin biological membranes

Despite their seemingly delicate appearance, thin biological membranes fulfill various crucial roles in the human body and can sustain substantial mechanical loads. Unlike engineering structures, biological membranes are able to grow and adapt to changes in their mechanical environment. Finite element modeling of biological growth holds the potential to better understand the interplay of membrane form and function and to reliably predict the effects of disease or medical intervention. However, standard continuum elements typically fail to represent thin biological membranes efficiently, accurately, and robustly. Moreover, continuum models are typically cumbersome to generate from surface-based medical imaging data. Here we propose a computational model for finite membrane growth using a classical midsurface representation compatible with standard shell elements. By assuming elastic incompressibility and membrane-only growth, the model a priori satisfies the zero-normal stress condition. To demonstrate its modular nature, we implement the membrane growth model into the general-purpose non-linear finite element package Abaqus/Standard using the concept of user subroutines. To probe efficiently and robustness, we simulate selected benchmark examples of growing biological membranes under different loading conditions. To demonstrate the clinical potential, we simulate the functional adaptation of a heart valve leaflet in ischemic cardiomyopathy. We believe that our novel approach will be widely applicable to simulate the adaptive chronic growth of thin biological structures including skin membranes, mucous membranes, fetal membranes, tympanic membranes, corneoscleral membranes, and heart valve membranes. Ultimately, our model can be used to identify diseased states, predict disease evolution, and guide the design of interventional or pharmaceutic therapies to arrest or revert disease progression.

On the Role of Mechanics in Chronic Lung Disease

Progressive airflow obstruction is a classical hallmark of chronic lung disease, affecting more than one fourth of the adult population. As the disease progresses, the inner layer of the airway wall grows, folds inwards, and narrows the lumen. The critical failure conditions for airway folding have been studied intensely for idealized circular cross-sections. However, the role of airway branching during this process is unknown. Here, we show that the geometry of the bronchial tree plays a crucial role in chronic airway obstruction and that critical failure conditions vary significantly along a branching airway segment. We perform systematic parametric studies for varying airway cross-sections using a computational model for mucosal thickening based on the theory of finite growth. Our simulations indicate that smaller airways are at a higher risk of narrowing than larger airways and that regions away from a branch narrow more drastically than regions close to a branch. These results agree with clinical observations and could help explain the underlying mechanisms of progressive airway obstruction. Understanding growth-induced instabilities in constrained geometries has immediate biomedical applications beyond asthma and chronic bronchitis in the diagnostics and treatment of chronic gastritis, obstructive sleep apnea and breast cancer.