Original Research PaperThe numerical analysis of particle-size distribution of clusters in shear flow at one-dimensional closed system and three-dimensional open system
Graphical abstract
Introduction
To improve the performance of a resin-based composite, such as a semiconductor-encapsulation material, electric-wire coating, and insulating composite for high-voltage equipment, the filler needs to become finer. As the filler becomes finer, the surface-area-to volume ratio becomes larger, so it easily aggregates due to the surface effect. This aggregation can greatly affect the performance of the resin composite material in terms of heat resistance, insulation, and thermal conductivity. Therefore, to maximize the performance, it is necessary to optimize the kneading process and the molding process at the design stage in a short time. To achieve that optimization, it is effective to numerically analyze the particle dispersion during the kneading process and the molding process.
Methods for numerical analysis of solid-liquid multiphase flow with the particle aggregation and dispersion can be roughly categorized in terms of the model they use, namely, a solid-liquid two-phase flow model or a mixture model. In the solid-liquid two-phase flow model, solid-phase particle motion and liquid-phase flow are calculated, and the interaction between the solid and liquid phases is evaluated. Various numerical methods, which differ in terms of the means used to calculate motion and interaction of each phase, have been proposed [1], [2], [3], [4], [5], [6], [7], [8], [9]. In the mixture model, which assumes a virtual fluid in which solid-phase particles are dispersed in the liquid phase, solid-phase particle motion is calculated from the stress applied by the flow of the virtual fluid [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25], [26]. However, the relative velocity between the solid and liquid phases is not calculated, so the interaction between the solid and liquid phases due to the relative velocity cannot be taken into account. In addition, it is necessary to accurately model the aggregation and dispersion kinetics of the solid-phase particles by using the flow of the virtual fluid. Although under these restrictions, the calculation cost of the mixture model is lower than that of the two-phase flow model, it is effective in estimating the particle dispersion in solid-liquid multiphase flow during molding process.
The two-dimensional numerical-simulation method proposed Masuda et al. [10], [11], [12] calculates the aggregation and dispersion behaviors of clusters in the microchannel flow on the basis of the thixotropy model established by Usui [13], [14]. In that thixotropy model, particle aggregation and dispersion kinetics are modeled by Brownian coagulation, shear coagulation, and shear breakage. In the model of Masuda et al., a non-uniform shear flow is calculated by computational fluid dynamics (CFD), and the particle aggregation and dispersion rate in that shear flow is evaluated by the thixotropy model. Then, the particle-size distribution (PSD) of the cluster at each element of the CFD model is calculated by solving the population-balance equation (PBE) for the clusters. The apparent-solid-volume fraction including the voids in the cluster is estimated from the obtained PSD, and the change in suspension viscosity is calculated. When applying the model of Masuda et al. to optimization of the molding process of resin composites, it is necessary to extend it to large-size clusters and three-dimensional flows with free surfaces.
In this study, which aimed to develop a numerical-analysis method for particle dispersion Masuda et al.'s method was extended to three-dimensional space. In three-dimensionalization, the PBE was discretized using the fixed-pivot technique of Kumar and Ramkrishna [27], [28] to solve the PBE at low calculation cost. The proposed method was incorporated into general-purpose CFD software FLOW-3D®, and its accuracy in a numerical simulation of backward-step flow was proved.
Section snippets
Numerical methods
The following PBE for calculating the PSD profile in a shear flow was developed by Masuda et al. [10], [11], [12]:where t [s] is time, i [-], j [-], and k [-] are numbers of primary particles in the cluster, ni [m−3] is
One-dimensional closed system (Case 1)
To evaluate whether the PBE results depend on the initial condion of the PSD, time-dependent changes of the average number of particles in cluster at various initial condition in Case 1-1 are shown in Fig. 4. The initial condition of PSD were set to ϕc,1 = 1, ϕc,5 = 1, ϕc,10 = 1, ϕc,15 = 1, and ϕc,20 = 1, respectively. In this case, the PBE results at steady state do not change when the PBE calcuration starts from well-dispersed smaller particles (ϕc,1 = 1) or larger aggregates (ϕc,20 = 1).
Conclusions
A method for analyzing the three-dimensional PSD in shear flow was proposed. The PSD is calculated by taking into account the aggregation and dispersion rate of the cluster by Brownian coagulation, shear coagulation, and shear breakage. The proposed method was incorporated into the general-purpose CFD software FLOW-3D®, and its accuracy was proved in a numerical simulation of backward-step flow. In the downstream region, the numerical and experiment result of distributions of average number of
References (31)
- et al.
Direct simulation of flowing colloidal dispersions by smoothed profile method
Adv. Powder Technol.
(2010) - et al.
Motion of particle agglomerate involving interparticle force in dilute suspension
Powder Technol.
(2008) Incorporation of smooth spherical bodies in the Lattice Boltzmann method
J. Comput. Phys.
(2011)- et al.
Large-scale CFD-DEM simulations of fluidized granular systems
Chem. Eng. Sci.
(2013) - et al.
Diffusion-based coarse graining in hybrid continuum–discrete solvers: applications in CFD–DEM
Int. J. Multiph. Flow.
(2015) - et al.
Model analysis on dispersion characteristics of fine particles in Newtonian molten polymer
Adv. Powder Technol.
(2009) - et al.
Comparison of numerical methods for solving population balance equations incorporating aggregation and breakage
Powder Technol.
(2009) - et al.
Discrete bivariate population balance modelling of heteroaggregation processes
J. Colloid Interface Sci.
(2009) - et al.
Population balance discretization for growth, attrition, aggregation, breakage and nucleation
Comput. Chem. Eng.
(2016) - et al.
Particle size distributions and viscosity of suspensions undergoing shear-induced coagulation and fragmentation
Chem. Eng. Sci.
(2003)
A modified aggregation based model for the accurate prediction of particle distribution and viscosity in magnetic nanofluids
Powder Technol.
Prediction of hydrodynamic and optical properties of TiO2/water suspension considering particle size distribution
Int. J. Heat Mass Transf.
Improved accuracy and convergence of discretized population balance for aggregation: the cell average technique
Chem. Eng. Sci.
On the solution of population balance equations by discretization—I. A fixed pivot technique
Chem. Eng. Sci.
Simulation method of colloidal suspensions with hydrodynamic interactions: fluid particle dynamics
Phys. Rev. Lett.
Cited by (2)
An analytical solution of the population balance equation for simultaneous Brownian and shear coagulation in the continuum regime
2020, Advanced Powder TechnologyCitation Excerpt :Guichard et al. [18] modelled the dynamics of nanoparticles under Brownian motion and turbulent effects using a differential-algebraic framework. Shimada et al. [19] analysed the three dimensional particle size distribution in a shear flow considering Brownian coagulation, shear coagulation and shear breakage. Besides these numerical studies, Swift and Friedlander [20] provided a theoretical analysis of the change in the particle number concentration for simultaneous Brownian coagulation and laminar shear in the continuum regime.