BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. Jonathan Chapman

S. Jonathan Chapman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynmath-phmath.MPBiological PhysicsQuantitative Methodsmath.APphysics.chem-phcond-mat.softmath.CAmath.NAphysics.app-phphysics.comp-phSubcellular Processes

Catalog footprint

What is connected

16works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2020arXiv

Modelling electrode heterogeneity in lithium-ion batteries: unimodal and bimodal particle-size distributions

In mathematical models of lithium-ion batteries, the highly heterogeneous porous electrodes are frequently approximated as comprising spherical particles of uniform size, leading to the commonly-used single-particle model (SPM) when transport in the electrolyte is assumed to be fast. Here electrode heterogeneity is modelled by extending this to a distribution of particle sizes. Unimodal and bimodal particle-size distributions (PSD) are considered. For a unimodal PSD, the effect of the spread of the distribution on the cell dynamics is investigated, and choice of effective particle radius when approximating by an SPM assessed. Asymptotic techniques are used to derive a correction to the SPM valid for narrow, but realistic, PSDs. In addition, it is shown that the heterogeneous internal states of all particles (relevant when modelling degradation, for example) can be efficiently computed after-the-fact. For a bimodal PSD, the results are well approximated by a double-particle model (DPM), with one size representing each mode. Results for lithium iron phosphate with a bimodal PSD show that the DPM captures an experimentally-observed double-plateau in the discharge curve, suggesting it is entirely due to bimodality.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Shear-induced instabilities of flows through submerged vegetation

We consider the instabilities of flows through a submerged canopy, and show how the full governing equations of the fluid-structure interactions can be reduced to a compact framework that captures many key features of vegetative flow. By modelling the canopy as a collection of homogeneous elastic beams, we predict the steady configuration of the plants in response to a unidirectional flow. This treatment couples the beam equations in the canopy to the fluid momentum equations. Our linear stability analysis suggests new insights into the development of instabilities at the surface of the vegetative region. In particular, we show that shear at the top of the canopy is a dominant factor in determining the onset of instabilities known as monami. Based on numerical and asymptotic analysis of the generalised eigenvalue problem, the system is shown to be stable if the canopy is sufficiently sparse or if the plants are sufficiently flexible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the boundary layer structure near a highly permeable porous interface

The method of matched asymptotic expansions is used to study the canonical problem of steady laminar flow through a narrow two-dimensional channel blocked by a tight-fitting finite-length highly permeable porous obstacle. We investigate the behaviour of the local flow close to the interface between the single-phase and porous regions (governed by the incompressible Navier--Stokes and Darcy flow equations, respectively). We solve for the flow in these inner regions in the limits of low and high Reynolds number, facilitating an understanding of the nature of the transition from Poiseuille to plug to Poiseuille flow in each of these limits. Significant analytic progress is made in the high-Reynolds-number limit, and we explore in detail the rich boundary layer structure that occurs. We consider the three-dimensional generalization to unsteady laminar flow through and around a tight-fitting highly permeable cylindrical porous obstacle within a Hele-Shaw cell. For the high-Reynolds-number limit, we give the coupling conditions and interfacial stress in terms of the outer flow variables, allowing information from a nonlinear three-dimensional problem to be obtained by solving a linear two-dimensional problem. Finally, we illustrate the utility of our analysis by considering the specific example of time-dependent forced far-field flow in a Hele-Shaw cell containing a porous cylinder with a circular cross-section. We determine the internal stress within the porous obstacle, which is key for tissue engineering applications, and the interfacial stress on the boundary of the porous obstacle, which has applications to biofilm erosion. In the high-Reynolds-number limit, we demonstrate that the fluid inertia can result in the cylinder experiencing a time-independent net force, even when the far-field forcing is periodic with zero mean.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reactive Boundary Conditions as Limits of Interaction Potentials for Brownian and Langevin Dynamics

A popular approach to modeling bimolecular reactions between diffusing molecules is through the use of reactive boundary conditions. One common model is the Smoluchowski partial absorption condition, which uses a Robin boundary condition in the separation coordinate between two possible reactants. This boundary condition can be interpreted as an idealization of a reactive interaction potential model, in which a potential barrier must be surmounted before reactions can occur. In this work we show how the reactive boundary condition arises as the limit of an interaction potential encoding a steep barrier within a shrinking region in the particle separation, where molecules react instantly upon reaching the peak of the barrier. The limiting boundary condition is derived by the method of matched asymptotic expansions, and shown to depend critically on the relative rate of increase of the barrier height as the width of the potential is decreased. Limiting boundary conditions for the same interaction potential in both the overdamped Fokker-Planck equation (Brownian Dynamics), and the Kramers equation (Langevin Dynamics) are investigated. It is shown that different scalings are required in the two models to recover reactive boundary conditions that are consistent in the high friction limit where the Kramers equation solution converges to the solution of the Fokker-Planck equation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Curvature suppresses the Rayleigh-Taylor instability

