Gareth Wyn Jones

Modeling Growth in Biological Materials

The biomechanical modeling of growing tissues has recently become an area of intense interest. In particular, the interplay between growth patterns and mechanical stress is of great importance, with possible applications to arterial mechanics, embryo morphogenesis, tumor development, and bone remodeling. This review aims to give an overview of the theories that have been used to model these phenomena, categorized according to whether the tissue is considered as a continuum object or a collection of cells. Among the continuum models discussed is the deformation gradient decomposition method, which allows a residual stress field to develop from an incompatible growth field. The cell-based models are further subdivided into cellular automata, center-dynamics, and vertex-dynamics models. Of these the second two are considered in more detail, especially with regard to their treatment of cell-cell interactions and cell division. The review concludes by assessing the prospects for reconciliation between these two fundamentally different approaches to tissue growth, and by identifying possible avenues for further research.

Design criteria for a printed tissue engineering construct: a mathematical homogenization approach

Cartilage tissue repair procedures currently under development aim to create a construct in which patient-derived cells are seeded and expanded ex vivo before implantation back into the body. The key challenge is producing physiologically realistic constructs that mimic real tissue structure and function. One option with vast potential is to print strands of material in a 3D structure called a scaffold that imitates the real tissue structure; the strands are composed of gel seeded with cells and so provide a template for cartilaginous tissue growth. The scaffold is placed in the construct and pumped with nutrient-rich culture medium to supply nutrients to the cells and remove waste products, thus promoting tissue growth. In this paper we use asymptotic homogenization to determine the effective flow and transport properties of such a printed scaffold system. These properties are used to predict the distribution of nutrient/waste products through the construct, and to specify design criteria for the scaffold that will optimize the growth of functional tissue.

The buckling of capillaries in solid tumours

We develop a model of the buckling (both planar and axial) of capillaries in cancer tumours, using nonlinear solid mechanics. The compressive stress in the tumour interstitium is modelled as a consequence of the rapid proliferation of the tumour cells, using a multiplicative decomposition of the deformation gradient. In turn, the tumour cell proliferation is determined by the oxygen concentration (which is governed by the diffusion equation) and the solid stress. We apply a linear stability analysis to determine the onset of mechanical instability, and the Liapunov–Schmidt reduction to determine the postbuckling behaviour. We find that planar modes usually go unstable before axial modes, so that our model can explain the buckling of capillaries, but not as easily their tortuosity. We also find that the inclusion of anisotropic growth in our model can substantially affect the onset of buckling. Anisotropic growth also results in a feedback effect that substantially affects the magnitude of the buckle.

Designing electronic properties of two-dimensional crystals through optimization of deformations

One of the enticing features common to most of the two-dimensional electronic systems that are currently at the forefront of materials science research is the ability to easily introduce a combination of planar deformations and bending in the system. Since the electronic properties are ultimately determined by the details of atomic orbital overlap, such mechanical manipulations translate into modified electronic properties. Here, we present a general-purpose optimization framework for tailoring physical properties of two-dimensional electronic systems by manipulating the state of local strain, allowing a one-step route from their design to experimental implementation. A definite example, chosen for its relevance in light of current experiments in graphene nanostructures, is the optimization of the experimental parameters that generate a prescribed spatial profile of pseudomagnetic fields in graphene. But the method is general enough to accommodate a multitude of possible experimental parameters and conditions whereby deformations can be imparted to the graphene lattice, and complies, by design, with graphenes elastic equilibrium and elastic compatibility constraints. As a result, it efficiently answers the inverse problem of determining the optimal values of a set of external or control parameters that result in a graphene deformation whose associated pseudomagnetic field profile best matches a prescribed target. The ability to address this inverse problem in an expedited way is one key step for practical implementations of the concept of two-dimensional systems with electronic properties strain-engineered to order. The general-purpose nature of this calculation strategy means that it can be easily applied to the optimization of other relevant physical quantities which directly depend on the local strain field, not just in graphene but in other two-dimensional electronic membranes.

Modelling apical constriction in epithelia using elastic shell theory

Apical constriction is one of the fundamental mechanisms by which embryonic tissue is deformed, giving rise to the shape and form of the fully-developed organism. The mechanism involves a contraction of fibres embedded in the apical side of epithelial tissues, leading to an invagination or folding of the cell sheet. In this article the phenomenon is modelled mechanically by describing the epithelial sheet as an elastic shell, which contains a surface representing the continuous mesh formed from the embedded fibres. Allowing this mesh to contract, an enhanced shell theory is developed in which the stiffness and bending tensors of the shell are modified to include the fibres’ stiffness, and in which the active effects of the contraction appear as body forces in the shell equilibrium equations. Numerical examples are presented at the end, including the bending of a plate and a cylindrical shell (modelling neurulation) and the invagination of a spherical shell (modelling simple gastrulation).

