Influence of bubble size, diffuser width and flow rate on the integral behaviour of bubble plumes

A large-eddy simulation based Eulerian-Lagrangian model is employed to quantify the impact of bubble size, diffuser diameter, and gas flow rate on integral properties of bubble plumes, such as the plume's width, centerline velocity, and mass flux. Calculated quantities are compared with experimental data and integral model predictions. Furthermore, the LES data were used to assess the behavior of the entrainment coefficient, the momentum amplification factor, and the bubble-to-momentum spread ratio. It is found that bubble plumes with constant bubble size and smaller diameter behave in accordance with integral plume models. Plumes comprising larger and non-uniform bubble sizes appear to deviate from past observations and model predictions. In multi-diameter bubble plumes, a bubble self-organisation takes place, i.e., small bubbles cluster in the center of the plume whilst large bubbles are found at the periphery of the plume. Multi-diameter bubble plumes also feature a greater entrainment rate than single-size bubble plumes, as well as a higher spread ratio and lower turbulent momentum rate. Once the plume is fully established, the size of the diffuser does not appear to affect integral properties of bubble plumes. However, plume development is affected by the diffuser width, as larger release areas lead to a delayed asymptotic behavior of the plume and consequently to a lower entrainment and higher spread ratio. Finally, the effect of the gas flow rate on the integral plume is studied and is deemed very relevant with regards to most integral plume properties and coefficients. This effect is already fairly well described by integral plume models

Secondary currents and turbulence over a non-uniformly roughened open-channel bed

Large-eddy simulations (LES) of the flow over a non-uniformly roughened channel bed are carried out to study the effect of non-uniform bed roughness on turbulence driven secondary currents and turbulence statistics. The channel bed is comprised of alternating rough and smooth strips, the width of which corresponds to the water depth. The Reynolds number based on hydraulic radius and bulk velocity is 34,000. The LES are successfully validated using experimental data. The secondary flow and bed roughness have a significant effect on the streamwise velocity and second order turbulence statistics. Turbulence is enhanced over rough strips and suppressed over smooth strips. Significant lateral momentum transfer takes place due to both advection and turbulence. The bed shear stresses over the smooth strips are approximately four times less than over the rough strips a result of near bed low momentum fluid being transported from the rough strips to the smooth strips and high momentum fluid being convected from the surface towards the bed. The most significant terms in the streamwise momentum equation are quantified and discussed with regard to momentum transfer

Scalability of an Eulerian-Lagrangian large-eddy simulation solver with hybrid MPI/OpenMP parallelisation

Eulerian-Lagrangian approaches capable of accurately reproducing complex fluid flows are becoming more and more popular due to the increasing availability and capacity of High Performance Computing facilities. However, the parallelisation of the Lagrangian part of such methods is challenging when a large number of Lagrangian markers are employed. In this study, a hybrid MPI/OpenMP parallelisation strategy is presented and implemented in a finite difference based large-eddy simulation code featuring the immersed boundary method which generally employs a large number of Lagrangian markers. A master-scattering-gathering strategy is used to deal with the handling of the Lagrangian markers and OpenMP is employed to distribute their computational load across several CPU threads. A classical domain-decomposition-based MPI approach is used to carry out the Eulerian, fixed-mesh fluid calculations. The results demonstrate that by using an effective combination of MPI and OpenMP the code can outperform a pure MPI parallelisation approach by up to 20%. Outcomes from this paper are of interest to various Eulerian-Lagrangian applications including the immersed boundary method, discrete element method or Lagrangian particle tracking