Corrected mean-field models for spatially-dependent advection-diffusion-reaction phenomena

In the exclusion-process literature, mean-field models are often derived by assuming that the occupancy status of lattice sites is independent. Although this assumption is questionable, it is the foundation of many mean-field models. In this work we develop methods to relax the independence assumption for a range of discrete exclusion process-based mechanisms motivated by applications from cell biology. Previous investigations that focussed on relaxing the independence assumption have been limited to studying initially-uniform populations and ignored any spatial variations. By ignoring spatial variations these previous studies were greatly simplified due to translational invariance of the lattice. These previous corrected mean-field models could not be applied to many important problems in cell biology such as invasion waves of cells that are characterised by moving fronts. Here we propose generalised methods that relax the independence assumption for spatially inhomogeneous problems, leading to corrected mean-field descriptions of a range of exclusion process-based models that incorporate (i) unbiased motility, (ii) biased motility, and (iii) unbiased motility with agent birth and death processes. The corrected mean-field models derived here are applicable to spatially variable processes including invasion wave type problems. We show that there can be large deviations between simulation data and traditional mean-field models based on invoking the independence assumption. Furthermore, we show that the corrected mean-field models give an improved match to the simulation data in all cases considered.

Solution methods for advection-diffusion-reaction equations on growing domains and subdomains, with application to modelling skin substitutes

Problems involving the solution of advection-diffusion-reaction equations on domains and subdomains whose growth affects and is affected by these equations, commonly arise in developmental biology. Here, a mathematical framework for these situations, together with methods for obtaining spatio-temporal solutions and steady states of models built from this framework, is presented. The framework and methods are applied to a recently published model of epidermal skin substitutes. Despite the use of Eulerian schemes, excellent agreement is obtained between the numerical spatio-temporal, numerical steady state, and analytical solutions of the model.

Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems on Anisotropically Refined Meshes

In this paper we consider the a posteriori and a priori error analysis of discontinuous Galerkin interior penalty methods for second-order partial differential equations with nonnegative characteristic form on anisotropically refined computational meshes. In particular, we discuss the question of error estimation for linear target functionals, such as the outflow flux and the local average of the solution. Based on our a posteriori error bound we design and implement the corresponding adaptive algorithm to ensure reliable and efficient control of the error in the prescribed functional to within a given tolerance. This involves exploiting both local isotropic and anisotropic mesh refinement. The theoretical results are illustrated by a series of numerical experiments.

A Pre-processing Moving Mesh Method for Discontinuous Galerkin Approximations of Advection-Diffusion-Reaction Problems

We propose a pre-processing mesh re-distribution algorithm based upon harmonic maps employed in conjunction with discontinuous Galerkin approximations of advection-diffusion-reaction problems. Extensive two-dimensional numerical experiments with different choices of monitor functions, including monitor functions derived from goal-oriented a posteriori error indicators are presented. The examples presented clearly demonstrate the capabilities and the benefits of combining our pre-processing mesh movement algorithm with both uniform, as well as, adaptive isotropic and anisotropic mesh refinement.

Asymptotic stability of a coupled Advection-Diffusion-Reaction system arising in bioreactor processes.

In this work, we perform an asymptotic analysis of a coupled system of two Advection-Diffusion-Reaction equations with Danckwerts boundary conditions, which models the interaction between a microbial population (e.g., bacterias), called biomass, and a diluted organic contaminant (e.g., nitrates), called substrate, in a continuous flow bioreactor. This system exhibits, under suitable conditions, two stable equilibrium states: one steady state in which the biomass becomes extinct and no reaction is produced, called washout, and another steady state, which corresponds to the partial elimination of the substrate. We use the method of linearization to give sufficient conditions for the asymptotic stability of the two stable equilibrium configurations. Finally, we compare our asymptotic analysis with the usual asymptotic analysis associated to the continuous bioreactor when it is modeled with ordinary differential equations.

A subgrid viscosity Lagrance-Galerkin method for convection-diffusion problems

We present and analyze a subgrid viscosity Lagrange-Galerk in method that combines the subgrid eddy viscosity method proposed in W. Layton, A connection between subgrid scale eddy viscosity and mixed methods. Appl. Math. Comp., 133: 14 7-157, 2002, and a conventional Lagrange-Galerkin method in the framework of P1⊕ cubic bubble finite elements. This results in an efficient and easy to implement stabilized method for convection dominated convection diffusion reaction problems. Numerical experiments support the numerical analysis results and show that the new method is more accurate than the conventional Lagrange-Galerkin one.

Critical time scales for advection–diffusion–reaction processes

