E ﬀ ects of Nanoparticle Additives on Spray Characteristics of Liquid Jets in Gaseous Crossﬂow

: Nanoﬂuids are attracting attention as future energy carriers owing to their high performance for improving combustion and heat transfer. In this study, the macroscopic characteristics of nanoﬂuid jets in a subsonic gaseous crossﬂow were investigated by focusing on the inﬂuence of nanoparticle additives on the breakup process. Based on a distribution map of the image grayscale standard deviation, we propose an improved method to process transverse injection shadowgraphs. A simpliﬁed model of the transition mechanism from column breakup to surface breakup at a small Weber number was established. The e ﬀ ects of nanoparticles on the jet trajectory and column fracture position were analyzed according to the deviations from the pure liquid. To interpret the e ﬀ ects of the nanoparticles, a new nondimensional parameter was introduced into the empirical correlation of the column fracture position. The results indicated that at low concentrations of nanoparticles, the surface tension of the nanoﬂuids increased slightly, while the viscosity increased signiﬁcantly (by up to 23%). These changes in the physical properties had little e ﬀ ect on the breakup regimes or jet trajectory. Moreover, the nanoparticles promoted cavitation inside the liquid column, resulting in an additional primary breakup mode for the nanoﬂuids. Consequently, the length of the column fracture was reduced by up to 20% compared with that of the basic ﬂuid.


Introduction
Liquid injection into a subsonic gaseous crossflow has been widely applied to energy and power systems, such as subsonic-combustion ramjets, lean premixed pre-evaporation low-emission burners, and diesel engines [1,2]. The spray of fuels significantly affects the combustion efficiency, flame stability, and emission characteristics; therefore, studies have been performed on the spray characteristics for decades [3,4]. Regarding the physical complexities of transverse injection, researchers have provided various useful correlations based on experiments to quantify the primary breakup regime, penetration characteristics, and atomization performance [5]. Several reports have indicated that the Weber number of the crossflow (We g ) plays a key role in defining nonturbulent primary breakup regimes [6,7]. The most widely used diagram for partitioning the different breakup modes of transverse injection is the We g -q map, which was proposed by Wu et al. in 1997 [1]. In their report, the primary breakup regimes were classified into four modes: enhanced capillary, bag, multi-mode, and shear breakup. Wu et al. described an interesting aspect of the breakup process-the surface breakup-and provided an empirical correlation of the transition line between surface breakup and onefold column breakup. Subsequently, many different correlations based on experimental data have been reported [8,9]. Broumand et al. summarized these correlations and proposed two boundary particles was reduced. Li et al. recently performed an experiment on the electrostatic atomization of an Al-ethanol nanofuel [35]. They found that as the concentration of nano-Al increased, the density, surface tension, and viscosity increased. Additionally, the nanofuel exhibited better atomization performance when an electrical field was applied. Similar to the addition of nanoparticles, Bolszo et al. added micron-sized water droplets to DF2 fuel and studied the atomization characteristics of the emulsion in air crossflow [36]. According to a dimensional analysis of current physical problems, it was verified that the liquid jet trajectory correlation reported by Wu et al. can successfully predict the penetration of an emulsion jet. An additional primary breakup mode for emulsion spray, which is called the interfacial tension breakup, was discovered to cause the column fracture in advance and influence the droplet size distribution. For the analysis and data fitting of the droplet size, a new nondimensional quantity was proposed to account for the effects of body forces and interfacial tension. This research method provides a good reference for the transverse injection of nanofluids.
It has been reported that the addition of nanoparticles to basic fluid can induce changes relevant to the atomization process, including changes in the physical properties and spray characteristics, such as the viscosity, surface tension, cavitation, and velocity distribution. Therefore, the detailed effects of these changes on fuel atomization deserve further study. In light of the meaningful physical problems and current research gaps, the effects of alumina nanoparticle dispersions on the transverse injection of water were examined in this study. A series of experiments were performed at room temperature and the atmospheric pressure. The macroscopic spray characteristics and primary breakup regimes were investigated, mainly via high-speed shadowgraphs and statistical image processing of the grayscale standard deviation (SD).

