Next Article in Journal
New Adaptive Control Strategy for a Wind Turbine Permanent Magnet Synchronous Generator (PMSG)
Next Article in Special Issue
Energy Efficiency for Social Buildings in Morocco, Comparative (2E) Study: Active VS. Passive Solutions Via TRNsys
Previous Article in Journal
Deterioration of the Mechanical Properties of FFF 3D-Printed PLA Structures
Previous Article in Special Issue
A Review of Energy Storage Participation for Ancillary Services in a Microgrid Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Project for a Much More Diverse Moroccan Strategic Version: The Generalization of Solar Water Heater

1
Ecole Supérieure de Technologie de Fes, Université Sidi Mohamed Ben Abdellah USMBA, Route d’Imouzzer, Fès BP 2427, Morocco
2
Laboratoire des Sciences de l’Ingénieur Appliquées à la Mécanique et au Génie Electrique (SIAME), Université de Pau et des Pays de l’Adour/E2S UPPA, EA4581, 64000 Pau, France
3
Laboratoire de Mathématiques, Modélisation en Physique Appliquée, Ecole Normale Supérieure de Fès, U.S.M.B.A, Route Bensouda, Fez BP 5206, Morocco
*
Author to whom correspondence should be addressed.
Submission received: 16 October 2020 / Revised: 21 December 2020 / Accepted: 22 December 2020 / Published: 24 December 2020
(This article belongs to the Special Issue Emerging Technologies for the Energy Systems of the Future)

Abstract

:
This paper presents a strategical project for the new version of the Moroccan energy policy. It highlights the technology of solar water heaters (SWH), studying energy, economic and environmental gains of SWH generalization to satisfy the total resident need proposing a new strategic version diversified in terms of adopted technologies (more than green electricity). A detailed analysis of thermal performances and economic profitability of direct thermosyphon solar water heaters (TSWH) for residential requirements in Morocco. The optimum design parameters were defined and investigated using the dynamic TRNSYS simulation program. The optimum system was simulated under the six climatic conditions of Morocco in order to assess the related performances in terms of the collector efficiency and solar fraction. The major finding of this work is that large-scale integration of TSWH into Moroccan residences could provide up to 70% of thermal energy loads. An economic study was also developed to predict the life-cycle savings generated by the generalization of this technology in Morocco for all residential building’s categories. Approximately 1250 million USD as national saving on the total energy bill can be achieved. The environmental effects were also assessed to achieve the aims of this work and to evaluate the CO2 emissions avoided due to this environmentally friendly solution.

1. Introduction

Sun is an abundant alternative source of clean, renewable, and sustainable energy [1]. However, it was found that solar energy remains until now not fully exploited because of the high initial cost of solar thermal technologies [2], particularly in Morocco, a member country of the Middle East and North Africa (MENA) and considered as a case study in this work. In fact, one of the most valuable new policies is the integration of solar energy in the building sector to satisfy the needs of domestic hot water in Morocco, principally through solar thermal technologies. Indeed, thermosyphon solar water heaters (TSWH) can reduce intensively the consumption of butane gas, develop the market of solar energy in Morocco, and as a result widen the overall area of installed collectors. The use of such technologies can also participate to reach national targets related to energy efficiency (EE) in the building sector [3]. Finally, solar thermal heaters can enhance the safety aspect during the utilization phase by avoiding the phenomenon of fatal accidents due to the exhaust of butane gas which was estimated at 60 deaths per year according to Gergaud et al. [4].
TSWH has become the most used technology around the world to heat water. It is defined as a passive solar technology based on natural convection, without the necessity of using mechanical or electrical power to pump and circulate the fluid neither its control. Fluid’s circulation is the result of the density difference caused by the solar fluid heating considered as a driving force. Thermosyphon is described as the oldest technology which was introduced into the world market. Its performance is equal to the active systems and remains efficient in some cases. Hence, it can be concluded that the thermosyphon technology is the simpler system that could be designed [5].
The performance of solar water heaters has been developed and assessed theoretically and experimentally in several works. For example, Michaelides et al. [6] studied, different configurations of solar water heating systems, especially the under meteorological and socio-economic conditions of Cyprus, using different load profiles where the results showed that the solar fraction balanced between 89 and 63% depending on the consumption pattern studied (low and high). Moreover, Carlsson et al. [7] studied the thermosyphon system by replacing the traditional materials with polymeric ones, in order to compare the thermal, and financial performances. It was found that the total energy costs could be reduced considerably using the polymeric material. Further detailed analyses on Thermosyphon performances were also investigated by Kalogirou et al. [8]. Authors studied the performance of a TSWH considering a series arrangement and formulated a novel heat exchanger model to evaluate the effects of many parameters such as the inclination angle, coupling geometries and aspect ratio on the global performance. The most influential parameter was found to be the aspect ratio which affected to temperature and efficiency according a linear trend [9,10].
Many researchers have focused on innovative techniques to improve the efficiency of TSWH; Vasiliev [11] discussed using of new two-phase termosyphons technologies with advanced performances in order to study identical long construction with different working fluids, to examine the thermal performances especial the resistance which can be lowed using the vapordynamic thermosiphon (VDT) thernosyphons and polymeric loop two-phase under different parameters such us the inclination. The innovation continuous to appear in many other works, such as Piotr Felinski’s [12] study which applies the phase change material (PCM) and the paraffin in evacuated tube collectors for domestic hot water application, to evaluate the effect of this materials and thermal performance. The results showed that the use of evacuated tube collector (ETC) with PCM and paraffin allowed height water temperatures and improved solar fractions during the peak loads and the periods of lowest solar radiation intensities under the typical meteorological year conditions, compared to the conventional types. While the conventional ETC showed in several works right performances in terms of water temperatures, yields, energy and exergy efficiencies for many applications; thermosyphonic and forced models [13,14,15]. Koffi et al. [16] studied theoretically and experimentally a TSWH with an internal exchanger in order to assess the effect of some operating parameters on the outlet temperature and collector thermal efficiency. This study was investigated for Ivory Coast climatic conditions in 2014. Zerrouki et al. [17] studied the outlook for the use of a thermosyphon system manufactured in Algeria to predict its efficiency under the Algerien climatic conditions. Another study, with the same purpose was carried out under clear nights conditions and using the thermosyphon system with at-plat collectors by Runsheng Tang et al. [18]. Many other works were developed recently to study the general performances TSWH. For instance, Zeghib et al. [19] in Algeria presented a new theoretical simulation model, to predict the performance and the system behaviors under a thermosiphonic operation. The model was carried out using a flat plat collector with an area of 2 m2 for a volume of 200 L. Consequently, the theoretical analysis results height heat energy output, and satisfying efficiencies of collector and auxiliary heater improved by the good stratification of the model designed compared to the fully mixed conventional models. Vieira et al. [20] performed a multi-parametric study of SWH system under the climatic conditions of Brisbane in Australia, using EnergyPlus 8.6 software. The model’s calculations and analysis showed that the split systems perform better than thermosiphon system in terms of service level and energy efficiency. Another parametric study based on the load profiles influence was carried out by ERICH HAHNE [21]. The analysis was performed using a solar combisystem using a storage tank with an internal thermosiphonally driven discharge unit, operating under realistic profiles of 1-min scales and a constant total yearly demand. The results showed that the load profile have a severe influence on the system performances, this is explained by the significant role of the duration and the flow rates of DHW (Domestic Hot Water) on the stratification and temperatures. S. Fung et al. [22] presented two existing systems in his work, drain water heat recovery (DWHR) system and two solar domestic water heaters (SDWH) intended for two houses, the study aims to evaluate the performances and the recovery potential of systems, using two collector technologies; flat plat collector and ETC as a result the DWHR present an effectiveness of 50%, and the productivity of sensors was evaluated. Several works were investigated to study other types of solar water heaters. The analysis of Hobbi et al. [23] could be cited as example. Authors simulated and optimized a forced circulation using the flat plat collectors for residential SWH application in Montréal Canada, the study aims to show the importance of solar energy which could provide between 83–97% of the hot water demand in summer. The aim is to improve the efficiency of systems integrated in buildings, in order to reach economical and efficient buildings. Something that can be achieved with a multitude of actions [24,25]. Thermosyphon effect has several other technological uses, Chen and Yang [26] solve the efficiency problem for concentrating solar cells due to their height temperature, by the integration of the loop thermosyphon in the heat dissipation system. Different fluid for the operating thermosyphon loop were used. The results of theoretical and experimental analyses showed that the acetone is better than water and ethanol in term of heat transfer. Solomon et al. [27] focused on the closed system presenting the comparison between the two phase closed thermosyphon (TPCT), porous copper coating and the uncoated, the results showed that difference in heat transfer coefficients between the two cases at an angle of 45° is 44% at a heat flux of 10 kW.m−2. Furthermore, more the copper coating is thin more the wall temperature of the evaporator is significant and the heat transfer coefficient is increased, which makes the thin coating a suitable for cooling high-density power. Kousksou et al. [28] presented a numerical study which is concerned with the integration of phase change materials (PCMs) in solar-based domestic hot water (DHW) systems, so as to enhance its overall performances.
In the literature, many works join to the technical study an economic analysis to present a significant result, especially for the strategically and policies works [29]. Recent economic studies regarding solar water heaters integration were published. The techno-economic benefits and reliabilities of solar water heaters using the Monte Carlo analysis were estimated by Rezvani et al. [30]. They focused on a product range manufactured by a local company in Australia. Sokka et al. [31] focused on the environmental impacts of solar energy. Greenhouse gas emissions and mitigation of climate change reduction were studied which is the main aim of the global climate policy. Bouhal et al. [32] investigated the impact of collector technology and load profile on the fractional savings of solar domestic water heaters under various climatic conditions in Morocco. While an energy analysis of solar domestic heating water systems was conducted by Allouhi et al. [33].
In the present work, the preliminary choice of the technology to be integrated into the Moroccan context was made based on the results of the different systems studied in the work exhibited in the literature review previously mentioned, the particularity of this study is to integrate the Moroccan climate context, economic and social in the choice and the technical elaboration of the adapted model in order to present a realistic study combining the technology and the current conjuncture of Morocco as part of an African area and fertile ground of the generalization project. For Morocco, solar energy is a crucial economic issue in coherence with the choice of sustainable development and the related energy policy, which aims promoting renewable energies to meet 20% of the country’s domestic energy needs, improving energy efficiency to achieve 12% reduction by 2020. Solar energy utilization is expected to play a pivotal role in environment protection and social development [34,35]. This orientation leaded many actions, laws and projects to integrate green solutions precisely the solar within different sectors [36], especially the building sector which presents 20.4% of the national total energy consumption [37], this sector is characterized, firstly by a quick development because of fast population growth and an increasing urbanization rate [38,39].
Through this study, large-scale integration of TSWH into various types of buildings in Morocco is discussed energetically and financially based on climatic data of six different climatic zones of Morocco [35]. The performance indicators describing the system in this work are the solar fraction and collector efficiency. The simulation outputs were then used to develop the economic study, which aims to quantify the possible financial gains of generalizing the utilization of TSWH in the domestic sector of Morocco. This assessment was made by comparing the solar system with conventional gas heaters which is the most used heating option in the country. Two scenarios of subsidies elimination for conventional energy sources are discussed and the national potentiality of CO2 emissions reduction resulting from generalizing this kind of water heaters is assessed.

