Fluid Dynamics & Materials Processing

DOI: 10.32604/fdmp.2021.012114

ARTICLE

Numerical Simulation of the Mixing and Hydrodynamics of Asphalt and Rubber in a Stirred Tank

Zechen Yao1, Renfeng Yang1,*and Haichao An2

1School of Construction Machinery, Chang’an University, Xi’an, 710064, China
2School of Transportation and Logistics, Taiyuan University of Science and Technology, Taiyuan, 030024, China
*Corresponding Author: Renfeng Yang. Email: yrfeng@chd.edu.cn
Received: 15 July 2020; Accepted: 24 December 2020

1  Introduction

At present, nearly 1.5 billion waste tires are generated around the world each year and the number is still increasing rapidly [1]. Waste tires are made of long polymer chains, most of which are crosslinked with sulfur bridges; it is very difficult to recover and recycle waste tires [2,3]. The disposal of used tires is a worldwide problem and burying or burning is the general treatment; hence, waste tires are called “black pollution” [4]. Asphalt-Rubber (AR) is a new material which is made by adding crumb waste tire rubber into base asphalt. This can reduce the amount of waste tires and it can also improve the road performance of asphalt by using natural rubber and carbon black from waste tires. Hence, AR technology is a green road construction and maintenance technology. With the adoption of the Paris climate agreement [5], a great number of countries are paying more attention to environmental protection, and AR technology has received more and more attention.

At present, there has been a lot of research on the modification effect of AR binders [6–10], which has resulted in remarkable achievements and promoted the development of AR technology. However, there have been few studies on the production process of AR binders, and uneven mixing phenomena always appear in AR production equipment. The AR binder is a solid-liquid material, and in recent years, there have been a lot of studies on the solid-liquid two-phase mixing process worldwide [11–18]. Nevertheless, a majority of the existing studies on the two-phase flow adopted water as the liquid phase and ignored the impact of the liquid’s viscosity, and the studies on the mixing process of AR binders with high viscosity are still basically missing.

CFD software has been proved to be very helpful for engineering design optimization [19–22]. Exploring the suspension quality of solid particles in the flow field and obtaining their solid-liquid suspension characteristics are very important for the design optimization and industrialization of stirred tanks. In the last two decades, CFD has been considered as a useful tool for simulating the fluid dynamics of two-phase flow in the stirred tank [14,16,23–25]. For instance, Hosseini et al. [26] predicted the normalized solid concentration profiles, cloud height, homogeneity based on CFD. Qi et al. [27] studied the impacts of particle density, particle size, liquid viscosity on mixing uniformity in a stirred tank with Smith turbine based on a CFD model. Gu et al. [25] simulated the effect of impeller type, impeller speed, impeller blades aperture ratio, particle size, and liquid viscosity on the suspension quality of solid particles in the stirred tank. Till et al. [28] explored the impact of the number of impeller blades on the mixing efficiency using a CFD model. Vilardi et al. [29] optimize the production of metallic iron nanoparticles in a baffled stirred tank reactor using CFD model. From this, it can be seen that CFD software is an effective method to simulate the suspension characteristics of solid particles.

In this paper, the suspension characteristics of asphalt and CR particles in a stirred tank with a Rushton impeller, a PBTD-6 impeller and a TXL impeller have been investigated by CFD simulation. The liquid velocity field, the axial CR concentration profiles, the contour plots of CR concentration and the cloud height in the stirred tank have been predicted. And the effects of impeller speeds, impeller types, CR particle sizes and initial CR particle loading on the suspension quality of CR particles in a stirred tank have been deeply investigated. The improved sample withdrawal method was used to verify the correctness of numerical modeling.

2  Physical Model and Experimental Methodology

2.1 Experimental Apparatus

As shown in Fig. 1, the experimental system consisted of a stirred tank with a diameter T = 0.3 m, equipped with four baffles with a width of W = T/10. The height (H) of liquid level was 0.3 m. The impeller’s off-bottom clearance was set at T/3. The impeller was driven by a motor (OuBang Company, China) with a power output of 250 W. The outer wall of the stirred tank was wrapped in a heating device with a power output of 7.2 KW, and a temperature controller (FengDa Company, China) was used to ensure a constant temperature of 180°C. The heating device was surrounded by a thermal insulation layer to reduce heat loss. A Rushton, a PBTD-6 and a TXL, as shown in Fig. 2, were adopted in the research; the diameter of these three impellers was d = T/3. The impeller’s power consumption was measured by a torque transducer (JinNuo Company, Model: JN-DN, China).

