Effect of Geometrical Parameters on Extraction Efﬁciency of the Annular Centrifugal Contactor

: The geometrical parameters of annular centrifugal contactors (ACCs) have an important inﬂuence on the extraction efﬁciency. The present work used a home-made 25 mm ACC constructed by 3D printing to investigate the effect of ﬁve geometrical parameters on the extraction efﬁciency. These parameters are annular width ( d ), clearance height ( H c ), rotor inlet diameter ( D in ), bottom vane number ( N ), and the bottom vane’s bending direction ( S ). Central composite design was employed to design the experiment, and the response surface methodology was used to analyze the data. The results show that H c and D in were positive for efﬁciency, while d and N were negative. When the bottom vane’s bending direction was the same as the liquid helical ﬂow direction, the efﬁciency improved compared to the straight vane. It is found that 3 mm d , 5 mm H c , 6 mm Din , and four clockwise covered vanes are the parameters where the efﬁciency reached the highest point of 94.5%. We analyzed the interactions between the parameters based on the coefﬁcients of the quadratic equation, and the interactions were not considered in previous studies. This work surprisingly reveals that the effects of the parameters on the extraction efﬁciency were not independent, and there were interactions between the parameters. The interaction between the rotor inlet diameter and annular width was signiﬁcant and could not be ignored. These results could serve as a reference for optimizing extraction processes and the design of ACCs.


Introduction
Extraction is one of the essential separation techniques in chemical processes. A component is separated from a solvent by transferring the component to an extractant which is immiscible to the solvent. This technique is conducted according to the different solubility of the components in the solvent and the extractant. Mixer-settlers, pulse extraction columns, and annular centrifugal contactors (ACCs) are common instruments used in liquid-liquid extraction [1,2]. The first two types of equipment utilize gravity for separation and result in a slow separation of the light phase from the heavy phase, while an ACC is an efficient device that utilizes centrifugal force. An ACC has the advantages of a small equipment size, a low liquid holdup, a short residence time, and good phase separation performance [3]. Therefore, ACCs have already been used in the nuclear industry, metallurgy, the chemical industry, and biological engineering [4][5][6][7].
The extraction efficiency is the main indicator to evaluate the mass transfer performance of ACCs. Different flow patterns in the ACC lead to different efficiencies. Available studies show that the extraction efficiency could be impacted by the operation parameters, geometrical parameters, and material system [8][9][10]. There have already been many works that studied the effect of operation parameters (total flow rate, phase ratio of two phases, and rotor speed) on efficiency. However, few researchers studied the effect of geometrical parameters, which include annular width (d), clearance height (H c ), rotor inlet Zhao et al. [11] found that the extraction efficiency increased with the annular width and clearance height. Xu et al. [12] showed that the smaller the annular width, the higher the efficiency. Birdwell et al. [13] obtained a lower efficiency when using a curved vane instead of a straight vane, while Leonard et al. [14] and Wardle et al. [15] illustrated that a curved vane could provide higher efficiency. Absolutely, these findings mentioned above have notable differences, and there are no studies that provide a reasonable explanation for the appearance of these differences.
We summarized the experimental design methods of various researchers. It was found that all existing works only study the effect of a single factor on the extraction efficiency once, while the interactions between parameters have never been considered. Therefore, it is necessary to consider the effects of multi-geometrical parameters on the efficiency simultaneously in the experimental design and to explore the influence pattern of geometrical parameters on efficiency. We employed the response surface methodology (RSM) to analyze the interactions between parameters. RSM uses a multivariate quadratic function (MQF) to represent the relationship between variables and response, while the interaction term coefficients characterize the degree of interaction between variables. RSM is widely applied to study extraction processes for its intuitive response surface plots [16,17]. The analysis steps of RSM can be referenced in Silva's work [18].
This work studied the effect of d, H c , and D in on the extraction efficiency and investigated the interaction between the three parameters. The experimental system was cyclohexane-isopropanol-water, and central composite design (CCD) was selected to optimize the experimental design. Moreover, this work compared the extraction efficiency obtained from the different bottom vanes (N and S). We analyzed how the geometrical parameters influence the extraction efficiency and provide a reasonable explanation for the discrepancies in the literature above. Our research is of great importance for the design and optimization of ACCs.

