Discrete and Continuous Dynamical Systems - B
June 2009 , Volume 11 , Issue 4
Select all articles
We consider a two-component Reaction-Diffusion system posed on non coincident spatial domains and featuring a reaction term involving an integral kernel. The question of global existence of componentwise nonnegative solutions is assessed. Then we investigate the stabilization of one of the solution components to zero via an internal control distributed on a small subdomain while preserving nonnegativity of both components. Our results apply to predator-prey systems.
A finite-volume scheme for a nonlocal three-component reaction-diffusion system modeling an epidemic disease with susceptible, infected, and recovered, individuals is analyzed. For this SIR model, the existence of solutions to the finite volume scheme and its convergence to a weak solution of the PDE is established. The convergence proof is based on deriving a series of apriori estimates and by using a general $L^p$ compactness criterion. Finally, numerical simulations from the finite volume scheme are given.
Within the context of multiscale computations, equation-free methods have been developed. In this approach, the evolution of a system is simulated on the macroscopic level while only a microscopic model is explicitly available. To this end, a coarse time stepper for the macroscopic variables can be constructed, based on appropriately initialized microscopic simulations. In this paper, we investigate the initialization of the microscopic simulator using the macroscopic variables only (called lifting in the equation-free framework) when the microscopic model is a molecular dynamics (MD) description of a mono-atomic dense fluid. We assume a macroscopic model to exist in terms of the lowest order velocity moments of the particle distribution (density, velocity and temperature). The major difficulty is to design a lifting operator that accurately reconstructs the physically correct state of the fluid (i.e., the higher order moments) at a reasonable computational cost. We construct a lifting operator, as well as a restriction operator for the reverse mapping. For a simple model problem, we perform a systematic numerical study to assess the time scales on which the lifting errors disappear after reinitialization (healing); we also examine the effects on the simulated macroscopic behavior. The results show that, although in some cases accurate initialization of the higher order moments is not crucial, in general a detailed study of the lifting operator is required.
In this paper we are concerned with the dynamic bifurcation of the complex Swift-Hohenberg equation on a closed interval in $\mathbb R$. We consider the equations under the Dirichlet and the periodic boundary conditions. It is shown that the equation bifurcates from the trivial solution to an attractor when the control parameter crosses the critical value.
It is known that some predator-prey system can possess a unique limit cycle which is globally asymptotically stable. For a prototypical predator-prey system, we show that the solution curve of the limit cycle exhibits temporal patterns of a relaxation oscillator, or a Heaviside function, when certain parameter is small.
We examine the stability and instability of solutions of a polynomial difference equation with state-dependent Gaussian perturbations, and describe a phenomenon that can only occur in discrete time. For a particular set of initial values, we find that solutions approach equilibrium asymptotically in a highly regulated fashion: monotonically and bounded above by a deterministic sequence. We observe this behaviour with a probability that can be made arbitrarily high by choosing the initial value sufficiently small.
However, for any fixed initial value, the probability of instability is nonzero, and in fact we can show that as the magnitude of the initial value increases, the probability of instability approaches $1$.
The work focuses on the behaviour at infinity of solutions to second order elliptic equation with first order terms in a semi-infinite cylinder. Neumann's boundary condition is imposed on the lateral boundary of the cylinder and Dirichlet condition on its base. Under the assumption that the coefficients stabilize to a periodic regime, we prove the existence of a bounded solution, its stabilization to a constant, and provide necessary and sufficient condition for the uniqueness.
We study the Hamiltonian system of two point vortices, embedded in external strain and rotation. This external deformation field mimics the influence of neighboring vortices or currents in complex flows. When the external field is stationary, the equilibria of the two vortices, symmetric with respect to the center of the plane, are determined. The stability analysis indicates that two saddle points lie at the crossing of separatrices, which bound streamfunction lobes having neutral centers.
When the external field varies periodically with time, resonance becomes possible between the forcing and the oscillation of vortices around the neutral centers. A multiple time-scale expansion provides the slow-time evolution equation for these vortices, which, for weak periodic deformation, oscillate within their original (steady) trajectory. These analytical results accurately compare with numerical integration of the complete equations of motion. As the periodic deformation field increases, this vortex oscillation migrates out of the original trajectories, towards the location of the separatrices. With a periodic external field, these separatrices have given way to heteroclinic trajectories with multiple self-intersections, as shown by the calculation of the Melnikov function.
Chaos appears in vortex trajectories as they enter the aperiodic domain around the heteroclinic curves. In fact, this chaotic domain progressively fills out the plane, replacing KAM tori and cantori, as the periodic deformation field reaches finite amplitude. The appearance of windows of periodicity is illustrated.
We consider differential equations coupled with the input-output memory relation defined by the Preisach operator. The differential equation relates an instant value of the rate of change of the output of the Preisach operator with an instant value of its input. We propose an algorithm for the linearisation of the evolution operator of the system and apply it to define the characteristic multiplier of periodic solutions, which determines their stability. Examples of the system considered include models of terrestrial hydrology and electronic oscillators with hysteresis.
In this paper, we develop a pseudospectral method for differential equations defined on unbounded domains. We first introduce Gauss-type interpolations using a family of generalized Laguerre functions, and establish basic approximation results. Then we propose a pseudospectral method for differential equations on unbounded domains, whose coefficients may degenerate or grow up. As examples, we consider two model problems. The proposed schemes match the underlying problems properly and exhibit spectral accuracy. Numerical results demonstrate the efficiency of this new approach.
The four dimensional Rössler system is investigated. For this system the Poincaré map exhibits chaotic dynamics with two expanding directions and one strongly contracting direction. It is shown that the 16th iterate of this Poincaré map has a nontrivial invariant set on which it is semiconjugated to the full shift on two symbols. Moreover, it is proven that there exist infinitely many homoclinic and heteroclinic solutions connecting periodic orbits of period two and four, respectively. The proof utilizes the method of covering relations with smooth tools (cone conditions).
The proof is computer assisted - interval arithmetic is used to obtain bounds of the Poincaré map and its derivative.
We consider in this paper the stabilized semi-implicit (in time) scheme and the splitting scheme for the Allen-Cahn equation $\phi_t-\Delta\phi+$ε$^-2f(\phi)=0$ arising from phase transitions in material science. For the stabilized first-order scheme, we show that it is unconditionally stable and the error bound depends on ε-1 in some lower polynomial order using the spectrum estimate of [2, 10, 11]. In addition, the first- and second-order operator splitting schemes are proposed and the accuracy are tested and compared with the semi-implicit schemes numerically.
In this paper, we consider the $\Gamma$-convergence of Landau-Lifshitz ferromagnetic model in the presence of Bloch wall in the disk of $\bb R^2$ with Dirichlet boundary condition.
Add your name and e-mail address to receive news of forthcoming issues of this journal:
[Back to Top]