Integration of Process Modeling , Design , and Optimization with an Experimental Study of a Solar-Driven Humidification and Dehumidification Desalination System

Solar energy is becoming a promising source of heat and power for electrical generation and desalination plants. In this work, an integrated study of modeling, optimization, and experimental work is undertaken for a parabolic trough concentrator combined with a humidification and dehumidification desalination unit. The objective is to study the design performance and economic feasibility of a solar-driven desalination system. The design involves the circulation of a closed loop of synthetic blend motor oil in the concentrators and the desalination unit heat input section. The air circulation in the humidification and dehumidification unit operates in a closed loop, where the circulating water runs during the daytime and requires only makeup feed water to maintain the humidifier water level. Energy losses are reduced by minimizing the waste of treated streams. The process is environmentally friendly, since no significant chemical treatment is required. Design, construction, and operation are performed, and the system is analyzed at different circulating oil and air flow rates to obtain the optimum operating conditions. A case study in Saudi Arabia is carried out. The study reveals unit capability of producing 24.31 kg/day at a circulating air rate of 0.0631 kg/s and oil circulation rate of 0.0983 kg/s. The tradeoff between productivity, gain output ratio, and production cost revealed a unit cost of 12.54 US$/m3. The impact of the circulating water temperature has been tracked and shown to positively influence the process productivity. At a high productivity rate, the humidifier efficiency was found to be 69.1%, and the thermal efficiency was determined to be 82.94%. The efficiency of the parabolic trough collectors improved with the closed loop oil circulation, and the highest performance was achieved from noon until 14:00 p.m.


Introduction
The vast majority of the water on Earth is salt water and surrounds available lands.The natural freshwater resources are limited, especially over a stressed water area in the Arabian Peninsula, where renewable energy sources, such as solar energy, are in the highest range.The most common desalination techniques are of a commercial size and designed for large cities, and small populations in villages experience a lack of freshwater availability [1].Saudi Arabia is located in the Middle East, which has various energy sources, including oil, natural gas, wind and solar radiation [2].The country is considered one of the largest oil producers in the world, and this industry requires further investment for environmental control in order to reduce the emissions, which contribute to the global warming issue.Therefore, renewable energy sources are more favorable research topics for the promotion of investment in clean technologies, and solar energy is a promising field of research in this connection.The average annual horizontal irradiation in the country is 2200 kwh/m 2 /year (251.14W/m 2 ) [2,3], which is in the high range compared to different locations around the world.Saudi Arabia's future expectation of solar power demand would be about 30% by the year 2032 [4].The average annual horizontal irradiance in Saudi Arabia is shown in Figure 1.
Processes 2018, 6, x FOR PEER REVIEW 2 of 35 investment for environmental control in order to reduce the emissions, which contribute to the global warming issue.Therefore, renewable energy sources are more favorable research topics for the promotion of investment in clean technologies, and solar energy is a promising field of research in this connection.The average annual horizontal irradiation in the country is 2200 kwh/m 2 /year (251.14W/m 2 ) [2,3], which is in the high range compared to different locations around the world.Saudi Arabia's future expectation of solar power demand would be about 30% by the year 2032 [4].The average annual horizontal irradiance in Saudi Arabia is shown in Figure 1.Solar energy can be collected using parabolic trough collectors (PTC), which have an efficiency of over 75% [6].The designed system in this work uses a type SAE 10W-30 synthetic blend engine oil, which circulates through the PTCs to transport the thermal energy to the humidifier heat exchanger; hence, oil flow rates are tested for system performance analysis.The PTCs show higher efficiencies when the circulating fluid rate is increased [7] up to an optimum value, at which a further increase can reduce the efficiency.The flow variation inside the absorber tube is performed [8] in order to relate the Reynolds number to the flow pattern, and the results were obtained in March, April, and May, with a laminar flow.In a review study [9], the PTC performance was found to be dependent on the optical efficiency, circulating fluid heat transfer coefficient, heat flux, reflectance of the reflector, the absorptance, the diameter of the absorber and the length of the collector.The PTC has been used in the desalination field due to its superior performance in solar light tracking and focal receiving [10].The parabolic concentrator tracking system is used [11] with the tubular solar still for the purpose of desalination, and the resultant production capacities are 0.28, 0.214 and 1.66 L/day, with a production cost of 0.033, 0.044 and 0.024 US$/L, respectively.Moreover, variable pressure humidification and dehumidification (HDH) is used in solar heat in different configurations [12], and the resultant operating costs were 0.034 and 0.041US$/L, respectively.The PTC is usually installed Solar energy can be collected using parabolic trough collectors (PTC), which have an efficiency of over 75% [6].The designed system in this work uses a type SAE 10W-30 synthetic blend engine oil, which circulates through the PTCs to transport the thermal energy to the humidifier heat exchanger; hence, oil flow rates are tested for system performance analysis.The PTCs show higher efficiencies when the circulating fluid rate is increased [7] up to an optimum value, at which a further increase can reduce the efficiency.The flow variation inside the absorber tube is performed [8] in order to relate the Reynolds number to the flow pattern, and the results were obtained in March, April, and May, with a laminar flow.In a review study [9], the PTC performance was found to be dependent on the optical efficiency, circulating fluid heat transfer coefficient, heat flux, reflectance of the reflector, the absorptance, the diameter of the absorber and the length of the collector.The PTC has been used in the desalination field due to its superior performance in solar light tracking and focal receiving [10].The parabolic concentrator tracking system is used [11] with the tubular solar still for the purpose of desalination, and the resultant production capacities are 0.28, 0.214 and 1.66 L/day, with a production cost of 0.033, 0.044 and 0.024 US$/L, respectively.Moreover, variable pressure humidification and dehumidification (HDH) is used in solar heat in different configurations [12], and the resultant operating costs were 0.034 and 0.041 US$/L, respectively.The PTC is usually installed where the tracking system is accommodated to enhance the performance of the collector.In a study of PTC, with and without a tracking system coupled with a double slope solar still, the tracking PTC is shown to have a higher still temperature and productivity compared with a non-tracked PTC and conventional solar still [13].The solar PTC performance varies with different latitude angles in different seasons; hence, the desalination performance shows an increase during the day [14].A solar driven desalination process can be implemented in arid areas, where large commercial technologies are not economically feasible for a small population; hence, the humidification and dehumidification process can be implemented for this purpose.The sophisticated commercial plants are inadequate in arid areas due to the infrastructure requirements, such as sophisticated switchgears, a power grid, fue transportation and pipelines for small production quantities.Moreover, they require a large number of employees, high technical level and maintenance programs; therefore, the solar HDH becomes a promising process for providing freshwater in arid areas with fewer complications associated with serving a small population village.The HDH process operates at lower temperatures and pressures compared to other desalination techniques, which is advantageous in solar desalination practice.Various arrangements of HDH processes can be implemented, and several studies have been performed in this field in order to improve the process outcome and economics [15].The air flow rate and temperature are crucial factors in the process, where heating the air supply can increase the specific humidity by 1.22 times compared to the ambient air [16].The gained output ratio (GOR) increases with an increasing air to water ratio.This increase can be achieved by changing either air or circulating water flow rates, and it may be effective up to an optimal value, at which a further increase in the air to water ratio can reduce the system GOR [17].Hamed et al. [18] concluded that the solar HDH desalination was strongly affected by ambient conditions and, therefore, that the system temperature rise positively affected the productivity from noon until 5 p.m., with the maximum productivity of 22 L/day and production cost of 0.0578 US$/L.Other parameters, such as the cooling water flow and temperature can alter the system productivity.Ahmed et al. [19] found that the HDH inlet cooling water temperature increases the unit productivity by 10 to 15 L/h, with a production cost of 0.01 US$/L.The solar HDH desalination could improve by understanding the effects of configuration and flow circulation rates on the optimum operating conditions.Combining the HDH desalination unit with a renewable-energy source is the target for sustainability and a clean development process.Solar desalination using collectors is still expensive compared to fossil fuels and needs an improvement in order to convert the solar energy to electrical or thermal energy [20].In order to investigate its process efficiency and economic viability, the process should be analyzed separately before hybridization with another energy source, such as the available fossil fuel in Saudi Arabia [21].Based on a review study by Kabeel et al. [15], further HDH design simulation, to understand the effects of air and water flow rates and optimal conditions, was recommended.The use of series PTCs, using the motor oil as an alternative to thermal fluids, and coupled to closed air with the recirculated water HDH system, is a promising design for the research field.This can enhance the desalination options for decentralized areas.
This project was implemented at King Abdul-Aziz University at Jeddah city to study the effect of various factors in the solar-driven HDH desalination unit using six PTCs in a series configuration.The system comprises multiple circulation loops in the solar heating cycle and the desalination air and water cycles.This study aimed to extend knowledge about solar energy utilization in the field of clean environment technologies and provide headway to overcome the freshwater demand in areas in which natural water resources are scarce.

Experimental Setup and Procedure
The experiments were performed at the King Abdul-Aziz university campus in an open area to allow for the required incidence from the sun at sunrise and sunset in order to achieve a better project analysis and evaluation of its performance.The location latitude is 21 degrees, 30 min north, and the longitude is 39 degrees, 15 min east.

