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

The following equations are considered:NEWLINENEWLINE- the Vlasov-Fokker-Planck equation NEWLINE\[NEWLINE\partial_t f+v\cdot\nabla_x f=\frac{1}{\varepsilon}L_{u,\theta}f+\eta_{P} \nabla_x \Phi \cdot\nabla_v f NEWLINE\]NEWLINE where \(L_{u,\theta}f=\mathrm{Div}_v ((v-u)f+\theta\nabla_v f),\)NEWLINENEWLINE- the Euler system NEWLINE\[NEWLINE\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}NEWLINE\]NEWLINE 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.NEWLINENEWLINEA numerical scheme for the simulation of these equations is introduced.NEWLINENEWLINEThis 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.NEWLINENEWLINEThe 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.NEWLINENEWLINENumerical 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.