A numerical study on the drag law of a gas bubble using dynamic body force method
Lingxin Zhang,Zecai Zhou,Jian Deng,Xueming Shao
DOI: https://doi.org/10.1063/5.0055646
IF: 4.6
2021-06-01
Physics of Fluids
Abstract:For a gas bubble rising in quiescent liquid, most previous drag force models correlate the drag coefficient <i>Cd</i> to dimensionless parameters such as the Reynolds number <i>Re</i>, the Eötvös number <i>Eo</i>, and the Morton number <i>Mo</i>. However, it is still an open question if the current models can be applied to non-buoyancy-driven bubbles in a wide range of parameters. In the present study, we investigate the effects of Reynolds number <i>Re</i> and Weber number <i>We</i> on the drag force coefficient <i>C<sub>d</sub></i> of bubbly flows. To obtain a converged drag force, a dynamic body force model is developed, which constraints the bubble at a fixed position. A large number of numerical simulations are carried out, spanning in the (<i>We</i>, <i>Re</i>) parametric space. The correlations of <i>C<sub>d</sub></i> to <i>We</i> and <i>Re</i> are studied, respectively. We find that for low <i>We</i>, as the bubble keeps spherical, its drag coefficients can be well predicted by the previous spherical bubble models. For medium–high <i>We</i>, when <i>Re</i> is low, the viscous force contributes most, and therefore, the bubble deformation or <i>We</i> plays an insignificant role on the drag. In contrast, when <i>Re</i> is high, the pressure force or shape drag becomes dominant, which is strongly dependent of <i>We</i>. Therefore, we can identify three different regions, spherical, <i>We</i>-dependent, and <i>We</i>-independent, in the parametric space. Finally, a new empirical drag force model for a non-buoyancy-driven bubble is proposed based on the <i>C<sub>d</sub></i>-<i>Re</i> and <i>C<sub>d</sub></i>-<i>We</i> correlations obtained by our numerical simulations. The current model is superior to the previous drag models, with predicting errors within 20%, by compared to the direct numerical simulation and experimental data. More profoundly, the current model can apply to a wider range of parameters, that is, <span class="equationTd inline-formula"><math> 1 ≤ R e ≤ 1000</math></span> and <span class="equationTd inline-formula"><math> 0 ≤ W e ≤ 20</math></span>, particularly for high <i>We</i>, when the bubble deformation cannot be neglected.
mechanics,physics, fluids & plasmas