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

Dynamics of polydisperse irreversible adsorption: a pharmacological example

Many drug delivery systems suffer from undesirable interactions with the host immune system. It has been experimentally established that covalent attachment (irreversible adsorption) of suitable macromolecules to the surface of the drug carrier can reduce such undesirable interactions. A fundamental understanding of the adsorption process is still lacking. In this paper, the classical random irreversible adsorption model is generalized to capture certain essential processes involved in pharmacological applications, allowing for macromolecules of different sizes, partial overlapping of the tails of macromolecules, and the influence of reactions with the solvent on the adsorption process. Working in one dimension, an integro-differential evolution equation for the adsorption process is derived, and the asymptotic behavior of the surface area covered and the number of molecules attached to the surface are studied. Finally, equation-free dynamic renormalization tools are applied to study the asymptotically self-similar behavior of the adsorption statistics

On the Lawrence–Doniach and Anisotropic Ginzburg–Landau Models for Layered Superconductors

The authors consider two models, the Lawrence-Doniach and the anisotropic Ginzburg-Landau models for layered superconductors such as the recently discovered high-temperature superconductors. A mathematical description of both models is given and existence results for their solution are derived. The authors then relate the two models in the sense that they show that as the layer spacing tends to zero, the Lawrence-Doniach model reduces to the anisotropic Ginzburg- Landau model. Finally, simplified versions of the models are derived that can be used to accurately simulate high-temperature superconductors

Modelling multiscale aspects of colorectal cancer

Colorectal cancer (CRC) is responsible for nearly half a million deaths annually world-wide [11]. We present a series of mathematical models describing the dynamics of the intestinal epithelium and the kinetics of the molecular pathway most commonly mutated in CRC, the Wnt signalling network. We also discuss how we are coupling such models to build a multiscale model of normal and aberrant guts. This will enable us to combine disparate experimental and clinical data, to investigate interactions between phenomena taking place at different levels of organisation and, eventually, to test the efficacy of new drugs on the system as a whole

Development of a Twin-Spool Turbofan Engine Simulation Using the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS)

The Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS) is a tool that has been developed to allow a user to build custom models of systems governed by thermodynamic principles using a template to model each basic process. Validation of this tool in an engine model application was performed through reconstruction of the Commercial Modular Aero-Propulsion System Simulation (C-MAPSS) (v2) using the building blocks from the T-MATS (v1) library. In order to match the two engine models, it was necessary to address differences in several assumptions made in the two modeling approaches. After these modifications were made, validation of the engine model continued by integrating both a steady-state and dynamic iterative solver with the engine plant and comparing results from steady-state and transient simulation of the T-MATS and C-MAPSS models. The results show that the T-MATS engine model was accurate within 3% of the C-MAPSS model, with inaccuracy attributed to the increased dimension of the iterative solver solution space required by the engine model constructed using the T-MATS library. This demonstrates that, given an understanding of the modeling assumptions made in T-MATS and a baseline model, the T-MATS tool provides a viable option for constructing a computational model of a twin-spool turbofan engine that may be used in simulation studies

The effect of weak inertia in rotating high-aspect-ratio vessel bioreactors

One method to grow artificial body tissue is to place a porous scaffold seeded with cells, known as a tissue construct, into a rotating bioreactor filled with a nutrient-rich fluid. The flow within the bioreactor is affected by the movement of the construct relative to the bioreactor which, in turn, is affected by the hydrodynamical and gravitational forces the construct experiences. The construct motion is thus coupled to the flow within the bioreactor. Over the timescale of a few hours, the construct appears to move in a periodic orbit but, over tens of hours, the construct drifts from periodicity. In the biological literature, this effect is often attributed to the change in density of the construct that occurs via tissue growth. In this paper, we show that weak inertia can cause the construct to drift from its periodic orbit over the same timescale as tissue growth. We consider the coupled flow and construct motion problem within a rotating high-aspect- ratio vessel bioreactor. Using an asymptotic analysis, we investigate the case where the Reynolds number is large but the geometry of the bioreactor yields a small reduced Reynolds number, resulting in a weak inertial effect. In particular, to accurately couple the bioreactor and porous flow regions, we extend the nested boundary layer analysis of Dalwadi et al. (J. Fluid Mech. vol. 798, pp. 88–139, 2016) to include moving walls and the thin region between the porous construct and the bioreactor wall. This allows us to derive a closed system of nonlinear ordinary differential equations for the construct trajectory, from which we show that neglecting inertia results in periodic orbits; we solve the inertia-free problem analytically, calculating the periodic orbits in terms of the system parameters. Using a multiple-scale analysis, we then systematically derive a simpler system of nonlinear ordinary differential equations that describe the long-time drift of the construct due to the effect of weak inertia. We investigate the bifurcations of the construct trajectory behaviour, and the limit cycles that appear when the construct is less dense than the surrounding fluid and the rotation rate is large enough. Thus, we are able to predict when the tissue construct will drift towards a stable limit cycle within the bioreactor and when it will drift out until it hits the bioreactor edg

LADEE Propulsion System Cold Flow Test

Lunar Atmosphere and Dust Environment Explorer (LADEE) is a NASA mission that will orbit the Moon. Its main objective is to characterize the atmosphere and lunar dust environment. The spacecraft development is being led by NASA Ames Research Center and scheduled for launch in 2013. The LADEE spacecraft will be operated with a bi-propellant hypergolic propulsion system using MMH and NTO as the fuel and oxidizer, respectively. The propulsion system utilizes flight-proven hardware on major components. The propulsion layout is composed of one 100-lbf main thruster and four 5-lbf RCS thrusters. The propellants are stored in four tanks (two parallel-connected tanks per propellant component). The propellants will be pressurized by regulated helium. A simulated propulsion system has been built for conducting cold flow test series to characterize the transient fluid flow of the propulsion system feed lines and to verify the critical operation modes, such as system priming, waterhammer, and crucial mission duty cycles. Propellant drainage differential between propellant tanks will also be assessed. Since the oxidizer feed line system has a higher flow demand than the fuel system does, the cold flow test focuses on the oxidizer system. The objective of the cold flow test is to simulate the LADEE propulsion fluid flow operation through water cold flow test and to obtain data for anchoring analytical models. The models will be used to predict the transient and steady state flow behaviors in the actual flight operations. The test activities, including the simulated propulsion test article, cold flow test, and analytical modeling, are being performed at NASA Marshall Space Flight Center. At the time of the abstract submission, the test article checkout is being performed. The test series will be completed by November, 201

Stochastic modelling of reaction-diffusion processes: algorithms for bimolecular reactions

Several stochastic simulation algorithms (SSAs) have been recently proposed for modelling reaction-diffusion processes in cellular and molecular biology. In this paper, two commonly used SSAs are studied. The first SSA is an on-lattice model described by the reaction-diffusion master equation. The second SSA is an off-lattice model based on the simulation of Brownian motion of individual molecules and their reactive collisions. In both cases, it is shown that the commonly used implementation of bimolecular reactions (i.e. the reactions of the form A + B -> C, or A + A -> C) might lead to incorrect results. Improvements of both SSAs are suggested which overcome the difficulties highlighted. In particular, a formula is presented for the smallest possible compartment size (lattice spacing) which can be correctly implemented in the first model. This implementation uses a new formula for the rate of bimolecular reactions per compartment (lattice site).Comment: 33 pages, submitted to Physical Biolog

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.Comment: 24 pages, submitted to Physical Biolog