Preparation and Properties of Nanofluids
Alumina nanoparticles of the γ crystal form (γ-Al 2 O 3 ) with high catalytic activity were selected as the additives. They had an average diameter of 20 nm and a true density of 3.97 g/cm 3 (Shanghai Aladdin Biochemical Technology, Inc., Shanghai, China). As repeated spraying was involved, deionized (DI) water (Shanghai Macklin Biochemical Co., Ltd., Shanghai, China), instead of fuel, was selected as the injection liquid. Sodium dodecylbenzene sulfonate (SDBS) was used as a surfactant to improve the suspension stability (Shanghai Aladdin Biochemical Technology Inc.). Al 2 O 3 -H 2 O nanofluids were prepared via a well-known two-step method. First, to simplify the comparison of fluids with different nanoparticle concentrations, SDBS at a fixed proportion of 0.2 wt.% was added to DI water (called the "basic fluid"), according to a dispersion stability test. Subsequently, the mixture was mechanically stirred for 20 min at a suitable rotation speed (approximately 600 rpm, to avoid a large amount of foam). Finally, a certain proportion of alumina nanoparticles were added in, and ultrasonic dispersion was conducted for 40 min (Scientz-2400F ultrasonic disperser, Ningbo Scientz Biotechnology Co., Ltd., Ningbo, China). The optimized dispersing conditions were an ultrasonic frequency of 20 kHz, power of 960 W, and a water bath temperature of 25 • C. In a turbidity test for dispersion stability evaluation, the prepared nanofluids exhibited no stratification or sedimentation over 4 h (as shown in Figure 1). With the addition of SDBS as an anionic surfactant, the alumina nanoparticles were well dispersed in water in a stable structure of the electron double layer.
The density of the nanofluids was calculated according to the mixing rule, i.e., ρ n f = ϕρ p + (1 − ϕ)ρ L , while the effect of the surfactant was disregarded. The viscosity was measured using a rotational viscometer (LVDV-1T, Shanghai FangRui Insurument Co., Ltd., Shanghai, China). The surface tension was measured via the pendant-drop method (Dataphysics OCA25, DataPhysics Instruments GmbH, Filderstadt, Germany). A single test was performed, with 3-5 repeated measurements for data reliability. The equation for converting the mass fraction of nanoparticles into the volume fraction is ϕ = ρ L / ρ p /w + ρ L − ρ p . The liquid properties are presented in Table 1 and Figure 2. The test results indicated that both the viscosity and surface tension of the nanofluids increased with an increase in the concentration of nanoparticles. With the addition of 0.51 vol.% nanoparticles, the surface tension of the Al2O3-H2O nanofluids increased with a maximum ratio of 2.6% relative to the basic fluid. The surface tension depends on the physical and chemical properties of dispersions, nanoparticles, and surfactants, as well as their mixed states. The increase in the surface tension in current experiment was attributed to the higher surface free energy of the nanofluids relative to the  The liquid properties are presented in Table 1 and Figure 2.  The liquid properties are presented in Table 1 and Figure 2. The test results indicated that both the viscosity and surface tension of the nanofluids increased with an increase in the concentration of nanoparticles. With the addition of 0.51 vol.% nanoparticles, the surface tension of the Al2O3-H2O nanofluids increased with a maximum ratio of 2.6% relative to the basic fluid. The surface tension depends on the physical and chemical properties of dispersions, nanoparticles, and surfactants, as well as their mixed states. The increase in the surface tension in current experiment was attributed to the higher surface free energy of the nanofluids relative to the  The test results indicated that both the viscosity and surface tension of the nanofluids increased with an increase in the concentration of nanoparticles. With the addition of 0.51 vol.% nanoparticles, the surface tension of the Al 2 O 3 -H 2 O nanofluids increased with a maximum ratio of 2.6% relative to the basic fluid.
The surface tension depends on the physical and chemical properties of dispersions, nanoparticles, and surfactants, as well as their mixed states. The increase in the surface tension in current experiment Energies 2020, 13, 1574 5 of 23 was attributed to the higher surface free energy of the nanofluids relative to the basic fluid. The reason for this at the microscale is that alumina nanoparticles tend to accumulate at the water-air interface, which increases the van der Waals attraction between them [37]. A small amount of surfactant significantly affected the surface tension. The surface tension of the basic fluid was 55% lower than DI water, which may significantly affect the liquid breakup.
The increase in viscosity is attributed to additional Brownian motion, stronger van der Waals forces, etc. As previously reported [38], even nanofluids with a small amount of nanoparticles are Newtonian fluids. In the present experiment, the viscosity of the Al 2 O 3 -H 2 O nanofluids increased significantly, by up to 23%. For a sparse suspension (ϕ < 2%), Einstein proposed a classical viscosity model [39]: µ n f /µ b f = 1 + 2.5ϕ. For nanoparticles, the actual volume fraction increases owing to the formation of a solvation shell (as shown in Figure 1). Thus, the calculation model needs to be partly modified. Combined with the experimental data, the new correlation is: where, d p0 represents the mean diameter of monodispersed nanoparticles (20 nm), and d m0 represents the equivalent diameter of water molecules (0.39 nm), which is calculated as follows: where M represents the molar mass of water molecules, N represents Avogadro' s constant, and ρ 0 represents the water density at 293 K. Compared with the classical equation, the equivalent volume fraction is changed to 14.9ϕ, and the thickness of the solvation layer is calculated as 2∆R/d p0 = 1.46. This parameter is used later in the breakup analysis.

Spray Experimental Facility
As shown in Figure 3, the experimental system mainly comprised three parts: a liquid spray unit, an airflow control unit, and a shadow imaging unit. The spray unit provided continuous liquid injection through a multifunctional pressure tank (5 L/1 MPa). The flow rate of the nanofluids was adjusted by a needle valve and measured using an electromagnetic flowmeter (RONK DN4, Chongqing Sichuan Instrument Co., Ltd., Chongqing, China). A high-precision pressure sensor was used to monitor the injection pressure. To inhibit extra interferences from cavitation and turbulence, a conical nozzle with a length/diameter ratio of 4 was used as described by Wu et al. [1] (Figure 3). With reference to the conclusion of Faeth et al. [40], a sufficiently long straight pipe section (length-to-diameter ratio of 40) was arranged upstream of the nozzle inlet. Therefore, the flow in the tube was fully developed, and the liquid at the nozzle inlet was uniform in all directions.
The airflow control unit uses a large-capacity gas tank (2000 L/4 MPa) for air supply. After filtration and drying using a cold dryer, the tank was filled with air by an uprated compressor. The airflow was precisely controlled using a decompressor group, buffer tank, and throttle valves. A multi-hole balance flowmeter (RONK DN65) was used to monitor the airflow rate, and a hot-wire anemometer (MODEL 8465, TSI Incorporated, Minnesota, USA) was installed in the test section for calibration. To obtain a uniform air crossflow, three rectifying nets were arranged at different connections of the variable cross-section upstream. The cross section of the test zone was rectangular, with a size of 100 × 75 mm 2 , and there was a transparent window made of quartz glass with a size of 100 × 200 mm 2 . A high-speed data collector was used to obtain the real-time flow status, including the pressure, temperature, flow rate of the airflow and liquid jetflow. To verify the stability of the air crossflow, a TSI 8465 sensor was used to measure the air velocities at two different positions under an air supply Energies 2020, 13, 1574 6 of 23 pressure of 0.6 MPa. One position was the center of the air-duct cross section, directly above the nozzle exit (Z = 50), and the other was 75 mm above the nozzle exit (Z = 75).  As shown in Figure 4, the average airflow velocities at these two locations were highly close (36.224 and 36.218 m/s), and they changed little within 5 s, indicating that the airflow had good spatial uniformity and time stability. The air velocity in the central area of the duct was relatively stable (with an SD of 0.03 m/s), while the velocity near the wall exhibited large fluctuations (with an SD of 0.26 m/s), which was mainly caused by disturbances due to the boundary layer. Additionally, the average velocity and its SD were calculated via three repetitive tests at different air supply pressures. The maximum deviation was < 4%, indicating the good reliability of the air supply system.  Experiments were performed at room temperature and the atmospheric pressure. The test conditions are presented in Table 2.   As shown in Figure 4, the average airflow velocities at these two locations were highly close (36.224 and 36.218 m/s), and they changed little within 5 s, indicating that the airflow had good spatial uniformity and time stability. The air velocity in the central area of the duct was relatively stable (with an SD of 0.03 m/s), while the velocity near the wall exhibited large fluctuations (with an SD of 0.26 m/s), which was mainly caused by disturbances due to the boundary layer. Additionally, the average velocity and its SD were calculated via three repetitive tests at different air supply pressures. The maximum deviation was <4%, indicating the good reliability of the air supply system.  As shown in Figure 4, the average airflow velocities at these two locations were highly close (36.224 and 36.218 m/s), and they changed little within 5 s, indicating that the airflow had good spatial uniformity and time stability. The air velocity in the central area of the duct was relatively stable (with an SD of 0.03 m/s), while the velocity near the wall exhibited large fluctuations (with an SD of 0.26 m/s), which was mainly caused by disturbances due to the boundary layer. Additionally, the average velocity and its SD were calculated via three repetitive tests at different air supply pressures. The maximum deviation was < 4%, indicating the good reliability of the air supply system.  Experiments were performed at room temperature and the atmospheric pressure. The test conditions are presented in Table 2.   Experiments were performed at room temperature and the atmospheric pressure. The test conditions are presented in Table 2. The shadow imaging unit consisted of a high-speed complementary metal-oxide semiconductor camera (Phantom V2640, Vision Research Inc., New Jersey, USA), a light-emitting diode (LED) (1000 W), and a homogenizing plate. A Nikon macro-lens (100 mm, f2.8D, Nikon Corporation, Tokyo, Japan) was installed on the camera. Spray shadowgraphs were obtained at a frame rate of 4000 fps and an exposure time of 5 µs, and the size of each image was 2048 × 1952. After calibration with a Vernier caliper, the size of a single pixel was 91.7 µm. To obtain well-focused images of the spray plume distributed across the entire cross section, the aperture was adjusted to obtain a suitable depth of field, while a sufficient light input was ensured during exposure.