Chemicals
Isopropanol (>98%) and cyclohexane (>98%) were obtained from Jiangtian Chemical Technology Co., Ltd., Tianjin, China. The solution to be extracted was ten wt.% isopropanolcyclohexane, and the extractant was water. During the process, the isopropanol in the cyclohexane phase was extracted into the aqueous phase. Figure 1 shows a schematic diagram of the ACC. The ACC was made by a UV-curable 3D printer (bene 3) from Nova Smart Technology Co., Ltd., Shenzhen, China. PX-88 resin (PMMA analogue crystalclear UV resin) from Zhemo Technology Co., Ltd., Shenzhen, China, was used as the print material. PX-88 resin has a good solvent resistance and could prevent the ACC from deforming during the experiment. In addition, its transparency enables us to measure the liquid height in the annular region. The dimensions of the homemade ACC are listed below: rotor's inner diameter was 25.0 mm; rotor's outer diameter was 28.2 mm; light/heavy phase inlet height was 44.5 mm; separation zone length was 60.0 mm; heavy phase weir diameter was 13.4 mm; light phase weir diameter was 12.0 mm.

Experimental Design
Design Expert V12 developed by Stat-Ease, Inc., was employed as the design tool. We designed a 3-factor, 5-level CCD experiment for continuous parameters (d, H c , and D in ). Table 1 shows the correspondence between the CCD's coding levels and the values of parameters. In addition, the discontinuous parameters (N and S) were experimented with individually. of parameters. In addition, the discontinuous parameters (N and S) were experimented with individually.    Figure 2 shows the schematic diagram of the experimental apparatus. Firstly, the heavy phase (water) was pumped into the ACC at 60 mL·min −1 . Secondly, the light phase (10 wt% isopropanol-cyclohexane) was pumped in with 30 mL·min −1 after the liquid flowed out of the heavy phase outlet. At 6 min after starting the experiment, the isopropanol concentration in the light phase outlet stream did not change, and we considered the process had reached a steady state. This measurement was consistent with the conclusion mentioned by Kadam et al. [19] that the extraction process would obtain the steady state after three times the residence time. Thirdly, after reaching the steady state, we recorded the annular liquid height (ALH) and analyzed the composition of the light phase outlet product. Finally, we kept the motor running and drained the annular liquid from the bottom discharge port, and its volume was the holdup volume (V m ) in the annular region. Each experiment was repeated 3 times. The motor speed was set to 2500 r·min −1 for all experiments and was rotated in a clockwise direction.

Methods
An Agilent 6820 gas chromatograph was used with a flame ionization detector at 120 °C and an Agilent SE-54 capillary column with dimensions of 30 m × 0.32 m × 0.25 μm. Split injections were made at 100 °C with a split ratio of 50:1.
The Taylor number (Ta) is a dimensionless number proposed by Taylor and defined as the ratio of the centrifugal force to viscous force, which reflects the flow pattern of the fluid [20]. Equation (1) can calculate the Ta of fluid in the annular region, and Equations   Figure 2 shows the schematic diagram of the experimental apparatus. Firstly, the heavy phase (water) was pumped into the ACC at 60 mL·min −1 . Secondly, the light phase (10 wt% isopropanol-cyclohexane) was pumped in with 30 mL·min −1 after the liquid flowed out of the heavy phase outlet. At 6 min after starting the experiment, the isopropanol concentration in the light phase outlet stream did not change, and we considered the process had reached a steady state. This measurement was consistent with the conclusion mentioned by Kadam et al. [19] that the extraction process would obtain the steady state after three times the residence time. Thirdly, after reaching the steady state, we recorded the annular liquid height (ALH) and analyzed the composition of the light phase outlet product. Finally, we kept the motor running and drained the annular liquid from the bottom discharge port, and its volume was the holdup volume (V m ) in the annular region. Each experiment was repeated 3 times. The motor speed was set to 2500 r·min −1 for all experiments and was rotated in a clockwise direction.

