Lattice Approximations for Stochastic Quasi-Linear Parabolic Partial Differential Equations Driven by Space-Time White Noise I

We approximate quasi-linear parabolic SPDEs substituting the derivatives in the space variable with finite differences. When the nonlinear terms in the equation are Lipschitz continuous we estimate the rate of L^p convergence of the approximations and we also prove their almost sure uniform convergence to the solution. When the nonlinear terms are not Lipschitz continuous we obtain this convergence in probability, if the pathwise uniqueness for the equation holds.     

Newton iterations in implicit time-stepping scheme for differential linear complementarity systems

We propose a generalized Newton method for solving the system of nonlinear equations with linear complementarity constraints in the implicit or semi-implicit time-stepping scheme for differential linear complementarity systems (DLCS). We choose a specific solution from the solution set of the linear complementarity constraints to define a locally Lipschitz continuous right-hand-side function in the differential equation. Moreover, we present a simple formula to compute an element in the Clarke generalized Jacobian of the solution function. We show that the implicit or semi-implicit time-stepping scheme using the generalized Newton method can be applied to a class of DLCS including the nondegenerate matrix DLCS and hidden Z-matrix DLCS, and has a superlinear convergence rate. To illustrate our approach, we show that choosing the least-element solution from the solution set of the Z-matrix linear complementarity constraints can define a Lipschitz continuous right-hand-side function with a computable Lipschitz constant. The Lipschitz constant helps us to choose the step size of the time-stepping scheme and guarantee the convergence.     

Product integration methods for solving a system of nonlinear Volterra integral equations

In this paper the technique of subtracting out singularities is used to derive explicit and implicit product Euler schemes with order one convergence and a product trapezoidal scheme with order two convergence for a system of Volterra integral equations with a weakly singular kernel. The convergence proofs of the numerical schemes are presented; these are nonstandard since the nonlinear function involved in the integral equation system does not satisfy a global Lipschitz condition.     

Convergence of the multigrid full approximation scheme for a class of elliptic mildly nonlinear boundary value problems

Summary The multigrid full approximation scheme (FAS MG) is a well-known solver for nonlinear boundary value problems. In this paper we restrict ourselves to a class of second order elliptic mildly nonlinear problems and we give local conditions, e.g. a local Lipschitz condition on the derivative of the continuous operator, under which the FAS MG with suitably chosen parameters locally converges. We prove quantitative convergence statements and deduce explicit bounds for important quantities such as the radius of a ball of guaranteed convergence, the number of smoothings needed, the number of coarse grid corrections needed and the number of FAS MG iterations needed in a nested iteration. These bounds show well-known features of the FAS MG scheme.     

Strong convergence of split-step backward Euler method for stochastic differential equations with non-smooth drift

In this paper, we are concerned with the numerical approximation of stochastic differential equations with discontinuous/nondifferentiable drifts. We show that under one-sided Lipschitz and general growth conditions on the drift and global Lipschitz condition on the diffusion, a variant of the implicit Euler method known as the split-step backward Euler (SSBE) method converges with strong order of one half to the true solution. Our analysis relies on the framework developed in [D. J. Higham, X. Mao and A. M. Stuart, Strong convergence of Euler-type methods for nonlinear stochastic differential equations, SIAM Journal on Numerical Analysis, 40 (2002) 1041-1063] and exploits the relationship which exists between explicit and implicit Euler methods to establish the convergence rate results.     

On stability of difference schemes. Central schemes for hyperbolic conservation laws with source terms

The stability of nonlinear explicit difference schemes with not, in general, open domains of the scheme operators are studied. For the case of path-connected, bounded, and Lipschitz domains, we establish the notion that a multi-level nonlinear explicit scheme is stable iff (if and only if) the corresponding scheme in variations is stable. A new modification of the central Lax–Friedrichs (LxF) scheme is developed to be of the second-order accuracy. The modified scheme is based on nonstaggered grids. A monotone piecewise cubic interpolation is used in the central scheme to give an accurate approximation for the model in question. The stability of the modified scheme is investigated. Some versions of the modified scheme are tested on several conservation laws, and the scheme is found to be accurate and robust. As applied to hyperbolic conservation laws with, in general, stiff source terms, it is constructed a second-order nonstaggered central scheme based on operator-splitting techniques.     

The tamed Milstein method for commutative stochastic differential equations with non-globally Lipschitz continuous coefficients