Image Processing
To counteract the interference of random jet oscillation, 1000 transient shadowgraphs of the stable jet flow with a time interval of 1 ms (one image selected every four frames) were captured and statistically processed.
First, an image-processing method for the jet boundary was established and verified, as shown in Figure 5. Referring to the method proposed by Kannaiyan et al. [33], we established a distribution map of the maximum grayscale SD to identify the boundary of a transverse jet (Max-SD method). Because the liquid column was unstable and oscillated, the zone near the boundary was occupied with liquid and air by turns. The corresponding gray value of the shadow image changed more sharply than that of other locations. Therefore, the maximum of the grayscale SD can be used to indicate the boundary position. The image-processing steps are as follows.
First, statistics were obtained for 1000 shadowgraph frames with an interval of 1 ms and a total duration of 1 s. The grayscale SD (i, j) corresponding to each spatial pixel cell (i, j) at 1000 different times were calculated. Additionally, the spatial distribution map of the SD (i, j) was drawn ( Figure 5a). Subsequently, by scanning point-by-point along the lines parallel to the X-axis, the maximum SD value and its location (i1, j1) were recorded for each line. Finally, the distribution map of the maximum SD positions was drawn in the form of a binary image (Figure 5b). Almost all the maximum-SD points were distributed in the upwind boundary of the transverse jet. This is attributed to two factors: first, the upwind boundary oscillated more intensely than the downwind boundary; second, the grayscale contrast in the downwind boundary was reduced by the nearby surface breakup.
To verify the accuracy of the Max-SD method, the classical method of Sobel edge detection was adopted (Sobel-ED method). Referring to the method reported by Bolszo et al. [36], the average grayscale image was first obtained using the statistics for 1000 shadowgraph frames. An optimal threshold was then determined via the maximum interclass variance method (using "graythresh" function in MATLAB), and the average image was converted into a binary image according to this threshold. Finally, the Sobel operator was used for edge detection to obtain a jet boundary image.
The jet boundaries obtained via the Max-SD and Sober-ED methods were plotted on the average jet shadowgraph (Figure 5c). The boundary curve obtained via the Max-SD method exhibited closer agreement with our expectations and the parabolic characteristics. Compared with the Sobel-ED Energies 2020, 13, 1574 8 of 23 method, the Max-SD method is usually more efficient. This new method does not require an additional process to filter out the interference due to the background and noise and is not affected by threshold selection. The jet boundaries obtained via the Max-SD and Sober-ED methods were plotted on the average jet shadowgraph (Figure 5c). The boundary curve obtained via the Max-SD method exhibited closer agreement with our expectations and the parabolic characteristics. Compared with the Sobel-ED method, the Max-SD method is usually more efficient. This new method does not require an additional process to filter out the interference due to the background and noise and is not affected by threshold selection.
Second, a new identification method for the liquid-column fracture location was established and verified, as shown in Figure 6. The fracture position was obtained using the contour map of the grayscale SD (SD-contour method). Because of the flow instability, the fracture position of the liquid column normally did not correspond to a precise point; statistically, it should correspond to a zone. The liquid in this area changed from a continuous cylinder to discrete ligaments or droplets. This morphological transformation affected the grayscale SD of the jet shadowgraphs. The objective of the SD-contour method is to determine the central position of the aforementioned area.
On the SD contour map, the liquid column is the main body in the area surrounded by the upwind and downwind boundary lines, which corresponds to a series of maximum SD values on the upper and lower sides, respectively. In the outer zone of the boundaries, the broken liquid ligaments or droplets are the main body. The contour line closest to the downwind boundary line connects the surface breakup zone on the right to the liquid-column fracture zone on the top. After the SD contour map is uniformly transformed into a 0-255 grayscale image, this contour line can be used to indicate the starting position of the liquid-column fracture. In this study, the nearest contour line outside the downwind boundary was tracked. The detection point was moved up to the intersection area of the upper and lower boundaries to determine the column fracture location, as shown in Figure 6a. In identifying the fracture position via the SD-contour method, the distance between the two adjacent contour lines affects the location precision. Thus, to reduce the deviation, the contour lines were drawn with a sufficient density, and the number of lines was kept constant (30 lines) in the processing. Second, a new identification method for the liquid-column fracture location was established and verified, as shown in Figure 6. The fracture position was obtained using the contour map of the grayscale SD (SD-contour method). Because of the flow instability, the fracture position of the liquid column normally did not correspond to a precise point; statistically, it should correspond to a zone. The liquid in this area changed from a continuous cylinder to discrete ligaments or droplets. This morphological transformation affected the grayscale SD of the jet shadowgraphs. The objective of the SD-contour method is to determine the central position of the aforementioned area.
On the SD contour map, the liquid column is the main body in the area surrounded by the upwind and downwind boundary lines, which corresponds to a series of maximum SD values on the upper and lower sides, respectively. In the outer zone of the boundaries, the broken liquid ligaments or droplets are the main body. The contour line closest to the downwind boundary line connects the surface breakup zone on the right to the liquid-column fracture zone on the top. After the SD contour map is uniformly transformed into a 0-255 grayscale image, this contour line can be used to indicate the starting position of the liquid-column fracture. In this study, the nearest contour line outside the downwind boundary was tracked. The detection point was moved up to the intersection area of the upper and lower boundaries to determine the column fracture location, as shown in Figure 6a. In identifying the fracture position via the SD-contour method, the distance between the two adjacent contour lines affects the location precision. Thus, to reduce the deviation, the contour lines were drawn with a sufficient density, and the number of lines was kept constant (30 lines) in the processing.
To verify the accuracy of the SD-contour method, a classical method was also adopted. Here, an image-processing program was used to search for the break point, replacing the human eye.
McDonell et al. [41] proposed searching for the fracture location by tracking the jet trajectory on transient shadowgraphs, as shown in Figure 6b. This is called the TT-Search method.
The processing results of the TT-Search method and the proposed SD-Contour method are plotted in Figure 6c. Under different We g values and nanoparticle concentrations, the trends of the two groups of data were consistent, and the differences between them were small. Therefore, the reliability of Energies 2020, 13, 1574 9 of 23 the SD-contour method is comparable to that of the TT-Search method. Furthermore, SD-Contour method does not need additional process to filter out the interference from background, as well as is not affected by threshold selection. To verify the accuracy of the SD-contour method, a classical method was also adopted. Here, an image-processing program was used to search for the break point, replacing the human eye.
McDonell et al. [41] proposed searching for the fracture location by tracking the jet trajectory on transient shadowgraphs, as shown in Figure 6b. This is called the TT-Search method.
The processing results of the TT-Search method and the proposed SD-Contour method are plotted in Figure 6c. Under different Weg values and nanoparticle concentrations, the trends of the two groups of data were consistent, and the differences between them were small. Therefore, the reliability of the SD-contour method is comparable to that of the TT-Search method. Furthermore, SD-Contour method does not need additional process to filter out the interference from background, as well as is not affected by threshold selection.
Third, to verify the convergence of the Max-SD and SD-Contour methods, the effects of the sample size of the shadowgraphs on the processing results were examined, as shown in Figure 7. 200-1000 samples were selected for image processing, and the results for the upwind jet boundary were essentially the same between the two methods when the number of images reached 200. The column fracture position Xb converged when the number of samples reached 600, which was consistent with the trend of the maximum SD. Therefore, in this study, 1000 shadowgraphs with 1ms intervals were selected to ensure good convergence. Third, to verify the convergence of the Max-SD and SD-Contour methods, the effects of the sample size of the shadowgraphs on the processing results were examined, as shown in Figure 7. 200-1000 samples were selected for image processing, and the results for the upwind jet boundary were essentially the same between the two methods when the number of images reached 200. The column fracture position X b converged when the number of samples reached 600, which was consistent with the trend of the maximum SD. Therefore, in this study, 1000 shadowgraphs with 1-ms intervals were selected to ensure good convergence. To verify the accuracy of the SD-contour method, a classical method was also adopted. Here, an image-processing program was used to search for the break point, replacing the human eye.
McDonell et al. [41] proposed searching for the fracture location by tracking the jet trajectory on transient shadowgraphs, as shown in Figure 6b. This is called the TT-Search method.
The processing results of the TT-Search method and the proposed SD-Contour method are plotted in Figure 6c. Under different Weg values and nanoparticle concentrations, the trends of the two groups of data were consistent, and the differences between them were small. Therefore, the reliability of the SD-contour method is comparable to that of the TT-Search method. Furthermore, SD-Contour method does not need additional process to filter out the interference from background, as well as is not affected by threshold selection.
Third, to verify the convergence of the Max-SD and SD-Contour methods, the effects of the sample size of the shadowgraphs on the processing results were examined, as shown in Figure 7. 200-1000 samples were selected for image processing, and the results for the upwind jet boundary were essentially the same between the two methods when the number of images reached 200. The column fracture position Xb converged when the number of samples reached 600, which was consistent with the trend of the maximum SD. Therefore, in this study, 1000 shadowgraphs with 1ms intervals were selected to ensure good convergence.

