Showing 1–2 of 2 results for author: Man, P L W
-
A Stochastic Algorithm for Parametric Sensitivity in Smoluchowski's Coagulation Equation
Authors:
I. Bailleul,
P. L. W. Man,
M. Kraft
Abstract:
In this article a stochastic particle system approximation to the parametric sensitivity in the Smoluchowski coagulation equation is introduced. The parametric sensitivity is the derivative of the solution to the equation with respect to some parameter, where the coagulation kernel depends on this parameter. It is proved that the particle system converges weakly to the sensitivity as the number…
▽ More
In this article a stochastic particle system approximation to the parametric sensitivity in the Smoluchowski coagulation equation is introduced. The parametric sensitivity is the derivative of the solution to the equation with respect to some parameter, where the coagulation kernel depends on this parameter. It is proved that the particle system converges weakly to the sensitivity as the number of particles N increases. A Monte Carlo algorithm is developed and variance reduction techniques are applied. Numerical experiments are conducted for two kernels: the additive kernel and one which has been used for studying soot formation in a free molecular regime. It is shown empirically that the techniques for variance reduction are indeed very effective and that the order of convergence is O(1/N). The algorithm is then compared to an algorithm based on a finite difference approximation to the sensitivity and it is found that the variance of the sensitivity estimators are considerably lower than that for the finite difference approach. Furthermore, two methods of establishing `efficiency' are considered and the new algorithm is found to be significantly more efficient.
△ Less
Submitted 21 October, 2009;
originally announced October 2009.
-
Coupling algorithm for calculating sensitivities of Smoluchowski's coagulation equation
Authors:
Peter L. W. Man,
James R. Norris,
Ismael F. Bailleul,
Markus Kraft
Abstract:
In this paper, two new stochastic algorithms for calculating parametric derivatives of the solution to the Smoluchowski coagulation equation are presented. It is assumed that the coagulation kernel is dependent on these parameters. The new algorithms (called `Single' and `Double') work by coupling two Marcus-Lushnikov processes in such a way as to reduce the difference between their trajectories…
▽ More
In this paper, two new stochastic algorithms for calculating parametric derivatives of the solution to the Smoluchowski coagulation equation are presented. It is assumed that the coagulation kernel is dependent on these parameters. The new algorithms (called `Single' and `Double') work by coupling two Marcus-Lushnikov processes in such a way as to reduce the difference between their trajectories, thereby significantly reducing the variance of central difference estimators of the parametric derivatives. In the numerical results, the algorithms are shown have have a O(1/N) order of convergence as expected, where N is the initial number of particles. It was also found that the Single and Double algorithms provide much smaller variances. Furthermore, a method for establishing `efficiency' is considered, which takes into account the variances as well as CPU run times, and the `Double' is significantly more `efficient' compared to the `Independent' algorithm in most cases.
△ Less
Submitted 20 October, 2009;
originally announced October 2009.