Next Article in Journal
Applicability of Composite Magnetic Membranes in Separation Processes of Gaseous and Liquid Mixtures—A Review
Previous Article in Journal
Experimental Investigation of Solar-Driven Hollow Fiber Membrane Liquid Dehumidification System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of CFD-DEM on the Impact of the Rolling Friction Coefficient on Deposition of Lignin Particles in a Single Ceramic Membrane Pore

1
College of Electromechanical Engineering, Qingdao University of Science and Technology, Qingdao 266061, China
2
Dongyue Group, Zibo 256401, China
*
Authors to whom correspondence should be addressed.
Membranes 2023, 13(4), 382; https://doi.org/10.3390/membranes13040382
Submission received: 6 January 2023 / Revised: 19 February 2023 / Accepted: 25 March 2023 / Published: 27 March 2023
(This article belongs to the Section Inorganic Membranes)

Abstract

:
The discrete element method coupled with the computational fluid dynamic (CFD-DEM) method is effective for studying the micro-flow process of lignin particles in ceramic membranes. Lignin particles may exhibit various shapes in industry, so it is difficult to model their real shapes in CFD-DEM coupled solutions. Meanwhile, the solution of non-spherical particles requires a very small time-step, which significantly lowers the computational efficiency. Based on this, we proposed a method to simplify the shape of lignin particles into spheres. However, the rolling friction coefficient during the replacement was hard to be obtained. Therefore, the CFD-DEM method was employed to simulate the deposition of lignin particles on a ceramic membrane. Impacts of the rolling friction coefficient on the deposition morphology of the lignin particles were analyzed. The coordination number and porosity of the lignin particles after deposition were calculated, based on which the rolling friction coefficient was calibrated. The results indicated that the deposition morphology, coordination number, and porosity of the lignin particles can be significantly affected by the rolling friction coefficient and slightly influenced by that between the lignin particles and membranes. When the rolling friction coefficient among different particles increased from 0.1 to 3.0, the average coordination number decreased from 3.96 to 2.73, and the porosity increased from 0.65 to 0.73. Besides, when the rolling friction coefficient among the lignin particles was set to 0.6–2.4, the spherical lignin particles could replace the non-spherical particles.

1. Introduction

With its rapid development, the papermaking industry has become a pillar in various countries [1]. However, it discharges a large amount of pulping waste liquor, resulting in environmental pollution and waste of resources. Pulping waste liquor is produced in the pulping and cooking stage and is the main pollution source of the papermaking industry. Due to its dark brown color, high viscosity, and stench, it is called black liquor [2,3,4,5]. As a representative organic matter in black liquor, lignin is not only a biomass energy with great utilization value, but also the only non-petroleum resource in nature that can provide renewable aryl compounds [6]. In the current context of a forest-based circular bioeconomy, lignin is of considerable commercial value, accounting for approximately 24–47% of the current pulp production revenue [7]. Therefore, an efficient and low-loss separation of lignin from black liquor is of great significance to alleviate environmental pollution and improve waste recycling.
The commonly used lignin separation methods include acid precipitation and membrane separation methods [8]. The acid precipitation method is characterized by a simple process and a low cost, but it requires a quantity of acid, which easily causes secondary pollution. Currently, the membrane separation technologies have made a figure in lignin separation due to its low pollution and high efficiency. Ceramic membranes present excellent thermal, chemical, and mechanical stabilities, long service life, high separation efficiency, and high mechanical strength [9,10,11]. Therefore, they have been well accepted as one of the most effective materials for lignin separation [12]. However, fouling in the filtration has become a key factor hindering the continuous industrial production of ceramic membranes, greatly reducing the lignin separation from black liquor. Therefore, the micro filtration mechanisms have attracted increasing attention, such as the movement of lignin particles, sedimentation of lignin particles on the membrane surface, and formation of filter cake. Although the visualization of micropore blockage has been rapidly developed [13], dynamic tracking of the movement and capture of lignin particles can’t be achieved due to the limitations of current experimental conditions and observation techniques. Meanwhile, the size of lignin particles is generally micron-level, and the experimental results cannot provide microscopic information on the interaction among lignin particles and that between lignin particles and fluid. Therefore, the discrete element method coupled with the computational fluid dynamic (CFD-DEM) has been proposed to investigate the dynamic characteristics of lignin particles and the formation mechanism of filter cake during filtration of ceramic membranes [14,15,16,17].
The CFD-DEM solution is a relatively new method that can dynamically characterize the particle trajectory, deposition morphology, interaction between particles and fluid, and collision among particles during filtration [18]. Therefore, it is extensively applied in the filtration and separation field. Deshpande et al. [19] adopted the DEM and CFD-DEM methods to investigate the precipitation and consolidation of packed bed/filter cake formed by monodisperse and bidisperse spherical particles under different flow conditions. Li and Marshall [20] simulated the precipitation of particles on a single fiber by the CFD-DEM method, and then discussed the impacts of adhesion on particle deposition. Hund et al. [21] employed the CFD-DEM method to simulate the formation of particle bridges in solid-liquid separation, studied the impacts of particle concentration and feed flow rate on the formation of particle bridges, preliminarily revealed the formation mechanism of particle bridges, and verified the accuracy of the simulation results by comparing them with the experimental results. Additionally, Qian et al. [22,23] and Cao et al. [24] established a multi-layer three-dimensional fiber medium filtration model and simulated the flow and precipitation of particles in the model by using the CFD-DEM method. The results revealed that it was convenient and feasible to study the flow and precipitation of small particles in filter media by the CFD-DEM method, which were consistent with the experimental results. However, in these studies, the shape of the particles was specially processed and simplified to be spherical. Most of the particles applied in industrial production are non-spherical, and the shape greatly influences the dynamic behavior of the particles [25]. The ellipsoid models [26], hyper ellipsoid models [27], bonded sphere models [28], multi-sphere models [29], and polyhedron models [30] have been developed in DEM to describe the shape of particles more accurately to achieve a more realistic simulation. The most important difference between spherical particles and non-spherical particles is that the formers are prone to rolling under the same conditions when the CFD-DEM method is adopted [25]. However, with more accurate descriptions of the shape of the particles, the amount of calculation in simulation is increasing, especially in solving the movement of micron-sized particles. Rolling motion of the spherical particles is controlled by the rolling friction coefficient among particles and that between the particles and the membrane, so the coefficient can be introduced into the DEM to generate the rolling resistance. Therefore, spherical particles can replace the non-spherical particles by controlling the numerical range of the rolling friction coefficient. In this context, people began to replace non-spherical particles with spherical particles, and the rolling friction coefficient has been artificially increased to show some properties of non-spherical particles [31], aiming to satisfy the simulation requirements and reduce the amount of calculation.
There are many studies on the range of the rolling friction coefficient for spherical particles to replace non-spherical particles. Wensrich and Katterfeld [32] revealed that the rolling friction coefficient could indeed bring some properties of non-spherical particles to spherical particles. Xie et al. [25] proposed a method of replacing non-spherical particles with spherical particles based on the rolling friction and then studied the applicability of this method under different processes and operating conditions. The results indicated that the non-spherical particles can be replaced, to some extent, to reduce the amount of calculation and speed up the solution if a reasonable rolling friction coefficient can be applied. Xiong et al. [31] simulated particle deposition on a single fiber, studied the impacts of the rolling friction coefficient on the deposition morphology, coordination number, and porosity of particles, and then calibrated the rolling friction coefficient among particles and that between particles and fibers in fiber filtration, making the simulation results closer to the real observations. In general, the rolling friction coefficient significantly affects the movement of spherical particles. Therefore, correctly defining the rolling friction coefficient is the key to ensure that the simulation process and results are in line with reality.
Many researchers have noticed that spherical particles can be substituted for non-spherical particles by setting a proper rolling friction coefficient. However, calibration of the rolling friction coefficient is always the key- and difficult-point during the simplification. In this study, non-spherical lignin particles were simplified to spherical particles by applying the rolling friction coefficient to the latter. Effects of the rolling friction coefficient on the deposition of lignin particles in ceramic membranes were studied by the CFD-DEM coupling method at the single pore level. Firstly, the capture of lignin particles was demonstrated dynamically. Then, effects of the rolling friction coefficient on the morphology, coordination number distribution, and porosity of the lignin particles were investigated and analyzed. Finally, the rolling friction coefficient was calibrated. Introducing the correct rolling friction coefficient could simplify the DEM model and improve its accuracy, providing a more reasonable rolling friction coefficient for subsequent research on the more complex deposition of lignin particles in porous ceramic membranes.