Flow Visualization
To verify the flow regime of injection from the nozzles used in the present study, a typical shadowgraph of the jet without a crossflow (u ∞ = 0, Re j = 15085) was examined, as shown in Figure 8.

Flow Visualization
To verify the flow regime of injection from the nozzles used in the present study, a typical shadowgraph of the jet without a crossflow (u∞ = 0, Rej = 15085) was examined, as shown in Figure 8. The surface of the liquid column was relatively smooth, with only a slight waveform disturbance. No droplet generation was observed before the liquid traveled a distance 30 times the jet diameter. Thus, the injection was characterized by a nonturbulent flow. In the experiments, as Weg increased from 16.3 to 100.7, the breakup mode of the basic fluid changed from bag breakup to multimode breakup and then to shear breakup. For nanofluids with a concentration of 2.0 wt.%, a critical state between the enhanced capillary breakup and the bag breakup was observed at Weg = 15.3. Overall, the spray morphology of the Al2O3-H2O nanofluids was almost identical to that of the basic fluid, although lower light transmittance was observed for the nanofluid jet. (For more details, a few typical videos of the jets can be seen in the Supplementary Materials.)

Breakup Regimes
To investigate the differences in the primary breakup regime due to the nanoparticles, a Weg-q regime map was drawn for all the liquids, as shown in Figure 9. The results were consistent with The surface of the liquid column was relatively smooth, with only a slight waveform disturbance. No droplet generation was observed before the liquid traveled a distance 30 times the jet diameter. Thus, the injection was characterized by a nonturbulent flow. In the experiments, as We g increased from 16.3 to 100.7, the breakup mode of the basic fluid changed from bag breakup to multi-mode breakup and then to shear breakup. For nanofluids with a concentration of 2.0 wt.%, a critical state between the enhanced capillary breakup and the bag breakup was observed at We g = 15.3. Overall, the spray morphology of the Al 2 O 3 -H 2 O nanofluids was almost identical to that of the basic fluid, although lower light transmittance was observed for the nanofluid jet. (For more details, a few typical videos of the jets can be seen in the Supplementary Materials.)

