Search

Understanding the generation of mechanical stress in drying, particle-laden films is important for a wide range of industrial processes. The cantilever experiment allows the stress in a drying film that has been deposited onto a thin plate to be quantified. Mechanical stresses in the film are transmitted to the plate and drive bending. Mathematical modelling enables the film stress to be inferred from measurements of the plate deflection. The aim of this paper is to present simplified models of the cantilever experiment that have been derived from the time-dependent equations of continuum mechanics using asymptotic methods. The film is described using nonlinear poroelasticity and the plate using nonlinear elasticity. In contrast to Stoney-like formulae, the simplified models account for films with non-uniform thickness and stress. The film model reduces to a single differential equation that can be solved independently of the plate equations. The plate model reduces to an extended form of the Foppl-von Karman (FvK) equations that accounts for gradients in the longitudinal traction acting on the plate surface. Consistent boundary conditions for the FvK equations are derived by resolving the Saint-Venant boundary layers at the free edges of the plate. The asymptotically reduced models are in excellent agreement with finite element solutions of the full governing equations. As the P\'eclet number increases, the time evolution of the plate deflection changes from $t$ to $t^{1/2}$, in agreement with experiments.

Drug-loaded hydrogels provide a means to deliver pharmaceutical agents to specific sites within the body at a controlled rate. The aim of this paper is to understand how controlled drug release can be achieved by tuning the initial distribution of drug molecules in a hydrogel. A mathematical model is presented for a spherical drug-loaded hydrogel. The model captures the nonlinear elasticity of the polymer network and thermodynamics of swelling. By assuming that the drug molecules are dilute, the equations for hydrogel swelling and drug transport partially decouple. A fast optimisation method is developed to accurately compute the optimal initial drug concentration by minimising the error between the numerical drug-release profile and a target profile. By taking the target drug efflux to be piecewise constant, the optimal initial configuration consists of a central drug-loaded core with isolated drug packets near the free boundary of the hydrogel. The optimal initial drug concentration is highly effective at mitigating the burst effect, where a large amount of drug is rapidly released into the environment. The hydrogel stiffness can be used to further tune the rate of drug release. Although stiffer gels lead to less swelling and hence reduce the drug diffusivity, the drug-release kinetics are faster than for soft gels due to the decreased distance that drug molecules must travel to reach the free surface.

The electric double layer (EDL) that forms at the interface between a polyelectrolyte gel and a salt bath is studied using asymptotic and numerical methods. Specifically, matched asymptotic expansions, based on the smallness of the Debye length relative to the typical gel dimensions, are used to construct solutions of the governing equations and derive electroneutral models with consistent jump conditions across the gel-bath interface. A general approach for solving the equations of incompressible nonlinear elasticity in a curved boundary layer is developed and used to resolve the gel mechanics in the EDL. A critical feature of the model is that it accounts for phase separation within the gel, which gives rise to diffuse interfaces with a characteristic thickness described by the Kuhn length. We show that the solutions of the electroneutral model can only be asymptotically matched to the solutions in the EDL, in general, when the Kuhn length greatly exceeds the Debye length. Conversely, if the Debye length is similar to or larger than the Kuhn length, then the entire gel can self-organise into periodic, electrically charged domains via phase separation. The breakdown of electroneutrality demonstrates that the commonly invoked electroneutral assumption must be used with caution, as it generally only applies when the Debye length is much smaller than the Kuhn length.

Various biological and chemical processes lead to the nucleation and growth of non-wetting fluid bubbles within the pore space of a granular medium, such as the formation of gas bubbles in liquid-saturated lake-bed sediments. In sufficiently soft porous materials, the non-wetting nature of these bubbles can result in the formation of open cavities within the granular solid skeleton. Here, we consider this process through the lens of phase separation, where thermomechanics govern the separation of the non-wetting phase from a fluid-fluid-solid mixture. We construct a phase-field model informed by large-deformation poromechanics, in which two immiscible fluids interact with a poroelastic solid skeleton. Our model captures the competing effects of elasticity and fluid-fluid-solid interactions. We use a phase-field damage model to capture the mechanics of the granular solid. As a model problem, we consider an initial distribution of non-wetting fluid in the pore space that separates into multiple cavities. We use simulations and linear-stability analysis to identify the key parameters that control phase separation, the conditions that favour the formation of cavities, and the characteristic size of the resulting cavities.