Figure 1: Schematic of the stirred tank

Figure 2: The structural dimensions of three impellers: (a) Rushton; (b) PBTD-6; (c) TXL

2.2 Materials and Methods

The dispersion of CR (ρs = 1170 kg/m3, particle diameter = 550 μm, particle shape is nearly spherical, Xi’an Zhongxuan Company, Xi’an, China) in asphalt (ρ1 = 900 kg/m3, μ = 1.5 Pa·s, SK Company, Seoul, Korea.) was investigated at 180°C. The CR particles were selected due to the difference of density. The asphalt is Newtonian fluid with constant viscosity at 180°C. The viscosity was tested by Rhein viscometer (Japan) following GB/T 22235-2008 [30]. The A mean CR loading of 20 vol%, which was always used in the preparation of the AR binder, was used in both the experiments and the simulations. The industrial production of AR binders was always divided into two steps: Pre-mixing and maturing mixing (solid-liquid suspension process). Hence, in order to be consistent with the reality, before the suspension experiment, the AR binders were pre-mixed for 5 min at 1500 rpm to ensure even mixing, and in order to avoid inadequate settlement caused by the lower density difference, mixing time of 500 s was used for the suspension experiments. At present, the sample withdrawal method [31,32] is always used to verify the accuracy of the CFD model. However, due to the sharp fall in the temperature of the AR binder when the samples were taken out and the density of asphalt being obviously changed, the conventional sample withdrawal method makes it difficult to accurately measure the content of the CR in the asphalt. In order to accurately measure the CR content, trichloroethylene and petrol were used to dissolve the asphalt in the samples. Then the CR particles were filtered out using a steel net with 300 mesh (50 μm) and dried in an oven, and the mass of the CR particles was then weighed. Finally, the volume fraction of the CR in the samples was obtained. The samples were withdrawn at seven axial positions of the stirred tank where z/H = 0.125, 0.25, 0.375, 0.5, 0.625, 0.75 and 0.875 in this paper. 100 g of AR binder was taken at every sampling position in this paper. It is worth mentioning that a large amount of trichloroethylene and petrol were needed to dissolve asphalt. Here we can refer to the literature Asphalt Solubility Test [33]. Moreover, due to trichloroethylene being a volatile carcinogen, gas masks (3M Company, Model: 3200, China) were worn during the experiments.

3  Numerical Model and Methods

3.1 CFD Model

In this paper, the solid-liquid flow was modeled based on Eulerian-Eulerian multifluid model. In this model, the phases were seen as continua, interpenetrating and interacting with each other in the flow field [15,18]. The governing equations of which are as follows:

Continuity equations [15]:

∂∂t(αiρi)+∇⋅(αiρivi→)=0,∑αi=1 (1)

Momentum equations [15]:

∂∂t(αiρivi→)+∇(αiρivi→vi→)=−αi∇p+∇τi¯¯+∑R→ij+αiρiFi→ (2)

where,

τi¯¯=αiμi(∇vi→+∇vi→T)+αi(λi−23μi)∇vi→I¯¯ (3)

where, i denotes the phase, i = s and i = l represent solid phase and liquid phase, respectively; αi is the volume fraction of each phase; ρi is the density of each phase; vi→ is the velocity vector of each phase; p is the pressure of all the phases; Rij→ is the interaction force between the phases; Fi→ is the external body force; τi¯¯ is the stress-strain tensor; μi is the shear viscosity of each phase; λi is the bulk viscosity of each phase and I¯¯ is the unit stress tensor.

The following formula is used to calculate the interphase force [18]:

R→ij=∑Kij(v→i−v→j);Ksl=αsρsfDτs;τs=ρsdp218μl (4)

