# A Three-Dimensional Simulation of Particle Distribution in a Separator and Structure Optimization with the Statistical Approach of Taguchi Method

Authors: Jinjin Liu; Kai Liu; Tong Zhao; Zhuofei Xu
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

A three-dimensional numerical simulation combining discrete phase method (DPM) and porous media based on the theory of Euler-Lagrange has been employed to investigate particles distribution in a separator. The DPM model is applied to monitor the movement of individual particles and calculate the contact force between them in the separator. The simulation results display the migration feature of dust particles over time and the distribution of particles on the surface element in porous region and reveal that the flow field influences the distribution uniformity of the particles in porous area directly. Based on the analysis, the structure of separator is optimized by the Taguchi method. An orthogonal relation motion has been established. The optimal solution is achieved by the calculation of the weight relationship. The calculated optimal structure is evaluated by the signal to noise (SNR). The result reveals that the values of SNR in case are eligible. As a result, the research of the separator points out a useful and improvable method for the parameter optimization of structure design.


---

## Body

## 1. Introduction

Coal, as the primal fossil fuel, has played an important role in national economy, such as the areas of electricity, heating, and transportation. Coal resource can be divided into high rank coal and low rank coal. Among them, low rank coal is not only characterized by the high moisture content but also the high volatility which can cause an explosion easily [1]. Meanwhile, high temperature gas containing dust is generated during burning, which causes serious pollution. Therefore, particles elimination is necessary before blowdown. Separator is important equipment in the filtration and removal process of low rank coal. Its core component is the filter core, which has the characteristics of a high temperature resistance and pressure resistance, strong corrosion resistance, and high filtration efficiency.In this paper, multiphase flow model is chosen for simulation analysis of the model combined with computational fluid dynamics (CFD). The commonly used mathematical models that have been raised so far can be grouped into two categories: one is the continuum-continuum approach, as represented by Euler-Euler model [2, 3], and the other is the continuum discrete approach [4] represented by the Euler-Lagrange model. The flow field of gas is treated as a continuum phase by Navies-Stokes equations in computational cell scale while the flow field of particles is treated as a discrete phase based on DPM. The physical status of individual particles, such as the trajectory in space, is obtained by the calculation of continuum phase of the air. A two-phase model of gas-solid flow is used to simulate the coupling of porous media module.Many scholars have studied the porous media model in the different structure. Wu [5] and Boysen [6] simulated the process of adsorption in porous. Yu Liming et al. [7] used CFD-DEM to simulate the process of pressure variation during filtration process. Wang et al. [8] researched the separation efficiency of particles with different sizes in the gas-solid flow. Zhao Tong et al. [9] researched the gas-solid flow whose particle movement and distribution are in a down-flow fluidized bed with CFD-DEM. M. Haim et al. [10] researched the particle-gas flows whose effect of the inlet conditions is on the numerical solutions. Xi Gao et al. [11] made a research on the evaluation of an enhanced filtered drag model compared with other eight different drag models in the three-dimensional simulations of gas-solid fluidized beds.However, minimum studies considered the more uniformity of particles in porous media by optimizing the structural parameters. Optimization design methods have used various metaheuristics such as nonlinear search optimization [12], the Taguchi method.The Taguchi method has been proposed by Genichi Taguchi [13, 14] and widely used in many fields. Rout [15] applied the Taguchi method to find the optimal setting of parameters in a manipulator and increase the positional accuracy. T. Mothilal et al. [16] applied the Taguchi method to optimize the operational parameters including the solid particle diameter, inlet air velocity, and inlet air temperature in a cyclone separator.In summary, the paper is organized as follows: Section2 introduces the basic theory of CFD, Euler-Lagrange Model, and the porous medium model as well; then the result of numerical simulation is discussed in a three-dimensional model of separator such as the migration distribution diagram of particles over time and the velocity distribution in filter area. Section 3 evaluates the model by using the statistics method of Taguchi to determine the weight of different structure parameters and get an effective improved plan.


## 2. The Approach of Modeling

### 2.1. Computational Fluid Theory

Computational fluid dynamics method is used in this study. Its basic approach is dividing the continuous physics region into discrete. Based on the Navier-Stokes equations, the fluid domain is calculated to obtain the key physical parameters of heated air such as speed, temperature, pressure, and flow positions.The flow condition can be judged from the Reynolds number. The turbulence model is selected which includes the conversation laws of mass, momentum, and energy. The turbulence equation should be considered. Thus k-ε model is taken as follows:(1)∂ρk∂t+∂ρkui∂xi=∂∂xjμ+μσk∂k∂xj+Gk+Gb-ρε-YM+Sk(2)∂ρε∂t+∂ρεui∂xi=∂∂xjμ+μσk∂k∂xj+G1εεk·Gk+C3εGb-G2ερε2k+SεHere k is turbulent kinetic energy in (1) and (2). ε means the dissipation rating. And we can know that Gk and Gb are the turbulent kinetic energy generation introduced by mean velocity gradient and buoyancy.YM is the contribution in the turbulent pulsating pressure expansion. C1ε,C2ε, and C3ε are empirical constants. σk and σε are Prandtl numbers to k and ε correspondingly. Sk and Sε are user defined options.


### 2.2. The Porous Medium Model

The filter area in separator needs to be calculated with model of porous media which is defined as laminar. A flow resistance value is introduced. The porous medium model based on superficial velocity can simulate pressure loss. The flow resistance of porous material in the fluid domain is calculated by adding a feature in momentum equation. The feature is composed of two parts, Darcy viscous resistance terms and inertial loss terms, which is expressed as(3)Sj=-∑j=13Dijμvj+∑j=13Cij12ρvvjHereDij andCij is the viscous resistance and inertia loss coefficient matrix, respectively. This negative result leads to a pressure drop in the unit of porous medium. μ is the fluid viscosity coefficient. The viscous resistance and inertial resistance is calculated by the following formula (4) and (5):(4)1α=150Dp21-ε2ε3(5)C2=3.5Dp1-ε2ε3whereDp is the particle diameter, ε is the porosity, and α is the permeability.


