A lattice Boltzmann study on the drag force in bubble swarms

2011 ◽  
Vol 679 ◽  
pp. 101-121 ◽  
Author(s):  
J. J. J. GILLISSEN ◽  
S. SUNDARESAN ◽  
H. E. A. VAN DEN AKKER

Lattice Boltzmann and immersed boundary methods are used to conduct direct numerical simulations of suspensions of massless, spherical gas bubbles driven by buoyancy in a three-dimensional periodic domain. The drag coefficient CD is computed as a function of the gas volume fraction φ and the Reynolds number Re = 2RUslip/ν for 0.03 φ 0.5 and 5 Re 2000. Here R, Uslip and ν denote the bubble radius, the slip velocity between the liquid and the gas phases and the kinematic viscosity of the liquid phase, respectively. The results are rationalized by assuming a similarity between the CD(Reeff)-relation of the suspension and the CD(Re)-relation of an individual bubble, where the effective Reynolds number Reeff = 2RUslip/νeff is based on the effective viscosity νeff which depends on the properties of the suspension. For Re ≲ 100, we find νeff ≈ ν/(1−0.6φ1/3), which is in qualitative agreement with previous proposed correlations for CD in bubble suspensions. For Re ≳ 100, on the other hand, we find νeff ≈ RUslipφ, which is explained by considering the turbulent kinetic energy levels in the liquid phase. Based on these findings, a correlation is constructed for CD(Re, φ). A modification of the drag correlation is proposed to account for effects of bubble deformation, by the inclusion of a correction factor based on the theory of Moore (J. Fluid Mech., vol. 23, 1995, p. 749).

Author(s):  
Orest Shardt ◽  
Jos Derksen

We describe the direct simulation of high-solids-fraction suspensions of non-spherical rigid particles that are slightly denser than the fluid. The lattice-Boltzmann method is used to solve the flow of the interstitial Newtonian fluid, and the immersed boundary method is used to enforce a no-slip boundary condition at the surface of each particle. The surface points for the immersed boundary method are also employed for collision handling by applying repulsive forces between nearby surface points. Due to the finite number of these points, the method simulates rough surface collisions. We also discuss methods for integrating the equations of particle motion at low density ratios and propose a method with improved accuracy. Rigid particles shaped like red blood cells were simulated. Simulations of a single particle showed that the particle settles in its original orientation when the Reynolds number is low (1.2) but flips to a higher drag orientation when the Reynolds number is higher (7.3). A simulation with a 45% solids volume fraction and a low solid over fluid density ratio showed the possibility of simulating blood as it is found in the body. A simulation at a lower solids volume fraction (35%) was used to compare the results with the erythrocyte sedimentation rate (ESR), a common blood test. The sedimentation rate was estimated as 0.2 mm/hr, which is an order of magnitude lower than a typical ESR of about 6 mm/hr for a healthy adult. The most likely reasons for the discrepancy are the omission of agglomeration-inducing inter-cellular forces from the simulations and the treatment of the red blood cells as rigid particles.


2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Wim-Paul Breugem ◽  
Vincent van Dijk ◽  
René Delfos

Two different direct-forcing immersed boundary methods (IBMs) were applied for the purpose of simulating slow flow through a real porous medium: the volume penalization IBM and the stress IBM. The porous medium was a random close packing of about 9000 glass beads in a round tube. The packing geometry was determined from an X-ray computed tomography (CT) scan in terms of the distribution of the truncated solid volume fraction (either 0 or 1) on a three-dimensional Cartesian grid. The scan resolution corresponded to 19.3 grid cells over the mean bead diameter. A facility was built to experimentally determine the permeability of the packing. Numerical simulations were performed for the same packing based on the CT scan data. For both IBMs the numerically determined permeability based on the Richardson extrapolation was just 10% lower than the experimentally found value. As expected, at finite grid resolution the stress IBM appeared to be the most accurate IBM.


10.30544/455 ◽  
2020 ◽  
Vol 26 (1) ◽  
pp. 103-119
Author(s):  
Kamel Chadi ◽  
Nourredine Belghar ◽  
Belhi Guerira ◽  
Aissam Messaoudi

In the present work, we have studied the thermal exchanges of different geometry sections of mini-channels of a cooler numerically. Particularly, we have chosen a mini channels cooler copper for cooling an electronic chip IGBT. In our simulation of three-dimensional (3D), we have compared the numerical results for the different forms of the proposed mini-channels and the three different types of nano-fluids by using the Cu-water, the Ag-water, and the Diamond-water with a volume fraction of 0.02%. The numerical results are obtained by choosing a Reynolds number (Re) between 100 and 900 and considering that the flow regime is stationary. The simulation was performed using commercial software, ANSYS-Fluent 15.0. The results obtained show that the increase of the exchange surface between the walls of the mini channels and the cooling fluid makes increases the heat exchange coefficient and the improvement of the maximum junction temperature of the electronic chip IGBT with the increase of the Reynolds number. The choice of nanoparticles has considerable effects on improving the heat transfer and the maximum junction temperature of the chip IGBT.


2020 ◽  
pp. 96-96
Author(s):  
Abchouyeh Asadi ◽  
Ganaoui El ◽  
Rasul Mohebbi ◽  
Mohammad Zarrabi ◽  
Omid Fard ◽  
...  

In this study, the Lattice Boltzmann Method (LBM) is employed in order to examine the fluid flow and forced convection heat transfer inside a two-dimensional horizontal channel with and without obstacles. In order to enhance the heat and thermal energy transfer within the channel, different obstacle arrangements are posed to the flow field and heat transfer with the purpose of studying their sensitivity to these changes. The results indicate that, when the value of the Reynolds number is maximum, the maximum average Nusselt numbers happens on the lower wall (Case 4). The paper extends the topic to the use of nanofluids to introduce a possibility to enhancement of the heat transfer in the channel with an array of the obstacles with forced convection. For this purpose, the AgMgO/water micropolar hybrid nanofluid is used, and the volume fraction of the nanoparticle (50% Ag and 50% MgO by volume) is set between 0 and 0.02. The results showed that, when the hybrid nanofluid is used instead of a typical nanofluid, the rate of the heat transfer inside the channel increases, especially for the high values of the Reynolds number, and the volume fraction of the nanoparticles. Increasing the volume fraction of the nanoparticles increase the local Nusselt number ( 1.17-fold). It is shown that the type of obstacle arrangement and the specific nanofluid can exerts significant effects on the characteristics of the flow field and heat transfer in the channel. This study provides a platform for using the LBM to examine fluid flow through discrete obstacles in offset positions.


1999 ◽  
Vol 385 ◽  
pp. 41-62 ◽  
Author(s):  
DEWEI QI

A lattice-Boltzmann method has been developed to simulate suspensions of both spherical and non-spherical particles in finite-Reynolds-number flows. The results for sedimentation of a single elliptical particle are shown to be in excellent agreement with the results of Huang, Hu & Joseph (1998) who used a finite-element method. Sedimentation of two-dimensional circular and rectangular particles in a two-dimensional channel and three-dimensional spherical particles in a tube with square cross-section is simulated. Computational results are consistent with experimentally observed phenomena, such as drafting, kissing and tumbling.


Sign in / Sign up

Export Citation Format

Share Document