Dynamical behaviors of a discrete HIV‐1 virus model with bilinear infective rate

We establish a discrete virus dynamic model by discretizing a continuous HIV‐1 virus model with bilinear infective rate using ‘hybrid’ Euler method. We discuss not only the existence and global stability of the uninfected equilibrium but also the existence and local stability of the infected equilibrium. We prove that there exists a crucial value similar to that of the continuous HIV‐1 virus dynamics, which is called the basic reproductive ratio of the virus. If the basic reproductive ratio of the virus is less than one, the uninfected equilibrium is globally asymptotically stable. If the basic reproductive ratio of the virus is larger than one, the infected equilibrium exists and is locally stable. Moreover, we consider the permanence for such a system by constructing a Lyapunov function v_n. Copyright © 2013 John Wiley & Sons, Ltd.     

A predator-prey model with Beddington-DeAngelis functional response: a non-standard finite-difference method

In this paper, we transform a continuous-time predator-prey model with Beddington–DeAngelis functional response into a discrete-time model by nonstandard finite difference scheme (NSFD). The NSFD model shows complete dynamic consistency with its continuous counterpart for any step size. However, the discrete model of same continuous system obtained by Euler forward method shows dynamic inconsistency for larger step size. Extensive numerical simulations have been done to compare the dynamics of NSFD system and Euler system. Our analysis reveals that dynamics of NSFD model is independent of the step-size, whereas the dynamics of the standard discrete model completely depends on the step-size and produces spurious dynamics like chaos.     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.     

Permanence of a discrete SIRS epidemic model with time delays

In this paper, we consider the permanence of a discrete SIRS epidemic model with time delays. This model is constructed from the discretization by the Euler method. Applying the technique to prove the existence of an eventual lower bound in a continuous epidemic model, a sufficient condition for the permanence of the above discrete SIRS epidemic model is obtained.     

Numerical analysis for generalized Forchheimer flows of slightly compressible fluids

In this article, we will consider the generalized Forchheimer flows for slightly compressible fluids. Using Muskat's and Ward's general form of Forchheimer equations, we describe the fluid dynamics by a nonlinear degenerate parabolic equation for density. The long‐time numerical approximation of the nonlinear degenerate parabolic equation with time dependent boundary conditions is studied. The stability for all time is established in a continuous time scheme and a discrete backward Euler scheme. A Gronwall's inequality‐type is used to study the asymptotic behavior of the solution. Error estimates for the solution are derived for both continuous and discrete time procedures. Numerical experiments confirm the theoretical analysis regarding convergence rates.     

Global stability for a discrete epidemic model for disease with immunity and latency spreading in a heterogeneous host population

In this paper, we propose a discrete epidemic model for disease with immunity and latency spreading in a heterogeneous host population, which is derived from the continuous case by using the well-known backward Euler method and by applying a Lyapunov function technique, which is a discrete version of that in the paper by Prüss et al. [J. Prüss, L. Pujo-Menjouet, G.F. Webb, R. Zacher, Analysis of a model for the dynamics of prions, Discrete Contin. Dyn. Syst. Ser. B 6 (2006) 225-235]. It is shown that the global dynamics of this discrete epidemic model with latency are fully determined by a single threshold parameter.     

On Kirchhoff's Model of Parabolic Type

In this article, the existence of a global strong solution for all finite time is derived for the Kirchhoff's model of parabolic type. Based on exponential weight function, some new regularity results which reflect the exponential decay property are obtained for the exact solution. For the related dynamics, the existence of a global attractor is shown to hold for the problem when the non-homogeneous forcing function is either independent of time or in L^∞(L²). With the finite element Galerkin method applied in spatial direction keeping time variable continuous, a semidiscrete scheme is analyzed, and it is also established that the semidiscrete system has a global discrete attractor. Optimal error estimates in L^∞(H¹) norm are derived which are valid uniformly in time. Further, based on a backward Euler method, a completely discrete scheme is analyzed and error estimates are derived. It is also further, observed that in cases where f = 0 or f = O(e^?γ₀t) with γ₀ > 0, the discrete solutions and error estimates decay exponentially in time. Finally, some numerical experiments are discussed which confirm our theoretical findings.     

Adaptive techniques for Landau–Lifshitz–Gilbert equation with magnetostriction

In this paper we propose a time–space adaptive method for micromagnetic problems with magnetostriction. The considered model consists of coupled Maxwell's, Landau–Lifshitz–Gilbert (LLG) and elastodynamic equations. The time discretization of Maxwell's equations and the elastodynamic equation is done by backward Euler method, the space discretization is based on Whitney edge elements and linear finite elements, respectively. The fully discrete LLG equation reduces to an ordinary differential equation, which is solved by an explicit method, that conserves the norm of the magnetization.     

Finite element approximation of nonlocal parabolic problem

In this article, we present a finite element scheme combined with backward Euler method to solve a nonlocal parabolic problem. An important issue in the numerical solution of nonlocal problems while using Newton's method is related to its structure. In fact differently from the local case where the Jacobian matrix is sparse and banded, in the nonlocal case the Jacobian matrix is dense and computations are much more onerous compared to that for differential equations. In order to avoid this difficulty, we use the technique given by Gudi (SIAM J Numer Anal 50 (2012), 657–668) for elliptic nonlocal problem of Kirchhoff type. We discuss the well‐posedness of the weak formulation at continuous as well as at discrete levels. We also derive a priori error estimates for semidiscrete and fully discrete formulations in L² and H¹ norms. Results based on the usual finite element method are provided to confirm the theoretical estimates. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 786–813, 2017     

Recursive Marginal Quantization of the Euler Scheme of a Diffusion Process

Abstract

We propose a new approach to quantize the marginals of the discrete Euler diffusion process. The method is built recursively and involves the conditional distribution of the marginals of the discrete Euler process. Analytically, the method raises several questions like the analysis of the induced quadratic quantization error between the marginals of the Euler process and the proposed quantizations. We show in particular that at every discretization step t_k of the Euler scheme, this error is bounded by the cumulative quantization errors induced by the Euler operator, from times t₀ = 0 to time t_k. For numerics, we restrict our analysis to the one-dimensional setting and show how to compute the optimal grids using a Newton–Raphson algorithm. We then propose a closed formula for the companion weights and the transition probabilities associated to the proposed quantizations. This allows us to quantize in particular diffusion processes in local volatility models by reducing dramatically the computational complexity of the search of optimal quantizers while increasing their computational precision with respect to the algorithms commonly proposed in this framework. Numerical tests are carried out for the Brownian motion and for the pricing of European options in a local volatility model. A comparison with the Monte Carlo simulations shows that the proposed method may sometimes be more efficient (w.r.t. both computational precision and time complexity) than the Monte Carlo method.     

An unconditionally stable uncoupled scheme for a triphasic Cahn–Hilliard/Navier–Stokes model

We propose an original scheme for the time discretization of a triphasic Cahn–Hilliard/Navier–Stokes model. This scheme allows an uncoupled resolution of the discrete Cahn–Hilliard and Navier‐Stokes system, which is unconditionally stable and preserves, at the discrete level, the main properties of the continuous model. The existence of discrete solutions is proved, and a convergence study is performed in the case where the densities of the three phases are the same. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq. 2013     

Finite element approximation of a forward and backward anisotropic diffusion model in image denoising and form generalization

A new forward–backward anisotropic diffusion model is introduced. The two limit cases are the Perona‐Malik equation and the Total Variation flow model. A fully discrete finite element scheme is studied using C⁰‐piecewise linear elements in space and the backward Euler difference scheme in time. A priori estimates are proven. Numerical results in image denoising and form generalization are presented.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Fully discrete T‐ψ finite element method to solve a nonlinear induction hardening problem

We study an induction hardening model described by Maxwell's equations coupled with a heat equation. The magnetic induction field is assumed a nonlinear constitutional relation and the electric conductivity is temperature‐dependent. The T ‐ψ method is to transform Maxwell's equations to the vector–scalar potential formulations and to solve the potentials by means of the finite element method. In this article, we present a fully discrete T ‐ψ finite element scheme for this nonlinear coupled problem and discuss its solvability. We prove that the discrete solution converges to a weak solution of the continuous problem. Finally, we conclude with two numerical experiments for the coupled system.     