2. Moroccan Context Related to the Study Goal

In 2008, Morocco faced many challenges. An energy system marked by an extreme dependence on the outside world, a predominance of petroleum products and hydrocarbons in general, sustained growth in demand, increased rural electrification, and high price and price volatility. Despite the relevance of the 2009 strategy’s vision and the break it has initiated, with the desire to increase the share of renewable energy. The strategy made it possible to secure the supply, to initiate the liberalization of the electricity market and position Morocco at the forefront of the climate agenda, however this first version of Moroccan energy strategy have several lacks and brakes, such as the focus on the photovoltaique renewable technology; that requires the adoption of new strategic orientations.
This energy strategy adopted in 2009 has consolidated many achievements and must now be revisited. It is necessary to focus on exploiting this deposit to position the Morocco as the energy transition leader. The exploitation of the Moroccan potential conducted in an integrated and inclusive way should benefit the citizen, the economy, and the State.
Many actions will be taken in the new energy strategy version such as training, research, development and innovation: invest more in human capital to build a pool of business skills and networks of researchers and engineers engaged in a global dynamic, both national and regional. Further, we should promote the initiatives undertaken by the ecosystem approach of “Institut de Recherche en Energie Solaire et Energies Nouvelles” (IRESEN) and its partners (universities). Besides the diversification of the energy technologies adopted (more than photovoltaiques and electricity production) is one of the necessaries actions. In this sense, Morocco lunched several projects to support the solar thermal technologies and specially for residential uses. The PROMASOL (programme de développement du marché marocain des chauffe-eau solaires) project was one of the first lunched since the first version which aims to develop the Moroccan solar thermal market, but it didn’t make significant results. For this reason, Morocco lunched a new project about solar water heaters which aims to present a new Moroccan SWH with an accessible price and adapted technology to the economic and the climatic context.
For that, an overview of the solar thermal market in morocco is necessary to introduce the present study, and situate this project proposition of (generalization of SWH).
Morocco’s equipment level has always remained around four times lower than the global average; this means, notably, its equipment level remained parallel to the world’s equipment level. Despite its good sunshine, Morocco is still not among the countries installed capacity in 2018 but, Figure 1 shows that it is listed in the third group of countries (total of five), having between 20 and 70 m2 for 1000 habitants [40].
Morocco’s share of installed solar thermal capacity is low but remained stable for three decades as showed in Figure 2, around 0.10%. This means that the growth in m2 installed in Morocco has been parallel to the world trend [41].
On a 30-year scale, Figure 3 shows that: prices in Dh/m2, cost insurance and freight (CIF, blue diamonds in Figure 3), fluctuated by 15% (vertical blue bars in Figure 3) around a prolonged decrease of nearly 0.34% per year (decreasing blue curve in Figure 3). Solar thermal is no longer in the maturation phase and no longer undergoes a “learning curve.” Thus, prices would have lost only 5 to 6% in 30 years, which is low but appreciable considering inflation and increased purchasing power. The thermal energy produced by solar collectors has always been more profitable than that produced by grid electricity. Even very modest, this drop in prices combined with the rise in electricity prices has made this profitability even more attractive. Although “fair play,” this competitiveness goes unnoticed because of the formidable competitor of thermal energy produced by butane gas, which is heavily subsidized. Indeed, the withdrawal of the compensation 12 kg bottle would place its price around 120–140 MAD (Moroccan dirham) and not 40 MAD. Because of this, in Morocco, a solar water heater depreciates in 10 years against its butane gas competitor, regardless of the advantages related to comfort and safety of use.
The installation of solar thermal systems would have contributed to creating a few dozen jobs until around 2007. The number of jobs created amounted to a few hundred until around 2018 and seemed to stagnate around a little over 200 [42].
The annual demand has almost multiplied by 15 in 25 years or average growth of nearly 11% per year, which is not very fast compared to the growth of other renewable energies such as photovoltaics [43]. It is even surprising that the demand is maintained over such a long period despite the harsh competition of butane gas.
Domestic Solar water heaters (DSWH) have always dominated the sales of the market as described in Figure 4, making a significant turnover presented in Figure 5. As showed in Figure 6 they tend to dominate it more and more to exceed 80% after having fallen around 75% in the middle of the decade 2000–2010, probably under the PROMASOL Program’s impetus, which was led by AMEE (Agency for Control and Energy Efficiency).
All the elements presented describe the Moroccan market as a fertile ground for the insertion of this technology and popularizing it. This study assesses the positive impacts of the integration and generalization of this technology at the national level in order to improve the previously presented results further and generate a significant share of the economic and environmental gains not only by the production of green electricity but also by the production of domestic hot water and the coverage of all residential needs by solar.

