Nathan Quadrio, "Multilevel Monte Carlo method for PDEs with fluid dynamic applications"
This thesis focuses on PDEs in which some of the parameters are not known exactly but affected by a certain amount of uncertainty, and hence described in terms of random variables/random fields. This situation is quite common in engineering applications.
A common goal in this framework is to compute statistical indices, like mean or variance, for some quantities of interest related to the solution of the equation at hand (“uncertainty quantification"). The main challenge in this task is represented by the fact that in many applications tens/hundreds of random variables may be necessary to obtain an accurate representation of the solution variability. The numerical schemes adopted to perform the uncertainty quantification should then be designed to reduce the degradation of their performance whenever the number of parameters increases, a phenomenon known as the “curse of dimensionality".
A method that acts in this direction is Monte Carlo sampling. Such a method is known to be dimension independent and very robust but it is also known for its very slow convergence rate. In this work, we describe Monte Carlo sampling together with solid error analysis and we provide a test for its robustness by integrating different functions with different regularities and by solving different PDE problems with random coefficients.
Later on, we introduce the technique of variance reduction and a further application of this idea that goes with the name of the Multilevel Monte Carlo (MLMC) method. The asymptotic cost of solving the stochastic problem with the multilevel method is proved to be significantly lower than that of the standard method and, in certain circumstances, grows only proportionally with respect to the cost of solving the deterministic problem. Numerical calculations demonstrating its effectiveness are presented and more complex problems such as elliptic PDEs in a domain with random geometry are also presented, this last task has been performed to test a code designed for the simulation of flow through porous media at the pore scale. The results are promising considering the very complex geometries that require extremely expensive discretizations.
This work is the final outcome of the participation in the Visiting Student Research Program at the King Abdullah University of Science and Technology, a project that allows students to conduct research with faculty mentors in selected areas of pure and applied sciences, and a Visiting Research Fellowship at the University of Texas in Austin.