2. Governing Equation of Fluid–Solid Two-Phase Flow

The coupling solution process based on CFD-DEM was completed under the Euler-Lagrange coupling framework, in which the Navier-Stokes (N-S) control equation was adopted to solve the fluid movement based on the continuum assumption. In addition, the Newton’s Second Law was adopted to solve the movement of each particle, and the coupling between the two is realized by the Newton’s Third Law. The flowchart of coupling is shown in Figure 1 [16].

2.1. Fluid Phase Control Equation

In the coupling solution process, the particle phase greatly affected the movement characteristics of the fluid. Therefore, the void fraction was necessary to be incorporated into the traditional N-S equation to characterize the volume of the fluid phase in a specific computational grid. The control equation of the fluid phase can be expressed as follows [33,34]:
ρ φ t + ρ φ u = 0
t ρ φ u + ρ φ ϑ u = ρ S + ϑ φ u + ρ φ g
S = i n F D / V
where, ϑ , u , and ρ refer to the viscosity, velocity, and density of the fluid with the units of Pa·s, m/s, and kg/m3, respectively; S is the sum of the fluid resistance F D acting on the volume of a grid cell, N; and V is the volume of the mesh unit, m3.

2.2. Discrete Model

During the movement of particles, the Newton’s Second Law was employed to solve the velocity and position of particles at each specific time based on the external force on particles. In this study, the external forces on particles included gravity, the interaction force between particles and fluid, the collision force among particles and that between particles and membranes, and the Van der Waals adhesion force among particles and that between particles and membranes. Therefore, the particle-phase control equation can be expressed as follows [14]:
m p d u p d t = m p g + F J K R + F f p + F c
I p d ω p d t = i = 1 k M i + M J K R , i
where, m p , u p , I p , ω p , and M i , are mass, velocity, inertia moment, rotational speed, and collision torque of the particle, respectively, with the units of kg, m/s, kg·m2, rad/s, and N·m, respectively. M J K R , i is the adhesive torque in collisions with other particles, N·m; F J K R refers to the Van der Waals adhesion force, N; F f p is the particle-fluid force, N; and F c denotes the particle–particle contact force, N.
F J K R is the Van der Waals adhesion, which can be calculated by Equation (6) [35,36]:
F J K R = 4 π γ E * α 3 / 2 + 4 E * 3 R * α 3
where, γ is the surface energy, J/m2; E * refers to the relative Young’s Modulus, Pa; α is the normal overlap, m; and R * represents the relative radius, m.
The equivalent Young’s Modulus E * and relative radius R * are defined as follows:
1 E * = 1 υ i 2 E i + 1 υ j 2 E j
1 R * = 1 R i + 1 R j
where, E i and E j are the Young’s Modulus of particle-i and particle-j, respectively, Pa; υ i and υ j are the Passion’s ratios of particle-i and particle-j, respectively; and R i and R j are radii of particle-i and particle-j, respectively, m.
F f p represents the force between the particle and fluid. In this study, the interaction force between liquid and solid phases only considers the drag force F d r a g and pressure gradient force F p , so F f p can be expressed as Equation (9).
F f p = F d r a g + F p
F d r a g represents the fluid drag force of the particle, as shown in Equation (10).
F d r a g = 0.5 C D ρ A u u p u u p
where, F d r a g is the drag force of the particle, N; A is the projection area of the particle, m2; C D refers to the drag coefficient, which depends on the Reynolds number R e , and can be calculated by Equation (11) [37].
C D = 24 R e                                                                     R e 0.5 24 1.0 + 0.25 R e 0.687 R e         0.5 < R e 1000 0.44                                                                     R e > 1000 ,   R e = φ ρ d p u u p μ
F p is the pressure gradient force caused by the pressure gradient of particles moving in the flow field, as shown in Equation (12) [38].
F p = d p / d x = ρ g ρ u d u / d x
F c is the particle–particle contact force, as shown in Equation (13) [39].
F c = F c n + F c t
where, F c n is the normal contact force acting on particles after particle collision, N, as expressed as Equation (14); and F c t is the tangential contact force acting on particles after particle collision, N, as written as Equation (15).
F c n , i j = k n α 3 / 2 c n u i j · n n
F c t , i j = k t δ c t v c t
where, u i j = u i u j refers to the velocity of particle-i relative to particle-j, m/s; n is the unit vector pointing from the centroid of particle-i to the centroid of particle-j; δ represents the tangential displacement of contact point, m; and v c t denotes the sliding velocity vector.