3. Methodology

The TSWH is simulated using the TRNSYS (transient system simulation tool) program. The main operating and design parameters were introduced; the design parameters include general specifications of the used thermal collector, storage tank and auxiliary system. Weather data files for the various Moroccan climatic zones are then integrated into the simulation process. A realistic load profile is as well considered for accurate prediction of thermal loads and system performances. Based on the hourly output data, namely water outlet in the solar collector and temperature inside the tank, energetic performance indexes were built. These indices are the monthly solar fractions and collector thermal efficiencies.
Based on the evaluated useful energy produced by the TSWH and the total thermal loads required, a large-scale economic study is developed to examine the potential of generalizing these systems into the Moroccan residential sector. This evaluation is based on a comparison of cumulative energy savings generated by switching from the conventional gas heaters to TSWH. Finally, it was possible to estimate the overall reductions in CO2 emissions.
The assessment is believed to be sufficiently accurate as it follows a rational methodology for quantifying the full-scale economic impact of generalizing these solar systems. First of all, it was necessary to evaluate the number of capita per climatic zone, this parameter is not yet available, it is available only in a regional database according to the High Commission of Plans. Therefore, it was essential to convert the available regional data to comply with the climatic zoning. A schematic diagram of the followed methodology is shown in Figure 7.

4. Simulation Process

4.1. Examined Configuration

The TSWH consists of a flat plate collector (FPC) connected to an auxiliary heater integrated inside a horizontal tank as presented in Figure 8. For this system, water flows through the system when warm water rises as cooler water sinks without the need for a circulation pump. The FPC collector is installed below the storage tank so that warm water will rise into the tank by the thermosiphon effect, its principle operates on the heated water’s characteristic, which becomes lighter than the cold water, ascends to the horizontal tank, replaces the cold water, which, heavier, descends downwards and passes through the solar panel.
First of all, optimum tilt angles per each zone were determined. These angles are summarized in Table 1.
Moreover, the technical characteristics of the solar collector used during the simulations under TRNSYS are presented in the Table 2.

4.2. Climatic Data and Geographical Specification

The Moroccan climatic conditions used to assess the thermal performance of the thermosyphon are presented in Table 3 and Figure 9. This new climatic zoning is established by the Moroccan Agency of Energy Efficiency (AMEE).
Because the Moroccan climate exhibits a great variability according to the geographical zone, it was necessary to perform a segmentation of the whole territory into areas with similar overall meteorological specifications. Accordingly, the Moroccan Agency of Energy Efficiency (AMEE) [44], proposed a climatic fragmentation of Morocco according to six climates, each climate is evidenced by a representative city. Table 2 shows general information about these zones. In addition, the main meteorological data (i.e., monthly average ambient temperature and solar irradiations) affecting the energetic performance of TSWHs is displayed in Figure 9. One can observe that the climate presents a remarkable variety which is expected to influence greatly the overall thermal behavior and performance of TSWHs.

4.3. Load Profile

The distribution of the hourly hot water consumption is affected by different factors. It varies according to the seasons, days, or the habits of each considered Moroccan family. Various hot water load profile has been studied on the literature review such as the constant, the early morning, the late morning, rand, or the late afternoon profiles. For the present study, the consumption of 202 L per day is described according to the realistic load profile presented in the Figure 10.

4.4. Physical Modeling

The type 45 operates according to a mathematical model [45] which presents the various phenomena describing its operation. The Bernoulli formula can be applied to define the pressure drop in the thermosyphon system’s nodes (Equation (1)):
Δ P i = ρ i g Δ h i + ρ i g h L i
where (i) is the index of the node, Δ h i the height of the ith node, ρ is the density, h L i is the frictional head loss in the piping, and g is the gravitational constant.
The flow rate on the thermosyphon system must be involved in order to satisfy the conditions describing the total pressure differences, which are determined using the equation above at any time by the next expression (2):
i = 1 i = N ρ i Δ h i = i = 1 i = N ρ i h L i
Many types of losses are present in the system; the most important one is the frictional head loses especially on the pipes, and the equation below describes this type of losses (3):
H p = f L v 2 2 d + K v 2 2
The parameter F is the friction factor, determined according to the value of the Reynolds coefficient Re which describe the nature of the flow and corrected to take in account the frictions on the connecting parts of pipes. The head loses pipe varies according to the position and the conditions in the system, these specifications are taken into account by the coefficient K, e.g., loses due to different types of bends, the parts of the tank which are connected to the collector.
The thermal performance of collector is modeled according to the Hottel–Whillier equation. Many parameters describe the system: the equation bellow (4) defines the F′UL parameter which is calculated by the FRUL and G parameters determined on the test conditions.
F U L = G t e s t   C P ln ( 1 F R U L G t e s t C p )
For the collector the temperature in the Kth node is calculated by the next expression (5):
T c k = T a + I T F R ( τ α ) F R U L + ( T c i T a I T F R ( τ α ) F R U L )   exp [ F U L G C p   ( k 1 / 2 ) N x ]
where Nx is the number of nodes over which the size of the collector deviates to obtain the weight of the fluid in the collector.
The changes in the fluid heat transfer coefficient results in F′ and UL which are considered negligible. The combination of the intercept efficiency at normal incidence FR (τα)n, and the incidence angle modifier (IAM).
((τα))n allows us to calculate the parameter FR(τα) (6):
( τ α ) ( τ α ) n = I b T ( τ α ) b ( τ α ) n + I d 1 + cos β 2 ( τ α ) s ( τ α ) n + I g 1 + cos β 2 ( τ α ) g ( τ α ) n I T
where the incidence angle modifiers for ground and sky diffuse radiation is calculated using the next Equation (7):
( τ α ) b ( τ α ) n = 1 b 0 ( 1 cos θ 1 )
The overall useful energy collected is determined using the Equations (8) and (9):
Q u = r   A c ( F R ( τ α )   I T F R U L ( T c i T a ) )
r = F R ( u s e ) F R ( t e s t ) = G ( 1 exp ( F U L G C P ) ) G t e s t ( 1 exp ( F U L G t e s t C P ) )
The outlet temperature from the flat plat collector is finally deduced as (10):
T C O = Q u m ˙ C p + T C i
The thermosyphon diagram executed under TRNSYS program is shown in Figure 11. The different module’s types used in this study are presented as follows:
  • Weather data (TYPE109-TMY2) this component reads meteorological data generated by meteonorm software and supply them in the TMY2 format. It calculates the necessary solar radiation for the calculation of TRNSYS software at any surface tilt.
  • Differential controller for temperature (TYPE 2). This component monitors and controls the average temperature of the tank, to control the operation of the auxiliary electric heater.
  • The general forcing function (TYPE 14) from TRNSYS library, this component describes and characterizes the hourly load profile corresponding to 202 L/day of hot water demand.

