Bifurcation analysis of a population model and the resulting SIS epidemic model with delay

This paper deals with the model for matured population growth proposed in Cooke et al. [Interaction of matiration delay and nonlinear birth in population and epidemic models, J. Math. Biol. 39 (1999) 332–352] and the resulting SIS epidemic model. The dynamics of these two models are still largely undetermined, and in this paper, we perform some bifurcation analysis to the models. By applying the global bifurcation theory for functional differential equations, we are able to show that the population model allows multiple periodic solutions. For the SIS model, we obtain some local bifurcation results and derive formulas for determining the bifurcation direction and the stability of the bifurcated periodic solution.     

A Series Solution of the Unsteady Von Kármán Swirling Viscous Flows

A new analytic technique is applied to solve the unsteady viscous flow due to an infinite rotating disk, governed by a set of two fully coupled nonlinear partial differential equations deduced directly from the exact Navier-Stokes equations. The system of coupled nonlinear partial differential equations is replaced by a sequence of uncoupled systems of linear ordinary differential equations. Different from all other previous analytic results, our series solution is accurate and valid for all time in the whole spatial region. Accurate expressions for skin friction coefficients are given, which are valid for all time. Such kind of series solutions have not been reported, to the best of our knowledge.      

Numerical approach to differential equations of fractional order

In this paper, the variational iteration method and the Adomian decomposition method are implemented to give approximate solutions for linear and nonlinear systems of differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper presents a numerical comparison between the two methods for solving systems of fractional differential equations. Numerical results show that the two approaches are easy to implement and accurate when applied to differential equations of fractional order.     

Exact analytic solutions of Schrödinger linear and nonlinear equations

Exact analytic solutions of Schrödinger linear partial differential equations are obtained. Moreover, the cubic nonlinear Schrödinger equation is treated with the use of a well-known functional analytic method and the existence of convergent power series solutions is proved. From these solutions, under certain initial conditions, similar results as those presented in the literature are obtained.     

一类具有非线性发生率的SIR传染病模型的全局稳定性

研究一类具有非线性发生率的SIR传染病模型.应用微分方程定性理论分别得到了该系统无病平衡点、地方病平衡点全局渐近稳定的充分条件,并进行了数值模拟.     

Asymptotic behaviors of stochastic epidemic models with jump-diffusion

In this paper, we classify the asymptotic behavior for a class of stochastic SIR epidemic models represented by stochastic differential systems where the Brownian motions and Lévy jumps perturb to the linear terms of each equation. We construct a threshold value between permanence and extinction and develop the ergodicity of the underlying system. It is shown that the transition probabilities converge in total variation norm to the invariant measure. Our results can be considered as a significant contribution in studying the long term behavior of stochastic differential models because there are no restrictions imposed to the parameters of models. Techniques used in proving results of this paper are new and suitable to deal with other stochastic models in biology where the noises may perturb to nonlinear terms of equations or with delay equations.     

An analytical solution for a nonlinear time-delay model in biology

In this paper, the homotopy analysis method is applied to develop a analytic approach for nonlinear differential equations with time-delay. A nonlinear model in biology is used as an example to show the basic ideas of this analytic approach. Different from other analytic techniques, the homotopy analysis method provides a simple way to ensure the convergence of the solution series, so that one can always get accurate approximations. A new discontinuous function is defined so as to express the piecewise continuous solutions of time-delay differential equations in a way convenient for symbolic computations. It is found that the time-delay has a great influence on the solution of the time-delay nonlinear differential equation. This approach has general meanings and can be applied to solve other nonlinear problems with time-delay.     

Laplace transform for the solution of higher order deformation equations arising in the homotopy analysis method

A new analytic approach for solving nonlinear ordinary differential equations with initial conditions is proposed. First, the homotopy analysis method is used to transform a nonlinear differential equation into a system of linear differential equations; then, the Laplace transform method is applied to solve the resulting linear initial value problems; finally, the solutions to the linear initial value problems are employed to form a convergent series solution to the given problem. The main advantage of the new approach is that it provides an effective way to solve the higher order deformation equations arising in the homotopy analysis method.     

Explicit series solutions of some linear and nonlinear Schrodinger equations via the homotopy analysis method

In this paper, by means of the homotopy analysis method (HAM), the solutions of some Schrodinger equations are exactly obtained in the form of convergent Taylor series. The HAM contains the auxiliary parameter ?, that provides a convenient way of controlling the convergent region of series solutions. This analytical method is employed to solve linear and nonlinear examples to obtain the exact solutions. HAM is a powerful and easy-to-use analytic tool for nonlinear problems.     

New explicit approximate solution of MHD viscoelastic boundary layer flow over stretching sheet

In this paper, the study the momentum and heat transfer characteristics in an incompressible electrically conducting non‐Newtonian boundary layer flow of a viscoelastic fluid over a stretching sheet. The partial differential equations governing the flow and heat transfer characteristics are converted into highly nonlinear coupled ordinary differential equations by similarity transformations. The resultant coupled highly nonlinear ordinary differential equations are solved by means of, homotopy analysis method (HAM) for constructing an approximate solution of heat transfer in magnetohydrodynamic (MHD) viscoelastic boundary layer flow over a stretching sheet with non‐uniform heat source. The proposed method is a strong and easy to use analytic tool for nonlinear problems and does not need small parameters in the equations. The HAM solutions contain an auxiry parameter, which provides a convenient way of controlling the convergence region of series solutions. The results obtained here reveal that the proposed method is very effective and simple for solving nonlinear evolution equations. The method is straightforward and concise, and it can also be applied to other nonlinear evolution equations in physics. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison between variational iteration method and homotopy perturbation method for linear and nonlinear partial differential equations with the nonhomogeneous initial conditions

