Numerical solution of the Dirichlet problem for nonlinear parabolic equations by a probabilistic approach

A number of new layer methods for solving the Dirichlet problemfor semilinear parabolic equations are constructed by usingprobabilistic representations of their solutions. The methodsexploit the ideas of weak sense numerical integration of stochasticdifferential equations in a bounded domain. Despite their probabilisticnature these methods are nevertheless deterministic. Some convergencetheorems are proved. Numerical tests are presented.     

The probability approach to numerical solution of nonlinear parabolic equations

A number of new layer methods for solving semilinear parabolic equations and reaction‐diffusion systems is derived by using probabilistic representations of their solutions. These methods exploit the ideas of weak sense numerical integration of stochastic differential equations. In spite of the probabilistic nature these methods are nevertheless deterministic. A convergence theorem is proved. Some numerical tests are presented. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 490–522, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10020     

Numerical algorithms for semilinear parabolic equations with small parameter based on approximation of stochastic equations

The probabilistic approach is used for constructing special layer methods to solve the Cauchy problem for semilinear parabolic equations with small parameter. Despite their probabilistic nature these methods are nevertheless deterministic. The algorithms are tested by simulating the Burgers equation with small viscosity and the generalized KPP-equation with a small parameter.

     

A fourth order ADI method for semidiscrete parabolic equations

A fourth order fourstep ADI method is described for solving the systems of ordinary differential equations which are obtained when a (nonlinear) parabolic initial-boundary value problem in two dimensions is semi-discretized. The local time-discretization error and the stability conditions are derived. By numerical experiments it is demonstrated that the (asymptotic) fourth order behaviour does not degenerate if the time step increases to relatively large values. Also a comparison is made with the classical ADI method of Peaceman and Rachford showing the superiority of the fourth order method in the higher accuracy region, particularly in nonlinear problems.     

A link of stochastic differential equations to nonlinear parabolic equations

Using Girsanov transformation,we derive a new link from stochastic differential equations of Markovian type to nonlinear parabolic equations of Burgers-KPZ type,in such a manner that the obtained BurgersKPZ equation characterizes the path-independence property of the density process of Girsanov transformation for the stochastic differential equation.Our assertion also holds for SDEs on a connected differential manifold.     

A comparative study of ADI splitting methods for parabolic equations in two space dimensions

The main purpose of the paper is a numerical comparison of three integration methods for semi-discrete parabolic partial differential equations in two space variables. Linear as well as nonlinear,equations are considered. The integration methods are the well-known ADI method of Peaceman and Rachford, a global extrapolation scheme of the classical ADI method to order four and a fourth order, four-step ADI splitting method.     

A pathwise solution for nonlinear parabolic equations with stochastic perturbations

We analyse here a semilinear stochastic partial differential equation of parabolic type where the diffusion vector fields are depending on both the unknown function and its gradient ∂ _xu with respect to the state variable, ∈ ℝⁿ. A local solution is constructed by reducing the original equation to a nonlinear parabolic one without stochastic perturbations and it is based on a finite dimensional Lie algebra generated by the given diffusion vector fields.     

Two‐Grid method for nonlinear parabolic equations by expanded mixed finite element methods

In this article, we develop a two‐grid algorithm for nonlinear reaction diffusion equation (with nonlinear compressibility coefficient) discretized by expanded mixed finite element method. The key point is to use two‐grid scheme to linearize the nonlinear term in the equations. The main procedure of the algorithm is solving a small‐scaled nonlinear equations on the coarse grid and dealing with a linearized system on the fine space using the Newton iteration with the coarse grid solution. Error estimation to the expanded mixed finite element solution is analyzed in detail. We also show that two‐grid solution achieves the same accuracy as long as the mesh sizes satisfy H = O(h^1/2). Two numerical experiments are given to verify the effectiveness of the algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Application of stochastic equivalence for solving parabolic partial differential equations

An approach to solving parabolic partial differential equations based on the method of stochastic characteristics is proposed. The method allows decomposition of the numerical procedure into separate unified blocks. The approximation error and the efficiency of the method are evaluated. An example is given.     

Application of wavelet transforms to the solution of boundary value problems for linear parabolic equations

A method based on wavelet transforms is proposed for finding weak solutions to initial-boundary value problems for linear parabolic equations with discontinuous coefficients and inexact data. In the framework of multiresolution analysis, the general scheme for finite-dimensional approximation in the regularization method is combined with the discrepancy principle. An error estimate is obtained for the stable approximate solution obtained by solving a set of linear algebraic equations for the wavelet coefficients of the desired solution.     

Existence and approximation of nonlocal optimal design problems driven by parabolic equations