Experimental Setup
The system is designed to utilize solar energy through multiple circulation procedures.The heat input section circulates the oil through the system of the parabolic concentrator, humidifier, bottom heat exchanger, and the oil expansion tank.The water circulates inside the humidifier in order to maintain maximum water temperature and prevent heat loss, which usually occurs in the open water cycle.On the humidifier-packing surface, the water vapor as well as the latent and sensible heat are transferred to the flowing air using a forced circulation technique, which helps to prevent the loss of air and to increase its temperature and humidity.For such a system, the circulating air mostly operates near the saturation curve and close to the equilibrium condition of air to water, therefore improving the humidifier efficiency.The experimental setup for testing the solar humidification and dehumidification desalination system is illustrated in Figure 2. The parabolic concentrators are used to increase the circulating oil temperature, in which the solar tracking system is used to maximize the beam radiation utilization by increasing the useful heat gain of the concentrators.The oil is driven by a gear oil pump, which provides a proper circulation of oil through the PTCs.Therefore, it delivers the gained energy to the heat exchanger on the bottom of the humidifier to circulate the water heating and allows for the oil thermal expansion in the oil tank.The oil circulation procedure allows the system to build up energy instead of using open circuit designs.The quantity of circulating oil in the tank would provide both stability and possible heat storage, since the circulating oil preserves the heat and requires time to complete the circulation cycle.In order to increase the humidifier inlet water temperature and to deliver the necessary sensible heat to the air and water interface, the water circulates through the packing of the humidifier.The water circulation procedure ensures both the heating and humidifying of air, gaining its energy from the heat exchanger on the bottom of the humidifier.The temperature of the circulating air is increased, and the air contains water vapor to deliver to the dehumidifier.In the dehumidifier, the cooling and condensing of the air and water vapor mixture occurs, and part of the water is extracted from the process as a distillate.

Experimental Setup
The system is designed to utilize solar energy through multiple circulation procedures.The heat input section circulates the oil through the system of the parabolic concentrator, humidifier, bottom heat exchanger, and the oil expansion tank.The water circulates inside the humidifier in order to maintain maximum water temperature and prevent heat loss, which usually occurs in the open water cycle.On the humidifier-packing surface, the water vapor as well as the latent and sensible heat are transferred to the flowing air using a forced circulation technique, which helps to prevent the loss of air and to increase its temperature and humidity.For such a system, the circulating air mostly operates near the saturation curve and close to the equilibrium condition of air to water, therefore improving the humidifier efficiency.The experimental setup for testing the solar humidification and dehumidification desalination system is illustrated in Figure 2. The parabolic concentrators are used to increase the circulating oil temperature, in which the solar tracking system is used to maximize the beam radiation utilization by increasing the useful heat gain of the concentrators.The oil is driven by a gear oil pump, which provides a proper circulation of oil through the PTCs.Therefore, it delivers the gained energy to the heat exchanger on the bottom of the humidifier to circulate the water heating and allows for the oil thermal expansion in the oil tank.The oil circulation procedure allows the system to build up energy instead of using open circuit designs.The quantity of circulating oil in the tank would provide both stability and possible heat storage, since the circulating oil preserves the heat and requires time to complete the circulation cycle.In order to increase the humidifier inlet water temperature and to deliver the necessary sensible heat to the air and water interface, the water circulates through the packing of the humidifier.The water circulation procedure ensures both the heating and humidifying of air, gaining its energy from the heat exchanger on the bottom of the humidifier.The temperature of the circulating air is increased, and the air contains water vapor to deliver to the dehumidifier.In the dehumidifier, the cooling and condensing of the air and water vapor mixture occurs, and part of the water is extracted from the process as a distillate.In the dehumidifier the air temperature and its humidity decreased to a point higher than the ambient conditions; therefore, the circulation of this stream allows for further process efficiency improvement.The humidifier circulating water is fed by making water to compensate for the amount In the dehumidifier the air temperature and its humidity decreased to a point higher than the ambient conditions; therefore, the circulation of this stream allows for further process efficiency improvement.The humidifier circulating water is fed by making water to compensate for the amount of water vapor being transferred to the dehumidifier and is extracted as a process yield.The fabricated and installed major components of the experiment are shown in Figure 3.

The Parabolic Concentrators
There are six parabolic concentrators aligned in a series.The solar beam radiation reflects from the concentrator on the lower part of the absorber tube.The six parabolic concentrators are made from glass fiber-reinforced polymer (GFRP), with a thickness of 2.3 mm.The GFRP material was selected due to its availability and cost-effectiveness and, moreover, because it is known to have an excellent thermal insulation.The GFRP are lined with a polymer-based mirror type 3M fused to a stainless steel sheet, with a thickness of 0.5 mm.The polymeric mirror allows the concentrator to have a high efficiency at the working incidence.The dimensions of one PTC used in this experiment are shown in Figure 4.The absorber tube is made from a copper material, which has a high thermal conductivity and is lightweight.The tube is coated with black polymeric paint, which allows the receiver tube to absorb the radiation more efficiently.
The design specification of the PTC and the fabrication material are listed in Table 1.The tracking system of the PTC is designed to work from east to west while the PTCs are oriented in the horizontal north-south direction.

The Parabolic Concentrators
There are six parabolic concentrators aligned in a series.The solar beam radiation reflects from the concentrator on the lower part of the absorber tube.The six parabolic concentrators are made from glass fiber-reinforced polymer (GFRP), with a thickness of 2.3 mm.The GFRP material was selected due to its availability and cost-effectiveness and, moreover, because it is known to have an excellent thermal insulation.The GFRP are lined with a polymer-based mirror type 3M fused to a stainless steel sheet, with a thickness of 0.5 mm.The polymeric mirror allows the concentrator to have a high efficiency at the working incidence.The dimensions of one PTC used in this experiment are shown in Figure 4.The absorber tube is made from a copper material, which has a high thermal conductivity and is lightweight.The tube is coated with black polymeric paint, which allows the receiver tube to absorb the radiation more efficiently.
The design specification of the PTC and the fabrication material are listed in Table 1.The tracking system of the PTC is designed to work from east to west while the PTCs are oriented in the horizontal north-south direction.
A thermally-treated protective glass covers the aperture plane, which aids in the reduction of the convective thermal loss due to winds flowing around the absorber tube.The solar beam radiation, measured as normal, on the PTC surface during the experiment is recorded by a precise solar power meter.The meter is an instrument, used in the solar research field, of type SM206, with an error range of ±5% of the measured value.A thermally-treated protective glass covers the aperture plane, which aids in the reduction of the convective thermal loss due to winds flowing around the absorber tube.The solar beam radiation, measured as normal, on the PTC surface during the experiment is recorded by a precise solar power meter.The meter is an instrument, used in the solar research field, of type SM206, with an error range of ±5% of the measured value.
Engine oil is used as the working fluid in this experiment due to its durability and lower price compared to commercial thermal fluids.The oil type is SAE 10W-30 synthetic blend engine oil.Engine oil has an advantage over water, since oil can sustain higher temperatures unlike water, which  Engine oil is used as the working fluid in this experiment due to its durability and lower price compared to commercial thermal fluids.The oil type is SAE 10W-30 synthetic blend engine oil.Engine oil has an advantage over water, since oil can sustain higher temperatures unlike water, which evaporates and develops hammering in the pipes, requiring high-quality piping and joint fabrication, and this increases the process capital cost.The thermo-physical properties of the circulating oil can change when its temperature increases during solar heating.The thermo-physical properties, shown in Table 2, are used to proceed with the system calculations.The humidifier tank is made from fiberglass-reinforced plastic (FRP) material that has the following advantages: it is lightweight, corrosion resistant, cost-effective and robust enough to handle the system working pressure and temperature.There are different sizes available for such a tank, which can be selected upon the operating capacities of the desalination unit.Moreover, the humidifier body is painted with a black polymer coating to enhance the sunbeam radiation absorption during the day, when the temperature of the system, designed to promote the circulating water, increases.The FRP material is known to have a low thermal conductivity; hence, the combination of both FRB and black paint will enhance the energy conservation in the humidifier.The opening on the top of the humidifier is equipped with a forced draft fan, driven by a 220 v 3-phase motor, to circulate the humid air through the system.A schematic diagram of the humidifier components, dimensions, and flow openings is shown in Figure 5. Several ports are fabricated for the returned air flow from the dehumidifier, makeup water line, circulating water inlet and oil circulation through the bottom heater.The humidifier is internally packed with canvas sheets to promote air-to-water contact by increasing the surface area due to its high porosity.This material is selected for its common availability, low cost and easy replacement, which reduces both capital and running costs.The packing material is supported in the middle of the humidifier tank using a transparent acrylic sheet manufactured from a synthetic polymer.The acrylic sheet is also used to seal and prevent air from passing around the packing in a shortcut, away from the packing.It is installed using adhesives around the middle edge of the humidifier, where the packing is placed.The circulating water is introduced into the humidifier through two top nozzles to ensure a conical shape of flow covering the packing surface area.Both nozzles are capable of delivering a circulating water pump capacity of 0.571 kg/s.The circulating air and water are flowing in a counter flow pattern inside the humidifier.The circulating water is a forced type, using a single-phase constant speed pump driven by a 220 v motor.The bottom heater of the humidifier is an unmixed counter flow heat exchanger, typically used in the internal combustion engine cooling system of automobiles.The heat exchanger manufacturing material is aluminum, with a total surface area of 11.074 m 2 .The heating oil flows inside the tubes are crosswise to the circulating water direction, where the water is maintained above the heat exchanger using a float valve in the makeup water inlet port.
The dehumidifier is a vertically mounted shell and tube heat exchanger, with one pass and no baffles.The humid air stream is introduced into the dehumidifier from the top and flows downward to the bottom of the heat exchanger, and it is then released from the outlet port and recycled for the humidifier.A schematic diagram of the dehumidifier shell and tube heat exchanger is shown in Figure 6.The design of the implemented heat exchanger uses 50 copper tubes, each with an outer diameter of 28.5 mm.The height of the copper tubes is 1300 mm, which together provide a 5.817 m 2 contact surface area between the humid air and the cooling water.There are two water boxes for the heat exchanger, which are made from galvanized steel, with a height of 50 mm and 470 mm in diameter, providing 0.0174 m 3 of cooling water for both water boxes.The cooling water flows from the bottom water box upward through the copper tubes and is then released from the top outlet pipe to ensure that the tubes are uniformly filled, even at lower water cooling flow rates.The humid air flows counter-current to the cooling water stream, in which part of the entrained water vapor condenses at the outer surface of the copper tubes and is then diverted to the distillate water port for collection.The dehumidifier is a vertically mounted shell and tube heat exchanger, with one pass and no baffles.The humid air stream is introduced into the dehumidifier from the top and flows downward to the bottom of the heat exchanger, and it is then released from the outlet port and recycled for the humidifier.A schematic diagram of the dehumidifier shell and tube heat exchanger is shown in Figure 6.The design of the implemented heat exchanger uses 50 copper tubes, each with an outer diameter of 28.5 mm.The height of the copper tubes is 1300 mm, which together provide a 5.817 m 2 contact surface area between the humid air and the cooling water.There are two water boxes for the heat exchanger, which are made from galvanized steel, with a height of 50 mm and 470 mm in diameter, providing 0.0174 m 3 of cooling water for both water boxes.The cooling water flows from the bottom water box upward through the copper tubes and is then released from the top outlet pipe to ensure that the tubes are uniformly filled, even at lower water cooling flow rates.The humid air flows counter-current to the cooling water stream, in which part of the entrained water vapor condenses at the outer surface of the copper tubes and is then diverted to the distillate water port for collection.