In this study, linear and nonlinear partial differential equations with the nonhomogeneous initial conditions are considered. We used Variational iteration method (VIM) and Homotopy perturbation method (HPM) for solving these equations. Both methods are used to obtain analytic solutions for different types of differential equations. Four examples are presented to show the application of the present techniques. In these schemes, the solution takes the form of a convergent series with easily computable components. The present methods perform extremely well in terms of efficiency and simplicity. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A general approach to get series solution of non-similarity boundary-layer flows

An analytic method for strongly non-linear problems, namely the homotopy analysis method (HAM), is applied to give convergent series solution of non-similarity boundary-layer flows. As an example, the non-similarity boundary-layer flows over a stretching flat sheet are used to show the validity of this general analytic approach. Without any assumptions of small/large quantities, the corresponding non-linear partial differential equation with variable coefficients is transferred into an infinite number of linear ordinary differential equations with constant coefficients. More importantly, an auxiliary artificial parameter is used to ensure the convergence of the series solution. Different from previous analytic results, our series solutions are convergent and valid for all physical variables in the whole domain of flows. This work illustrates that, by means of the homotopy analysis method, the non-similarity boundary-layer flows can be solved in a similar way like similarity boundary-layer flows. Mathematically, this analytic approach is rather general in principle and can be applied to solve different types of non-linear partial differential equations with variable coefficients in science and engineering.     

Solving nonlinear fractional partial differential equations using the homotopy analysis method

In this article, the homotopy analysis method is applied to solve nonlinear fractional partial differential equations. On the basis of the homotopy analysis method, a scheme is developed to obtain the approximate solution of the fractional KdV, K(2,2), Burgers, BBM‐Burgers, cubic Boussinesq, coupled KdV, and Boussinesq‐like B(m,n) equations with initial conditions, which are introduced by replacing some integer‐order time derivatives by fractional derivatives. The homotopy analysis method for partial differential equations of integer‐order is directly extended to derive explicit and numerical solutions of the fractional partial differential equations. The solutions of the studied models are calculated in the form of convergent series with easily computable components. The results of applying this procedure to the studied cases show the high accuracy and efficiency of the new technique. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Solving systems of fractional differential equations using differential transform method

This paper presents approximate analytical solutions for systems of fractional differential equations using the differential transform method. The fractional derivatives are described in the Caputo sense. The application of differential transform method, developed for differential equations of integer order, is extended to derive approximate analytical solutions of systems of fractional differential equations. The solutions of our model equations are calculated in the form of convergent series with easily computable components. Some examples are solved as illustrations, using symbolic computation. The numerical results show that the approach is easy to implement and accurate when applied to systems of fractional differential equations. The method introduces a promising tool for solving many linear and nonlinear fractional differential equations.     

Melnikov analysis of chaos in a general epidemiological model

The purpose of this paper is to study a SIR model of epidemic dynamics with a periodically modulated nonlinear incidence rate. We must go, for the first time, through a series of coordinate transformations to bring the equations into amenable to Melnikov analysis. This analysis establishes mathematically the existence of chaotic motion of the models by Melnikov's method. The numerical simulations are made for the conclusions in this paper.     

加权无标度网络中的疾病传播

提出具有加权传播率和非线性传染能力的SIR模型和SIS模型,通过平均场方法证明了这两个模型在加权无标度网络中可以存在非零的传播阈值,从而传播率需要跨越更大的传播阈值才能流行.并且得到的结果在特殊情况下可退化为已有的一些经典结论.     

Analytic approximate solution of three‐dimensional Navier‐Stokes equations of flow between two stretchable disks

The similarity transform for the steady three‐dimensional Navier‐Stokes equations of flow between two stretchable disks gives a system of nonlinear ordinary differential equations which is analytically solved by applying a newly developed method, namely, the homotopy analysis method. The analytic solutions of the system of nonlinear ordinary differential equations are constructed in the series form. The convergence of the obtained series solutions is analyzed. The validity of our solutions is verified by the numerical results. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

On the homotopy analysis method for backward/forward-backward stochastic differential equations

In this paper, an analytic approximation method for highly nonlinear equations, namely the homotopy analysis method (HAM), is employed to solve some backward stochastic differential equations (BSDEs) and forward-backward stochastic differential equations (FBSDEs), including one with high dimensionality (up to 12 dimensions). By means of the HAM, convergent series solutions can be quickly obtained with high accuracy for a FBSDE in a 6-dimensional case, within less than 1 % CPU time used by a currently reported numerical method for the same case [34]. Especially, as dimensionality enlarges, the increase of computational complexity for the HAM is not as dramatic as this numerical method. All of these demonstrate the validity and high efficiency of the HAM for the backward/forward-backward stochastic differential equations in science, engineering, and finance.     

A class of Lyapunov functions and the global stability of some epidemic models with nonlinear incidence

In this paper, by investigating an SIR epidemic model with nonlinear incidence, we present a new technique for proving the global stability of the endemic equilibrium, which consists of introducing a variable transformation and constructing a more general Lyapunov function. For the model we obtain the following results. The disease-free equilibrium is globally stable in the feasible region as the basic reproduction number is less than or equal to unity, and the endemic equilibrium is globally stable in the feasible region as the basic reproduction number is greater than unity.The generality of the technique is illustrated by considering certain nonlinear incidences and SIS and SIRS epidemic models.     

Numerical comparison of methods for solving linear differential equations of fractional order

In this article, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving linear differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper will present a numerical comparison between the two methods and a conventional method such as the fractional difference method for solving linear differential equations of fractional order. The numerical results demonstrates that the new methods are quite accurate and readily implemented.