The Effect of Geometrical, Operational, Mixing Methods, and Rheological Parameters on Discharge Coefficients of Internal ‐ Mixing Twin ‐ Fluid Atomizers

: Accurate prediction of the discharge coefficient (C D ) for internal ‐ mixing twin ‐ fluid (IMTF) atomizers is challenging, the effect of control factors remains inadequately understood, and comparative data on the C D of IMTF atomizers are unavailable. This work presents an experimental study on C D for different IMTF atomizers with a wide range of factors, including the gas ‐ to ‐ liquid ratio (GLR), the inlet ‐ overpressure ratio ( ∆ p mix /p amb ), the orifice length ‐ to ‐ diameter ratio (L o /d o ), and the liquid viscosity ( μ L ). Five atomizers with different internal ‐ mixing principles were probed on a cold test rig, including the frequently studied outside ‐ in ‐ gas (OIG) and inside ‐ out ‐ gas (IOG) effervescent types, the recently ‐ introduced outside ‐ in ‐ liquid (OIL) and air ‐ core ‐ liquid ‐ ring (ACLR) atomizers, and our new design named the swirling ‐ air ‐ core ‐ liquid ‐ ring (SACLR) atomizer. The results demonstrate that C D is governed mainly by GLR, and reduces if GLR, L o /d o , or μ L is increased. An increase in ∆ p mix /p amb causes a C D reduction up to ∆ p mix /p amb = 0.98, and C D increases for a higher ∆ p mix /p amb . Surprisingly, differences in C D amid examined atomizers were found negligible, although the flow visualization inside the orifice showed a significantly different flow character for each one of the atomizers. Finally, a general C D correlation fitting with an R 2 ≥ 0.99 for all the tested nozzles was proposed. The results amend the present knowledge, allow design optimization, and provide flow rate prediction for a variety of IMTF atomizers.

