Sebastian Aland

Wetting Resistance at Its Topographical Limit: The Benefit of Mushroom and Serif T Structures

Springtails (Collembola) are wingless arthropods adapted to cutaneous respiration in temporarily rain-flooded habitats. They immediately form a plastron, protecting them against suffocation upon immersion into water and even low-surface-tension liquids such as alkanes. Recent experimental studies revealed a high-pressure resistance of such plastrons against collapse. In this work, skin sections of Orthonychiurus stachianus are studied by transmission electron microscopy. The micrographs reveal cavity side-wall profiles with characteristic overhangs. These were fitted by polynomials to allow access for analytical and numerical calculations of the breakthrough pressure, that is, the barrier against plastron collapse. Furthermore, model profiles with well-defined geometries were used to set the obtained results into context and to develop a general design principle for the most robust surface structures. Our results indicate the decisive role of the sectional profile of overhanging structures to form a robust heterogeneous wetting state for low-surface-tension liquids that enables the omniphobicity. Furthermore, the design principles of mushroom and serif T structures pave the way for omniphobic surfaces with a high-pressure resistance irrespective of solid surface chemistry.

Extracting Cell Stiffness from Real-Time Deformability Cytometry: Theory and Experiment

Cell stiffness is a sensitive indicator of physiological and pathological changes in cells, with many potential applications in biology and medicine. A new method, real-time deformability cytometry, probes cell stiffness at high throughput by exposing cells to a shear flow in a microfluidic channel, allowing for mechanical phenotyping based on single-cell deformability. However, observed deformations of cells in the channel not only are determined by cell stiffness, but also depend on cell size relative to channel size. Here, we disentangle mutual contributions of cell size and cell stiffness to cell deformation by a theoretical analysis in terms of hydrodynamics and linear elasticity theory. Performing real-time deformability cytometry experiments on both model spheres of known elasticity and biological cells, we demonstrate that our analytical model not only predicts deformed shapes inside the channel but also allows for quantification of cell mechanical parameters. Thereby, fast and quantitative mechanical sampling of large cell populations becomes feasible.

A continuum model of colloid-stabilized interfaces

Colloids that are partially wetted by two immiscible fluids can become confined to fluid-fluidinterfaces. At sufficiently high volume fractions, the colloids may jam and the interface may crystallize. Examples include bicontinuous interfacially jammed emulsion gels (bijels), which were proposed in this study by Stratford et al. [Science 309, 2198 (2005)] as a hypothetical new class of soft materials in which interpenetrating, continuous domains of two immiscible viscous fluids are maintained in a rigid state by a jammed layer of colloidal particles at their interface. We develop a continuum model for such a system that is capable of simulating the long-time evolution. A Navier-Stokes-Cahn-Hilliard model for the macroscopic two-phase flow system is combined with a surface phase-field-crystal model for the microscopic colloidal system along the interface. The presence of colloids introduces elastic forces at the interface between the two immiscible fluid phases. An adaptive finite element method is used to solve the model numerically. Using a variety of flow configurations in two dimensions, we demonstrate that as colloids jam on the interface and the interfacecrystallizes, the elastic force may be strong enough to make the interface sufficiently rigid to resist external forces, such as an applied shear flow, as well as surface tension induced coarsening in bicontinuous structures.

Diffuse interface models of locally inextensible vesicles in a viscous fluid

We present a new diffuse interface model for the dynamics of inextensible vesicles in a viscous fluid with inertial forces. A new feature of this work is the implementation of the local inextensibility condition in the diffuse interface context. Local inextensibility is enforced by using a local Lagrange multiplier, which provides the necessary tension force at the interface. We introduce a new equation for the local Lagrange multiplier whose solution essentially provides a harmonic extension of the multiplier off the interface while maintaining the local inextensibility constraint near the interface. We also develop a local relaxation scheme that dynamically corrects local stretching/compression errors thereby preventing their accumulation. Asymptotic analysis is presented that shows that our new system converges to a relaxed version of the inextensible sharp interface model. This is also verified numerically. To solve the equations, we use an adaptive finite element method with implicit coupling between the Navier-Stokes and the diffuse interface inextensibility equations. Numerical simulations of a single vesicle in a shear flow at different Reynolds numbers demonstrate that errors in enforcing local inextensibility may accumulate and lead to large differences in the dynamics in the tumbling regime and smaller differences in the inclination angle of vesicles in the tank-treading regime. The local relaxation algorithm is shown to prevent the accumulation of stretching and compression errors very effectively. Simulations of two vesicles in an extensional flow show that local inextensibility plays an important role when vesicles are in close proximity by inhibiting fluid drainage in the near contact region.

Time integration for diffuse interface models for two-phase flow

