Stochastic analysis for lagrangian particles simulation : application to colloidal particle collision
Analyse stochastique pour la simulation de particules lagrangiennes : application aux collisions de particules colloïdes
Résumé
This thesis broadly concerns colloidal particle simulation which plays an important role in understanding two-phase flows. More specifically, we track the particles inside a turbulent flow and model their dynamics as a stochastic process, their interactions as perfectly elastic collisions where the influence of the flow is modelled by a drift on the velocity term. By coupling each particle and considering their relative position and velocity, the perfectly elastic collision becomes a specular reflection condition. We put forward a time discretisation scheme for the resulting Lagrange system with specular boundary conditions and prove that the convergence rate of the weak error decreases at most linearly in the time discretisation step. The evidence is based on regularity results of the Feynman-Kac PDE and requires some regularity on the drift. We numerically experiment a series of conjectures, amongst which the weak error linearly decreasing for drifts that do not comply with the theorem conditions. We test the weak error convergence rate for a Richardson Romberg extrapolation. We finally deal with Lagrangian/Brownian approximations by considering a Lagrangian system where the velocity component behaves as a fast process. We control the weak error between the position of the Lagrangian system and an appropriately chosen uniformly elliptic diffusion process and subsequently prove a similar control by introducing a specular reflecting boundary on the Lagrangian and an appropriate reflection on the elliptic diffusion.
Cette thèse s'inscrit dans le cadre de la simulation de particules colloïdales. Plus précisément, nous nous intéressons aux particules dans un écoulement turbulent et modélisons leur dynamique par un processus lagrangien, leurs interactions comme des collisions parfaitement élastiques où l'influence de l'écoulement est modélisée par un terme de force sur la composante vitesse du système. En couplant les particules deux par deux et considérant leurs position et vitesse relatives, la collision parfaitement élastique devient une condition de réflexion spéculaire. Nous proposons un schéma de discrétisation en temps pour le système Lagrangien résultant avec des conditions aux bords spéculaires et prouvons que l'erreur faible diminue au plus linéairement dans le pas de discrétisation temporelle. La démonstration s’appuie sur des résultats de régularité de l'EDP Feynman-Kac et requiert une certaine régularité sur le terme de force. Nous expérimentons numériquement certaines conjectures, dont l’erreur faible diminuant linéairement pour des termes de force qui ne respectent pas les conditions du théorème. Nous testons le taux de convergence de l’erreur faible pour l’extrapolation de Romberg. Enfin, nous nous intéressons aux approximations Lagrangiennes/Browniennes en considérant un système Lagrangien où la composante vitesse se comporte comme un processus rapide. Nous contrôlons l'erreur faible entre la composante position du modèle Lagrangien et un processus de diffusion uniformément elliptique. Nous démontrons ensuite un contrôle similaire en introduisant une limite réfléchissante spéculaire sur le système Lagrangien et une réflexion appropriée sur la diffusion elliptique.
Origine : Version validée par le jury (STAR)
Loading...