Search

• Development and study of a synergistic co-infection model for crops. • Basic and Invasion Reproduction numbers drive the dynamics of the system. • Application to the Maize Lethal Necrosis case. • Among all possible control strategies, early roguing is the best one. Crops are often subject to intense attacks by pests and diseases. Among them, Maize Lethal Necrosis (MLN) is a serious disease that impact maize crops in many Southern countries. It results from the synergistic interaction of two plant viruses, transmitted by two vectors. In this paper, we develop a general deterministic epidemic model of crop-vector-borne disease for synergistic co-infection, with a particular focus on MLN disease. A theoretical analysis shows that different thresholds exists that drive the dynamics of the system: the well known basic reproduction numbers and also invasion reproduction numbers. The latter are essential for the emergence or not of the MLN disease. After a global sensitivity analysis, we illustrate our results through numerical simulations and discuss potential control methods such as vector control and roguing. [ABSTRACT FROM AUTHOR]

In this work, we consider numerical solutions of the FitzHugh–Nagumo system of equations describing the propagation of electrical signals in nerve axons. The system consists of two coupled equations: a nonlinear partial differential equation and a linear ordinary differential equation. We begin with a review of the qualitative properties of the nonlinear space independent system of equations. The subequation approach is applied to derive dynamically consistent schemes for the submodels. This is followed by a consistent and systematic merging of the subschemes to give three explicit nonstandard finite difference schemes in the limit of fast extinction and slow recovery. A qualitative study of the schemes together with the error analysis is presented. Numerical simulations are given to support the theoretical results and verify the efficiency of the proposed schemes. [ABSTRACT FROM AUTHOR]

A mathematical model presented in Berge T, Lubuma JM-S, Moremedi GM, Morris N Shava RK, A simple mathematical model for Ebola in Africa, J Biol Dyn 11(1): 42-74 (2016) for the transmission dynamics of Ebola virus is extended to incorporate vaccination and change of behavior for self-protection of susceptible individuals. In the new setting, it is shown that the disease-free equilibrium is globally asymptotically stable when the basic reproduction number is less than or equal to unity and unstable when . In the latter case, the model system admits at least one endemic equilibrium point, which is locally asymptotically stable. Using the parameters relevant to the transmission dynamics of the Ebola virus disease, we give sensitivity analysis of the model. We show that the number of infectious individuals is much smaller than that obtained in the absence of any intervention. In the case of the mass action formulation with vaccination and education, we establish that the number of infectious individuals decreases as the intervention efforts increase. In the new formulation, apart from supporting the theory, numerical simulations of a nonstandard finite difference scheme that we have constructed suggests that the results on the decrease of the number of infectious individuals is valid. [ABSTRACT FROM AUTHOR]

The Schrödinger equation is a model for many physical processes in quantum physics. It is a singularly perturbed differential equation where the presence of the small reduced Planck's constant makes the classical numerical methods very costly and inefficient. We design two new schemes. The first scheme is the nonstandard finite volume method, whereby the perturbation term is approximated by nonstandard technique, the potential is approximated by its mean value on the cell and the complex dependent boundary conditions are handled by exact schemes. In the second scheme, the deficiency of classical schemes is corrected by the inner expansion in the boundary layer region. Numerical simulations supporting the performance of the schemes are presented. [ABSTRACT FROM PUBLISHER]

We design explicit nonstandard finite difference schemes for the nonlinear Allen–Cahn reaction diffusion equation in the limit of very small interaction length. In the proposed scheme, the perturbation parameter is part of the argument of the functional step size, thereby minimizing the restrictions normally associated with standard explicit finite difference schemes. The derivation involves splitting the equation into the space-independent and the time-independent different models. An exact nonstandard scheme is proposed for the space-independent model and energy conservative schemes are proposed for the time-independent model. We show the power of the derived scheme over the existing schemes through several numerical examples. [ABSTRACT FROM AUTHOR]

We present a model of biofilm growth in a long channel where the biomass is assumed to have the rheology of a viscous polymer solution. We examine the competition between growth and erosion-like surface detachment due to the flow. A particular focus of our investigation is the effect of the biofilm growth on the fluid flow in the pores, and the issue of whether biomass can grow sufficiently to shut off fluid flow through the pores, thus clogging the pore space. Net biofilm growth is coupled along the pore length via flow rate and nutrient transport in the pore flow. Our 2D model extends existing results on stability of 1D steady state biofilm thicknesses to show that, in the case of flows driven by a fixed pressure drop, full clogging of the pore can indeed happen in certain cases dependent on the functional form of the detachment term. [ABSTRACT FROM AUTHOR]

Mathematical modelling of transport phenomena in food processes is vital to understand the process dynamics. In this work, we study the process of double sided cooking of meat by developing a mathematical model for the simultaneous heat and mass transfer. The constitutive equations for the heat and mass transport are based on Fourier conduction, and the Flory–Huggins theory respectively, formulated for a two-phase transport inside a porous medium. We investigate a reduced one-dimensional case to verify the model, by applying appropriate boundary conditions. The results of the simulation agree well with experimental findings reported in literature. Finally, we comment upon the sensitivity of the model to the porosity of meat. [ABSTRACT FROM AUTHOR]