We develop a theory for drying-induced stresses in sessile, poroelastic drops undergoing evaporation on rigid surfaces. Using a lubrication-like approximation, the governing equations of three-dimensional nonlinear poroelasticity are reduced to a single thin-film equation for the drop thickness. We find that thin drops experience compressive elastic stresses but the total in-plane stresses are tensile. The mechanical response of the drop is dictated by the initial profile of the solid skeleton, which controls the in-plane deformation, the dominant components of elastic stress, and sets a limit on the depth of delamination that can potentially occur. Our theory suggests that the alignment of desiccation fractures in colloidal drops is selected by the shape of the drop at the point of gelation. We propose that the emergence of three distinct fracture patterns in dried blood drops is a consequence of a non-monotonic drop profile at gelation. We also show that depletion fronts, which separate wet and dry solid, can invade the drop from the contact line and localise the generation of mechanical stress during drying. Finally, the finite element method is used to explore the stress profiles in drops with large contact angles.

In this study we use non-equilibrium thermodynamics to systematically derive a phase-field model of a polyelectrolyte gel coupled to a hydrodynamic model for a salt solution surrounding the gel. The governing equations for the gel account for the free energy of the internal interfaces which form upon phase separation, the nonlinear elasticity of the polyelectrolyte network, and multi-component diffusive transport following a Stefan--Maxwell approach. The time-dependent model describes the evolution of the gel across multiple time and spatial scales and so is able to capture the large-scale solvent flux and the emergence of long-time pattern formation in the system. We explore the model for the case of a constrained gel undergoing uni-axial deformations. Numerical simulations show that rapid changes in the gel volume occur once the volume phase transition sets in, as well as the triggering of spinodal decomposition that leads to strong inhomogeneities in the lateral stresses, potentially leading to experimentally visible patterns.

We analyse the dynamics of different routes to collapse of a constrained polyelectrolyte gel in contact with an ionic bath. The evolution of the gel is described by a model that incorporates non-linear elasticity, Stefan-Maxwell diffusion and interfacial gradient free energy to account for phase separation of the gel. A bifurcation analysis of the homogeneous equilibrium states reveals three solution branches at low ion concentrations in the bath, giving way to only one above a critical ion concentration. We present numerical solutions that capture both the spatial heterogeneity and the multiple time-scales involved in the process of collapse. These solutions are complemented by two analytical studies. Firstly, a phase-plane analysis that reveals the existence of a depletion front for the transition from the highly swollen to the new collapsed equilibrium state. This depletion front is initiated after the fast ionic diffusion has set the initial condition for this time regime. Secondly, we perform a linear stability analysis about the homogeneous states that show that for a range of ion concentrations in the bath, spinodal decomposition of the swollen state gives rise to localized solvent-rich(poor) and, due to the electro-neutrality condition, ion-poor(rich) phases that coarsen on the route to collapse. This dynamics of a collapsing polyelectrolyte gel has not been described before.

In this work we have investigated the evolutionary dynamics of a generalist pathogen, e.g. a virus population, that evolves towards specialisation in an environment with multiple host types. We have particularly explored under which conditions generalist viral strains may rise in frequency and coexist with specialist strains or even dominate the population. By means of a nonlinear mathematical model and bifurcation analysis, we have determined the theoretical conditions for stability of nine identified equilibria and provided biological interpretation in terms of the infection rates for the viral specialist and generalist strains. By means of a stability diagram we identified stable fixed points and stable periodic orbits, as well as regions of bistability. For arbitrary biologically feasible initial population sizes, the probability of evolving towards stable solutions is obtained for each point of the analyzed parameter space. This probability map shows combinations of infection rates of the generalist and specialist strains that might lead to equal chances for each type becoming the dominant strategy. Furthermore, we have identified infection rates for which the model predicts the onset of chaotic dynamics. Several degenerate Bogdanov-Takens and zero-Hopf bifurcations are detected along with generalized Hopf and zero-Hopf bifurcations. This manuscript provides additional insights into the dynamical complexity of host-pathogen evolution towards different infection strategies.

We study two thermo-electrochemical models for lithium-ion batteries. The first is based on volume averaging the electrode microstructure whereas the second is based on the pseudo-two-dimensional (P2D) approach which treats the electrode as a collection of spherical particles. A scaling analysis is used to reduce the volume-averaged model and show that the electrochemical reactions are the dominant source of heat. Matched asymptotic expansions are used to compute solutions of the volume-averaged model for the cases of constant applied current, oscillating applied current, and constant cell potential. The asymptotic and numerical solutions of the volume-averaged model are in remarkable agreement with numerical solutions of the thermal P2D model for (dis)charge rates up to 2C, and reasonable agreement is found at 4C. Homogenisation is then used to derive a thermal model for a battery consisting of several connected lithium-ion cells. Despite accounting for the Arrhenius dependence of the reaction coefficients, we show that thermal runaway does not occur in the model. Instead, the cell potential is simply pushed closer to the open-circuit potential. We also show that in many cases, the homogenised battery model can be solved analytically, making it ideal for use in on-board thermal management systems.