TF atomizers use different approaches to produce two-phase mixtures for the intensification of liquid break-up. The TF atomizers in the present study significantly differ, for example, from the flash-boiling atomizers, even though both these categories work with a two-phase flow at the exit orifice. Flash-boiling atomizers use a single supply line with the bulk liquid phase only, and the twophase flow inside the atomizer is a product of a phase-change process (i.e., an intense gas phase nucleation or flashing boiling process). This effect can be achieved by exposing the bulk liquid to a sudden and sharp pressure drop (during an isothermal expansion process) or by isobaric heating. As a result, the liquid changes from a sub-cooled state to vapor or a superheated liquid state [19]. The typical flashing boiling atomizer contains two sharp-edged orifices (0.1 and 5 mm in diameter, respectively) with an expansion chamber that separates them. The flashing of the bulk liquid, which is usually pre-heated, takes place through the first small orifice. Due to the specific design, liquid physical properties, and the working conditions, a small portion from the liquid vaporizes to create a two-phase flow inside this atomizer. On the other hand, it is difficult for the phase change phenomena inside the internal-mixing twin-fluid (IMTF) atomizers to occur for reasons that are discussed in Section 2.2. In the TF atomizers used here, both the gas and liquid phases are separately introduced into the atomizer and are mixed into a two-phase fluid. The gas-liquid mixing can be carried out internally or externally. Internal-mixing twin-fluid (IMTF) nozzles feature lower atomizing gas consumption and higher energy transfer compared to their external-mixing twin-fluid (EMTF) counterparts (e.g., airblast atomizers [20]) [5,[21][22][23], which reduce operating costs.
The gas and liquid enter the IMTF nozzle through separate ports and interact inside a mixing chamber forming a two-phase mixture with a character varying depending on the mixing conditions (mixing chamber geometry, fluid pressures, and rheological properties). The mixture discharges from the exit orifice forming a spray. The discharge coefficient (CD) measures the extent to which the liquid flow uses the available exit orifice area [22].
The purpose of an atomizer is to deliver a specified amount of a liquid into a target place (defined with spray cone angle, droplet concentration, and liquid mass flux) and disperse it with the required quality (the spray characteristics such as droplet size and velocity, and droplet size distribution, are considered here). The amount of liquid discharged per time is represented with the flow rate or more conveniently with the CD of the atomizer. The importance of each factor category depends on the application. Precise liquid dosing and metering in a whole range of operating conditions are crucial in power generation units, such as combustion engines, turbines, burners, and furnaces; the throughput of the liquid governs mixing ratios of fuel with air in reacting sprays and directly affects combustion efficiency as well as pollutant formation. The amount of discharged liquid plays an essential role in process engineering and medical applications, particularly in mixing and dosing devices, chemical reactors, or columns for pollutant removal and capture. It is also considered in agriculture, food industry, or spray coating devices.
The atomized liquid is delivered to the target application in a specified amount within a range of operating conditions. Therefore, the discharge coefficient (CD) knowledge allows for predicting the actual liquid flow rate for given operating conditions; thus, playing an important role in the nozzle design process [22,[24][25][26]. It is advantageous to have a prediction tool for discharge characteristics of the atomizer before its realization. Iterative or trial-error techniques of fabrication and testing of flow characteristics in the whole range of operating conditions are costly and time-consuming. Such a calculation procedure giving correct data simplifies the atomizer design, and the final experimental verification at only several operation points is required. The discharge coefficient is also a suitable parameter for the validation of numerical simulations [27] that are challenging due to the two-phase flow nature as pointed out below. The CD is defined as a ratio of the actual liquid mass flow rate to the theoretical flow rate as can be described by Equation (1) [28,29]: where m . , A , ρ , and ∆p are the actual liquid mass flow rate, the exit orifice area, the liquid density, and the inlet-overpressure, respectively.
The CD determination of the IMTF nozzles is not such an easy task because both phases flow through the same exit orifice area [26]. Additionally, the gas-liquid flow is complex, random, and challenging in the study since this two-phase flow combines the characteristics of a deformable interface and the compressibility of one of the phases. A notable effort was paid to study the factors that mark the CD of the IMTF nozzles in order to obtain the CD predictive correlations found in the literature. Previous works are summarized in Table 1. All previous investigations are consistent in the effect of gas-to-liquid ratio (GLR) on CD (i.e., an increase in GLR leads to a non-linear decrease in the CD, while the decrement rate is rapid at low GLRs). However, there is some contradiction about the effect of liquid properties, the inlet-overpressure ratio, and the exit orifice structure on the CD. As a result, a comparison between the previously proposed models with the same conditions and geometries showed a notable difference in predictions of about 25%. Ochowiak [24] concluded that the full model for the CD calculation is unavailable yet, since there are complex relations between dependent and independent factors that affect the nozzle discharge.
Applications of IMTF nozzles cover a wide range of operating conditions and liquid rheological properties. The literature review (i.e., Table 1) suggests that the investigated factor ranges are in many cases very limited and narrow. Some factors were not included in the investigations at all, and a general and precise CD model has not yet been presented. Each of the studies was carried out on a single IMTF nozzle type, and only two old types (outside-in-gas (OIG) and inside-out-gas (IOG)) were examined. Therefore, each of the presented correlations is valid only for the type used, and no discharge characteristics are available for the newly introduced and promising IMTF nozzles. Additionally, the different IMTF nozzles possess a different way and intensity of mixing. This results in forming different internal flow regimes and characteristics [30,31]. Models based on the physical analysis of the two-phase discharge (i.e., [28]) suggest that, at the same operating conditions, the internal flow form affects the discharge characteristics; homogeneous gas-liquid flow (e.g., the homogenous bubbly flow) results in a CD decrease, and it is vice versa for separated flow (e.g., slug flow). Therefore, the internal-mixing method is expected to affect CD values, and the literature shows that this factor has not yet been investigated.
Consequently, in the present work, we tried to address these gaps by conducting a detailed experimental study on the discharge characteristics (i.e., CD) of different types of IMTF nozzles. Several factors were considered, including operation conditions, liquid types, and exit orifice geometry. The covered ranges of these factors were widened in comparison with the available literature. To study the effect of the internal-mixing method factor, for the first attempt, five types of IMTF nozzles with widely different internal-mixing methods were tested under the same conditions, liquid types, and orifice geometry. Two of the selected nozzles (outside-in-gas (OIG ) and inside-outgas (IOG)) are well-known and frequently used [1,5,10,13,24,28,[30][31][32][33][34][35][36]. The following two designs (outside-in-liquid (OIL) and air-core-liquid-ring (ACLR)) were introduced recently; they are promising due to a better atomization performance [10,30,31,34], [37,38]. The last (swirling-air-coreliquid-ring (SACLR)) nozzle represents our new design, inspired by the nozzles described by García et al. [7] Stähle and colleagues [34]. Finally, we proposed a precise model for CD prediction as a function of all the investigated factors.
-c is a constant, depends on the liquid properties, and is equal to 0.385 for water. CD decreases with increasing gasto-liquid ratio (GLR), while its value increases when the inlet-overpressure increases. The effect of viscosity is negligible. [39] -For horizontal flow. This model is valid for 2000 < Re < 20,000, with a coefficient of determination (R 2 ) of 0.982. The inlet-overpressure has a small effect on CD; increasing the gas pressure leads to decrease the CD. [28] 0.04 0.02 The terms that reflect the effect of physical properties were adopted from Chin and Lefebvre [29]. G can be obtained by combining two models (the homogeneous flow model (HFM) and the separated flow model (SFM)) of two-phase discharge, which can be obtained via the online multiphase flow calculator [40]. CD increases when liquid viscosity, surface tension, or GLR increases, while it is vice versa for the orifice length-to-diameter ratio.
[ 35] 0.75 The atomizer construction shows a strong effect on CD, in which CD increases when the orifice diameter do decreases. [24] 0.0475 -CD decreases with increasing liquid viscosity or GLR. In addition, CD is affected by the flow character and the atomization mechanism.

Materials and Methods
Experiments for measurement of CD and other relevant data were conducted by a cold atmospheric test rig at room temperature of 20 °C. A schematic drawing for the test rig is shown in Figure 1. The air was used as the gas phase, while water and aqueous glycerol solutions simulated the liquid phase. A wide range of operation factors was covered during the nozzle testing (i.e., GLR ranged between 1% and 65% and ∆pmix/pamb varied between 0.29 and 5.9).
An illustration of the operating conditions used for all the nozzles testing such as mass flow rate of the gas and the liquid, and the corresponding GLRs with their uncertainty values are documented in Table S1. Note that the ∆pmix/pamb was kept constant, while the GLR varied. Increasing the liquid flow rate resulted in a gradual decrease in the gas flow rate. The μL ranged between 0.001 and 0.2132 Pa•s (see Table 2), and the Lo/do ratio of the exit orifice varied between 0.5 and 4.2, while the convergence angle upstream the orifice (2β) was kept constant at 90° (see Table 3). The rheological properties of the liquids were found as follows. The density was calculated by dividing the known liquid volume by its weight, while the viscosity was measured by a controlled stress-rotary rheometer (DHR-2, TA Instruments, New Castle, DE, USA). Viscosity measurements were conducted at a constant shear rate of 10 s −1 .

Fluid Supply System
The liquid was driven by a commercial plunger pump (AR-RC 14.16 N, Modena, Italy) through a cylindrical accumulator used for damping the flow fluctuations. Sets of devices such as a needle valve for volume flow adjustment and a rotameter (DK800, KROHNE, Duisburg, Germany) were used. Rotameters with measuring ranges of 0.25-2.5 L/h, 1.2-12 L/h, 4-40 L/h, and 16-160 L/h were used for measuring the liquid volumetric flow rate with an accuracy of ±2.5% of the measuring value.
The pressure gauge was fixed to measure the liquid pressure, as shown in Figure 1. The flow directed to the nozzle was controlled by a bypass valve. The pressurized air was delivered from a highpressure air cylinder, passed through a gate valve, a pressure regulator used to control the ∆pmix/pamb. A Bourdon pressure gauge (ZYIA, Zhejiang, China, accuracy ±1.6% of the full scale) was fixed very close to the nozzle in order to measure the gas pressure. A non-return valve was used for the prevention of the water backflow to the gas line. A thermal mass flow controller (YJ-700CF-Air-10-500 SCCM, and YJ-700CV-Air-0.5-25 SLPM, Nanning, KONXIN, China, accuracy ±1% of FS) measured the gas mass flow rate. The funnel collected the sprayed liquid. Two K-type temperature sensors were located after the pressure gauges and used for measuring the gas and liquid temperatures.
When conducting the experiments of high viscosity mixtures, the plunger pump was replaced with a gear one. Additionally, before these experiments, the liquid flow meters were calibrated for metering the mixtures with different viscosities. The uncertainties for liquid density, liquid flow rate, and gas flow rate were used in the error analysis calculations for both the GLR and CD according to Holman [44]. The application of the mass flow controller for gas precluded the issue of gas density change when varying the pressure that would affect the true flow rate value. It also allowed us to achieve a large number of flow rate measurements (see Table S1) required for accurate and precise regression analysis.
The inlet conditions such as the flow fluctuations of the supplied fluids, eventual presence of dissolved gas in the liquid or inlet turbulence intensity can affect the performance of the atomizers and harm their comparison. The inlet fluid fluctuations were prevented using several precautions. The liquid supply line was equipped with a cylindrical accumulator, which provided a damping of liquid fluctuations eventually caused by the plunger pump. Additionally, to avoid the pulsations in the gas supply line, a high-pressure air cylinder was used instead of a reciprocating compressor to supply the pressurized air to the atomizer. Moreover, to avoid the effect of the possible presence of dissolved gas in the liquid, a treated water was used instead of the tap water. The treated water also prevents the metal components of the test rig to get oxidized and rusted.

Tested Nozzles
All tested IMTF nozzles were designed to work nearly under the same operation conditions; they feature comparable dimensions, as displayed in Figure 2, but their internal-mixing method for the two phases differs.
The OIG nozzle ( Figure 2a) was investigated many times [1,5,10,13,28,[30][31][32][33][34]36]. The mixing method is usually as follows: The liquid flows inside the main perforated chamber (aerator), while the gas is injected into the flowing liquid stream from outside through the aerator holes, perpendicularly to the main axis, resulting in a formation of a two-phase flow inside the mixing chamber that is located in the space between the last row of injection holes and the top of the exit orifice. The aerator contains 20 injection holes, each 1 mm in diameter, arranged in five rows. These rows are separated by 10 mm and turned by 45° to the neighboring one. The internal flow pattern is mainly influenced by the GLR and varies among bubbly, slug, intermittent, and annular flow [13,36].
The OIL nozzle (Figure 2b) was first introduced by Mlkvik et al. [30] in 2015 and studied later [31,42]. Its design is similar to the OIG nozzle, but the mixing method is totally different; the gas flows through the main perforated chamber, while the liquid enters into the flowing gas stream from outside through the holes at the chamber wall. This mixing method is most similar to the Y-Jet nozzle mixing. The internal flow pattern for the OIL nozzle was theoretically predicted to be annular at any GLR [30,31,42], but this hypothesis has not yet been experimentally proven.
The IOG nozzle (Figure 2c) has been frequently investigated as well [5,24,28,35,45,46]. In this configuration, the liquid flows inside the main chamber, while the gas enters into the aerator tube and bubbles outward to mix with the liquid. The aerator, as in the previous designs, contains 20 injection holes, each 1 mm in diameter, arranged in five rows. These rows are separated by 10 mm and turned by 45° to the neighboring one. The aerator internal diameter is 4 mm, and its wall is 1 mm thick. The internal flow pattern is heavily influenced by the GLR and thus similar to the OIG configuration.
Correspondingly, the ACLR nozzle (Figure 2d) was designed and introduced by Stähle et al. [34] and studied later by other researchers [10][11][12]. In contrast to the OIG and IOG nozzles, the air in the ACLR nozzle is not injected into the liquid through multiple holes perpendicularly to the main axis, but through a single capillary positioned coaxially to the main axis and it exits very close to the discharge orifice. This arrangement produces a high-velocity air core surrounded by a liquid annulus in the mixing position, which provides an annular flow pattern inside the orifice at any GLR [34].
Finally, the new SACLR nozzle was designed to introduce an additional type of internal mixing which involves a swirling motion. Figure 2e shows that, in this case, the air enters through a tangential inlet at the top and swirls inside the conical conduct towards the exit hole at the base. This hole is coaxial to the main atomizer axis and ends very close to the top of the discharge orifice. The air exits this hole to enter into the mixing chamber in the form of a swirling air core, while the liquid is introduced parallel to this core, which forms a liquid ring around the swirling core directly upstream the exit orifice. Considering the phase change phenomena, it is difficult for the liquid inside the IMTF atomizers to be exposed to an intense nucleation at any region. Neither the internal design nor the operating conditions help to initiate such a process. As shown in the literature review, and in the present study, these atomizers usually work under very low inlet-overpressure (∆pmix/pamb <6) and the mixing of the two phases carries out at a constant pressure. This means that the liquid pressure outside the main perforated chamber (i.e., in the OIL design case) is slightly higher than that inside the chamber. In addition, the number and the diameter of the holes in the chamber are large (i.e., 1 mm, see Figure 2) in comparison with the first, small orifice of the flash-boiling atomizer (i.e., 0.1 mm in diameter). In addition, the liquid temperature is low and is usually equal to the room temperature. Consequently, the liquid flows through these holes with a negligible pressure drop, and therefore without any phase transition. After the mixing process, the two-phase (water-air) mixture flows through the convergent/exit orifice region (the convergence angle upstream the orifice is 90°). During this, the average mixture pressure gradually drops down until it equalizes with the atmospheric pressure. Furthermore, it was observed in other studies [47][48][49] that the nucleation phenomena considerably reduces and delays if a small amount of dissolved gas is present in the pure liquid. In other words, the gas bubble presence in the liquid inhibits the phase change process if compared to the pure liquid case. Therefore, the cold water-air mixture flowing through the convergent/exit orifice region in the IMTF atomizer is not likely to be exposed to the phase change process. It is expected that this phenomenon does not affect the discharge characteristics of the IMTF atomizers.

Internal Flow Visualization
For the flow visualization inside the exit orifice, the segment containing the contraction zone and exit orifice was machined (with do = 0.85 mm, and Lo/do 1.1) in a square cross-section acrylic (Plexiglas) block to reduce the optical distortion caused by the curved shapes. This part was backlit illuminated by a diffused 500 W light. A high-speed camera (Fastcam SA 1.1, Photron, Tokyo, Japan) with a long distance microscope (Questar-QM1) was used to capture the flow images with 2000 frames/s rate and 1024 × 1024 pixel image resolution.

Methodology of Data Analysis
The literature review revealed that CD is affected by working conditions [26], the exit office size and geometry [28,35], and the physical properties of the phases [24,28]. Here a general function applies: After performing the dimensional analysis, CD can be expressed by Equation (3): The form of CD correlation can be obtained by applying a non-linear regression analysis with the least-squares method to the experimental data. Note that the gas-to-liquid mass flow ratio is defined as shown by Equation (4): The effect of the investigated factors on CD was evaluated via an analysis of variance (ANOVA) test. ANOVA is a form of statistics hypothesis testing that determines the effect of one or more factors on an outcome variable by comparing the means of different groups. ANOVA is widely used in the analysis of experimental data [50]. The test result, which is usually presented by the probability (pvalue) that is calculated from the null hypothesis and the experimental data sample, is considered statistically significant when the p-value is less than a pre-specified threshold. A predictor that keeps a low p-value has a significant effect on the model accuracy as changing that predictor value causes a considerable change in the model response. On the other hand, high p-value predictors do not have a significant effect on the model response. For the identification of significant differences, the alpha value was set equal to 0.05.
For each case, 250 flow visualization images were captured and processed, and this number was found sufficient for the minimization of bias errors in the measurement. Thereafter, MATLAB software was used for determining the gas phase boundaries at the entrance of the exit orifice and subsequently determining the diameter of the gas core (see Figure 3). The gas core diameter was used to measure the time variation of the gas phase inside the exit orifice, which corresponds with the internal flow regime [51,52]. Additionally, this diameter can be used to indicate some parameters for the gas-liquid flow (i.e., the gas void fraction and the corresponding slip ratio), as represented by Equations (5)-(8).
αi is the time-varying void fraction, AG is the cross-sectional area occupied by the gas phase, Ao is the orifice cross-sectional area, CVαi is the coefficient of variation of αi, which reflects its fluctuation, σα is the standard deviation of αi, αavg is the time-averaged void fraction, and N is the image number. The CVαi increase indicates a wide distribution range of αi and its high fluctuations. S is the slip ratio (also called the velocity ratio), which represents the actual velocities ratio of both gas and liquid phases, respectively. UG and UL are phase superficial velocities of the gas and liquid, respectively.

Exit orifice walls (blue)
Gas phase boundaries (red) Time Figure 3. Processing procedure for determining the gas core diameter at the entrance of the exit orifice.

Analysis of the Experimental Errors
The uncertainties in the experimental results are caused by the errors in the primary measurements. Experimental error analysis was provided according to Holman [44]. Let R be a function given as R = f (X1, X2, X3… Xn), where X1, X2, X3 ,..., Xn are independent variables with uncertainties of WX1, WX2, WX3 ,..., WXn, respectively. Therefore, the uncertainty in the R value (WR) can be calculated as follows: and the relative error is The relative error of the dependent variables is determined as follows: Density: Mass flow rate of liquid and gas: Gas-to-liquid ratio: Discharge coefficient: The relative errors for the CD are displayed as error bars at each measuring point in the following result graphs. For instance, the uncertainties for mass flow and GLR are documented in Table S1.

Results and Discussion
In this section, the results of the OIL nozzle are presented, while the results of the other nozzles are presented in Figures S1-S4. The comparisons between the obtained results of the tested nozzles will be illustrated in Section 3.4.

The Effect of GLR and Orifice Structure on CD
The effect of both GLR and Lo/do on the CD of the OIL nozzle is shown in Figure 4. It is clear that CD features a sharp decrease when increasing the GLR to ≈30%, while a further GLR increase has little effect on CD (see Figure 4a where the absolute values for the CD are documented). This behavior results from the link between the GLR and the flow area occupied by the liquid inside the exit orifice [26]. The actual area available for the liquid flow decreases with any increase in the gas flow rate (i.e., with increasing GLR). It is apparent that the discharge coefficients of the twin-fluid atomizers are, due to the gas blockage, lower than that reported for plain-orifice atomizers and single-fluid nozzles [53], simplex [54], and even spill-return pressure-swirl atomizers [55]. It moreover drops down below 0.1 for GLRs increasing above 10%. Figure 4b documents that, with a constant GLR, CD decreases with an increasing Lo/do ratio. Here the friction force on the liquid flow plays a role, as it is proportional to the exit orifice length. Additionally, the sensitivity of CD to the Lo/do variation is high at low GLRs, while it gradually decreases with the GLR increase. The obtained results show that, when Lo/do increases from 0.5 to 4.2, there is a pronounced drop of CD by 35.6% and 18% at GLRs of 3.7% and 50.2%, respectively. The reason is that, at low GLRs, the liquid occupies the largest portion of the exit orifice and is exposed to the largest friction force. Therefore, a small increase in the exit orifice length results in significant resistance to the liquid flow, and a notable CD reduction. At high GLRs, the amount of liquid passing through the orifice is very small relative to the gas amount, and the friction force acting on this liquid is also small.
It was found that an exponential decay function CD = A × EXP (B/ (GLR + C)) fits the experimental data accurately as can be seen in Figure 4. The results of other nozzles (OIG, IOG, ACLR, and SACLR) show similar trends as the OIL nozzle (see Figure S1). Regarding the results of the IOG nozzle, it is consistent with Jedelsky and colleauges [28] and Chin and Lefebvre [29] and inconsistent with Ochowiak [45]. CD was found to increase to a maximum value as Lo/do increases to 2.04 and then to decrease with a further increase in Lo/do. Finally, results from the ANOVA show a statistically significant difference in CD with varying GLRs and Lo/do (i.e., a p-value equal to 6.94 × 10 −27 and 5.49 × 10 −20 , respectively), so the effect of both GLRs and Lo/do factors on CD cannot be neglected.

The Effect of the Inlet-Overpressure Ratio on CD
The effect of ∆pmix/pamb (pamb was constant in all the experiments, i.e., atmospheric pressure) on the CD of the OIL nozzle is shown in Figure 5. Figure 5a documents the notable effect of GLRs on CD at different ∆pmix/pamb, while the effect of the ∆pmix/pamb is shown in Figure 5b. Figure 5b shows that an increase in ∆pmix/pamb at constant GLRs results in the CD reduction to its minimum value, which occurred at ∆pmix/pamb = 0.98. Then CD increased for higher ∆pmix/pamb (i.e., > 0.98). Figure 6 illustrates the reasons for this non-monotonic phenomenon. Increasing the ∆pmix at a constant GLR causes a decrease in the gas volume flux at the exit orifice due to the high gas compressibility compared to liquids. As a result, the liquid mass flux (G) increases. G features a significant change at ∆pmix ≥ 0.1 MPa in comparison with its change at lower ∆pmix. According to the known relation CD  G/(∆pmix) 0.5 [39], at ∆pmix < 0.1 MPa, the effect of the pressure increase is predominant and overcomes the effect of the G increase, so CD decreases; it is vice versa for ∆pmix > 0.1 MPa (see Figures 5b and 6). Additionally, as aforementioned, Figure 5b shows that the increase in GLR suppresses the effect of ∆pmix/pamb on CD because of a decrease in the liquid amount. This phenomenon is of significant importance for IMTF atomizers; firstly, it is related to CD, which is important and an influential factor for atomizers. Secondly, this phenomenon occurs at regimes of low pressures and low GLRs (see Figure 5b), where the nozzles feature the highest atomization efficiency [56]. GLR  Additionally, we found that the other IMTF nozzles (OIG, IOG, ACLR, and SACLR) perform similarly with respect to the OIL nozzle (see Figure S2). It is worth mentioning that we observed this non-monotonic phenomenon in a wider range for the ∆pmix/pamb in this study than those in the literature (see Table 1). This finally explains and unifies the previous contradicting results by Chen and Lefebvre [26] and Ramamurthi et al. [39], who found that CD monotonically increased and decreased with the ∆pmix/pamb in narrower operation ranges, respectively. Lastly, results from ANOVA show that the influence of ∆pmix/pamb on CD is statistically significant and cannot be neglected (i.e., a p-value is equal to 1.88 × 10 −10 ).

The Effect of Liquid Viscosity on CD
The effect of viscosity on the CD of the OIL nozzle is shown in Figure 7. CD decreases as the μL/μG increases at constant GLRs, as shown in the figures. This effect can be interpreted via the relation between CD and the Reynolds number for the liquid flow (ReL) (Equation (15)): (15) where UL is the cross-sectional average velocity of the liquid phase. CD depends on ReL in a linear manner [24,35,57], as confirmed in Figure 8a. The value of ReL drops with rising liquid viscosity at constant GLRs (Figure 8b). By increasing the liquid viscosity, the internal friction of the liquid increases, which leads to a decrease in the discharge velocity of the liquid as well as the suppression of the flow turbulence [58,59], and CD consequently decreases. The same effect on ReL was found with the increase in GLR (Figure 8b). The sensitivity of CD to the variation in μL/μG is high at low GLRs and it gradually decreases with increasing GLR; the increase of μL/μG from 55 to 11,779 results in a reduction in CD by 53.7% at GLR = 3.7%, while it is 44% at GLR = 45.8%.
The results of other IMTF nozzles show similar tendencies as the OIL nozzle (see Figure S3). Our findings on the viscosity effect are in good agreement with the results reported by Ochowiak [24] for the IOG nozzles, while they are inconsistent with those reported by Chen and Lefebvre [26] for the OIG type. Once more, ANOVA results show that the p-value of the effect of μL/μG on the CD is equal to 4.60 × 10 −13 , so this effect cannot be neglected.

The Effect of the Mixing Method on CD
As mentioned, the nozzle type outlines the mixing process. The influence of the mixing on CD is shown in Figure 9 at different conditions. It is clear in Figure 9a-d that the differences in the CD values are very small and may even approach zero. This finding was confirmed by ANOVA results that showed insignificant differences in CD values among different nozzles (i.e., p-value is greater than 0.15). Overall, the effect of the mixing method factor is negligible. It is interesting that every nozzle type possesses its own internal flow pattern and characteristics, because the mixture formation and its character are strongly influenced by the mixing process [30,31], especially at low GLRs. The literature sources (see Table 4) show that the internal flow regime varies from one nozzle to another between bubbly, slug, intermittent, and annular. Shepard [60] found that the flow patterns can change during the residence time of the mixture inside the atomizer, so the regimes at the exit orifice can differ from those at the gas injection area. Our visualization results for the actual flow inside the exit orifice for the tested nozzles showed significant differences in the flow nature between nozzles.  Figure 10 shows the character of αi at the entrance of the exit orifice for the tested nozzles at the given condition (i.e., GLR = 1.2%, ∆pmix/pamb = 0.98, and μL/μG = 55). The variation in the character of αi between the tested atomizers is certainly related to the nature of the internal flow regime inside the mixing chamber. It is clear in Figure 10a-e that αi is not stable and features a different character for each nozzle, which means the internal flow regime inside each nozzle is quite different from the others. This also indicates that the internal flow of these nozzles is inherently unsteady, but with varying levels. Only the exit orifice for the OIL nozzle continuously possesses an annular flow, while the other nozzles possess different flow regimes.
The corresponding CVαi values for αi in Figure 10, while CD values and slip ratios for each nozzle are shown in Figure 11. Figure 11a shows that the CVαi features substantial differences between nozzles. In contrast, the differences in CD between nozzles were not significant (i.e., p-value is equal to 0.24). These results reveal that the CD of the IMTF nozzles is independent of the form of the internal flow regime.
It is worth mentioning that the CD prediction model based on the analytical physical models of two-phase discharge (i.e., homogeneous flow model (HFM) and separated flow model (SFM), which correspond to the two-phase flow discharge without slip between phases and the discharge with maximum slip, respectively) predict that the CD for the SFM should be higher than that for the HFM by three times [62]. The present experimental data show that there is no difference ( 0) in CD values for the nozzles that feature the highest and lowest slip ratio (i.e., SACLR and OIL, respectively) (see Figure 11a,b). Therefore, CD prediction models should not rely solely on the analytical models of twophase discharge. Additionally, a comparison between the inlet turbulence intensity (Tu1) values of the continuous single-phase flow entering the mixing chamber of all the atomizers is displayed in Figure 12. The gas represents continuous single phase for the OIL atomizer, while it is the liquid in the other configurations (OIG, IOG, ACLR, and SACLR), as shown in Figure 2. Tu1 can be calculated as follows [63,64]: Figure 12 compares the differences in Tu1 among the tested nozzles. It suggests that the mixing process of the two phases, and consequently the nature of the produced two-phase mixture were possibly affected by the inlet conditions. The maximum difference in the Tu1 is about 7.4% and applies between the OIL atomizer (i.e., the highest Tu1) and the SACLR atomizer (i.e., the lowest Tu1). While the difference in the CD between the same atomizers (under the same operating condition) is only about 0.4%. It is implied that the CD values were negligibly affected by the differences in Tu1.
Finally, Table 4 indicates that the internal flow inside the ACLR was annular flow, while our visualization showed a different flow character. This may be due to the change in the ratio between exit orifice diameter and gas capillary diameter; it was about 0.2 and 1 in the present and previous study, respectively. We set these dimensions so as to test all the nozzles under the same operating conditions and exit orifice geometry. Additionally, it is interesting that the flow nature inside the exit orifice of the OIL nozzle was more stable than the other nozzles. This behavior gives an advantage to this type, given that this stability will be reflected in the quality of the spray. Thus, the OIL nozzle can be useful in applications that require a more stable spray with high quality (combustion, surface coating, and powder generation).

The Proposed CD Model and a Comparison with Published Models
Since the effect of the mixing method is negligible, a general CD prediction model can be developed for different IMTF nozzles. The non-linear least-squares regression technique was applied to the experimental data of CD for each nozzle, and it was found that the CD can be correlated as follows: This general equation is valid for all the IMTF nozzles (i.e., internal-mixing types). Table 5 summarizes the values of the coefficient of determination (R 2 ) and the sum of the squared residuals (SS). A comparison of the experimentally found CD with the CD predicted by Equation (17) for the OIL nozzle is shown in Figure 13, and the others (OIG, IOG, ACLR, and SACLR) are documented in Figure S4. The figures and the table confirm that the proposed CD model fits the experimental data for each mixing type very closely with R 2 ≥ 0.99.
A comparison between the present model with the other models published for the OIG and IOG nozzles in Figure 14 shows that the CD predicted by the published models for both nozzles features a significant dispersion along the expected trend. This could be because of the neglection of some important factors in these equations, or because the ranges of the investigated factors are too narrow and inconsistent with present experiments (see Table 1), and/or because the effect of some factors is inconsistent with the findings of the present work. Note that the SFM (a = 0, b = 1) with polytropic expansion was selected during the calculations of Jedelsky's model in Figure 14a.  Figure 13. Comparison of the experimental and predicted CD from the present model in Equation (17) for the OIL nozzle.

Conclusions
In this paper, a detailed experimental study was carried out to investigate and analyze the effect of many factors upon the discharge coefficient of different types of internal-mixing twin-fluid atomizers. These factors mainly are the gas-to-liquid mass flow ratio, the inlet-overpressure ratio, the orifice length-to-diameter ratio, the internal mixing method of the two phases, and the liquid viscosity. Five different types of IMTF nozzles-OIL, OIG, IOG, ACLR, and SACLR-were tested under the same working conditions, liquid types, and exit orifice geometry. Additionally, the internal flow at the entrance of the exit orifice was visualized using a high-speed camera at a certain operating condition for each nozzle. Based on the results, the main conclusions are as follows: 1. The GLR, Lo/do, μL/μG, and ∆pmix/pamb have a significant effect on the CD of IMTF nozzles, while the effect of the mixing-method factor is negligible. 2. CD decreases with the increase of GLR, Lo/do, or μL/μG, and the same trend was observed with the increase of ∆pmix/pamb up to a value of 0.98, but CD increases with any further increase in ∆pmix/pamb. 3. The sensitivity of CD to changes in Lo/do, ∆pmix/pamb, or μL/μG is high at a low GLR but declines at a high GLR. 4. Flow visualization inside the exit orifice showed a diverse internal flow nature for each nozzle, which indicates that the CD of the IMTF nozzles is independent of the form of the internal flow regime at any given condition. 5. The internal flow inside the exit orifice of the tested IMTF nozzles is inherently unsteady, but with varying levels; the OIL nozzle features a more stable internal flow behavior than the other nozzles. 6. Finally, a general well-fitting CD prediction model (Equation (17)) with an R 2 ≥ 0.99 was proposed, which is valid for any type of IMTF nozzle and with wider validity ranges for the included factors.

Conflicts of Interest:
The authors declare no conflicts of interest.