Breakup Regimes
To investigate the differences in the primary breakup regime due to the nanoparticles, a We g -q regime map was drawn for all the liquids, as shown in Figure 9. The results were consistent with those Energies 2020, 13, 1574 11 of 23 reported by Wu et al. [1]. In the map, the critical We g of water, basic fluid, and nanofluid overlapped, indicating that the nanoparticles had no effect on the transition of breakup modes.
Energies 2020, 13, x FOR PEER REVIEW 11 of 22 those reported by Wu et al. [1]. In the map, the critical Weg of water, basic fluid, and nanofluid overlapped, indicating that the nanoparticles had no effect on the transition of breakup modes. The transition from column breakup to surface breakup is an important topic. In the surface breakup mode, the early stripping of liquid ligaments significantly affects the formation and distribution of fine droplets on the downstream side, as well as the development of surface waves at the upwind boundary.
An auxiliary method was also proposed to distinguish the surface breakup from the column breakup, as shown in Figure 9. According to the SD contour map, all images with similar q values (i.e., the jets with similar trajectories) were combined to compare the breakup starting positions near the downwind boundary. The starting positions could be essentially divided into two types: the positions at a small Weg were significantly higher (corresponding to column breakup), and the positions at a large Weg were significantly lower (corresponding to surface breakup). In some cases, the starting positions obtained from 3-4 repetitive tests were not uniform. Some of them were slightly higher and some were lower, and the cases corresponded to a mixed breakup mode.
Several different empirical correlations of the transition curve have been proposed. However, no relevant analysis of the internal mechanism or mathematical description has been reported thus far. The physical mechanism underlying this transition is discussed below.
Whether surface breakup occurs before liquid-column fracture is determined by the relative speeds of the two processes. When the speed at which ligaments were stripped from the sideways boundary layer exceeded the speed of disturbance development along the liquid jet, a surface breakup mode was observed. When transverse injection was conducted at a small Weg (15 < Weg < 100), the spray exhibited bag breakup or multi-mode breakup. In this case, the surface wave generated by aerodynamic instability gradually developed into a column wave, which affected the entire shape of the jet. Consequently, the jet trajectory changed from a smooth parabola to a U-type wave form with a wavelength of λc (as shown in Figure 10a). The transition from column breakup to surface breakup is an important topic. In the surface breakup mode, the early stripping of liquid ligaments significantly affects the formation and distribution of fine droplets on the downstream side, as well as the development of surface waves at the upwind boundary.
An auxiliary method was also proposed to distinguish the surface breakup from the column breakup, as shown in Figure 9. According to the SD contour map, all images with similar q values (i.e., the jets with similar trajectories) were combined to compare the breakup starting positions near the downwind boundary. The starting positions could be essentially divided into two types: the positions at a small We g were significantly higher (corresponding to column breakup), and the positions at a large We g were significantly lower (corresponding to surface breakup). In some cases, the starting positions obtained from 3-4 repetitive tests were not uniform. Some of them were slightly higher and some were lower, and the cases corresponded to a mixed breakup mode.
Several different empirical correlations of the transition curve have been proposed. However, no relevant analysis of the internal mechanism or mathematical description has been reported thus far. The physical mechanism underlying this transition is discussed below.
Whether surface breakup occurs before liquid-column fracture is determined by the relative speeds of the two processes. When the speed at which ligaments were stripped from the sideways boundary layer exceeded the speed of disturbance development along the liquid jet, a surface breakup mode was observed. When transverse injection was conducted at a small We g (15 < We g < 100), the spray exhibited bag breakup or multi-mode breakup. In this case, the surface wave generated by aerodynamic instability gradually developed into a column wave, which affected the entire shape of the jet. Consequently, the jet trajectory changed from a smooth parabola to a U-type wave form with a wavelength of λ c (as shown in Figure 10a). The liquid column was significantly deflected under the drag force of a crossflow. An additional force caused by relative motion cannot then be neglected, and was relevant to the injection initial velocity uj. Moreover, the windward area of the liquid column gradually increased as the cross section was compressed from a circular shape to an elliptical shape. Therefore, the total aerodynamic force continuously increased along the jet flow direction. A liquid column unit with a single wavelength was selected for force analysis.
When the aerodynamic force (Fa) exceeded the surface tension of both ends (2Fσ), the liquidcolumn fractured. Therefore, the critical state of column breakup should satisfy the following condition: Here, S1 represents the windward area of the liquid column: According to an analysis of the Rayleigh-Taylor instability performed by Ng et al. [42], λc is expressed as: Because the variation range of the current air velocity was relatively small, the drag coefficient CD-round was considered as a constant, with reference to the results of Wu et al. [1].
The drag coefficient CD changes with the deformation of the liquid column. Considering the similarity between the primary breakup of a transverse jet and the secondary breakup of a droplet [7], a simplified model of the droplet drag coefficient was cited, which was proposed by Reitz et al. [43]: Here, CD-round represents the drag coefficient of the liquid column with a circular cross section. It was treated as a constant within the Reynolds-number range of this study.
The equivalent velocity of air crossflow ug at the position of column fracture is given as: The liquid column was significantly deflected under the drag force of a crossflow. An additional force caused by relative motion cannot then be neglected, and was relevant to the injection initial velocity u j . Moreover, the windward area of the liquid column gradually increased as the cross section was compressed from a circular shape to an elliptical shape. Therefore, the total aerodynamic force continuously increased along the jet flow direction. A liquid column unit with a single wavelength was selected for force analysis.
When the aerodynamic force (F a ) exceeded the surface tension of both ends (2F σ ), the liquid-column fractured. Therefore, the critical state of column breakup should satisfy the following condition: (4) Here, S 1 represents the windward area of the liquid column: According to an analysis of the Rayleigh-Taylor instability performed by Ng et al. [42], λc is expressed as: Because the variation range of the current air velocity was relatively small, the drag coefficient C D-round was considered as a constant, with reference to the results of Wu et al. [1].
The drag coefficient C D changes with the deformation of the liquid column. Considering the similarity between the primary breakup of a transverse jet and the secondary breakup of a droplet [7], a simplified model of the droplet drag coefficient was cited, which was proposed by Reitz et al. [43]: Here, C D-round represents the drag coefficient of the liquid column with a circular cross section. It was treated as a constant within the Reynolds-number range of this study. The equivalent velocity of air crossflow u g at the position of column fracture is given as: Here, θ represents the angle between the X-axis and the tangent of the jet trajectory at the fracture point.
According to a mechanical analysis performed by Wu et al. [1], the theoretical jet trajectory (simplified as air flowing around a circular cylinder) can be described as follows: The X-position of the column fracture is X b = C X . Here, both A and C X are constants. Thus: where the constant B is given by Because the variation range of the physical parameters for the current fluids was relatively small, to simplify the analysis, ρ L /ρ g , µ L /µ g , and σ were considered as constant.
By substituting S 1 and u g for the aforementioned expressions into Equation (4), F a can be approximated as follows: Moreover, F σ can be expressed as: Here, l c represents the cross-sectional perimeter of the liquid column. Under significant deformation, the cross section can be simplified as a rectangle: In contrast, surface breakup was regarded as a process of ligaments being stripped from the boundary layer, and the mechanism was similar to Rayleigh breakup. The boundary layer on the front and rear sides developed and accelerated continuously, accompanied by deflection and deformation of the liquid column. When the momentum of the liquid in the boundary layer exceeded the surface tension of a ligament with a size of δ, surface breakup occurred. According to an analysis performed by Sallam et al. [7], the mechanism can be expressed as follows: The surface breakup time t li can be expressed as: Energies 2020, 13, 1574 14 of 23 where B 2 , C l , C t , d, ρ L , and µ g are constants.
The development distance of the boundary layer (u l ·t li ) should be equivalent to the short-axis length d i of the liquid-column cross section. Therefore: Hence, a constant B 1 and B 2 can be introduced. We obtain: Moreover, according to the mass conservation law, the cross-sectional area of the liquid column should change slightly. Subsequently: By introducing Equations (12)- (19) into Equation (3), a final correlation is obtained: where C 0 and C 1 are constants. C 1 is defined as According to the experimental data, the constants can be determined via data fitting of the dimensionless parameter (We g , q) corresponding to the mixed breakup mode. Finally, a semiempirical correlation of the transition line between the column breakup and surface breakup at a small We g is proposed, as follows: ,R 2 = 0.95, We g = 15 − 100 (21) As shown in Figure 9, the current correlation (solid line) is consistent with the correction of Wu et al. [1] (dotted line) for the intersection test cases (We g = 57-100) [1]. Additionally, the critical value of q increases exponentially at low We g , according to the present experiment results and correlation. The jet is close to the enhanced capillary breakup mode when We g < 20. The breakup is mainly dominated by surface tension rather than crossflow, and surface breakup hardly occurs.
In conclusion, nanoparticle additives at low concentrations do not significantly affect the transformation of primary breakup modes.