4.5. Performance Indices

4.5.1. Collector Efficiency

The collector efficiency describes the ratio of the total useful energy gain to the solar energy absorbed by the collector. It is calculated using the equation below (11):
E c o l l = Q u c o l l / ( A c   I c o l l )
Qucoll is the energy rate from heat source, and Icoll is the total radiation on the tilted surface.

4.5.2. Solar Fraction

The solar fraction Fsol represents the solar energy available to meet the needs. It presents a clear contribution indication of the solar system in meeting the thermal load. The solar fraction is calculated using the equation below (12):
F s o l = 1 ( Q a u x / Q D H W )
where QDHW is the energy rate delivered to the load describing the necessary energy to meet all domestic hot water needs, and Qaux is the auxiliary energy. It is important to state that the solar fraction is evaluated on a monthly basis.

5. Model Validation

The developed model was simulated with firstly the meteorological data of Cyprus in order to carry out a comparison with the study performed by Michaelides et al. study [6]. A typical meteorological year (TMY) for Nosica was built for this end.
The system simulation was lunched for similar conditions as indicated in [6]. Table 4 shows the technical specifications of the studied TSWH while Table 4 and Table 5 indicates the general assumptions concerning technical parameters of validation case, the load profile and water consumption.
Figure 12 presents a comparison between monthly solar fractions of the TSWH investigated by Michaelides et al. study [6] and those predicted by the introduced TRNSYS model.
The average relative error is about 5%, which is an acceptable value. The deviation between the results can be attributed to the update in climatic conditions caused when using the Meteonorm database that corresponded to the period 1986–1992 for the study of Michaelides et al. [6] while the prepared metorological data file used in our calculations corresponds to the period 1991–2010.

6. Results and Discussion

6.1. Energetic Evaluation

As presented in Figure 13a, the thermosyphon system performs differently in the six cases. The average solar fraction ranges between 46% to 71%, and the maximum solar fraction value was noticed in Errachidia 71%. This result is explained by the importance of solar irradiations in this zone. The annual average solar fraction is about 69% in the first zone represented by Agadir which is a sunny coastal zone. It is in general the most favorable zone for this system, especially in winter.
Tangier (Zone 2) and Marrakech (Zone 5) exhibit approximately the same solar fraction variation, but (Zone 5) has a higher annual average value of 66% versus 61% in (Zone 2). This difference is due to the favorable conditions of Marrakech which is considered as one of the sunniest areas in Morocco. Zone 3 presents a total annual average solar fraction of 52%, which is lower compared to the previously discussed sites. This result can be explained by the fact that Fez is characterized by a Mediterranean climate, dominated by continental and Atlantic effects, which accounts for the remarkable fall of the solar fraction to 17.4% during the coldest month of winter. Finally, the coldest zone, Ifrane has an averaged solar fraction value of only 46%. It is concluded that Ifrane’s weather is the most unfavorable climate for installing solar TSWHs because Ifrane remains the coldest zone in Morocco with an average annual temperature of 10.8 °C. In addition, one can observe that most of the hot water requirements in the Moroccan zones, previously cited, are totally satisfied during the sunniest three months of the year (July, August, and September) with a solar fraction achieving 100%.
Figure 13b presents the dynamic behavior of the collector efficiency for the six zones. The most favorable zone, where the collector performs with high efficiency is Zone 6 represented by Errachidia, with a maximum value of 86% in July, and an average annual value of 74%. The collector operates in Agadir with 70% of efficiency, in Marrakech with 66%, Fez and Tangier with 54% and 59% respectively and finally in Ifrane with 51%. These results confirm the explanations previously presented for solar fraction. In fact, when the thermal efficiency of the collector is maximized, the useful thermal energy transmitted to the working fluid is as well maximized and thus less auxiliary energy will be required.

6.2. Economic Assessment

This section aims to assess the total financial gains generated by the integration and the generalization of solar water heaters in Morocco into the residential sector, according to the six climatic zones and for different house’s types. The solar system will be compared to a reference system which is the conventional gas boiler according to two scenarios:
  • The first scenario takes into account the state subsidies on gas, which presents a purchasing power support of households. In this case, the consumer benefits on each bottle of gas that purchased a reduction of 67% on its price.
  • The second scenario does not consider the subsidies.
The number of thermosyphon systems required to be installed was determined according to the situation of the building in Morocco. The number of houses counted on 2014 is 5.8 million, distributed as follow: occupied houses, secondary, vacant and professional use houses as presented in Table 6.
The number of occupied houses in Morocco is 4.582 million houses: this value is officially distributed over 12 regional zones according to (the High Commissariat of the Moroccan Plan) [34].
A conversion of this distribution according to the climatic zones was investigated as presented in Table 7. It is found that this number of houses per climatic zones can be also fragmented into many types such as: modern Moroccan houses which represent the main construction type in Morocco with a percentage of 63%, followed by apartments (25%), villas 4%, traditional Moroccan houses 4%, and finally anarchic construction 4%.
These statistics will enable the determination of the thermosyphon’s number to be integrated at a national scale.
The study will focus on apartments, modern Moroccan houses, and villas. It should be highlighted that, because of the presence of these variants in the construction sector, it is recommended to consider different capacities related to the type of house for an appropriate analysis. Table 8 presents the system capacity concerning the construction type. The generalization of this system into the building sector will generate an important total investment of 5876.42 million USD. This total investment is found to be 1145.5 million USD, 4144.42 million USD, and 586.50 million USD for modern Moroccan houses (150 L), apartments (200 L), and villas (300 L) respectively.

6.3. Total Annual Gains

The total gains or the life cycle savings of the thermosyphon generalization in Morocco by using the subsidies and non-subsidies gas butane prices is presented in Figure 7. This indicator represents the global cost savings resulting from the integration of the thermosyphon system instead of buying gas butane cylinders.
It is concluded, from the comparison of these two scenarios, that, obviously, the gains generated considering subsidies are less than the second case for the six zones. In fact, the government subsidies (on fuels) policy is representing a barrier towards development solar energy projects in Morocco. Moreover, the first zone represents the most important gains by 197.11 million USD in the first scenario against only 44.32 million USD in the second zone. This difference can be explained by the highest solar fraction (68.8%), together with the high intensity of the population concentration with a percentage of 46%. In the zone represented by Agadir, this optimal combination of solar fraction and the number of homes for the first climatic zone gives optimized results and very encouraging gains as it is previously represented. This value becomes significant and would reach 448 million USD per year by canceling government grants.
The third zone is holding the second place with total gains of 69.34 million USD per year. It can be improved by removing the state subsidies to reach 214.63 million USD per year. The other climatic zones achieve lower gains than the previous ones, despite their significant solar fractions which reach 66% (Errachidia) representing the most favorable climate conditions for the system. This result is explained by the minimal percentage of houses.
The energy gains generated by the generalization of thermosyphon systems are presented in Figure 14b. The first zone achieve an energy annual gain of 5943.46 GWh per year, which is distributed according to the different studied configurations: apartment with 1214.60 GWh per year, modern Moroccan houses by 4340.19 GWh annually, and 388.67 GWh generated for the villas type. The integration of solar water heaters in the modern Moroccan houses generates the maximum gains in the six climate zones compared to other types due to the remarkable presence of this type in the typological distribution of housing in the Moroccan territory. In fact, it represents 63% of the total occupied houses, followed by the buildings and finally villas.
Figure 15a–c present the evolution of the annual gains in an interval of 12 years, and under the condition of (non-subsidies) for several type of houses treated in this paper.
The payback period for the proposed sizes and under the six climatic data condition can be determined based on Figure 14. The return on investment (ROI) varies between 3 and 6 years according to the performances of the system related to the meteorological context, the investment and the gains generated which are related also to the performance and the total expenditures produced by the traditional gas butane system. In general, after a period of 6 years the thermosyphons installed in the climatic zones will cover the initial investment and in general they will enable 100% of the life cycle savings.

