Design, Energy, Environmental and Cost Analysis of an Integrated Collector Storage Solar Water Heater Based on Multi-Criteria Methodology

: The paper presents a design and operation analysis of an Integrated Collector Storage (ICS) solar water heater, which consists of an asymmetric Compound Parabolic Concentrating (CPC) reﬂector trough, while the water tank comprises two concentric cylinders. The annulus between these vessels is partially depressurized and contains a small amount of water in the bottom of the outer vessel which dominantly contributes to the heat transfer from the outer to the inner cylinder. A multi-criteria optimization algorithm is applied to re-evaluate the design speciﬁcations of the parabolic surface, thus modifying the design of the entire ICS system and predict the necessary number of units for achieving the highest possible effectiveness with minimized fabrication costs and environmental impacts. The environmental footprint of the device is assessed through Life Cycle Assessment (LCA). The produced thermal energy in conjunction with the environmental and economic results are evaluated as a function of different conﬁguration parameters regarding the water storage conditions, the solar radiation and the total pressure inside the annulus. The ultimate aim of the evaluation process is to offer new perspectives on the design principles of environmentally friendly and cost-effective devices with improved thermal performance.


Introduction
Solar Water Heating (SWH) systems contribute to the reduction of greenhouse emissions during lifetime operation. Nowadays, different types of SWH systems are used, such as flat plate and evacuated tube thermosiphonic units (Flat Plate Thermosiphonic Units: FPTU; Evacuated Tube Thermosiphonic Units: ETTU) and forced circulating systems (Flat Plate Forced Circulating Units: FPFCU; Evacuated Tube Forced Circulating Units: ETFCU) and Compound Parabolic Concentrating (CPC) reflectors in connection with storage tanks for domestic applications. Usually, CPC reflectors are adopted in Integrated Collector Storage (ICS) water heating systems [1,2]. Solar water heaters are considered "clean" and low-cost thermal energy systems. However, the cost of specific types of solar water heaters with high construction complexity has increased rapidly, i.e., the combination of evacuated tube collectors with CPC-type reflectors could be such a reason for increased cost systems [3]. Souliotis et al. [4] and Singh et al. [5] analysed recent developments in the ICS solar systems and Papaefthimiou et al. [6] presented the developments in energy modelling and management of solar water heating systems during the last decade.
Several scientific studies have presented results on the environmental viability of such systems. Crawford et al. [7] focused on energy consumption and CO 2 emissions during the operation of FPTU water heaters, while Kalogirou [8] presented the environmental and economic analysis of a thermosiphon SWH used for domestic applications in Cyprus with two types of backups (electricity and diesel). Focusing on SWH systems, Yin Hang et al. [9] compared their economic and environmental performance with those conventional for typical residential buildings in the United States, while Colle et al. [10] studied both SWH and PV systems.
Several studies have employed Life Cycle Assessment (LCA) methodology to evaluate the environmental profile of solar thermal systems throughout their whole life cycle. The energy and environmental balance throughout the entire life cycle of flat plate solar collectors were initially presented by Ardente et al. [11], while Crawford and Treloar [12] estimated the embodied energy, the greenhouse gas (GHG) emissions, and the payback period during the whole life cycle of installed solar systems in two Australian cities. The environmental impact assessment of the entire life cycle of an innovative prototype vacuum tube solar collector using Phase Change Material (PCM) was presented by Jachura and Sekret [13]. Furthermore, Karasu and Dincer [14] used a cradle-to-grave LCA to analyse integrated borehole type thermal energy storage (TES) systems in Canadian buildings. Lamnatou et al. in [15][16][17] presented a comprehensive LCA studies analysing the environmental profile of solar devices, either stand alone or Building Integrated Solar Thermal (BIST) systems. The results revealed a high dependence on the materials used for the concentrator/storage tank, and the need for improved recycling processes.
The optimal design of SWH systems aims at low-cost systems with improved thermal behaviour and environmental footprint. An optimization process can focus on the design or/and operation of solar thermal systems in various applications. In this direction, Hobson and Norton [18] proposed a design algorithm for direct thermosiphon SWHs by presenting a characteristic curve for domestic SWHs in order to predict the annual solar fraction for a validated numerical simulation model. On the other hand, Gordon and Rabl [19] proposed an optimization in the design of industrial solar plants without heat storage, with a focus on the lowest fabrication cost. Frei and Vogelsanger [20] presented an optimum design methodology for solar thermal systems for domestic hot water and space heating, including comparative results to more traditional designs. In the work of Lund and Peltola [21], the predesign and the optimization of solar heating systems with seasonable storage were studied, demonstrating a simulation design tool which evaluates the energy system's size and dimensions. In addition, Abdel-Malek [22] described the optimum design of SWHs through the analysis of the design parameters depending on the solar incidence for various areas and the demand for hot water. De Winter [23] presented optimum designs for SWH with different tank designs for a single-family home.
The majority of existing literature regarding the optimization of solar thermal systems does not take into account the anticipated environmental footprint of the proposed systems. However, Ko [24] presented an optimization method for the design of a SWH system on the basis of the Life Cycle Cost (LCC) methodology. A genetic algorithm was employed to optimize the configuration and sizing of three constraint conditions, including energy balance, solar fraction and available space to install solar collectors. Furthermore, Patrčević et al. [25] presented a methodology for conducting a cost-optimal energy performance calculation of a SWH system, used for space heating and domestic hot water needs. The calculation was based on dynamic hourly methods dealing with storage-tank water temperature calculations and provided proposals for modifications in order to overcome the observed ambiguities and shortcomings. The present work presents a detailed study of an ICS system comprising two concentric cylindrical tanks placed horizontally and an asymmetric CPC reflector. The primary objective of the paper is to identify the optimal sizing of the ICS system and to evaluate the optimum dimensions of each part of the reflector. The proposed optimization technique aims at fulfilling four basic criteria: (i). the highest mean daily efficiency, (ii). the lowest water tank's thermal losses at night, (iii). the minimization of the fabrication cost, and (iv). the reduction of environmental impacts throughout the entire life cycle of the device, for different geographic locations of Greece. The second objective of this work is to provide a detailed and comparative energy, cost, and environmental analysis in various operational modes for the proposed ICS units. The economic profile of the models is studied using the Net Present Value (NPV) and the simple payback period, in order to estimate the economic viability of such systems, especially when installed on islands.
The adopted optimization process is not only novel but also offers new perspectives in the design principles of environmentally friendly and cost-effective devices. Parameters such as the size of the reflector (and its geometry) and the required number of units for optimum operation, are all taken into consideration in the pursuit of achieving the combined targets of minimum environmental footprint and maximum performance. To the best of our knowledge, the existing literature on domestic SWH systems mainly focuses on either specific design aspects of devices (noting that the optical analysis is not gaining the interest of most researchers) or on the LCA analysis. The combination of both these evaluation aspects in a single optimization procedure has not been previously reported.
The layout of the paper is as follows. Section 2 presents all technical details about the design, the operational concept, the environmental and economic profile of the of the proposed ICS system. Additionally, the multi-criteria design problem is presented alongside the proposed optimization technique. Section 3 contains the optimum dimensions for the parabolic parts and the comparison of the results, while in the last section, the conclusions of the study are discussed.