We investigate the one-dimensional growth of a solid into a liquid bath, starting from a small crystal, using the Guyer-Krumhansl and Maxwell-Cattaneo models of heat conduction. By breaking the solidification process into the relevant time regimes we are able to reduce the problem to a system of two coupled ordinary differential equations describing the evolution of the solid-liquid interface and the heat flux. The reduced formulation is in good agreement with numerical simulations. In the case of silicon, differences between classical and non-classical solidification kinetics are relatively small, but larger deviations can be observed in the evolution in time of the heat flux through the growing solid. From this study we conclude that the heat flux provides more information about the presence of non-classical modes of heat transport during phase-change processes., Comment: 29 pages, 6 figures, 2 tables + Supplementary Material

In this paper we formulate a Stefan problem appropriate when the thermophysical properties are distinct in each phase and the phase-change temperature is size or velocity dependent. Thermophysical properties invariably take different values in different material phases but this is often ignored for mathematical simplicity. Size and velocity dependent phase change temperatures are often found at very short length scales, such as nanoparticle melting or dendrite formation; velocity dependence occurs in the solidification of supercooled melts. To illustrate the method we show how the governing equations may be applied to a standard one-dimensional problem and also the melting of a spherically symmetric nanoparticle. Errors which have propagated through the literature are highlighted. By writing the system in non-dimensional form we are able to study the large Stefan number formulation and an energy-conserving one-phase reduction. The results from the various simplifications and assumptions are compared with those from a finite difference numerical scheme. Finally, we briefly discuss the failure of Fourier's law at very small length and time-scales and provide an alternative formulation which takes into account the finite time of travel of heat carriers (phonons) and the mean free distance between collisions., Comment: 39 pages, 5 figures