Global Results: Cumulative Gains for Different House’s Types

The overall average solar fraction presents a national value which is estimated at 62%. Table 9 presents the total gains for different type of house according to the two scenarios (with and without subsidies): the generalization of solar water heaters have an important role to play in the achievement of the development program’s goals of the Moroccan market for solar water heaters by (PROMASOL).
This program aims to provide the evolution continuity of thermosyphon installed capacity since 2000 obtained by the analysis of imports conducted by the (AMEE) which is a result presented in Figure 16. This figure is presenting the cumulative area of the SWH installed in Morocco since 2000 to the expectations desired to achieve.

6.4. Environmental Assessment

Using solar water heaters will also make it possible, in contrast to the conventional water heater, to reduce emissions of greenhouse gases. It should be recalled that recently, Morocco announced its engagement to reduce 13% of these gases by 2030, and ratified the Kyoto Protocol, which entered into force in 2005 Table 10 presents the annual CO2 emissions avoided by the integration and the generalization of the TSWH.
The maximum CO2 emissions avoided are remarkably visualized on the first climatic zone which presents the highest number of houses, and as a result the highest number of SWH installed, and energy savings.
The calculation of the CO2 emissions avoided is made based on the quantity of the butane gas consumed in the conventional system; this consumption of the butane gas will be eliminated by implementing the generalization project. This step required adopting the butane gas characteristics where it generates 230 g of CO2 emissions for each kWh produced. Furthermore, for bottles of 6 kg, the discharges correspond to 157 tons of CO2 equivalent. The bottles of 10 kg, under the same conditions, give 160 kg of CO2 equivalent. They finally amount to a bottle of 13 kg to 150 kg of CO2 equivalent.

7. SWH Generalization a Bridge Idea to a Promoter Economic Project (Manufacturing Moroccan SWH)

Nowadays, the global perspective focuses on the economy, the human, energy, and the environment, achieving optimum energy consumption and reducing greenhouse gas emissions. Besides, with this 2020 pandemic, the economic sector’s growth has become more urgent than before.
In Morocco, the project of this paper can be realized currently only by the import of this technology; this requires reflection to a local production unit of solar water heaters. Here, the idea of an innovative project was born (SOL’R SHEMSY project). We are working through the development of a Moroccan solar water heater ready to be industrialized locally. Hence, the project integrates all these goals and offers a balanced solution; technology, economy, and social and environment.
The local production will present solar water heaters with affordable prices adapted to the social context, breaking the most significant barrier of this technology’s extensive integration, which is the high prices resulting from the importation effect.
The manufacturing project will bring various high-quality products intended for multiple customer categories to meet all social and industrial clients’ needs, presenting the high technology with the best efficiencies, yields, and reasonable prices adapted to the Moroccan purchasing power. In particular, the main obstacle to its adoption is the high-security and weak support of African states. In economic and social terms, this project will create permanent jobs and achieve significant energy gains by reducing the energy bills costs for a range of sectors: housing sector, hotels, hospitals, schools and universities, sports sector (heated swimming water), industrial and agricultural sector. Strengthening the economic Moroccan–African industry through the manufacturing and marketing of solar water heaters and why not export. Finally, generalizing the use of these products would significantly reduce carbon dioxide emissions and combat the dangers of the gas heaters that kill dozens a year.

8. Conclusions

In this paper the thermosyphon solar water heater was presented, modeled, simulated and assessed according to six Moroccan climatic zones which represent the recent zoning carried out by the Moroccan Agency of Energy Efficiency (AMEE), in order to predict the system performances, and the dynamic behaviors of different parameters which are describing the thermosyphon. The system was modeled and simulations where launched under TRNSYS simulation program.
According to the simulation results the thermosyphon’s performances in terms of the annual average solar fraction varies between 50% and 70%. The results confirmed that the climatic conditions are the key of the simulation. They have to be taken into consideration to optimize the overall system efficiency. Moreover, an economic study which is about the generalization of this technology was conducted to analyze the potential benefits offered by this solar solution in different zones in Morocco from a weather point of view, and to predict the total life-cycle savings. In fact, life-cycle savings were optimal in the first zone represented by Agadir, due to the optimal combination of performances and high needs that are concentrated in this zone. To complete this study, an environmental section was developed to exhibit the amount of CO2 emissions that can be avoided by the implantation of this project. The results confirmed that a considerable amount of greenhouse gas emissions is reduced.
The penetration of large-scale solar water heaters in Morocco is prevented by many barriers, such as the understanding’s lack among the population of the financial unattractiveness to poor households and the benefits they provide. However, thermosyphon will be accepted easier when the population understands all these technology advantages. In the case of SWH systems, the benefits are not immediately visible for a payback period of six years: the energy bill will be reduced, but not fast enough in order to justify the cost of installation. Thus, the governmental subsidies will be an important incentive solution for beneficiaries. The generalization of this type of installation will not be possible without the introduction of incentive mechanisms that involves the public authorities, subsidizing the purchase of solar water heaters which will enable citizens to equip themselves at a lower cost.
The advantages linked to the generalization of the thermosyphon in Morocco are multiple, both economically and environmentally. Citizens could save their energy bill in the long term, since after an average of six years they would have profitable purchase price. Moreover, manufacturing the solar water heater in Morocco will lower the import bill cost. Furthermore, it will also help the development of subcontracting and involve SMEs (Small and Medium Enterprises) in the manufacture of this product. Hence, entrepreneurship in this solar field will grow and be encouraged. Therefore, jobs will be created which will absorb unemployment.

Author Contributions

Conceptualization, F.Z.G. and A.A.; Methodology, F.Z.G. and A.A.; Software, F.Z.G. and A.A.; Supervision, T.K., A.J. and A.B.; Visualization, H.E.-H.; Writing—original draft, F.Z.G. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

SOL’R SHEMSY project is financed by IRESEN(Institut de Recherche en Energie Solaire et Energie Nouvelle).

Data Availability Statement

Data available in a publicly accessible repository.

Acknowledgments