### 2.3. The Coupling of Discrete Phase and Porous Medium

The particles trajectory is tracked using the DPM model. The particle motion formula is on a Lagrange coordinate system. The particles motion formula of discrete phase is shown in (6) and the inertia force is considered. (6)dupdt=FDu-up+gxρp-ρρp+FxHere up is the velocity of particle.ρ and ρp are the densities of the continuous and discrete phase, respectively. gx is the acceleration of gravity in x direction; g(ρp-ρ)/ρp represents the joint force of gravity and buoyancy when particle is moving;FD(u-up) is the air resistance acting on the particle in formula (6); andFD represents the resistance impacting on particles in unit-mass:(7)FD=18μρpdp2CDRe24(8)Re=ρdpup-uμ(9)CD=a1+a2Re+a3ReR e is the relative Reynolds number andCD is the product ratio of fluid dynamic pressure.FD has impact on the projection which is the movement direction of the particles, namely, the drag force coefficient.dp is the particle diameter here. The control equations of coupling model in the continuous phase, discrete phase, and porous media [17, 18] are as follows:(10)ερCpf∂Tf∂t+ρCpεv∇Tf=ε∇kf∇Tf+1-ε6hfdTs-Tf(11)ρCs∂Ts∂t=∇ks∇Ts+6hfdTs-Tf+ρs∑i-∇Hi∂qi∂tC p and Cs are the specific heat at constant pressure of the continuous and discrete phase, respectively; hƒ is the thermal coefficient of gas and solid phase; kƒ is axial heat transfer coefficient for the gas phase; ks is thermal conductivity coefficient for solid phase, unit for; Tƒ and Ts are the temperature of gas and solid phase, respectively. Here transfer rate of mass equation is used:(12)∂qi∂t=kiqi∗-qiAdsorption balance equations [19, 20] are as (13)qi∗=Kpi1+∑k=1NbkPk,K=kiexp⁡k2Tf,k3exp⁡k4Tfk i is the mass transfer coefficient of the component i; qi is the adsorption quantity for solid phase of componenti; q∗ is the maximum adsorption amount for solid phase components.


### 2.4. Principle and Simulation Conditions

#### 2.4.1. Separator Principle

The computational domain and grid arrangements of the separator are shown in Figure1. The simulation domain consists of inlet, outlet, six Venturi tubes, and inferior vena including six filter parts. The diameters of inlet and outside are both 400 mm. The angle between the centerlines of the inlet and the horizontal direction is 45°. There are six Venturi tubes in the separator. What is more, a back pulse made by vibration wave is used to reduce particles accumulated in filters elements. There are six filters in inferior vena which can be seen as porous mediums. Each filter has a diameter of 840 mm and a length of 1980 mm.Figure 1
Calculation domain and grid arrangement.The semi-implicit method for pressure-link equation (SIMPLE) algorithm is used as the numerical method. In computational fluid dynamics, the SIMPLE algorithm is widely used to solve the Navier-Stokes equations, and it has been extensively used to solve fluid flow and heat transfer problems [21, 22]. A no-slip condition is set to the air phase on the walls. The simulation is started with the random generation of particles combined with gas phase in inlet. With the setting of gravitational force, lift force, and drag force settling, the particles enter into the inferior vena and pervade in the area of porous mediums. Finally, the part of the particles will sink into bottom; others will pass filter holes and reach the cavity Venturi tube and the export. In general we want to keep less particles away from passing through the porous area. Thus the filter will work for a longer time and not be broken easily. We hope more particles are blown down to the bottom of separator directly.


#### 2.4.2. Simulation Conditions

Calculation domain is meshed by Gambit. The total number of meshes is 6,197,815; therefore, the time of calculation of the meshed model needs about 48 hours on a workstation (System is Windows 8. The processor is Xeon 1231v3. The memory is 16G). The calculation is done in Fluent 14.0 and the result is extracted by Tecplot 360 21013R1. The type of grid is chosen as tetrahedron structure. The grid density is increased appropriately in the six filter parts. The Jacobian value is positive.In this study, we set medium volatile coal (Coal-Mv) as the particles of discrete phase in calculation. Compared with highly volatile coal (Coal-Hv) and low volatile coal (Coal-Lv), Coal-Mv is suitable for simulating most coal powders in working conditions [23]. The physical parameters are set according to the real property in Table 1.Table 1
The parameter settings of calculation.
Gas phase Density(kg/m3) 1.225 Viscosity(kg/m-s) 1.7894e-05 Reference Temperature(k) 298.15 Velocity Magnitude(m/s) 20 Turbulent Intensity(%) 10 Hydraulic Diameter(m) 0.4  Particles phase Cp(Specific Heat )(j/kg-k) 1550 Molecular Weight (kg/kg-mol) 17.237 Standard State Enthalpy(j/kg-mol) -5.601e+07 Reference Temperature(k) 298.15 Velocity (m/s) 20 Turbulent Kinetic Energy 0.8 Momentum 0.2  Porous area Porosity 0.4 face permeability(m2) 1.19e-11 Porous permeability 0.5 Viscous Resistance 8.4e+10 Inertial Resistance 5.1e+04The Rosin-Rammler model is used to describe the distribution of particle diameters in this work. The Rosin-Rammler model was first applied by Rosin-Rammler to describe a particle size distribution. The expression function of the Rosin-Rammler is described as(14)Fd=1-exp⁡-dcmwhere F is the distribution function of particles with a diameter of d. c is the characteristic diameter which is always the average diameter of particles. The exponent m is uniformity constant [24, 25].Inlet and outlet of the mode are mentioned in Figure1. The area of filters is defined as the laminar zone porous. Other parameters in Ansys are chosen in Table 1. In Table 1, air is simulated as gas phase while particles phase is simulated as discrete phase. Filters are represented as porous area in separator. The physical parameters are set according to the real property. As to the velocity, diameter, energy, and momentum, average values are chosen to initialize and calculate. Momentum is the subrelaxation factor in Ansys 14.0 to reduce appropriately between the calculation results of the previous layer and the later one, which could avoid the nonlinear divergence induced by excessive difference during the iterative process.In order to make result closer to reality, other options are set based on the real actual working conditions.