2.3. The Directional Constant Torque Model

As mentioned in the first section, the most important difference between spherical and non-spherical particles during the movement is that the formers are easy to roll under the same conditions. Therefore, the principle of simplifying non-spherical particles into spherical particles can be described as follows. The rolling resistance is applied to the simplified spherical particles, so that they show the movement characteristics of non-spherical particles. In this case, they can replace the non-spherical particles. However, the rolling resistance is difficult to be determined in various industries. In the DEM, the directional constant torque model is the most widely accepted rolling resistance model, which can be calculated with Equation (16).
M r = μ F c n R c ω c
where, M r is the rolling friction resistance, N; μ represents the rolling friction coefficient; R c refers to the distance between the contact point and the center of the sphere, m; and ω c is the angular velocity of the contact point, rad/s. Therefore, the determination of the rolling resistance can be converted to the determination of the rolling friction coefficient. Changing the rolling friction coefficient can alter the size of the rolling resistance, and the rolling friction coefficient can be calibrated through studying a series of physical quantities. In this way, the non-spherical particles can be replaced with the spherical particles.

3. Computational Set-Up

3.1. Geometry and Computational Domain

According to the microscopic characterizations of the structure and morphology, the filter channel in the ceramic membrane was formed by the accumulation of ceramic particles with a shape similar to a sphere. Such a porous structure for a ceramic membrane was difficult to model, and many membrane pores would negatively affect the ability to characterize the deposition morphology of particles in each membrane pore, making it difficult to calibrate the rolling friction coefficient. Therefore, the porous structure of the ceramic membrane was simplified in this study by analogy with the method introduced by Tao et al. [17] and Xiong et al. [31]. Figure 2 illustrated the simplified pore structure of a single ceramic membrane and the size of the computational domain. Among them, the X and Y directions of the computational domain size were twice the diameter of the ceramic particle, which exhibits similar results with those obtained by Xiong et al. [31] in the single fiber filtration model. The computational domain size along the direction of fluid flow is set to six times the diameter of the ceramic particle, according to the reports of Xiong et al. [31] and Qian et al. [23], to ensure the uniform flow of the black liquor and lignin particles in the inlet and outlet. Moreover, the shape of the ceramic particles was simplified accordingly. The ceramic particle spheres, similar to spheres, were simplified as spheres, and the rolling friction coefficient between the lignin particles and ceramic membranes was introduced in the solution setting. Calibrating the rolling friction coefficient between the lignin particles and ceramic membranes can make the physical model closer to reality and greatly accelerate the calculation.

3.2. Boundary Conditions and Parameter Settings

In this study, the EDEM 2020 and Fluent 2020R2 were coupled to calculate the particles’ movement and fluid fields through User-Defined Function (UDF). The coupling between CFD and DEM is implemented using the Eulerian coupling method, which considers the momentum exchange between liquid phase and particle phase, and the influence of particles relative to the liquid phase. Particles in the DEM are translated and rotated by the explicit time integral method, while the governing equation of fluid in the CFD is solved by the SIMPLE algorithm under the pressure base solver. The second order scheme is employed to disperse the pressure term, and the first order scheme discretizes the other terms. Meanwhile, the drag model of Ergun and Wen & Yu was selected. In addition, the black liquor adopted the speed import and pressure export. The lignin particles entered the flow field under the driving force of the black liquor. To ensure that the lignin particles entered the flow fields smoothly, the lignin particle generating surface was behind the velocity entrance. In consideration of the structural characteristics of the ceramic membrane, surfaces of the ceramic particles were set as the non-slip boundary condition, and the other boundaries of the computational domain were set as the symmetric boundary condition. The total calculation time was 2.2 ms, and the lignin particles were generated only within 0–2 ms, which ensured the complete deposition or flowing out of all the lignin particles in the computational domain at the end of the simulation. The main parameters used for simulation in this study are shown in Table 1 and Table 2.
A semi-analytical CFD-DEM coupling interface was employed in this study. Instead of precisely analyzing the flow state around each particle, the semi-analytical interface redistributed the particle-phase volume spatially by introducing kernel functions and capturing the physical information of the background flow fields in the expanding region. The physical information in the background flow fields could be obtained reasonably and accurately no matter how many grids were covered. Therefore, the semi-analytical interface could simulate the movement of particles in a fluid grid that is equivalent to or even slightly smaller than the size of particles. It solved the problem that the unresolved CFD-DEM coupling interface cannot be applied to the field of ceramic membrane filtration, and solved the problem that the resolved CFD-DEM coupling interface had a large amount of calculation. Therefore, the size of the fluid domain grid was set within a range equal to or slightly smaller than the diameter of the lignin particle in this study. In this case, the grid independence was verified based on the clean filtration stage, as displayed in Table 3. The semi-analytical CFD-DEM coupling interface exhibited good stability, and the simulation results could be correctly represented when the number of grids was 46,400. Figure 3 demonstrated the fluid domain grid model when the number of grids was 46,400. As it illustrated, the fluid domain was meshed by ICEM CFD, and the hexahedral mesh was generated by the block topology function under the ICEM CFD. The minimum and maximum mesh sizes were 0.2 μm and 0.7 μm, respectively, and the maximum mesh mass was 1. It should be noted that the grid quality here was measured by the determinant with a value range of 0–1. The more the value was close to 1, the more perfect the grid was.