Design and Operation Concept of the ICS System
Souliotis et al. [26,27] proposed the reference ICS model, illustrated in Figure 1. The water storage consists of two cylindrical vessels, one inside the other. In particular, the annulus between the two cylinders is partially depressurized containing a small amount of water in the bottom (4 at 15 • C) which acts as PCM. PCM materials can improve the thermal behaviour of solar devices by reducing thermal losses [28][29][30][31][32][33][34]. The vapor that is produced creates a heat diode transfer mechanism from the outer vessel to the inner water storage tank ( Figure 1).

Environmental and Energy Analysis during Fabrication and Installation Phases
The ICS device ( Figure 1) is mainly composed of five parts: the casing, the insulation, the asymmetric parabolic reflector, the solar glass cover, and two concentric water tanks. The design and construction details of a similar system are presented in previous work by Arnaoutakis et al. [35]. Table 1 presents all quantities (in kg) of the employed materials for the system in Figure 1.

Environmental and Energy Analysis during Fabrication and Installation Phases
The ICS device ( Figure 1) is mainly composed of five parts: the casing, the insulation, the asymmetric parabolic reflector, the solar glass cover, and two concentric water tanks. The design and construction details of a similar system are presented in previous work by Arnaoutakis et al. [35]. Table 1 presents all quantities (in kg) of the employed materials for the system in Figure 1.