Methods
An Agilent 6820 gas chromatograph was used with a flame ionization detector at 120 • C and an Agilent SE-54 capillary column with dimensions of 30 m × 0.32 m × 0.25 µm. Split injections were made at 100 • C with a split ratio of 50:1.
The Taylor number (Ta) is a dimensionless number proposed by Taylor and defined as the ratio of the centrifugal force to viscous force, which reflects the flow pattern of the fluid [20]. Equation (1) can calculate the Ta of fluid in the annular region, and Equations (2) and (3) provide the critical Taylor number (Ta Cr ) with the axial flow [21]. In this experiment, the Taylor number was 1.28 × 10 7 -5.91 × 10 7 , and the corresponding critical Taylor number was 4.83 × 10 3 -4.46 × 10 3 . Thus, the Ta is more significant than 1000 times the Ta Cr , indicating that the flow patterns in the annular region are in turbulent Taylor flow.
where Ta is the Taylor number of the fluid; Ta Cr is the critical Taylor number; r io is the outer radius of the rotor; d is the annular width; ω is the angular velocity of the rotor; ν is the kinematic viscosity of the fluid; Re z is the axial flow Reynolds number in the annular region; u is the liquid flow velocity along the axial direction in the annular region.
(2) and (3) provide the critical Taylor number (Ta Cr ) with the axial flow [21]. In this experiment, the Taylor number was 1.28 × 10 7 -5.91 × 10 7 , and the corresponding critical Taylor number was 4.83 × 10 3 -4.46 × 10 3 . Thus, the Ta is more significant than 1000 times the Ta Cr , indicating that the flow patterns in the annular region are in turbulent Taylor flow.
Ta = riod 2 ω 2 /ν 2 (1) where Ta is the Taylor number of the fluid; Ta Cr is the critical Taylor number; r io is the outer radius of the rotor; d is the annular width; ω is the angular velocity of the rotor; ν is the kinematic viscosity of the fluid; Re z is the axial flow Reynolds number in the annular region; u is the liquid flow velocity along the axial direction in the annular region.

Extraction Efficiency
The extraction efficiency is the evaluation index of the mass transfer performance of ACCs, and the higher the efficiency, the better the mass transfer performance. The extraction efficiency defined by the Murphree equation is shown in Equation (4) [22]. We took the average of three experiments as the extraction efficiency under a certain condition.
where E O is the extraction efficiency of the organic phase (%); w O,in and w O,out are the isopropanol mass fractions in the light phase inlet stream and outlet stream, respectively; w O,eq is the mass fraction of isopropanol in the organic phase when the extraction process reaches equilibrium.

Results and Discussion
Error analysis of the data obtained through three replicate experiments showed that the relative error of the extraction efficiency value was 0.920%.
We obtained Equation (5) using the multivariate quadratic function (MQF) in the response surface methodology (RSM) to fit the data of extraction efficiency. Table 2 shows the ANOVA for the Equation (5), where the p-value of the equation is less than 0.0001, which implies the equation is significant, and there is only a 0.01% chance that such an Fvalue (19.0) could occur due to noise. Furthermore, the regression coefficient of Equation (5) is 0.945, indicating that the equation fits the data well. Figure 3 shows the response surface of equation (5), where the x-y plane represents the experimental interval for any

Extraction Efficiency
The extraction efficiency is the evaluation index of the mass transfer performance of ACCs, and the higher the efficiency, the better the mass transfer performance. The extraction efficiency defined by the Murphree equation is shown in Equation (4) [22]. We took the average of three experiments as the extraction efficiency under a certain condition.
where E O is the extraction efficiency of the organic phase (%); w O,in and w O,out are the isopropanol mass fractions in the light phase inlet stream and outlet stream, respectively; w O,eq is the mass fraction of isopropanol in the organic phase when the extraction process reaches equilibrium.

Results and Discussion
Error analysis of the data obtained through three replicate experiments showed that the relative error of the extraction efficiency value was 0.920%.
We obtained Equation (5) using the multivariate quadratic function (MQF) in the response surface methodology (RSM) to fit the data of extraction efficiency. Table 2 shows the ANOVA for the Equation (5), where the p-value of the equation is less than 0.0001, which implies the equation is significant, and there is only a 0.01% chance that such an Fvalue (19.0) could occur due to noise. Furthermore, the regression coefficient of Equation (5) is 0.945, indicating that the equation fits the data well. Figure 3 shows the response surface of equation (5), where the x-y plane represents the experimental interval for any two of the three parameters, and the z-axis means the E O . The fluctuation in the response surface indicates the variation in the efficiency.
ALH (i.e., annular liquid height) and V m (i.e., holdup volume) were measured to explain the effect of geometrical parameters on E O , and the data are shown in Figures 4  and 5, respectively. The error of ALH and V m mainly comes from the measurement. In this study, the absolute error of ALH is 0.3 mm, and V m is 0.1 mL. The change in the geometrical  surface indicates the variation in the efficiency. ALH (i.e., annular liquid height) and V m (i.e., holdup volume) were measured to explain the effect of geometrical parameters on E O , and the data are shown in Figures 4 and  5, respectively. The error of ALH and V m mainly comes from the measurement. In this study, the absolute error of ALH is 0.3 mm, and V m is 0.1 mL. The change in the geometrical parameters caused a shift in the phases' flow pattern, consequently changing the ALH and V m , which caused the variability of the E O .    Figure 4a shows that ALH decreases significantly from 14.8 to 6.9 mm, and Figure 5a shows V m increases slightly from 9.8 to 10.0 mL as the annular width increases.     The reason for the change in E O with the annular width can be summarized as the following two points: (1) As the annular width increased, the V m increased. The mixing time of the annular liquid was proportional to the V m . Therefore, the mixing time became longer, and the extractant extracted more components, which improved mass transfer performance. (2) The ALH decreased when the annular width increased, leading to unit volume power (P/V) reduction. Thus, the mass transfer coefficient decreased, which was detrimental to mass transfer.

