Research on the Utilization of Saline Alkali Water Resources Based on Two-Phase Flow

: In order to reasonably use solar energy to solve problems such as land desertiﬁcation and soil salinization in Southern Xinjiang, this paper proposes a system combining photothermal and ﬂash evaporation technologies, which use local saline water for desalination treatment to achieve secondary utilization of water resources. Firstly, we introduce the whole system of the photovoltaic desalination plant. As an important heat-collecting element of the system, the solar tube is the key to whether this plant can work e ﬃ ciently. Then, we carry out the detection and theoretical derivation of data along the tube. We establish a two-phase ﬂow model of saline water in the tube, considering convective heat transfer, and deﬁne the formula of the heat collecting e ﬃ ciency factor. Finally, based on iterative calculation, the temperature trend of the tube and the change law of the two-phase ﬂow are obtained, and the ecological and economic beneﬁts and energy e ﬃ ciency of the system are analyzed.


Introduction
Southern Xinjiang is an area in China that is blessed with abundant solar energy, but has the issue of severe saline alkaline water pollution. As such, it is of great significance to make full use of the available solar energy to solve water pollution problems. The Xinjiang region has the characteristics of long days and sufficient solar radiation. As shown by the annual radiation of the Aksu meteorological station (Figure 1), there is 5467 MJ/m 2 solar radiation per year on average. (There are many monitoring stations in various places, but the Aksu meteorological station is in the Tarim Basin and located in an area that can best represent the local solar radiation). At the same time, the local solar radiation intensity changes with the seasons as shown in Figure 2.  The direct use of solar energy to heat salt water for desalination can meet the growing demand for water resources. Therefore, solar water desalination technologies may become one of the basic water resources for future survival and sustainable development. According to a survey, by 2025, 14% of the global population will be forced to use or drink desalinated water [1]. Solar thermal desalination is a medium-low temperature technology, which poses frontier research problems that have been recognized by the world's new energy and renewable energy departments [2]. Some common desalination devices already exist, such as humidification and dehumidification (HDH), membrane distillation (MD), multistage distillation (MSD), and multieffect evaporation (MED) devices [3][4][5][6]. The use of solar energy to achieve the desalination of saline water involves two-phase flow and boiling heat transfer phenomena. Two-phase flow has been widely used in various industrial fields, such as air conditioning [7], humidification [8], water cooling [9], power engineering applications [10], and industrial hot water heating [11]. Furthermore, the number of calculation formulas for flow heat transfer coefficients in the tube has been constantly increasing [12][13][14][15]. Badar et al. [16][17][18] used a steady-state heat transfer model to analyze the performance of single-fluid and dual-fluid pairs of collector tubes under the same conditions. Thermal distillation in conventional desalination plants often releases large amounts of greenhouse gases due to the consumption of large amounts of energy, leading to global warming [19]. At present, many countries are changing the development direction of seawater desalination technology, focusing on environmentally friendly water treatment processes. Elsarrag et al. [20] studied the advantages of combining solar pools with solar evaporators, recycling the high heat waste brine (brine is a waste product in other desalination systems) as a steam heat source, and using the surface water of the solar pool as cooling water to be used in the condensing system. Their study pointed out that the desalination cost of the solar pond was smaller than that of other desalination methods, and the processing capacity is large. Due to its low efficiency, a solar still cannot be operated on a large scale. The efficiency of solar stills can be sharply increased by adding solar collectors, while most of the energy losses occur as a dissipation of heat, lacking effective energy recycling. Despite the fact that multistage flash distills saline water by flashing a portion of the water into steam in multiple stages, essentially reusing the dissipation energy, its cost performance is still poor.
Mamouri et al. [21][22][23] developed a salt water desalination device that combines a solar collector tube with an inclination angle of 35 • C and an evaporator, and studied the influence of heat collection temperature and solar radiation on the performance of the collector tube. On this basis, Mosleh et al. [24] added a trough-type reflector to increase the absorption rate of the collector tube; Li et al. [25] ventilated the bottom of the brine evaporation tank to expand the gas-liquid contact area to increase the amount of water vapor generated. Al-harahsheh et al. [26] added a phase change heat storage tank to improve the salt water desalination device to effectively use heat energy.
These technologies differ greatly in terms of energy consumption, cost, efficiency, and output. The solar distillation process is simple, and the water production efficiency is too low due to the mode of only relying on the natural heating, evaporation, and condensation of the waterbody and the large area. In order to achieve stable working operation, multistage flash evaporation requires additional heat sources, resulting in higher overall energy consumption. The high cost of humidification and dehumidification equipment, the complicated process flow leads to large heat loss, and a relatively high relative humidity is required. The output of water will increase with the increase of humidity. Moreover, the energy transfer of these desalination systems is mostly one-way, and the remaining heat energy after partial heating is not utilized.
Most of Xinjiang is remote with insufficient rainfall throughout the year. Considering the local natural conditions and ensuring a certain amount of water supply, the energy utilization is maximized as much as possible without adding additional heat sources. The system combines the advantages of each process, under the premise of controlling the floor space, input cost, and adapting to the local environment, it achieves a stable water output and provides a method for solving the salt water desalination in remote and arid areas. The whole system mainly includes a water storage tank, condenser, heater, and evaporator. The heater consists of a series of horizontally arranged solar vacuum glass tubes, arranged in parallel. The heat exchanger consists of a condenser, with a heat exchange tube located in the condenser. The cost of the system mainly includes equipment, transportation, construction, and maintenance. The total is about 160,000 USD. The entire system station covers about 400 m 2 .
Based on the current situation in Xinjiang and the traditional desalination methods, this paper proposes a new desalination system adapted to the local area, described in Section 1. Section 2 introduces the components of the desalination system, and introduces in detail the theoretical analysis method for the core heating element of the system and the establishment of the two-phase flow model. Through field experiments, the temperature changes along the tube were obtained under different flow rates and described in Section 3. The comparison with the efficiency curve of the established two-phase flow model proved the applicability of the model. Through numerical simulation, plots of heat collection efficiency and temperature change under different mass flow rates are drawn. Section 4 explains the relationship between the heat collecting efficiency of the tube and the change of temperature and flow in combination with the analysis of the internal heat transfer mechanism. In the comparison of the advantages and limitations of the traditional desalination process, the practical value of the system is reflected, and the benefits of the system are analyzed from three aspects. Section 5 describes the feasibility of the system for solving the problem of local saline-alkali water desalination, and summarizes the expression of the heat collection efficiency factor and the gas-liquid phase transformation formula established in the theoretical analysis of the two-phase flow of the tube. This article aims to combine the existing desalination technologies to propose a desalination system suitable for Southern Xinjiang, and provide new ideas for solving water resources problems in remote areas.