4. Results and Discussion

4.1. Particle Deposition Process

Figure 4 illustrates the deposition process of the lignin particles with the rolling friction coefficient among the lignin particles of μ p p = 1.2 , and that between the lignin particles and membranes of μ p m = 1.0 . The deposition process of the lignin particles could be roughly divided into two stages: the capture of the ceramic membrane (A–B), and the capture of the deposited lignin particles (B–F). In the initial stage of filtration (A–B), the lignin particles were mainly captured and approximately evenly distributed on the upstream side surface of the ceramic membrane. At this stage, the number of lignin particles captured was not large, and the pressure drop increased slightly. As indicated by the pressure drop curve in the A–B stage (Figure 5), the pressure drop increased linearly but the growth rate was relatively slow. As the filtration process proceeded (B–F), the lignin particles began to be captured by the lignin particles deposited, leading to a significant change in the deposition morphology. The deposited lignin particles began to grow outward, forming a dendritic structure, which became more obvious and independent as the filtration process progressed. Finally, the deposition morphology of the lignin particles showed the shape of a “forest.” At this stage, most of the lignin particles were captured by the dendritic structure, greatly increasing the number of captured particles and the pressure drop. As demonstrated in the pressure drop curve in the B–F stage in Figure 5, the slope of the pressure drop curve with time gradually increased, indicating that the efficiency in capturing lignin particles was constantly increasing. However, it should be noted that the lignin particles were only generated within 0–2 ms, so the pressure drop decreased gradually when the filtration time was longer than 2 ms (E–F stage in Figure 5).
The cumulative number of penetrating particles and the efficiency in capturing lignin particles during filtration were counted and calculated to further investigate the impacts of dendritic structure growth on the efficiency in capturing lignin particles. As demonstrated in Figure 6, the cumulative number of penetrating particles and the efficiency in capturing lignin particles both increased as the filtration time increased. However, the cumulative number of penetrated particles grew slowly, and the efficiency in capturing lignin particles exhibited a faster growth rate. This is because the deposited lignin particles formed dendritic structures, exerting a secondary trapping effect on the lignin particles, so that an increasing number of lignin particles were captured, and the retention rate was increased. In other words, formation of the dendritic structure significantly affected the efficiency in capturing lignin particles.

4.2. Impacts of the Rolling Friction Coefficient on Deposition Morphology of Lignin Particles

The actual lignin particles had greater rolling friction resistance than the ideal spherical lignin particles due to their shapes, thereby limiting the movement among lignin particles, and that between lignin particles and membranes [31]. Previous studies have shown that the deposition morphology of lignin particles was correlated with the external force. When the rolling friction resistance was not enough to resist the external force, the lignin particles would move until they reached a new equilibrium state [40], which significantly affected the deposition morphology of the lignin particles. Therefore, effects of the rolling friction coefficient on the deposition morphology of lignin particles were quantitatively characterized when the rolling friction coefficient among the lignin particles was μ p p = 0 3 , and that between the lignin particles and membranes was μ p m = 0 3 .
Figure 7 shows the deposition morphology of lignin particles under different rolling friction coefficients. It revealed that regardless of μ p m , the deposition morphology of lignin particles showed a strong regularity with increasing μ p p . When μ p p 0.6 , the particles were deposited on the surface of the ceramic membrane in the form of agglomeration, and no dendritic structure was observed, which was inconsistent with the experimental results of Huang et al. [41]. The reason is that the lignin particles were easier to roll because of a too-small rolling resistance, and the dendritic structure formed at the early stage lodged down to the ceramic membrane surface so that the deposition morphology of lignin particles exhibited an agglomeration. As the rolling friction coefficient among lignin particles increased, the corresponding rolling resistance increased, and the lignin particles were difficult to rotate around the contact point. At this moment, an obvious dendritic structure was observed. As there was an increase in the rolling friction coefficient among the lignin particles, the dendritic structure became increasingly obvious and independent, making the deposition morphology of the lignin particles look like a “forest.” This finding is consistent with the deposition morphology of the particles on the single fiber observed by Huang et al. [41].
Moreover, Figure 7 discloses that the rolling friction coefficient between the lignin particles and membranes had no obvious effect on the deposition morphology of the lignin particles. This might be because the number of contacts among the lignin particles was much larger than that between the lignin particles and membranes [31]. Table 4 showed the proportion of the number of contacts between the lignin particles and membranes to the total number of contacts under different rolling friction coefficients. The number of contacts between the lignin particles and membranes accounted for only 18.4% at most, weakening the effects of the rolling friction coefficient between them on the deposition morphology of lignin particles.

