An Euler–Lagrange method considering bubble radial dynamics for modeling sonochemical reactors

Unsteady numerical computations are performed to investigate the flow field, wave propagation and the structure of bubbles in sonochemical reactors. The turbulent flow field is simulated using a two-equation Reynolds-Averaged Navier–Stokes (RANS) model. The distribution of the acoustic pressure is solved based on the Helmholtz equation using a finite volume method (FVM). The radial dynamics of a single bubble are considered by applying the Keller–Miksis equation to consider the compressibility of the liquid to the first order of acoustical Mach number. To investigate the structure of bubbles, a one-way coupling Euler–Lagrange approach is used to simulate the bulk medium and the bubbles as the dispersed phase. Drag, gravity, buoyancy, added mass, volume change and first Bjerknes forces are considered and their orders of magnitude are compared. To verify the implemented numerical algorithms, results for one- and two-dimensional simplified test cases are compared with analytical solutions. The results show good agreement with experimental results for the relationship between the acoustic pressure amplitude and the volume fraction of the bubbles. The two-dimensional axi-symmetric results are in good agreement with experimentally observed structure of bubbles close to sonotrode.