where, Kij is the momentum exchange coefficient between the phases; fD is the function of the drag force; τs is the particles’ relaxation time and dp is the solid particles’ diameter. In the commonly used drag model, the Gidaspow model is combined with the Ergun formula on the basis of the Wen-Yu model, and the ensuing calculation is more accurate. The exchange coefficient was calculated as follows [11,34]:

Ksl=34CDαsαlρl|v→s−v→l|dpαl−2.65 (5)

where, CD is the drag coefficient [11,34]:

CD=24αlRes[1+0.15(αlRes)0.687] (6)

The relative Reynolds number is defined as [18,24]:

Res=ρldp|v→s−v→l|μl (7)

According to the definition of the Re number (Re = ρNd2/(60μ)), the maximum value was Re = 42 in the mixing system, which indicates that laminar flow occurred in the stirred tank.

3.2 Methodology and Boundary Conditions

In this paper, the MRF technology was employed to simulate the rotation zone with impeller. As shown in Fig. 3, the unstructured tetrahedral grid was used in rotation zone due to the complexity of impeller structure, while the structured hexahedral grid was used in non-rotation zone to improve the accuracy and convergence of calculation. Some models are partitioned using unstructured grids due to the convergence. A grid independence test was carried out to avoid the influence of the grid number on the accuracy of the calculation. Taking the Rushton impeller system as an example, the variation of the velocity magnitude with z/H (z being the Z axis coordinate) has been shown in Fig. 4 for the division of three grid numbers. When the number of grids was increased from 682301 to 1325864, the relative deviation was less than 3%. Within the allowable range of the project, the number of grids used for the Rushton impeller system was 682301. Similarly, the number of grids used for the PBTD-6 system and the TXL system were 623578 and 685264, respectively. Fig. 5 has shown the comparison of the local axial CR concentration profiles both with the temperature boundary condition (180°C) and without; this has shown that the temperature had little impact on the distribution of particles. For the sake of simplicity, the temperature boundary conditions were not considered in the CFD model. Considering the impact of steady-state MRF simulations of the two-phase system on the model [35], the average CR concentration computed by the model is 20 vol%.

Figure 3: Grid model of the stirred tank: (a) Rushton; (b) PBTD-6; (c) TXL; (d) Stirred tank

Figure 4: Grid independence verification (particle diameter = 550 μm, r/R = 0.4, N = 60 rpm)

Figure 5: Temperature boundary condition verification (particle diameter = 550 μm, r/R = 0.8, N = 60 rpm)

The ANSYS FLUENT 18.0 was used to calculate simulation model. Asphalt is the incompressible fluid. All terms of governing equations were discretized by a first-order upwind and the coupling of the pressure and the velocity was carried out using the SIMPLEC algorithm [24]. The steady-state model was used for calculation. The relative residual was set at 1 × 10-4, which was considered to be the convergence index.

4  Results and Discussion

4.1 Verification of the Model

In order to validate the correctness of numerical simulation, the results of computation were validated by experimental measurements. The improved sample withdrawal method was employed to measure the local CR concentration. Fig. 6 has shown the variation trend of the local axial CR concentration profiles between the simulation and the experimental results were similar although there are some deviations. Therefore, the correctness of CFD model is proved.

Figure 6: Verification of the model: Comparison of the local axial CR concentration profiles between the simulation and the experiment (particle diameter = 550 μm, r/R = 0.8, N = 300 rpm)

4.2 Distribution of the Liquid Velocity Field

Fig. 7 has illustrated the two-dimensional liquid velocity contour for the Rushton impeller, the PBTD-6 impeller and the TXL impeller. It can be found that a typical two-loop flow pattern was obtained when the Rushton impeller and the PBTD-6 impeller were adopted in the stirred tank. When these two impellers were employed, the fluid flowed from the impeller to the wall and divided into upper and lower streams, converged at the center of the stirred shaft, which is consistent with Liu’s research [18]. Meanwhile, it can be seen that the strength of the flow field near the bottom of the stirred tank with the Rushton and the TXL was greater than the PBTD-6 impeller system, which is beneficial to the CR particles’ suspension.

Figure 7: Liquid velocity field distribution corresponding to the different impellers (N = 60 rpm, particle diameter = 550 μm): (a) Rushton; (b) PBTD-6; (c) TXL