4.3. Impacts of the Rolling Friction Coefficient on the Deposition Structure of Lignin Particles

Generally, coordination number (Q) and porosity ( ε ) are two important parameters describing the structural properties of particle packing [31]. The former refers to the number of particles in contact with the central particles, and its size reflects the agglomeration properties of the stacking structure. The latter reflects the compactness of the accumulation structure. Therefore, this section discusses the effects of the rolling friction coefficient on the average coordination number, coordination number distribution, and porosity of lignin particles.
Figure 8 displays the curve of the average coordination number of lignin particles with the rolling friction coefficient. It illustrated that the average coordination number decreased with the increasing rolling friction coefficient among the lignin particles, and the change law was consistent with the results obtained by Xiong et al. [31]. When the rolling friction coefficient among the lignin particles increased from 0.1 to 3.0, the average coordination number decreased from 3.96 to 2.73; it indicated that the number of contacts among the lignin particles gradually decreased, and the dendritic structure became more independent as the rolling friction coefficient among the lignin particles increased. Such a conclusion was consistent with the change law of the deposition morphology of lignin particles mentioned in Section 4.2. Besides, the rolling friction coefficient between the lignin particles and membranes slightly affected the average coordination number, which might be because the number of contacts between the lignin particles and membranes was too small (as described in Section 4.2). Yang et al. analyzed the stacking structure of settling particles, only under the action of gravity, and revealed that the average coordination number of particles at 1 μm was 2.13 [42]. However, the fluid force would increase the average coordination number, so the reasonable and acceptable minimum average coordination number was determined as 2.13 in this study; that is, the rolling friction coefficient in this simulation range could meet the requirements of the coordination number.
Figure 9 describes the impacts of the rolling friction coefficient among the lignin particles on the coordination number distribution. The increase in the rolling friction coefficient made the dendritic structure formed by the deposition of lignin particles more independent, decreasing the distribution range of the coordination number with an increase in the rolling friction coefficient among the lignin particles. The research conducted by Yang et al. suggested that particles at 1 μm were only subject to the action of gravity, and the distribution range of the coordination number was mainly 1–4 [42]. However, the contact among different particles became closer due to the action of fluid force, which expanded the distribution range of the coordination number.
Figure 10 exhibits the impacts of the rolling friction coefficient between the lignin particles and membranes on the coordination number distribution. It illustrates that the rolling friction coefficient between the lignin particles and membranes did not affect the coordination number distribution. Therefore, it can be concluded that the rolling friction coefficient among the lignin particles was the primary factor affecting the coordination number distribution.
The porosity of the particle packing structure can be calculated by Equation (17) [31]:
ε = 1 Q / Q 0 1 / m n 1 / 4
where, Q is the mean coordination number; Q 0 = 2.02; m = 87.38; and n = 25.81.
Figure 11 displays the variation curve of porosity with the rolling friction coefficient. As the rolling friction coefficient between the lignin particles and membranes did not affect the average coordination number, it also did not affect the porosity [31]. We observed that the porosity increased from 0.65 to 0.73 when the rolling friction coefficient among the lignin particles increased from 0.1 to 3.0. This was because the increase in the rolling friction coefficient between the lignin particles improved the rolling resistance among the lignin particles, which enhanced the interaction among the lignin particles, thus increasing the porosity. According to the experimental results of Dingwell et al., porosity of the filter cake formed by the lignin particles was 0.61–0.71 [43]. When the rolling friction coefficient among the lignin particles was 0.1–2.4, the porosity was 0.65–0.72. Therefore, it was more in line with the actual porosity when the rolling friction coefficient among the lignin particles was set to 0.1–2.4 in the simulation, according to the experimental results of Dingwell et al.

5. Conclusions

In this study, a method of simplifying the shape of lignin particles into spheres was proposed to solve the difficult CFD-DEM solution and low computational efficiency of non-spherical lignin particles. Then, the CFD-DEM method was adopted to simulate the deposition of spherical lignin particles in the pores of a single ceramic membrane. Later, impacts of the rolling friction coefficient on the deposition morphology, average coordination number, coordination number distribution, and porosity of lignin particles were investigated and analyzed. In addition, the rolling friction coefficient was calibrated. Finally, the following conclusions were drawn:
(1)
The deposition of lignin particles on ceramic membranes was dynamic, which mainly included capturing ceramic membranes in the initial filtration and deposited lignin particles. Formation of a dendritic structure not only made the deposition morphology of lignin particles look like a “forest,” but also greatly improved the efficiency in capturing the lignin particles.
(2)
The rolling friction coefficient among the lignin particles crucially affected the deposition morphology, average coordination number, coordination number distribution, and porosity of the particles; the average coordination number decreased from 3.96 to 2.73, and the porosity increased from 0.65 to 0.73, when it increased from 0.1 to 3.0.
(3)
Reasonably providing a rolling friction coefficient among the lignin particles could replace spherical lignin particles with non-spherical particles. Impacts of the rolling friction coefficient on the deposition morphology, coordination number, and porosity of lignin particles enabled the simulation to be closer to the real lignin filtration by setting the rolling friction coefficient among the lignin particles as 0.6–2.4.

Author Contributions