Measuring and Controls
Both the desalination unit and the solar tracking system of the PTCs are controlled and monitored using Laboratory Virtual Instrument Engineering Workbench (LabVIEW) software.For solar system tracking, the designed graphical interface using a visual programming language is enabled to control the tracking DC motor motion through the amplified signal, sent by a Galil motion control card.A control loop with a feedback signal allows the motion control system to continuously track throughout the day, as shown in Figure 7.

Measuring and Controls
Both the desalination unit and the solar tracking system of the PTCs are controlled and monitored using Laboratory Virtual Instrument Engineering Workbench (LabVIEW) software.For solar system tracking, the designed graphical interface using a visual programming language is enabled to control the tracking DC motor motion through the amplified signal, sent by a Galil motion control card.A control loop with a feedback signal allows the motion control system to continuously track throughout the day, as shown in Figure 7.The east and west rotations of PTCs are protected by limit switches, installed in both directions for controlling the motor operation at sunset, to return the PTCs to the home position and stop the motor when the home position is reached as well as prepare the collector for the sunrise position of the next day.An absolute type encoder, mounted on the DC motor shaft, controls the angular position of the PTCs to monitor the motion activity.The desalination unit drives and the circulating oil pump drive speed are controlled by an AC frequency inverter, which provides a wide range of speed control; hence, the flow capacities are controlled based on the desired testing value.The circulating water pump of the humidifier is kept at a constant speed, while the air circulating fan speed is altered.The system oil/air/water temperatures are measured using K-type thermocouples, which have a wide The east and west rotations of PTCs are protected by limit switches, installed in both directions for controlling the motor operation at sunset, to return the PTCs to the home position and stop the motor when the home position is reached as well as prepare the collector for the sunrise position of the next day.An absolute type encoder, mounted on the DC motor shaft, controls the angular position of the PTCs to monitor the motion activity.The desalination unit drives and the circulating oil pump drive speed are controlled by an AC frequency inverter, which provides a wide range of speed control; hence, the flow capacities are controlled based on the desired testing value.The circulating water pump of the humidifier is kept at a constant speed, while the air circulating fan speed is altered.The system oil/air/water temperatures are measured using K-type thermocouples, which have a wide range of sensing temperatures, with an accepted accuracy at error limits of ±0.75% [22].The thermocouples are connected to the computer through a field point thermocouple input model.Daytime relative humidity, wind speed, and ambient temperatures are recorded.The circulating air conditions are recorded, and its temperature and humidity are measured by a humidity sensor mounted on the air ducts.The product water from the desalination unit is collected from the distillate port under the dehumidifier and measured in a scaled water storage tank.The measured experimental data uncertainty analysis is described in Appendix A, however, the maximum standard uncertainty is found to be 0.94.

Experimental Procedure
Each experimental operation and recording is performed on a daily basis.The oil expansion tank level is checked in an adequate position.The oil, air, humidifier water, dehumidifier cooling water circulation are started at 06:20 a.m.daily to ensure the system functionality.The circulation provides a uniform temperature over the cycles and complete humidifier packing constitutes a wet condition.The desired oil and air circulating flow rates are established using the AC frequency inverters.The PTCs are checked in the home position, preparing for the sunrise and tracking the solar motion throughout the day.The computer graphical interface of LabVIEW is initiated to start the tracking system using a Galil motion control card.Similarly, the system temperatures and humidity recordings are stored in an excel sheet for the whole experiment run.The solar beam radiation is measured as normal on the horizontal level of the PTCs and is recorded every hour.The distillate obtained from the desalination unit is collected and measured hourly.Sun tracking and system recording are continued until the end of the day, when the return limit switch is activated to rotate the PTCs back to the home position.For the next testing conditions, the circulating oil and air flow rates are changed based on the test plan to record several operation modes in order to analyze the designed system performance.

Experiment Analysis
The desalination system can be split into three major control volumes in order to perform analyses and calculations for each one.The control volumes are classified as follows: (1) Solar parabolic concentrator, (2) Humidifier, (3) Dehumidifier.
The direction of the heat flow starts from the solar collector through the humidifier to the dehumidifier, where part of the heat is retained by recirculating the air back to the humidifier, and the insignificant heat loss is dumped into the cooling reservoir.

Solar Irradiance and Tracking Angles
The beam from the sun can generally be considered as the direct beam and diffused beam; hence, the solar irradiance to the parabolic concentrator would be combined beams.Due to the geometrical design of the parabolic concentrator, the direct beam is considered the useful part of the solar irradiance, and it is affected by the incidence angle and called the tilted beam radiation [23], which can be given by the following relation: where I and G bn are the solar irradiance and the normal beam radiation, respectively.The incidence angle θ must be calculated in order to obtain the useful solar irradiance to the parabolic concentrator and can be found by calculating other solar angles.

• Latitude angle (L)
This depends on the location of the site on the Earth and the site of the experiment found above the equator; hence, the latitude of the experiment site is positively valued and found to be 21.5 • .

• Hour angle (h)
Solar noon is defined as that time of day when the sun appears due south (north), and the solar time is the time of day measured from solar noon [24].The reference day hour is 12:00 p.m., where the mid daylight hours are calculated to have the maximum expected radiation.The hour angles are negatively valued for the first half of the day and positively valued in the afternoon.The following relation is employed to calculate the hour angle [17,19]: where t s is the solar time based on the local time and substituted in a 24-h system.

• Declination angle (δ)
The declination angel is influenced by the location of the Earth with respect to the sun.The value of this angle varies positively or negatively depending on the season and specifically on the number of days in the year due to the rotation of the Earth.The following relation gives the declination angle [23,25]: where N is the number of the day in the year.For the experiment site, the maximum positive declination angle for the earth is 23.45 in the summer, and the maximum negative value is −23.45 in the winter.

•
Altitude angel (α) The altitude angle is calculated for the horizontal surface of the experiment site.It depends on the elevation of the sun in the sky; hence, the projection of the sunbeam on the horizontal surface forms this angle.Therefore, the experimental site latitude, hour angle and declination angle of the earth with respect to the sun have a great influence on this angle.The following relation is used to calculate it [23,26]: • Azimuth angle (φ) The azimuth angle is concerned with the south direction of the experimental site.It is the angle between the south direction of the experimental site and the sunbeam projection on the horizontal site and can be calculated as follows [23]: • Surface azimuth angle (φ s ) The surface azimuth angle is the deviation of the projection, on a normal horizontal plane on the surface, from the local meridian [27].For the present case, where the tracking system operates from east to west, this angle will be either −90 • , when the concentrator is facing east, or 90 • , when it is facing west [23].

•
Zenith angle (θ z ) The magnitude of the zenith angle is obtained from the difference between the altitude angles and the vertical position on the surface, which depends on the sun's position in the sky.The zenith angle is calculated by the following relation [28]: • The incident angle (θ) This is the sunbeam projection angle and is measured from the vertical axis to the surface.It is obtained by Kreith and Kreider, and Duffie and Beckman [23,27,29,30], and is shown as follows:

Concentrator Thermal Analysis
The optical efficiency of the solar parabolic concentrator is required in order to proceed with the calculation of the thermal efficiency of the system.It represents the percentage of the absorbed solar irradiance "energy" that falls on the absorber tube.This efficiency accounts for the system's random and non-random errors.The equation used to calculate the optical efficiency is given by [23,27]: where A f is the geometric ratio between the aperture lost area A l of the concentrator and the total aperture area A a and can be obtained by the following relation: The geometric ratio of the design parabolic concentrator is calculated and found to be A f = 0.511.The intercept factor γ for the concentrator is obtained by Abdel-Hady F. et al. [31] and found to be 0.84, when the absorber tube outer diameter is 0.42 mm, and still 16.3% of rays dissipated away from the absorber tube.Therefore, the absorber tube diameter is enlarged to 0.54 mm to allow for maximum sunbeam coverage; hence, the intercept factor becomes 1, and the collector concentration ratio C = W/πD o is decreased to 9.984.The lost area, due to the end effect, is not an effective area and is significantly affected by the solar incidence angle θ, which can be calculated by the following relation [23]: The concentrator instantaneous efficiency is known to be the ratio of the amount of heat gained to the amount of solar irradiance that falls on the absorber tube within its effective surface area.The following relation is used to calculate this efficiency [23,27,28]: The heat removal factor F R , in Equation ( 13) is employed to calculate the collector instantaneous efficiency, which is based on the concentrator fluid inlet temperature T i and ambient temperature [28,30].
From the experimental data, obtained daily at different operating conditions, the required variables are measured and analyzed.ASHRAE Standard 93-2003, "Methods of testing to determine the thermal performance of solar collectors" (ASHRAE 2003), is used.The concentrator instantaneous efficiency η is calculated from Equation (12), where the term (T i − T a )/I in Equation ( 13) is measured experimentally.The plot of the instantaneous efficiency against (T i − T a )/I reveals the value of F R η o , which represents the intercept of the line, and (F R U L /C) is the slope of this line; hence, the values of F R and U L are obtained.The overall heat transfer coefficient of the absorber tube can be calculated from the following relation [23,28]: where the convective heat transfer coefficient h ab inside the absorbers tube is calculated from the following relation [32]: Nu = 0.023Re 0.8 Pr 0.4 (14) For the fully developed laminar flow pattern, the Nu number value is 4.36 [32].The Re and Pr numbers are obtained from the following relations, using the fluid properties inside the absorber tube.

Humidifier Bottom Heater Section
A cross flow type heat exchanger, in which the oil flows inside the tubes in a cross direction to the outlet water from the humidifier packing, which flows downward through the heat exchanger, is incorporated in the humidifier.The heat exchanger design parameters and the experimental data are used to calculate the heat exchanger overall heat transfer coefficient in the following manner.The effectiveness relation of the cross-flow heat exchanger, with a single pass and both fluids unmixed, are given below [32,33]: where C r is the heat capacity ratio, which is given as C r = C min /C max .The heat capacities are calculated for the cold stream C w = m w c w , and the hot fluid C oil = m oil c oil is then categorized based on their values as maximum or minimum.The value of the heat exchanger effectiveness is calculated from the experiment heat analysis as follows: From the previously obtained parameters, Equation ( 18) is solved for the NTU value at the operating condition.Therefore, NTU is used to calculate the overall heat transfer coefficient of the humidifier bottom heat exchanger as shown below.
The mass and energy balance across the heater is obtained, based on the configuration of the streams, as shown in Figure 7.
Mass balance: Energy balance: where

Humidifier Packing Section
According to the designed humidification and dehumidification desalination unit, shown in Figure 8, the heat and mass balance equations are obtained from the control volume 1, as follows: Mass balance across the packing: Energy balance for air and water streams:

Humidifier Packing Section
According to the designed humidification and dehumidification desalination unit, shown in Figure 8, the heat and mass balance equations are obtained from the control volume 1, as follows: Mass balance across the packing: Energy balance for air and water streams: The mass transfer coefficient of the humidifier is obtained by using the log mean enthalpy difference (LMED) method.This method showed a good consistency with ε-NTU [34,35]; hence, the formulation of the balance is obtained as follows: where ΔHm is the log mean enthalpy difference: where f is the developed analytically-based correction factor, which is calculated from:

Dehumidifier
The flow streams of the dehumidifier (CV2) are shown in Figure 9.The heat and mass balance of the dehumidifier are formulated as follows: Mass balance: Energy balance: The mass transfer coefficient of the humidifier is obtained by using the log mean enthalpy difference (LMED) method.This method showed a good consistency with ε-NTU [34,35]; hence, the formulation of the balance is obtained as follows: where ∆H m is the log mean enthalpy difference: where f is the developed analytically-based correction factor, which is calculated from:

Dehumidifier
The flow streams of the dehumidifier (CV2) are shown in Figure 9.The heat and mass balance of the dehumidifier are formulated as follows: Mass balance: Energy balance: where where In order to calculate the overall heat transfer coefficient in the condenser, the log mean enthalpy difference LMED method [32] is used.This method is adequately used in the present case, where the dehumidifier is considered as a shell and tube heat exchanger with a single pass.The applied balance across the condenser is shown below.
where ΔTlm is calculated as:

Humidifier Efficiency
This factor shows the capability of the humidifier to transfer the amount of water vapor, compared to that in the full saturation condition.It indicates the performance of the design, along with the packing, and can be calculated as follows: In order to calculate the overall heat transfer coefficient in the condenser, the log mean enthalpy difference LMED method [32] is used.This method is adequately used in the present case, where the dehumidifier is considered as a shell and tube heat exchanger with a single pass.The applied balance across the condenser is shown below.
where ∆T lm is calculated as: 3.4.System Key Performance Parameters

Humidifier Efficiency
This factor shows the capability of the humidifier to transfer the amount of water vapor, compared to that in the full saturation condition.It indicates the performance of the design, along with the packing, and can be calculated as follows: The energy-based effectiveness of the humidifier is a non-dimensional value, which can evaluate the humidifier performance based on the maximum possible change in the total enthalpy rate.The energy effectiveness is calculated based on the following [36]: 3.4.2.System Gain Output Ratio (GOR) The solar energy is used in this system as the heating source.The yield from the desalination is the amount of distillate produced; hence, the yield over the input energy to the system represents the gain output ratio of the system:

Humidifier Thermal Efficiency
In order to measure the energy conversion in the humidifier, the input energy to the system is related to the output from the humidifier.It indicates how much thermal energy is lost in the humidifier and not being used for the purpose of humidification:

.4. Convective Heat and Mass Transfer Coefficients
The ratio of the convective heat transfer coefficient to the convective mass transfer coefficient represents the amount of heat per unit volume required to raise the temperature by one degree.This ratio is affected by the type of process components, insulation and air/water stream flow variations.The ratio is calculated from the Chilton-Colburn analogy [26]: where Sc = µ a ρ a D water/air (38) The diffusivity of water in the air is calculated by the following relation [37] in the temperature range of 273 K < T < 373 K: D water/air = 1.97 × 10 −5 T a avg 256 (39)

Thermodynamic Relations
The properties of the fluids are calculated using the thermodynamic relations.The following relation calculates the absolute humidity [38]: The partial pressure P v of water vapor is calculated from Antoine equation [39]: The humid air enthalpy is calculated by [38]: The specific heat of humid air c a [39] is obtained from the following relation:

Process Economics
The preliminary economic analysis is based on the experiment cost information obtained during the construction process.The capital investment data, including the installation, are collected and shown in Table 3, which can be used for the production cost analysis.The summation of fixed and variable operation costs gives the cost of the production, where the total operation cost will include the operation, maintenance, and the annual capital cost for the process, which can be obtained by the following relation [40]: The maintenance cost is taken as 3% [41,42] of the annual capital charges, which can be less than 2%, as estimated by El-Dessouky et al. [43].The uniform series of payment method is used to calculate the annual capital cost (C A ).The capital recovery of the project, invested at a 5% interest rate for 20 years, is obtained from the following relation, as an annualized form [44]: The calculations accounted for the variable operating cost, which is influenced by the price of the energy consumed by the process in the drives and control system.The local energy price is 0.048 US$/kWh [45], as the standard price for the industrial sector in Saudi Arabia.The unit is assumed to operate 10 h per day, with the availability of 90% of the year.The following relation is used to calculate the production cost:

Results and Discussion
The system was operated and tested to evaluate the key performance factors and the productivity at different operating conditions.