Materials and Methods
The benefits of this system are as follows. The heater is made of quartz glass and does not react with saline-alkaline water. Therefore, no more pollution is produced. The glass tube has inner tubes and outer tubes. The inner tube has a heat absorbent coating to absorb most of the solar heat, and the vacuum treatment between the outer and inner tubes can prevent heat loss to the utmost extent. Moreover, the glass tubes are connected in series, which can increase the area for absorbing sunlight and imp heating efficiency. The inside of the evaporator is evacuated, with a flashing function, which can lower the boiling point of saline-alkaline water and evaporate it into steam in a short time. The hot saline-alkaline water at the outlet of the vacuum glass tube enters the evaporator in the form of a jet stream. The heat exchanger consists of a condenser containing heated steam and heat exchange tubes containing cold saline-alkaline water. The cold saline-alkaline water and steam exchange heat at the surface of the exchange tubes, such that the energy is not easily lost, result from the large heat exchange area. Making full use of the energy of hot steam to heat the cold saline water leads to higher energy efficiency. The schematic diagram of the system as shown in Figure 3. Evaporator: The evaporator consists of an evaporation chamber, where the chamber is in a low vacuum state at 0.7-0.8 bar to achieve the purpose of lowering the boiling point of salt water. There is a spin nozzle inside the evaporator, which uses turbulent diffusion to improve the efficiency of gas-liquid separation. The inlet of the evaporator is connected to the outlet of the vacuum glass tube, while the top of the evaporation chamber is connected to the inlet of the condenser. As evaporation continues to occur, fresh water is continuously produced and the salty aqueous solution in the lower part of the flash tank continuously becomes more concentrated, finally reaching the critical salinity (25%) for meeting the salt conditions. At this point, the circulation pump and the return valve are closed, the discharge valve is opened, and the concentrated saline solution in the lower part of the flash tank is discharged. The structure diagram of the evaporator and salt pond are shown in Figure 4. Condensation device: Cold salt water enters the device from the lower part and flows out from the upper part. The cold salt water acts as a continuous cold source for cooling the hot steam into fresh water, while hot steam enters from the upper inlet as a continuous heat source for preheating the cold salt water; finally, the fresh water is collected and removed from the lower outlet. The design of the condenser is shown in Figure 5. Solar vacuum tube: The heat-drying device is made up of many horizontal vacuum tubes, arranged in parallel. The solar vacuum tube used in this study consists of two coaxial quartz tubes, which are divided into an inner tube and an outer tube. The inner tube is coated with a heat absorbing layer, and a vacuum is applied in the gap between the inner and outer tubes. The solar vacuum tube can continuously absorb solar heat and reduce heat loss. By connecting in series, the area of received sunlight can be increased, such that the salty water temperature in the tube increases with distance. The inner and outer tube structure of the solar vacuum tube is shown in Figure 6 and the arrangement of the tubes is shown in Figure 7.

