Long Time Asymptotic Behavior of Solution of Difference Scheme for a Semilinear Parabolic Equation

In this paper we prove that the solution of implicit difference scheme for a semilinear parabolic equation converges to the solution of difference scheme for the corresponding nonlinear stationary problem as $t\rightarrow\infty$. For the discrete solution of nonlinear parabolic problem, we get its long time asymptotic behavior which is similar to that of the continuous solution. For simplicity, we consider one-dimensional problem.     

Numerical simulation of a Finite Moment Log Stable model for a European call option

Compared to the classical Black-Scholes model for pricing options, the Finite Moment Log Stable (FMLS) model can more accurately capture the dynamics of the stock prices including large movements or jumps over small time steps. In this paper, the FMLS model is written as a fractional partial differential equation and we will present a new numerical scheme for solving this model. We construct an implicit numerical scheme with second order accuracy for the FMLS and consider the stability and convergence of the scheme. In order to reduce the storage space and computational cost, we use a fast bi-conjugate gradient stabilized method (FBi-CGSTAB) to solve the discrete scheme. A numerical example is presented to show the efficiency of the numerical method and to demonstrate the order of convergence of the implicit numerical scheme. Finally, as an application, we use the above numerical technique to price a European call option. Furthermore, by comparing the FMLS model with the classical B-S model, the characteristics of the FMLS model are also analyzed.     

The solution of fractional order epidemic model by implicit Adams methods

We consider the numerical solution of the fractional order epidemic model on long time intervals of a non-fatal disease in a population. Under real-life initial conditions the problem needs to be treated by means of an implicit numerical scheme. Here we consider the use of implicit fractional linear multistep methods of Adams type. Numerical results are presented.     

A Fourier method and an extrapolation technique for Stokes’ first problem for a heated generalized second grade fluid with fractional derivative

In this article, we consider Stokes’ first problem for a heated generalized second grade fluid with fractional derivative (SFP-HGSGF). Implicit and explicit numerical approximation schemes for the SFP-HGSGF are presented. The stability and convergence of the numerical schemes are discussed using a Fourier method. In addition, the solvability of the implicit numerical approximation scheme is also analyzed. A Richardson extrapolation technique for improving the order of convergence of the implicit scheme is proposed. Finally, a numerical test is given. The numerical results demonstrate the good performance of our theoretical analysis.     

A numerical approximation of parabolic stochastic partial differential equations driven by a Poisson random measure

The paper deals with SPDEs driven by Poisson random measures in Banach spaces and its numerical approximation. We investigate the accuracy of space and time approximation. As the space approximation we consider spectral methods and as time approximation the implicit Euler scheme and the explicit Euler scheme. AMS subject classification (2000) 60H15, 35R30     

Two-dimensional inverse quasilinear parabolic problem with periodic boundary condition

In this study, we consider a coefficient problem of a quasi-linear two-dimensional parabolic inverse problem with periodic boundary and integral over determination conditions. We prove the existence, uniqueness and continuously dependence upon the data of the solution by iteration method. Also, we consider numerical solution for this inverse problem by using linearization and the implicit finite-difference scheme.     

An Analysis of the Blended Three-Step Backward Differentiation Formula Time-Stepping Scheme for the Navier-Stokes-Type System Related to Soret Convection

In this article, we investigate the stability and convergence of a new class of blended three-step Backward Differentiation Formula (BDF) time-stepping scheme for spatially discretized Navier-Stokes-type system modeling Soret driven convective flows. A Galerkin mixed finite element spatial discretization is assumed, and the temporal discretization is by the implicit blended three-step BDF scheme. The blended BDF scheme is more accurate than the classical second order accurate two-step BDF (BDF2) scheme, yet strongly A-stable. We consider an implicit, linearly extrapolated version of the scheme to improve its efficiency. We present optimal finite element error estimates and prove the scheme is unconditionally stable and convergent. Numerical experiments are presented that compare the scheme to the classical BDF2 scheme.     

Evolution of phase transitions in methane hydrate

We consider a simplified model of methane hydrates which we cast as a nonlinear evolution problem. For its well-posedness we extend the existing theory to cover the case in which the problem involves a measurable family of graphs. We represent the nonlinearity as a subgradient and prove a useful comparison principle, thus optimal regularity results follow. For the numerical solution we apply a fully implicit scheme without regularization and use the semismooth Newton algorithm for a solver, and the graph is realized as a complementarity constraint (CC). The algorithm is very robust and we extend it to define an easy and superlinearly convergent fully implicit scheme for the Stefan problem and other multivalued examples.     

Implicit iteration scheme with perturbed mapping for common fixed points of a finite family of nonexpansive mappings

In this paper we propose a new implicit iteration scheme with perturbed mapping for approximation of common fixed points of a finite family of nonexpansive mappings. We establish some convergence theorems for this implicit iteration scheme. In particular, necessary and sufficient conditions for strong convergence of this implicit iteration scheme were obtained.     

Temperature fields in a semiconductor film heated by a laser beam

We consider the heating of a thin semiconductor film on a dielectric substrate by a laser beam. The problem is solved using an implicit symmetric difference scheme. The results of a numerical experiment using real data are presented.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 60, pp. 53–59, 1986.     

An implicit finite difference scheme for the nonlinear size-structured population model

In this paper we consider the "fully nonlinear" size structured population model. We develop an implicit finite difference scheme to approximate the solution of this nonlinear partial differential equation. The convergence of this approximation to a unique bounded variation solution of this model is obtained. Numerical results to an example problem are presented.     

