Mathematics > Probability
[Submitted on 21 Oct 2009]
Title:A Stochastic Algorithm for Parametric Sensitivity in Smoluchowski's Coagulation Equation
View PDFAbstract: 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.
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.