Conceptualization, Methodology, software, writing—original draft, H.W. and J.W.; data curation, investigation, resources, J.W. and P.F.; writing—review & editing, X.W. and Y.L.; Writing—original draft, S.W. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors give sincere gratitude to Wu Junfei from Qingdao University of Science and Technology for his continued support of research activities.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Song, P.; Ping, L.; Gao, J.; Li, X.; Zhu, M.; Wang, J. Ecotoxicological effects of fertilizers made from pulping waste liquor on earthworm Eisenia fetida. Ecotoxicol. Environ. Saf. 2018, 166, 237–241. [Google Scholar] [CrossRef] [PubMed]
  2. Sun, M.; Wang, Y.; Shi, L.; Klemeš, J. Uncovering energy use, carbon emissions and environmental burdens of pulp and paper industry: A systematic review and meta-analysis. Renew. Sust. Energ. Rev. 2018, 92, 823–833. [Google Scholar] [CrossRef]
  3. Chen, T.; Li, Y.; Lei, L.; Hong, M.; Sun, Q.; Hou, Y. Influence of residual black liquor in pulp on wastewater pollution after bleaching process. Bioresources 2017, 12, 2031–2039. [Google Scholar] [CrossRef] [Green Version]
  4. Chen, T.; Li, Y.; Lei, L.; Hong, M.; Sun, Q.; Hou, Y. The influence of stock consistency on the pollution load in washing process. Bioresources 2016, 11, 2214–2223. [Google Scholar] [CrossRef] [Green Version]
  5. Ji, X. Experimental Study and Numerical Simulation on Straw Pulp Black Liquor Combustion in Fluidized Bed for Direct Causticization. Ph.D. Thesis, Harbin Institute of Technology, Harbin, China, 2017. [Google Scholar]
  6. Zhang, Y. Study on extraction and antioxidation of lignin in papermaking black liquor. Master’s Thesis, Fudan University, Shanghai, China, 2009. [Google Scholar]
  7. Wenger, J.; Haas, V.; Stern, T. Why can we make anything from lignin except money? Towards a broader economic perspective in lignin research. Curr. For. Rep. 2020, 6, 294–308. [Google Scholar] [CrossRef]
  8. Vives, M.; Thuvander, J.; Arkell, A.; Lipnizki, F. Low-molecular-weight lignin recovery with nanofiltration in the kraft pulping process. Membranes 2022, 12, 310. [Google Scholar] [CrossRef]
  9. Singh, P.; Manikandan, N.; Purnima, M.; Pakshirajan, K.; Pugazhenthi, G. Recovery of lignin from water and methanol using low-cost kaolin based tubular ceramic membrane. J. Water Process Eng. 2020, 38, 101615. [Google Scholar] [CrossRef]
  10. Humpert, D.; Ebrahimi, M.; Stroh, A.; Czermak, P. Recovery of Lignosulfonates from Spent Sulfite Liquor Using Ceramic Hollow-Fiber Membranes. Membranes 2019, 9, 45. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, W. The Preparation of Inorganic Microporous Membrane and Application in Lignin Separation. Master’s Thesis, South China University of Technology, Guangzhou, China, 2018. [Google Scholar]
  12. Colyar, K.; Pellegrino, J.; Kadam, K. Fractionation of pre-hydrolysis products from lignocellulosic biomass by an ultrafiltration ceramic tubular membrane. Sep. Sci. Technol. 2018, 43, 447–476. [Google Scholar] [CrossRef]
  13. Agbangla, G.; Climent, É.; Bacchin, P. Experimental investigation of pore clogging by microparticles: Evidence for a critical flux density of particle yielding arches and deposits. Sep. Purif. Technol. 2012, 101, 42–48. [Google Scholar] [CrossRef] [Green Version]
  14. Zhu, H.; Zhou, Z.; Yang, R.; Yu, A. Discrete particle simulation of particulate systems: Theoretical developments. Chem. Eng. Sci. 2007, 62, 3378–3396. [Google Scholar] [CrossRef]
  15. Zhu, H.; Zhou, Z.; Yang, R.; Yu, A. Discrete particle simulation of particulate systems: A review of major applications and findings. Chem. Eng. Sci. 2008, 63, 5728–5770. [Google Scholar] [CrossRef]
  16. Wang, H.; Wu, J.; Fu, P.; Qu, Z.; Zhao, W.; Song, Y. CFD-DEM study of bridging mechanism of particles in ceramic membrane pores under surface filtration conditions. Processes 2022, 10, 475. [Google Scholar] [CrossRef]
  17. Tao, R.; Yang, M.; Li, S. Effect of adhesion on clogging of microparticles in fiber filtration by DEM-CFD simulation. Powder Technol. 2020, 360, 289–300. [Google Scholar] [CrossRef]
  18. Hu, Z.; Zeng, H.; Ge, Y.; Wang, W.; Wang, J. Simulation and experiment of gas-solid flow in a safflower sorting device based on the CFD-DEM coupling method. Processes 2021, 9, 1239. [Google Scholar] [CrossRef]
  19. Deshpande, R.; Antonyuk, S.; Iliev, O. Study of the filter cake formed due to the sedimentation of monodispersed and bidispersed particles using discrete element method–computational fluid dynamics simulations. AIChE J. 2019, 65, 1294–1303. [Google Scholar] [CrossRef]
  20. Li, S.; Marsha, J. Discrete element simulation of micro-particle deposition on a cylindrical fiber in an array. J. Aerosol Sci. 2007, 38, 1031–1046. [Google Scholar] [CrossRef]
  21. Hund, D.; Lösch, P.; Kerner, M.; Ripperger, S.; Antonyuk, S. CFD-DEM study of bridging mechanisms at the static solid-liquid surface filtration. Powder Technol. 2020, 361, 600–609. [Google Scholar] [CrossRef]
  22. Qian, F.; Huang, N.; Zhu, X.; Lu, J. Numerical study of the gas–solid flow characteristic of fibrous media based on SEM using CFD–DEM. Powder Technol. 2013, 249, 63–70. [Google Scholar] [CrossRef]
  23. Qian, F.; Huang, N.; Lu, J.; Han, Y. CFD–DEM simulation of the filtration performance for fibrous media based on the mimic structure. Comput. Chem. Eng. 2014, 71, 478–488. [Google Scholar] [CrossRef]
  24. Cao, B.; Wang, S.; Dong, W.; Zhu, J.; Qian, F.; Lu, J.; Han, Y. Investigation of the filtration performance for fibrous media: Coupling of a semi-analytical model with CFD on Voronoi-based microstructure. Sep. Purif. Technol. 2020, 251, 117364. [Google Scholar] [CrossRef]
  25. Xie, C.; Ma, H.; Zhao, Y. Investigation of modeling non-spherical particles by using spherical discrete element model with rolling friction. Eng. Anal. Bound. Elem. 2019, 105, 207–220. [Google Scholar] [CrossRef]
  26. Xu, W.; Chen, H. Mesostructural characterization of particulate composites via a contact detection algorithm of ellipsoidal particles. Powder Technol. 2012, 211, 296–305. [Google Scholar] [CrossRef]
  27. Ma, H.; Zhao, Y. Modelling of the flow of ellipsoidal particles in a horizontal rotating drum based on DEM simulation. Chem. Eng. Sci. 2017, 172, 636–651. [Google Scholar] [CrossRef]
  28. Guo, Y.; Wassgren, C.; Hancock, B.; Ketterhagen, W.; Curtis, J. Granular shear flows of flat disks and elongated rods without and with friction. Phys. Fluids 2013, 25, 063304. [Google Scholar] [CrossRef]
  29. Shamsi, M.; Mirghasemi, A. Numerical simulation of 3D semi-real-shaped granular particle assembly. Powder Technol. 2012, 221, 431–446. [Google Scholar] [CrossRef]
  30. Wachs, A.; Girolami, L.; Vinay, G.; Ferrer, G. Grains3D, a flexible DEM approach for particles of arbitrary convex shape—Part I: Numerical model and validations. Powder Technol. 2012, 224, 374–389. [Google Scholar] [CrossRef]
  31. Xiong, G.; Gao, Z.; Hong, C.; Qiu, B.; Li, S. Effect of the rolling friction coefficient on particles’ deposition morphology on single fibre. Comput. Geotech. 2020, 121, 103450. [Google Scholar] [CrossRef]
  32. Wensrich, C.; Katterfeld, K. Rolling friction as a technique for modelling particle shape in DEM. Powder Technol. 2012, 217, 409–417. [Google Scholar] [CrossRef]
  33. Chu, K.; Wang, B.; Yu, A.; Vince, A. CFD-DEM modelling of multiphase flow in dense medium cyclones. Powder Technol. 2009, 193, 235–247. [Google Scholar] [CrossRef]
  34. Zhou, Z.; Kuang, S.; Chu, K.; Yu, A. Discrete particle simulation of particle–fluid flow: Model formulations and their applicability. J. Fluid Mech. 2010, 661, 482–510. [Google Scholar] [CrossRef]
  35. Johnson, K.; Kendall, K.; Roberts, A. Surface energy and the contact of elastic solids. Proc. R. Soc. Lond. A Math. Phys. Sci. 1971, 324, 301–303. [Google Scholar]
  36. Liu, G.; Li, S.; Yao, Q. A JKR-based dynamic model for the impact of micro-particle with a flat surface. Powder Technol. 2011, 207, 215–223. [Google Scholar] [CrossRef]
  37. Magalhães, H.; Gomez, R.; Leite, B.; Nascimento, J.; Brito, M.; Araújo, M.; Cavalccante, D.; Lima, E.; Lima, A.; Neto, S. Investigating the dialysis treatment using hollow fiber membrane: A new approach by CFD. Membranes 2022, 12, 710. [Google Scholar] [CrossRef] [PubMed]
  38. Anderson, T.; Jackson, R. Fluid-mechanical description of fluidized beds. Equations of motion. Ind. Eng. Chem. Fundamentals 1967, 6, 527–539. [Google Scholar] [CrossRef]
  39. Wang, G.; Hao, W.; Wang, J. Discrete Element Method and Its Practice on EDEM, 1st ed.; Northwestern Polytechnical University Press: Xi’an, China, 2010; pp. 21–24. [Google Scholar]
  40. Mohamed, A.; Gutierrez, M. Comprehensive study of the effects of rolling resistance on the stress-strain and strain localization behavior of granular materials. Granul. Matter 2010, 12, 527–541. [Google Scholar] [CrossRef]
  41. Huang, B.; Yao, Q.; Li, S.; Zhao, H.; Song, Q.; You, C. Experimental investigation on the particle capture by a single fiber using microscopic image technique. Powder Technol. 2006, 168, 1125–1133. [Google Scholar] [CrossRef]
  42. Yang, R.; Zou, R.; Yu, A. Computer simulation of the packing of fine particles. Phys. REV. E 2000, 62, 3900–3908. [Google Scholar] [CrossRef] [PubMed]
  43. Dingwell, K.; Sedin, M.; Theliander, H. Filtration of lignin from a lignocellulose-based ethanol pilot plant. Environ. Eng. Sci. 2011, 28, 775–779. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart of the CFD-DEM coupling.
