Simulation of fluid-particles flows: heavy particles, flowing regime, and asymptotic-preserving schemes (Q357515)

The following equations are considered: - the Vlasov-Fokker-Planck equation \[ \partial_t f+v\cdot\nabla_x f=\frac{1}{\varepsilon}L_{u,\theta}f+\eta_{P} \nabla_x \Phi \cdot\nabla_v f \] where \(L_{u,\theta}f=\mathrm{Div}_v ((v-u)f+\theta\nabla_v f),\) - the Euler system \[ \begin{cases} \partial_t \rho+\mathrm{Div}_x (\rho u)&=0, \\ \partial_t (\rho u)+\mathrm{Div}_x (\rho u\otimes u)+\nabla_x p&=\frac{1}{\varepsilon}\frac{\rho_P}{\rho_F}\mathcal{{F}}-\eta_F \rho\nabla_x \Phi, \\ \partial_t (\rho E)+\mathrm{Div}_x ((\rho E+p)u)&=\frac{1}{\varepsilon}\frac{\rho_P}{\rho_F}\mathcal{{E}}-\eta_F \rho u\cdot\nabla_x \Phi, \end{cases} \] in the case when the parameter \(\varepsilon\) can become small. In this regime, the conservation and dissipation properties are the basis to describe the asymptotic behavior of the model. A numerical scheme for the simulation of these equations is introduced. This scheme is able to capture the correct asymptotic behavior without suffering the prohibitive scaling-parameter-dependent numerical constraints, so this method belongs to the class of the asymptotic-preserving schemes. The method is based on the implicit treatment of the potentially stiff terms and relies on the possibility of updating the macroscopic unknowns by solving simple linear systems. Numerical examples are given. The ability of the scheme in preserving mass and energy and in dissipating entropy, the influence of the external force, asymptotic preserving properties, the temperature dependent viscosity, and effects of the state law and different kinetic approximations are considered.