System Productivity
The HDH unit productivity is affected by the circulating air to water ratio as well as by the circulating oil flow rate for the designed solar-driven process.The process has been tested at different air flow rates with a constant oil flow over a certain time period; then, the oil flow rate is changed, and the air flow is set to the next desired testing flow rates.System productivity, at different air and oil flow rates, is illustrated in Figure 10.The objective of increasing the air flow rate is to obtain the highest entrained water in the air at the maximum efficiency of the humidifier in order to separate them by condensation in the dehumidifier section.The maximum obtained system productivity is 24.31 kg/day at the air flow of 0.0631 kg/s.This productivity is enhanced by the oil flow of 0.0983 kg/s, which delivered the thermal energy to the humidifier at the maximum balance between the incoming solar irradiance and the air cooling effect in the humidifier.Moreover, the effect of increasing air flow rate on the productivity, which shows a reduction in the process yield when the air flow rate is increased, is presented in Figure 10.
Processes 2018, 6, x FOR PEER REVIEW 18 of 35 highest entrained water in the air at the maximum efficiency of the humidifier in order to separate them by condensation in the dehumidifier section.The maximum obtained system productivity is 24.31 kg/day at the air flow of 0.0631 kg/s.This productivity is enhanced by the oil flow of 0.0983 kg/s, which delivered the thermal energy to the humidifier at the maximum balance between the incoming solar irradiance and the air cooling effect in the humidifier.Moreover, the effect of increasing air flow rate on the productivity, which shows a reduction in the process yield when the air flow rate is increased, is presented in Figure 10.At a circulating oil flow rate of 0.0983 kg/s, the highest yield is obtained.The productivity of the system is reduced by 26% when the air flow increased to 0.1315 kg/s and further decreased by 52.2% and 70.1% at air flow rates of 0.1998 kg/s and 0.2682 kg/s, respectively.The circulating water temperature has a significant role in the humidifier performance, where the system's highest productivity is obtained at a circulating water temperature of 49.23 °C.The circulating water temperature is decreased when the circulating air flow rate increases, and the maximum reduction was 21.7% at the highest air flow rate.The experimental output is shown to have some fluctuations due to many factors that interfere with the process stability, such as the climate conditions and characteristics of the independent components.In general, each tested condition shows a trend in the process behavior despite the fluctuation rate, which indicates the system outcome.The heat input on the bottom heater of the humidifier provokes the productivity of the system by increasing the humidifier circulating water temperature, which is shown in Figure 11.In the humidifier unit, the highest oil temperature is achieved when the unit performance is combined with the oil flow rate, achieving the required energy input for the maximum circulating water temperature.Increasing the oil flow positively affects the process productivity until the point of the tradeoff between either the installation of an enormous humidifier heat exchanger or a larger solar concentrator aperture area.The reduction of heating at higher oil flow rates is attributed to many factors, such as the absorber heat removal coefficient and less resident time of the oil flows inside the absorber tube and humidifier heat exchanger.At a circulating oil flow rate of 0.0983 kg/s, the highest yield is obtained.The productivity of the system is reduced by 26% when the air flow increased to 0.1315 kg/s and further decreased by 52.2% and 70.1% at air flow rates of 0.1998 kg/s and 0.2682 kg/s, respectively.The circulating water temperature has a significant role in the humidifier performance, where the system's highest productivity is obtained at a circulating water temperature of 49.23 • C. The circulating water temperature is decreased when the circulating air flow rate increases, and the maximum reduction was 21.7% at the highest air flow rate.The experimental output is shown to have some fluctuations due to many factors that interfere with the process stability, such as the climate conditions and characteristics of the independent components.In general, each tested condition shows a trend in the process behavior despite the fluctuation rate, which indicates the system outcome.The heat input on the bottom heater of the humidifier provokes the productivity of the system by increasing the humidifier circulating water temperature, which is shown in Figure 11.In the humidifier unit, the highest oil temperature is achieved when the unit performance is combined with the oil flow rate, achieving the required energy input for the maximum circulating water temperature.Increasing the oil flow positively affects the process productivity until the point of the tradeoff between either the installation of an enormous humidifier heat exchanger or a larger solar concentrator aperture area.The reduction of heating at higher oil flow rates is attributed to many factors, such as the absorber heat removal coefficient and less resident time of the oil flows inside the absorber tube and humidifier heat exchanger.

Energy of the System
Based on the oil circulation flow rate, the solar energy gain varies; hence, the input energy to the heat exchanger on the bottom of the humidifier is affected by this variation.The circulating air entering the humidifier has the effect of cooling the system, which has lost the sensible heat in the dehumidifier.The energy input to the heat exchanger on the bottom of the humidifier is illustrated in Figure 12.In the present system, the energy input increased, in most cases, as the oil circulation rate increased up to the desired heating condition, which ranges from 2.39 to 3.07 kWh at an oil flow of 0.0983 kg/s.At a higher humidifier air flow rate, the bottom heat exchanger did not yield the required input energy, and the energy supplied was reduced by 39.3% at the maximum oil flow rate compared to the highest obtained air flow rate of 0.0631 kg/s.

Energy of the System
Based on the oil circulation flow rate, the solar energy gain varies; hence, the input energy to the heat exchanger on the bottom of the humidifier is affected by this variation.The circulating air entering the humidifier has the effect of cooling the system, which has lost the sensible heat in the dehumidifier.The energy input to the heat exchanger on the bottom of the humidifier is illustrated in Figure 12.In the present system, the energy input increased, in most cases, as the oil circulation rate increased up to the desired heating condition, which ranges from 2.39 to 3.07 kWh at an oil flow of 0.0983 kg/s.At a higher humidifier air flow rate, the bottom heat exchanger did not yield the required input energy, and the energy supplied was reduced by 39.3% at the maximum oil flow rate compared to the highest obtained air flow rate of 0.0631 kg/s.

Energy of the System
Based on the oil circulation flow rate, the solar energy gain varies; hence, the input energy to the heat exchanger on the bottom of the humidifier is affected by this variation.The circulating air entering the humidifier has the effect of cooling the system, which has lost the sensible heat in the dehumidifier.The energy input to the heat exchanger on the bottom of the humidifier is illustrated in Figure 12.In the present system, the energy input increased, in most cases, as the oil circulation rate increased up to the desired heating condition, which ranges from 2.39 to 3.07 kWh at an oil flow of 0.0983 kg/s.At a higher humidifier air flow rate, the bottom heat exchanger did not yield the required input energy, and the energy supplied was reduced by 39.3% at the maximum oil flow rate compared to the highest obtained air flow rate of 0.0631 kg/s.

The GOR
At the low oil flow rate, the process revealed no significance in the unit productivity compared to the high oil flow rates.The variation of the averaged GOR versus oil flow rate of the desalination unit, under the tested conditions, is illustrated in Figure 13.The GOR was increasing at a higher oil flow rate due to the increase in the supply of the thermal energy from the solar concentrator.Additionally, the GOR decreases at a higher circulating air flow rate due to the enhanced cooling effect on the process.Moreover, at oil circulations of 0.984 kg/s and lower, the unit GOR and productivity was significantly reduced, and the process output was adversely affected.

The GOR
At the low oil flow rate, the process revealed no significance in the unit productivity compared to the high oil flow rates.The variation of the averaged GOR versus oil flow rate of the desalination unit, under the tested conditions, is illustrated in Figure 13.The GOR was increasing at a higher oil flow rate due to the increase in the supply of the thermal energy from the solar concentrator.Additionally, the GOR decreases at a higher circulating air flow rate due to the enhanced cooling effect on the process.Moreover, at oil circulations of 0.984 kg/s and lower, the unit GOR and productivity was significantly reduced, and the process output was adversely affected.The circulation rates of the oil and air are shown to have higher GORs at an oil flow rate of 0.1202 kg/s, and the process GOR is higher at the lowest circulating air flow rate of 0.0631 kg/s compared with other tested conditions.

Humidifier Efficiency and Thermal Efficiency
The humidifier efficiency is linked to the design, packing type and air-to-water ratio.All these parameters are reflected on the mass transfer coefficient and therefore shown in the performance efficiency of the humidifier.The variation of the humidifier efficiency at different air to water flow ratios and different heat input conditions, represented by oil flow rates, are shown in Figure 14.The amount of vapor transferred is above 60%, compared to the ideal saturated condition, in order to humidify the air.The thermal efficiency of the humidifier, where the performance appears to be lower when the oil flow rate increases to 0.0983 kg/s and above, is shown in Figure 15.The humidifier has the efficiency of 78% and above to deliver the energy from the heat input to the coldest part, which is the humid air in this condition.
The humidifier average energy-based effectiveness for all cases is shown in Figure 16.At higher circulating air flow rates, the energy-based effectiveness is shown to be decreased.Air flow rate has an influence on the cooling of the system and reduces the actual heat transfer rate.Therefore, a further energy supply is required.For all cases, increasing the oil flow rate has an insignificance effect of energy-based effectiveness reduction compared to the increase to the circulating air flow rate.At the highest test air flow of 0.2682 kg/s, the average energy-based effectiveness is shown to be from 0.4 to 0.476.The circulation rates of the oil and air are shown to have higher GORs at an oil flow rate of 0.1202 kg/s, and the process GOR is higher at the lowest circulating air flow rate of 0.0631 kg/s compared with other tested conditions.

Humidifier Efficiency and Thermal Efficiency
The humidifier efficiency is linked to the design, packing type and air-to-water ratio.All these parameters are reflected on the mass transfer coefficient and therefore shown in the performance efficiency of the humidifier.The variation of the humidifier efficiency at different air to water flow ratios and different heat input conditions, represented by oil flow rates, are shown in Figure 14.The amount of vapor transferred is above 60%, compared to the ideal saturated condition, in order to humidify the air.The thermal efficiency of the humidifier, where the performance appears to be lower when the oil flow rate increases to 0.0983 kg/s and above, is shown in Figure 15.The humidifier has the efficiency of 78% and above to deliver the energy from the heat input to the coldest part, which is the humid air in this condition.
The humidifier average energy-based effectiveness for all cases is shown in Figure 16.At higher circulating air flow rates, the energy-based effectiveness is shown to be decreased.Air flow rate has an influence on the cooling of the system and reduces the actual heat transfer rate.Therefore, a further energy supply is required.For all cases, increasing the oil flow rate has an insignificance effect of energy-based effectiveness reduction compared to the increase to the circulating air flow rate.At the highest test air flow of 0.2682 kg/s, the average energy-based effectiveness is shown to be from 0.4 to 0.476.

Enthalpy Difference
The enthalpy of the humidifier outlet is the combination of the enthalpy of the dry air and the enthalpy of the entrapped water vapor within the stream.The outlet enthalpy is measured relative to the inlet air enthalpy in the humidifier.The calculated enthalpies are related to the specific and latent heats at the operating temperature and the amount of vapor in the streams, which have been computed by Equation (42).The enthalpy difference of the humidifier, at different operating air flow rates and different circulating oil flow rates, is illustrated in Figure 17.

