Discrete & Continuous Dynamical Systems - B
January 2021 , Volume 26 , Issue 1
20 years of DCDS-B
Select all articles
We study local bifurcations of periodic solutions to time-periodic (systems of) integrodifference equations over compact habitats. Such infinite-dimensional discrete dynamical systems arise in theoretical ecology as models to describe the spatial dispersal of species having nonoverlapping generations. Our explicit criteria allow us to identify branchings of fold- and crossing curve-type, which include the classical transcritical-, pitchfork- and flip-scenario as special cases. Indeed, not only tools to detect qualitative changes in models from e.g. spatial ecology and related simulations are provided, but these critical transitions are also classified. In addition, the bifurcation behavior of various time-periodic integrodifference equations is investigated and illustrated. This requires a combination of analytical methods and numerical tools based on Nyström discretization of the integral operators involved.
This paper is concerned with a numerical solution of the acoustic scattering by a bounded impenetrable obstacle in three dimensions. The obstacle scattering problem is formulated as a boundary value problem in a bounded domain by using a Dirichlet-to-Neumann (DtN) operator. An a posteriori error estimate is derived for the finite element method with the truncated DtN operator. The a posteriori error estimate consists of the finite element approximation error and the truncation error of the DtN operator, where the latter is shown to decay exponentially with respect to the truncation parameter. Based on the a posteriori error estimate, an adaptive finite element method is developed for the obstacle scattering problem. The truncation parameter is determined by the truncation error of the DtN operator and the mesh elements for local refinement are marked through the finite element approximation error. Numerical experiments are presented to demonstrate the effectiveness of the proposed method.
In this paper, we present and study
In this paper, we discretize and analyze a stoichiometric optimal foraging model where the grazer's feeding effort depends on the producer's nutrient quality. We systematically make comparisons of the dynamical behaviors between the discrete-time model and the continuous-time model to study the robustness of model predictions to time discretization. When the maximum growth rate of producer is low, both model types admit similar dynamics including bistability and deterministic extinction of the grazer caused by low nutrient quality of the producer. Especially, the grazer is benefited from optimal foraging similarly in both discrete-time and continuous-time models. When the maximum growth rate of producer is high, dynamics of the discrete-time model are more complex including chaos. A phenomenal observation is that under extremely high light intensities, the grazer in the continuous-time model tends to perish due to poor food quality, however, the grazer in the discrete-time model persists in regular or irregular oscillatory ways. This significant difference indicates the necessity of studying discrete-time models which naturally include species' generations and are thus more popular in theoretical biology. Finally, we discuss how the shape of the quality-based feeding function regulates the beneficial or restraint effect of optimal foraging on the grazer population.
In the article, some bilinear evolution equations in Hilbert space driven by paths of low regularity are considered and solved explicitly. The driving paths are scalar-valued and continuous, and they are assumed to have a finite
We present a spatial food chain model on a bounded domain coupled with optimal control theory to examine harvesting strategies. Motivated by the fishery industry in the Black Sea, the anchovy stock and two more trophic levels are modeled using nonlinear parabolic partial differential equations with logistic growth, movement by diffusion and advection, and Neumann boundary conditions. Necessary conditions for the optimal harvesting control are established. The objective for the problem is to find the spatial optimal harvesting strategy that maximizes the discounted net value of the anchovy population. Numerical simulations using data from the Black Sea are presented. We discuss spatial features of harvesting and the effects of diffusion and advection (migration speed) on the anchovy population. We also present the landing of anchovy and its net profit by applying two different harvesting strategies.
In this paper, we analyze the interaction of localized patterns such as traveling wave solutions for reaction-diffusion systems with nonlocal effect in one space dimension. We consider the case that a nonlocal effect is given by the convolution with a suitable integral kernel. At first, we deduce the equation describing the movement of interacting localized patterns in a mathematically rigorous way, assuming that there exists a linearly stable localized solution for general reaction-diffusion systems with nonlocal effect. When the distances between localized patterns are sufficiently large, the motion of localized patterns can be reduced to the equation for the distances between them. Finally, using this equation, we analyze the interaction of front solutions to some nonlocal scalar equation. Under some assumptions, we can show that the front solutions are interacting attractively for a large class of integral kernels.
A discrete delay is included to model the time between the capture of the prey and its conversion to viable biomass in the simplest classical Gause type predator-prey model that has equilibrium dynamics without delay. As the delay increases from zero, the coexistence equilibrium undergoes a supercritical Hopf bifurcation, two saddle-node bifurcations of limit cycles, and a cascade of period doublings, eventually leading to chaos. The resulting periodic orbits and the strange attractor resemble their counterparts for the Mackey-Glass equation. Due to the global stability of the system without delay, this complicated dynamics can be solely attributed to the introduction of the delay. Since many models include predator-prey like interactions as submodels, this study emphasizes the importance of understanding the implications of overlooking delay in such models on the reliability of the model-based predictions, especially since temperature is known to have an effect on the length of certain delays.
To underpin the concern that environmental change can flip an ecosystem from stable persistence to sudden total collapse, we consider a class of so-called ecoepidemic models, predator – prey/host – parasite systems, in which a base species is prey to a predator species and host to a micro-parasite species. Our model uses generalized frequency-dependent incidence for the disease transmission and mass action kinetics for predation.
We show that a large variety of dynamics can arise, ranging from dynamic persistence of all three species to either total ecosystem collapse caused by high transmissibility of the parasite on the one hand or to parasite extinction and prey-predator survival due to low parasite transmissibility on the other hand. We identify a threshold parameter (tipping number) for the transition of the ecosystem from uniform prey/host persistence to total extinction under suitable initial conditions.
The question of whether and how two competing consumers can coexist on a single limiting resource has a long tradition in ecological theory. We build on a recent seasonal (hybrid) model for one consumer and one resource, and we extend it by introducing a second consumer. Consumers reproduce only once per year, the resource reproduces throughout the"summer" season. When we use linear consumer reproduction between years, we find explicit expressions for the trivial and semi-trivial equilibria, and we prove that there is no positive equilibrium generically. When we use non-linear consumer reproduction, we determine conditions for which both semi-trivial equilibria are unstable. We prove that a unique positive equilibrium exists in this case, and we find an explicit analytical expression for it. By linear analysis and numerical simulation, we find bifurcations from the stable equilibrium to population cycles that may appear through period-doubling or Hopf bifurcations. We interpret our results in terms of climate change that changes the length of the"summer" season.
A complete Lyapunov function characterizes the behaviour of a general discrete-time dynamical system. In particular, it divides the state space into the chain-recurrent set where the complete Lyapunov function is constant along trajectories and the part where the flow is gradient-like and the complete Lyapunov function is strictly decreasing along solutions. Moreover, the level sets of a complete Lyapunov function provide information about attractors, repellers, and basins of attraction.
We propose two novel classes of methods to compute complete Lyapunov functions for a general discrete-time dynamical system given by an iteration. The first class of methods computes a complete Lyapunov function by approximating the solution of an ill-posed equation for its discrete orbital derivative using meshfree collocation. The second class of methods computes a complete Lyapunov function as solution of a minimization problem in a reproducing kernel Hilbert space. We apply both classes of methods to several examples.
We study the well-posedness for the mildly compressible Navier-Stokes-Cahn-Hilliard system with non-constant viscosity and Landau potential in two and three dimensional domains.
We consider the ecological and evolutionary dynamics of a reaction-diffusion-advection model for populations residing in a one-dimensional advective homogeneous environment, with emphasis on the effects of boundary conditions and domain size. We assume that there is a net loss of individuals at the downstream end with rate
For the two-species competition model, we show that the infinite dispersal rate is evolutionarily stable for
This paper is concerned with a semi-linear elliptic problem with Robin boundary condition:
The coevolution or coexistence of multiple viruses with multiple hosts has been an important issue in viral ecology. This paper is to study the mathematical properties of the solutions of a chemostat model for two host species and two virus species. By virtue of the global dynamics of its submodels and the theories of uniform persistence and Hopf bifurcation, we derive sufficient conditions for the coexistence of two hosts with two viruses and coexistence of two hosts with one virus, as well as occurrence of Hopf bifurcation.
We provide a short review of existing models with multiple taxis performed by (at least) one species and consider a new mathematical model for tumor invasion featuring two mutually exclusive cell phenotypes (migrating and proliferating). The migrating cells perform nonlinear diffusion and two types of taxis in response to non-diffusing cues: away from proliferating cells and up the gradient of surrounding tissue. Transitions between the two cell subpopulations are influenced by subcellular (receptor binding) dynamics, thus conferring the setting a multiscale character.
We prove global existence of weak solutions to a simplified model version and perform numerical simulations for the full setting under several phenotype switching and motility scenarios. We also compare (via simulations) this model with the corresponding haptotaxis-chemotaxis one featuring indirect chemorepellent production and provide a discussion about possible model extensions and mathematical challenges.
We review and discuss various aspects that the modeling of pharmacometric properties has on the structure of optimal solutions in mathematical models for cancer treatment. These include (i) the changes in the interpretation of the solutions as pharmacokinetic (PK) models are added, respectively deleted from the modeling and (ii) qualitative changes in the structures of optimal controls that occur as pharmacodynamic (PD) models are varied. The results will be illustrated with a sample of models for cancer treatment.
This paper is aimed to provide a new theory for the formation of the solar surface eruptions and sunspots. The key ingredient of the study is the new anti-diffusive effect of heat, based on the recently developed statistical theory of heat by the authors [
Viral dynamics within plant hosts can be important for understanding plant disease prevalence and impacts. However, few mathematical modeling efforts aim to characterize within-plant viral dynamics. In this paper, we derive a simple system of delay differential equations that describes the spread of infection throughout the plant by barley and cereal yellow dwarf viruses via the cell-to-cell mechanism. By incorporating ratio-dependent incidence function and logistic growth of the healthy cells, the model can capture a wide range of biologically relevant phenomena via the disease-free, endemic, mutual extinction steady states, and a stable periodic orbit. We show that when the basic reproduction number is less than
In this paper, we review some recent results on the nonlinear dynamics of delayed differential equation models describing the interaction between tumor cells and effector cells of the immune system, in which the delays represent times necessary for molecule production, proliferation, differentiation of cells, transport, etc. First we consider a tumor-immune system interaction model with a single delay and present results on the existence and local stability of equilibria as well as the existence of Hopf bifurcation in the model when the delay varies. Second we investigate a tumor-immune system interaction model with two delays and show that the model undergoes various possible bifurcations including Hopf, Bautin, Fold-Hopf (zero-Hopf), and Hopf-Hopf bifurcations. Finally we discuss a tumor-immune system interaction model with three delays and demonstrate that the model exhibits more complex behaviors including chaos. Numerical simulations are provided to illustrate the nonlinear dynamics of the delayed tumor-immune system interaction models. More interesting issues and questions on modeling and analyzing tumor-immune dynamics are given in the discussion section.
In the study of equilibrium solutions for partial differential equations there are so many equilibria that one cannot hope to find them all. Therefore one usually concentrates on finding individual branches of equilibrium solutions. On the one hand, a rigorous theoretical understanding of these branches is ideal but not generally tractable. On the other hand, numerical bifurcation searches are useful but not guaranteed to give an accurate structure, in that they could miss a portion of a branch or find a spurious branch where none exists. In a series of recent papers, we have aimed for a third option. Namely, we have developed a method of computer-assisted proofs to prove both existence and isolation of branches of equilibrium solutions. In the current paper, we extend these techniques to the Ohta-Kawasaki model for the dynamics of diblock copolymers in dimensions one, two, and three, by giving a detailed description of the analytical underpinnings of the method. Although the paper concentrates on applying the method to the Ohta-Kawasaki model, the functional analytic approach and techniques can be generalized to other parabolic partial differential equations.
In 1993, Holt and Lawton introduced a stochastic model of two host species parasitized by a common parasitoid species. We introduce and analyze a generalization of these stochastic difference equations with any number of host species, stochastically varying parasitism rates, stochastically varying host intrinsic fitnesses, and stochastic immigration of parasitoids. Despite the lack of direct, host density-dependence, we show that this system is dissipative i.e. enters a compact set in finite time for all initial conditions. When there is a single host species, stochastic persistence and extinction of the host is characterized using external Lyapunov exponents corresponding to the average per-capita growth rates of the host when rare. When a single host persists, say species
We develop in this paper efficient and accurate numerical schemes based on the scalar auxiliary variable (SAV) approach for the generalized Zakharov system and generalized vector Zakharov system. These schemes are second-order in time, linear, unconditionally stable, only require solving linear systems with constant coefficients at each time step, and preserve exactly a modified Hamiltonian. Ample numerical results are presented to demonstrate the accuracy and robustness of the schemes.
In this paper, we present a rigorous mathematical analysis of a free boundary problem modeling the growth of a vascular solid tumor with a necrotic core. If the vascular system supplies the nutrient concentration
The asymptotic behavior of a model for 2D incompressible stochastic micropolar fluid flows with rough noise on a Poincaré domain is investigated. First, the existence and uniqueness of solutions to an evolution equation arising from the underlying stochastic micropolar fluid model is established via the Galerkin method and energy method. Then the existence of a random attractor is studied by using the theory of random dynamical systems for which the noise is dealt with by appropriate reproducing kernel Hilbert space.
We present a review of the different techniques available to study a special kind of fractal basins of attraction known as Wada basins, which have the intriguing property of having a single boundary separating three or more basins. We expose several approaches to identify this topological property that rely on different, but not exclusive, definitions of the Wada property.
Searching saddle points on the potential energy surface is a challenging problem in the rare event. When there exist multiple saddle points, sampling different initial guesses are needed in local search methods in order to find distinct saddle points. In this paper, we present a novel global optimization-based dimer method (GOD) to efficiently search saddle points by coupling ant colony optimization (ACO) algorithm with optimization-based shrinking dimer (OSD) method. In particular, we apply OSD method as a local search algorithm for saddle points and construct a pheromone function in ACO to update the global population. By applying a two-dimensional example and a benchmark problem of seven-atom island on the (111) surface of an FCC crystal, we demonstrate that GOD shows a significant improvement in computational efficiency compared with OSD method. Our algorithm is the first try to apply the global optimization technique in searching saddle points and it offers a new framework to open up possibilities of adopting other global optimization methods.
Add your name and e-mail address to receive news of forthcoming issues of this journal:
[Back to Top]