Random vortex and expansion-rate model for Oberbeck-Boussinesq fluid flows
Zihao Guo, Zhongmin Qian, Zihao Shen
TL;DR
This work develops a vorticity–expansion-rate framework to simulate compressible Oberbeck–Boussinesq-type flows using a random vortex method. It combines a new parabolic integral representation with Biot–Savart constructions to obtain stochastic representations for velocity, temperature, and density, and extends the approach to wall-bounded domains via reflection and cut-off techniques. The method yields implicit McKean–Vlasov–type representations and decoupled SDEs in a constant-density limit, enabling Monte Carlo simulations that capture Benard-like convection and density–temperature coupling. Numerical experiments in a half-space demonstrate boundary-driven vortex dynamics and small density variations consistent with near-incompressible OB physics, supported by error-term analysis in the Appendix.
Abstract
By using a formulation of a class of compressible viscous flows with a heat source via vorticity and expansion-rate, we study the Oberbeck-Boussinesq flows. To this end we establish a new integral representation for solutions of parabolic equations subject to certain boundary condition, which allows us to develop a random vortex method for certain compressible flows and to compute numerically solutions of their dynamical models. Numerical experiments are carried out, which not only capture detailed Bénard convection but also are capable of providing additional information on the fluid density and the dynamics of expansion-rate of the flow.