### 2.5. Simulation Result and Discussion

This section discusses the results of the simulation. The result positions of the three cross-sectional in porous region in Figure2. Particle distribution in the three sections is shown in Table 2. We can find the following: (1) The particles enter into the inferior vena of separator, combined with continuous phase. Most particles cannot penetrate the porous media because there are viscous resistance and inertial resistance. (2) With the increase of time, particles slowly goes from import to inferior vena area. The number of particles gradually increases from bottom to top. Meanwhile, a large number of particles accumulated in the middle area. (3) Due to the fact that the particle phase is loaded continuously and the pressure is increased in the inferior vena, particles slowly enter into the filter which is close to inlet about the time of 23.4 s. (4) The color of the particles indicates that the particles sink slowly to bottom of separator under the effect of gravity with time increases.Table 2
Particle distribution at different times.
position Simulation result T:126s T:162s T:198s T:234s Y=-1100mm  Y=-1500mm  Y=-1900mmFigure 2
Three-dimensional particle distribution (t=286s).Table3 shows the velocity distribution images in three cross-sectional. These sections illustrate the flow field of continuous phase under the influence of solid phase, and solid phase affects velocity distribution of the continuous phase. Take 60 speed values in each face to calculate the mean square error for each surface and it is shown in Table 4. Yi represents the mean square error of sample values in each cross section which is calculated as (15)Yi=1N∑i=1Nxi-μ2Table 3
The image of velocity distribution.
position T:126s T:162s T:198s T:234s Y=-1100mm  Y=-1500mm  Y=-1900mmTable 4
The results of mean square error.
Position of section 126s 162s 198s 234s -1100mm 10.877 9.957 9.073 8.838 -1500mm 7.886 7.852 8.083 5.441 -1900mm 6.777 6.777 5.685 6.230N is the number of data point and xi is the value of each point.μ is the average value of all data points. The mean square error of each section is relatively large.In engineering, the replacement of filter has a high cost and complicated operation. Therefore, we hope that most particles will not penetrate filter. A large number of particles accumulated in the middle area will cause filter to be partially stressed resulting in the aging rate of the filter, so particle should distribute uniformly in the filter to reduce the accumulation of particles partially in the filter. The particles are mainly affected by the continuous phase. Therefore, the goal of the separator structure optimization is to get a more homogeneous flow field in porous media zone. In our paper, the uniformity of flow field is evaluated by velocity and pressure.


## 2.1. Computational Fluid Theory

Computational fluid dynamics method is used in this study. Its basic approach is dividing the continuous physics region into discrete. Based on the Navier-Stokes equations, the fluid domain is calculated to obtain the key physical parameters of heated air such as speed, temperature, pressure, and flow positions.The flow condition can be judged from the Reynolds number. The turbulence model is selected which includes the conversation laws of mass, momentum, and energy. The turbulence equation should be considered. Thus k-ε model is taken as follows:(1)∂ρk∂t+∂ρkui∂xi=∂∂xjμ+μσk∂k∂xj+Gk+Gb-ρε-YM+Sk(2)∂ρε∂t+∂ρεui∂xi=∂∂xjμ+μσk∂k∂xj+G1εεk·Gk+C3εGb-G2ερε2k+SεHere k is turbulent kinetic energy in (1) and (2). ε means the dissipation rating. And we can know that Gk and Gb are the turbulent kinetic energy generation introduced by mean velocity gradient and buoyancy.YM is the contribution in the turbulent pulsating pressure expansion. C1ε,C2ε, and C3ε are empirical constants. σk and σε are Prandtl numbers to k and ε correspondingly. Sk and Sε are user defined options.


## 2.2. The Porous Medium Model

The filter area in separator needs to be calculated with model of porous media which is defined as laminar. A flow resistance value is introduced. The porous medium model based on superficial velocity can simulate pressure loss. The flow resistance of porous material in the fluid domain is calculated by adding a feature in momentum equation. The feature is composed of two parts, Darcy viscous resistance terms and inertial loss terms, which is expressed as(3)Sj=-∑j=13Dijμvj+∑j=13Cij12ρvvjHereDij andCij is the viscous resistance and inertia loss coefficient matrix, respectively. This negative result leads to a pressure drop in the unit of porous medium. μ is the fluid viscosity coefficient. The viscous resistance and inertial resistance is calculated by the following formula (4) and (5):(4)1α=150Dp21-ε2ε3(5)C2=3.5Dp1-ε2ε3whereDp is the particle diameter, ε is the porosity, and α is the permeability.


## 2.3. The Coupling of Discrete Phase and Porous Medium

The particles trajectory is tracked using the DPM model. The particle motion formula is on a Lagrange coordinate system. The particles motion formula of discrete phase is shown in (6) and the inertia force is considered. (6)dupdt=FDu-up+gxρp-ρρp+FxHere up is the velocity of particle.ρ and ρp are the densities of the continuous and discrete phase, respectively. gx is the acceleration of gravity in x direction; g(ρp-ρ)/ρp represents the joint force of gravity and buoyancy when particle is moving;FD(u-up) is the air resistance acting on the particle in formula (6); andFD represents the resistance impacting on particles in unit-mass:(7)FD=18μρpdp2CDRe24(8)Re=ρdpup-uμ(9)CD=a1+a2Re+a3ReR e is the relative Reynolds number andCD is the product ratio of fluid dynamic pressure.FD has impact on the projection which is the movement direction of the particles, namely, the drag force coefficient.dp is the particle diameter here. The control equations of coupling model in the continuous phase, discrete phase, and porous media [17, 18] are as follows:(10)ερCpf∂Tf∂t+ρCpεv∇Tf=ε∇kf∇Tf+1-ε6hfdTs-Tf(11)ρCs∂Ts∂t=∇ks∇Ts+6hfdTs-Tf+ρs∑i-∇Hi∂qi∂tC p and Cs are the specific heat at constant pressure of the continuous and discrete phase, respectively; hƒ is the thermal coefficient of gas and solid phase; kƒ is axial heat transfer coefficient for the gas phase; ks is thermal conductivity coefficient for solid phase, unit for; Tƒ and Ts are the temperature of gas and solid phase, respectively. Here transfer rate of mass equation is used:(12)∂qi∂t=kiqi∗-qiAdsorption balance equations [19, 20] are as (13)qi∗=Kpi1+∑k=1NbkPk,K=kiexp⁡k2Tf,k3exp⁡k4Tfk i is the mass transfer coefficient of the component i; qi is the adsorption quantity for solid phase of componenti; q∗ is the maximum adsorption amount for solid phase components.


