Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)

 
 

 

Rate of convergence of a stochastic particle method for the Kolmogorov equation with variable coefficients


Authors: Pierre Bernard, Denis Talay and Luciano Tubaro
Journal: Math. Comp. 63 (1994), 555-587, S11
MSC: Primary 65M12; Secondary 35K57, 60J15, 60J60
DOI: https://doi.org/10.1090/S0025-5718-1994-1250770-3
MathSciNet review: 1250770
Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: In a recent paper, E. G. Puckett proposed a stochastic particle method for the nonlinear diffusion-reaction PDE in $ [0,T] \times \mathbb{R}$ (the so-called "KPP" (Kolmogorov-Petrovskii-Piskunov) equation):

$\displaystyle \left\{ \begin{array}{*{20}{c}} \frac{{\partial u}}{{\partial t}}... ... + f(u), \hfill \\ u(0, \cdot ) = {u_0}( \cdot ), \hfill \\ \end{array} \right.$

where $ 1 - {u_0}$ is the cumulative function, supposed to be smooth enough, of a probability distribution, and f is a function describing the reaction. His justification of the method and his analysis of the error were based on a splitting of the operator A. He proved that, if h is the time discretization step and N the number of particles used in the algorithm, one can obtain an upper bound of the norm of the random error on $ u(T,x)$ in $ {L^1}(\Omega \times \mathbb{R})$ of order $ 1/{N^{1/4}}$, provided $ h = \mathcal{O}(1/{N^{1/4}})$, but conjectured, from numerical experiments, that it should be of order $ \mathcal{O}h + \mathcal{O}(1/\sqrt N )$, without any relation between h and N.

We prove that conjecture. We also construct a similar stochastic particle method for more general nonlinear diffusion-reaction-convection PDEs

$\displaystyle \left\{ \begin{array}{*{20}{c}} \frac{{\partial u}}{{\partial t}} = Lu + f(u), \hfill \\ u(0,\cdot) = {u_0}(\cdot), \hfill \\ \end{array} \right.$

where L is a strongly elliptic second-order operator with smooth coefficients, and prove that the preceding rate of convergence still holds when the coefficients of L are constant, and in the other case is $ \mathcal{O}(\sqrt h ) + \mathcal{O}(1/\sqrt N )$.

The construction of the method and the analysis of the error are based on a stochastic representation formula of the exact solution u.


References [Enhancements On Off] (What's this?)

  • [1] A. Bensoussan and J. L. Lions, Applications des inéquations variationnelles en contrôle stochastique, Dunod, Paris, 1978. MR 0513618 (58:23923)
  • [2] P. Cannarsa and V. Vespri, Generation of analytic semigroups by elliptic operators with unbounded coefficients, SIAM J. Math. Anal. 18 (1987), 857-872. MR 883572 (88d:47053)
  • [3] B. Chauvin and A. Rouault, KPP equation and supercritical branching Brownian motion in the subcritical speed area. Application to spatial trees, Probab. Theory Related Fields 80 (1988), 299-314. MR 968823 (90b:60113)
  • [4] -, A stochastic simulation for solving scalar reaction-diffusion equations, Adv. Appl. Probab. 12 (1990), 88-100. MR 1039378 (91d:65009)
  • [5] -, Supercritical branching Brownian motion and K-P-P equation in the critical speed-area, Math. Nachr. 149 (1990), 41-59. MR 1124793 (92m:60075)
  • [6] A. Friedman, Stochastic differential equations and applications, Vol. 1, Academic Press, New York, 1975. MR 0494490 (58:13350a)
  • [7] T. Hida, Brownian motion, Springer-Verlag, Berlin and New York, 1980. MR 562914 (81a:60089)
  • [8] H. Kunita, Stochastic differential equations and stochastic flows of diffeomorphisms, Ecole d'Eté de Saint-Flour XII, Lecture Notes in Math., vol. 1097, Springer-Verlag, Berlin and New York, 1984. MR 876080 (87m:60127)
  • [9] G. N. Milshtein, Approximate integration of stochastic differential equations, Theory Probab. Appl. 19 (1974), 557-562. MR 0356225 (50:8696)
  • [10] A. Pazy, Semigroups of linear operators and applications to partial differential equations, Springer-Verlag, Berlin and New York, 1983. MR 710486 (85g:47061)
  • [11] G. Da Prato and E. Sinestrari, Differential operators with non dense domain, Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4) 14 (1987), 285-344. MR 939631 (89f:47062)
  • [12] E. G. Puckett, Convergence of a random particle method to solutions of the Kolmogorov equation, Math. Comp. 52 (1989), 615-645. MR 964006 (90h:65008)
  • [13] F. Rothe, Global solutions of reaction-diffusion systems, Lecture Notes in Math., vol. 1072, Springer-Verlag, Berlin and New York, 1984. MR 755878 (86d:35071)
  • [14] A. S. Sherman and C. S. Peskin, A Monte Carlo method for scalar reaction diffusion equations, SIAM J. Sci. Statist. Comput. 7 (1986), 1360-1372. MR 857799
  • [15] B. Stewart, Generation of analytic semigroups by strongly elliptic operators, Trans. Amer. Math. Soc. 199 (1974), 141-162. MR 0358067 (50:10532)
  • [16] D. Talay, Simulation and numerical analysis of stochastic differential systems: a review, Rapport de Recherche INRIA, vol. 1313, 1990 (and to appear in Effective Stochastic Analysis (P. Kree and W. Wedig, eds.), Springer-Verlag).

Similar Articles

Retrieve articles in Mathematics of Computation with MSC: 65M12, 35K57, 60J15, 60J60

Retrieve articles in all journals with MSC: 65M12, 35K57, 60J15, 60J60


Additional Information

DOI: https://doi.org/10.1090/S0025-5718-1994-1250770-3
Article copyright: © Copyright 1994 American Mathematical Society

American Mathematical Society