73 research outputs found
Analysis of the accuracy and convergence of equation-free projection to a slow manifold
In [C.W. Gear, T.J. Kaper, I.G. Kevrekidis, and A. Zagaris, Projecting to a
Slow Manifold: Singularly Perturbed Systems and Legacy Codes, SIAM J. Appl.
Dyn. Syst. 4 (2005) 711-732], we developed a class of iterative algorithms
within the context of equation-free methods to approximate low-dimensional,
attracting, slow manifolds in systems of differential equations with multiple
time scales. For user-specified values of a finite number of the observables,
the m-th member of the class of algorithms (m = 0, 1, ...) finds iteratively an
approximation of the appropriate zero of the (m+1)-st time derivative of the
remaining variables and uses this root to approximate the location of the point
on the slow manifold corresponding to these values of the observables. This
article is the first of two articles in which the accuracy and convergence of
the iterative algorithms are analyzed. Here, we work directly with explicit
fast--slow systems, in which there is an explicit small parameter, epsilon,
measuring the separation of time scales. We show that, for each m = 0, 1, ...,
the fixed point of the iterative algorithm approximates the slow manifold up to
and including terms of O(epsilon^m). Moreover, for each m, we identify
explicitly the conditions under which the m-th iterative algorithm converges to
this fixed point. Finally, we show that when the iteration is unstable (or
converges slowly) it may be stabilized (or its convergence may be accelerated)
by application of the Recursive Projection Method. Alternatively, the
Newton-Krylov Generalized Minimal Residual Method may be used. In the
subsequent article, we will consider the accuracy and convergence of the
iterative algorithms for a broader class of systems-in which there need not be
an explicit small parameter-to which the algorithms also apply
Reduction for Michaelis-Menten-Henri kinetics in the presence of diffusion
The Michaelis-Menten-Henri (MMH) mechanism is one of the paradigm reaction mechanisms in biology and chemistry. In its simplest form, it involves a substrate that reacts (reversibly) with an enzyme, forming a complex which is transformed (irreversibly) into a product and the enzyme. Given these basic kinetics, a dimension reduction has traditionally been achieved in two steps, by using conservation relations to reduce the number of species and by exploiting the inherent fast-slow structure of the resulting equations. In the present article, we investigate how the dynamics change if the species are additionally allowed to diffuse. We study the two extreme regimes of large diffusivities and of small diffusivities, as well as an intermediate regime in which the time scale of diffusion is comparable to that of the fast reaction kinetics. We show that reduction is possible in each of these regimes, with the nature of the reduction being regime dependent. Our analysis relies on the classical method of matched asymptotic expansions to derive approximations for the solutions that are uniformly valid in space and time
Equation-Free Analysis of Macroscopic Behavior in Traffic and Pedestrian Flow
Equation-free methods make possible an analysis of the evolution of a few
coarse-grained or macroscopic quantities for a detailed and realistic model
with a large number of fine-grained or microscopic variables, even though no
equations are explicitly given on the macroscopic level. This will facilitate a
study of how the model behavior depends on parameter values including an
understanding of transitions between different types of qualitative behavior.
These methods are introduced and explained for traffic jam formation and
emergence of oscillatory pedestrian counter flow in a corridor with a narrow
door
A mathematical framework for critical transitions: normal forms, variance and applications
Critical transitions occur in a wide variety of applications including
mathematical biology, climate change, human physiology and economics. Therefore
it is highly desirable to find early-warning signs. We show that it is possible
to classify critical transitions by using bifurcation theory and normal forms
in the singular limit. Based on this elementary classification, we analyze
stochastic fluctuations and calculate scaling laws of the variance of
stochastic sample paths near critical transitions for fast subsystem
bifurcations up to codimension two. The theory is applied to several models:
the Stommel-Cessi box model for the thermohaline circulation from geoscience,
an epidemic-spreading model on an adaptive network, an activator-inhibitor
switch from systems biology, a predator-prey system from ecology and to the
Euler buckling problem from classical mechanics. For the Stommel-Cessi model we
compare different detrending techniques to calculate early-warning signs. In
the epidemics model we show that link densities could be better variables for
prediction than population densities. The activator-inhibitor switch
demonstrates effects in three time-scale systems and points out that excitable
cells and molecular units have information for subthreshold prediction. In the
predator-prey model explosive population growth near a codimension two
bifurcation is investigated and we show that early-warnings from normal forms
can be misleading in this context. In the biomechanical model we demonstrate
that early-warning signs for buckling depend crucially on the control strategy
near the instability which illustrates the effect of multiplicative noise.Comment: minor corrections to previous versio
Equation-Free Analysis of Two-Component System Signalling Model Reveals the Emergence of Co-Existing Phenotypes in the Absence of Multistationarity
Phenotypic differences of genetically identical cells under the same environmental conditions have been attributed to the inherent stochasticity of biochemical processes. Various mechanisms have been suggested, including the existence of alternative steady states in regulatory networks that are reached by means of stochastic fluctuations, long transient excursions from a stable state to an unstable excited state, and the switching on and off of a reaction network according to the availability of a constituent chemical species. Here we analyse a detailed stochastic kinetic model of two-component system signalling in bacteria, and show that alternative phenotypes emerge in the absence of these features. We perform a bifurcation analysis of deterministic reaction rate equations derived from the model, and find that they cannot reproduce the whole range of qualitative responses to external signals demonstrated by direct stochastic simulations. In particular, the mixed mode, where stochastic switching and a graded response are seen simultaneously, is absent. However, probabilistic and equation-free analyses of the stochastic model that calculate stationary states for the mean of an ensemble of stochastic trajectories reveal that slow transcription of either response regulator or histidine kinase leads to the coexistence of an approximate basal solution and a graded response that combine to produce the mixed mode, thus establishing its essential stochastic nature. The same techniques also show that stochasticity results in the observation of an all-or-none bistable response over a much wider range of external signals than would be expected on deterministic grounds. Thus we demonstrate the application of numerical equation-free methods to a detailed biochemical reaction network model, and show that it can provide new insight into the role of stochasticity in the emergence of phenotypic diversity
- …