Fluid-structure interaction for the flexible filament's propulsion hanging in the free stream

Transferring Bubble Breakage Models Tailored for Euler-Euler Approaches to Euler-Lagrange Simulations

Most bubble breakage models have been developed for multiphase simulations using Euler-Euler (EE) approaches. Commonly, they are linked with population balance models (PBM) and are validated by making use of Reynolds-averaged Navier-Stokes (RANS) turbulence models. The latter, however, may be replaced by alternate approaches such as Large Eddy simulations (LES) that play a pivotal role in current developments based on lattice Boltzmann (LBM) technologies. Consequently, this study investigates the possibility of transferring promising bubble breakage models from the EE framework into Euler-Lagrange (EL) settings aiming to perform LES. Using our own model, it was possible to reproduce similar bubble size distributions (BSDs) for EL and EE simulations. Therefore, the critical Weber (Wecrit) number served as a threshold value for the occurrence of bubble breakage events. Wecrit depended on the bubble daughter size distribution (DSD) and a set minimum time between two consecutive bubble breakage events. The commercial frameworks Ansys Fluent and M-Star were applied for EE and EL simulations, respectively. The latter enabled the implementation of LES, i.e., the use of a turbulence model with non-time averaged entities. By properly choosing Wecrit, it was possible to successfully transfer two commonly applied bubble breakage models from EE to EL. Based on the mechanism of bubble breakage, Wecrit values of 7 and 11 were determined, respectively. Optimum Wecrit were identified as fitting the shape of DSDs, as this turned out to be a key criterion for reaching optimum prediction quality. Optimum Wecrit values hold true for commonly applied operational conditions in aerated bioreactors, considering water as the matrix.

Suppressing Viscous Fingering in Porous Media with Wetting Gradient

The viscous fingering phenomenon often occurs when a low-viscosity fluid displaces a high-viscosity fluid in a homogeneous porous media, which is an undesirable displacement process in many engineering applications. The influence of wetting gradient on this process has been studied over a wide range of capillary numbers (7.5 × 10^-6 to 1.8 × 10^-4), viscosity ratios (0.0025 to 0.04), and porosities (0.48 to 0.68), employing the lattice Boltzmann method. Our results demonstrate that the flow front stability can be improved by the gradual increase in wettability of the porous media. When the capillary number is less than 3.5 × 10^-5, the viscous fingering can be successfully suppressed and the transition from unstable to stable displacement can be achieved by the wetting gradient. Moreover, under the conditions of high viscosity ratio (M > 0.01) and large porosity (Φ > 0.58), wetting gradient improves the stability of the flow front more significantly.

An Explicit-Correction-Force Scheme of IB-LBM Based on Interpolated Particle Distribution Function

In order to obtain a better numerical simulation method for fluid-structure interaction (FSI), the IB-LBM combining the lattice Boltzmann method (LBM) and immersed boundary method (IBM) has been studied more than a decade. For this purpose, an explicit correction force scheme of IB-LBM was proposed in this paper. Different from the current IB-LBMs, this paper introduced the particle distribution function to the interpolation process from the fluid grids to the immersed boundary at the mesoscopic level and directly applied the LBM force models to obtain the interface force with a simple form and explicit process. Then, in order to ensure the mass conservation in the local area of the interface, this paper corrected the obtained interface force with the correction matrix, forming the total explicit-correction-force (ECP) scheme of IB-LBM. The results of four numerical tests were used to verify the order of accuracy and effectiveness of the present method. The streamline penetration is limited and the numerical simulation with certain application significance is successful for complex boundary conditions such as the movable rigid bodies (free oscillation of the flapping foil) and flexible deformable bodies (free deformation of cylinders). In summary, we obtained a simple and alternative simulation method that can achieve good simulation results for engineering reference models with complex boundary problems.