Effect of Geometrical Parameters on Extraction Efficiency
Kadam et al. [17] had given the mass transfer coefficient (k c a) of the centrifugal extraction as a function of P/V, as shown in Equation (6).
The reason for the change in E O with the annular width can be summarized as the following two points: (1) As the annular width increased, the V m increased. The mixing time of the annular liquid was proportional to the V m . Therefore, the mixing time became longer, and the extractant extracted more components, which improved mass transfer performance. (2) The ALH decreased when the annular width increased, leading to unit volume power (P/V) reduction. Thus, the mass transfer coefficient decreased, which was detrimental to mass transfer.
Kadam et al. [17] had given the mass transfer coefficient (k c a) of the centrifugal extraction as a function of P/V, as shown in Equation (6). k c a = 0.64(P/V) 0.22 (6) Arafat et al. [23] calculated the useful work (P) conducted by the rotor in Equation (7). Equation (8) could be used to calculate the rotor worked area's volume (V).
The adjustment factor j in Equation (7) can be obtained by Equations (9) and (10). j = 0.0554(lgN Re ) + 1.368 (3.0 × 10 3 < N Re < 1.0 × 10 6 ) (9) where k c a is the mass transfer coefficient; P is the useful work conducted by the rotor on the liquid in the annular region; V is the volume of fluid in the area of direct rotor work, i.e., the volume of liquid in the annular region after removing the clearance liquid at the bottom of the rotor; H is the liquid height of V in the annular region; r io is the outer radius of the rotor; j is the adjustment factor; ρ is the density of the mixed fluid; µ is the viscosity of the mixed fluid; d is the width of the annular region; N Re is the Reynolds number of the fluid in the annular region; ω is the angular velocity of the rotor.
With the annular width increased from 3.0 to 5.0 mm, the mixing time increased by 2.04%, while the k c a decreased by 12.4%, which noted that the negative impact caused by the increase in the annular width dominated the efficiency. Therefore, the E O dropped when the annular width increased. We predict through the MQF that the positive impact of the mixing time would dominate the efficiency when the d is more than 5.16 mm. However, due to the size limitation of the 3D printing equipment, the E O at larger annular widths was not investigated.
In this subsection, we explain how the geometrical parameters influence the E O from the perspective of mass transfer theory. The introduction of ALH and V m into the E O 's study could help us to understand the mass transfer process in the ACC. Figure 3a,b show that the E O reaches the maximum when the clearance height is 5 mm. The growth rate of E O with the clearance height is the same for different annular widths and rotor inlet diameters. As shown in Figures 4b and 5b, with the clearance height increasing from 3.0 to 5.0 mm, the increment of the ALH has the same value as the increment of the clearance height, and the V m increases obviously by 22.3%.

Clearance Height
As the clearance height increased, the ALH increased, but the H remained unchanged. The mass transfer coefficient remained the same according to Equations (6)- (10). Notably, the increase in V m significantly prolonged the mixing time of the two-phase liquid in the annular region and intensified the mass transfer process. Therefore, the E O increased with the clearance height.
By analyzing the V m and ALH, we found that the change in the clearance height had no effect on k c a, but could affect E O by changing V m . Therefore, a larger clearance height is beneficial to E O .

