Abstract
The Regularised Inertial Dean–Kawasaki model (RIDK) – introduced by the authors and J. Zimmer in earlier works – is a nonlinear stochastic PDE capturing fluctuations around the meanfield limit for large-scale particle systems in both particle density and momentum density. We focus on the following two aspects. Firstly, we set up a Discontinuous Galerkin (DG) discretisation scheme for the RIDK model: we provide suitable definitions of numerical fluxes at the interface of the mesh elements which are consistent with the wave-type nature of the RIDK model and grant stability of the simulations, and we quantify the rate of convergence in mean square to the continuous RIDK model. Secondly, we introduce modifications of the RIDK model in order to preserve positivity of the density (such a feature only holds in a “high-probability sense” for the original RIDK model). By means of numerical simulations, we show that the modifications lead to physically realistic and positive density profiles. In one case, subject to additional regularity constraints, we also prove positivity. Finally, we present an application of our methodology to a system of diffusing and reacting particles. Our Python code is available in open-source format.
Reference50 articles.
1. Abramowitz M. and Stegun I.A., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Vol. 55. US Government Printing Office (1964).
2. Particle approximation of the Wasserstein diffusion
3. Arnold D.N., Brezzi F., Cockburn B. and Marini D., Discontinuous Galerkin methods for elliptic problems, in Discontinuous Galerkin Methods. Springer, Berlin Heidelberg (2000) 89–101.
4. Multigrid in H (div) and H (curl)
5. Baňas L., Gess B. and Vieth C., Numerical approximation of singular-degenerate parabolic stochastic PDEs. Preprint arXiv:2012.12150 (2020).
Cited by
2 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献