Concentrating Heat Mechanism of Tubes
Based on the analysis method of Duffie and Beckman [27], we propose an improved model to meet the requirements of the coaxial effect and two-phase flow research of solar vacuum tubes.
Take a section of tubes for analysis, as shown in Figure 8. The solar radiation energy absorbed per unit area of the quartz inner tube is S = ατ·I, where I is the radiation intensity, α is the absorption rate, and τ is the transmittance of the quartz outer tube. When salty water flows through the tubes, it absorbs most of the heat, while a small heat loss is recorded as q loss . The heat absorbed by the tubes per unit area is: q abs = ατ·I − q loss .
(1) The temperature of the endothermic coating gradually rises under the solar radiation, forming a temperature difference with the ambient environment. Heat flows from the endothermic coating through the gas between the inner and outer quartz tubes and then through the outer quartz tube to the surrounding environment. This heat loss is recorded as q loss .
In order to calculate the magnitude of the heat loss, the value of the total thermal resistance U needs to be determined. Therefore, assuming that the ambient temperature of the left and right semicircles is the same, the expression of U can be analyzed by considering the tube as an axisymmetric body. The heat loss mechanism of the tube is shown in Figure 9: This heat loss is caused by radiative heat transfer and convective heat transfer. We obtain the heat loss coefficient of the semicircular tube U as: where h rig is the radiation heat transfer coefficient of the endothermic coating of the outer quartz tube (the subscript r indicates radiation), h cig is the sum of the convective heat transfer coefficient and thermal conductivity of the endothermic coating of the outer quartz tube (the subscript c indicates convection and heat conduction, and the subscript ig indicates inner tube to outer tube), h rga is the radiation heat transfer coefficient of the outer quartz tube to the air (the subscript r, again, indicates radiation), and h cga is the convective heat transfer coefficient of the outer quartz tube to the air (the subscript c, again, indicates convection and heat conduction, and the subscript ga indicates outer tube to the air). When the vacuum between the inner and outer quartz tubes is high enough, convective heat transfer and heat conduction can be isolated (i.e., h cig ≈ 0). Radiative heat transfer is only considered in the heat transfer process of the inner quartz tube to the outer quartz tube. The radiative heat transfer coefficient between the inner quartz tube and the outer quartz tube is calculated according to the formula of Duffie and Beckman [28]: where σ is the Boltzmann constant W/ m 2 ·K 4 ε p is the emissivity of the endothermic coating, ε g is the emissivity of the outer quartz tube, A g is the surface area of the outer quartz tube, F ig is the angle of view factor of the endothermic coating of the outer quartz tube, T a is the ambient temperature (the subscript a represents air), T abs is the temperature of the endothermic coating, and T g is the temperature of the outer quartz tube. The radiative heat transfer coefficient of the outer quartz tube to the ambient environment is: The convective heat transfer coefficient of the outer quartz tube to the ambient environment is: where V air is the local wind velocity (m/s) and the units of h cga are W/ m 2 ·K . Under steady-state conditions, the heat transfer per unit area from the endothermic coating to the outer quartz tube is the same as the heat loss per unit area from the endothermic coating to the ambient environment. The net heat increase Q f net of salty water in the vacuum tube can be expressed as: where F c is the heat collecting efficiency factor of the tubes, U is the heat loss coefficient of the whole tubes, and A abs is the heat absorbing area of the tube. In the end, we get the water temperature at the outlet of the tube as: where T f i is the water temperature at the inlet, T f o is the water temperature at the outlet (subscript i represents the inlet, subscript o represents the outlet), and C f is the specific heat capacity of salty water (subscript f represents the waterbody). Salty water enters the inner tubes with mass flow rate . m (for convenience, hereinafter referred to as the flow rate) and temperature T i , with the temperature of salty water leaving the tube section being T o .