Soft condensed matter (SCM) physics has recently gained importance for a large class of engineering materials. The treatment of food materials from a soft matter perspective, however, is only at the surface and is gaining importance for understanding the complex phenomena and structure of foods. In this work, we present a theoretical treatment of navy beans from a SCM perspective to describe the hydration kinetics. We solve the transport equations within a porous matrix and employ the Flory–Huggin’s equation for polymer–solvent mixture to balance the osmotic pressure. The swelling of the legume seed is modelled as a moving boundary with an explicit transient equation. The model exhibits a good agreement with the experimental observations and is capable of explaining the stages of hydration. Sensitivity analysis indicated that the degree of hydration is dependent on the bean size and is also sensitive to the selection of the intrinsic permeability of the bean. [ABSTRACT FROM AUTHOR]

We present a theory for the coupled flow of ice, subglacial water and subglacial sediment, which is designed to represent the processes which occur at the bed of an ice sheet. The ice is assumed to flow as a Newtonian viscous fluid, the water can flow between the till and the ice as a thin film, which may thicken to form streams or cavities, and the till is assumed to be transported, either through shearing by the ice, squeezing by pressure gradients in the till, or by fluvial sediment transport processes in streams or cavities. In previous studies, it was shown that the dependence of ice sliding velocity on effective pressure provided a mechanism for the generation of bedforms resembling ribbed moraine, while the dependence of fluvial sediment transport on water film depth provides a mechanism for the generation of bedforms resembling mega-scale glacial lineations. Here, we combine these two processes in a single model, and show that, depending largely on the granulometry of the till, instability can occur in a range of types which range from ribbed moraine through three-dimensional drumlins to mega-scale glacial lineations. [ABSTRACT FROM AUTHOR]

Purpose – For a partial differential equation with a fourth-order derivative such as the Cahn-Hilliard equation, it is always a challenge to design numerical schemes that can handle the restrictive time step introduced by this higher order term. The purpose of this paper is to employ a fractional splitting method to isolate the convective, the nonlinear second-order and the fourth-order differential terms. Design/methodology/approach – The full equation is then solved by consistent schemes for each differential term independently. In addition to validating the second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution. Findings – The scheme is second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution. Originality/value – The authors believe that this is the first time the equation is handled numerically using the fractional step method. Apart from the fact that the fractional step method substantially reduces computational time, it has the advantage of simplifying a complex process efficiently. This method permits the treatment of each segment of the original equation separately and piece them together, in a way that will be explained shortly, without destroying the properties of the equation. [ABSTRACT FROM AUTHOR]

We provide and analyse a model for the growth of bacterial biofilms based on the concept of an extracellular polymeric substance as a polymer solution, whose viscoelastic rheology is described by the classical Flory-Huggins theory. We show that one-dimensional solutions exist, which take the form at large times of travelling waves, and we characterize their form and speed in terms of the describing parameters of the problem. Numerical solutions of the time-dependent problem converge to the travelling wave solutions. [ABSTRACT FROM AUTHOR]

Abstract: We consider the penetration of liquid into a porous medium which is, in part, impervious to the liquid. Corners arise at points where the medium becomes impervious and numerical difficulties can arise due to a conflict in the boundary conditions. In addition there are start-up problems associated with numerically modelling the flow problem because the flow domain is initially empty. We consider some techniques for exploiting similarity techniques to develop accurate numerical algorithms which avoid the above difficulties. [Copyright &y& Elsevier]

We provide effective and practical guidelines on the choice of the complex denominator function of the discrete derivative as well as on the choice of the nonlocal approximation of nonlinear terms in the construction of nonstandard finite difference (NSFD) schemes. Firstly, we construct nonstandard one-stage and two-stage theta methods for a general dynamical system defined by a system of autonomous ordinary differential equations. We provide a sharp condition, which captures the dynamics of the continuous model. We discuss at length how this condition is pivotal in the construction of the complex denominator function. We show that the nonstandard theta methods are elementary stable in the sense that they have exactly the same fixed-points as the continuous model and they preserve their stability, irrespective of the value of the step size. For more complex dynamical systems that are dissipative, we identify a class of nonstandard theta methods that replicate this property. We apply the first part by considering a dynamical system that models the Ebola Virus Disease (EVD). The formulation of the model involves both the fast/direct and slow/indirect transmission routes. Using the specific structure of the EVD model, we show that, apart from the guidelines in the first part, the nonlocal approximation of nonlinear terms is guided by the productive-destructive structure of the model, whereas the choice of the denominator function is based on the conservation laws and the sub-equations that are associated with the model. We construct a NSFD scheme that is dynamically consistent with respect to the properties of the continuous model such as: positivity and boundedness of solutions; local and/or global asymptotic stability of disease-free and endemic equilibrium points; dependence of the severity of the infection on self-protection measures. Throughout the paper, we provide numerical simulations that support the theory. [ABSTRACT FROM AUTHOR]