Rotor Inlet Diameter
From Figure 3c, we find that the annular width governs the effect of the rotor inlet diameter on the E O . When the annular width is 3.0 mm, the E O increases by 6.18% significantly with the rotor inlet diameter increases. However, the E O only increases by 2.99% when the annular width is 4.0 mm. The increasing speed slows down compared with 3.0 mm. E O maintains a value of 86.0% as the annular width reaches 5.0 mm. ALH and V m data are shown in Figures 4c and 5c, respectively.
The ACC draws the fluid into the rotor, relying on the rotor's self-priming force caused by the pressure differential between the rotor inlet and the cavity inside the rotor. The greater the self-priming force, the more easily annular fluid is pulled into the rotor. Hence, the ALH and V m reduce. Bernstein et al. [24] found the self-priming force of the rotor was negatively related to the rotor inlet diameter by theoretical derivation.
When the annular width was 3.0 mm, the rotor's self-priming force decreased as the rotor inlet diameter increased, which led to difficulty in drawing the two-phase liquid into the rotor. The accumulation of V m prolonged the mixing time of the two-phase liquid, and the extraction efficiency rose. However, an increase in the annular width would cause the ALH to decrease. The rotor inhaled more gas because of a lower liquid level, weakening the effect of the rotor inlet diameter on the efficiency. Therefore, the rise in E O when the annular width is 4.0 mm was less than 3.0 mm. If the annular width reached 5.0 mm, the ALH was only 6.9 mm. The rotor inhaled a large amount of gas. At this time, the pressure differential had no changes with the rotor inlet diameter increases. Thus, the V m did not change, and the E O was maintained at 86.0%.
Although the increase in the rotor inlet diameter is beneficial for the E O , the level of the annular width affects the degree of the E O rise. This is due to the interaction between the annular width and the rotor inlet diameter. We will explain this interaction in Section 3.2 using a response surface model.  ALH was only 6.9 mm. The rotor inhaled a large amount of gas. At this time, the pressure differential had no changes with the rotor inlet diameter increases. Thus, the V m did not change, and the E O was maintained at 86.0%.
Although the increase in the rotor inlet diameter is beneficial for the E O , the level of the annular width affects the degree of the E O rise. This is due to the interaction between the annular width and the rotor inlet diameter. We will explain this interaction in Section 3.2 using a response surface model. Figure 6 describes the five bottom vanes employed to investigate the effects of the vane's shape and number on the extraction efficiency. The result (Table 3) shows 4CV has the highest E O of 94.5%. Comparing with 6V and 8V, the E O of 4V increases by 2.3% and 5.2%, respectively. Wardle et al. [25] found that the bottom vane significantly affected the V m through computational fluid dynamics (CFD) simulation and experimental study. The function of the bottom vane is to break the liquid helical flow in the annular region and direct the liquid flow to the rotor inlet through vortices. 4V is the most commonly used type of vane in centrifugal contactor studies [26][27][28]. Compared with 4V, the bending direction of 4CV is the same as the liquid helical flow, which leads 4CV to enhance the fluid resistance flow to the rotor inlet. Therefore, the guiding role of 4CV was weakened, and the fluid was unfavorable to flow to the rotor inlet. 4CV had a higher V m and a longer mixing time than 4V. Thus, when 4CV was used instead of 4V, the E O rose.

Bottom Vane
4ACV had an opposite direction to the liquid spiral flow, which weakened the flow resistance at the bottom of the ACC and enhanced the bottom vane's guiding effect. 4ACV promoted the fluid flow to the rotor inlet. Thus, the V m was decreased, and the mixing time of the two-phase liquid was shortened. 4ACV lead to a reduction in E O .
The variation in E O with the bottom vane number is shown in Table 3. The addition of the bottom vane number enhanced vanes' guiding effect. The liquid was favorable to the flow toward the rotor inlet. Hence, 8V had a lower V m than 4V and 6V, and the mixing time of 8V was shorter. Consequently, the E O decreased with the increase in the bottom vane number.
We compared the E O of five different structures of bottom vanes in this subsection and interpreted the results by analyzing the flow of the fluid inside the ACC. For the straight vanes discussed in the current work, the E O is higher with a lower number of  Wardle et al. [25] found that the bottom vane significantly affected the V m through computational fluid dynamics (CFD) simulation and experimental study. The function of the bottom vane is to break the liquid helical flow in the annular region and direct the liquid flow to the rotor inlet through vortices. 4V is the most commonly used type of vane in centrifugal contactor studies [26][27][28]. Compared with 4V, the bending direction of 4CV is the same as the liquid helical flow, which leads 4CV to enhance the fluid resistance flow to the rotor inlet. Therefore, the guiding role of 4CV was weakened, and the fluid was unfavorable to flow to the rotor inlet. 4CV had a higher V m and a longer mixing time than 4V. Thus, when 4CV was used instead of 4V, the E O rose.
4ACV had an opposite direction to the liquid spiral flow, which weakened the flow resistance at the bottom of the ACC and enhanced the bottom vane's guiding effect. 4ACV promoted the fluid flow to the rotor inlet. Thus, the V m was decreased, and the mixing time of the two-phase liquid was shortened. 4ACV lead to a reduction in E O .
The variation in E O with the bottom vane number is shown in Table 3. The addition of the bottom vane number enhanced vanes' guiding effect. The liquid was favorable to the flow toward the rotor inlet. Hence, 8V had a lower V m than 4V and 6V, and the mixing time of 8V was shorter. Consequently, the E O decreased with the increase in the bottom vane number.
We compared the E O of five different structures of bottom vanes in this subsection and interpreted the results by analyzing the flow of the fluid inside the ACC. For the straight vanes discussed in the current work, the E O is higher with a lower number of bottom vanes. We find that 4CV could raise the E O better than 4V, while 4ACV, which is opposite to the flow direction, was not conducive to the E O . This conclusion could explain why the results