Establishment of the Mathematical Model of Two-Phase Flow in Tubes
Solar radiant energy is continuously absorbed by the tubes, and the waterbody is constantly heated. If the tubes are long enough, the average water temperature at the local position of the tubes will reach the saturation temperature under the corresponding pressure, and flow boiling will occur. In this case, there are two states in the fluid in the tubes, and the ordinary single-phase flow mathematical model is not applicable. It is better to use a two-phase flow model to describe the motion of fluid in the tubes. The heat transfer rate of a two-phase flow is usually much higher than that of a normal single-phase flow. However, a two-phase flow is much more complicated, in terms of physical phenomena, than a single-phase flow, and multiple coefficients need to be considered to predict the heat transfer. It is very difficult to directly study vaporization and heat transfer in the tubes by using a two-phase differential equation. It is, therefore, necessary to introduce the approximate hydrodynamic assumptions: (1) The heat flux in the microcell is the same.
(2) The heat radiation characteristics of the endothermic coating are independent of the initial water temperature.  7) When the water flowing in the tubes reaches the boiling point, the water temperature in the inner tubes is equal to the saturation temperature, and the heat transfer coefficient is equal to the two-phase boiling heat transfer coefficient.
Taking a microbody of salty water in the flow direction for thermal energy balance analysis (as shown in Figure 9), the average temperature is T f . The microbody absorbs heat from the coating and upstream water, and transfers heat downstream. The heat transferred from the endothermic coating to the salt water is calculated based on the heat transfer coefficient, h f , of the coating to the waterbody. According to Fourier's law, the heat passing through the cross-section of the waterbody per unit time is proportional to the rate of temperature change and the cross-sectional area perpendicular to the cross-sectional direction. The direction of heat transfer is opposite to the direction of temperature rise. The heat transfer of water in the flow direction can be expressed as −k f dT f dx , where k f is the thermal conductivity of water. If vaporization occurs during the flow to produce water vapor, the change in thermal energy in the gas phase is represented by the difference in enthalpy between the outlet and the inlet of the microelement. As shown in Figure 10. The input heat minus the lost heat is equal to the net heat absorption of salty water. According to the heat absorption and heat transfer of the coating, the thermal energy balance equation on the unit coating area is established as: where m abs is the mass of the endothermic coating, c abs is the specific heat capacity of the endothermic coating (the subscript abs represents the endothermic coating), h 1 is the total heat transfer coefficient of the endothermic coating to the outer tubes, T g is the temperature of the outer tube, U is the total heat transfer coefficient of the heat absorbing coating to the environment, and T a is the ambient temperature. When salty water flows through the inner tube, the heat of the tubes is absorbed, the water temperature rises, and the heat energy balance equation of salty water per unit area is: where .
m∆H is the change of enthalpy in salt water. When the fluid in the tubes is at rest, the constant input of heat causes the water temperature to continuously increase, causing it to be in a nonconstant state. When the water flows in the tubes, the continuously incoming heat is continuously taken away by the water flow, such that the temperature on the local section does not rise with time, showing a near-constant state. The heat absorption of the static fluid in the tubes belongs to the research of time category. If the fluid flows, the heat in the section is taken away, the section temperature is maintained at a stable value, the heat absorbed by the fluid in the next section causes the temperature to be higher than the previous one, and spatial temperature change occurs. This also converts the change of static water temperature with time into the change of flowing water temperature with space, greatly shortening the heating time. As the temperature along the process increases, it is convenient for us to observe the phenomena of vaporization, boiling, heat accumulation, heat transfer, and heat loss in the waterbody, which are also convenient for gas-liquid separation and condensate collection.

Heat Collecting Efficiency Factor of Tubes
At both ends of the tubes, there are elbows that are symmetric both upward and downward, with no vacuum insulation. Suppose the two elbows are 2L apart, the thickness of quartz inner tubes is δ, and the thermal conductivity is k i . As the endothermic coating is a good conductor, the temperature gradient along the cross-section of the endothermic coating is negligible. As the elbows are exposed to the environment, causing a large heat loss, the temperature of the endothermic coating at the bare place is less than that in the center of the tube section, and part of the heat is transferred from the center to both sides. According to the principle of the heat sink, the temperature at the center of the tubes is the highest, both sides are symmetrically dissipated, and the tubes are symmetric about the center. Therefore, half of the tube is selected for thermodynamic analysis, as shown in Figure 11. The length of the tube section in Figure 11a is L and the temperature of the endothermic coating at the junction of the straight tubes and elbows is recorded as T L . The microbody of the endothermic coating in the flow direction is shown in Figure 11b.
According to Fourier's law, the heat passing through a section in a unit time is proportional to the rate of temperature change, which is perpendicular to the direction of the section. The direction of heat transfer is opposite to that of temperature rise. As a result, heat conduction of the endothermic coating in the flow direction can be expressed as −k i δ dT dx , in which k i is the thermal conductivity of the endothermic coating. The thermal energy conservation equation of the microbody is recorded as: By solving this equation, we can define the heat collecting efficiency factor F c of the endothermic coating as: The tube system is made up of a series of vacuum quartz glass tubes connected in series and the plane layout is in the form of up and down symmetry. Therefore, the heat collecting efficiency factor of each tube section is the same:

Iterative Calculation of Salt Water Temperature
When the salty water flows through the tubes, the total heat gain of the tube system and the water temperature distribution along the path can be calculated iteratively. We calculate each unit according to the conservation of thermal energy.
The elbow is the n+1 unit, and there is only heat loss at the elbow due to the absence of the endothermic coating and the vacuum insulation. The inlet and outlet water temperatures of unit 1 are T f i1 and T f o1 , the inlet and outlet water temperatures of unit 2 are T f i2 and T f o2 , and so on (i.e., the inlet and outlet water temperatures of unit n + 1 are T f i,n+1 and T f o,n+1 , respectively). As shown in Figure 12.
In the iterative calculation, when the salty water has not reached the boiling temperature, T sat , under the pressure of the tubes, it is considered that the tube is a single-phase flow. We use the single-phase flow heat transfer coefficient h f to substitute Formula (12) in order to calculate the heat energy absorbed by the unit waterbody. When the unit average water temperature T f m exceeds the saturation temperature T sat , the T f m of this unit is set to T sat . The net heat increment of gas and water is calculated according to the two-phase flow boiling heat transfer coefficient h f ,boil . We use the method of Shah [29] to estimate the two-phase heat transfer coefficient: (1) Calculate the dimensionless parameter N: (2a) For the dimensionless parameter N > 1.0, we calculate the heat transfer coefficients h NcB and h c , according to the following formula (the subscript NcB stands for nucleate boiling, subscript c stands for convection) and choose the larger of the two as h f ,boil : The heat flux of the inner tube wall is: where h f is the heat transfer coefficient of the liquid phase, the flow rate of the liquid phase is G(1 − α), and α is the mass fraction of the gas. The heat transfer coefficient of the liquid phase is: in which h f ,in is obtained from the Nusselt number of a fully developed laminar flow of single-phase fluid in tubes: where k f is the thermal conductivity of the fluid and N u is the Nusselt number of the liquid laminar flow. (2b) For the dimensionless parameter 1.0 > N > 0.1, the calculation of h c still uses Formula (21). In the bubble suppression mechanism, it was determined that h NcB depends on the following formula: where F is the efficiency factor of the heat sink. We choose the larger of h c and h NcB as h f ,boil .
(2c) For the dimensionless parameter N < 0.1, we use Formula (21) to calculate h c and determine h NcB , depending on the following formula: and select the larger of h c and h NcB as h f ,boil . The constant F is determined by the following relationship:

Iterative Calculation of Gas Mass Fraction
One of the main variables of forced convection boiling is the mass fraction α of steam, or vaporization amount, which is used to calculate the latent heat increment, Q L , of the fluid between the two positions in the tube: where the subscript L represents latent heat; H 1 and H 2 are the enthalpies at positions 1 and 2, respectively; and α 1 and α 2 are the mass fractions of vaporization at positions 1 and 2, respectively. The steam mass of each unit section can be calculated using Equation (27). Assuming that the water in a unit in the tube reaches or approaches the saturation temperature, the mass of steam at the inlet of unit j is zero. After the saturated waterbody further absorbs the heat of the tubes, steam is generated at the outlet of unit j. The mass of steam flowing from unit j is equal to that flowing into j+1 unit, as shown in Figure 13. For any straight tube unit, the steam mass fraction is: For units at elbows, the steam mass fraction is: where Q c,n+1 is the heat released at unit n+1 and the subscript c represents condensation.

Calculation of Theoretical Efficiency of a Single Tube
When using the iterative method to calculate the temperature of the endothermic coating and the heat absorbed by the tubes, we assume that the inlet water temperature of the first unit is T f i1 , the temperature of the quartz outer tube is T g1 , and the temperature of the endothermic coating is iteratively calculated as T abs1 . Then, we determine the value of the total heat transfer coefficient U, and then calculate the net heat gain Q f net1 and the outlet water temperature T f o1 of this unit. The outlet water temperature T f o1 of the first unit is equal to the inlet water temperature T f i2 of the second unit, and the process is repeated continuously until the inlet and outlet water temperatures of the nth unit are obtained. The heat gain at the elbow unit n + 1 is negative, satisfying T f o,n = T f i,n+1 . The average water temperature in the tubes is approximately taken as the temperature of the endothermic coating, and the total heat loss at unit n + 1 is obtained as: where U n+1 is the total heat loss coefficient at unit n + 1 and T f m,n+1 is the average water temperature at unit n + 1. The theoretical efficiency η of a single solar vacuum tube system can be determined as follows: where T f o and T f i are the outlet water temperature and inlet water temperature of the entire tube, respectively.

Analysis of the Final Temperature of the Tube
The tube outlet is connected to the evaporator inlet and continuously provides an external heat source for the evaporation system. The inner wall is insulated from the outside. The salt water enters the evaporator after fully absorbing heat in the pipe and gas-liquid separation occurs. The water vapor is continuously pumped away by the exhaust fan at the top of the evaporator, and the salinity of the salty waterbody increases after the evaporation of water, remaining at the bottom of the evaporation tank. In this system, high-temperature salt water enters the evaporator at different flow rates, and a variable-speed axial fan with a maximum capacity of 273 m 3 /h is used to extract the water vapor generated in the evaporation tank at different steady-state flow rates.
In order to measure the temperature and humidity, a set of monitoring instruments with two humidity sensors and three temperature sensors with a humidity error range of ±3% and a temperature error range of ±0.3 K, respectively, are placed in the evaporator. The sensors monitor the temperature of the high-saline waterbody remaining at the bottom of the evaporator, the humidity and temperature of the gas-liquid separation site at the water inlet, and the gas temperature at the top of the evaporator.