## 2.4. Principle and Simulation Conditions

### 2.4.1. Separator Principle

The computational domain and grid arrangements of the separator are shown in Figure1. The simulation domain consists of inlet, outlet, six Venturi tubes, and inferior vena including six filter parts. The diameters of inlet and outside are both 400 mm. The angle between the centerlines of the inlet and the horizontal direction is 45°. There are six Venturi tubes in the separator. What is more, a back pulse made by vibration wave is used to reduce particles accumulated in filters elements. There are six filters in inferior vena which can be seen as porous mediums. Each filter has a diameter of 840 mm and a length of 1980 mm.Figure 1
Calculation domain and grid arrangement.The semi-implicit method for pressure-link equation (SIMPLE) algorithm is used as the numerical method. In computational fluid dynamics, the SIMPLE algorithm is widely used to solve the Navier-Stokes equations, and it has been extensively used to solve fluid flow and heat transfer problems [21, 22]. A no-slip condition is set to the air phase on the walls. The simulation is started with the random generation of particles combined with gas phase in inlet. With the setting of gravitational force, lift force, and drag force settling, the particles enter into the inferior vena and pervade in the area of porous mediums. Finally, the part of the particles will sink into bottom; others will pass filter holes and reach the cavity Venturi tube and the export. In general we want to keep less particles away from passing through the porous area. Thus the filter will work for a longer time and not be broken easily. We hope more particles are blown down to the bottom of separator directly.


### 2.4.2. Simulation Conditions

Calculation domain is meshed by Gambit. The total number of meshes is 6,197,815; therefore, the time of calculation of the meshed model needs about 48 hours on a workstation (System is Windows 8. The processor is Xeon 1231v3. The memory is 16G). The calculation is done in Fluent 14.0 and the result is extracted by Tecplot 360 21013R1. The type of grid is chosen as tetrahedron structure. The grid density is increased appropriately in the six filter parts. The Jacobian value is positive.In this study, we set medium volatile coal (Coal-Mv) as the particles of discrete phase in calculation. Compared with highly volatile coal (Coal-Hv) and low volatile coal (Coal-Lv), Coal-Mv is suitable for simulating most coal powders in working conditions [23]. The physical parameters are set according to the real property in Table 1.Table 1
The parameter settings of calculation.
Gas phase Density(kg/m3) 1.225 Viscosity(kg/m-s) 1.7894e-05 Reference Temperature(k) 298.15 Velocity Magnitude(m/s) 20 Turbulent Intensity(%) 10 Hydraulic Diameter(m) 0.4  Particles phase Cp(Specific Heat )(j/kg-k) 1550 Molecular Weight (kg/kg-mol) 17.237 Standard State Enthalpy(j/kg-mol) -5.601e+07 Reference Temperature(k) 298.15 Velocity (m/s) 20 Turbulent Kinetic Energy 0.8 Momentum 0.2  Porous area Porosity 0.4 face permeability(m2) 1.19e-11 Porous permeability 0.5 Viscous Resistance 8.4e+10 Inertial Resistance 5.1e+04The Rosin-Rammler model is used to describe the distribution of particle diameters in this work. The Rosin-Rammler model was first applied by Rosin-Rammler to describe a particle size distribution. The expression function of the Rosin-Rammler is described as(14)Fd=1-exp⁡-dcmwhere F is the distribution function of particles with a diameter of d. c is the characteristic diameter which is always the average diameter of particles. The exponent m is uniformity constant [24, 25].Inlet and outlet of the mode are mentioned in Figure1. The area of filters is defined as the laminar zone porous. Other parameters in Ansys are chosen in Table 1. In Table 1, air is simulated as gas phase while particles phase is simulated as discrete phase. Filters are represented as porous area in separator. The physical parameters are set according to the real property. As to the velocity, diameter, energy, and momentum, average values are chosen to initialize and calculate. Momentum is the subrelaxation factor in Ansys 14.0 to reduce appropriately between the calculation results of the previous layer and the later one, which could avoid the nonlinear divergence induced by excessive difference during the iterative process.In order to make result closer to reality, other options are set based on the real actual working conditions.


## 2.4.1. Separator Principle

The computational domain and grid arrangements of the separator are shown in Figure1. The simulation domain consists of inlet, outlet, six Venturi tubes, and inferior vena including six filter parts. The diameters of inlet and outside are both 400 mm. The angle between the centerlines of the inlet and the horizontal direction is 45°. There are six Venturi tubes in the separator. What is more, a back pulse made by vibration wave is used to reduce particles accumulated in filters elements. There are six filters in inferior vena which can be seen as porous mediums. Each filter has a diameter of 840 mm and a length of 1980 mm.Figure 1
Calculation domain and grid arrangement.The semi-implicit method for pressure-link equation (SIMPLE) algorithm is used as the numerical method. In computational fluid dynamics, the SIMPLE algorithm is widely used to solve the Navier-Stokes equations, and it has been extensively used to solve fluid flow and heat transfer problems [21, 22]. A no-slip condition is set to the air phase on the walls. The simulation is started with the random generation of particles combined with gas phase in inlet. With the setting of gravitational force, lift force, and drag force settling, the particles enter into the inferior vena and pervade in the area of porous mediums. Finally, the part of the particles will sink into bottom; others will pass filter holes and reach the cavity Venturi tube and the export. In general we want to keep less particles away from passing through the porous area. Thus the filter will work for a longer time and not be broken easily. We hope more particles are blown down to the bottom of separator directly.