Numerical analysis of semilinear stochastic evolution equations in Banach spaces

The solution of stochastic evolution equations generally relies on numerical computation. Here, usually the main idea is to discretize the SPDE spatially obtaining a system of SDEs that can be solved by e.g., the Euler scheme. In this paper, we investigate the discretization error of semilinear stochastic evolution equations in L^p-spaces, resp. Banach spaces. The space discretization may be done by Galerkin approximation, for the time discretization we consider the implicit Euler, the explicit Euler scheme and the Crank–Nicholson scheme. In the last section, we give some examples, i.e., we consider an SPDEs driven by nuclear Wiener noise approximated by wavelets and delay equation approximated by finite differences.     

单物种人口模型指数隐式Euler方法的振动性

本文研究了用以描述单物种人口模型的延迟Logistic方程的数值振动性.对方程应用隐式Euler方法进行求解,针对离散格式定义了指数隐式Euler方法,证明了该方法的收敛阶为1.根据线性振动性理论获得了数值解振动的充分条件.进而还对非振动数值解的性质作了讨论.最后用数值算例对理论结果进行了验证.     

On the time discrete approximation of the Brinkman–Forchheimer equations

In this work, we study the structural stability of the fully implicit Euler scheme for the Brinkman‐Forchheimer equations. More precisely, we consider the time discretization scheme of the unsteady Brinkman–Forchheimer equations, and we prove the existence of solutions. Moreover, we derive somebia priori estimates of the discrete in time solutions. Next, with the aid of the discrete Gronwall lemma, we show that the numerical solutions depend continuously on the Brinkman and the Forchheimer coefficient. Copyright © 2011 John Wiley & Sons, Ltd.     

A uniformly convergent scheme for a system of reaction–diffusion equations

In this work a system of two parabolic singularly perturbed equations of reaction–diffusion type is considered. The asymptotic behaviour of the solution and its partial derivatives is given. A decomposition of the solution in its regular and singular parts has been used for the asymptotic analysis of the spatial derivatives. To approximate the solution we consider the implicit Euler method for time stepping and the central difference scheme for spatial discretization on a special piecewise uniform Shishkin mesh. We prove that this scheme is uniformly convergent, with respect to the diffusion parameters, having first-order convergence in time and almost second-order convergence in space, in the discrete maximum norm. Numerical experiments illustrate the order of convergence proved theoretically.     

Enumeration Approach for Linear Complementarity Problems Based on a Reformulation-Linearization Technique

In this paper, we consider the linear complementarity problem (LCP) and present a global optimization algorithm based on an application of the reformulation-linearization technique (RLT). The matrix M associated with the LCP is not assumed to possess any special structure. In this approach, the LCP is formulated first as a mixed-integer 0–1 bilinear programming problem. The RLT scheme is then used to derive a new equivalent mixed-integer linear programming formulation of the LCP. An implicit enumeration scheme is developed that uses Lagrangian relaxation, strongest surrogate and strengthened cutting planes, and a heuristic, designed to exploit the strength of the resulting linearization. Computational experience on various test problems is presented.     

Mean curvature flow with obstacles

We consider the evolution of fronts by mean curvature in the presence of obstacles. We construct a weak solution to the flow by means of a variational method, corresponding to an implicit time-discretization scheme. Assuming the regularity of the obstacles, in the two-dimensional case we show existence and uniqueness of a regular solution before the onset of singularities. Finally, we discuss an application of this result to the positive mean curvature flow.     

Approximation for Semilinear Stochastic Evolution Equations

We investigate the approximation by space and time discretization of quasi linear evolution equations driven by nuclear or space time white noise. An error bound for the implicit Euler, the explicit Euler, and the Crank–Nicholson scheme is given and the stability of the schemes are considered. Lastly we give some examples of different space approximation, i.e., we consider approximation by eigenfunction, finite differences and wavelets.     

On the order optimality of the regularization via inexact Newton iterations

Inexact Newton regularization methods have been proposed by Hanke and Rieder for solving nonlinear ill-posed inverse problems. Every such a method consists of two components: an outer Newton iteration and an inner scheme providing increments by regularizing local linearized equations. The method is terminated by a discrepancy principle. In this paper we consider the inexact Newton regularization methods with the inner scheme defined by Landweber iteration, the implicit iteration, the asymptotic regularization and Tikhonov regularization. Under certain conditions we obtain the order optimal convergence rate result which improves the suboptimal one of Rieder. We in fact obtain a more general order optimality result by considering these inexact Newton methods in Hilbert scales.     

Error estimates for a mixed finite element discretization of some degenerate parabolic equations

We consider a numerical scheme for a class of degenerate parabolic equations, including both slow and fast diffusion cases. A particular example in this sense is the Richards equation modeling the flow in porous media. The numerical scheme is based on the mixed finite element method (MFEM) in space, and is of one step implicit in time. The lowest order Raviart–Thomas elements are used. Here we extend the results in Radu et al. (SIAM J Numer Anal 42:1452–1478, 2004), Schneid et al. (Numer Math 98:353–370, 2004) to a more general framework, by allowing for both types of degeneracies. We derive error estimates in terms of the discretization parameters and show the convergence of the scheme. The features of the MFEM, especially of the lowest order Raviart–Thomas elements, are now fully exploited in the proof of the convergence. The paper is concluded by numerical examples.