Numerical Simulation of Two-Phase Flow in the Tube
A two-phase flow model was used to numerically simulate the saltwater warming effect after the solar tube absorbs sunlight and heat. We selected the following parameters as input conditions: solar radiation intensity I = 1000 W/m 2 , water temperature at the entrance of the drying tube T f i = 30°C, ambient temperature T a = 38°C, wind speed V air = 4 m/s, and flow rates of  Table 1.

Experimental Verification
In order to verify the accuracy of the established theoretical model of two-phase flow, a flow pump was added in the experiment, and three sets of flow rate conditions were selected under the control of the pump. Three flow rates of 8.8, 4.4, and 2.6 kg/h were selected for the experiment, the inlet water temperature of the tube was 30 • C, and the effective solar radiation monitored on site was 1000 W/m 2 . There were 11 temperature probes along the drying tube, the radiation intensity and initial water temperature of these three working conditions were consistent with the numerical values selected in the simulation, and the variable between each working condition was only the flow rate. Setting the heat transfer in the tube to a constant state based on the heat transfer equilibrium equation, we iteratively solved for the temperature of the water, and numerically simulated the forced convection boiling after the water in the tube reached the saturation temperature and compared it with the experimental data, as shown in Table 2. When the flow rate is 8.8 kg/h, the theoretical efficiency of single-phase flow deviates significantly from the experimental efficiency as the temperature of the heat-absorbing coating increases. When(T abs − T a )/I = 0.115, the deviation can reach up to 12%. The reason is that the single-phase flow model considers the complete vacuum (P = 10 −5 mb), and the convective heat transfer coefficient was ignored in the calculation. However, it is difficult to achieve a complete vacuum during the experiment. The two-phase flow model can consider partial vacuum, which can avoid the deviation between the single-phase flow model and the measured value. At this time, the theoretical efficiency and the experimental efficiency curve are basically in line with the error within 2%. When the flow rate is 4.4 kg/h, the waterbody reaches the boiling point near the outlet, and a smaller amount of steam begins to be produced in the tube. The deviation between the single-phase flow model and the experimental model becomes larger as the gas mass fraction increases. When the flow rate is 2.6 kg/h, the waterbody reaches the boiling point at 1500 m, and vapor begins to be continuously generated in the tube. The heat collection efficiency curve calculated by the two-phase flow model agrees better with the experiment, indicating that the consideration of partial vacuum is also applicable to the gas-liquid two-phase flow state. When the water reaches a boiling state, the temperature of the heat-absorbing coating is at a high temperature, which causes the heat loss to increase, so that the influence of the vacuum factor on the heat collection efficiency of the drying tube decreases in the proportion of the total influencing factors, and the accuracy of the single-phase flow model is somewhat improved.
The heat collection efficiency of the tube is calculated according to the temperature measured under different flow rates. The comparison between theoretical efficiency and actual efficiency of each working condition is shown in Figure 14.

Analysis of Heat Collecting Efficiency of the Tubes
The outlet temperature of the tubes is mainly related to the amount of solar radiation and the flow rate of the water in the tube. When the amount of solar radiation is constant, the flow of salty water determines the temperature at the outlet. In order to analyze the relationship between the water flow rate and the outlet temperature of the tube, we adjusted the different flow conditions while keeping the solar radiation intensity and initial temperature almost the same as the simulation conditions, and recorded the outlet temperature during steady-state heat transfer under each condition, as shown in Figure 15. When the salty water temperature T f changed with the flow rate m, the correlation coefficient of the trend line reached 0.99 or higher, and the water temperature can be directly calculated by the power function: Figure 15. Curve of outlet temperature and mass flow rate.
According to the numerical changes obtained by simulation, the relationship between heat collecting efficiency, mass flow rate, and water temperature is established, as shown in Figure 16. When the maximum water temperature is less than 50 • C, the actual heat transfer efficiency of the tube is low, and called a low efficiency zone. The latent heat of vaporization is generated when the water temperature exceeds 100 • C. The theoretical efficiency of the tube is no longer based solely on the temperature change of the waterbody. Therefore, when the maximum water temperature is higher than 100 • C, it is called an overheating zone.
When the flow rate is reduced from 18 to 4.4 kg/h, the heat collecting efficiency of the liquid flow drops by 6.7%. When the flow rate is lower than 4.4 kg/h, the slower the salt water velocity is, the easier it is to reach the boiling condition, but the heat loss is also greater. In the two-phase flow zone, as the share of the gas phase continues to increase, local deceleration, backflow, and stagnation will occur, so that the water temperature no longer rises but remains at the level of 120 • C; the heat collecting efficiency is reduced to 0.79 nearby, and slightly fluctuates. As the flow continues to decrease, the liquid phase can almost completely vaporize before the water flows out of the tube. In order to prevent accidents due to the excessive share of the gas phase in the tube, it is necessary to set up a pressure relief valve. In addition, after the two-phase flow leaves the tube, the gas will take away the stored latent heat, thereby reducing the heat collecting efficiency. To sum up, with the change of flow, the maximum temperature of the flowing water and the heat collecting efficiency of the tube show the opposite trend. That is, the heat collecting efficiency decreases with the decrease of the flow rate, and the water temperature increases with the decrease of the flow rate. The theoretical analysis of two-phase flow can reasonably describe this phenomenon.