Enthalpy Difference
The enthalpy of the humidifier outlet is the combination of the enthalpy of the dry air and the enthalpy of the entrapped water vapor within the stream.The outlet enthalpy is measured relative to the inlet air enthalpy in the humidifier.The calculated enthalpies are related to the specific and latent heats at the operating temperature and the amount of vapor in the streams, which have been computed by Equation (42).The enthalpy difference of the humidifier, at different operating air flow rates and different circulating oil flow rates, is illustrated in Figure 17.At lower air flow rates, the enthalpy difference of the system is shown to be greater, and the oil circulation rate of 0.0983 kg/s is shown to increase the enthalpy difference for most of the air streams.

Humidifier Mass Transfer Coefficient
The average mass transfer coefficient of the experiment is plotted against the oil flow rate and is shown in Figure 18.The effect of both heating and air/water circulation on the coefficient is shown in this figure.The air-water flow rate significantly affects the value of the mass transfer coefficient during all experimental runs.Increasing the air flow rate reveals an increase of the mass transfer coefficient up to the air flow of 0.2682 kg/s, at which point the coefficient declined due to concentration and heat potential complications.The process is known to have a mass transfer due to the concentration difference on the airside of the interface and latent heat potentials.The mass transfer coefficient is calculated from Equation (26), in which the process driving force is expressed in terms of the enthalpy difference, where the enthalpy is a thermodynamic property calculated using Equation (42).Hence, the driving force is the function of both temperature and humidity.The variation of the circulating oil flow changes the heating input to the humidifier, which is reflected in the mass transfer behavior.The air-to-water ratio the humidity in the system and dominantly controls the mass transfer coefficient values, where the temperature influences the system efficiency when further heat is introduced.At air flow rates of 0.0632 and 0.2682 kg/s, the system could not gain the required heating due to an insufficient air flow, low air flow rate or excessive air flow, which led to a cooling effect at higher air flow rates.At lower air flow rates, the enthalpy difference of the system is shown to be greater, and the oil circulation rate of 0.0983 kg/s is shown to increase the enthalpy difference for most of the air streams.

Humidifier Mass Transfer Coefficient
The average mass transfer coefficient of the experiment is plotted against the oil flow rate and is shown in Figure 18.The effect of both heating and air/water circulation on the coefficient is shown in this figure.The air-water flow rate significantly affects the value of the mass transfer coefficient during all experimental runs.Increasing the air flow rate reveals an increase of the mass transfer coefficient up to the air flow of 0.2682 kg/s, at which point the coefficient declined due to concentration and heat potential complications.The process is known to have a mass transfer due to the concentration difference on the airside of the interface and latent heat potentials.The mass transfer coefficient is calculated from Equation (26), in which the process driving force is expressed in terms of the enthalpy difference, where the enthalpy is a thermodynamic property calculated using Equation (42).Hence, the driving force is the function of both temperature and humidity.The variation of the circulating oil flow changes the heating input to the humidifier, which is reflected in the mass transfer behavior.The air-to-water ratio affects the humidity in the system and dominantly controls the mass transfer coefficient values, where the temperature influences the system efficiency when further heat is introduced.At air flow rates of 0.0632 and 0.2682 kg/s, the system could not gain the required heating due to an insufficient air flow, low air flow rate or excessive air flow, which led to a cooling effect at higher air flow rates.

Convective Mass Transfer Coefficient
The convective heat to mass transfer coefficient plot, where the coefficient increases as the air flow rate increases, is shown in Figure 19.At lower air flow rates, this ratio decreases when the oil flow rate reaches 0.0983 kg/s, at which point maximum productivity is obtained.This indicates the reduction in the thermal loss and the usefulness of the heat for enhancing the humidifier-driving force through the mass transfer controlled by the humidity difference.At a lower oil flow rate, the process either loses heat or has an insufficient heat input due to higher air flow rates.The air flow rate at 0.2682 shows a higher ratio trend, which can be further improved by using a higher-grade insulation material to enhance the process productivity at higher air flow rates and avoid heat loss.For the present system, the productivity was shown to have declined when the air increased to this limit; hence, optimization is considered a practical case, which could be implemented in order to obtain a tradeoff between the cost, productivity, and size for a good performance.The best unit performance was obtained at the lowest tested air flow and the point at which the oil flow reached 0.0983 kg/s, where the lower loss was obtained per unit volume.

Convective Mass Transfer Coefficient
The convective heat to mass transfer coefficient plot, where the coefficient increases as the air flow rate increases, is shown in Figure 19.At lower air flow rates, this ratio decreases when the oil flow rate reaches 0.0983 kg/s, at which point maximum productivity is obtained.This indicates the reduction in the thermal loss and the usefulness of the heat for enhancing the humidifier-driving force through the mass transfer controlled by the humidity difference.At a lower oil flow rate, the process either loses heat or has an insufficient heat input due to higher air flow rates.The air flow rate at 0.2682 shows a higher ratio trend, which can be further improved by using a higher-grade insulation material to enhance the process productivity at higher air flow rates and avoid heat loss.For the present system, the productivity was shown to have declined when the air increased to this limit; hence, optimization is considered a practical case, which could be implemented in order to obtain a tradeoff between the cost, productivity, and size for a good performance.The best unit performance was obtained at the lowest tested air flow and the point at which the oil flow reached 0.0983 kg/s, where the lower loss was obtained per unit volume.

Convective Mass Transfer Coefficient
The convective heat to mass transfer coefficient plot, where the coefficient increases as the air flow rate increases, is shown in Figure 19.At lower air flow rates, this ratio decreases when the oil flow rate reaches 0.0983 kg/s, at which point maximum productivity is obtained.This indicates the reduction in the thermal loss and the usefulness of the heat for enhancing the humidifier-driving force through the mass transfer controlled by the humidity difference.At a lower oil flow rate, the process either loses heat or has an insufficient heat input due to higher air flow rates.The air flow rate at 0.2682 shows a higher ratio trend, which can be further improved by using a higher-grade insulation material to enhance the process productivity at higher air flow rates and avoid heat loss.For the present system, the productivity was shown to have declined when the air increased to this limit; hence, optimization is considered a practical case, which could be implemented in order to obtain a tradeoff between the cost, productivity, and size for a good performance.The best unit performance was obtained at the lowest tested air flow and the point at which the oil flow reached 0.0983 kg/s, where the lower loss was obtained per unit volume.

Optimum Condition
The results of the experiment are analyzed in terms of productivity, production cost and GOR in order to find the best operating conditions that meet the process economics with a high yield.The behavior of the process at different circulating oil and air flow rates, obtaining a tradeoff between the maximum productivity, production cost, and the GOR, is illustrated in Figure 20.Under all operating conditions, the process reveals a high productivity and GOR, with the lowest operating cost, at a circulating oil flow of 0.0983 kg/s.The maximum productivity obtained is 24.31 kg/day at a circulating air flow rate of 0.0631 kg/s, and the GOR is found to be 0.51.This condition has the minimum operating cost, where the production cost attained is 12.54 US$/m 3 , which is 0.0125 US$/liter.Since the process used solar energy for heating, the heating energy cost would not affect the unit production cost, although the capital investment would.The process economics can be improved by the reduction of the system capital investment, such as automation and controls, since the system works at low pressures and temperatures, or by increasing the productivity by enhancing the power utilization using an auxiliary heating system operated by solar photovoltaic (PV) panels.The general output of the process reduces as the circulating air flow rate increases.At the highest air flow rate, the economics of the process is found to be discouraging, with a production cost of 38.92 US$/m 3 under the optimal circulating oil condition.Similarly, the GOR of the process has a tendency toward further reduction when the circulating air flow rate is increased.There is also uncertainty in the operating variables due to the variation during different seasons throughout the year, which could affect the overall system results, but the process is expected to behave in a manner similar to the results obtained.The more expensive design may allow for a higher circulating oil temperature of the parabolic concentrator, which requires a larger heat exchanger surface area or humidifier air and water contact area; hence, the benefits could require a higher system capital cost.The maximum output of the experiment, compared with the previous work in the field of solar desalination using the HDH method, is presented in Table 4.Despite the design and capacities of the processes in other research work, the process design in the present study has economic improvements, which enhance the commercialization potential of the process, as demonstrated in Table 4.

Collector Efficiency
The behavior of the averaged efficiency of the PTC under the operating conditions is illustrated in Figure 21.A decrease of the instantaneous efficiency of the PTC, following the increase of the circulating oil flow rate, is shown in this figure .The low efficiencies at high circulating oil flow rates, attributed to the shorter residence time of oil in the absorber tubes, reduce the heat gain of the PTC.