## 2.4.2. Simulation Conditions

Calculation domain is meshed by Gambit. The total number of meshes is 6,197,815; therefore, the time of calculation of the meshed model needs about 48 hours on a workstation (System is Windows 8. The processor is Xeon 1231v3. The memory is 16G). The calculation is done in Fluent 14.0 and the result is extracted by Tecplot 360 21013R1. The type of grid is chosen as tetrahedron structure. The grid density is increased appropriately in the six filter parts. The Jacobian value is positive.In this study, we set medium volatile coal (Coal-Mv) as the particles of discrete phase in calculation. Compared with highly volatile coal (Coal-Hv) and low volatile coal (Coal-Lv), Coal-Mv is suitable for simulating most coal powders in working conditions [23]. The physical parameters are set according to the real property in Table 1.Table 1
The parameter settings of calculation.
Gas phase Density(kg/m3) 1.225 Viscosity(kg/m-s) 1.7894e-05 Reference Temperature(k) 298.15 Velocity Magnitude(m/s) 20 Turbulent Intensity(%) 10 Hydraulic Diameter(m) 0.4  Particles phase Cp(Specific Heat )(j/kg-k) 1550 Molecular Weight (kg/kg-mol) 17.237 Standard State Enthalpy(j/kg-mol) -5.601e+07 Reference Temperature(k) 298.15 Velocity (m/s) 20 Turbulent Kinetic Energy 0.8 Momentum 0.2  Porous area Porosity 0.4 face permeability(m2) 1.19e-11 Porous permeability 0.5 Viscous Resistance 8.4e+10 Inertial Resistance 5.1e+04The Rosin-Rammler model is used to describe the distribution of particle diameters in this work. The Rosin-Rammler model was first applied by Rosin-Rammler to describe a particle size distribution. The expression function of the Rosin-Rammler is described as(14)Fd=1-exp⁡-dcmwhere F is the distribution function of particles with a diameter of d. c is the characteristic diameter which is always the average diameter of particles. The exponent m is uniformity constant [24, 25].Inlet and outlet of the mode are mentioned in Figure1. The area of filters is defined as the laminar zone porous. Other parameters in Ansys are chosen in Table 1. In Table 1, air is simulated as gas phase while particles phase is simulated as discrete phase. Filters are represented as porous area in separator. The physical parameters are set according to the real property. As to the velocity, diameter, energy, and momentum, average values are chosen to initialize and calculate. Momentum is the subrelaxation factor in Ansys 14.0 to reduce appropriately between the calculation results of the previous layer and the later one, which could avoid the nonlinear divergence induced by excessive difference during the iterative process.In order to make result closer to reality, other options are set based on the real actual working conditions.


## 2.5. Simulation Result and Discussion

This section discusses the results of the simulation. The result positions of the three cross-sectional in porous region in Figure2. Particle distribution in the three sections is shown in Table 2. We can find the following: (1) The particles enter into the inferior vena of separator, combined with continuous phase. Most particles cannot penetrate the porous media because there are viscous resistance and inertial resistance. (2) With the increase of time, particles slowly goes from import to inferior vena area. The number of particles gradually increases from bottom to top. Meanwhile, a large number of particles accumulated in the middle area. (3) Due to the fact that the particle phase is loaded continuously and the pressure is increased in the inferior vena, particles slowly enter into the filter which is close to inlet about the time of 23.4 s. (4) The color of the particles indicates that the particles sink slowly to bottom of separator under the effect of gravity with time increases.Table 2
Particle distribution at different times.
position Simulation result T:126s T:162s T:198s T:234s Y=-1100mm  Y=-1500mm  Y=-1900mmFigure 2
Three-dimensional particle distribution (t=286s).Table3 shows the velocity distribution images in three cross-sectional. These sections illustrate the flow field of continuous phase under the influence of solid phase, and solid phase affects velocity distribution of the continuous phase. Take 60 speed values in each face to calculate the mean square error for each surface and it is shown in Table 4. Yi represents the mean square error of sample values in each cross section which is calculated as (15)Yi=1N∑i=1Nxi-μ2Table 3
The image of velocity distribution.
position T:126s T:162s T:198s T:234s Y=-1100mm  Y=-1500mm  Y=-1900mmTable 4
The results of mean square error.
Position of section 126s 162s 198s 234s -1100mm 10.877 9.957 9.073 8.838 -1500mm 7.886 7.852 8.083 5.441 -1900mm 6.777 6.777 5.685 6.230N is the number of data point and xi is the value of each point.μ is the average value of all data points. The mean square error of each section is relatively large.In engineering, the replacement of filter has a high cost and complicated operation. Therefore, we hope that most particles will not penetrate filter. A large number of particles accumulated in the middle area will cause filter to be partially stressed resulting in the aging rate of the filter, so particle should distribute uniformly in the filter to reduce the accumulation of particles partially in the filter. The particles are mainly affected by the continuous phase. Therefore, the goal of the separator structure optimization is to get a more homogeneous flow field in porous media zone. In our paper, the uniformity of flow field is evaluated by velocity and pressure.


## 3. The Optimizations of Structure for Separator