For stochastic differential equations (SDEs) with a superlinearly growing and globally one-sided Lipschitz continuous drift coefficient, the classical explicit Euler scheme fails to converge strongly to the exact solution. Recently, an explicit strongly convergent numerical scheme, called the tamed Euler method, has been proposed in [8] for such SDEs. Motivated by their work, we here introduce a tamed version of the Milstein scheme for SDEs with commutative noise. The proposed method is also explicit and easily implementable, but achieves higher strong convergence order than the tamed Euler method does. In recovering the strong convergence order one of the new method, new difficulties arise and kind of a bootstrap argument is developed to overcome them. Finally, an illustrative example confirms the computational efficiency of the tamed Milstein method compared to the tamed Euler method.     

An adaptive scheme for the approximation of dissipative systems

We propose a new scheme for the long time approximation of a diffusion when the drift vector field is not globally Lipschitz. Under this assumption, a regular explicit Euler scheme–with constant or decreasing step–may explode and implicit Euler schemes are CPU-time expensive. The algorithm we introduce is explicit and we prove that any weak limit of the weighted empirical measures of this scheme is a stationary distribution of the stochastic differential equation. Several examples are presented including gradient dissipative systems and Hamiltonian dissipative systems.     

On convergence of finite volume schemes for one-dimensional two-phase flow in porous media

This paper deals with development and analysis of finite volume schemes for a one-dimensional nonlinear, degenerate, convection-diffusion equation having application in petroleum reservoir and groundwater aquifer simulation. The main difficulty is that the solution typically lacks regularity due to the degenerate nonlinear diffusion term. We analyze and compare three families of numerical schemes corresponding to explicit, semi-implicit, and implicit discretization of the diffusion term and a Godunov scheme for the advection term. L^∞ stability under appropriate CFL conditions and BV estimates are obtained. It is shown that the schemes satisfy a discrete maximum principle. Then we prove convergence of the approximate solution to the weak solution of the problem. Results of numerical experiments using the present approach are reported.     

非线性波动方程的弱隐式与显式差分方法

广泛出现于物理、化学、机械动力学、生物、几何学等领域的非线性波动方程已经有很多的研究工作,Sine-Gordon方程和非线性受迫振动方程就是典型的例子.周毓麟教授在[1]中研究了非线性波动方程组     

A memory‐efficient finite volume method for advection‐diffusion‐reaction systems with nonsmooth sources

We present a parallel matrix‐free implicit finite volume scheme for the solution of unsteady three‐dimensional advection‐diffusion‐reaction equations with smooth and Dirac‐Delta source terms. The scheme is formally second order in space and a Newton–Krylov method is employed for the appearing nonlinear systems in the implicit time integration. The matrix‐vector product required is hardcoded without any approximations, obtaining a matrix‐free method that needs little storage and is well‐suited for parallel implementation. We describe the matrix‐free implementation of the method in detail and give numerical evidence of its second‐order convergence in the presence of smooth source terms. For nonsmooth source terms, the convergence order drops to one half. Furthermore, we demonstrate the method's applicability for the long‐time simulation of calcium flow in heart cells and show its parallel scaling. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq31: 143–167, 2015     

Second-order characteristic schemes in time and age for a nonlinear age-structured population model

In this paper, we analyze two new second-order characteristic schemes in time and age for an age-structured population model with nonlinear diffusion and reaction. By using the characteristic difference to approximate the transport term and the average along the characteristics to treat the nonlinear spatial diffusion and reaction terms, an implicit second-order characteristic scheme is proposed. To compute the nonlinear approximation system, an explicit second-order characteristic scheme in time and age is further proposed by using the extrapolation technique. The global existence and uniqueness of the solution of the nonlinear approximation scheme are established by using the theory of variation methods, Schauder’s fixed point theorem, and the technique of prior estimates. The optimal error estimates of second order in time and age are strictly proved for both the implicit and the explicit characteristic schemes. Numerical examples are given to illustrate the performance of the methods.     

The Secant method for nondifferentiable operators

In this paper, we use the Secant method to find a solution of a nonlinear operator equation in Banach spaces. A semilocal convergence result is obtained. For that, we consider a condition for divided differences which generalizes the usual ones, i.e., Lipschitz continuous or Hölder continuous conditions. Besides, we apply our results to approximate the solution of a nonlinear equation.     

A mixed finite element method with Lagrange multipliers for nonlinear exterior transmission problems

Summary. We apply a mixed finite element method to numerically solve a class of nonlinear exterior transmission problems in R ² with inhomogeneous interface conditions. Besides the usual unknowns required for the dual-mixed method, which include the gradient of the temperature in this nonlinear case, our approach makes use of the trace of the outer solution on the transmission boundary as a suitable Lagrange multiplier. In addition, we use a boundary integral operator to reduce the original transmission problem on the unbounded region into a nonlocal one on a bounded domain. In this way, we are lead to a two-fold saddle point operator equation as the resulting variational formulation. We prove that the continuous formulation and the associated Galerkin scheme defined with Raviart-Thomas spaces are well posed, and derive the a-priori estimates and the corresponding rate of convergence. Then, we introduce suitable local problems and deduce first an implicit reliable and quasi-efficient a-posteriori error estimate, and then a fully explicit reliable one. Finally, several numerical results illustrate the effectivity of the explicit estimate for the adaptive computation of the discrete solutions. Mathematics Subject Classification (2000): 65N30, 65N38, 65N22, 65F10This research was partially supported by CONICYT-Chile through the FONDAP Program in Applied Mathematics, and by the Dirección de Investigación of the Universidad de Concepción through the Advanced Research Groups Program.     

Strong convergence rates for backward Euler on a class of nonlinear jump-diffusion problems

We generalise the current theory of optimal strong convergence rates for implicit Euler-based methods by allowing for Poisson-driven jumps in a stochastic differential equation (SDE). More precisely, we show that under one-sided Lipschitz and polynomial growth conditions on the drift coefficient and global Lipschitz conditions on the diffusion and jump coefficients, three variants of backward Euler converge with strong order of one half. The analysis exploits a relation between the backward and explicit Euler methods.     

SINGLE PROJECTION ALGORITHM FOR VARIATIONAL INEQUALITIES IN BANACH SPACES WITH APPLICATION TO CONTACT PROBLEM

We study the single projection algorithm of Tseng for solving a variational inequality problem in a 2-uniformly convex Banach space. The underline cost function of the variational inequality is assumed to be monotone and Lipschitz continuous. A weak convergence result is obtained under reasonable assumptions on the variable step-sizes. We also give the strong convergence result for when the underline cost function is strongly monotone and Lipchitz continuous. For this strong convergence case, the proposed method does not require prior knowledge of the modulus of strong monotonicity and the Lipschitz constant of the cost function as input parameters, rather, the variable step-sizes are diminishing and non-summable. The asymptotic estimate of the convergence rate for the strong convergence case is also given. For completeness, we give another strong convergence result using the idea of Halpern's iteration when the cost function is monotone and Lipschitz continuous and the variable step-sizes are bounded by the inverse of the Lipschitz constant of the cost function.Finally, we give an example of a contact problem where our proposed method can be applied.     

Strong Convergence of a Fully Discrete Scheme for Multiplicative Noise Driving SPDEs with Non-Globally Lipschitz Continuous Coefficients

This work investigates strong convergence of numerical schemes for nonlinear multiplicative noise driving stochastic partial differential equations under some weaker conditions imposed on the coefficients avoiding the commonly used global Lipschitz assumption in the literature. Space-time fully discrete scheme is proposed, which is performed by the finite element method in space and the implicit Euler method in time. Based on some technical lemmas including regularity properties for the exact solution of the considered problem, strong convergence analysis with sharp convergence rates for the proposed fully discrete scheme is rigorously established.     

Implicit Scheme for Stochastic Parabolic Partial Diferential Equations Driven by Space-Time White Noise

We extend Rothe's method of solving linear parabolic PDEs to the case of nonlinear SPDEs driven by space-time white noise. When the nonlinear terms are Lipschitz functions we prove almost sure convergence of the approximations uniformly in time and space. When the nonlinear drift term is only measurable we obtain the convergence in probability, by using Malliavin calculus.     

A modification of the Chebyshev method

In this paper we use a one-parametric family of second-orderiterations to solve a nonlinear operator equation in a Banachspace. Two different analyses of convergence are shown. First,under standard Newton-Kantorovich conditions, we establish aKantorovich-type convergence theorem. Second, another Kantorovich-typeconvergence theorem is proved, when the first Frchet-derivativeof the operator satisfies a Lipschitz condition. We also givean explicit expression for the error bound of the family ofmethods in terms of a real parameter 0.     

Convergence of the Tamed Euler Method for Stochastic Differential Equations with Piecewise Continuous Arguments Under Non-global Lipschitz Continuous Coefficients

In this paper, we deal with the strong convergence of numerical methods for stochastic differential equations with piecewise continuous arguments (SEPCAs) with at most polynomially growing drift coefficients and global Lipschitz continuous diffusion coefficients. An explicit and time-saving tamed Euler method is used to solve this type of SEPCAs. We show that the tamed Euler method is bounded in pth moment. And then the convergence of the tamed Euler method is proved. Moreover, the convergence order is one-half. Several numerical simulations are shown to verify the convergence of this method.