4462 articles – 13150 references  [version française]
HAL: hal-00677529, version 1

See detailed view  BibTeX,EndNote,...
Monte Carlo approximations of the Neumann problem
Sylvain Maire 1, 2, Etienne Tanré ( ) 2
(2012-03-08)

We introduce Monte Carlo methods to compute the solution of elliptic equations with pure Neumann boundary conditions. We first prove that the solution obtained by the stochastic representation has a zero mean value with respect to the invariant measure of the stochastic process associated to the equation. Pointwise approximations are computed by means of standard and new simulation schemes especially devised for local time approximation on the boundary of the domain. Global approximations are computed thanks to a stochastic spectral formulation taking into account the property of zero mean value of the solution. This stochastic formulation is asymptotically perfect in terms of conditioning. Numerical examples are given on the Laplace operator on a square domain with both pure Neumann and mixed Dirichlet-Neumann boundary conditions.
1:  Laboratoire des Sciences de l'Information et des Systèmes (LSIS)
CNRS : UMR6168 – Arts et Métiers ParisTech – Université Paul Cézanne - Aix-Marseille III – Université de la Méditerranée - Aix-Marseille II – Université de Provence - Aix-Marseille I – Université Sud Toulon Var
2:  TOSCA (INRIA Sophia Antipolis / INRIA Nancy - Grand Est/ IECN)
INRIA – CNRS : UMR7502 – Université de Lorraine
Mathematics/Probability
Attached file list to this document: 
PDF
MCNeumann2012.pdf(269.9 KB)