A separator model is given by an enterprise in China. Our purpose is to improve parameters about inlet and filters based on CFD. As to a separator, performance is always decided by the distribution of velocity and pressure. In ideal status, a uniform distribution of flow field means that the particle will contact the filters equivalently; therefore, all the filters paly the same role. Otherwise, particles will focus on some parts of few filters which will be damaged easily. These filters may be blocked by too much concentrated particles, while an unfair distribution also means that some area of filters does not work at all. For these reasons, mean square errors of velocity and pressure are used to evaluate the result of optimization, and a smaller value means a better performance.As to a given separator in specific working conditions, the mainframe cannot be changed, and the filter is already made according to standard; therefore, their sizes always cannot be changed. Inlet is a key part to decide the distribution of fluid filed; thus both the position and inclination angle of inlet are chosen to optimize. In the inner of separator, the utilization rate of filters is decided by the installation site directly, since it is so, it cannot be neglected that the distance between two filters in the direction-Z in Figure1. The distance-X is neglected here to simplify the process since its result is similar to distance in direction-Z.Based on the reasons above, three parameters are selected to optimize the model according to the mean square error from distribution of both velocity and pressure. Taguchi method is used to realize the target.


### 3.1. The Theory of Taguchi Method

The basic principle of the Taguchi method, avoiding the blindness of traverse type analysis, is the inspection of the effectiveness from the experimental variables which is the synthesis procedure conducted by orthogonal array with minimum number of experiments [15, 26]. This approach can obtain effectively the contribution rate and weight relations from various elements of target. The method can be used to realize the optimization and reduce test times compared to exhaustion method. The main steps of structure optimization for separator are as follows: (1) Determining the key elements of the research object, select all elements of the plan. (2) Design the orthogonal table for all combination schemes according to the relationship between elements. (3) Select the orthogonal schemes and calculate each of them. (4) According to the results, get the order of weight for all selected elements based on the contribution rate.


### 3.2. The Design of Orthogonal Table

In the study based on Taguchi method, to represent three parameters, we take the angle between the inlet axis and the horizontal direction as factor A, the distance between the inlet and bottom of filter as factor B, and the distance of two filters in direction-Z as factor C. Control factors and its levels are designed and shown in Table5.Table 5
The statistics of element.
Factor Level 1 Level 2 Level 3 A (degree) 45 22.5 0 B (m) 1.267 0.725 0.200 C (m) 1.220 1.020 0.820Factor A is often between 0° to 45° in actual equipment, thus 0°, 45°, and 22.5° are chosen. As to factor B and factor C, the extreme values and their medium value are taken into consideration, ranging from the highest and lowest position. Different values are noted as different levels in Taguchi method as shown in Table5.Table6 is the orthogonal table for Taguchi analysis with factors A, B, and C. A correct orthogonal table needs to fulfill with following two conditions: (1) In a column, the numbers of different values should be the same. (2) The same levels have the same frequency in every two columns.Table 6
The proof of the orthogonal.
case A B C AB BC CA 1 + + + + + + 2 + - - - + - 3 + 0 0 0 0 0 4 - + - - - + 5 - - 0 + 0 0 6 - 0 + 0 0 - 7 0 + 0 0 0 0 8 0 - + 0 - 0 9 0 0 - 0 0 0 Total 0 0 0 0 0 0L9 orthogonal array is used here; it contains nine tests, three levels, and three variables (L means Latin square and 9 means the number of cases). The orthogonal table is recorded in Table 7 and it can be used as experiment scheme.Table 7
The orthogonal table of separator.
case A B C 1 A1 B1 C1 2 A1 B2 C2 3 A1 B3 C3 4 A2 B1 C2 5 A2 B2 C3 6 A2 B3 C1 7 A3 B1 C3 8 A3 B2 C2 9 A3 B3 C1 For example, A1 stands for factor A with Level 1.


### 3.3. The Analysis of the SNR

In Taguchi method, the ratio of signal to noise (SNR) is always used as the key value to evaluate the performance of each experiment in Taguchi analysis. In this study, the target is to find smaller mean square errors; thus SNR can be expressed as (16). A bigger value of SNR is needed since it is negatively correlated with Y [27].(16)η=-10lg1n∑i=1nYi2η represents the SNR.n is the number of replications. Since it is expected to be smaller that both the mean square errors of velocity and pressure, a smaller value is what we needed. In Taguchi method, the ratio of signal to noise (SNR) is used instead of mean square errors as SNR can expand the range of smaller values. In (16), Yi2 is used to eliminate the negative values and an average is calculated. A logarithmic relation can amplify the average of Yi2 when it is under 1; this means the contrast ratio for smaller values is increased and it is benefit to our optimization based on mean square errors. Then the result is multiplied by 10 to observe easily. Finally, an opposite number is chosen and that means that a smaller mean square error has a bigger SNR as (16). We can find that SNR can expand the range of smaller values with logarithmic relation and describe the result more clearly. In another way, SNR is more sensitive to a small value and that is what we care; thus we take SNR instead of Y. From the explanation above, a better result has a bigger SNR as to (16).The SNR, a function of noise, is a good performance measure for robustness, is a good performance measure for robustness. Thus, it is expected that the results of statistically significant should represent the optimal combination of kinematic parameters [28]. The η is used to evaluate the fluctuation of mean square error about combination parameters in the article.Values ofY and η from 180 data points on three sections of each model are calculated. 9 cases in Table 7 are calculated and the results are recorded in Tables 8 and 9. Table 8 shows the result with velocity while Table 9 with pressure.Table 8
The results about velocity.
case A(degree) B(m) C(m) Y η 1 45 1.267 1.220 3.000 -9.542 2 45 0.725 1.020 1.052 -0.440 3 45 0.200 0.820 1.895 -5.552 4 22.5 1.267 1.020 1.338 -2.529 5 22.5 0.725 0.820 0.677 3.388 6 22.5 0.200 1.220 0.121 18.344 7 0 1.267 0.820 1.687 -4.542 8 0 0.725 1.020 0.181 14.846 9 0 0.200 1.220 0.790 2.047Table 9
The results about pressure.
case A(degree) B(m) C(m) Y η 1 45 1.267 1.220 38.132 -31.626 2 45 0.725 1.020 16.057 -24.113 3 45 0.200 0.820 43.186 -32.707 4 22.5 1.267 1.020 3.895 -11.810 5 22.5 0.725 0.820 40.566 -32.163 6 22.5 0.200 1.220 2.042 -6.201 7 0 1.267 0.820 55.801 -34.933 8 0 0.725 1.020 3.405 -10.642 9 0 0.200 1.220 31.813 -30.052Then the contribution rate is expressed as (17). It reflects the importance of different factors to η. z is the zth control factors, namely, the factors in Table 4. ηmaxz-ηminz is the range of different factors. m is the number of levels of a factor. Contribution rates are recorded in Tables 10 and 11.(17)CRz=ηmax⁡z-ηmin⁡z∑z=1mηmax⁡z-ηmin⁡z×100%Table 10
The contribution of factors based on velocity.
A B C SNR 1 -5.178 -5.528 3.616 2 6.401 5.931 3.959 3 4.117 4.947 -2.235  Range 11.579 11.459 6.194  Contribution Rate 39.611% 39.200% 21.189%Table 11
The contribution of factors based on pressure.
A B C SNR 1 -29.482 -26.123 -22.626 2 -16.725 -22.306 -15.522 3 -25.209 -22.987 -33.267  Range -12.757 -0.681 -17.745  Contribution Rate 40.910% 2.184% 56.906%From Table10, A and B play important roles in determine the uniform of velocity, because they both have a contribution over 39%. The influence of C is much smaller. From Table 11, it can be found that C, over 56%, is considered to be the most important factor to distribution of pressure and A is less important with 41%, while factor B has almost no impact on pressure with 2.184%. In the work, we must select the factors with contribution rate.In order to obtain better parameters of the separator, both distributions are considered. The values ofη according to factors are shown in Figures 3 and 4. Target values of η are bigger ones, with contribution rate; the best scheme is A2B2C2 according to Figure 3. As to the effective of pressure, A2B2C2 is still the best scheme. Although it comes from different physical quantities, the result is the same.Figure 3
SNR according to velocity.Figure 4
SNR according to pressure.It can be seen that there is a positive correlation between the distribution of velocity and pressure. A better scheme can be obtained with Taguchi method to get a uniform distribution of fluid filed in inner part.