Quantized transport, strain-induced perfectly conducting modes and valley filtering on shape-optimized graphene Corbino devices

The extreme mechanical resilience of graphene and the peculiar coupling it hosts between lattice and electronic degrees of freedom have spawned a strong impetus toward strain-engineered graphene where, on the one hand, strain augments the richness of its phenomenology and makes possible new concepts for electronic devices, and on the other hand, new and extreme physics might take place. Here, we demonstrate that the shape of substrates supporting graphene sheets can be optimized for approachable experiments where strain-induced pseudomagnetic fields (PMF) can be tailored by pressure for directionally selective electronic transmission and pinching-off of current flow down to the quantum channel limit. The Corbino-type layout explored here furthermore allows filtering of charge carriers according to valley and current direction, which can be used to inject or collect valley-polarized currents, thus realizing one of the basic elements required for valleytronics. Our results are based on a framework developed to realistically determine the combination of strain, external parameters, and geometry optimally compatible with the target spatial profile of a desired physical property-the PMF in this case. Characteristic conductance profiles are analyzed through quantum transport calculations on large graphene devices having the optimal shape.

Asymptotic analysis of a buckling problem for an embedded spherical shell

The axisymmetric buckling of a spherical shell embedded in an elastic medium with uniaxial compression at infinity is examined in the limit of small shell thickness ratio. An asymptotic method is developed by considering the paradigm problem of a beam attached to a Winkler substrate of variable stiffness, which in the small aspect ratio limit displays the same behavior as the shell. The asymptotic method is then applied to the Euler–Lagrange equations corresponding to shell buckling. The system is analyzed in two distinguished limits, displaying good agreement with the full numerical results.

Grasping with a soft glove: intrinsic impedance control in pneumatic actuators

The interaction of a robotic manipulator with unknown soft objects represents a significant challenge for traditional robotic platforms because of the difficulty in controlling the grasping force between a soft object and a stiff manipulator. Soft robotic actuators inspired by elephant trunks, octopus limbs and muscular hydrostats are suggestive of ways to overcome this fundamental difficulty. In particular, the large intrinsic compliance of soft manipulators such as ‘pneu-nets’—pneumatically actuated elastomeric structures—makes them ideal for applications that require interactions with an uncertain mechanical and geometrical environment. Using a simple theoretical model, we show how the geometric and material nonlinearities inherent in the passive mechanical response of such devices can be used to grasp soft objects using force control, and stiff objects using position control, without any need for active sensing or feedback control. Our study is suggestive of a general principle for designing actuators with autonomous intrinsic impedance control.

Planar morphometry, shear and optimal quasi-conformal mappings

To characterize the diversity of planar shapes in such instances as insect wings and plant leaves, we present a method for the generation of a smooth morphometric mapping between two planar domains which matches a number of homologous points. Our approach tries to balance the competing requirements of a descriptive theory which may not reflect mechanism and a multi-parameter predictive theory that may not be well constrained by experimental data. Specifically, we focus on aspects of shape as characterized by local rotation and shear, quantified using quasi-conformal maps that are defined precisely in terms of these fields. To make our choice optimal, we impose the condition that the maps vary as slowly as possible across the domain, minimizing their integrated squared-gradient. We implement this algorithm numerically using a variational principle that optimizes the coefficients of the quasi-conformal map between the two regions and show results for the recreation of a sample historical grid deformation mapping of D’Arcy Thompson. We also deploy our method to compare a variety of Drosophila wing shapes and show that our approach allows us to recover aspects of phylogeny as marked by morphology.

Optimal control of plates using incompatible strains

A flat plate will bend into a curved shell if it experiences an inhomogeneous growth field or if constrained appropriately at a boundary. While the forward problem associated with this process is well studied, the inverse problem of designing the boundary conditions or growth fields to achieve a particular shape is much less understood. We use ideas from variational optimization theory to formulate a well posed version of this inverse problem to determine the optimal growth field or boundary condition that will give rise to an arbitrary target shape, optimizing for both closeness to the target shape and for smoothness of the growth field. We solve the resulting system of PDE numerically using finite element methods with examples for both the fully non-symmetric case as well as for simplified one-dimensional and axisymmetric geometries. We also show that the system can also be solved semi-analytically by positing an ansatz for the deformation and growth fields in a circular disk with given thickness profile, leading to paraboloidal, cylindrical and saddle-shaped target shapes, and show how a soft mode can arise from a non-axisymmetric deformation of a structure with axisymmetric material properties.