Penetration Behaviors
The jet penetration behavior is an important factor in combustion chamber design, as it affects the distribution of the spray plume, mixing rate with O 2 , and evaporation rate. To investigate the effects of nanoparticle additives, a comparison of the upwind jet boundary between nanofluids and the basic fluid was performed while q was kept constant. The results are presented in Figure 11

Penetration Behaviors
The jet penetration behavior is an important factor in combustion chamber design, as it affects the distribution of the spray plume, mixing rate with O2, and evaporation rate. To investigate the effects of nanoparticle additives, a comparison of the upwind jet boundary between nanofluids and the basic fluid was performed while q was kept constant. The results are presented in Figure 11. The results indicated that adding a small amount of nanoparticles had almost no effect on the jet trajectory. Thus, for the jet trajectory of nanofluids, the nondimensional ordinate Z is an The results indicated that adding a small amount of nanoparticles had almost no effect on the jet trajectory. Thus, for the jet trajectory of nanofluids, the nondimensional ordinate Z is an exponential function of X and q. Related constants were determined using the MATLAB Curve Fitting Toolbox.
The best-fit correlation for the present experiments was as follows: The correlation coefficient of the fit was 0.98. This fitting result is approximate to the previously reported expression, i.e., Z = A Xq, which is based on mechanical analysis. Both indices of q are 0.5, while the present index of X is slightly smaller owing to a stronger drag force caused by the column deformation. As shown in Figure 12, the fitting curves are consistent with the experimental data of the jet trajectory over a wide q range.
Energies 2020, 13, x FOR PEER REVIEW 15 of 22 exponential function of X and q. Related constants were determined using the MATLAB Curve Fitting Toolbox. The best-fit correlation for the present experiments was as follows: The correlation coefficient of the fit was 0.98. This fitting result is approximate to the previously reported expression, i.e., = , which is based on mechanical analysis. Both indices of q are 0.5, while the present index of X is slightly smaller owing to a stronger drag force caused by the column deformation. As shown in Figure 12, the fitting curves are consistent with the experimental data of the jet trajectory over a wide q range.

Column Fracture Regularity
To confirm the facticity of the additional physical breakup mode for the nanofluids, the column fracture positions were plotted in different graphs. As shown in Figure 13, the X-positions of column fracture (Xb) were scattered within the range of 9.9 ± 2.4 under the present conditions of Weg = 25-92 and q = 32-186. Compared with the Xb value of 8.06 ± 1.46 reported by Wu et al. [1], there was a deviation, which was mainly attributed to two factors: the different breakup modes at a small Weg and the different methods of selecting the fracture location.

Column Fracture Regularity
To confirm the facticity of the additional physical breakup mode for the nanofluids, the column fracture positions were plotted in different graphs. As shown in Figure 13, the X-positions of column fracture (X b ) were scattered within the range of 9.9 ± 2.4 under the present conditions of We g = 25-92 and q = 32-186. Compared with the X b value of 8.06 ± 1.46 reported by Wu et al. [1], there was a deviation, which was mainly attributed to two factors: the different breakup modes at a small We g and the different methods of selecting the fracture location.

Column Fracture Regularity
To confirm the facticity of the additional physical breakup mode for the nanofluids, the column fracture positions were plotted in different graphs. As shown in Figure 13, the X-positions of column fracture (Xb) were scattered within the range of 9.9 ± 2.4 under the present conditions of Weg = 25-92 and q = 32-186. Compared with the Xb value of 8.06 ± 1.46 reported by Wu et al. [1], there was a deviation, which was mainly attributed to two factors: the different breakup modes at a small Weg and the different methods of selecting the fracture location. We g = 25.5 ± 0.5 We g = 48.5 ± 0.6 We g = 86. 8