The concept of local accumulation time (LAT) was introduced by Berezhkovskii and coworkers in 2010–2011 to give a finite measure of the time required for the transient solution of a reaction–diffusion equation to approach the steady–state solution (Biophys J. 99, L59 (2010); Phys Rev E. 83, 051906 (2011)). Such a measure is referred to as a critical time. Here, we show that LAT is, in fact, identical to the concept of mean action time (MAT) that was first introduced by McNabb in 1991 (IMA J Appl Math. 47, 193 (1991)). Although McNabb’s initial argument was motivated by considering the mean particle lifetime (MPLT) for a linear death process, he applied the ideas to study diffusion. We extend the work of these authors by deriving expressions for the MAT for a general one–dimensional linear advection–diffusion–reaction problem. Using a combination of continuum and discrete approaches, we show that MAT and MPLT are equivalent for certain uniform–to-uniform transitions; these results provide a practical interpretation for MAT, by directly linking the stochastic microscopic processes to a meaningful macroscopic timescale. We find that for more general transitions, the equivalence between MAT and MPLT does not hold. Unlike other critical time definitions, we show that it is possible to evaluate the MAT without solving the underlying partial differential equation (pde). This makes MAT a simple and attractive quantity for practical situations. Finally, our work explores the accuracy of certain approximations derived using the MAT, showing that useful approximations for nonlinear kinetic processes can be obtained, again without treating the governing pde directly.

Transient analysis of mammalian cell growth in hollow fibre bioreactor

A mathematical model describing the dynamics of mammalian cell growth in hollow fibre bioreactor operated in closed shell mode is developed. Mammalian cells are assumed to grow as an expanding biofilm in the extra-capillary space surrounding the fibre. Diffusion is assumed to be the dominant process in the radial direction while axial convection dominates in the lumen of the bioreactor. The transient simulation results show that steep gradients in the cell number are possible under the condition of substrate limitation. The precise conditions which result in nonuniform growth of cells along the length of the bioreactor are delineated. The effect of various operating conditions, such as substrate feed rate, length of the bioreactor and diffusivity of substrate in different regions of the bioreactor, on the bioreactor performance are evaluated in terms of time required to attain the steady-state. The rime of growth is introduced as a measure of effectiveness factor for the bioreactor and is found to be dependent on two parameters, a modified Peclet number and a Thiele modulus. Diffusion, reaction and/or convection control regimes are identified based on these two parameters. The model is further extended to include dual substrate growth limitations, and the relative growth limiting characteristics of two substrates are evaluated. (C) 1997 Elsevier Science Ltd.

A robust layer adapted difference method for singularly perturbed two-parameter parabolic problems

A finite difference method for a time-dependent singularly perturbed convection-diffusion-reaction problem involving two small parameters in one space dimension is considered. We use the classical implicit Euler method for time discretization and upwind scheme on the Shishkin-Bakhvalov mesh for spatial discretization. The method is analysed for convergence and is shown to be uniform with respect to both the perturbation parameters. The use of the Shishkin-Bakhvalov mesh gives first-order convergence unlike the Shishkin mesh where convergence is deteriorated due to the presence of a logarithmic factor. Numerical results are presented to validate the theoretical estimates obtained.

Concentration distribution in solution crystal growth: effect of moving interface conditions

The linear diffusion-reaction theory with finite interface kinetics is employed to describe the dissolution and the growth processes. The results show that it is imperative to consider the effect of the moving interfaces on the concentration distribution at the growth interface for some cases. For small aspect ratio and small gravity magnitude, the dissolution and the growth interfaces must be treated as the moving boundaries within an angle range of 0 degrees < gamma < 50 degrees in this work. For large aspect ratio or large gravity magnitude, the effect of the moving interfaces on the concentration distribution at the growth interface can be neglected except for gamma < - 50 degrees.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

Concentration distribution in crystallization from solution under microgravity

Concentration distribution in crystallization from solution under microgravity is numerically studied. A quasi-steady state growth and dissolution in a 2D rectangular enclosure filled with sodium chlorate (NaClO3) aqueous solution, in which one wall is the growth surface of the crystal and the opposite one is the dissolution surface, is considered. The solute transport process at the growth surface is described by the diffusion-reaction theory with finite interface kinetics coefficient. The results show that the concentration at the growth surface is supersaturated and the supersaturation distribution is of non-uniformity, i.e. the supersaturation in a region facing an incoming flow is high. On the other hand, the non-uniformity of supersaturation at the growth surface is closely related to the gravity level even under microgravity, it exponentially increases as the thermal Rayleigh number on behalf of the gravity level rises.