We propose a variant of the @q-scheme for diffuse interface models for two-phase flow, together with three new linearization techniques for the surface tension. These involve either additional stabilizing force terms, or a fully implicit coupling of the Navier-Stokes and Cahn-Hilliard equation. In the common case that the equations for interface and flow are coupled explicitly, we find a time step restriction which is very different to other two-phase flow models and in particular is independent of the grid size. We also show that the proposed stabilization techniques can lift this time step restriction. Even more pronounced is the performance of the proposed fully implicit scheme which is stable for arbitrarily large time steps. We demonstrate in a Taylor-flow application that this superior coupling between flow and interface equation can decrease the computation time by several orders of magnitude.

Buckling Instability of Viral Capsids—A Continuum Approach

The crystallographic structure of spherical viruses is modeled using a multiscale approach combining a macroscopic Helfrich model for morphology evolution with a microscopic approximation of a classical density functional theory for the protein interactions. The derivation of the model is based on energy dissipation and conservation of protein number density. The resulting set of equations is solved within a diffuse domain approach using finite elements and shows buckling transitions of spherical shapes into faceted viral shapes.

A Mechanistic Collective Cell Model for Epithelial Colony Growth and Contact Inhibition

We present a mechanistic hybrid continuum-discrete model to simulate the dynamics of epithelial cell colonies. Collective cell dynamics are modeled using continuum equations that capture plastic, viscoelastic, and elastic deformations in the clusters while providing single-cell resolution. The continuum equations can be viewed as a coarse-grained version of previously developed discrete models that treat epithelial clusters as a two-dimensional network of vertices or stochastic interacting particles and follow the framework of dynamic density functional theory appropriately modified to account for cell size and shape variability. The discrete component of the model implements cell division and thus influences cell size and shape that couple to the continuum component. The model is validated against recent in vitro studies of epithelial cell colonies using Madin-Darby canine kidney cells. In good agreement with experiments, we find that mechanical interactions and constraints on the local expansion of cell size cause inhibition of cell motion and reductive cell division. This leads to successively smaller cells and a transition from exponential to quadratic growth of the colony that is associated with a constant-thickness rim of growing cells at the cluster edge, as well as the emergence of short-range ordering and solid-like behavior. A detailed analysis of the model reveals a scale invariance of the growth and provides insight into the generation of stresses and their influence on the dynamics of the colonies. Compared to previous models, our approach has several advantages: it is independent of dimension, it can be parameterized using classical elastic properties (Poissons ratio and Youngs modulus), and it can easily be extended to incorporate multiple cell types and general substrate geometries.

Adaptive diffuse domain approach for calculating mechanically induced deformation of trabecular bone

Remodelling of trabecular bone is essentially affected by the mechanical load of the trabeculae. Mathematical modelling and simulation of the remodelling process have to include time-consuming calculations of the displacement field within the complex trabecular structure under loading. We present an adaptive diffuse domain approach for calculating the elastic bone deformation based on micro computer tomogram data of real trabecular bone structures and compared it with a conventional voxel-based finite element method. In addition to allowing for higher computational efficiency, the adaptive approach is characterised by a very smooth representation of the bone surface, which suggests that this approach would be suitable as a basis for future simulations of bone resorption and formation processes within the trabecular structure.

Phase Field Models for Two-Phase Flow with Surfactants and Biomembranes

We give an overview on recent developments of phase field models for two-phase flows with surfactants and lipid bilayer membranes. Starting from the two-phase flow model of a clean fluid-fluid interface we discuss the time discretization and boundary conditions for dynamic and static contact angles. Using the adsorption models of Henry and Langmuir, soluble surfactants are included in the diffuse interface formulation. To consider lipid bilayer membranes the model is extended by membrane bending stiffness and membrane inextensibility. We present phase field models to include these elastic effects, with a particular focus on the inextensibility constraint for which we discuss different phase field variants from the literature and present numerical tests.

A phase-field model for fluid–structure interaction

Abstract In this paper, we develop a novel phase-field model for fluid–structure interaction (FSI), that is capable to handle very large deformations as well as topology changes like contact of the solid to a wall. The model is based on a fully Eulerian description of the velocity field in both, the fluid and the elastic domain. Viscous and elastic stresses in the Navier–Stokes equations are restricted to the corresponding domains by multiplication with their characteristic functions. The solid is described as a hyperelastic neo-Hookean material and the elastic stress is obtained by solving an additional Oldroyd-B – like equation. Thermodynamically consistent forces are derived by energy variation. The convergence of the derived equations to the traditional sharp interface formulation of fluid–structure interaction is shown by matched asymptotic analysis. The model is evaluated in a challenging benchmark scenario of an elastic body traversing a fluid channel. A comparison to reference values from Arbitrary Lagrangian Eulerian (ALE) simulations shows very good agreement. We highlight some distinct advantages of the new model, like the avoidance of re-triangulations and the stable inclusion of surface tension. Further, we demonstrate how simple it is to include contact dynamics into the model, by simulating a ball bouncing off a wall. We extend this scenario to include adhesion of the ball, which to our knowledge, cannot be simulated with any other FSI model. While we have restricted simulations to fluid–structure interaction, the model is capable to simulate any combination of viscous fluids, visco-elastic fluids and elastic solids.