The chaos and control of a food chain model supplying additional food to top-predator

The control and management of chaotic population is one of the main objectives for constructing mathematical model in ecology today. In this paper, we apply a technique of controlling chaotic predator–prey population dynamics by supplying additional food to top-predator. We formulate a three species predator–prey model supplying additional food to top-predator. Existence conditions and local stability criteria of equilibrium points are determined analytically. Persistence conditions for the system are derived. Global stability conditions of interior equilibrium point is calculated. Theoretical results are verified through numerical simulations. Phase diagram is presented for various quality and quantity of additional food. One parameter bifurcation analysis is done with respect to quality and quantity of additional food separately keeping one of them fixed. Using MATCONT package, we derive the bifurcation scenarios when both the parameters quality and quantity of additional food vary together. We predict the existence of Hopf point (H), limit point (LP) and branch point (BP) in the model for suitable supply of additional food. We have computed the regions of different dynamical behaviour in the quantity–quality parametric plane. From our study we conclude that chaotic population dynamics of predator prey system can be controlled to obtain regular population dynamics only by supplying additional food to top predator. This study is aimed to introduce a new non-chemical chaos control mechanism in a predator–prey system with the applications in fishery management and biological conservation of prey predator species.     

Applications of He’s principles to partial differential equations

This paper applies the variational iteration method (VIM) and semi-inverse variational principle to obtain solutions of linear and nonlinear partial differential equations. The nonlinear model is considered from gas dynamics, fluid dynamics and Burgers equation. The linear model is the heat transfer (diffusion) equation. Results show that variational iteration method is a powerful mathematical tool for solving linear and nonlinear partial differential equations, and therefore, can be widely applied to engineering problems.     

Integration of mussel in fish farm: Mathematical model and analysis

The paper deals with the dynamical behavior of fish and mussel population in a fish farm where external food is supplied. The ecosystem of the fish farm is represented by a set of nonlinear differential equations involving the nutrient (food), fish and mussels. We have studied the boundedness, local stability and global stability of the model system. We have incorporated the discrete type gestational delay of fish and analyze effect of the delay on the dynamical behavior of the model system. The delay parameter complicates the dynamics depending on the external food from changing the stable state to unstable damped periodic trajectories leading to a limit cycle oscillation. We have studied the Hopf-bifurcation of the model system in the neighborhood of the coexisting equilibrium point considering delay as a variable bifurcation parameter. We have performed numerical simulation to verify the analytical results. The entire study reveals that the external food supply controls the dynamics of the system.     

Evolving to the edge of chaos: Chance or necessity?

We show that ecological systems evolve to edges of chaos (EOC). This has been demonstrated by analyzing three diverse model ecosystems using numerical simulations in combination with analytical procedures. It has been found that all these systems reside on EOC and display short-term recurrent chaos (strc). The first two are non-linear food chains and the third one is a linear food chain. The dynamics of first two is dictated by deterministic changes in system parameters. In contrast to this, dynamics of the third model system (the linear food chain) is governed by both deterministic changes in system parameters as well as exogenous stochastic perturbations (unforeseen changes in initial conditions) of these dynamical systems.     

Evolving to the edge of chaos: Chance or necessity?

We show that ecological systems evolve to edges of chaos (EOC). This has been demonstrated by analyzing three diverse model ecosystems using numerical simulations in combination with analytical procedures. It has been found that all these systems reside on EOC and display short-term recurrent chaos (strc). The first two are non-linear food chains and the third one is a linear food chain. The dynamics of first two is dictated by deterministic changes in system parameters. In contrast to this, dynamics of the third model system (the linear food chain) is governed by both deterministic changes in system parameters as well as exogenous stochastic perturbations (unforeseen changes in initial conditions) of these dynamical systems.     

Development and Analysis of a Mathematics Aptitude Test for Gifted Elementary School Students

This paper describes the development and preliminary analysis of a mathematical test targeted for high mathematical ability elementary school students, the Stanford Education Program for Gifted Youth (EPGY) Mathematical Aptitude Test (SEMAT). A version was administered to 248 students, 9–11 years old, in EPGY. The SEMAT was developed because no other satisfactory test was designed or normed for this population. Most standardized tests assess mathematics proficiency for the general population so that gifted students' scores cluster in the few top percentiles. The SEMAT discriminated among this extreme upper end. Item response theory determined proficiency estimates, which were then used as scores to predict various outcomes in EPGY. The SEMAT proved to be a strong predictor of acceleration in EPGY.     

Using dynamic models to support inferences of insider threat risk

Two modeling approaches were integrated to address the problem of predicting the risk of an attack by a particular insider. We present a system dynamics model that incorporates psychological factors including personality, attitude and counterproductive behaviors to simulate the pathway to insider attack. Multiple runs of the model that sampled the population of possible personalities under different conditions resulted in simulated cases representing a wide range of employees of an organization. We then structured a Bayesian belief network to predict attack risk, incorporating important variables from the system dynamics model and learning the conditional probabilities from the simulated cases. Three scenarios were considered for comparison of risk indicators: An average employee (i.e., one who scores at the mean of a number of personality variables), an openly disgruntled malicious insider, and a disgruntled malicious insider who decides to conceal bad behaviors. The counterintuitive result is that employees who act out less than expected, given their particular level of disgruntlement, can present a greater risk of being malicious than other employees who exhibit a higher level of counterproductive behavior. This result should be tempered, however, considering the limited grounding of some of the model parameters. Nevertheless, this approach to integrating system dynamics modeling and Bayesian belief networks to address an insider threat problem demonstrates the potential for powerful prediction and detection capability in support of insider threat risk mitigation.     

Vortex dynamics for the nonlinear wave equation

Vortex dynamics for the nonlinear wave equation is a typical model of the “particle and field” theories of classical physics. The formal derivation of the dynamical law was done by J.Neu. He also made an interesting connection between vortex dynamics and the Dirac theory of electrons. Here we give a rigorous mathematical proof of this natural dynamical law. © 1999 John Wiley & Sons, Inc.     

Challenges of living in the harsh environments: A mathematical modeling study

We propose and analyze a mathematical model, which mimics community dynamics of plants and animals in harsh environments. The mathematical model exploits type IV functional responses whose idiosyncrasies have been recognized only in recent years. The interaction of the middle predator with the top predator is cast into Leslie-Gower scheme. Linear and non-linear stability analyses are performed to get an idea of the stability behavior of the model food chain. It turns out that carrying capacity of the prey and the immunity parameter of the middle predator are two crucial parameters governing the model. Availability of alternative food options to the generalist predator also plays a key role in deciding the model dynamics.Simulation runs performed on this model provide insight into population dynamics of monkeys of macaque family found in northern Japan. These monkeys are social animals which reproduce sexually. The characteristic feature of the model dynamics is that the generalist predator (macaque monkeys) is able to avoid impending extinction frequently and recovers at a rate which falsify threats from exogenous external forces; extreme weather conditions, etc.     

Mathematical modelling and computational analysis of enzymatic degradation of xenobiotic polymers

We study the enzymatic degradation of xenobiotic polymers mathematically. As a mathematical model, we derive a linear second-order hyperbolic partial differential equation which governs the evolution of the weight distribution with respect to the molecular weight. Given an initial weight distribution and a final weight distribution, we formulate a problem to determine a degradation rate. We establish a necessary and sufficient condition for which the problem has a local solution. We also introduce a numerical technique based on our analysis, and present a numerical result that we obtained applying weight distributions before and after enzymatic degradation of polyvinyl alcohol.     

Dynamics of an HIV-1 infection model with cell mediated immunity

In this paper, we study the dynamics of an improved mathematical model on HIV-1 virus with cell mediated immunity. This new 5-dimensional model is based on the combination of a basic 3-dimensional HIV-1 model and a 4-dimensional immunity response model, which more realistically describes dynamics between the uninfected cells, infected cells, virus, the CTL response cells and CTL effector cells. Our 5-dimensional model may be reduced to the 4-dimensional model by applying a quasi-steady state assumption on the variable of virus. However, it is shown in this paper that virus is necessary to be involved in the modeling, and that a quasi-steady state assumption should be applied carefully, which may miss some important dynamical behavior of the system. Detailed bifurcation analysis is given to show that the system has three equilibrium solutions, namely the infection-free equilibrium, the infectious equilibrium without CTL, and the infectious equilibrium with CTL, and a series of bifurcations including two transcritical bifurcations and one or two possible Hopf bifurcations occur from these three equilibria as the basic reproduction number is varied. The mathematical methods applied in this paper include characteristic equations, Routh–Hurwitz condition, fluctuation lemma, Lyapunov function and computation of normal forms. Numerical simulation is also presented to demonstrate the applicability of the theoretical predictions.     

Chaotic vibrations in nonlinear problems of bar structures

In the response of nonlinear mathematical models which describe vibrations of structural elements one could observe an irregular behaviour which is called chaos. Loss of the information on initial states in deterministic dynamical systems after a short time of theirs evolution, increasing amplitudes of displacements, velocities and accelerations, sensitive dependency on initial conditions makes chaos dangerous phenomenon in mechanics of construction. In this article quantitative (bifurcation diagrams, Poincare sections and Fourier power spectrum analysis) identication methods of the chaotic dynamics in geometrically nonlinear model of one DOF Mises truss are shown. Main goal of this article is to show and verify dangerous influence of chaos (in the engineering sense) on the analyzed structure. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On noise modeling in a nerve fibre

We present a novel mathematical approach to model noise in dynamical systems. We do so by considering the dynamics of a chain of diffusively coupled Nagumo cells affected by noise. We show that the noise in a variable representing the transmembrane current can be effectively modeled as fluctuations in the model parameters corresponding to electric resistance and capacitance of the membrane. These fluctuations may account for the interactions between the membrane and the surrounding (physiological) solution as well as for the thermal effects. The proposed approach to model noise in a nerve fibre is an alternative to the standard technique based on the consideration of additive stochastic current perturbation (the Langevin type equations) and differs from it in important mathematical aspects, particularly, it points out to the non-Markov dynamics of transmembrane potential. Our scheme relates to a time scale which is shorter than the relaxation times of involved physiological processes.     

SCALE AND TIME EFFECTS ON MATHEMATICAL MODELS FOR TRANSPORT IN THE ENVIRONMENT

The purpose of this paper is to analyse mathematical models used in environmental modelling.Following a brief survey of the development in modelling scale-and time-dependent dispersion processes in the environment,this paper compares three similarity solutions,one of which is a solution of the generalized Feller equation(GF)with fractal parameters,and the other two for the newly-developed generalized Fokker-Planck equation(GFP).The three solutions are derived with parameters having physical significance.Data from field experiments are used to verify the solutions.The analyses indicate that the solutions of both GF and GFP represent the physically meaningful natural processes,and simulate the realistic shapes of tracer breakthrough curves.     

Pseudospectra localizations and their applications

Usefulness of spectral localizations in analysis of various matrix properties, such as stability of dynamical systems, has led us to derive a pseudospectra localization technique using the ideas that come from diagonally dominant matrices. In such way, many theoretical and practical applications of pseudospectra (robust stability, transient behavior, nonnormal dynamics, etc.) can be linked with specific relationships between matrix entries. This allows one to understand certain phenomena that occur in practice better, as we show for the realistic model of soil energetic food web. The novelty of the presented results, therefore, lies not only in new mathematical formulations but also in the conceptual sense because it links stability with empirical data and their uncertainty limitations. Copyright © 2015 John Wiley & Sons, Ltd.     

Modeling treatment of cancer using virotherapy with generalized logistic growth of tumor cells

Abstract

Virotherapy is an effective strategy in cancer treatment. It eliminates tumor cells without harming the healthy cells. In this article, a deterministic mathematical model to understand the dynamics of tumor cells in response to virotherapy is formulated and analyzed by incorporating cytotoxic T lymphocytes (CTLs). The basic reproduction number and the immune response reproduction number are computed and different equilibria of the proposed model are found. The local stability of different equilibria is discussed in detail. Further, the proposed model is extended to stochastic model. Numerical simulation is performed for both deterministic and stochastic models. It is observed that when both the reproduction numbers are greater than one, which corresponds to existence of unique nontrivial equilibrium point, dynamics of deterministic and stochastic models are almost same. The deterministic model shows a very complex dynamics when one or both the reproduction numbers are below one. The system exhibits both backward bifurcation and Hopf-bifurcation for suitable sets of parameters and in this situation it is not easy to predict the dynamics of cancer cells and virus particles. The existence of backward bifurcation demonstrates the fact that partial success of virotherapy can be achieved even if the immune response reproduction number is less than one.     

非线性粘弹性梁的动力学行为

建立了描述受周期荷载作用的均匀粘弹性梁动力学行为的非线性偏微分-积分方程,梁的材料满足Leaderman非线性本构关系,对于两端简支的情形用Galerkin方法进行了2阶截断后,简化为常微分-积分方程,进一步简化为便于进行数值实验的常微分方程,最后用数值方法比较了1阶和2阶截断系统的动力学行为。     

Optimal harvesting policy of an inland fishery resource under incomplete information

A new mathematical model for finding the optimal harvesting policy of an inland fishery resource under incomplete information is proposed in this paper. The model is based on a stochastic control formalism in a regime‐switching environment. The incompleteness of information is due to uncertainties involved in the body growth rate of the fishery resource: a key biological parameter. Finding the most cost‐effective harvesting policy of the fishery resource ultimately reduces to solving a terminal and boundary value problem of a Hamilton‐Jacobi‐Bellman equation: a nonlinear and degenerate parabolic partial differential equation. A simple finite difference scheme for solving the equation is then presented, which turns out to be convergent and generates numerical solutions that comply with certain theoretical upper and lower bounds. The model is finally applied to the management of Plecoglossus altivelis, a major inland fishery resource in Japan. The regime switching in this case is due to the temporal dynamics of benthic algae, the main food of the fish. Model parameter values are identified from field measurement results in 2017. Our computational results clearly show the dependence of the optimal harvesting policy on the river environmental and biological conditions. The proposed model would serve as a mathematical tool for fishery resource management under uncertainties.     

中医方药量化研究中“相对药量”的数学模型体系

建立中医方药量华研究的"相对药量"概念模型体系;方法:运用微分方程理论;提出五种情况下,中药常用量范围内相对药量概念的数学模型,分别为直线模型、指数函数模型、对数函数模型、二次函数模型(开口向上和向下两种),并说明该模型体系的和理性与适应性;结论:相对药量概念核心体系的建立,增加了中医方药"相对药量"可比性的全面性,这对进一步研究单位药乃至方剂中各药在性、味、归经等方面的影响程度及其规律性,将起到至关重要的作用.