A Stochastic Filling and Modeling Algorithm of Non-Equal Diameter Particles with Specified Probability Density for Porous Permeable Materials
Abstract
:1. Introduction
2. Generation Algorithm for Meso-Distribution Particles
2.1. Compactness Algorithm for Particles
2.1.1. Formation of Non-Equal Diameter Particles
2.1.2. Determine Initial Position
Algorithm 1 Algorithm of initial position determination |
for k in [1, fov*Nz] z_max = −∞ for j in [1, fov*Ny] y_max = −∞ prev_x_max = −∞ for i in [1, fov*Nx] generate particle with random size r if i = 1 then x_min = r else x_min = prev_x_max + r if j = 1 then y_min = r else y_min = prev_y_max + r if k = 1 then z_min = r else z_min = prev_z_max + r x_new = u(x_min, x_min+init_noise) y_new = u(y_min, y_min+init_noise) y_new = u(z_min, z_min+init_noise) place the newly generated particle at (x_new, y_new, z_new) prev_x_max = max(prev_x_max, x_new+r) y_max = max(y_max, y_new+r) z_max = max(z_max, z_new+r) prev_y_max = y_max prev_z_max = z_max |
2.1.3. Compactness of Particles
Algorithm 2 Compactness algorithm |
fork in [0, Kmax−1] for j in [0, Jmax−1] for i in [0, Imax−1] xboundary = max (x|x is the x coordinate of points on the particles at M [i−1,j,k]) yboundary = max (y|y is the y coordinate of points on the particles at M [i,jn,k], where jn in [j−nneighbor, j+nneighbor]) zboundary = max (z|z is the z coordinate of points on the particles at M [i,jn,kn], where jn in [j−nneighbor, j+nneighbor] and kn in [k−nneighbor, k+nneighbor] Assume: Plane Px parallel to the yz-plane intersects the x-axis at xboundary. Plane Py parallel to the xz-plane intersects the y-axis at yboundary. Plane Pz parallel to the xy-plane intersects the z-axis at zboundary. dx = distance between M [i,j,k] and Px dy = distance between M [i,j,k] and Py dz = distance between M [i,j,k] and Pz while max (dx, dy, dz) ⩾ if dx = max (dx, dy, dz) subtract dx/3 from the x coordinate of M [i,j,k] else if dy = max (dx, dy, dz) subtract dy/3 from the y coordinate of M [i,j,k] else dz = max (dx, dy, dz) subtract dz/3 from the z coordinate of M [i,j,k] Recalculate dx,dy,dz if loop has executed more than 1000 times: Report exception and quit |
2.2. Filtering Algorithm
Algorithm 3 Filtering algorithm |
fork in [0, Kmax−1] for j in [0, Jmax−1] for I in [0, Imax−1] Let x,y,z be the coordinates of the object M [i,j,k] if s > v2/vg∙∙fx,y,z Delete M [i,j,k] |
3. Influence of Parameters on the Algorithm
3.1. Distribution of Particles with Different Size Ranges with Same Probability Density
3.2. Particles Distribution with Same Probability Density and Different Volume Fraction
4. Algorithm Analysis
4.1. Computational Efficiency
4.2. The Randomness of Generated Samples
4.2.1. Random Test of Particles Size
4.2.2. Random Test of Particle Position
5. Permeability of Porous Media
5.1. Lattice Boltzmann Method
5.2. Darcy’s Law
5.3. Lattice Boltzmann Method Simulation Procedures
- (1)
- Setting of initial conditions;
- (2)
- Execute collision at time ;
- (3)
- Boundary processing;
- (4)
- Calculate macroscopic quantities;
- (5)
- Check whether convergence exists. If not, return to Step 2. Otherwise, go to the next step;
- (6)
- Output the result.
5.4. Simulation Results
5.4.1. The Influence of the Distribution Law
5.4.2. The Influence of the Numbers of Particles and Grids
6. Conclusions
- The filling process avoids the huge calculation burden that is caused by the continuous iterative interference detection;
- In the process of generating the initial position of the particles or realizing the compact stacking, the changes of each particle are only related to its small neighborhood, so it has a high compactness efficiency;
- The computational complexity of the algorithm is first order, while that of the traditional algorithm is much larger than second order, which illustrates the computational efficiency of the proposed algorithm;
- With this algorithm, the size and position of the particles can be distributed according to the arbitrary probability density based on the requirements, and the specified volume fraction can also be realized according to the requirements.
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
Symbols and Abbreviations
probability density | |
3-dimensional array | |
a small positive value | |
a small positive integer | |
volume of the cube | |
probability density | |
total volume of the particles | |
volume fraction | |
particle number | |
number of repeated generation position | |
a positive constant | |
number of all random numbers | |
discrete distribution function | |
velocity space | |
discrete time step | |
current time step | |
change caused by collision | |
lattice velocity | |
weight coefficient | |
fluid velocity | |
Darcy velocity | |
permeability of porous media | |
fluid pressure | |
Hamiltonian operator | |
predicted value of permeability | |
rate parameter of the exponential distribution |
References
- Noubactep, C.; Caré, S.; Togue-Kamga, F.; Schöner, A.; Woafo, P. Extending Service Life of Household Water Filters by Mixing Metallic Iron with Sand. Clean Soil Air Water 2010, 38, 951–959. [Google Scholar] [CrossRef] [Green Version]
- Wang, W.; Parker, K.H. Effect of deformable porous surface layers on the motion of a sphere in a narrow cylindrical tube. J. Fluid Mech. 1995, 283, 287–305. [Google Scholar] [CrossRef] [Green Version]
- Sobieski, W. Waterfall Algorithm as a tool of investigation the geometrical features of granular porous media. Comput. Part. Mech. 2021, 9, 551–567. [Google Scholar] [CrossRef]
- Primera, J.; Hasmy, A.; Woignier, T. Numerical study of pore sizes distribution in gels. J. Sol-Gel Sci. Technol. 2003, 26, 671–675. [Google Scholar] [CrossRef]
- Tory, E.M.; Cochrane, N.A.; Waddell, S.R. Anisotropy in Simulated Random Packing of Equal Spheres. Nature 1986, 220, 1023–1024. [Google Scholar] [CrossRef]
- Jodrey, W.S.; Tory, E.M. Computer Simulation of Isotropic, Homogeneous, Dense Random Packing of Equal Spheres. Powder Technol. 1981, 30, 111–118. [Google Scholar] [CrossRef]
- Cundall, P.; Strack, O. Discussion: A discrete numerical model for granular assemblies. Geotechnique 1980, 30, 331–336. [Google Scholar] [CrossRef] [Green Version]
- Bezrukov, A.; Bargiel, M.; Stoyan, D. Statistical analysis of simulated random packings of spheres. Part. Part. Syst. Charact. Meas. Descr. Part. Prop. Behav. Powders Other Disperse Syst. 2002, 19, 111–118. [Google Scholar] [CrossRef]
- Bailakanavar, M.; Liu, Y.; Fish, J.; Zheng, Y. Automated modeling of random inclusion composites. Eng. Comput. 2012, 30, 609–625. [Google Scholar] [CrossRef]
- Xin, Z.; Miao, W.; Wang, Y.; Chen, H. Simulation of Tensile Fracture of ZrO2 Toughened Al2O3 Particles Reinforced Fe45 Composites by Finite-Discrete Element Coupling Method. Acta Mater. Sin. 2019, 36, 1471–1479. (In Chinese) [Google Scholar]
- Stevenl, B.; Peter, R.K.; David, M. Network Model Evaluation of Permeability and Spatial Correlation in a Real Random Sphere Packing. Transp. Porous Media 1993, 11, 53–70. [Google Scholar]
- Heijs, A.W.J. Numerical evaluation of the permeability and the Kozeny constant for two types of porous media. Phys. Rev. E 1995, 51, 4346–4353. [Google Scholar] [CrossRef] [PubMed]
- Bryant, S.; Blunt, M. Prediction of relative permeability in simple porous media. Phys. Rev. A 1992, 46, 2004–2011. [Google Scholar] [CrossRef] [PubMed]
- Raats, P.A.C. Dynamics of Fluids in Porous Media. Soil Sci. Soc. Am. J. 1973, 37, vi. [Google Scholar] [CrossRef]
- Zhu, J.; Xi, Z.; Tang, H.; Tan, P. Characterization of porous structures and fractal theory. Rare Met. Mater. Eng. 2006, 35, 452–456. [Google Scholar]
- Hosmane, B.S. Improved likelihood ratio tests and Pearson chi-squared tests for independence in two dimensional contingency tables. Commun. Stat. 1986, 16, 1875–1888. [Google Scholar] [CrossRef]
- Kline, R.; Kline, R.B.; Kline, R. Principles and Practice of Structural Equation Modelling. J. Am. Stat. Assoc. 2005, 4, 1941–1947. [Google Scholar]
- Drezner, Z.; Zerom, O.T.D. A Modified Kolmogorov-Smirnov Test for Normality. Commun. Stat.–Simul. Comput. 2010, 39, 693–704. [Google Scholar] [CrossRef] [Green Version]
- Lilliefors, H.W. On The Kolmogorov-Smirnov Test for Normality with Mean and Variance Unknown. J. Am. Stat. Assoc. 1967, 62, 399–402. [Google Scholar] [CrossRef]
- Bhatnagar, P.L.; Gross, E.P.; Krook, M. A model for collision processes in gases: I. Small amplitude processes in charged and neutral one-component systems. Phys. Rev. 1953, 94, 511–525. [Google Scholar] [CrossRef]
- Qian, Y.H.; D‘Humières, D.; Lallemand, P. Lattice BGK Models for Navier-Stokes Equation. EPL (Europhys. Lett.) 1991, 17, 479–484. [Google Scholar] [CrossRef]
- Maier, R.S.; Kroll, D.M.; Kutsovsky, Y.E.; Davis, H.T.; Bernard, R.S. Simulation of flow through bead packs using the lattice Boltzmann method. Phys. Fluids 1998, 10, 60–74. [Google Scholar] [CrossRef]
- Neumann, M.; Stenzel, O.; Willot, F.; Holzer, L.; Schmidt, V. Quantifying the influence of microstructure on effective conductivity and permeability: Virtual materials testing. Int. J. Solid Struct. 2020, 184, 211–220. [Google Scholar] [CrossRef]
Range of Particle Sizes | The Distribution of Particles |
---|---|
R = 0.05–0.15 mm | |
R = 0.15–0.25 mm | |
R = 0.25–0.35 mm |
Volume Fraction | The Distribution of Particles | |
---|---|---|
10% | 12.550 | |
20% | 25.504 | |
30% | 37.743 |
Test Statistics | |
---|---|
Statistic | Value |
11.720 | |
Degrees of freedom | 9 |
Asymptotic significance | 0.230 |
Sample Size | Average | Standard Deviation | Partial Degrees | Kurtosis | Kolmogorov–Smirnov Test | Shapiro–Wilk Test | ||
---|---|---|---|---|---|---|---|---|
The Statistic D Value | p | The Statistic W Value | p | |||||
301 | 3.014 | 0.511 | −0.028 | −0.219 | 0.029 | 0.742 | 0.996 | 0.610 |
Distribution | Average Velocity | |
---|---|---|
Uniform | 0.586 × 10−3 | 5.204 × 10−10 |
Normal | 0.533 × 10−3 | 5.335 × 10−10 |
Exponential | 0.681 × 10−3 | 6.757 × 10−10 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 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
Zhang, W.; He, L.; Wang, F.; Zhang, G. A Stochastic Filling and Modeling Algorithm of Non-Equal Diameter Particles with Specified Probability Density for Porous Permeable Materials. Materials 2022, 15, 4733. https://doi.org/10.3390/ma15144733
Zhang W, He L, Wang F, Zhang G. A Stochastic Filling and Modeling Algorithm of Non-Equal Diameter Particles with Specified Probability Density for Porous Permeable Materials. Materials. 2022; 15(14):4733. https://doi.org/10.3390/ma15144733
Chicago/Turabian StyleZhang, Wei, Lile He, Fazhan Wang, and Guangyong Zhang. 2022. "A Stochastic Filling and Modeling Algorithm of Non-Equal Diameter Particles with Specified Probability Density for Porous Permeable Materials" Materials 15, no. 14: 4733. https://doi.org/10.3390/ma15144733
APA StyleZhang, W., He, L., Wang, F., & Zhang, G. (2022). A Stochastic Filling and Modeling Algorithm of Non-Equal Diameter Particles with Specified Probability Density for Porous Permeable Materials. Materials, 15(14), 4733. https://doi.org/10.3390/ma15144733