Volume Fraction of Nanoparticles, ϕ
We g = 25.5 ± 0.5 We g = 48.9 ± 0.4 We g = 87.8 ± 0.6 Figure 13. X-positions of column fracture in different cases: (a) Different We g and q; (b) Different We g and ϕ.
The fracture location is mostly selected on the upwind boundary of the jet in the literature, while the position was selected in the middle area of the upper and lower jet boundaries in this study. From the viewpoint of applications, it is sometimes more meaningful to detect the location in the middle area, which can indicate the central position with the maximum fracture probability of a random-oscillation liquid column.
The results indicate that X b increased with a decrease in We g but had almost no correlation with q ( Figure 13a). Z b increased with q but had almost no correlation with We g (Figure 14a). The fracture distance of the liquid column decreased as the nanoparticle concentration increased. Relative to the basic fluid, the fracture location X b of nanofluids descended more rapidly at a larger We g number, and the location Z b descended more rapidly at a smaller q number. This is consistent with the conclusions of Kannaiyan et al. [33] and Li et al. [35]. When the mass concentration of nano-alumina was 2.0 wt.% (corresponding volume fraction was 0.51%), the fracture distance X b in the present test was shortened by approximately 20% (We g = 92, Re j = 19074) compared with that of the basic fluid, and the breakup length was reduced by approximately 17% (Re p = 27450) reported by Kannaiyan et al. [33], which was on the same order of magnitude. and the location Zb descended more rapidly at a smaller q number. This is consistent with the conclusions of Kannaiyan et al. [33] and Li et al. [35]. When the mass concentration of nano-alumina was 2.0 wt.% (corresponding volume fraction was 0.51%), the fracture distance Xb in the present test was shortened by approximately 20% (Weg = 92, Rej = 19074) compared with that of the basic fluid, and the breakup length was reduced by approximately 17% (Rep = 27450) reported by Kannaiyan et al. [33], which was on the same order of magnitude. In contrast to the present study, Kannaiyan et al. used a swirling spray nozzle, and the reduction of the breakup length was attributed to the decrease in the nanofluid surface tension. Possibly because of different basic fluids, we found that the surface tension of water was slightly increased by adding nanoparticles, which is consistent with the results of Saad and Li [37]. Therefore, there should be other reasons why the nanoparticles change the fracture positions in transverse injection. According to [29,36], it is likely that the cavitation in the jet core is changed by the nanoparticles. It is necessary to determine whether the cavitation is enhanced by the nanoparticles, which would cause internal interface breakup. Therefore, a round jet test (u∞ = 0) was performed using a twicemagnification transparent nozzle made of plexiglass.
As shown in Figure 15, in the absence of visible bubbles upstream the nozzle exit, no bubbles were observed inside the basic fluid jet. However, for nanofluids with a concentration of 0.5 wt.%, a stream of bubbles with sizes of 0.2-0.4 mm appeared in the jet core region. The bubbles emerged from the jet inner zone with enhanced disturbance near the nozzle (blue rectangle in Figure 15b). Thus, the addition of nanoparticles increased the amount of nucleation sites and promoted cavitation in the jet. In contrast to the present study, Kannaiyan et al. used a swirling spray nozzle, and the reduction of the breakup length was attributed to the decrease in the nanofluid surface tension. Possibly because of different basic fluids, we found that the surface tension of water was slightly increased by adding nanoparticles, which is consistent with the results of Saad and Li [37]. Therefore, there should be other reasons why the nanoparticles change the fracture positions in transverse injection. According to [29,36], it is likely that the cavitation in the jet core is changed by the nanoparticles. It is necessary to determine whether the cavitation is enhanced by the nanoparticles, which would cause internal interface breakup. Therefore, a round jet test (u ∞ = 0) was performed using a twice-magnification transparent nozzle made of plexiglass.
As shown in Figure 15, in the absence of visible bubbles upstream the nozzle exit, no bubbles were observed inside the basic fluid jet. However, for nanofluids with a concentration of 0.5 wt.%, a stream of bubbles with sizes of 0.2-0.4 mm appeared in the jet core region. The bubbles emerged from the jet inner zone with enhanced disturbance near the nozzle (blue rectangle in Figure 15b). Thus, the addition of nanoparticles increased the amount of nucleation sites and promoted cavitation in the jet. Additionally, the generation of bubbles caused interface breakup inside the jet, accelerating the fracture of the liquid column.
Energies 2020, 13, x FOR PEER REVIEW 17 of 22 Additionally, the generation of bubbles caused interface breakup inside the jet, accelerating the fracture of the liquid column. Accordingly, a new nondimensional quantity associated with nanoparticles should be introduced for the correction of the column fracture position. In this paper, it is defined as PV, which is mainly determined by two factors. The first consists of basic characteristics of particle phase, including the particle concentration, true density, size distribution, and thickness of the solvation layer. The second consists of characteristics relevant to the jet flow field, mainly including the Peclet number (Pep) and Reynolds number of the particles (Rep), which represents the dispersion uniformity and following feature.
At the nanoscale, both the particle dispersion inhomogeneity and moving relaxation are low. Thus, Pep and Rep can be ignored. Among the remaining factors, the particle concentration was the most influential. Referring to the corrected form of the nanofluid viscidity based on the equivalent volume concept, PV is expressed as follows: Hence, PV represents the extent to which the nanofluids deviate from the pure liquid. For the present Al2O3-H2O nanofluids, C was calculated as = 1.46 × 3 × 1.019 = 9.5.
According to the foregoing analysis, the X and Z coordinates of the column fracture position can be expressed as follows: According to the experimental data, the following corrections are established: where the correlation coefficients of the fits were both 0.96. As shown in Figure 16, the fitting curves agreed well with the experimental data. For the addition of nanoparticles, the X-position of the liquid-column fracture was proportional to PV 3 , and the Z position was proportional to PV 2 . Accordingly, a new nondimensional quantity associated with nanoparticles should be introduced for the correction of the column fracture position. In this paper, it is defined as P V , which is mainly determined by two factors. The first consists of basic characteristics of particle phase, including the particle concentration, true density, size distribution, and thickness of the solvation layer. The second consists of characteristics relevant to the jet flow field, mainly including the Peclet number (Pe p ) and Reynolds number of the particles (Re p ), which represents the dispersion uniformity and following feature.
At the nanoscale, both the particle dispersion inhomogeneity and moving relaxation are low. Thus, Pe p and Re p can be ignored. Among the remaining factors, the particle concentration was the most influential. Referring to the corrected form of the nanofluid viscidity based on the equivalent volume concept, P V is expressed as follows: Hence, P V represents the extent to which the nanofluids deviate from the pure liquid. For the present Al 2 O 3 -H 2 O nanofluids, C was calculated as C = 1.46 3 × 3 × 1.019 = 9.5.
According to the foregoing analysis, the X and Z coordinates of the column fracture position can be expressed as follows: According to the experimental data, the following corrections are established: where the correlation coefficients of the fits were both 0.96. As shown in Figure 16, the fitting curves agreed well with the experimental data. For the addition of nanoparticles, the X-position of the liquid-column fracture was proportional to P V 3 , and the Z position was proportional to P V 2 . Energies 2020, 13

