Gradient-based Monte Carlo methods for hyperbolic and kinetic equations

Lorenzo Pareschi
Heriot-Watt University
Mathematics

Particle methods based on evolving the spatial derivatives of the solution were originally introduced to simulate reaction-diffusion processes, inspired by vortex methods for the Navier--Stokes equations. Such methods, referred to as gradient random walk methods, were extensively studied in the '90s and have several interesting features, such as being grid free, automatically adapting to the solution by concentrating elements where the gradient is large and significantly reducing the variance of the standard random walk approach. In this talk, we revive these ideas by showing how to generalize the approach to a larger class of partial differential equations, including hyperbolic systems of conservation laws and kinetic equations.
This is a joint work with G. Bertaglia and R. Caflisch

Presentation (PDF File)

Back to International Conference on Multiscale Modeling and Simulation based on Physics and Data