### 3.4. Result

From the result above, A2B2C2 is considered as the best scheme. The mean square errors of velocity and pressure are shown in Figures5 and 6, respectively.Figure 5
The mean square error of velocity in all cases.Figure 6
The mean square error of pressure in all cases.The mean square error of the 10th scheme is close to the optimal, as well as the case of 6th scheme and the case of 8th scheme. The most obvious characteristic of the three cases is the smaller angle between import and horizontal. Therefore, control factor A has great influence for the flow field uniformity in the filter sections. We can draw the conclusion that the Taguchi method can also evaluate the importance of a single factor.The Taguchi methods can not only evaluate the importance of elements and also evaluate the rate of contribution of the single factor. This method applied in structure design can be refined and extended in a follow-up study such as increasing the number of parameters and the degree of level factors. A better scheme can be got to reduce the cost of design and manufacture through the analysis and calculation.


## 3.1. The Theory of Taguchi Method

The basic principle of the Taguchi method, avoiding the blindness of traverse type analysis, is the inspection of the effectiveness from the experimental variables which is the synthesis procedure conducted by orthogonal array with minimum number of experiments [15, 26]. This approach can obtain effectively the contribution rate and weight relations from various elements of target. The method can be used to realize the optimization and reduce test times compared to exhaustion method. The main steps of structure optimization for separator are as follows: (1) Determining the key elements of the research object, select all elements of the plan. (2) Design the orthogonal table for all combination schemes according to the relationship between elements. (3) Select the orthogonal schemes and calculate each of them. (4) According to the results, get the order of weight for all selected elements based on the contribution rate.


## 3.2. The Design of Orthogonal Table

In the study based on Taguchi method, to represent three parameters, we take the angle between the inlet axis and the horizontal direction as factor A, the distance between the inlet and bottom of filter as factor B, and the distance of two filters in direction-Z as factor C. Control factors and its levels are designed and shown in Table5.Table 5
The statistics of element.
Factor Level 1 Level 2 Level 3 A (degree) 45 22.5 0 B (m) 1.267 0.725 0.200 C (m) 1.220 1.020 0.820Factor A is often between 0° to 45° in actual equipment, thus 0°, 45°, and 22.5° are chosen. As to factor B and factor C, the extreme values and their medium value are taken into consideration, ranging from the highest and lowest position. Different values are noted as different levels in Taguchi method as shown in Table5.Table6 is the orthogonal table for Taguchi analysis with factors A, B, and C. A correct orthogonal table needs to fulfill with following two conditions: (1) In a column, the numbers of different values should be the same. (2) The same levels have the same frequency in every two columns.Table 6
The proof of the orthogonal.
case A B C AB BC CA 1 + + + + + + 2 + - - - + - 3 + 0 0 0 0 0 4 - + - - - + 5 - - 0 + 0 0 6 - 0 + 0 0 - 7 0 + 0 0 0 0 8 0 - + 0 - 0 9 0 0 - 0 0 0 Total 0 0 0 0 0 0L9 orthogonal array is used here; it contains nine tests, three levels, and three variables (L means Latin square and 9 means the number of cases). The orthogonal table is recorded in Table 7 and it can be used as experiment scheme.Table 7
The orthogonal table of separator.
case A B C 1 A1 B1 C1 2 A1 B2 C2 3 A1 B3 C3 4 A2 B1 C2 5 A2 B2 C3 6 A2 B3 C1 7 A3 B1 C3 8 A3 B2 C2 9 A3 B3 C1 For example, A1 stands for factor A with Level 1.


## 3.3. The Analysis of the SNR