Collector Efficiency
The behavior of the averaged efficiency of the PTC under the operating conditions is illustrated in Figure 21.A decrease of the instantaneous efficiency of the PTC, following the increase of the circulating oil flow rate, is shown in this figure .The low efficiencies at high circulating oil flow rates, attributed to the shorter residence time of oil in the absorber tubes, reduce the heat gain of the PTC.The circulating water temperature lines of the humidifier, shown in Figure 19, illustrate the condition of the humidification process within the operating oil flow, where the lower oil flow revealed a lower productivity and high PTC efficiency.Moreover, the temperatures of the circulating water in the humidifier were found to be higher at the oil flow rate of 0.0983 kg/s than with other settings, which shows the optimum oil flow rate for the heat exchange in the humidifier.The average PTC efficiency at the lowest oil flow rate is 27.8% and is reduced by 49.5% at the maximum flow rate.The circulating water temperature lines of the humidifier, shown in Figure 19, illustrate the condition of the humidification process within the operating oil flow, where the lower oil flow revealed a lower productivity and high PTC efficiency.Moreover, the temperatures of the circulating water in the humidifier were found to be higher at the oil flow rate of 0.0983 kg/s than with other settings, which shows the optimum oil flow rate for the heat exchange in the humidifier.The average PTC efficiency at the lowest oil flow rate is 27.8% and is reduced by 49.5% at the maximum flow rate.In the designed system, the ratio of the heat gained to the amount of beam radiation concentrated in the absorber tube is revealed to be higher, especially from the middle to the end of the solar day due to closed system circulation.Moreover, for all of the cases in which the incidence angle increased throughout the day, the losses due to the end effect decreased; hence, the PTC efficiency increased.When the beam radiation reduced at the end of the day, the system still can deliver the energy to the humidifier at a higher rate due to the stored heat in the oil expansion tank and the circulating air and water of the humidifier.The behavior of the circulating oil inlet and outlet temperatures in the parabolic concentrators is shown in Figure 22.Generally, at higher oil flow rates, the system heat storage capabilities increased due to multiple circulation procedures.This positively affects the value of the parabolic concentrator efficiency; hence, the minimum beam radiation utilization is obtained.In the designed system, the ratio of the heat gained to the amount of beam radiation concentrated in the absorber tube is revealed to be higher, especially from the middle to the end of the solar day due to closed system circulation.Moreover, for all of the cases in which the incidence angle increased throughout the day, the losses due to the end effect decreased; hence, the PTC efficiency increased.When the beam radiation reduced at the end of the day, the system still can deliver the energy to the humidifier at a higher rate due to the stored heat in the oil expansion tank and the circulating air and water of the humidifier.The behavior of the circulating oil inlet and outlet temperatures in the parabolic concentrators is shown in Figure 22.Generally, at higher oil flow rates, the system heat storage capabilities increased due to multiple circulation procedures.This positively affects the value of the parabolic concentrator efficiency; hence, the minimum beam radiation utilization is obtained.The circulating water temperature lines of the humidifier, shown in Figure 19, illustrate the condition of the humidification process within the operating oil flow, where the lower oil flow revealed a lower productivity and high PTC efficiency.Moreover, the temperatures of the circulating water in the humidifier were found to be higher at the oil flow rate of 0.0983 kg/s than with other settings, which shows the optimum oil flow rate for the heat exchange in the humidifier.The average PTC efficiency at the lowest oil flow rate is 27.8% and is reduced by 49.5% at the maximum flow rate.In the designed system, the ratio of the heat gained to the amount of beam radiation concentrated in the absorber tube is revealed to be higher, especially from the middle to the end of the solar day due to closed system circulation.Moreover, for all of the cases in which the incidence angle increased throughout the day, the losses due to the end effect decreased; hence, the PTC efficiency increased.When the beam radiation reduced at the end of the day, the system still can deliver the energy to the humidifier at a higher rate due to the stored heat in the oil expansion tank and the circulating air and water of the humidifier.The behavior of the circulating oil inlet and outlet temperatures in the parabolic concentrators is shown in Figure 22.Generally, at higher oil flow rates, the system heat storage capabilities increased due to multiple circulation procedures.This positively affects the value of the parabolic concentrator efficiency; hence, the minimum beam radiation utilization is obtained.The heat removal factor F R is known to be the measure of the thermal resistance by the absorbed radiation in the circulating oil [27].The heat removal factor changes at different oil and air flow rates, as shown in Figure 23.At very high oil flow rates, it possible that the oil outlet temperature becomes closer to that of the inlet; therefore, the factor is reduced.From the figure, at the highest oil flow rate of 0.12 kg/s, the oil temperature difference decreases and reduces the heat removal factor.In all the tested cases, the oil flow rate of 0.0983 kg/s was shown to have high system productivity and maintained oil temperature difference, and the higher air flow rates revealed higher values of the heat removal factor.The heat removal factor FR is known to be the measure of the thermal resistance by the absorbed radiation in the circulating oil [27].The heat removal factor changes at different oil and air flow rates, as shown in Figure 23.At very high oil flow rates, it possible that the oil outlet temperature becomes closer to that of the inlet; therefore, the factor is reduced.From the figure, at the highest oil flow rate of 0.12 kg/s, the oil temperature difference decreases and reduces the heat removal factor.In all the tested cases, the oil flow rate of 0.0983 kg/s was shown to have high system productivity and maintained oil temperature difference, and the higher air flow rates revealed higher values of the heat removal factor.

Absorber Tube Overall Heat Transfer Coefficient
The overall heat transfer coefficient, including the absorber tube wall, calculated using the heat removal factor, the heat removal coefficient and the convective heat transfer coefficient inside the absorber tube, is obtained from the recorded experimental data.The behavior of the overall heat transfer coefficient, at different circulating oil flow rates and humidifier air flow rates, is illustrated in Figure 24.The circulating oil inside the tube is found to be laminar in all the tested flows, and the corresponding Nusselt number was 4.364.The overall heat transfer coefficient depends on the conduction through the tube wall, the convection inside the tube and the heat removal loss coefficient, and the ambient conditions, such as the wind speed and temperature, influence its value.The convection coefficient for the laminar flow has less significance for the value of the overall heat transfer coefficient, and the wall conduction coefficient does not have an impact on changing the overall heat transfer coefficient.Therefore, the calculated heat removal coefficient profoundly influences the process affected by the ambient conditions.Increasing the circulating oil flow rate will increase the tendency toward a higher heat loss; hence, the overall heat transfer coefficient will increase.At a lower air flow rate, such as 0.0631 kg/s, the absorber tube will gain more heat at a lower oil residence time inside the absorber tube; hence, the system will have a tendency to increase the heat removal coefficient.The heat removal value could be decreased using a glass vacuum tube instead of a PTC glass cover.

Absorber Tube Overall Heat Transfer Coefficient
The overall heat transfer coefficient, including the absorber tube wall, calculated using the heat removal factor, the heat removal coefficient and the convective heat transfer coefficient inside the absorber tube, is obtained from the recorded experimental data.The behavior of the overall heat transfer coefficient, at different circulating oil flow rates and humidifier air flow rates, is illustrated in Figure 24.The circulating oil inside the tube is found to be in all the tested flows, and the corresponding Nusselt number was 4.364.The overall heat transfer coefficient depends on the conduction through the tube wall, the convection inside the tube and the heat removal loss coefficient, and the ambient conditions, such as the wind speed and temperature, influence its value.The convection coefficient for the laminar flow has less significance for the value of the overall heat transfer coefficient, and the wall conduction coefficient does not have an impact on changing the overall heat transfer coefficient.Therefore, the calculated heat removal coefficient profoundly influences the process affected by the ambient conditions.Increasing the circulating oil flow rate will increase the tendency toward a higher heat loss; hence, the overall heat transfer coefficient will increase.At a lower air flow rate, such as 0.0631 kg/s, the absorber tube will gain more heat at a lower oil residence time inside the absorber tube; hence, the system will have a tendency to increase the heat removal coefficient.The heat removal value could be decreased using a glass vacuum tube instead of a PTC glass cover.

Conclusions
This work has integrated modeling, techno-economic analysis, design, and experimental work for the assessment and optimization of a solar-driven humidification and dehumidification desalination unit for arid areas.The system is characterized by high efficiency, relatively simple operation and low maintenance.The GOR of the process, considering the primary purpose of the system construction, would have less priority compared to the other studied factors.Solar energy, concentrated by a PTC system, and an associated HDH desalination unit have been studied in order to evaluate the system productivity, efficiency, and economic feasibility.The design, construction, and testing have been implemented, and the following results have been obtained for a case study in Saudi Arabia:

•
The HDH system productivity is significantly influenced by the circulating air-to-water flow ratio, and the highest system productivity obtained is 24.31 kg/day.The best productivity curve for all the tested air flow rates is obtained at a 0.0983 kg/s circulating oil flow rate.The system

Conclusions
This work has integrated modeling, techno-economic analysis, design, and experimental work for the assessment and optimization of a solar-driven humidification and dehumidification desalination unit for arid areas.The system is characterized by high efficiency, relatively simple operation and low maintenance.The GOR of the process, considering the primary purpose of the system construction, would have less priority compared to the other studied factors.Solar energy, concentrated by a PTC system, and an associated HDH desalination unit have been studied in order to evaluate the system productivity, efficiency, and economic feasibility.The design, construction, and testing have been implemented, and the following results have been obtained for a case study in Saudi Arabia:

•
The HDH system productivity is significantly influenced by the circulating air-to-water flow ratio, and the highest system productivity obtained is 24.31 kg/day.The best productivity curve for all the tested air flow rates is obtained at a 0.0983 kg/s circulating oil flow rate.The system showed a drastic reduction in productivity of 70.1% when the air flow rate increased to 0.2683 kg/s.

•
Increasing the humidifier circulating water temperature positively affects the humidifier productivity.This has been achieved by modulating the heat input to the humidifier by altering the oil circulation through the PTCs.The maximum humidifier circulating water temperature with the highest productivity is 49.23 • C. The circulating water temperature is reduced by 21.7% when the highest circulating air flow rate of 0.2683 kg/s is achieved.

•
Increasing the circulating oil flow rate positively influences the system's productivity until the point when this starts to decline; hence, a tradeoff between either a larger size humidifier heat exchanger or larger PTC surface area is required to enhance productivity.

•
The energy consumption of the process under the best conditions varies between 2.39 and 3.07 kWh, and the energy delivered to the humidifier is decreased, when further circulating oil flow is increased, and reached a 39.3% reduction at the highest oil flow rate.

