Numerical Study of the Movement of Single Fine Particles in Porous Media Based on LBM-DEM
Abstract
:1. Introduction
2. Mathematic Model
2.1. Model of LBM
2.2. Solid–Solid Coupled Model
2.3. Liquid–Solid Coupled Model
2.4. Dynamics Model of Discrete Particle
3. Physical Model and Numerical Simulation Conditions
4. Model Validation
5. Numerical Simulation Results and Analysis
5.1. Migrations and Sedimentations of Fine Particles in Porous Media under Different Particle Densities
5.2. Migration and Sedimentation of a Fine Particle in Porous Media under Different Diameters
5.3. Effect of Fine Particles on the Flow Field in Porous Media
5.4. Migration and Sedimentation of Two Fine Particles in Porous Media
5.5. Migration and Sedimentation of Three Fine Particles in Porous Media
6. Summary and Conclusions
- (1)
- A fine particle with a density of 2.1–2.5 g/cm3 can travel longer distance under fluid drag, gravity and collision force, although it is more likely to deviate from the channel and collide with the internal wall of porous media; however, a fine particle with a density of 3.0 g/cm3 will be deposited on an internal wall after a shorter route because of its greater gravity.
- (2)
- The migration distance of a fine particle is not linearly correlated with its density. A fine particle with a density less than 2.0 g/cm3 is liable to stagnate in a low-infiltration area in porous media.
- (3)
- A fine particle with a smaller diameter is also likely to be stranded in a low-infiltration area even though its Stokes number is small.
- (4)
- A fine particle blocking the pore throat in porous media results in an evident pressure fluctuation of the flow field, intensifying the momentum transfer between the flow field and fine particles. An increase in diameter of a fine particle causes a greater increase in the pressure drop of the liquid–solid flow compared to an increase in the fine particle’s density.
- (5)
- If the horizontal distance between two non-contact fine particles is less than 0.25 cm, the presence of one fine particle may improve the migration of the other, and the other will travel a longer distance in the flow field in porous media.
- (6)
- When the pressure drop of the liquid–solid flow approaches a stable state, the values of the pressure drop of the liquid–solid flow containing two or three fine particles in porous media are slightly greater than that of the liquid–solid flow containing only one fine particle.
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
Nomenclature
distribution function | |
equilibrium distribution function | |
velocity vector | |
external force term | |
transition function | |
lattice sound speed | |
lattice spacing | |
time step | |
pressure of the fluid | |
macroscopic fluid velocity | |
normal contact force (N) | |
tangential contact force (N) | |
normal stiffness coefficient | |
unit vector from particle i to particle j | |
restitution coefficient determined by experiment | |
tangential stiffness coefficient | |
slip velocity of contact point (m/s) | |
, | radii of particles i and j (m) |
tangential unit vector | |
contact force (N) | |
contact torque (N m) | |
local force acting on fluid | |
Lagrange force density | |
velocity of the particle boundary discrete point | |
velocity pre-collision | |
matrix with interpolation items | |
matrix with extensions | |
, | fluid Euler coordinates |
, | solid Lagrange coordinates |
dr | thickness of the solid forcing shell |
ds | distance between two adjacent Lagrangian points on the solid boundary |
velocity of the Euler point of the fluid | |
fluid force of particle | |
fluid torque of particle | |
number of Lagrange points on the particle boundary | |
number of local Euler meshes around the particle | |
Lagrange point coordinate | |
center coordinate of the particle mass | |
mass of solid particles (kg) | |
rotational inertia of solid particles (kg·m2) | |
gravity (N) | |
υ | kinematic viscosity (m2/s) |
Cd | resistance coefficient |
Fd | resistance of the particle (N) |
pressure drop (Pa) | |
t | time (s) |
k | permeability |
L | length of the porous media |
W | width of the porous media |
ki | initial permeability of porous media. |
Greek Symbols | |
discrete velocity direction | |
dimensionless relaxation time | |
fluid density | |
, , | model parameters |
normal overlap of the contact particles | |
, | normal and tangential damping coefficients |
tangential displacement | |
tangential friction coefficient of the particles | |
weight coefficient | |
, | angular velocities of particle i and j (rad/s) |
Dirac function | |
porosity | |
µ | viscosity of fluid |
References
- Conley, G.; Beck, N.; Riihimaki, C.A.; Tanner, M. Quantifying clogging patterns of infiltration systems to improve urban stormwater pollution reduction estimates. Water Res. X 2020, 7, 100049. [Google Scholar] [CrossRef]
- Verma, S.; Daverey, A.; Sharma, A. Slow sand filtration for water and wastewater treatment—A review. Environ. Technol. Rev. 2017, 6, 47–58. [Google Scholar] [CrossRef]
- Zhang, B.L.; Zhou, S.; Li, C.L.; Zhang, C.K.; Zhang, C. A research on waste water treatment method of marine desulfurization washing system. Adv. Mater. Res. 2013, 803, 43–46. [Google Scholar] [CrossRef]
- Zhang, B.L.; Zhou, S.; Li, C.L.; Zhang, C.K.; Zhang, C. The status of exploitation techniques of natural gas hydrate. Chin. J. Chem. Eng. 2019, 27, 2133–2147. [Google Scholar]
- Jiménez, S.; Micó, M.; Arnaldos, M.; Medina, F.; Contreras, S. State of the art of produced water treatment. Chemosphere 2018, 192, 186–208. [Google Scholar] [CrossRef]
- Gitis, V.; Rubinstein, I.; Livshits, M.; Ziskind, G. Deep-bed filtration model with multistage deposition kinetics. Chem. Eng. J. 2010, 163, 78–85. [Google Scholar] [CrossRef]
- Khuzhayorov, B.; Fayziev, B.; Ibragimov, G.; Arifin, N.M. A deep bed filtration model of two-component suspension in dual-zone porous medium. Appl. Sci. 2020, 10, 2793. [Google Scholar] [CrossRef]
- Alem, A.; Elkawafi, A.; Ahfir, N.-D.; Wang, H. Filtration of kaolinite particles in a saturated porous medium: Hydrodynamic effects. Hydrogeol. J. 2013, 21, 573–586. [Google Scholar] [CrossRef]
- Zhang, P.; Bai, B.; Jiang, S. Coupled effects of hydrodynamic forces and pore structure on suspended particle transport and deposition in a saturated porous medium. Rock Soil Mech. 2016, 37, 1307–1316. [Google Scholar]
- Du, X.; Ye, X.; Zhang, X. Clogging of saturated porous media by silt-sized suspended solids under varying physical conditions during managed aquifer recharge. Hydrol. Process. 2018, 32, 2254–2262. [Google Scholar] [CrossRef]
- Abbar, B.; Alem, A.; Pantet, A.; Marcotte, S.; Ahfir, N.-D.; Duriatti, D. Experimental investigation on removal of suspended particles from water using flax fibre geotextiles. Environ. Technol. 2017, 38, 2964–2978. [Google Scholar] [CrossRef]
- Bennacer, L.; Ahfir, N.; Alem, A.; Wang, H. Coupled effects of ionic strength, particle size, and flow velocity on transport and deposition of suspended particles in saturated porous media. Transp. Porous Media 2017, 118, 251–269. [Google Scholar] [CrossRef]
- Park, H.; Pan, X.; Lee, C.; Choi, J.-I. A pre-conditioned implicit direct forcing based immersed boundary method for incompressible viscous flows. J. Comput. Phys. 2016, 314, 774–799. [Google Scholar] [CrossRef]
- Wang, S.; Li, H.; Wang, R.; Tian, R.; Sun, Q.; Ma, Y. Numerical simulation of flow behavior of particles in a porous media based on CFD-DEM. J. Pet. Sci. Eng. 2018, 171, 140–152. [Google Scholar] [CrossRef]
- Mirzaee, H.; Henn, T.; Krause, M.J.; Goubergrits, L.; Schumann, C.; Neugebauer, M.; Kuehne, T.; Preusser, T.; Hennemuth, A. MRI-based computational hemodynamics in patients with aortic coarctation using the lattice Boltzmann methods: Clinical validation study. J. Magn. Reson. Imaging JMRI 2017, 45, 139–146. [Google Scholar] [CrossRef]
- Agarwal, G.; Singh, S. Coupling of first collision source method with Lattice Boltzmann Method for the solution of neutron transport equation. Ann. Nucl. Energy 2022, 175, 109242. [Google Scholar] [CrossRef]
- Fu, S.; Su, W.; Zhang, H.; Wang, L. An immersed moving boundary for fast discrete particle simulation with complex geometry. Chem. Eng. Sci. 2024, 283, 119407. [Google Scholar] [CrossRef]
- Kempe, T.; Fröhlich, J. An improved immersed boundary method with direct forcing for the simulation of particle laden flows. J. Comput. Phys. 2012, 231, 3663–3684. [Google Scholar] [CrossRef]
- Tavanashad, V.; Subramaniam, S. Fully resolved simulation of dense suspensions of freely evolving buoyant particles using an improved immersed boundary method. Int. J. Multiph. Flow 2020, 132, 103396. [Google Scholar] [CrossRef]
- Alzetta, G. Multiscale modeling of fiber reinforced materials via non-matching immersed methods. Comput. Struct. 2020, 239, 106334. [Google Scholar] [CrossRef]
- Jiang, F.; Liu, H.; Chen, X.; Tsuji, T. A coupled LBM-DEM method for simulating the multiphase fluid-solid interaction problem. J. Comput. Phys. 2022, 454, 110963. [Google Scholar] [CrossRef]
- Ma, H.; Zhao, Y. An approach to distribute the marker points on non-spherical particle/boundary surface within the IBM-LBM framework. Eng. Anal. Bound. Elem. 2019, 108, 254–266. [Google Scholar] [CrossRef]
- Farahani, M.V.; Foroughi, S.; Norouzi, S.; Jamshidi, S. Mechanistic study of fines migration in porous media using lattice Boltzmann method coupled with rigid body physics engine. J. Energy Resour. Technol. 2019, 141, 58. [Google Scholar]
- Li, Y.-X.; Alshbool, M.H.; Lv, Y.-P.; Khan, I.; Khan, M.R.; Issakhov, A. Heat and mass transfer in MHD Williamson nanofluid flow over an exponentially porous stretching surface. Case Stud. Therm. Eng. 2021, 26, 100975. [Google Scholar] [CrossRef]
- Khan, M.R.; Elkotb, M.A.; Matoog, R.T.; Alshehri, N.A.; Abdelmohimen, M.A.H. Thermal features and heat transfer enhancement of a Casson fluid across a porous stretching/shrinking sheet: Analysis of dual solutions. Case Stud. Therm. Eng. 2021, 28, 101594. [Google Scholar] [CrossRef]
- Alshehri, N.A.; Abidi, A.; Khan, M.R.; Reddy, Y.D.; Rasheed, S.; Alali, E.; Galal, A.M. Unsteady Convective MHD Flow and Heat Transfer of a Viscous Nanofluid across a Porous Stretching/Shrinking Surface: Existence of Multiple Solutions. Crystals 2021, 11, 1359. [Google Scholar] [CrossRef]
- Wu, J.; Aidun, C.K. Simulating 3D deformable particle suspensions using lattice Boltzmann method with discrete external boundary force. Int. J. Numer. Methods Fluids 2010, 62, 765–783. [Google Scholar] [CrossRef]
- Jiang, M.; Liu, Z. A boundary thickening-based direct forcing immersed boundary method for fully resolved simulation of particle-laden flows. J. Comput. Phys. 2019, 390, 203–231. [Google Scholar] [CrossRef]
- Fo, B.; Xu, R.; Xi, J.; Lu, Y.; Song, X.; Cai, J.; Gu, Z. Numerical simulation of fine particle liquid—Solid flow in porous media based on LBM-IBM-DEM. Can. J. Chem. Eng. 2022, 101, 3576–3591. [Google Scholar] [CrossRef]
- Guo, Z.; Shi, B.; Wang, N. Lattice BGK model for incompressible Navier-Stokes equation. J. Comput. Phys. 2000, 165, 288–306. [Google Scholar] [CrossRef]
- Zou, Q.; Hou, S.; Chen, S.; Doolen, G.D. An improved incompressible lattice Boltzmann model for time-independent flow. J. Stat. Phys. 1995, 81, 35–48. [Google Scholar] [CrossRef]
- Lin, Z.; Fang, H.; Tao, R. Improved lattice Boltzmann model for incompressible two-dimensional steady flow. Phys. Rev. E 1996, 54, 6323–6330. [Google Scholar] [CrossRef]
- He, X.; Luo, L.S. Lattice Boltzmann model for the incompressible Navier-Stokes equation. J. Stat. Phys. 1997, 88, 927–944. [Google Scholar] [CrossRef]
- Glowinski, R.; Pan, T.-W.; Hesla, T.; Joseph, D. A distributed Lagrange multiplier/fictitious domain method for particulate flows. Int. J. Multiph. Flow 1999, 25, 755–794. [Google Scholar] [CrossRef]
Physical Property Parameter | Values |
---|---|
Normal elastic coefficient of particles (g/s2) | 8 × 105 |
Coefficient of restitution of particles | 0.9 |
Friction coefficient of particles | 0.3 |
Time step (s) | 0.0000012963 |
Parameter | Value |
---|---|
Density of fluid (g/cm3) | 1.0 |
Inlet velocity of fluid Vin (cm/s) | −2.0 |
Length in X-axis L (cm) | 2 |
Height in Y-axis W (cm) | 4 |
Diameter of large particles (cm) | 0.16 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Zhou, Y.; Fo, B.; Xu, R.; Xi, J.; Cai, J. Numerical Study of the Movement of Single Fine Particles in Porous Media Based on LBM-DEM. Sustainability 2024, 16, 7346. https://doi.org/10.3390/su16177346
Zhou Y, Fo B, Xu R, Xi J, Cai J. Numerical Study of the Movement of Single Fine Particles in Porous Media Based on LBM-DEM. Sustainability. 2024; 16(17):7346. https://doi.org/10.3390/su16177346
Chicago/Turabian StyleZhou, Yinggui, Bin Fo, Ruifu Xu, Jianfei Xi, and Jie Cai. 2024. "Numerical Study of the Movement of Single Fine Particles in Porous Media Based on LBM-DEM" Sustainability 16, no. 17: 7346. https://doi.org/10.3390/su16177346
APA StyleZhou, Y., Fo, B., Xu, R., Xi, J., & Cai, J. (2024). Numerical Study of the Movement of Single Fine Particles in Porous Media Based on LBM-DEM. Sustainability, 16(17), 7346. https://doi.org/10.3390/su16177346