Theoretical Analysis
There is pressurized fluid in the tube, which has a certain gas phase mass fraction. After the salty water enters the evaporation system, the space suddenly increases and the exhaust fan at the top continuously draws air, causing the pressure in the waterbody environment to fall, the boiling point of the waterbody to decrease, and a part of the waterbody to become overheated. As the saturation temperature at the overheated state is lower than the original boiling point, the evaporation rate of the waterbody is accelerated (or even boiling) and the vaporization needs to absorb a large amount of latent heat. As the external heat source of the evaporation system, the tube directly determines the efficiency of gas-liquid separation in the evaporation system. The higher the inlet water temperature, the more obvious the gas-liquid separation effect. The heat transfer in the tube is set to a constant state, such that the temperature change of the static waterbody with time is transformed into the temperature change of the moving waterbody with space, which is conducive to observing the phenomena of the waterbody starting to vaporize and boil, heat accumulation, heat transfer, and heat loss. Thermodynamic analysis of the phenomenon of heat dissipation from the endothermic coating at the center to the two ends gives the relationship between the heat collection efficiency factor F c of the tube and the tube length L and the total heat transfer coefficient U of the endothermic coating to the environment. The heat transfer mechanism of the waterbody in the tube is obtained theoretically, and a two-phase flow model on the water temperature and gas phase mass fraction in the tube-that is, the heat and mass transfer balance equations-is used to analyze the critical conditions of water vaporization and boiling, as well as to calculate the thermal efficiency of the tube. Based on the heat transfer equilibrium equation, the temperature of the endothermic coating and the waterbody in the tube can be iteratively solved. Considering that overpressure occurs during the actual operation of the tube, the boiling point of the waterbody becomes larger as the pressure increases during the numerical simulation of the two-phase flow. When the flow rate is reduced from 8.81 to 4.4 kg/h, the heat collection efficiency in the liquid flow zone is reduced by 6.68%. When the flow rate is lower than 4.4 kg/h, such as 2.62 kg/h, the fluid in the tubes is a gas-liquid two-phase flow, indicating that the slower the flow rate of the salt water, the easier it is to reach the boiling condition; however, the heat loss of the tubes is also greater. In the two-phase flow zone, as the proportion of the gas phase increases, the phenomena of deceleration, recirculation, and stagnation will occur in the tubes, such that the water temperature in the drying pipe will not rise and instead remain at 120 • C. Therefore, the heat collection efficiency is reduced to 0.79 and only slightly fluctuates. As the flow rate continues to decrease, the liquid phase is almost completely vaporized before the water flows out of the tubes. In order to prevent accidents caused by an excessive proportion of gas phase in the tubes, it is necessary to install a pressure relief valve. In addition, after the two-phase flow leaves the tubes, the gas will carry away the stored latent heat, thereby reducing the heat collecting efficiency of the tubes. As a result, we introduce this latent heat into the evaporation system, along with the salty water, and release it to drive the evaporation system for gas and liquid two-phase separation.
In summary, as the flow rate changes, the maximum temperature of the flowing water and the heat collection efficiency of the tubes are opposite; that is, the heat collection efficiency of the tubes decreases as the flow rate decreases and the water temperature increases as the flow rate decreases. The theoretical analysis of two-phase flow can reasonably describe this phenomenon.