•
The humidifier efficiency ranges between 69.1% and 75.6%, and the thermal efficiencies are in the range of 76.6% and 82.94% at a circulating oil flow rate of 0.0983 kg/s.

•
The cost analysis of the system indicates a reduction of the production cost under the highest productivity condition.A tradeoff between the cost, productivity, and GOR is necessary.The production cost is 12.54 US$/m 3 under a condition where the highest GOR of 0.51 is obtained, and the maximum productivity is achieved.

•
The efficiency of the PTCs is calculated, and the results showed a significant increase in the performance at a low circulating oil flow rate.The oil circulation procedure in a closed loop enhances the system efficiency by reducing the losses and enhancing the energy storage in the oil expansion tank.
Finally, the PTC coupled with the HDH desalination system reveals a promising desalination technique in arid areas, where no commercial desalination plants are available for a small population.The system outcome can be improved by the further development of a hybrid system of PV, PTC, and auxiliary heaters in combination with an HDH desalination unit.

Recommendations for Future Work
The solar-driven HDH desalination process would be more effective, with a higher performance, cost-efficient and positive environmental impact in the desalination technology, if the following recommendations were implemented.

•
The solar-driven HDH system can be coupled with other thermal desalination technologies to increase the desalination process yield.The dual-purpose plants for power and desalination can be significantly improved.The optimization of the dual-purpose system can benefit from and be enhanced using a shortcut method to address the need for conceptual design studies [51].Moreover, decoupling the mass balances and the heat balances and heat-transfer sizing equations using the linear programming (LP) optimization formulation simplifies the computations and significantly reduces the model size and complexity [52].A hybrid design can include both thermal and reverse osmosis plants [53]; therefore, economic optimization can be carried out in the overall system to identify the process configuration tradeoffs in various scenarios.Mixed-integer nonlinear programming (MINLP) is used for 20 RO membrane modules, and 58% of savings in freshwater use is achieved as compared to the existing base case operations [54].The solution to this program provides the optimal arrangement, types, and sizes of the reverse-osmosis units and the booster pumps [55].Such an application showed an improvement in the optimization, determining the optimal mix of solar energy, thermal storage, and fossil fuel to attain the maximum annual profit of the shale gas production system studied [56], which can be applied to the solar-driven HDH system studied.

•
Heat integration can be performed to reduce the operating cost by coupling the HDH desalination system to industrial facilities.This approach has been implemented along with the optimization of multi-effect distillation (MED) and membrane distillation (MD) and revealed a successful cost reduction in different scenarios [57].Moreover, coupling multi-stage flash (MSF) with MED desalination processes was shown to improve the system output due to the use of the waste heat recovery [58].Therefore, the HDH desalination can be integrated with solar collectors, absorption refrigeration cycles, and organic Rankine cycles.

•
The existing desalination systems can consider the HDH as a modular desalination technology in the retrofitting option, which increases the capacity of the output [59].

•
The economic factors are the most important factors relating to project improvement; hence, the low cost-integrated solar HDH desalination system is an option that incorporates sustainability in the design and selection of projects.The metric approach can provide an evaluation of project profitability and promote sustainability related to the capital cost [60].

Figure 1 .
Figure 1.Average annual horizontal irradiation in Saudi Arabia.Reproduced with permission from the solar resource map © 2018 Solargis [5].

Figure 1 .
Figure 1.Average annual horizontal irradiation in Saudi Arabia.Reproduced with permission from the solar resource map © 2018 Solargis [5].

Figure 2 .
Figure 2. Schematic diagram of the humidifier.

Figure 2 .
Figure 2. Schematic diagram of the humidifier.

Processes 2018, 6 ,
x FOR PEER REVIEW 5 of 35 of water vapor being transferred to the dehumidifier and is extracted as a process yield.The fabricated and installed major components of the experiment are shown in Figure3.

Figure 3 .
Figure 3. Photos of the experimental setup and main components: (a) parabolic trough collector (PTC) with absorber tube; (b) humidification and dehumidification (HDH) desalination unit; and (c) oil expansion tank.

Figure 3 .
Figure 3. Photos of the experimental setup and main components: (a) parabolic trough collector (PTC) with absorber tube; (b) humidification and dehumidification (HDH) desalination unit; and (c) oil expansion tank.

Figure 4 .
Figure 4. Schematic diagram of the parabolic concentrator.

Figure 4 .
Figure 4. Schematic diagram of the parabolic concentrator.

Figure 7 .
Figure 7. Schematic diagram of the experiment controls.

Figure 7 .
Figure 7. Schematic diagram of the experiment controls.

Figure 8 .
Figure 8. Schematic of the humidifier balance.

Figure 8 .
Figure 8. Schematic of the humidifier balance.

Figure 9 .
Figure 9. Schematic of the dehumidifier balance.

Figure 9 .
Figure 9. Schematic of the dehumidifier balance.

Figure 10 .
Figure 10.The productivity of the system at different air and oil flow rates.

Figure 10 .
Figure 10.The productivity of the system at different air and oil flow rates.

Figure 11 .
Figure 11.Average HD water temperature at different air and oil flow rates.

Figure 12 .
Figure 12.HD average input energy at different air and oil flow rates.

Figure 11 .
Figure 11.Average HD water temperature at different air and oil flow rates.

Processes 2018, 6 , 35 Figure 11 .
Figure 11.Average HD water temperature at different air and oil flow rates.

Figure 12 .
Figure 12.HD average input energy at different air and oil flow rates.

Figure 12 .
Figure 12.HD average input energy at different air and oil flow rates.

Figure 13 .
Figure 13.HD average input energy at different air and oil flow rates.

Figure 13 .
Figure 13.HD average input energy at different air and oil flow rates.

Figure 14 .
Figure 14.Average humidifier efficiency at different air and oil flow rates.

Figure 15 .
Figure 15.Average humidifier thermal efficiency at different air and oil flow rates.

Figure 16 .
Figure 16.Average energy-based effectiveness ℇ at different air and oil flow rates.

Figure 14 . 35 Figure 14 .
Figure 14.Average humidifier efficiency at different air and oil flow rates.

Figure 15 .
Figure 15.Average humidifier thermal efficiency at different air and oil flow rates.

Figure 16 .
Figure 16.Average energy-based effectiveness ℇ at different air and oil flow rates.

Figure 15 . 35 Figure 14 .
Figure 15.Average humidifier thermal efficiency at different air and oil flow rates.

Figure 15 .
Figure 15.Average humidifier thermal efficiency at different air and oil flow rates.

Figure 16 .
Figure 16.Average energy-based effectiveness ℇ at different air and oil flow rates.

Figure 16 .
Figure 16.Average energy-based effectiveness ε at different air and oil flow rates.

Figure 17 .
Figure 17.Average enthalpy difference at different air and oil flow rates.

Figure 17 .
Figure 17.Average enthalpy difference at different air and oil flow rates.

Figure 18 .
Figure 18.Average HD mass transfer coefficient at different air and oil flow rates.

Figure 19 .Figure 18 .
Figure19.Average HD convective heat and mass transfer coefficient ratio at different air and oil flow rates.

Figure 18 .
Figure 18.Average HD mass transfer coefficient at different air and oil flow rates.

Figure 19 .Figure 19 .
Figure 19.Average HD convective heat and mass transfer coefficient ratio at different air and oil flow rates.

Figure 20 .
Figure 20.System optimization at different air and oil flow rates.

Figure 20 .
Figure 20.System optimization at different air and oil flow rates.

Figure 21 .
Figure 21.Average collector efficiency at different air and oil flow rates.

Figure 22 .
Figure 22.Average PTC oil temperature at different air and oil flow rates.

Figure 21 .
Figure 21.Average collector efficiency at different air and oil flow rates.

Figure 21 .
Figure 21.Average collector efficiency at different air and oil flow rates.

Figure 22 .
Figure 22.Average PTC oil temperature at different air and oil flow rates.

Figure 22 .
Figure 22.Average PTC oil temperature at different air and oil flow rates.

Figure 23 .
Figure 23.Heat removal factor FR of the PTCs with the change of system productivity and oil ΔT at different air and oil flow rates.

Figure 23 .
Figure 23.Heat removal factor F R of the PTCs with the change of system productivity and oil ∆T at different air and oil flow rates.

Figure 24 .
Figure 24.Average absorber tube overall heat transfer coefficient.

Figure 24 .
Figure 24.Average absorber tube overall heat transfer coefficient.

Table 1 .
The PTC geometrical parameters.

Table 1 .
The PTC geometrical parameters.

Table 2 .
Thermo-physical properties of the SAE 10W-30 synthetic blend engine oil. T

Table 3 .
Capital cost of the HDH unit.

Table 4 .
Comparison between the current result and previous work.
Water temperature, • C U Overall heat transfer coefficient, W/m 2 • C U c Overall heat transfer coefficient of the absorber tube, W/m 2 • C U L Heat removal coefficient of the absorber tube, W/m 2 • C V Volume of the humidifier, m 3 W Aperture width, m a Interface area, m −1 c a Specific heat of humid air, kJ/kg • C c w Specific heat of water, kJ/kg • C f c Enthalpy correction factor, kJ/kg h Hour angle, degrees h ab Oil convective heat transfer coefficient, kJ/s m 2 • C h c Convective heat transfer coefficient, kJ/s m 2 • C h fg Latent heat of evaporation of water at ambient conditions, kJ/kg h m Latent heat of evaporation of water at base temperature, kJ/kg