The authors are grateful for the funding support provided by the Institut de Recherche en Energie Solaire et Energies Nouvelles (IRESEN) under the project SOL’R SHEMSY.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AcCollector area (m2).
FRCollector heat removal efficiency factor (-).
FRULSlope of the collector efficiency curve (kJ h−1 K−1 m−2).
F R ( τ α ) Intercept efficiency corrected for non-normal incidence (-).
FFriction factor (-).
F′Collector efficiency factor (-).
ULHeat loss coefficient (kJ h−1 K−1 m−2).
GtestCollector mass flux at test conditions (kg h−1 m−2).
GCollector flow rate per unit area (kg h−1 m−2).
ρ i Density of ith node (kg m−3).
Δ h i Height of the ith node (mm).
Δ P i Change in pressure across the ith node (bar).
h Li Frictional head loss in the ith node on the pipe (-).
HpThe friction head losses on the pipes (-).
r Ratio of collector heat removal efficiency factor, FR, to the test conditions (-).
N x Number of equal sized collector nodes (-).
KCorrection coefficient for additional friction due to developing flow in the pipe (-).
kNode number (-).
ReReynolds coefficient (-).
f Friction factor for flow in pipes (-).
v Velocity of the fluid (m s−1).
m ˙ Thermo siphon flow rate (kg h−1).
C p Specific heat of working fluid (kJ kg−1 °C−1).
g Gravitational constant (N m2 kg−2).
( τ α )Product of the cover transmittance and the absorber absorptance (-).
( τ α ) b ( τ α ) for beam radiation depends on the incidence angle (-).
( τ α ) n ( τ α ) at normal incidence (-).
( τ α ) s ( τ α ) for sky diffuse radiation (-).
( τ α ) g ( τ α ) for ground reflected radiation (-).
IbTBeam radiation incident of the solar collector (kJ h−1 m−2).
IdDiffuse horizontal radiation (kJ h−1 m−2).
ITGlobal radiation incident on the solar collector (kJ h−1 m−2).
I T Total incident radiation per unit area (kJ h−1 m−2).
β Collector slope above the horizontal plane (°).
θ Incident angle for beam radiation (°).
b0Negative of the 2nd order coefficient in the IAM curve fit equation (-).
IAMIncident angle modifier (-).
T a Ambient temperature (°C).
T ci Collector inlet temperature (°C).
T CO Collector outlet temperature (°C).
LhLength of collector headers (mm).
LcoLength of collector outlet (mm).
LciLength of collector inlet (mm).
di,doDiameter of collector inlet and outlet pipes (mm).
Li, L0Length on inlet and outlet piping (mm).
HcVertical distance between outlet and inlet of collectors (mm).
HOVertical distance between outlet of tank and inlet to collector (mm).
DrRiser tank diameter (mm).
DhHeader tank diameter (mm).
VTank volume (L).
LPipe length (mm).
dPipe diameter (mm).
NcNumber of parallel collector risers (-).
UapiConductance for heat losses from collector inlet pipe (kJ h−1 k−1).
UApoConductance for heat losses from collector outlet pipe (kJ h−1 k−1).
(UA)sOverall UA value for storage tank (kJ h−1 k−1).
Qaux,maxMaximum heating rate (kW).
PmaxMaximum pressure (bar).
KTThermal losses of the tank ((kJ h−1 m−2).
CscCapacity of solar circuit (L).
EcollEfficiency of the collector (%).
FsolSolar fraction of the system (%).
QauxAuxiliary heating load (kJ h−1).
QDHWEnergy rate to load (kJ h−1).
QucollEnergy rate from heat source (kJ h−1).
IcollTotal radiation on tilted surface (kJ h−1 m−2).
ROIReturn on investment (years).

References

  1. Williams, J.R. Solar energy: Technology and applications. STIA 1974, 75, 12425. [Google Scholar]
  2. Ong, K. An improved computer program for the thermal performance of a solar water heater. Sol. Energy 1976, 18, 183–191. [Google Scholar] [CrossRef]
  3. Bleu, P. Stratégie méditerranéenne pour le développement durable. In Mise à Jour des Indicateurs de Suivi; 2013. [Google Scholar]
  4. Gergaud, O. Modélisation Énergétique et Optimisation Économique d’un Système de Production Éolien et Photovoltaïque Couplé au Réseau et Associé à un Accumulateur; 2002. [Google Scholar]
  5. Anderson, R.; Kreith, F. Natural convection in active and passive solar thermal systems. In Advances in Heat Transfer; Elsevier: Amsterdam, The Netherlands, 1987; Volume 18, pp. 1–86. [Google Scholar]
  6. Michaelides, I.; Lee, W.; Wilson, D.; Votsis, P. Computer simulation of the performance of a thermosyphon solar water-heater. Appl. Energy 1992, 41, 149–163. [Google Scholar] [CrossRef]
  7. Carlsson, B.; Meir, M.; Rekstad, J.; Preiß, D.; Ramschak, T. Replacing traditional materials with polymeric materials in solar thermosiphon systems–Case study on pros and cons based on a total cost accounting approach. Sol. Energy 2016, 125, 294–306. [Google Scholar] [CrossRef]
  8. Kalogirou, S.A.; Papamarcou, C. Modelling of a thermosyphon solar water heating system and simple model validation. Renew. Energy 2000, 21, 471–493. [Google Scholar] [CrossRef]
  9. Tecchio, C.; Oliveira, J.; Paiva, K.; Mantelli, M.; Gandolfi, R.; Ribeiro, L. Thermal performance of thermosyphons in series connected by thermal plugs. Exp. Therm. Fluid Sci. 2017, 88, 409–422. [Google Scholar] [CrossRef]
  10. Xu, Z.; Zhang, Y.; Li, B.; Wang, C.-C.; Ma, Q. Heat performances of a thermosyphon as affected by evaporator wettability and filling ratio. Appl. Therm. Eng. 2018, 129, 665–673. [Google Scholar] [CrossRef]
  11. Vasiliev, L.L.; Grakovich, L.; Rabetsky, M.; Vassiliev, L.L., Jr.; Zhuravlyov, A. Thermosyphons with innovative technologies. Appl. Therm. Eng. 2017, 111, 1647–1654. [Google Scholar] [CrossRef]
  12. Feliński, P.; Sekret, R. Effect of PCM application inside an evacuated tube collector on the thermal performance of a domestic hot water system. Energy Build. 2017, 152, 558–567. [Google Scholar] [CrossRef]
  13. Kumar, S.; Dubey, A.; Tiwari, G. A solar still augmented with an evacuated tube collector in forced mode. Desalination 2014, 347, 15–24. [Google Scholar] [CrossRef]
  14. Hassanien, R.H.E.; Li, M.; Tang, Y. The evacuated tube solar collector assisted heat pump for heating greenhouses. Energy Build. 2018, 169, 305–318. [Google Scholar] [CrossRef]
  15. Essa, M.A.; Mostafa, N.H. Theoretical and experimental study for temperature distribution and flow profile in all water evacuated tube solar collector considering solar radiation boundary condition. Sol. Energy 2017, 142, 267–277. [Google Scholar] [CrossRef]
  16. Koffi, P.; Andoh, H.; Gbaha, P.; Touré, S.; Ado, G. Theoretical and experimental study of solar water heater with internal exchanger using thermosiphon system. Energy Convers. Manag. 2008, 49, 2279–2290. [Google Scholar] [CrossRef]
  17. Zerrouki, A.; Boumedien, A.; Said, N.; Tedjiza, B. Input/output test results and long-term performance prediction of a domestic thermosiphon solar water heater in Algiers, Algeria. Renew. Energy 2002, 25, 153–161. [Google Scholar] [CrossRef]
  18. Tang, R.; Cheng, Y.; Wu, M.; Li, Z.; Yu, Y. Experimental and modeling studies on thermosiphon domestic solar water heaters with flat-plate collectors at clear nights. Energy Convers. Manag. 2010, 51, 2548–2556. [Google Scholar] [CrossRef]
  19. Zeghib, I.; Chaker, A. Simulation of a solar domestic water heating system. Energy Procedia 2011, 6, 292–301. [Google Scholar] [CrossRef] [Green Version]
  20. Shrivastava, R.; Kumar, V.; Untawale, S. Modeling and simulation of solar water heater: A TRNSYS perspective. Renew. Sustain. Energy Rev. 2017, 67, 126–143. [Google Scholar] [CrossRef]
  21. Jordan, U.; Vajen, K. Influence of the DHW load profile on the fractional energy savings: A case study of A solar combi-system with TRNSYS simulations. Sol. Energy 2001, 69, 197–208. [Google Scholar] [CrossRef]
  22. Tanha, K.; Fung, A.S.; Kumar, R. Performance of two domestic solar water heaters with drain water heat recovery units: Simulation and experimental investigation. Appl. Therm. Eng. 2015, 90, 444–459. [Google Scholar] [CrossRef]
  23. Hobbi, A.; Siddiqui, K. Optimal design of a forced circulation solar water heating system for a residential unit in cold climate using TRNSYS. Sol. Energy 2009, 83, 700–714. [Google Scholar] [CrossRef]
  24. Ge, J.; Wu, J.; Chen, S.; Wu, J. Energy efficiency optimization strategies for university research buildings with hot summer and cold winter climate of China based on the adaptive thermal comfort. J. Build. Eng. 2018, 18, 321–330. [Google Scholar] [CrossRef]
  25. Bellia, L.; Borrelli, M.; De Masi, R.F.; Ruggiero, S.; Vanoli, G.P. University building: Energy diagnosis and refurbishment design with cost-optimal approach. Discussion about the effect of numerical modelling assumptions. J. Build. Eng. 2018, 18, 1–18. [Google Scholar] [CrossRef]
  26. Chen, S.; Yang, J. Loop thermosyphon performance study for solar cells cooling. Energy Convers. Manag. 2016, 121, 297–304. [Google Scholar] [CrossRef]
  27. Solomon, A.B.; Daniel, V.A.; Ramachandran, K.; Pillai, B.; Singh, R.R.; Sharifpur, M.; Meyer, J. Performance enhancement of a two-phase closed thermosiphon with a thin porous copper coating. Int. Commun. Heat Mass Transf. 2017, 82, 9–19. [Google Scholar] [CrossRef]
  28. Kousksou, T.; Bruel, P.; Cherreau, G.; Leoussoff, V.; El Rhafiki, T. PCM storage for solar DHW: From an unfulfilled promise to a real benefit. Sol. Energy 2011, 85, 2033–2040. [Google Scholar] [CrossRef]
  29. Dwaikat, L.N.; Ali, K.N. The economic benefits of a green building–Evidence from Malaysia. J. Build. Eng. 2018, 18, 448–453. [Google Scholar] [CrossRef]
  30. Rezvani, S.; Bahri, P.; Urmee, T.; Baverstock, G.; Moore, A. Techno-economic and reliability assessment of solar water heaters in Australia based on Monte Carlo analysis. Renew. Energy 2017, 105, 774–785. [Google Scholar] [CrossRef] [Green Version]
  31. Sokka, L.; Sinkko, T.; Holma, A.; Manninen, K.; Pasanen, K.; Rantala, M.; Leskinen, P. Environmental impacts of the national renewable energy targets–A case study from Finland. Renew. Sustain. Energy Rev. 2016, 59, 1599–1610. [Google Scholar] [CrossRef]
  32. Bouhal, T.; Agrouaz, Y.; Allouhi, A.; Kousksou, T.; Jamil, A.; El Rhafiki, T.; Zeraouli, Y. Impact of load profile and collector technology on the fractional savings of solar domestic water heaters under various climatic conditions. Int. J. Hydrog. Energy 2017, 42, 13245–13258. [Google Scholar] [CrossRef]
  33. Allouhi, A.; Jamil, A.; Kousksou, T.; El Rhafiki, T.; Mourad, Y.; Zeraouli, Y. Solar domestic heating water systems in Morocco: An energy analysis. Energy Convers. Manag. 2015, 92, 105–113. [Google Scholar] [CrossRef]
  34. Arouri, M.E.H.; Youssef, A.B.; M’henni, H.; Rault, C. Energy consumption, economic growth and CO2 emissions in Middle East and North African countries. Energy Policy 2012, 45, 342–349. [Google Scholar] [CrossRef] [Green Version]
  35. Kousksou, T.; Allouhi, A.; Belattar, M.; Jamil, A.; El Rhafiki, T.; Arid, A.; Zeraouli, Y. Renewable energy potential and national policy directions for sustainable development in Morocco. Renew. Sustain. Energy Rev. 2015, 47, 46–57. [Google Scholar] [CrossRef]
  36. Ministry of Energy, Water and Environment, Department of the Environment. Oasis in Morocco: Adaptation to Climate Change. Resilient Oasis: Adaptation to Climate Change; 2008. [Google Scholar]
  37. Goldemberg, J.; Coelho, S.T. Renewable energy—Traditional biomass vs. modern biomass. Energy Policy 2004, 32, 711–714. [Google Scholar] [CrossRef]
  38. Morocco HPCI. Available online: https://www.hpci-office.jp/cgi-bin/hpcidatabase/app/summary.cgi?view=list&col=f70&date_from=&date_to=&doclang=all&invisiblecol=0&outerhpci=1&prjcls=field&resrc=65551&row=C&s=hp160268&select=32&sop=and&lang=en (accessed on 16 December 2020).
  39. El-Katiri, L. OCP Policy Center; 2016. [Google Scholar]
  40. Gonçalves, H. Solar Heating and Cooling in Portugal. In Solar Heating and Cooling Programme-International Energy Agency, Energy Day; 2018. [Google Scholar]
  41. REN 21. Renewables 2019–Global Status Report; 2019. [Google Scholar]
  42. Bennouna, A. Trois Décennies D’énergie Solaire Thermique au Maroc; EcoActu: Casablanca, Morocco, 2019. [Google Scholar]
  43. Bennouna, A. Trois Décennies D’énergie Solaire Photovoltaïque au Maroc; EcoActu: Casablanca, Morocco, 2019. [Google Scholar]
  44. Moroccan Agency of Energy Efficiency A.M.E.E. Available online: http://www.amee.ma (accessed on 23 December 2020).
  45. TRNSYS Studio Documentation, Mathematical. Available online: https://sel.me.wisc.edu/trnsys/ (accessed on 16 December 2020).
Figure 1. Evolution of the total installed area per inhabitant in Morocco compared to the world.
Figure 1. Evolution of the total installed area per inhabitant in Morocco compared to the world.
Inventions 06 00002 g001
Figure 2. Evolution of the accumulation of the installed area in Morocco (% from the total in world).
Figure 2. Evolution of the accumulation of the installed area in Morocco (% from the total in world).
Inventions 06 00002 g002
Figure 3. Evolution of the unit wholesale price of thermal collectors (MAD per m2) [42].
Figure 3. Evolution of the unit wholesale price of thermal collectors (MAD per m2) [42].
Inventions 06 00002 g003
Figure 4. Evolution of the solar thermal collector’s cumulative area installed annually in Morocco vs. the national market (imported: installed and not installed).
Figure 4. Evolution of the solar thermal collector’s cumulative area installed annually in Morocco vs. the national market (imported: installed and not installed).
Inventions 06 00002 g004
Figure 5. Turnover evolution (solar thermal Moroccan market) for 3500 MAD/m2. [42].
Figure 5. Turnover evolution (solar thermal Moroccan market) for 3500 MAD/m2. [42].
Inventions 06 00002 g005
Figure 6. Domestic solar water heaters sales in Morocco.
Figure 6. Domestic solar water heaters sales in Morocco.
Inventions 06 00002 g006
Figure 7. Methodology.
Figure 7. Methodology.
Inventions 06 00002 g007
Figure 8. Thermosyphon system schematic.
Figure 8. Thermosyphon system schematic.
Inventions 06 00002 g008
Figure 9. Weather data.
Figure 9. Weather data.
Inventions 06 00002 g009
Figure 10. Realistic Moroccan load profile for hot water consumption.
Figure 10. Realistic Moroccan load profile for hot water consumption.
Inventions 06 00002 g010
Figure 11. Thermosyphon diagram executed under the TRNSYS program.
Figure 11. Thermosyphon diagram executed under the TRNSYS program.
Inventions 06 00002 g011
Figure 12. Solar fraction of Michaelides study versus validation study.
Figure 12. Solar fraction of Michaelides study versus validation study.
Inventions 06 00002 g012
Figure 13. (a) Solar fraction. (b) Collector’s efficiency.
Figure 13. (a) Solar fraction. (b) Collector’s efficiency.
Inventions 06 00002 g013
Figure 14. (a) Total annual gains (with state subsidies). (b) Total annual gains (without state subsidies).
Figure 14. (a) Total annual gains (with state subsidies). (b) Total annual gains (without state subsidies).
Inventions 06 00002 g014
Figure 15. (a) Cumulative gains for buildings (non-subsidies scenario). (b) Cumulative gains for modern Moroccan houses (non-subsidies scenario). (c) Cumulative gains for villas (non-subsidies scenario).
Figure 15. (a) Cumulative gains for buildings (non-subsidies scenario). (b) Cumulative gains for modern Moroccan houses (non-subsidies scenario). (c) Cumulative gains for villas (non-subsidies scenario).
Inventions 06 00002 g015
Figure 16. Cumulative area of the SWH installed in Morocco since 2000 and expectations.
Figure 16. Cumulative area of the SWH installed in Morocco since 2000 and expectations.
Inventions 06 00002 g016
Table 1. Optimal tilt angle.
Table 1. Optimal tilt angle.
ZoneAgadirMarrakechErrachidiaFESIfraneTangier
Optimal angle30°31°31°32°30°32°
Table 2. Thermosyphon solar water heater characteristics.
Table 2. Thermosyphon solar water heater characteristics.
ParameterValueUnit
Collector area2.2m2
Intercept efficiency0.74-
Efficiency slope12.33kJ·h−1·m−2·K−1
Tested flow rate30kg·h−1·m−2
Collector slope30–32°
Number of parallel collector risers10-
Riser diameter0.2mm
Header diameter0.2mm
Header length1.0mm
Collector inlet to outlet distance1.273mm
Collector inlet to tank outlet distance0.25mm
Collector inlet diameter0.2mm
Length of collector inlet1.0mm
Collector outlet diameter0.2mm
Length of collector outlet1.0mm
Outlet pipe losses coefficient15kJ·h−1·m−2·K−1
Tank volume202L
Fluid specific heat4.190kJ·kg−1·K−1
Fluid density1000.0kg·m−3
Tank configurationHORIZENTAL-
Maximum heating rate2kW
Maximum pressure8Bars
Thermal losses of the tank1.8kW·h−1·m−2
Capacity of solar circuit5L
Table 3. General Information of the climatic zones under investigation.
Table 3. General Information of the climatic zones under investigation.
Climatic ZoneRepresentative CityAltitude (m)CoordinatesClimateTemperature °C
Z1AGADIR3130°25′ N 9°36′ WSubtropical-semiarid14.3–23.2
Z2TANGER2035°46′ N 5°48′ WMediterranean hot12.6–24.2
Z3FES40334°03′ N 4°58′ WMediterranean/continental9.6–26
Z4IFRANE201933°32′ N 4°58′ WHumid temperate climate6.2–25.4
Z5MARRAKECH45731°37′ N 8°00′ WSemiarid12.4–28
Z6ER-RACHIDIA14131°55′ N 4°25′ WHot desert9.1–33.2
Table 4. Consumption profile of the validation study.
Table 4. Consumption profile of the validation study.
Profile CodeDaily HW Consumption L/PersonTemperature (°C)Pattern
DOM3050Michaelides
Table 5. Technical parameters of the validation study.
Table 5. Technical parameters of the validation study.
ParameterValue
Ac2.72 m2
FR0.791
FRUL24 kJ·h−1·m−2·k−1
Gtest96 kg·h−1·m−2
β40°
NR20
Dr15 mm
di,do22 mm, 22 mm
Hc1000 mm
HO1150 mm
Hr250 mm
Ha300 mm
Hth370 mm
Li, Lo2000 mm, 520 mm
Uapi5.1 kJ·h−1·k−1
UAPo5.1 kJ·h−1·k−1
V162 L
(UA)s6.45 kJ·h−1·k−1
Qaux3 kW
Table 6. Distribution of houses in Morocco according to the 2014 census.
Table 6. Distribution of houses in Morocco according to the 2014 census.
TypePercentage
Housing occupied (habitat)79%
Secondary housing6%
Vacant9%
Professional use3%
Table 7. Distribution of occupied houses per climatic zone in Morocco.
Table 7. Distribution of occupied houses per climatic zone in Morocco.
ZonePercentage of Houses %Number of Houses
1462,109,220
212.5572,790
329.31,340,970
43.1144,960
54.1189,740
65224,320
Table 8. The number of thermosyphon solar water heaters (TSWH) to be installed.
Table 8. The number of thermosyphon solar water heaters (TSWH) to be installed.
Size150 L200 L300 L
ZoneNumberArea m2NumberArea m2NumberArea m2
Zone 1527,4001,160,1001,413,2003,109,00084,400371,300
Zone 2143,200315,1003,837,800844,30023,000100,800
Zone 3335,300737,600898,5001,976,60053,700236,000
Zone 436,30079,80097,200213,700580025,600
Zone 547,500104,400127,200279,700760033,400
Zone 656,000123,400150,300330,700890039,500
Table 9. The total gains for different type of houses according to the two scenarios (with and without subsidies).
Table 9. The total gains for different type of houses according to the two scenarios (with and without subsidies).
Type of HouseApartmentsModern and Traditional Moroccan HousesVillas
Total Investment1145.54144.42586.4
Subsidies
Annual expenditure of gas Water heater (USD)115.473009.4536.95
Gains (USD)71.591865.8622.91
No subsidies
Annual expenditure of gas Water heater (USD)357.40957.83114.36
Gains (USD)221.59593.8670.90
The annual gains (the number of gas bottles)17.05 million45.68 million5.46 million
Table 10. The annual CO2 emissions avoided by the integration and the generalization of the TSWH.
Table 10. The annual CO2 emissions avoided by the integration and the generalization of the TSWH.
ZoneGains (Number of Gas Cylinder) (million)Gains (million m3 of Butane Gas)CO2 Emissions Avoided (tones/year)
Zone144,610.75220,3771,856,000
Zone 28.3641.29446,900
Zone 316.5181.56891,800
Zone 41.597.8685,300
Zone 53.114.83160,200
Zone 63.561.756203,700
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gargab, F.Z.; Allouhi, A.; Kousksou, T.; El-Houari, H.; Jamil, A.; Benbassou, A. A New Project for a Much More Diverse Moroccan Strategic Version: The Generalization of Solar Water Heater. Inventions 2021, 6, 2. https://0-doi-org.brum.beds.ac.uk/10.3390/inventions6010002

AMA Style

Gargab FZ, Allouhi A, Kousksou T, El-Houari H, Jamil A, Benbassou A. A New Project for a Much More Diverse Moroccan Strategic Version: The Generalization of Solar Water Heater. Inventions. 2021; 6(1):2. https://0-doi-org.brum.beds.ac.uk/10.3390/inventions6010002

Chicago/Turabian Style

Gargab, Fatima Zohra, Amine Allouhi, Tarik Kousksou, Haytham El-Houari, Abdelmajid Jamil, and Ali Benbassou. 2021. "A New Project for a Much More Diverse Moroccan Strategic Version: The Generalization of Solar Water Heater" Inventions 6, no. 1: 2. https://0-doi-org.brum.beds.ac.uk/10.3390/inventions6010002

Article Metrics

Back to TopTop