Comparison
The water production cost of desalination depends on capital cost and operation: equipment and installation, energy, maintenance, and replacement. The energy cost and capital cost occupy nearly 81% [30] of total water production cost of conventional membrane distillation. By comparing the existing conventional desalination treatment, the technologies differ greatly in terms of energy consumption, cost, efficiency, and output due to their different processes. It can be seen from Table 3 that the energy consumption of membrane treatment is far lower than other methods. At the same time, the high cost of the core component of the reverse osmosis membrane and the high frequency of replacement have led to an increase in the overall water production cost. The solar distillation process is simple but has a large area and low efficiency. The output of multistage flash evaporation is low, and the energy consumption is huge due to the need for an external stable heat source. The humidification and dehumidification equipment are expensive, the process flow is complicated, and the heat loss is large. Combining the characteristics of population distribution and water resources in Xinjiang, the system not only controls the overall installation and operating costs, but also ensures that an appropriate amount of water can be produced without chemical treatment, and the local wastewater is reasonably recycled, thereby appropriately solving the local agricultural water shortage and saline-alkali land management problems. The cost of the system mainly includes equipment, transportation, construction, and maintenance. The total is about 160,000 USD. Suitable for inland and rural areas Recycle waste saline water Improve agricultural environment Low water production Affected by local sandstorms Water production efficiency of the system is unstable

Benefit Analysis
Economic benefits: Southern Xinjiang is located in the southern foothills of the Tianshan Mountains, with an area of about 10% of China. The climate is dry and the annual evaporation is 1877.5-2337.4 mm. Due to the generally high soil salinity in the Tarim Basin, the salt content of the 0-30 cm soil layer is 50-300 mg/kg, reaching a maximum of 600-800 mg/kg at the surface, where 47% of the area is abandoned due to secondary salinization; the actual reserved area is only 18,666.7 km 2 [36]. According to a survey [37], soil salinization reduces Xinjiang's grain output by about 720 million kg/year, accounting for about 8.6% of the total grain output of the whole year, causing economic losses of about 3.5 billion RMB. For conventional desalination equipment [38,39], the site construction is restricted by many geographical conditions, and it covers a large area and has a high investment cost. The vast area of Xinjiang, the uneven population distribution, and the characteristics of a large number of remote rural areas all need to be taken into account [40]. There are many drainage canals in various regions, combined with the advantages of the system's total area of less than 400 m 2 and the low construction cost, it is a good choice to improve local economic development.
Ecological benefits: The survey showed that in the Aksu irrigation area, the current water consumption rate of cotton is 0.65-0.75 kg/m 3 , and the water consumption rate for wheat is 1.2 kg/m 3 . This shows that there is still room for improvement in the water productivity of wheat and cotton in this area, while the economic benefits of cotton are better than those of wheat. There are a lot of cotton fields near the system. Due to the greater water consumption of salt field cultivation, the system can improve the water supply by purifying the waste saline water in the drainage canals without sacrificing the original ecological stability. The water is used to grow local cotton fields, effectively improving the local ecological planting environment.
Energy benefits: The energy utilization efficiency of the system is an important parameter to measure its feasibility and benefit. Through the observation records of the inlet and outlet water temperatures, radiation, gasification, and other parameters of the pipe section, the total heat absorption of the heat-absorbing coating is calculated from the heat collection efficiency factor as F c :Q abs = A abs F c [ατI − U(T abs − T a )] = 172, 251.43 KJ, compared with the actual heat gain calculated by the temperature change of the inlet and outlet and the latent heat of vaporization: Q net = Q f + Q f g = 145, 896.97 KJ, where Q f is the sensible heat and Q f g is the latent heat of vaporization. The energy utilization efficiency η can, then, be obtained as: η = Q net Q abs ×100% = 84.7%%, which reflects the feasibility of the energy efficiency of the tubes.
In the future, mainly the heat loss of the tube will be studied. The heat loss at the elbow is relatively small (compared to the straight pipe). The experimental and numerical simulation results show that the accuracy is high enough when only analyzing the heat loss of the straight. However, if the heat loss of the elbow can be considered in the research, the accuracy will be further improved.

Conclusions
Conclusion 1: Considering the problems of strong sunshine, dry climate, and shortage of water resources in Xinjiang, this paper combines solar evaporator and flash evaporation technologies to propose a method for the separation of brackish water based on concentrated heat. Taking into account the characteristics of high temperature resistance and corrosion resistance of quartz tubes and evaporators, sufficient daily routines ensure the stable operation of the system, maximize the utilization rate of the system, and minimize the loss rate.
Conclusion 2: In the study of the system, through the monitoring, numerical simulation, and iterative calculation of the tube, the phase transition formula of the gas-liquid two-phase was obtained as ρ v c v dT v dt = ρ f duH dx . Through thermodynamic analysis of the phenomenon of heat dissipation of the endothermic coating at the center along the pipe section towards both ends, the relationship between the heat collection efficiency factor F c in the tube and the length L of the tube was obtained as F c = 1 − tan mL mL . The total heat transfer coefficient U of the coating to the environment could then be obtained, which revealed the heat transfer mechanism of water. This provides a theoretical basis for ensuring the stable and effective operation of the system.