Calculating the Malliavin derivative of some stochastic mechanics problems.
| Autor: | Paul Hauseux, Jack S Hale, Stéphane P A Bordas |
|---|---|
| Jazyk: | angličtina |
| Rok vydání: | 2017 |
| Předmět: | |
| Zdroj: | PLoS ONE, Vol 12, Iss 12, p e0189994 (2017) |
| Druh dokumentu: | article |
| ISSN: | 1932-6203 16173341 |
| DOI: | 10.1371/journal.pone.0189994 |
| Popis: | The Malliavin calculus is an extension of the classical calculus of variations from deterministic functions to stochastic processes. In this paper we aim to show in a practical and didactic way how to calculate the Malliavin derivative, the derivative of the expectation of a quantity of interest of a model with respect to its underlying stochastic parameters, for four problems found in mechanics. The non-intrusive approach uses the Malliavin Weight Sampling (MWS) method in conjunction with a standard Monte Carlo method. The models are expressed as ODEs or PDEs and discretised using the finite difference or finite element methods. Specifically, we consider stochastic extensions of; a 1D Kelvin-Voigt viscoelastic model discretised with finite differences, a 1D linear elastic bar, a hyperelastic bar undergoing buckling, and incompressible Navier-Stokes flow around a cylinder, all discretised with finite elements. A further contribution of this paper is an extension of the MWS method to the more difficult case of non-Gaussian random variables and the calculation of second-order derivatives. We provide open-source code for the numerical examples in this paper. |
| Databáze: | Directory of Open Access Journals |
| Externí odkaz: |
|