The dynamics of a thin liquid film on the underside of a curved cylindrical substrate is studied. The evolution of the liquid layer is investigated as the film thickness and the radius of curvature of the substrate are varied. A dimensionless parameter (a modified Bond number) that incorporates both geometric parameters, gravity, and surface tension is identified, and allows the observations to be classified according to three different flow regimes: stable films, films with transient growth of perturbations followed by decay, and unstable films. Experiments and theory confirm that, below a critical value of the Bond number, curvature of the substrate suppresses the Rayleigh-Taylor instability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exponential asymptotics and problems with coalescing singularities

Problems in exponential asymptotics are typically characterized by divergence of the associated asymptotic expansion in the form of a factorial divided by a power. In this paper, we demonstrate that in certain classes of problems that involve coalescing singularities, a more general type of exponential-over-power divergence must be used. As a model example, we study the water waves produced by flow past an obstruction such as a surface-piercing ship. In the low speed or low Froude limit, the resultant water waves are exponentially small, and their formation is attributed to the singularities in the geometry of the obstruction. However, in cases where the singularities are closely spaced, the usual asymptotic theory fails. We present both a general asymptotic framework for handling such problems of coalescing singularities, and provide numerical and asymptotic results for particular examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Model reduction for slow-fast stochastic systems with metastable behaviour

The quasi-steady-state approximation (or stochastic averaging principle) is a useful tool in the study of multiscale stochastic systems, giving a practical method by which to reduce the number of degrees of freedom in a model. The method is extended here to slow-fast systems in which the fast variables exhibit metastable behaviour. The key parameter that determines the form of the reduced model is the ratio of the timescale for the switching of the fast variables between metastable states to the timescale for the evolution of the slow variables. The method is illustrated with two examples: one from biochemistry (a fast-species-mediated chemical switch coupled to a slower-varying species), and one from ecology (a predator-prey system). Numerical simulations of each model reduction are compared with those of the full system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analysis of the two-regime method on square meshes

The two-regime method (TRM) has been recently developed for optimizing stochastic reaction-diffusion simulations. It is a multiscale (hybrid) algorithm which uses stochastic reaction-diffusion models with different levels of detail in different parts of the computational domain. The coupling condition on the interface between different modelling regimes of the TRM was previously derived for one-dimensional models. In this paper, the TRM is generalized to higher dimensional reaction-diffusion systems. Coupling Brownian dynamics models with compartment-based models on regular (square) two-dimensional lattices is studied in detail. In this case, the interface between different modelling regimes contain either flat parts or right-angled corners. Both cases are studied in the paper. For flat interfaces, it is shown that the one-dimensional theory can be used along the line perpendicular to the TRM interface. In the direction tangential to the interface, two choices of the TRM parameters are presented. Their applicability depends on the compartment size and the time step used in the molecular-based regime. The two-dimensional generalization of the TRM is also discussed in the case of corners.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New gravity-capillary waves at low speeds. Part 1: Linear geometries

When traditional linearised theory is used to study gravity-capillary waves produced by flow past an obstruction, the geometry of the object is assumed to be small in one or several of its dimensions. In order to preserve the nonlinear nature of the obstruction, asymptotic expansions in the low-Froude or low-Bond number limits can be derived, but here, the solutions invariably predict a waveless surface at every order. This is because the waves are in fact, exponentially small, and thus beyond-all-orders of regular asymptotics; their formation is a consequence of the divergence of the asymptotic series and the associated Stokes Phenomenon. By applying techniques in exponential asymptotics to this problem, we have discovered the existence of new classes of gravity-capillary waves, from which the usual linear solutions form but a special case. In this paper, we present the initial theory for deriving these waves through a study of gravity-capillary flow over a linearised step; this will be done using two approaches: in the first, we derive the surface waves using the the standard method of Fourier transforms; in the second, we derive the same result using exponential asymptotics. Ultimately, these two methods give the same result, but conceptually, they offer different insights into the study of the low-Froude, low-Bond number problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New gravity-capillary waves at low speeds. Part 2: Nonlinear geometries

When traditional linearised theory is used to study gravity-capillary waves produced by flow past an obstruction, the geometry of the object is assumed to be small in one or several of its dimensions. In order to preserve the nonlinear nature of the obstruction, asymptotic expansions in the low-Froude or low-Bond number limits can be derived, but here, the solutions are waveless to every order. This is because the waves are in fact, exponentially small, and thus beyond-all-orders of regular asymptotics; their formation is a consequence of the divergence of the asymptotic series and the associated Stokes Phenomenon. In Part 1, we showed how exponential asymptotics could be used to study the problem when the size of the obstruction is first linearised. In this paper, we extend the analysis to the nonlinear problem, thus allowing the full geometry to be considered at leading order. When applied to the classic problem of flow over a step, our analysis reveals the existence of six classes of gravity-capillary waves, from which two share a connection with the usual linearised solutions first discovered by Lord Rayleigh. The new solutions arise due to the availability of multiple singularities in the geometry, coupled with the interplay of gravitational and cohesive effects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Waveless ships in the low speed limit: Results for multi-cornered hulls

