An adaptive pseudo-spectral method for reaction diffusion problems

The spectral interpolation error was considered for both the Chebyshev pseudo-spectral and Galerkin approximations. A family of functionals I sub r (u), with the property that the maximum norm of the error is bounded by I sub r (u)/J sub r, where r is an integer and J is the degree of the polynomial approximation, was developed. These functionals are used in the adaptive procedure whereby the problem is dynamically transformed to minimize I sub r (u). The number of collocation points is then chosen to maintain a prescribed error bound. The method is illustrated by various examples from combustion problems in one and two dimensions

On the Birth of Isolas

Isolas are isolated, closed curves of solution branches of nonlinear problems. They have been observed to occur in the buckling of elastic shells, the equilibrium states of chemical reactors and other problems. In this paper we present a theory to describe analytically the structure of a class of isolas. Specifically, we consider isolas that shrink to a point as a parameter τ of the problem, approaches a critical value τ_0. The point is referred to as an isola center. Equations that characterize the isola centers are given. Then solutions are constructed in a neighborhood of the isola centers by perturbation expansions in a small parameter ε that is proportional to (τ-τo), with a appropriately determined. The theory is applied to a chemical reactor problem

Mean Field Effects for Counterpropagating Traveling Wave Solutions of Reaction-Diffusion Systems

In many problems, e.g., in combustion or solidification, one observes traveling waves that propagate with constant velocity and shape in the x direction, say, are independent of y and z and describe transitions between two equilibrium states, e.g., the burned and the unburned reactants. As parameters of the system are varied, these traveling waves can become unstable and give rise to waves having additional structure, such as traveling waves in the y and z directions, which can themselves be subject to instabilities as parameters are further varied. To investigate this scenario we consider a system of reaction-diffusion equations with a traveling wave solution as a basic state. We determine solutions bifurcating from the basic state that describe counterpropagating traveling waves in directions orthogonal to the direction of propagation of the basic state and determine their stability. Specifically, we derive long wave modulation equations for the amplitudes of the counterpropagating traveling waves that are coupled to an equation for a mean field, generated by the translation of the basic state in the direction of its propagation. The modulation equations are then employed to determine stability boundaries to long wave perturbations for both unidirectional and counterpropagating traveling waves. The stability analysis is delicate because the results depend on the order in which transverse and longitudinal perturbation wavenumbers are taken to zero. For the unidirectional wave we demonstrate that it is sufficient to consider the cases of (i) purely transverse perturbations, (ii) purely longitudinal perturbations, and (iii) longitudinal perturbations with a small transverse component. These yield Eckhaus type, zigzag type, and skew type instabilities, respectively. The latter arise as a specific result of interaction with the mean field. We also consider the degenerate case of very small group velocity, as well as other degenerate cases, which yield several additional instability boundaries. The stability analysis is then extended to the case of counterpropagating traveling waves

Stochastic Chemical Reactions in Micro-domains

Traditional chemical kinetics may be inappropriate to describe chemical reactions in micro-domains involving only a small number of substrate and reactant molecules. Starting with the stochastic dynamics of the molecules, we derive a master-diffusion equation for the joint probability density of a mobile reactant and the number of bound substrate in a confined domain. We use the equation to calculate the fluctuations in the number of bound substrate molecules as a function of initial reactant distribution. A second model is presented based on a Markov description of the binding and unbinding and on the mean first passage time of a molecule to a small portion of the boundary. These models can be used for the description of noise due to gating of ionic channels by random binding and unbinding of ligands in biological sensor cells, such as olfactory cilia, photo-receptors, hair cells in the cochlea.Comment: 33 pages, Journal Chemical Physic

Reactive-diffuse System with Arrhenius Kinetics: Peculiarities of the Spherical Goemetry

The steady reactive-diffusive problem for a non isothermal permeable pellet with first-order Arrhenius kinetics is studied. In the large activation-energy limit, asymptotic solutions are derived for the spherical geometry. The solutions exhibit multiplicity and it is shown that a suitable choice of parameters can lead to an arbitrarily large number of solutions, thereby confirming a conjecture based upon past computational experiments. Explicit analytical expressions are given for the multiplicity bounds (ignition and extinction limits). The asymptotic results compare very well with those obtained numerically, even for moderate values of the activation energy

Thermodynamic Limit Of The Ginzburg-Landau Equations

We investigate the existence of a global semiflow for the complex Ginzburg-Landau equation on the space of bounded functions in unbounded domain. This semiflow is proven to exist in dimension 1 and 2 for any parameter values of the standard cubic Ginzburg-Landau equation. In dimension 3 we need some restrictions on the parameters but cover nevertheless some part of the Benjamin-Feijer unstable domain.Comment: uuencoded dvi file (email: [email protected]

Spatial and spatio-temporal patterns in a cell-haptotaxis model

We investigate a cell-haptotaxis model for the generation of spatial and spatio-temporal patterns in one dimension. We analyse the steady state problem for specific boundary conditions and show the existence of spatially hetero-geneous steady states. A linear analysis shows that stability is lost through a Hopf bifurcation. We carry out a nonlinear multi-time scale perturbation procedure to study the evolution of the resulting spatio-temporal patterns. We also analyse the model in a parameter domain wherein it exhibits a singular dispersion relation

Thermodynamic Limit of a Nonequilibrium Steady-State: Maxwell-Type Construction for a Bistable Biochemical System

We show that the thermodynamic limit of a bistable phosphorylation-dephosphorylation cycle has a selection rule for the "more stable" macroscopic steady state. The analysis is akin to the Maxwell construction. Based on the chemical master equation approach, it is shown that, except at a critical point, bistability disappears in the stochastic model when fluctuation is sufficiently low but unneglectable. Onsager's Gaussian fluctuation theory applies to the unique macroscopic steady state. With initial state in the basin of attraction of the "less stable" steady state, the deterministic dynamics obtained by the Law of Mass Action is a metastable phenomenon. Stability and robustness in cell biology are stochastic concepts.Comment: 12 pages, 2 figure

The Escape Problem for Irreversible Systems

The problem of noise-induced escape from a metastable state arises in physics, chemistry, biology, systems engineering, and other areas. The problem is well understood when the underlying dynamics of the system obey detailed balance. When this assumption fails many of the results of classical transition-rate theory no longer apply, and no general method exists for computing the weak-noise asymptotics of fundamental quantities such as the mean escape time. In this paper we present a general technique for analysing the weak-noise limit of a wide range of stochastically perturbed continuous-time nonlinear dynamical systems. We simplify the original problem, which involves solving a partial differential equation, into one in which only ordinary differential equations need be solved. This allows us to resolve some old issues for the case when detailed balance holds. When it does not hold, we show how the formula for the mean escape time asymptotics depends on the dynamics of the system along the most probable escape path. We also present new results on short-time behavior and discuss the possibility of focusing along the escape path.Comment: 24 pages, APS revtex macros (version 2.1) now available from PBB via `get oldrevtex.sty

On the validity of the geometrical theory of diffraction by convex cylinders

In this paper we consider the scattering of a wave from an infinite line source by an infinitely long cylinder C. The line source is parallel to the axis of C , and the cross section C of this cylinder is smooth, closed and convex. C is formed by joining a pair of smooth convex arcs to a circle C 0 , one on the illuminated side, and one on the dark side, so that C is circular near the points of diffraction. By a rigorous argument we establish the asymptotic behavior of the field at high frequencies, in a certain portion of the shadow S that is determined by the geometry of C in S. The leading term of our asymptotic expansion is the field predicted by the geometrical theory of diffraction.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/46181/1/205_2004_Article_BF00248157.pd