Summary and Conclusions
For the first time, a nanofluid jet injected into a subsonic gaseous crossflow was investigated with a focus on three aspects: the primary breakup regimes, penetration behaviors, and column fracture features. A series of experiments were performed at room temperature and the atmospheric pressure. The experimental setup involved one conical nozzle, three sets of injection pressures, five different concentrations of nanoparticles, and seven crossflow rates. First, an improved imageprocessing method was developed for statistical analysis of transverse injection. According to the distribution map of the grayscale SD, the jet trajectory, column fracture position, and morphology of the column/surface breakup were obtained. Second, for the first time, a simplified mathematical description of the transition mechanism from column breakup to surface breakup at a small Weg was established. The correlation of the transition line, which was obtained via morphological modeling and force analysis, agreed well with the experimental correlation reported by Wu et al. [1]. Finally, the effects of adding alumina nanoparticles on the upwind jet trajectory and column fracture position were investigated through comparisons between nanofluids and the basic fluid. On the basis of research method for emulsion injection and the equivalent volume concept for fitting the nanofluid viscidity, a new nondimensional parameter called PV was introduced. The following conclusions are drawn.
In the case of a small Weg, the fracture of the liquid column was caused by column wave expansion with aerodynamic instability, while the surface breakup was dominated by the boundarylayer stripping process caused by the crossflow. When the speed at which ligaments were stripped from the boundary layer exceeded the speed of disturbance development along the liquid jet, the surface breakup mode was observed. When the momentum of the liquid in the boundary layer exceeded the surface tension of the ligament with a size equivalent to δ, droplets were stripped sideways from the column. The following semiempirical correlation of the transition line between the column breakup and surface breakup was proposed: With the addition of a small amount of nanoparticles (φ < 2.0 wt.%), the surface tension increased slightly, while the viscosity increased significantly (by up to 23%). The changes in these physical properties hardly affected the breakup regimes and jet trajectory. The trajectory correlation of the transverse injection for nanofluids can be expressed as follows:

Summary and Conclusions
For the first time, a nanofluid jet injected into a subsonic gaseous crossflow was investigated with a focus on three aspects: the primary breakup regimes, penetration behaviors, and column fracture features. A series of experiments were performed at room temperature and the atmospheric pressure. The experimental setup involved one conical nozzle, three sets of injection pressures, five different concentrations of nanoparticles, and seven crossflow rates. First, an improved image-processing method was developed for statistical analysis of transverse injection. According to the distribution map of the grayscale SD, the jet trajectory, column fracture position, and morphology of the column/surface breakup were obtained. Second, for the first time, a simplified mathematical description of the transition mechanism from column breakup to surface breakup at a small We g was established. The correlation of the transition line, which was obtained via morphological modeling and force analysis, agreed well with the experimental correlation reported by Wu et al. [1]. Finally, the effects of adding alumina nanoparticles on the upwind jet trajectory and column fracture position were investigated through comparisons between nanofluids and the basic fluid. On the basis of research method for emulsion injection and the equivalent volume concept for fitting the nanofluid viscidity, a new nondimensional parameter called P V was introduced. The following conclusions are drawn.
In the case of a small We g , the fracture of the liquid column was caused by column wave expansion with aerodynamic instability, while the surface breakup was dominated by the boundary-layer stripping process caused by the crossflow. When the speed at which ligaments were stripped from the boundary layer exceeded the speed of disturbance development along the liquid jet, the surface breakup mode was observed. When the momentum of the liquid in the boundary layer exceeded the surface tension of the ligament with a size equivalent to δ, droplets were stripped sideways from the column. The following semiempirical correlation of the transition line between the column breakup and surface breakup was proposed: , R 2 = 0.95, We g = 15 − 100 With the addition of a small amount of nanoparticles (ϕ < 2.0 wt.%), the surface tension increased slightly, while the viscosity increased significantly (by up to 23%). The changes in these physical properties hardly affected the breakup regimes and jet trajectory. The trajectory correlation of the transverse injection for nanofluids can be expressed as follows: The addition of nanoparticles promotes cavitation and bubble generation inside the liquid jet, resulting in an additional primary breakup mode called internal interface breakup. Consequently, the liquid column fractures, and its horizontal length is reduced by up to 20%. The X and Z coordinates of the column fracture location are expressed as follows: The primary breakup of nanofluid jets in a crossflow at a small We g is depicted in Figure 17. After being ejected from the nozzle, the liquid column is deflected and compressed by the crossflow. Two types of breakups occur, accompanied by an increase in instability: (i) the formation of a column wave when the threshold of the amplitude of the surface wave is reached, which causes the liquid column to fracture, and (ii) cavitation caused by disturbance in the nanofluid jet, which is followed by internal interface breakup. In some cases of a high We g or q, the surface breakup occurs in the form of boundary-layer stripping, similar to Rayleigh breakup, which results in smaller droplets. The addition of nanoparticles promotes cavitation and bubble generation inside the liquid jet, resulting in an additional primary breakup mode called internal interface breakup. Consequently, the liquid column fractures, and its horizontal length is reduced by up to 20%. The X and Z coordinates of the column fracture location are expressed as follows: The primary breakup of nanofluid jets in a crossflow at a small Weg is depicted in Figure 17. After being ejected from the nozzle, the liquid column is deflected and compressed by the crossflow. Two types of breakups occur, accompanied by an increase in instability: i) the formation of a column wave when the threshold of the amplitude of the surface wave is reached, which causes the liquid column to fracture, and ii) cavitation caused by disturbance in the nanofluid jet, which is followed by internal interface breakup. In some cases of a high Weg or q, the surface breakup occurs in the form of boundary-layer stripping, similar to Rayleigh breakup, which results in smaller droplets.    Acknowledgments: A special thanks to Wei Fu and Liufang Yu, for their suggestions on test method and supports in image device.

Conflicts of Interest:
The authors declare no conflict of interest. liquid column average drag coefficient width of liquid column cross-section, m liquid jet diameter at the nozzle exit, m equivalent diameter of water molecules, m length of liquid column cross-section, m mean diameter of monodispersed nanoparticles, m thickness of solvation layer, m aerodynamic drag force, N surface tension, N characteristic height of column cell, m cross-section perimeter of liquid column, m Ohnesorge number, µ L / ρ L d j σ particle Peclet number, d p u L − u p /α p nanofluids-jet parameter, 1 − C·ϕ jet momentum flux ratio, ρ L u 2 j /ρ g u 2 ∞ jet Reynolds number, ρ L u j d j /µ L particle Reynolds number, ρ L d p u L − u p /µ L windward area of liquid column, m 2 rounded cross-section area of liquid column, m 2 elliptic cross-section area of liquid column, m 2 onset of surface breakup time, s air crossflow velocity, m/s liquid velocity at the nozzle exit, m/s mean velocity of ligament in boundary layer, m/s mass fraction of nanoparticles, wt.% gas phase Weber number, ρ g u 2 ∞ d j /σ non-dimensional x-coordinate, x/d j non-dimensional z-coordinate, z/d j thickness of liquid boundary layer, m angle between X axis and tangent of jet trajectory column wave length, m dynamic viscosity, Pa·s kinematic viscosity, m 2 /s density, kg/m 3 surface tension, N/m volume fraction of nanoparticles, vol.% Subscripts b bf g j L nf p ∞ location of liquid column fracture basic fluid gas jet liquid nanofluids nanoparticle ambient gas