In Taguchi method, the ratio of signal to noise (SNR) is always used as the key value to evaluate the performance of each experiment in Taguchi analysis. In this study, the target is to find smaller mean square errors; thus SNR can be expressed as (16). A bigger value of SNR is needed since it is negatively correlated with Y [27].(16)η=-10lg1n∑i=1nYi2η represents the SNR.n is the number of replications. Since it is expected to be smaller that both the mean square errors of velocity and pressure, a smaller value is what we needed. In Taguchi method, the ratio of signal to noise (SNR) is used instead of mean square errors as SNR can expand the range of smaller values. In (16), Yi2 is used to eliminate the negative values and an average is calculated. A logarithmic relation can amplify the average of Yi2 when it is under 1; this means the contrast ratio for smaller values is increased and it is benefit to our optimization based on mean square errors. Then the result is multiplied by 10 to observe easily. Finally, an opposite number is chosen and that means that a smaller mean square error has a bigger SNR as (16). We can find that SNR can expand the range of smaller values with logarithmic relation and describe the result more clearly. In another way, SNR is more sensitive to a small value and that is what we care; thus we take SNR instead of Y. From the explanation above, a better result has a bigger SNR as to (16).The SNR, a function of noise, is a good performance measure for robustness, is a good performance measure for robustness. Thus, it is expected that the results of statistically significant should represent the optimal combination of kinematic parameters [28]. The η is used to evaluate the fluctuation of mean square error about combination parameters in the article.Values ofY and η from 180 data points on three sections of each model are calculated. 9 cases in Table 7 are calculated and the results are recorded in Tables 8 and 9. Table 8 shows the result with velocity while Table 9 with pressure.Table 8
The results about velocity.
case A(degree) B(m) C(m) Y η 1 45 1.267 1.220 3.000 -9.542 2 45 0.725 1.020 1.052 -0.440 3 45 0.200 0.820 1.895 -5.552 4 22.5 1.267 1.020 1.338 -2.529 5 22.5 0.725 0.820 0.677 3.388 6 22.5 0.200 1.220 0.121 18.344 7 0 1.267 0.820 1.687 -4.542 8 0 0.725 1.020 0.181 14.846 9 0 0.200 1.220 0.790 2.047Table 9
The results about pressure.
case A(degree) B(m) C(m) Y η 1 45 1.267 1.220 38.132 -31.626 2 45 0.725 1.020 16.057 -24.113 3 45 0.200 0.820 43.186 -32.707 4 22.5 1.267 1.020 3.895 -11.810 5 22.5 0.725 0.820 40.566 -32.163 6 22.5 0.200 1.220 2.042 -6.201 7 0 1.267 0.820 55.801 -34.933 8 0 0.725 1.020 3.405 -10.642 9 0 0.200 1.220 31.813 -30.052Then the contribution rate is expressed as (17). It reflects the importance of different factors to η. z is the zth control factors, namely, the factors in Table 4. ηmaxz-ηminz is the range of different factors. m is the number of levels of a factor. Contribution rates are recorded in Tables 10 and 11.(17)CRz=ηmax⁡z-ηmin⁡z∑z=1mηmax⁡z-ηmin⁡z×100%Table 10
The contribution of factors based on velocity.
A B C SNR 1 -5.178 -5.528 3.616 2 6.401 5.931 3.959 3 4.117 4.947 -2.235  Range 11.579 11.459 6.194  Contribution Rate 39.611% 39.200% 21.189%Table 11
The contribution of factors based on pressure.
A B C SNR 1 -29.482 -26.123 -22.626 2 -16.725 -22.306 -15.522 3 -25.209 -22.987 -33.267  Range -12.757 -0.681 -17.745  Contribution Rate 40.910% 2.184% 56.906%From Table10, A and B play important roles in determine the uniform of velocity, because they both have a contribution over 39%. The influence of C is much smaller. From Table 11, it can be found that C, over 56%, is considered to be the most important factor to distribution of pressure and A is less important with 41%, while factor B has almost no impact on pressure with 2.184%. In the work, we must select the factors with contribution rate.In order to obtain better parameters of the separator, both distributions are considered. The values ofη according to factors are shown in Figures 3 and 4. Target values of η are bigger ones, with contribution rate; the best scheme is A2B2C2 according to Figure 3. As to the effective of pressure, A2B2C2 is still the best scheme. Although it comes from different physical quantities, the result is the same.Figure 3
SNR according to velocity.Figure 4
SNR according to pressure.It can be seen that there is a positive correlation between the distribution of velocity and pressure. A better scheme can be obtained with Taguchi method to get a uniform distribution of fluid filed in inner part.


## 3.4. Result

From the result above, A2B2C2 is considered as the best scheme. The mean square errors of velocity and pressure are shown in Figures5 and 6, respectively.Figure 5
The mean square error of velocity in all cases.Figure 6
The mean square error of pressure in all cases.The mean square error of the 10th scheme is close to the optimal, as well as the case of 6th scheme and the case of 8th scheme. The most obvious characteristic of the three cases is the smaller angle between import and horizontal. Therefore, control factor A has great influence for the flow field uniformity in the filter sections. We can draw the conclusion that the Taguchi method can also evaluate the importance of a single factor.The Taguchi methods can not only evaluate the importance of elements and also evaluate the rate of contribution of the single factor. This method applied in structure design can be refined and extended in a follow-up study such as increasing the number of parameters and the degree of level factors. A better scheme can be got to reduce the cost of design and manufacture through the analysis and calculation.


## 4. Conclusion

In this paper, a separator is simulated with CFD, and key parameters of a separator are optimized with Taguchi method. The results in this paper are as follows:( 1 ) The process of movement and diffusion of particles is obtained. The accumulation areas of particle are found, and the causes of accumulation are analyzed through the distribution figures. An evaluation for the structure of the original separator is taken.( 2 ) Taguchi method is introduced to take orthogonal experiment. It can reduce 27 schemes to 9. Then the weights of key parameters are determined through the contribution rate and it can be known which parameter has a bigger influence on separator( 3 ) The optimal scheme, namely, A2B2C2, is determined according to the results of a bigger SNR value. It is verified by numerical simulation and found that mean square error of distributions from both velocity and pressure is close to the minimum. The optimization of the structure is realized.In this paper a theoretical foundation and effective means has been provided and help in extending the knowledge on the research of optimization design in separator.



---

*Source: 7439245-2018-11-27.xml*