A note on the convergence of discretized dynamic iteration

In this note we prove convergence results, including error estimates, for the dynamic iteration scheme where the forward Euler and backward Euler method are used to compute the iterates. The proofs are interesting in that they are exact analogues of the proof for the continuous case, using discrete versions of Gronwall's inequality.     

Mean‐square stability of the backward Euler–Maruyama method for neutral stochastic delay differential equations with jumps

This paper is mainly considered whether the mean‐square stability of neutral stochastic delay differential equations (NSDDEs) with jumps is shared with that of the backward Euler–Maruyama method. Under the one‐sided Lipschitz condition and the linear growth condition, the trivial solution of NSDDEs with jumps is proved to be mean‐square stable by using the functional comparison principle and the Barbalat's lemma. It is shown that the backward Euler–Maruyama method can reproduce the mean‐square stability of the trivial solution under the same conditions. The implicit backward Euler–Maruyama method shows better characteristic than the explicit Euler–Maruyama method for the reason that it works without the linear growth condition on the drift coefficient. Compared with some existing results, our results do not need to add extra condition on the neutral part. The conclusions can be applied to NSDDEs and SDDEs with jumps. The effectiveness of the theoretical results is illustrated by an example. Copyright © 2016 John Wiley & Sons, Ltd.     

Diffusion approximation for nonparametric autoregression

A nonparametric statistical model of small diffusion type is compared with its discretization by a stochastic Euler difference scheme. It is shown that the discrete and continuous models are asymptotically equivalent in the sense of Le Cam's deficiency distance for statistical experiments, when the discretization step decreases with the noise intensity ε. Received: 12 April 1996 / Revised version: 29 October 1997     

Convergence of the Point Vortex Methods for Euler Equation on Half Plane

1.IntroductionThevortexmethodsareefficientnumericalmethodofsimulatingincompressibleflowathighReynold'snumber-TheconvergenceofthevortexmethodsfortheinitialvalueproblemsofEulerequationwasfirstobtainedbyHald[41,thentheresultswereimprovedbyseveralauthors[1'2'3'5].Butinfact,manypracticalproblemsareconsideredinaboundeddomainoranexteriordomain,andthenumericalboundaryconditionhasanimportanteffectonnumericalresult.Theparticletrajectoriesofexactsolutionwillnotgooutfromthedomain,butitisnotthecaseinprat…     

Multiple grid and Osher's scheme for the efficient solution of the steady Euler equations

An iterative method is developed for the solution of the steady Euler equations for inviscid flow. The system of hyperbolic conservation laws is discretized by a finite-volume Osher-discretization. The iterative method is a multiple grid (FAS) iteration with symmetric Gauss-Seidel (SGS) as a relaxation method. Initial estimates are obtained by full multigrid (FMG). In the pointwise relaxation the equations are kept in block-coupled form and local linearization of the equations and the boundary conditions is considered. The efficient formulation of Osher's discretization of the 2-D non-isentropic steady Euler equations and its linearization is presented. The efficiency of FAS-SGS iteration is shown for a transonic model problem. It appears that, for the problem considered, the rate of convergence is almost independent of the gridsize and that for all meshsizes the discrete system is solved up to truncation error accuracy in only a few (2 or 3) iteration cycles.     

Dynamical behaviors of a discrete two-dimensional competitive system exactly driven by the large centre

In this paper, a new discrete large-sub-center system is obtained by using the Euler and nonstandard discretization methods for the corresponding continuous system. It is surprised that all dynamic behaviors of the discrete system are exactly driven by the large-center equation, for example, the stabilities, the bifurcations, the period-doubling orbits, and the chaotic dynamics, etc. Additionally, the global asymptotical stability, the existence of exact 2-periodic solutions, the flip bifurcation theorem, and the invariant set of the sub-center equation is also given. These results reveal far richer dynamics of the discrete model compared with the continuous model. Through numerical simulation, we can observe some complex dynamic behaviors, such as period-doubling cascade, periodic windows, chaotic dynamics, etc. Especially, our theoretical results are also showed by those numerical simulations.