The role of thermal relaxation in nanoparticle melting is studied using a mathematical model based on the Maxwell--Cattaneo equation for heat conduction. The model is formulated in terms of a two-phase Stefan problem. We consider the cases of the temperature profile being continuous or having a jump across the solid-liquid interface. The jump conditions are derived from the sharp-interface limit of a phase-field model that accounts for variations in the thermal properties between the solid and liquid. The Stefan problem is solved using asymptotic and numerical methods. The analysis reveals that the Fourier-based solution can be recovered from the classical limit of zero relaxation time when either boundary condition is used. However, only the jump condition avoids the onset of unphysical `supersonic' melting, where the speed of the melt front exceeds the finite speed of heat propagation. These results conclusively demonstrate that the jump condition, not the continuity condition, is the most suitable for use in models of phase change based on the Maxwell--Cattaneo equation. Numerical investigations show that thermal relaxation can increase the time required to melt a nanoparticle by more than a factor of ten. Thus, thermal relaxation is an important process to include in models of nanoparticle melting and is expected to be relevant in other rapid phase-change processes.

A mathematical model is presented for thermal transport in nanowires with rectangular cross sections. Expressions for the effective thermal conductivity of the nanowire across a range of temperatures and cross-sectional aspect ratios are obtained by solving the Guyer--Krumhansl hydrodynamic equation for the thermal flux with a slip boundary condition. Our results show that square nanowires transport thermal energy more efficiently than rectangular nanowires due to optimal separation between the boundaries. However, circular nanowires are found to be even more efficient than square nanowires due to the lack of corners in the cross section, which locally reduce the thermal flux and inhibit the conduction of heat. By using a temperature-dependent slip coefficient, we show that the model is able to accurately capture experimental data of the effective thermal conductivity obtained from Si nanowires, demonstrating that phonon hydrodynamics is a powerful framework that can be applied to nanosystems even at room temperature.

We present a porous electrode model for lithium-ion batteries using Butler--Volmer reaction kinetics. We model lithium concentration in both the solid and fluid phase along with solid and liquid electric potential. Through asymptotic reduction, we show that the electric potentials are spatially homogeneous which decouples the problem into a series of time-dependent problems. These problems can be solved on three distinguished time scales, an early time scale where capacitance effects in the electrode dominate, a mid-range time scale where a spatial concentration gradient forms in the electrolyte, and a long-time scale where each of the electrodes saturate and deplete with lithium respectively. The solid-phase concentration profiles are linear functions of time and the electrolyte potential is everywhere zero, which allows the model to be reduced to a system of two uncoupled ordinary differential equations. Analytic and numerical results are compared with full numerical simulations and experimental discharge curves demonstrating excellent agreement., Comment: Accepted in SIAM Journal on Applied Mathematics

Nanoscale solidification is becoming increasingly relevant in applications involving ultra-fast freezing processes and nanotechnology. However, thermal transport on the nanoscale is driven by infrequent collisions between thermal energy carriers known as phonons and is not well described by Fourier's law. In this paper, the role of non-Fourier heat conduction in nanoscale solidification is studied by coupling the Stefan condition to the Guyer--Krumhansl (GK) equation, which is an extension of Fourier's law, valid on the nanoscale, that includes memory and non-local effects. A systematic asymptotic analysis reveals that the solidification process can be decomposed into multiple time regimes, each characterised by a non-classical mode of thermal transport and unique solidification kinetics. For sufficiently large times, Fourier's law is recovered. The model is able to capture the change in the effective thermal conductivity of the solid during its growth, consistent with experimental observations. The results from this study provide key quantitative insights that can be used to control nanoscale solidification processes.

Communicating transactions is a form of distributed, non-isolated transactions which provides a simple construct for building concurrent systems. In this paper we develop a logical framework to express properties of the observable behaviour of such systems. This comprises three nominal modal logics which share standard communication modalities but have distinct past and future modalities involving transactional commits. All three logics have the same distinguishing power over systems because their associated weak bisimulations coincide with contextual equivalence. Furthermore, they are equally expressive because there are semantics-preserving translations between their formulae. Using the logics we can clearly exhibit subtle example inequivalences. This work presents the first property logics for non-isolated transactions., Comment: Published in Information and Computation. This is author's final version

In the standard testing theory of DeNicola-Hennessy one process is considered to be a refinement of another if every test guaranteed by the former is also guaranteed by the latter. In the domain of web services this has been recast, with processes viewed as servers and tests as clients. In this way the standard refinement preorder between servers is determined by their ability to satisfy clients. But in this setting there is also a natural refinement preorder between clients, determined by their ability to be satisfied by servers. In more general settings where there is no distinction between clients and servers, but all processes are peers, there is a further refinement preorder based on the mutual satisfaction of peers. We give a uniform account of these three preorders. In particular we give two characterisations. The first is behavioural, in terms of traces and ready sets. The second, for finite processes, is equational.

The superheating that usually occurs when a solid is melted by volumetric heating can produce irregular solid-liquid interfaces. Such interfaces can be visualised in ice, where they are sometimes known as Tyndall stars. This paper describes some of the experimental observations of Tyndall stars and a mathematical model for the early stages of their evolution. The modelling is complicated by the strong crystalline anisotropy, which results in an anisotropic kinetic undercooling at the interface; it leads to an interesting class of free boundary problems that treat the melt region as infinitesimally thin.

We present a timed process calculus for modelling wireless networks in which individual stations broadcast and receive messages; moreover the broadcasts are subject to collisions. Based on a reduction semantics for the calculus we define a contextual equivalence to compare the external behaviour of such wireless networks. Further, we construct an extensional LTS (labelled transition system) which models the activities of stations that can be directly observed by the external environment. Standard bisimulations in this LTS provide a sound proof method for proving systems contextually equivalence. We illustrate the usefulness of the proof methodology by a series of examples. Finally we show that this proof method is also complete, for a large class of systems.

We define a pi-calculus variant with a costed semantics where channels are treated as resources that must explicitly be allocated before they are used and can be deallocated when no longer required. We use a substructural type system tracking permission transfer to construct coinductive proof techniques for comparing behaviour and resource usage efficiency of concurrent processes. We establish full abstraction results between our coinductive definitions and a contextual behavioural preorder describing a notion of process efficiency w.r.t. its management of resources. We also justify these definitions and respective proof techniques through numerous examples and a case study comparing two concurrent implementations of an extensible buffer., Comment: 51 pages, 7 figures

In this paper the evolution of a binary mixture in a thin-film geometry with a wall at the top and bottom is considered. By bringing the mixture into its miscibility gap so that no spinodal decomposition occurs in the bulk, a slight energetic bias of the walls towards each one of the constituents ensures the nucleation of thin boundary layers that grow until the constituents have moved into one of the two layers. These layers are separated by an interfacial region where the composition changes rapidly. Conditions that ensure the separation into two layers with a thin interfacial region are investigated based on a phase-field model. Using matched asymptotic expansions a corresponding sharp-interface problem for the location of the interface is established. It is then argued that this newly created two-layer system is not at its energetic minimum but destabilizes into a controlled self-replicating pattern of trapezoidal vertical stripes by minimizing the interfacial energy between the phases while conserving their area. A quantitative analysis of this mechanism is carried out via a thin-film model for the free interfaces, which is derived asymptotically from the sharp-interface model., Comment: Submitted 23/12/2013

Session types are used to describe and structure interactions between independent processes in distributed systems. Higher-order types are needed in order to properly structure delegation of responsibility between processes. In this paper we show that higher-order web-service contracts can be used to provide a fully-abstract model of recursive higher-order session types. The model is set-theoretic, in the sense that the meaning of a contract is given in terms of the set of contracts with which it complies. The proof of full-abstraction depends on a novel notion of the complement of a contract. This in turn gives rise to an alternative to the type duality commonly used in systems for type-checking session types. We believe that the notion of complement captures more faithfully the behavioural intuition underlying type duality., Comment: Added definitions of m-closed terms, of 'dual', and a discussion to show the problems of the complement function

We propose a process calculus to model high level wireless systems, where the topology of a network is described by a digraph. The calculus enjoys features which are proper of wireless networks, namely broadcast communication and probabilistic behaviour. We first focus on the problem of composing wireless networks, then we present a compositional theory based on a probabilistic generalisation of the well known may-testing and must-testing pre- orders. Also, we define an extensional semantics for our calculus, which will be used to define both simulation and deadlock simulation preorders for wireless networks. We prove that our simulation preorder is sound with respect to the may-testing preorder; similarly, the deadlock simulation pre- order is sound with respect to the must-testing preorder, for a large class of networks. We also provide a counterexample showing that completeness of the simulation preorder, with respect to the may testing one, does not hold. We conclude the paper with an application of our theory to probabilistic routing protocols.

We introduce a notion of real-valued reward testing for probabilistic processes by extending the traditional nonnegative-reward testing with negative rewards. In this richer testing framework, the may and must preorders turn out to be inverses. We show that for convergent processes with finitely many states and transitions, but not in the presence of divergence, the real-reward must-testing preorder coincides with the nonnegative-reward must-testing preorder. To prove this coincidence we characterise the usual resolution-based testing in terms of the weak transitions of processes, without having to involve policies, adversaries, schedulers, resolutions, or similar structures that are external to the process under investigation. This requires establishing the continuity of our function for calculating testing outcomes., Comment: In Proceedings QAPL 2011, arXiv:1107.0746

Process behaviour is often defined either in terms of the tests they satisfy, or in terms of the logical properties they enjoy. Here we compare these two approaches, using extensional testing in the style of DeNicola, Hennessy, and a recursive version of the property logic HML. We first characterise subsets of this property logic which can be captured by tests. Then we show that those subsets of the property logic capture precisely the power of tests., Comment: In Proceedings EXPRESS'10, arXiv:1011.6012

We develop a version of the pi-calculus, picost, where channels are interpreted as resources which have costs associated with them. Code runs under the financial responsibility of owners; they must pay to use resources, but may profit by providing them. We provide a proof methodology for processes described in picost based on bisimulations. The underlying behavioural theory is justified via a contextual characterisation. We also demonstrate its usefulness via examples.

We view channels as the main form of resources in a message-passing programming paradigm. These channels need to be carefully managed in settings where resources are scarce. To study this problem, we extend the pi-calculus with primitives for channel allocation and deallocation and allow channels to be reused to communicate values of different types. Inevitably, the added expressiveness increases the possibilities for runtime errors. We define a substructural type system which combines uniqueness typing and affine typing to reject these ill-behaved programs.

In 1992 Wang & Larsen extended the may- and must preorders of De Nicola and Hennessy to processes featuring probabilistic as well as nondeterministic choice. They concluded with two problems that have remained open throughout the years, namely to find complete axiomatisations and alternative characterisations for these preorders. This paper solves both problems for finite processes with silent moves. It characterises the may preorder in terms of simulation, and the must preorder in terms of failure simulation. It also gives a characterisation of both preorders using a modal logic. Finally it axiomatises both preorders over a probabilistic version of CSP., Comment: 33 pages

We propose a simple global computing framework, whose main concern is code migration. Systems are structured in sites, and each site is divided into two parts: a computing body, and a membrane, which regulates the interactions between the computing body and the external environment. More precisely, membranes are filters which control access to the associated site, and they also rely on the well-established notion of trust between sites. We develop a basic theory to express and enforce security policies via membranes. Initially, these only control the actions incoming agents intend to perform locally. We then adapt the basic theory to encompass more sophisticated policies, where the number of actions an agent wants to perform, and also their order, are considered., Comment: 23 pages; to appear in Logical Methods in Computer Science