Figure 1. Flowchart of the CFD-DEM coupling.
Membranes 13 00382 g001
Figure 2. Pore structure of the ceramic membrane and the size of the computational domain. (a) Isometric view; (b) YZ plan view.
Figure 2. Pore structure of the ceramic membrane and the size of the computational domain. (a) Isometric view; (b) YZ plan view.
Membranes 13 00382 g002
Figure 3. Mesh model of the fluid domain.
Figure 3. Mesh model of the fluid domain.
Membranes 13 00382 g003
Figure 4. Change in the deposition morphology of the lignin particles with time. (a) A (T = 0.1 ms); (b) B (T = 0.5 ms); (c) C (T = 1.0 ms); (d) D (T = 1.5 ms); (e) E (T = 2.0 ms); (f) F (T = 2.2 ms).
Figure 4. Change in the deposition morphology of the lignin particles with time. (a) A (T = 0.1 ms); (b) B (T = 0.5 ms); (c) C (T = 1.0 ms); (d) D (T = 1.5 ms); (e) E (T = 2.0 ms); (f) F (T = 2.2 ms).
Membranes 13 00382 g004
Figure 5. Pressure drop curve with time.
Figure 5. Pressure drop curve with time.
Membranes 13 00382 g005
Figure 6. Curves of the cumulative number of penetrated particles and efficiency in capturing lignin particles with time.
Figure 6. Curves of the cumulative number of penetrated particles and efficiency in capturing lignin particles with time.
Membranes 13 00382 g006
Figure 7. Impacts of the rolling friction coefficient on the deposition morphology of lignin particles. (a) μ p m = 0.1 ; (b) μ p m = 1.0 ; (c) μ p m = 2.0 ; (d) μ p m = 3.0 .
Figure 7. Impacts of the rolling friction coefficient on the deposition morphology of lignin particles. (a) μ p m = 0.1 ; (b) μ p m = 1.0 ; (c) μ p m = 2.0 ; (d) μ p m = 3.0 .
Membranes 13 00382 g007aMembranes 13 00382 g007b
Figure 8. Impacts of the rolling friction coefficient on average coordination number.
Figure 8. Impacts of the rolling friction coefficient on average coordination number.
Membranes 13 00382 g008
Figure 9. Impacts of the rolling friction coefficient among the lignin particles on coordination number distribution. (a) μ p m = 0.1 ; (b) μ p m = 1.0 ; (c) μ p m = 2.0 ; (d) μ p m = 3.0 .
Figure 9. Impacts of the rolling friction coefficient among the lignin particles on coordination number distribution. (a) μ p m = 0.1 ; (b) μ p m = 1.0 ; (c) μ p m = 2.0 ; (d) μ p m = 3.0 .
Membranes 13 00382 g009
Figure 10. Impacts of the rolling friction coefficient between the lignin particles and membranes on the coordination number distribution (a) μ p p = 0.6 ; (b) μ p p = 1.8 ; (c) μ p p = 3.0 .
Figure 10. Impacts of the rolling friction coefficient between the lignin particles and membranes on the coordination number distribution (a) μ p p = 0.6 ; (b) μ p p = 1.8 ; (c) μ p p = 3.0 .
Membranes 13 00382 g010
Figure 11. Impacts of the rolling friction coefficient on porosity.
Figure 11. Impacts of the rolling friction coefficient on porosity.
Membranes 13 00382 g011
Table 1. Physical parameters.
Table 1. Physical parameters.
MaterialParticleMembraneBlack Liquor
Diameter (μm)110-
Density (kg/m3)145131001004
Shear modulus (Pa)2 × 1077 × 1010-
Poisson’ ratio0.250.2-
Viscosity (Pa·s)--1.467
Velocity (m/s)0.5-0.5
Table 2. Collision and simulation parameters [17,23,31].
Table 2. Collision and simulation parameters [17,23,31].
Collision ParametersCoefficient of RestitutionCoefficient of Static FrictionCoefficient of Rolling FrictionSurface Energy (J/m2)
Particle–particle0.12.00.1–30.6
Particle–membrane0.12.00.1–31
Simulation ParametersParticle Generation Rate/sTotal Number of ParticlesTime Step of DEM/sTime Step of CFD/s
5 × 10510001 × 10−101 × 10−8
Table 3. Mesh independence verification.
Table 3. Mesh independence verification.
GroupMesh QuantityPressure Drop (Pa)
128,900527.34924
237,544527.9762
346,400528.85345
450,270529.08069
559,048529.62501
Table 4. Proportion of the number of contacts between the lignin particles and membranes to the total.
Table 4. Proportion of the number of contacts between the lignin particles and membranes to the total.
μ p m 0.11.02.03.0
μ p p
0.113.8%17.3%13.7%17.3%
0.612.1%11.1%10.4%11.3%
1.216.7%11.2%14.6%15.5%
1.818.4%12.8%13.0%14.2%
2.417.5%15.1%14.4%16.3%
3.017.2%14.2%16.4%13.1%
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.