In the low-speed limit, a blunt ship modeled as two-dimensional semi-infinite body with a single corner can never be made waveless. This was the conclusion of the previous part of our work in Trinh et al. (2011), which focused on the Dagan & Tulin (1972) model of ship waves in the low speed limit. In this accompanying paper, we continue our investigations with the study of more general, piecewise-linear, or multi-cornered ships. The low-speed or low-Froude limit, coupled with techniques in exponential asymptotics allows us to derive explicit formulae relating the geometry of the hull to the form of the waves. Configurations with closely spaced corners present a non-trivial extension of the theory, and we present the general methodology for their study. Lastly, numerical computations of the nonlinear ship-wave problem are presented in order to confirm the analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diffusion of multiple species with excluded-volume effects

Stochastic models of diffusion with excluded-volume effects are used to model many biological and physical systems at a discrete level. The average properties of the population may be described by a continuum model based on partial differential equations. In this paper we consider multiple interacting subpopulations/species and study how the inter-species competition emerges at the population level. Each individual is described as a finite-size hard core interacting particle undergoing Brownian motion. The link between the discrete stochastic equations of motion and the continuum model is considered systematically using the method of matched asymptotic expansions. The system for two species leads to a nonlinear cross-diffusion system for each subpopulation, which captures the enhancement of the effective diffusion rate due to excluded-volume interactions between particles of the same species, and the diminishment due to particles of the other species. This model can explain two alternative notions of the diffusion coefficient that are often confounded, namely collective diffusion and self-diffusion. Simulations of the discrete system show good agreement with the analytic results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Excluded-volume effects in the diffusion of hard spheres

Excluded-volume effects can play an important role in determining transport properties in diffusion of particles. Here, the diffusion of finite-sized hard-core interacting particles in two or three dimensions is considered systematically using the method of matched asymptotic expansions. The result is a nonlinear diffusion equation for the one-particle distribution function, with excluded-volume effects enhancing the overall collective diffusion rate. An expression for the effective (collective) diffusion coefficient is obtained. Stochastic simulations of the full particle system are shown to compare well with the solution of this equation for two examples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From Brownian Dynamics to Markov Chain: an Ion Channel Example

A discrete rate theory for general multi-ion channels is presented, in which the continuous dynamics of ion diffusion is reduced to transitions between Markovian discrete states. In an open channel, the ion permeation process involves three types of events: an ion entering the channel, an ion escaping from the channel, or an ion hopping between different energy minima in the channel. The continuous dynamics leads to a hierarchy of Fokker-Planck equations, indexed by channel occupancy. From these the mean escape times and splitting probabilities (denoting from which side an ion has escaped) can be calculated. By equating these with the corresponding expressions from the Markov model the Markovian transition rates can be determined. The theory is illustrated with a two-ion one-well channel. The stationary probability of states is compared with that from both Brownian dynamics simulation and the hierarchical Fokker-Planck equations. The conductivity of the channel is also studied, and the optimal geometry maximizing ion flux is computed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multiscale reaction-diffusion algorithms: PDE-assisted Brownian dynamics

Two algorithms that combine Brownian dynamics (BD) simulations with mean-field partial differential equations (PDEs) are presented. This PDE-assisted Brownian dynamics (PBD) methodology provides exact particle tracking data in parts of the domain, whilst making use of a mean-field reaction-diffusion PDE description elsewhere. The first PBD algorithm couples BD simulations with PDEs by randomly creating new particles close to the interface which partitions the domain and by reincorporating particles into the continuum PDE-description when they cross the interface. The second PBD algorithm introduces an overlap region, where both descriptions exist in parallel. It is shown that to accurately compute variances using the PBD simulation requires the overlap region. Advantages of both PBD approaches are discussed and illustrative numerical examples are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Realistic boundary conditions for stochastic simulations of reaction-diffusion processes

Many cellular and subcellular biological processes can be described in terms of diffusing and chemically reacting species (e.g. enzymes). Such reaction-diffusion processes can be mathematically modelled using either deterministic partial-differential equations or stochastic simulation algorithms. The latter provide a more detailed and precise picture, and several stochastic simulation algorithms have been proposed in recent years. Such models typically give the same description of the reaction-diffusion processes far from the boundary of the simulated domain, but the behaviour close to a reactive boundary (e.g. a membrane with receptors) is unfortunately model-dependent. In this paper, we study four different approaches to stochastic modelling of reaction-diffusion problems and show the correct choice of the boundary condition for each model. The reactive boundary is treated as partially reflective, which means that some molecules hitting the boundary are adsorbed (e.g. bound to the receptor) and some molecules are reflected. The probability that the molecule is adsorbed rather than reflected depends on the reactivity of the boundary (e.g. on the rate constant of the adsorbing chemical reaction and on the number of available receptors), and on the stochastic model used. This dependence is derived for each model.

0 citations0 reviewsNo code linkedOpen access