This work is a follow‐up to a series of articles by the authors where the same topic for the elliptic case is analyzed. In this article, a class of nonlocal optimal design problem driven by parabolic equations is examined. After a review of results concerning existence and uniqueness for the state equation, a detailed formulation of the nonlocal optimal design is given. The state equation is of nonlocal parabolic type, and the associated cost functional belongs to a broad class of nonlocal integrals. In the first part of the work, a general result on the existence of nonlocal optimal design is proved. The second part is devoted to analyzing the convergence of nonlocal optimal design problems toward the corresponding classical problem of optimal design. After a slight modification of the problem, either on the cost functional or by considering a new set of admissibility, the G‐convergence for the state equation and, consequently, the convergence of the nonlocal optimal design problem are proved.     

A nonlinear iteration method for solving a two-dimensional nonlinear coupled system of parabolic and hyperbolic equations

A nonlinear iteration method for solving a class of two-dimensional nonlinear coupled systems of parabolic and hyperbolic equations is studied. A simple iterative finite difference scheme is designed; the calculation complexity is reduced by decoupling the nonlinear system, and the precision is assured by timely evaluation updating. A strict theoretical analysis is carried out as regards the convergence and approximation properties of the iterative scheme, and the related stability and approximation properties of the nonlinear fully implicit finite difference (FIFD) scheme. The iterative algorithm has a linear constringent ratio; its solution gives a second-order spatial approximation and first-order temporal approximation to the real solution. The corresponding nonlinear FIFD scheme is stable and gives the same order of approximation. Numerical tests verify the results of the theoretical analysis. The discrete functional analysis and inductive hypothesis reasoning techniques used in this paper are helpful for overcoming difficulties arising from the nonlinearity and coupling and lead to a related theoretical analysis for nonlinear FI schemes.     

On global behavior of solutions to an inverse problem for nonlinear parabolic equations

We find conditions on data guaranteeing global nonexistence of solutions to an inverse source problem for a class of nonlinear parabolic equations. We also establish a stability result on a bounded domain for a problem with the opposite sign on the power type nonlinearity.     

A new approach to numerical solution of fixed-point problems and its application to delay differential equations

In this paper we consider a certain approximation of fixed-points of a continuous operator A mapping the metric space into itself by means of finite dimensional ε(h)-fixed-points of A. These finite dimensional functions are obtained from functions defined on discrete space grid points (related to a parameter h→0) by applying suitably chosen extension operators p_h. A theorem specifying necessary and sufficient conditions for existence of fixed-points of A in terms of ε(h)-fixed-points of A is given. A corollary which follows the theorem yields an approximate method for a fixed-point problem and determines conditions for its convergence. An example of application of the obtained general results to numerical solving of boundary value problems for delay differential equations is provided.Numerical experiments carried out on three examples of boundary value problems for second order delay differential equations show that the proposed approach produces much more accurate results than many other numerical methods when applied to the same examples.     

WELL-POSEDNESS OF A PARABOLIC INVERSE PROBLEM

1.IntroductionItiswellknownthatinverseproblemsinpartialdifferentialequations,mostofwhichhavenotyetbeensolveduptonow,remainasachallengeinappliedmathematics.Therefore,manymathematiciansstudiedvariousinverseproblemsforparabolicequa-tions.FOrasimplesurveywereferto[1,2,8,9]foridentifyingcoefficients,[7]foridentifyingboundaryvalues,[4,10]foridentifyingsourcetermsofparabolicequations.Wehavenotincludedalotofpapersconcerningthecomputationalmethodsusedforsolvinginverseparabolicproblems.Inthispaperthein…     

倒向随机方程期权定价模型的一类随机算法

期权定价问题可以转化为对倒向随机微分方程的求解,进而转化为对相应抛物型偏微分方程的求解.为了求解与倒向随机微分方程相应的二阶拟线性抛物型微分方程初值问题,引入一类新的随机算法-分层方法取代传统的确定性数值算法.这种数值方法理论上是通过弱显式欧拉法,离散其相应随机系统解的概率表示而得到.该随机算法的收敛性在文中得到证明,其稳定性是自然的.并构造了易于数值实现的基于插值的算法,实证研究说明这种算法能很好地提供期权定价模型的数值模拟.     

Linearly implicit methods for nonlinear parabolic equations

We construct and analyze combinations of rational implicit and explicit multistep methods for nonlinear parabolic equations. The resulting schemes are linearly implicit and include as particular cases implicit-explicit multistep schemes as well as the combination of implicit Runge-Kutta schemes and extrapolation. An optimal condition for the stability constant is derived under which the schemes are locally stable. We establish optimal order error estimates.

     

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     

Quasilinear degenerate parabolic equations of Kirchhoff type

We investigate the evolution problem where H is a Hilbert space, A is a self‐adjoint linear non‐negative operator on H with domain D(A), and is a continuous function. We prove that if , and , then there exists at least one global solution, which is unique if either m never vanishes, or m is locally Lipschitz continuous. Moreover, we prove that if for all , then this problem is well posed in H. On the contrary, if for some it happens that for all , then this problem has no solution if with β small enough. We apply these results to degenerate parabolic PDEs with non‐local non‐linearities. Copyright © 1999 John Wiley & Sons, Ltd.