Share and Cite

MDPI and ACS Style

Wang, H.; Wang, X.; Wu, Y.; Wang, S.; Wu, J.; Fu, P.; Li, Y. Study of CFD-DEM on the Impact of the Rolling Friction Coefficient on Deposition of Lignin Particles in a Single Ceramic Membrane Pore. Membranes 2023, 13, 382. https://doi.org/10.3390/membranes13040382

AMA Style

Wang H, Wang X, Wu Y, Wang S, Wu J, Fu P, Li Y. Study of CFD-DEM on the Impact of the Rolling Friction Coefficient on Deposition of Lignin Particles in a Single Ceramic Membrane Pore. Membranes. 2023; 13(4):382. https://doi.org/10.3390/membranes13040382

Chicago/Turabian Style

Wang, Hao, Xinyuanrui Wang, Yongping Wu, Song Wang, Junfei Wu, Ping Fu, and Yang Li. 2023. "Study of CFD-DEM on the Impact of the Rolling Friction Coefficient on Deposition of Lignin Particles in a Single Ceramic Membrane Pore" Membranes 13, no. 4: 382. https://doi.org/10.3390/membranes13040382

APA Style

Wang, H., Wang, X., Wu, Y., Wang, S., Wu, J., Fu, P., & Li, Y. (2023). Study of CFD-DEM on the Impact of the Rolling Friction Coefficient on Deposition of Lignin Particles in a Single Ceramic Membrane Pore. Membranes, 13(4), 382. https://doi.org/10.3390/membranes13040382

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop