Analysis of incomplete factorizations for a nine-point approximation to a convection–diffusion model problem

We study the stability of zero-fill incomplete LU factorizations of a nine-point coefficient matrix arising from a high-order compact discretisation of a two-dimensional constant-coefficient convection–diffusion problem. Nonlinear recurrences for computing entries of the lower and upper triangular matrices are derived and we show that the sequence of diagonal entries of the lower triangular factor is unconditionally convergent. A theoretical estimate of the limiting value is derived and we show that this estimate is a good predictor of the computed value. The unconditional convergence of the diagonal sequence of the lower triangular factor to a positive limit implies that the incomplete factorization process never encounters a zero pivot and that the other diagonal sequences are also convergent. The characteristic polynomials associated with the lower and upper triangular solves that occur during the preconditioning step are studied and conditions for the stability of the triangular solves are derived in terms of the entries of the tridiagonal matrices appearing in the lower and upper subdiagonals of the block triangular system matrix and a triplet of parameters which completely determines the solution of the nonlinear recursions. Results of ILU-preconditioned GMRES iterations and the effects of orderings on their convergence are also described.     

Analysis of a delayed stochastic predator–prey model in a polluted environment

In this paper, we investigate two-species Lotka–Volterra delayed stochastic predator–prey systems, with and without pollution, denoted by (M)$(M)$ and (_M0)$(M_0)$, respectively. We show that there exists a unique non-negative solution in each system that is permanent in time average under certain conditions. Moreover, the non-permanence of model (M)$(M)$ is studied. Finally, computer simulations are carried out to verify our results.     

Analysis of reinforced concrete structures subjected to dynamic loads with a viscoplastic Drucker–Prager model

The behavior of reinforced concrete structures subjected to dynamic loads is analyzed. The concrete material is modelled by an elasto-viscoplastic law, whose inviscid counterpart is the Drucker–Prager model. A viscous regularization is introduced in order to avoid the mesh dependency effects that usually appear when strain softening occurs. The model is implemented in a general finite element computer code for fast transient analysis of fluid-structure systems, based on an explicit central difference scheme. The model is activated to both continuum elements and layered shell elements. So, realistic numerical analyses of complex 3-D engineering problems are simple and efficient. Three examples, two of which are modelled with layered shell elements, are presented below.     

Analysis of a Cahn–Hilliard–Brinkman model for tumour growth with chemotaxis

Phase field models recently gained a lot of interest in the context of tumour growth models. Typically Darcy-type flow models are coupled to Cahn–Hilliard equations. However, often Stokes or Brinkman flows are more appropriate flow models. We introduce and mathematically analyse a new Cahn–Hilliard–Brinkman model for tumour growth allowing for chemotaxis. Outflow boundary conditions are considered in order not to influence tumour growth by artificial boundary conditions. Existence of global-in-time weak solutions is shown in a very general setting.     

Analysis of a diffusive Leslie–Gower predator–prey model with nonmonotonic functional response

In this paper, a diffusive Leslie–Gower predator–prey system with nonmonotonic functional respond is studied. We obtain the persistence of this model and show the local asymptotic stability of positive constant equilibrium by linearized analysis and the global stability by constructing Liapunov function. Besides, Turing instability of this equilibrium is obtained. The existence and nonexistence of positive nonconstant steady states of this model are established. Furthermore, by numerical simulations we illustrate the patterns of prey and predator.     

Analysis on dynamics of a population model with predator–prey-dependent functional response

We consider a reaction–diffusion population model with predator–prey-dependent functional response. Firstly, we discuss the conditions which ensure the model has a unique positive constant solution. Secondly, we investigate the dynamical properties of the model, including the large time behaviors of the nonconstant solutions and the local and the global asymptotic stability of the positive constant solution.     

Analysis of a plankton–fish model with external toxicity and nonlinear harvesting

Ricerche di Matematica - In this paper, we consider a three species plankton–fish system that incorporates external toxicity and nonlinear harvesting. We consider that the growth of species...     

Harvesting of a phytoplankton–zooplankton model

Considering that some phytoplankton and zooplankton are harvested for food, a phytoplankton–zooplankton model with harvesting is proposed and investigated. First, stability conditions of equilibria and existence conditions of a Hopf-bifurcation are established. Our results indicate that over exploitation would result in the extinction of the population and an appropriate harvesting strategy should ensure the sustainability of the population which is in line with reality. Furthermore, the existence of bionomic equilibria and the optimal harvesting policy are discussed. The present value of revenues is maximized by using Pontryagin’s maximum principle subject to the state equations and the control constraints. We discussed the case of optimal equilibrium solution. It is found that the shadow prices remain constant over time in optimal equilibrium when they satisfy the transversality condition. It is established that the zero discounting leads to the maximization of economic revenue and that an infinite discount rate leads to complete dissipation of economic rent. Finally, some numerical simulations are given to illustrate our results.     

Uniform stabilization to equilibrium of a nonlinear fluid–structure interaction model

