Discrete and Continuous Dynamical Systems - B
October 2021 , Volume 26 , Issue 10
Select all articles
In this paper, based on the transmission mechanism of echinococcosis in China, we propose a discrete-time dynamical model for the transmission of echinococcosis. The research results indicate that transmission dynamics of this discrete-time model are determined by basic reproduction number
A reliable and efficient a posteriori error estimator is presented for a weak Galerkin finite element method without stabilizer for the second order elliptic equation with mixed boundary conditions. The upper bound of the estimator is proved by Helmholtz decomposition technique and lower bound is hold naturally. The performance of the estimator is illustrated by numerical experiments.
In this paper, we propose a simple tumor-immune system model with time delay of tumor action, where two kinds of effects of the tumor cells (i.e. stimulation and neutralization) on the effector cells are considered. The local stability of the model is obtained by analyzing the characteristic equations of the model at the corresponding equilibria, the sufficient conditions on the global stability are found by applying the Fluctuation Lemma and constructing the different convergent sequences. The obtained results show that, compared to the results for the model without time delay, the time delay of tumor action can affect the stability of tumor equilibrium of the model as the stimulation effect of the tumor cells is strong enough, while the delay is harmless for the stability of tumor equilibrium under the neutralization of tumor cells. For the appropriate neutralization of tumor cells on effector cells, the bistability of the tumor free equilibrium and the stronger tumor equilibrium can appear. In the case of stimulation of tumor cells, the sufficiently large time delay can lead to the appearance of a stable periodic solution by Hopf bifurcation, and the numerical simulation illustrates that the amplitude of the periodic orbit increases with time delay. We also discuss the dependence of the tumor equilibrium and the time delay threshold, determining the stability of the tumor equilibrium, on tumor action. The related conditions determining dynamics of the model are expressed by certain formulae with biological meanings.
We consider an ecological model consisting of two species of predators competing for their common prey with explicit interference competition. With a proper rescaling, the model is portrayed as a singularly perturbed system with one fast (prey dynamics) and two slow variables (dynamics of the predators). The model exhibits a variety of rich and interesting dynamics, including, but not limited to mixed-mode oscillations (MMOs), featuring concatenation of small and large amplitude oscillations, relaxation oscillations and bistability between a semi-trivial equilibrium state and a coexistent oscillatory state. More interestingly, in a neighborhood of singular Hopf bifurcation, long lasting transient dynamics in the form of chaotic MMOs or relaxation oscillations are observed as the system approaches the periodic attractor born out of supercritical Hopf bifurcation or a semi-trivial equilibrium state respectively. The transient dynamics could persist for hundreds or thousands of generations before the ecosystem experiences a regime shift. The time series of population cycles with different types of irregular oscillations arising in this model stem from a biological realistic feature, namely, by the variation in the intraspecific competition amongst the predators. To explain these oscillations, we use bifurcation analysis and methods from geometric singular perturbation theory. The numerical continuation study reveals the rich bifurcation structure in the system, including the existence of codimension-two bifurcations such as fold-Hopf and generalized Hopf bifurcations.
Anti-angiogenesis therapy has been an emerging cancer treatment which may be further combined with chemotherapy to enhance overall survival of cancer patients. In this paper, we investigate a system of nonlinear ordinary differential equations describing a microenvironment consisting of host cells, tumor cells, immune cells and endothelial cells while incorporating treatment combination with chemotherapy and anti-angiogenesis therapy. We perform a dynamical systems analysis demonstrating that our model is able to capture the three phases of cancer immunoediting: elimination, equilibrium, and escape. In addition, we present transcritical bifurcations for relevant parameter values that correspond to the progression from the elimination phase to the equilibrium phase. A range of medically useful tumor doubling times were simulated to determine how combined therapy affects the tumor microenvironment over the course of a 250 day treatment. This analysis found two additional bifurcation parameters that move the system of equations from the equilibrium phase to the elimination phase. We determine that the most important aspect of an effective therapy is the activation of the anti-tumor immune response.
Layek and Pati (Phys. Lett. A, 2017) studied a nonlinear system of five coupled equations, which describe thermal relaxation in Rayleigh-Benard convection of a Boussinesq fluid layer, heated from below. Here we return to that paper and use techniques from dynamical systems theory to analyse the codimension-one Hopf bifurcation and codimension-two double-zero Bogdanov-Takens bifurcation. We determine the stability of the bifurcating limit cycle, and produce an unfolding of the normal form for codimension-two bifurcation for the Layek and Pati's model.
The singular support of the global attractor is introduced. It is shown that the singular support of the global attractor for a damped BBM equation equals to the singular support of the force term. This gives a delicate description of the local regularity, which roughly says that the attractor is smooth exactly where the force is smooth.
In this paper, we analyze the qualitative dynamics of a generalized Nosé-Hoover oscillator with two parameters varying in certain scope. We show that if a solution of this oscillator will not tend to the invariant manifold
In this paper, we establish a new online game addiction model with low and high risk exposure. With the help of the next generation matrix, the basic reproduction number
We study the Cauchy problem for the equations describing a viscous compressible and heat-conductive fluid in two dimensions. By imposing a weight function to initial density to deal with Sobolev embedding in critical space, and constructing an ad-hoc truncation to control the quadratic nonlinearity appeared in energy equation, we establish the local in time existence of unique strong solution with large initial data. The vacuum state at infinity or the compactly supported density is permitted. Moreover, we provide a different approach and slightly improve the weighted
We investigate the dynamics of a three-dimensional system modeling a molecular mechanism for the circadian rhythm in Drosophila. We first prove the existence of a compact attractor in the region with biological meaning. Under the assumption that the dimerization reactions are fast, in this attractor we reduce the three-dimensional system to a simpler two-dimensional system on the persistent normally hyperbolic slow manifold.
In this paper we prove that the stochastic Navier-Stokes equations with stable Lévy noise generate a random dynamical systems. Then we prove the existence of random attractor for the Navier-Stokes equations on 2D spheres under stable Lévy noise (finite dimensional). We also deduce the existence of a Feller Markov Invariant Measure.
We prove the norm inflation phenomena for the Boussinesq system on
We study for nonlinear Kirchhoff's model of pseudo parabolic type by considering its two different problems.
In this paper, by combining of fractional centered difference approach with alternating direction implicit method, we introduce a mixed difference method for solving two-dimensional Riesz space fractional advection-dispersion equation. The proposed method is a fourth order centered difference operator in spatial directions and second order Crank-Nicolson method in temporal direction. By reviewing the consistency and stability of the method, the convergence of the proposed method is achieved. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed technique.
We study positive solutions to steady state reaction diffusion equations of the form:
In this paper, we consider a class of spatially heterogeneous reaction diffusion rabies model which was used to describe population dynamics of the rabies epidemic disease observed in Europe. The dynamics of both the original non-degenerate reaction-diffusion system and its corresponding shadow system are investigated in great details. Firstly, we prove that under certain conditions, the in-time solutions of both the original non-degenerate reaction-diffusion system and its shadow system exist globally and remain uniformly bounded. Secondly, we are capable of showing that the shadow system is the nice approximations for the original non-degenerate reaction-diffusion system when the diffusion rate
We investigate numerically the advancing motion of 3D droplets spreading on physically flat chemically heterogeneous surfaces with periodic structures. We use the Navier-Stokes-Cahn-Hilliard equations with the generalized Navier boundary conditions to model the motion of droplets. Based on a convex splitting scheme, we have done numerical simulations and compared the results between different surface patterns quantitatively. We study the effect of pattern property on the advancing motion of three phase contact lines, the critical volume at the contact line jump and the effective advancing angles. By increasing the volume of droplet slowly on heterogeneous surfaces with different pattern property, we find that the advancing contact line approaches an equiangular octagon for the patterned surface with periodic squares separated by channels and approaches a regular hexagon for the patterned surface with periodic circles in regular hexagonal arrays. The shape of three-phase contact line is much more determined by the macro structure of the pattern than the micro structure of the pattern in each period.
Many fractional processes can be represented as an integral over a family of Ornstein–Uhlenbeck processes. This representation naturally lends itself to numerical discretizations, which are shown in this paper to have strong convergence rates of arbitrarily high polynomial order. This explains the potential, but also some limitations of such representations as the basis of Monte Carlo schemes for fractional volatility models such as the rough Bergomi model.
This paper deals with the number of limit cycles for planar piecewise smooth near-Hamiltonian or near-integrable systems with a switching curve. The main task is to establish a so-called first order Melnikov function which plays a crucial role in the study of the number of limit cycles bifurcated from a periodic annulus. We use the function to study Hopf bifurcation when the periodic annulus has an elementary center as its boundary. As applications, using the first order Melnikov function, we consider the number of limit cycles bifurcated from the periodic annulus of a linear center under piecewise linear polynomial perturbations with three kinds of quadratic switching curves. And we obtain three limit cycles for each case.
A mechanistic model describing the anaerobic mineralization of chlorophenol in a three-step food-web is investigated. The model is a six-dimensional system of ordinary differential equations. In our study, the phenol and the hydrogen inflowing concentrations are taken into account as well as the maintenance terms. The case of a large class of growth kinetics is considered, instead of specific kinetics. We show that the system can have up to eight types of steady states and we analytically determine the necessary and sufficient conditions for their existence according to the operating parameters. In the particular case without maintenance, the local stability conditions of all steady states are determined. The bifurcation diagram shows the behavior of the process by varying the concentration of influent chlorophenol as the bifurcating parameter. It shows that the system exhibits a bi-stability where the positive steady state can lose stability undergoing a supercritical Hopf bifurcation with the emergence of a stable limit cycle.
In this paper we are concerned with the existence of stable stationary solutions for the problem
This paper considers a stochastic single-species model with Lévy noises and time periodic coefficients. By Lyapunov functions and stochastic estimates, the threshold conditions between the time-average persistence in probability and extinction for the model are derived where Lévy noises play an important role in persistence and extinction of populations. It is shown that the time-average persistence in probability of the model implies the existence and uniqueness of positive periodic solution and the existence and uniqueness of periodic measure of the model. An example and its numerical simulations are given to verify the effectiveness of the theoretical results.
Chebyshev criterion is a powerful tool in the study of limit cycle bifurcations in dynamical systems based on Abelian integrals, but it is difficult when the Abelian integrals involve parameters. In this paper, we consider the Abelian integrals on the periodic annuli of a Hamiltonian with one parameter, arising from the generalized Liénard system, and identify the parameter values such that the Abelian integrals have Chebyshev property. In particular, the bounds on the number of zeros of the Abelian integrals are derived for different parameter intervals. The main mathematical tools are transformations and polynomial boundary theory, which overcome the difficulties in symbolic computations and analysis, arising from large parametric-semi-algebraic systems.
We examine the asymptotic behavior of the non-autonomous non-local fractional stochastic reaction-diffusion equations on
In this paper, we investigate the dynamical behavior of extinction and survival in tumor growth model with immunization under stochastic perturbation. Firstly, the model describing the growth of cancer cells monitored by immune cells is established. Then, the steady probability distribution of tumor cells for different noise intensities and immune parameter intensities, and necessary conditions for extinction and different survival of cancer cells are obtained by numerical and theoretical method. Besides, it is found that the extinction and survival of cancer cells rely on the state of immunization and noise. Finally, stochastic simulations are taken to test the theoretical analytical results. The results of our work are beneficial to discover the evolution mechanism and design effective immunotherapy of tumor.
Add your name and e-mail address to receive news of forthcoming issues of this journal:
[Back to Top]