Environmental and Energy Analysis during Fabrication and Installation Phases
The ICS device ( Figure 1) is mainly composed of five parts: the casing, the insulation, the asymmetric parabolic reflector, the solar glass cover, and two concentric water tanks. The design and construction details of a similar system are presented in previous work by Arnaoutakis et al. [35]. Table 1 presents all quantities (in kg) of the employed materials for the system in Figure 1. The environmental and energy analysis of the studied system is investigated through the LCA method, which follows the ISO standards 14040, 14044 [36,37] and is implemented with SimaPro 8.5 [38] and ecoinvent 3.4 database [39]. The environmental impact assessment methodology used is Recipe Endpoint 2016 Hierarchist [40]. Figure 3 summarizes the environmental footprint of the studied system. The environmental and energy analysis of the studied system is investigated through the LCA method, which follows the ISO standards 14040, 14044 [36,37] and is implemented with SimaPro 8.5 [38] and ecoinvent 3.4 database [39]. The environmental impact assessment methodology used is Recipe Endpoint 2016 Hierarchist [40]. Figure 3 summarizes the environmental footprint of the studied system. As presented in Figure 3, the highest environmental impacts are caused by the extensive use of steel for the construction of the outer cover and the water tank and are mainly assigned to the following impact categories (with decreasing magnitude): climate change, human health, climate change ecosystems, particulate matter formation and human toxicity [39]. The first two categories increase the greenhouse gas emissions that affect climate change, while the last two contribute to potential lung and respiratory problems in human health [39]. Similarly, the necessity for using large quantities of polyurethane (PUR) during the fabrication process mainly influences the climate change category, while the generated wastes from the formulation of the thermal insulation inside the casing affect the particulate matter formation and human toxicity categories [39]. The rest of the materials are assigned with a relatively low impact score. The total energy demands during the fab- As presented in Figure 3, the highest environmental impacts are caused by the extensive use of steel for the construction of the outer cover and the water tank and are mainly assigned to the following impact categories (with decreasing magnitude): climate change, human health, climate change ecosystems, particulate matter formation and human toxicity [39]. The first two categories increase the greenhouse gas emissions that affect climate change, while the last two contribute to potential lung and respiratory problems in human health [39]. Similarly, the necessity for using large quantities of polyurethane (PUR) during the fabrication process mainly influences the climate change category, while the generated wastes from the formulation of the thermal insulation inside the casing affect the particulate matter formation and human toxicity categories [39]. The rest of the materials are assigned with a relatively low impact score. The total energy demands during the fabrication and installation of the ICS system is 43.83 MJ (~12.2 kWh), focusing on the collector with 10.08 MJ (~2.8 kWh), the water tank with 27.73 MJ (~7.7 kWh), the support system with 5.5 MJ (~1.5 kWh) and the installation with 0.52 MJ (~0.15 kWh) [35]. All the above form the production cost during the fabrication and installation phases, with a total value of approximately EUR 292, the production electricity cost EUR 2, the material cost EUR 190 and the salary cost EUR 100 [35].

Energy, Environmental and Economic Analysis in Operation Mode
The energy analysis, throughout the operation phase, can provide the produced thermal quantity of the SWH during a yearly operation period [41]. The annual produced thermal energy (Q THE⁄Year ) of the studied system ( Figure 1) was evaluated by the following equation [42]: where:  [26,27], Souliotis and Tripanagnostopoulos [41], Rabl [43] and Arnaoutakis et al. [44].
Apart from the energy analysis, the environmental and economic studies can provide information about the impact and the viability of a solar thermal energy system during the entire operation period. To approach and investigate both, RETScreen Expert, a dedicated energy management software [45] is utilised. RETScreen Expert is a software that helps decision-makers to determine the technical and financial feasibility of potential renewable energy, energy efficiency, and cogeneration (combined heat and power) projects. Each project goes through a five-step analysis that includes energy, cost, emission, financial, and sensitivity/risk sub-elements.
The net annual greenhouse emission reduction for the ICS system of Figure 1 is approximated by using the test results from various cases of different FPTU systems in the RETScreen computational tool. All these tests were performed based on the weather data of Greece.
The main economic parameters are the Net Present Value (NPV) and the payback period. In finance, the NPV or the Net Present Worth (NPW) [46] is an indicator of the profitability of a project or product that is calculated by subtracting the Present Values (PV) of cash outflows (including initial cost) from the present values of cash inflows over a period of time [47]. Incoming and outgoing cash flows can also be described as benefit and cost cash flows, respectively [47]. The following equation is the formula for the NPV calculation: where: The payback period is the time required to recover the cost of a given investment and is an important determinant of whether to undertake the project [48]. Typically, the payback period is calculated via the expression:

Design Analysis
The main objective of this subsection is the fabrication of totally new water heaters with higher efficiency without significantly increasing the production cost of the system. The design details of the studied system have already been described in the literature [26,49]. The reflector is composed of two asymmetric parabolic parts ((AB), (C A )) and there is only one involute part, the (BC), respectively. Taking into account the geometry of Figure 1, the optimization problem is described by the following set of equations: A detailed description of a similar design problem was presented by Arnaoutakis et al. [44]. Here, the main difference is the existence of three objective functions instead of two: mean daily efficiency n d (ψ m , ϕ m , B, C), thermal loss coefficient U S (D, F) and production cost. In addition, this similar design problem has been considered for the case of only one ICS unit, with the values of the involved parameters (B, C) being kept fixed [44]. In the present work, the parameters B, C, D, and F involved in Equations (1) and (4) are not constant but vary between lower and upper limits, which were defined by previous experimental studies in ICS systems [26,49]. These parameters, along with the angular parameters ψ m and ϕ m (maximum reflector's angles as indicated in Figure 1) and the parameter z (number of the proposing installed ICS units in the building), represent the decision variables for the optimization problem of Equation (4). The objective is to find the size of the system of Figure 1, along with the number of proposed ICS units, which produces the maximum possible amount of thermal energy, whilst at the same time keeping the production cost to a minimum value. The first two objective functions are similar with the coefficients n d (ψ m , ϕ m , B, C) and U S (D, F), which have been described in Section 2.3. In addition, the quantity (−C P (ψ m , ϕ m )) of the third objective function have already been analysed in the existing literature [43].
Based on the geometry of Figure 1 and the analysis of the studies of Souliotis et al. [26], Arnaoutakis et al. [44] and Tripanagnostopoulos and Souliotis [49], the optimum horizontal and vertical dimensions for the parabolic parts (AB) and (C A ), respectively, are expressed using the following formulas: Similarly, [26,44,[49][50][51] the dimensions (lengths) for the parabolic (AB), (C A ) and the involute part (BC), respectively, are calculated as: where ψ m , ϕ m are the angular parameters, constituting also two of the decision variables of the optimization problem, R T is the radius of the water storage and ω m is the maximum angle for the involute part (BC) of the reflector. The angle ω m does not belong to the decision variables of the problem and it takes the specific value ω m = 135 • , as provided by the existing literature [26,49]. The Equation (4) are solved using a multi-criteria design analysis method that optimizes a group of nonlinear objective functions and is based on an existing algorithm used in problems with several decision variables [52,53].
The entire optimization procedure is performed with an inhouse FORTRAN and Visual Basic code. The program code is interfaced with a database containing parameters for various materials (e.g., optical and thermal attributes and prices per preference unit) [26,44,[54][55][56][57]. In addition, the database includes weather data of different Greek locations [58] and estimations of the hot water consumption's dependency on the number of occupants in the installation building. The code supports different reflector types (symmetric or asymmetric) and water tank capacities. Other sub-functions of the algorithm are the evaluation of the produced thermal energy using Equation (1), the greenhouse gas emissions in operation mode for the designed ICS system, and the techno-economical procedure, using Equations (2) and (3) for the calculation of the NPV and Payback period parameters. The annual emission reduction of the redesigned ICS devices and the consumption's dependency with the number of occupants in the building are approached using the RETScreen software and various test results of the systems, which are considered to be installed in Greek territory. The materials, which are derived from the database, are classified through the following environmental criteria: The above categorization considers impacts during fabrication and installation phases and is depicted taking into account the environmental footprint of Figure 3 and the Recipe Endpoint 2016 Hierarchist Methodology.
Additional operations of the program refer to the prediction of the optimum number of the solar thermal water units that will be installed in the building and the evaluation of the produced thermal energy and the greenhouse gas emissions during operation phase for all the designed devices. The lower and upper limits, especially for the parameter e in the 7th constraint of the Equation (4), referring to the prediction of the installed units, are based on the following assumption: the water consumption depends on the domestic hot water needs per occupant (50-60 /day). If that energy reaches some very low values, the designed "one" unit cannot cover the hot water needs of that selected by the input menu occupants, and it is necessary to fabricate more than one device which will be installed in the studied residence (declared in the optimization problem as the coefficient e in the upper limit of the 7th constraint). The same limits in the last constraint (which are expressed in MWh/Year) are based on similar assumptions (using RETScreen's program existing case studies for solar thermal systems) such as the parameter f depending on hot water consumption per occupant. In addition, the quantity: Q THE⁄Year (ψ m , ϕ m , B, C, D, F) of the same constraint is influenced by the desired water temperature, which is also related to the number of people that will use the ICS system.
Thus, in the optimization procedure, if the suggested unit is highly efficient, the specific device can produce the appropriate amount of thermal energy for all the occupants in the building and the algorithm does not allow the making of other water heaters, except the one which has already been designated by the code.
Using a Graphical User Interface (GUI environment), the program extracts the optimum reflector's geometry, the environmental impacts along the entire life cycle of the device, the tank's geometry, the mean daily efficiency, the thermal loss coefficient, the production cost of the system, and the appropriate number of the ICS units. The output form of the algorithm supports optical design plans that describe the optimum geometry for the entire solar water heater. These plans can help the fabrication of specific devices in the industry. The flow charts that summarize the above procedures are depicted below (Figures 4 and 5).     Figure 4).
In the above problem, three original "complexes" of ≥ + 1 points each are created. N is the number of decision variables of the problem (in the existing case N = 7), consisting of a feasible starting point. The variables , are common in two of the A similar design approach has already been referred to in the existing literature [42,50], with significant differences in various points that are described analytically in the following paragraphs.
In the above problem, three original "complexes" of K ≥ N + 1 points each are created. N is the number of decision variables of the problem (in the existing case N = 7), consisting of a feasible starting point. The variables ψ m , ϕ m are common in two of the three "complexes" (these two complexes correspond to the functions n d (ψ m , ϕ m , B, C, z) and −C p (ψ m , ϕ m , z)) and the variable z is common in all three "complexes". K − 1 additional point for each "complex", generated from random numbers. Constraints for all the independent variables are given by the following formulas: where i = 1, 2, . . . , N and j = 1, 2, . . . , K − 1 and r ij are random numbers varying between 0-0.4 for Equations (10) and (11), 0.9-1.0 for Equations (12)- (15), and 3-10 for the last formula. In addition, h 1 to h 7 and g 1 to g 7 are the upper and lower limits of the nonlinear inequality constraints of the optimization problem, as presented in Equation (4). Based on the study of Arnaoutakis et al. [44], the feasible starting point was related to the constraints 58.24 • ≤ ψ m < 90 • and 57.95 • ≤ ϕ m < 90 • . In the specific case, there is an additional point (which similarly was related to constraint 1 ≤ z ≤ e). Thus, in order to locate this point, random numbers (between 0-0.4: first two equations, 0.9-1.0: Equations (12)-(15), and 1-10: the last formula) are used for the parameter r ij . The optimum solution may be lower or equal to the starting point, but the final decision will be provided by the optimization methodology of Figure 5. In addition, it is necessary that the random numbers r ij in Equation (16) be in integer form, such as the result of the formula: r ij · h 7 −g 7 10 (Int: the first three letters of the word "Integer"). If the result is decimal, the integer part of the above calculation is selected.
The algorithm evaluates, except for the variables ψ m , ϕ m (optimum reflector's angles and optimum number of the installed ICS units), the parameters B, C, D, and F, which define the mean daily efficiency and the thermal loss coefficient of the studied ICS system.
The selected configuration of the reflector and the water tank is defined by the input menu using the GUI environment (Figure 4), which also contains information about the materials of the system's components, the climate conditions of the geographical area, and the number of people that will use the specific solar thermal system. Next, two choices are available from the program:

•
The energy analysis option, which calculates the produced thermal energy during the entire life cycle of the system and is accompanied by the evaluation of the annual emission reduction during its operation.

•
The economic analysis choice, which can provide information about the profit of the investment.
The total produced thermal energy from all the designed ICS units for a house with four occupants is exhibited (using the program) by the following equation: Year (17) where z is the number of the ICS units in the building and Q THE⁄Year is the produced thermal energy (MWh/Year) for a single designed unit.
The input parameters for the designed ICS systems are the type of the reflector, the capacity of the water tank, the type of the constructive materials, the meteorological data of the installation area, the number of occupants in the building, and the desired temperature of the water storage. The algorithm approaches the annually produced thermal energy of the geographic location in "Athens" that is set as a reference value for the designed ICS units (the experimental study of the prototype ICS device [26] took place in the city of Patra's which belongs in the "Athens" climate zone [58]). Thus, the designed ICS unit will be able to meet the domestic hot water needs of a house with four occupants located in the "Athens" area. Using the above assumption [27], the code considers five cases for the initial total pressure of partial vacuum within the annulus of the water storage, applying them in the existing multi-criteria design analysis methodology:

Results
The main goal of the entire procedure is to find the ICS design with the lowest cost that presents almost the same mean daily efficiency and thermal loss coefficient with conventional FPTU devices, which are competitive systems to the ICS ones. The water storage volume was set at 48.18 [27] for all the redesigned devices and the materials for the reflector and the tank were chosen to be aluminium and steel (materials which are used in Souliotis et al. studies [26,27]), respectively. After the completion of the simulation procedure, the optimum dimensions for the parabolic parts (AB) and (C A ) of the reflector (which were evaluated through the maximum angles ψ m and ϕ m , respectively, as indicated in Figure 1 and Equations (7) and (9)), are depicted in Figure 6. The input parameters for the designed ICS systems are the type of the reflector, the capacity of the water tank, the type of the constructive materials, the meteorological data of the installation area, the number of occupants in the building, and the desired temperature of the water storage. The algorithm approaches the annually produced thermal energy of the geographic location in "Athens" that is set as a reference value for the designed ICS units (the experimental study of the prototype ICS device [26] took place in the city of Patra's which belongs in the "Athens" climate zone [58]). Thus, the designed ICS unit will be able to meet the domestic hot water needs of a house with four occupants located in the "Athens" area. Using the above assumption [27], the code considers five cases for the initial total pressure of partial vacuum within the annulus of the water storage, applying them in the existing multi-criteria design analysis methodology:

Results
The main goal of the entire procedure is to find the ICS design with the lowest cost that presents almost the same mean daily efficiency and thermal loss coefficient with conventional FPTU devices, which are competitive systems to the ICS ones. The water storage volume was set at 48.18ℓ [27] for all the redesigned devices and the materials for the reflector and the tank were chosen to be aluminium and steel (materials which are used in Souliotis et al. studies [26,27]), respectively. After the completion of the simulation procedure, the optimum dimensions for the parabolic parts (AB) and (C'A') of the reflector (which were evaluated through the maximum angles and , respectively, as indicated in Figure 1 and Equations (7) and (9)), are depicted in Figure 6. Next, using CAD (Computer Aided Design) software, two additional rectangular surfaces are evaluated. These geometrical shapes will be located (hypothetically) in the two parts of the reflector in order to formulate the scale for the parabolic part of the CPC surface in the CAD environment (Figures 7 and 8, respectively). It is considered that there is not any modification in the dimensions of the involute part [44], as discussed above for Next, using CAD (Computer Aided Design) software, two additional rectangular surfaces are evaluated. These geometrical shapes will be located (hypothetically) in the two parts of the reflector in order to formulate the scale for the parabolic part of the CPC surface in the CAD environment (Figures 7 and 8, respectively). It is considered that there is not any modification in the dimensions of the involute part [44], as discussed above for the ω m angular parameter.  The above methodology (applied in industry) can produce different reflectors in a  The above methodology (applied in industry) can produce different reflectors in a very short time with larger or smaller parabolic parts, compared with a prototype parabolic shape which was derived from the studies by Souliotis and Tripanagnostopoulos  The above methodology (applied in industry) can produce different reflectors in a very short time with larger or smaller parabolic parts, compared with a prototype parabolic shape which was derived from the studies by Souliotis and Tripanagnostopoulos [26,27]. The derived designed drawings by this method can cooperate with a Computerized Numerical Control (CNC) machine in order to cut and bend the metal sheeting, which will be used in the fabrication of the studied reflector.
The aggregate results (of the programming code for a single ICS unit) for the mean daily efficiency and the thermal loss coefficient are depicted in Figures 9 and 10, respectively, compared with similar configurations (for the initial total pressure of partial vacuum within the annulus of the water storage) of the reference ICS model [27]. The aggregate results (of the programming code for a single ICS unit) for the mean daily efficiency and the thermal loss coefficient are depicted in Figures 9 and 10, respectively, compared with similar configurations (for the initial total pressure of partial vacuum within the annulus of the water storage) of the reference ICS model [27].  In both Figures 9 and 10, the solid lines represent the variation of the mean daily efficiency and thermal losses coefficient for the original models, while the dash lines represent the corresponding variations for the optimized models. It is important to note that Figure 9. Variation of the mean daily efficiency for the optimized and reference ICS models.

2, 15, x FOR PEER REVIEW 15 of 21
The aggregate results (of the programming code for a single ICS unit) for the mean daily efficiency and the thermal loss coefficient are depicted in Figures 9 and 10, respectively, compared with similar configurations (for the initial total pressure of partial vacuum within the annulus of the water storage) of the reference ICS model [27].  In both Figures 9 and 10, the solid lines represent the variation of the mean daily efficiency and thermal losses coefficient for the original models, while the dash lines represent the corresponding variations for the optimized models. It is important to note that In both Figures 9 and 10, the solid lines represent the variation of the mean daily efficiency and thermal losses coefficient for the original models, while the dash lines represent the corresponding variations for the optimized models. It is important to note that the variations in both mean daily efficiencies and thermal losses coefficients for the original models, depicted with the solid lines, were obtained based on experimental test results. The experimental results in both mean daily efficiencies and thermal losses coefficients are additionally pointed out with several markers in the above figures. Figure 9 shows that the systems with initial pressure "670 mbar" and "990 mbar" perform better than the rest of them. The systems with the lowest values in the mean daily efficiency and the thermal loss coefficient are the "245 mbar" and "86 mbar", respectively. Compared with the original models, as depicted in Figures 9 and 10, the designed systems appear to have higher mean daily efficiency and a slightly lower thermal loss coefficient. For all the above, it can be considered that ICS systems of initial pressure "670 mbar" and "990 mbar" approach the optimal thermal behaviour, thus taking into account the daily and the night operation of them. These findings are consistent with the results obtained by the experimental study of the original ICS devices [26,27]. Table 2 presents the fabrication costs of the five cases studied through the algorithm. The specific table summarizes the fabrication cost that was derived by the code for each separate case, the fabrication cost of the reference ICS model [35], and the differences between the optimized cases and the original ones. It can be seen from Table 2 that the "990 mbar" case exhibits the minimum difference (between the cost of the optimized ICS unit and the cost of the original system), while the "86 mbar" depicts the maximum difference. The "670 mbar" case presents the second minimum difference after the "990 mbar" option. Since the method simulates the domestic hot water needs of a house with four occupants, the most suitable cases correspond to five ICS systems "670 mbar" and "990 mbar" with fabrication costs of EUR 210.74 and 208.98 per unit, respectively. When taking into account the analysis of the previous paragraphs and the information provided from Table 2, the "990" and "670 mbar" cases represent more accurately the optimum solution of the problem, because they are a compromise between the maximum efficiency and the variation in the fabrication cost for the designed ICS unit. The above cases mostly approximate the optimum solution of Equation (4) thus, they will be selected.
Next, the energy and cost analysis (for systems installed in a house with four occupants), including the possible number of the designed units for those selected by the algorithm's best cases (670 mbar and 990) will be implemented. The redesigned cases 670 mbar and 990 mbar are renamed to ICS 670 and ICS 990, respectively, while all cases for the original models are renamed to ICS original. A summary of the results for the redesigned systems, depicting information about the optimum dimensions of the parabolic concentrator, the number of ICS units which will be installed in the building, the produced annual thermal energy, the annual emission reduction, and the production cost of the designed systems, are presented in Table 3. The water storage volume, for the proposed ICS models was considered to be 48.18 , following the design construction details [26,27,49]. In addition, in order to validate the dependency of results with the variations in climate conditions, three different locations in southern, central, and northern Greece were considered (Chania, Athens, and Thessaloniki, respectively). For simplicity, the columns of Table 3 are classified using the keywords C: Chania, A: Athens, T: Thessaloniki. In Table 4, the environmental impacts in the production and installation phase are presented. Both redesigned best cases systems behave the same for all components.  Table 4. Environmental impacts in fabrication and installation phases for the two selected, by the algorithm, best cases, following the description of the + and-rating sequence in Section 2.4.

Outer Cover Pollutants
According to Table 3, the designed models try to reach the producing energy of the "Athens" area [27], thus modifying the reflector's geometry. Based on that assumption, the program decreases the size of the reflector, the production cost, and the number of installed units in the building that will be installed in southern locations. On the other hand, the algorithm predicts more units for northern locations with a larger reflector width and therefore higher production costs. The next step is the evaluation of the profitability for one unit of these studied cases, considering only the geographic location of Chania (Crete), which appears to be the most cost-effective and energy-efficient proposed solution by the algorithm, which will be attempted using the economic analysis procedure of the GUI environment. Table 5 presents the economical results for the redesigned systems. The aggregate results for the produced thermal energy and the net annual emission reduction for the original and best-redesigned two models are presented in Figure 11. The results refer to one unit for each model. It is evident from Figure 11 that system ICS 990 demonstrates very high energy and environmental performance compared to the rest of the ICS systems and constitutes the best choice for the design of specific devices, if taking into account only the maximization of the producing thermal energy. Although the original ICS system presents comparable ranges to the optimum ICS systems (ICS 670 and ICS 990) for both produced thermal energy and net annual emission reduction, the latter operate marginally better. The optimum configuration of these units contributes to relatively low environmental impacts during operation mode, depending on the existing thermal energy production. The system ICS 670 follows and corresponds to the optimum configuration for the designed ICS vessels because it combines energy and environmental efficiency with the minimization of the fabrication cost and the prediction of the minimum number of units. According to Table 3, the algorithm decreases the number of units for systems that will be installed in geographic locations with higher solar radiation, while it increases the corresponding units for systems installed in northern regions. In addition, it should be stressed that the redesigned systems are valuable investments against the original ICS unit, taking into account the NPV and simple payback period values, as presented in Table 5.

Conclusions
A detailed energy, environmental, economic, and design analysis has been presented in two stages for ICS solar thermal systems, focusing on domestic hot water applications in Greece. The prototype ICS unit is a specially designed device with an asymmetric CPC reflector that comprises two concentric cylindrical vessels in the water storage tank. The annulus between the cylinders is partially depressurized containing a small amount of water that serves as PCM material, creating a thermal transfer mechanism from the environment to the inner vessels.
In order to improve the energy, economic and environmental efficiency of the prototype ICS system, a multi-criteria optimization method has been employed. In the first stage, the design of the parabolic parts of the CPC reflector's surface (for varying climate conditions) and the evaluation of the minimum number of the units (in house with four occupants) are implemented. This process maximizes the energy and environmental performance, whilst at the same time minimising the production cost of the studied systems. Based on this optimization, new ICS systems are designed comprising a totally new geometry of the reflector's surface with optimized horizontal and vertical dimensions, as It is evident from Figure 11 that system ICS 990 demonstrates very high energy and environmental performance compared to the rest of the ICS systems and constitutes the best choice for the design of specific devices, if taking into account only the maximization of the producing thermal energy. Although the original ICS system presents comparable ranges to the optimum ICS systems (ICS 670 and ICS 990) for both produced thermal energy and net annual emission reduction, the latter operate marginally better. The optimum configuration of these units contributes to relatively low environmental impacts during operation mode, depending on the existing thermal energy production. The system ICS 670 follows and corresponds to the optimum configuration for the designed ICS vessels because it combines energy and environmental efficiency with the minimization of the fabrication cost and the prediction of the minimum number of units. According to Table 3, the algorithm decreases the number of units for systems that will be installed in geographic locations with higher solar radiation, while it increases the corresponding units for systems installed in northern regions. In addition, it should be stressed that the redesigned systems are valuable investments against the original ICS unit, taking into account the NPV and simple payback period values, as presented in Table 5.

Conclusions
A detailed energy, environmental, economic, and design analysis has been presented in two stages for ICS solar thermal systems, focusing on domestic hot water applications in Greece. The prototype ICS unit is a specially designed device with an asymmetric CPC reflector that comprises two concentric cylindrical vessels in the water storage tank. The annulus between the cylinders is partially depressurized containing a small amount of water that serves as PCM material, creating a thermal transfer mechanism from the environment to the inner vessels.
In order to improve the energy, economic and environmental efficiency of the prototype ICS system, a multi-criteria optimization method has been employed. In the first stage, the design of the parabolic parts of the CPC reflector's surface (for varying climate conditions) and the evaluation of the minimum number of the units (in house with four occupants) are implemented. This process maximizes the energy and environmental performance, whilst at the same time minimising the production cost of the studied systems. Based on this optimization, new ICS systems are designed comprising a totally new geometry of the reflector's surface with optimized horizontal and vertical dimensions, as these are evaluated by the optimum angles ψ m and ϕ m of the parabolic curves. The adopted procedure includes the design of the entire ICS system's components. In the second stage, a methodology for the formulation of the specific reflector's geometry, using the above optimization technique and CAD environment, is proposed. The results indicate that the ICS devices are cost-effective, compared to the prototype ICS units. Several tests depicted the variation of the mean daily efficiency and the thermal loss coefficient of the improved ICS devices, which differ in the initial total pressure inside the annulus. In particular, the "670 mbar" choice produces the minimum number of the designed vessels, which can operate with higher performance compared to the original system and represents the optimum solution of the design analysis.
The conducted energy, environmental and cost analysis for varying climate conditions, combined with the design analysis methodology in the second stage, indicate that the cases "670 and 990 mbar" improve the energy performance of the prototype device. The analysis suggests the installation of a lower number of units for southern regions. Based on the energy and environmental analysis, it is evident that the above cases are much more efficient and environmentally friendly during operation mode. In addition, the economic study proves that the proposed ICS devices exhibit better economical behaviour compared to the original one. The present study might serve as a guidance for the design and selection of improved ICS devices, which will be installed in buildings for domestic applications.