We consider uniform stability to a nontrivial equilibrium of a nonlinear fluid–structure interaction (FSI) defined on a two or three dimensional bounded domain. Stabilization is achieved via boundary and/or interior feedback controls implemented on both the fluid and the structure. The interior damping on the fluid combining with the viscosity effect stabilizes the dynamics of fluid. However, this dissipation propagated from the fluid alone is not sufficient to drive uniformly to equilibrium the entire coupled system. Therefore, additional interior damping on the wave component or boundary porous like damping on the interface is considered. A geometric condition on the interface is needed if only boundary damping on the wave is active. The main technical difficulty is the mismatch of regularity of hyperbolic and parabolic component of the coupled system. This is overcome by considering special multipliers constructed from Stokes solvers. The uniform stabilization result obtained in this article is global for the fully coupled FSI model.     

Analysis of a predator–prey model with modified Holling–Tanner functional response and time delay

In paper, a predator–prey model with modified Holling–Tanner functional response and time delay is discussed. It is proved that the system is permanent under some appropriate conditions. The local stability of the equilibria is investigated. By constructing a suitable Lyapunov functional, sufficient conditions are derived for the global stability of the positive equilibrium of the model.     

Dynamics of a generalized Ricker–Beverton–Holt competition model subject to Allee effects

In this article, we propose and study a generalized Ricker–Beverton–Holt competition model subject to Allee effects to obtain insights on how the interplay of Allee effects and contest competition affects the persistence and the extinction of two competing species. By using the theory of monotone dynamics and the properties of critical curves for non-invertible maps, our analysis show that our model has relatively simple dynamics, i.e. almost every trajectory converges to a locally asymptotically stable equilibrium if the intensity of intra-specific competition intensity exceeds that of inter-specific competition. This equilibrium dynamics is also possible when the intensity of intra-specific competition intensity is less than that of inter-specific competition but under conditions that the maximum intrinsic growth rate of one species is not too large. The coexistence of two competing species occurs only if the system has four interior equilibria. We provide an approximation to the basins of the boundary attractors (i.e. the extinction of one or both species) where our results suggests that contest species are more prone to extinction than scramble ones are at low densities. In addition, in comparison to the dynamics of two species scramble competition models subject to Allee effects, our study suggests that (i) Both contest and scramble competition models can have only three boundary attractors without the coexistence equilibria, or four attractors among which only one is the persistent attractor, whereas scramble competition models may have the extinction of both species as its only attractor under certain conditions, i.e. the essential extinction of two species due to strong Allee effects; (ii) Scramble competition models like Ricker type models can have much more complicated dynamical structure of interior attractors than contest ones like Beverton–Holt type models have; and (iii) Scramble competition models like Ricker type competition models may be more likely to promote the coexistence of two species at low and high densities under certain conditions: At low densities, weak Allee effects decrease the fitness of resident species so that the other species is able to invade at its low densities; While at high densities, scramble competition can bring the current high population density to a lower population density but is above the Allee threshold in the next season, which may rescue a species that has essential extinction caused by strong Allee effects. Our results may have potential to be useful for conservation biology: For example, if one endangered species is facing essential extinction due to strong Allee effects, then we may rescue this species by bringing another competing species subject to scramble competition and Allee effects under certain conditions.     

On solutions to a FitzHugh–Rinzel type model

Ricerche di Matematica - A ternary autonomous dynamical system of FitzHugh–Rinzel type is analyzed. The system, at start, is reduced to a nonlinear integro differential equation. The...     

Analysis of spatiotemporal patterns in a single species reaction–diffusion model with spatiotemporal delay

Employing the theories of Turing bifurcation in the partial differential equations, we investigate the dynamical behavior of a single species reaction–diffusion model with spatiotemporal delay. The linear stability and the conditions for the occurrence of Turing bifurcation in this model are obtained. Moreover, the amplitude equations which represent different spatiotemporal patterns are also obtained near the Turing bifurcation point by using multiple scale method. In Turing space, it is found that the spatiotemporal distributions of the density of this researched species have spots pattern and stripes pattern. Finally, some numerical simulations corresponding to the different spatiotemporal patterns are given to verify our theoretical analysis.     

A new approach to the stabilization of a fluid–structure interaction model

This article is a continuation of our work on a linear fluid–structure interaction model [Grobbelaar-Van Dalsen, On a fluid–structure model in which the dynamics of the structure involves the shear stress due to the fluid, J. Math. Fluid Mech. 10(3) (2008), pp. 388–401; Grobbelaar-Van Dalsen, Strong stability for a fluid––structure model, Math. Methods Appl. Sci., 32(2009) pp. 1452–1466]. The model describes the interaction between a 3-D incompressible fluid and a 2-D plate, the interface, which coincides with a flat flexible part of the surface of the vessel containing the fluid. The mathematical model comprises the Stokes equations and the equations for the longitudinal deflections of the plate with the inclusion of the shear stress that the fluid exerts on the plate. A dissipative damping mechanism of Kelvin–Voigt type is applied to the interior of the plate. While our earlier work shows that weak solutions in a space of finite energy are strongly asymptotically stable under no-slip transmission conditions at the interface with uniform exponential stability only attainable under an additional domination condition, the present research is directed at achieving uniform exponential stability of weak solutions without imposing the domination condition. Using energy methods we establish uniform exponential decay under a modified transmission condition at the interface. This condition entails that the fluid velocity at the interface is coupled to a linear combination of the plate velocity and displacement.