Fig. 8 has illustrated the axial, radial, and tangential liquid velocities on the sampling line L (shown in Fig. 1) for the Rushton impeller, the PBTD-6 impeller and the TXL impeller. These three velocities in the stirred tank with the Rushton impeller were higher than those of the other two impeller systems, which indicated that CR particles in Rushton impeller system were easy to suspend at same impeller speed. The radial velocities of the TXL impeller were lower than those of PBTD-6 impeller, while its tangential velocities were higher than those of the PBTD-6 impeller.

Figure 8: Comparison of the liquid velocity along sampling line L in Fig. 1, corresponding to the different impellers (N = 60 rpm, particle diameter = 550 μm): (a) Axial velocity; (b) Radial velocity; (c) Tangential velocity

4.3 Effects of Impeller Types and Impeller Speeds

Fig. 9 has illustrated the local axial CR concentration profiles at different impeller speeds for the Rushton impeller, the PBTD-6 impeller and the TXL impeller. Fig. 10 has shown the computed contour plots of the CR volumetric fraction on a horizontal plane at different vessel heights (i.e., z = 0H, 0.3H, 0.5H, 0.7H, 0.9H) for different impeller speeds. As displayed in Fig. 9, at a low impeller speed (i.e., N = 60 rpm), it can be observed that the Ch/Cavg values were very small at the upper part of the stirred tank, which shown that few CR particles suspended in the upper part of the stirred tank, and the similar phenomenon can be also seen in Fig. 10. With the increase of the impeller speed, plenty of CR particles was dispersed from the lower part of the stirred tank and the degree of homogeneity was improved. It can be seen that the Ch/Cavg values of the Rushton system at the upper part were larger than that of the PBTD-6 impeller system and the TXL impeller system at N = 180, 300, 420 rpm. As shown in Fig. 10, the CR particles’ suspension quality for the PBTD-6 impeller system was the worst of these three impeller systems. The cloud height was defined as the distance from the highest point of the iso-surface of the average CR concentration to the bottom of the stirred tank [13], which can reflect the quality of the suspension. The cloud heights of the CR particles in the stirred tank at different impeller speeds for the Rushton impeller, the PBTD-6 impeller and the TXL impeller have been illustrated in Fig. 11. The cloud heights increased with the impeller speed and it can be observed that the cloud height of the PBTD-6 system was lower than those of the other two impeller systems, indicating that its mixing quality was not good, which is consistent with the CR concentration distribution in Figs. 9, and 10. There was little difference between the two cloud heights of the Rushton system and the TXL system, and the cloud heights of the Rushton impeller system was slightly lower than those of the TXL system at N = 300, 420 rpm. The reason that the suspension quality of CR particles in a stirred tank with Rushton impeller is better may be the normal direction of the blade is the same as the linear velocity direction of the impeller, and the Rushton impeller has the higher fluid velocity, which leads to better particles suspension under same impeller speed.

Figure 9: Effects of impeller speeds on the local axial CR concentration profiles (particle diameter = 550 μm, r/R = 0.8): (a) N = 60 rpm; (b) N = 180 rpm; (c) N = 300 rpm; (d) N = 420 rpm

Figure 10: Contour plots of the CR volumetric fraction on a horizontal plane at different vessel heights (i.e., z = 0H, 0.3H, 0.5H, 0.7H, 0.9H) for different impeller speeds (particle diameter = 550 μm): (a) Rushton; (b) PBTD-6; (c) TXL

Figure 11: Cloud heights with Cavg = 0.2 at different impeller speeds (particle diameter = 550 μm)

4.4 Effects of CR Particle Diameter

In practical pavement engineering, #20 mesh (830 μm), #30 mesh (550 μm) and #40 mesh (350 μm) CR is often used to prepare AR binders, and particle size is one of the key factors affecting suspension quality. Fig. 12 shows the local axial CR concentration profiles effected by the CR particle diameter. It can be observed that the suspension quality decreased with the increase of the CR’s particle diameter. This can be interpreted by the free setting velocity (Vt) [15]:

Vt=[4gdp(ρs−ρl)3CDρl]0.5 (8)