Response Surface Model Explanation of Each Factor Affecting the E O
As shown in Equation (5), the RSM (i.e., response surface methodology) represented the relationship between the geometrical parameters and E O by an MQF (i.e., multivariate quadratic function). This method reflects the direction and magnitude of geometrical parameters' effect on E O in the form of regression coefficients from statistical principles. Thus, MQF provides a new way of thinking by studying the change rate of geometrical parameters relative to E O and then giving a statistical interpretation for the experimental results. This idea is also one of the advantages of applying the RSM to assist the experimental design in this study.
The first derivatives of the E O in Equation (5) concerning three geometrical parameters were obtained separately. The first derivatives represented the slope of the curves. As shown in Equations (11)-(13), we used the second derivative of its first derivative with respect to another geometric parameter to indicate how another parameter affected the slopes between curves. The closer the second derivative was to zero, the more independent the effect of the parameter on the E O .
None of the second derivatives in the equations above are zero, which shows that the effect of each geometric parameter on the E O is not independent. There is an interaction between the parameters. Single-factor experiments cannot analyze the interactions between parameters, while fitting the interactions into the correlation equation is one of the advantages of RSM [29].
∂E O /∂D in = 1.04D in − 1.46d − 0.0554H c − 2.93 ∂ 2 E O /(∂D in ∂d) = −1.46 (13) Figure 4 displays that the E O varies with each factor. The slope change between the curves is not significant in Figure 4a,b, which correspond to the second derivatives close to zero in Equations (11) and (12). Therefore, the interactions of d-H c and H c -D in are weak. However, the second derivative of E O on D in -d is −1.48, indicating the variation in d interferes with the effect of D in on E O significantly. This result consists of the mass transfer theory analysis in Section 3.1.3.
We used the regression coefficients of the MQF to explain the interactions between geometrical parameters, demonstrating the applicability of the RSM in the study of the multiparametric effects of E O .

Conclusions
The present work studied the E O of an ACC with different geometrical parameters. We dissected the reasons for the effects of these parameters based on mass transfer theory and analyzed the interactions between the parameters according to the MQF.
We found that the E O is proportional to the clearance height and rotor inlet diameter. However, the effect of the rotor inlet diameter on E O is limited by the annular width. The E O is inversely proportional to the annular width and bottom vane's number. Moreover, the MQF predicted that the E O increases with the annular width when the annular width is more than 5.16 mm. The E O of 4CV is higher than 4V because 4CV has the same bending direction as the fluid flow. In contrast, the bending direction of 4ACV is opposite to the fluid flow, and the E O is lower than 4V. Analysis of the MQF indicates that each parameter's effect on the E O is not independent. The interaction of D in -d is more substantial, and we must consider it. This suggests that we should consider the interactions between the parameters comprehensively when conducting the experimental design or equipment design to determine a reasonable range of parameter values.
Due to the limitation of the existing 3D printing equipment in the laboratory, we did not investigate a wider range of geometrical parameters. However, the results in this work still show that the ACC's geometrical parameters influence the E O significantly, and there are interactions between parameters. These results could provide a reference for the design of ACCs and the optimization of the extraction process.