where, g is the acceleration due to gravity; dp is the particle’s diameter; ρs is the CR’s density; ρl is the liquid’s density and CD is the drag coefficient. The settling velocity increases with the particle diameter, which means that larger particles need more energy to suspend in the stirred tank. Hence, this suggests that the impeller speed should be increased to ensure mixing uniformity when AR binders with a lower mesh number are produced.

Figure 12: Effects of particle diameter on axial CR concentration profiles (Rushton impeller, N = 60 rpm, r/R = 0.8)

4.5 Effects of the Initial CR Particle Loading

Fig. 13 has shown the effects of the initial CR particle loading on the local axial CR concentration profiles. The suspension quality of CR particles in the stirred tank with a lower initial CR loading of 15 vol% was better than that of 20 vol% and 25 vol% because CR particles in the mixing system with higher content have to gain more energy to overcome the gravity to suspend.

Figure 13: Effects of the initial CR particle loading on the axial CR concentration profiles (Rushton impeller, N = 60 rpm, particle diameter = 550 μm, r/R = 0.8)

4.6 Evaluation of the CR Particle Suspension Quality

In this paper, the relative standard deviation (RSD) [24] was employed to illustrate the variation of suspension quality with power consumption in the stirred tank. The RSD was calculated by the local CR concentration at different positions in the stirred tank, where z/H = 0.125, 0.25, 0.375, 0.5, 0.625, 0.75 and 0.875, and the smaller RSD value represented the better suspension quality in stirred tank. The RSD was defined [24]:

RSD=1Cavg1n−1∑h=1n(Ch−Cavg)2 (9)

where, n is the number of sampling points, Ch is the local axial CR concentration at height h, and Cavg is the average CR concentration. Fig. 14 has illustrated the variation of RSD with the power consumption (P) for the Rushton impeller, the PBTD-6 impeller and the TXL impeller for experiments. As can be seen that the RSD decreased as the power consumption increased, which indicated that the increase of the power consumption was beneficial to improve the uniformity of CR particles suspension. It can also be found that the suspension quality of the TXL system was better than that of the Rushton and the PBTD-6 impeller system, and the mixing effect of the PBTD-6 impeller was the worst. In other words, the power consumption of the TXL impeller with regard to the just off-bottom suspension state was the lowest of these three impellers. The impeller speeds at each constant power consumption in experiments for Fig. 14 are shown in Tab. 1.

Figure 14: Comparison of the RSD at different power consumptions in experiments

Table 1: Impeller speeds at each constant power consumption in experiments

5  Conclusions

The CR particles suspension characteristic in asphalt in a stirred tank with a Rushton impeller, a PBTD-6 impeller and a TXL impeller were explored by numerical simulation in this paper. The effects of impeller types, impeller speeds, CR particle diameter and initial CR particle loading on the flow field characteristics and the suspension quality were deeply studied. Results showed that the local axial CR concentration profiles calculated by the CFD model were similar with the experimental results. The average axial, radial, and tangential liquid velocities in the lower part of the stirred tank with the Rushton impeller were the highest of these three impeller systems at the same impeller speed. As expected, the suspension quality of CR particles was increased with impeller speed. At the same speed, the CR particles’ suspension quality in the stirred tank with the PBTD-6 was the worst of the three impeller systems. A larger particle size and higher initial CR particle loading were not conducive to the suspension of CR particles. Moreover, the TXL impeller was able to enhance the degree of homogeneity compared with the PBTD-6 impeller and the Rushton impeller at the same power consumption. This paper may help to improve the production quality of AR binders and extend the implication of two-phase fluid dynamics on mixing system with high temperature and high viscosity.

However, there are some limitations in the experiments. Due to the high temperature of experimental process may easily lead to the breakdown of motor and torque transducer and the impact of the reaction of asphalt and CR particles on accuracy of extracting solid particles in asphalt, finding the substitute whose physical property is similar with asphalt to carry out experiments at room temperature may be good for the research of AR experiment in the future.

Acknowledgement: The authors gratefully acknowledge the helpful comments and suggestions from several reviewers, which significantly improve the earlier versions of this paper.

Funding Statement: This work was supported by the Fundamental Research Funds for the Central Universities [300102250112]; the Fundamental Research Funds for the Central Universities, CHD [300102251711].

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.