Experimental and Numerical Thermal Analysis of Multi-Layered Microchannel Heat Sink for Concentrating Photovoltaic Application

Concentrating photovoltaic has a major challenge due to the high temperature raised during the process which reduces the efficiency of the solar cell. A multi-layered microchannel heat sink technique is considered more efficient in terms of heat removal and pumping power among many other cooling techniques. Thus, in the current work, multi-layered microchannel heat sink is used for concentrating photovoltaic cooling. The thermal behavior of the system is experimentally and numerically investigated. The results show that in extreme heating load of 30 W/cm2 with heat transfer fluid flow rate of 30 mL/min, increasing the number of layers from one to four reduces the heat source temperature from 88.55 to 73.57 ◦C. In addition, the single layered MLM heat sink suffers from the highest non-uniformity in the heat source temperature compared to the heat sinks with the higher number of layers. Additionally, the results show that increasing the number of layers from one to four reduces the pressure drop from 162.79 to 32.75 Pa.


Introduction
Concentrating photovoltaic (CPV) technology has gained much more attention than the conventional photovoltaic due its competitive cost [1].Eighty percent of the concentrated sunlight is absorbed by the solar cell, only a small portion is converted to electrical energy, and the remaining is converted to thermal energy which causes an increase of the solar cell temperature [2,3].A study showed that the solar cell temperature could reach 1400 • C under 500 concentration ratio if it is fully insulated [4].Studies show that the temperature effect on the solar cell efficiency depends on its type for instance the output power declines in an order of 0.65%/ • C, 0.4%/ • C and 0.248%/ • C for crystalline solar cells, thin-film cell and triple-junction solar cells respectively [5].Additionally, the excessive temperature causes a long-term degradation of the solar cell if it exceeds the maximum limit [6].Temperature distribution on the surface of a solar cell due to the non-uniform illumination on the surface is another parameter that affects the solar cell efficiency [7].This implies the importance of removing the generated heat efficiently.
Solar cells are classified according to arrangement as single-cell, linear and densely packed cells [8].In a single solar cell receiver, the heat is concentrated at the center of the receiver.This implies that the heat needs to be transferred to the base layers with lowest thermal resistance and equally to achieve the reduction of the solar cell maximum temperature and the temperature uniformity.Theristis and O'Donovan [9] found that CPV single solar cells can be cooled using a passive approach having thermal resistance of 1.7 K/W but in extreme outdoor conditions, lower thermal resistance (1.4 K/W) is needed which means an active cooling system is required.On the other hand, Aldossary et al. [10] have examined experimentally and numerically different cooling techniques for CPV in a harsh environment.The simulations show that the round pin and straight fins heat sinks are not sufficient to maintain the single solar cell below the operational temperature limit of 80 • C. Additionally, the study experimentally investigated the thermal behavior of a 1 cm 2 single solar cell receiver under high concentration attached to a water cooled heat sink using a rectangular cooling channel with 10 mm height.The results show that the solar cell maximum surface temperature was maintained at 60 • C for water velocity of 0.01 m/s with a heat transfer coefficient of 1645 W/(m 2 K).
A multi-layered microchannel (MLM) heat sink technique is considered more efficient in terms of pumping power and heat removal capability among many other cooling techniques [11].Vafai and Zhus [12] were the first to propose the concept of multi-layered microchannel heat sink (MLM) in electronic cooling applications.They numerically studied the performance of two layers for a counter flow and a uniform heat source and found, interestingly, that the pressure drop and temperature rise on the base surface were reduced compared to the conventional one layered heat sink.In later years, many studies have been published investigating the multi-layered techniques of different channel shapes, materials, flow and fluids.The optimum number of layers was investigated by Wei and Joshi [13] where the number of layers varied between one and six.The analysis shows a significant effect of the number of layers on the thermal resistance and pumping power, and the optimum number of layers was found to be three.In addition, studies have proposed different materials of layers such as silicon [14][15][16], copper [17,18], aluminum [19], steel and diamond [11].Additionally, the studies have investigated different fluid flow patterns such as parallel flow [20][21][22], counter flow [23][24][25][26] and both flows comparison [27][28][29].Different coolants have been tested for the MLM heat sink applications such as water [30], air [31,32], nanofluids [33,34], deionized water and ethylene glycol [35].The MLM heat sink investigations have also studied the effect of channel shape such as rectangle [36][37][38], square [39,40], triangle [19] and wavy [41].Table 1 summarizes the MLM heat sink studies and their important findings.
The majority of the studies conducted for MLM heat sink have considered a uniform heat source and very few have considered the non-uniform heat source.Wei et al. [42] have compared the behavior of MLM heat sink under uniform and partial heat source numerically and experimentally.The experiments have used several heaters to simulate the partial heating scenario.The results show that the heat sink for partial heating has higher total thermal resistance than the full heating due to localized heating.Additionally, the partial heating effect increases as flow rate increases.The authors have pointed out the importance of the location of the heat source to minimize the localized heating effect.In addition, the flow directions (counter flow and parallel flow) have been explored.The results show that double layered heat sinks experience similar thermal resistance of 0.09 • C/(W/cm 2 ) for both flow patterns.However, the heat sink temperature in the counter flow pattern has 40% less temperature non-uniformity than parallel flow.Recently, Ansari and Kim [43] have explored, numerically, the performance of a double-layered heat sink under non-uniform heating conditions with a random hotspot for microprocessor cooling applications.However, the study showed the effect of the flow pattern without any comparison with a single layered heat sink.• For a single layer counter flow: Heat sink performance (heat transfer and hydrodynamic) in laminar flow can be similar of that in turbulent flow.

•
For a double layer counter flow: Heat sink performance in laminar flow is better than that of the turbulent.

•
Optimum fins width, channel width and aspect ratio of two-layered microchannel heat sink.

•
The optimum number of layers for lower thermal resistance and pumping power.

•
Thermal resistance is inversely proportional to the pumping power.

•
Thermal resistance depends on the channel length.The shorter length, the small thermal resistance.

•
In this study the optimum number of layers is three.

•
To determine the lowest thermal resistance by varying the height sink dimensions such as height and width of channels, height and width of fins.

•
For materials with low thermal conductivity, the thermal resistance increases with increase of layers.

•
The materials thermal conductivity is not important for low coolant flow rate.

•
To determine the effect of the fluid flow on the overall heat transfer performance.

•
The performance of the heat sink is dependent on the flow characteristics in the channels.

•
To investigate the effect of presence of passive microstructure in the channel

•
The increase of the height of passive microstructure compare to microchannel height leads to lower thermal resistance.

•
The comparison of performance of single and multiple layers of minichannel heat sink.

•
The increase of the minichannel layers leads to reduce the thermal resistance and pumping power.

•
Investigation of non-uniform allocation of microchannel for more temperature uniformity and less pumping power.

•
The study was able to achieve 50% less pumping power with sufficient cooling.

•
Analysis of the effect of individual geometric parameters in the performance of the heat sink and to find out its optimum parameters.

•
Optimal thermal resistance decreases with pumping power.

•
As the coolant pumping power increases, the optimal channel number, lower and upper aspect ratio increase and the optimal width ratio decreases.

•
Analysis of enhancement double-layered heat sink by using a nanofluid and geometric parameters variation.

•
Significant thermal performance improvement of heat sink when using the nanofluid compared to water.

•
The heat sink effectiveness declines significantly under high pumping power. [

•
Investigation of the effect of wavy on the heat sink performance.

•
The effects of channel number, aspect ratio and velocity ratio on the overall thermal performance.
Energies 2019, 12, 122 6 of 25 Applying an MLM heat sink in a CPV cooling application requires attention due to the rapid temperature increase in the CPV module and the importance of solar cell temperature uniformity.Few studies have been found that apply this cooling technique in the CPV applications [45,46].Al Siyabi et al. [45,47] have numerically conducted a thermal analysis using an MLM heat sink in a single solar cell receiver.They considered a non-uniform heat source for the MLM heat sink.The results show a significant temperature reduction of the solar cell when the number of layers increased from single layer to four-layers by 15 • C. In addition, the pressure drop decreased by more than 50% when using a four-layers heat sink and reduced the temperature non-uniformity of the solar cell.
As stated above, few studies investigated the MLM heat sink technique in CPV applications.The main objective of this study is to experimentally examine the behavior of the MLM heat sink integrated with CPV assembly.Therefore, experiments have been conducted to evaluate the thermal performance of an MLM heat sink for different numbers of layers, different heat transfer fluid (HTF) flow rates and heating power rates.A heat source of electrical resistance heater (ERH) has been considered to represent the CPV assembly.In addition, a numerical model has been developed to further analyze the heat sink behavior for measurements that cannot be obtained using the experimental approach.

Test Loop
A schematic and a photograph of the experimental setup is shown in Figure 1.The experimental setup consists of a heat sink testing module, HTF circulation bath with a built-in pump (MultiTemp ® III, Pharmacia Biotech AB, Uppsala, Sweden), a variable area glass flowmeter (NFX-1/4"-25-B, Nixon flowmeters, Cheltenham, UK), differential pressure meter (Digitron 2080P, Rototherm group, Swansea, UK), K-type thermocouples, data logger (Keithley 2700, Keithley Instruments, Inc., Cleveland, OH, USA) and a desktop computer.An HTF with constant flow rate and temperature is supplied using the circulation bath at a temperature of 24 • C. The HTF volume flow rate is measured using the variable area flowmeter and controlled using a built-in needle valve.The cold HTF picks up the heat from the MLM heat sink and passes to the circulation bath.K-type thermocouples are placed in the inlet, outlet and within the test module in different locations.The temperature is recorded with an interval of 10 s until the steady state is reached.

Test Module
The test module, illustrated in Figure 2, was designed to stack layers between one and four.The case was made of a transparent acrylic and consists of two pieces, top and bottom.Bolts were used to join the case pieces and to stand any forces.A rubber sealing was placed between the top and bottom parts to prevent any leakage issues.The HTF inlet and outlet ports were at the side of the case using an 8 mm diameter copper connector.The area in-between the inlet and the MLM entrance were kept to develop a uniform flow of the HTF.The HTF inlet and outlet temperatures were measured using a 1 mm K type thermocouple stainless steel, which was inserted just after the inlet and just before the outlet ports.The pressure drop across the MLM heat sink was captured using two pressure ports using the 8 mm copper connectors.A 0.5 mm K type thermocouple stainless steel was inserted in the middle channel of the first layer middle channel to measure the temperature at this location.Other thermocouple locations were distributed in the middle of the top surface of the heat source to capture the temperature uniformity.

Test Module
The test module, illustrated in Figure 2, was designed to stack layers between one and four.The case was made of a transparent acrylic and consists of two pieces, top and bottom.Bolts were used to join the case pieces and to stand any forces.A rubber sealing was placed between the top and bottom parts to prevent any leakage issues.The HTF inlet and outlet ports were at the side of the case using an 8 mm diameter copper connector.The area in-between the inlet and the MLM entrance were kept to develop a uniform flow of the HTF.The HTF inlet and outlet temperatures were measured using a 1 mm K type thermocouple stainless steel, which was inserted just after the inlet and just before the outlet ports.The pressure drop across the MLM heat sink was captured using two pressure ports using the 8 mm copper connectors.A 0.5 mm K type thermocouple stainless steel was inserted in the middle channel of the first layer middle channel to measure the temperature at this location.Other thermocouple locations were distributed in the middle of the top surface of the heat source to capture the temperature uniformity.A CPV assembly (3C44A, AZUR SPACE Solar Power GmbH, Heilbronn, Germany) is used in this study.The assembly consists of a triple junction (TP) solar cell of size 10 mm × 10 mm as shown in Figure 3a.As shown in Figure 3b, the CPV assembly is composed of a solar cell located under a conductive layer made of copper.For the electrical insulation, a ceramic layer was placed under the copper layer.Finally, the full arrangement was placed in a second copper layer for thermal conduction purposes.The assembly also consisted of two by-pass diodes and two electrical terminals.As pointed out, this study was dealing with thermal behavior of the MLM heat sink in cooling the CPV and, for simplicity, the solar cell in the CPV assembly was replaced by the electrical resistance heater (ERH) (PWR263S-35, Bourns, Inc, UK) to simulate the heat load generated by the CPV.Therefore, the solar cell was removed by placing the CPV assembly on a preheated hot plate.Then, the solar cell was pulled and the ERH was placed in the same location.The ERH power rating was 30 W rating of a size of 10 × 10 mm and was operated using a DC power supply.A CPV assembly (3C44A, AZUR SPACE Solar Power GmbH, Heilbronn, Germany) is used in this study.The assembly consists of a triple junction (TP) solar cell of size 10 mm × 10 mm as shown in Figure 3a.As shown in Figure 3b, the CPV assembly is composed of a solar cell located under a conductive layer made of copper.For the electrical insulation, a ceramic layer was placed under the copper layer.Finally, the full arrangement was placed in a second copper layer for thermal conduction purposes.The assembly also consisted of two by-pass diodes and two electrical terminals.As pointed out, this study was dealing with thermal behavior of the MLM heat sink in cooling the CPV and, for simplicity, the solar cell in the CPV assembly was replaced by the electrical resistance heater (ERH) (PWR263S-35, Bourns, Inc., Bedford, UK) to simulate the heat load generated by the CPV.Therefore, the solar cell was removed by placing the CPV assembly on a preheated hot plate.Then, the solar cell was pulled and the ERH was placed in the same location.The ERH power rating was 30 W rating of a size of 10 × 10 mm and was operated using a DC power supply.
The microchannel plates were made from aluminum and the width (W) and length (L) were 32 × 30 mm, respectively.The other microchannel dimensions are shown in Figure 4a and Table 2.
Micro-sawing technique was used to fabricate the rectangular channels machining.Micro-sawing is a low cost microchannel fabrication method and offers a better surface finish compared to other techniques [48].A scan microscope was used to inspect the channels' shape and to ensure that the sides of the channels are straight and within the accuracy range.As can be seen in Figure 4b, the verticality of the walls was properly achieved, and no significant inclination is visible.Moreover, the bottom surfaces of the channels are relatively smooth.A CPV assembly (3C44A, AZUR SPACE Solar Power GmbH, Heilbronn, Germany) is used in this study.The assembly consists of a triple junction (TP) solar cell of size 10 mm × 10 mm as shown in Figure 3a.As shown in Figure 3b, the CPV assembly is composed of a solar cell located under a conductive layer made of copper.For the electrical insulation, a ceramic layer was placed under the copper layer.Finally, the full arrangement was placed in a second copper layer for thermal conduction purposes.The assembly also consisted of two by-pass diodes and two electrical terminals.As pointed out, this study was dealing with thermal behavior of the MLM heat sink in cooling the CPV and, for simplicity, the solar cell in the CPV assembly was replaced by the electrical resistance heater (ERH) (PWR263S-35, Bourns, Inc, UK) to simulate the heat load generated by the CPV.Therefore, the solar cell was removed by placing the CPV assembly on a preheated hot plate.Then, the solar cell was pulled and the ERH was placed in the same location.The ERH power rating was 30 W rating of a size of 10 × 10 mm and was operated using a DC power supply.The microchannel plates were made from aluminum and the width (W) and length (L) were 32 × 30 mm, respectively.The other microchannel dimensions are shown in Figure 4a and Table 2. Microsawing technique was used to fabricate the rectangular channels machining.Micro-sawing is a low cost microchannel fabrication method and offers a better surface finish compared to other techniques [48].A scan microscope was used to inspect the channels' shape and to ensure that the sides of the channels are straight and within the accuracy range.As can be seen in Figure 4b, the verticality of the walls was properly achieved, and no significant inclination is visible.Moreover, the bottom surfaces of the channels are relatively smooth.Table 2.The dimensions of the microchannel.

Test Procedures and Uncertainty Analysis
A typical test run was started when the test module was assembled and all instruments were connected.The modified CPV assembly was attached to the top layer of the heat sink and a thin layer of silicone thermal transfer paste was applied for better thermal conductivity between the two plates.Water was used here as HTF throughout all experiments.Before starting any test, water was circulated from the right to left direction at a temperature of 24 °C and to the desired flow rate for 1 h to ensure the temperature of the heat sink was uniform and at its initial conditions, and also to remove the trapped air bubbles from the loop.Throughout the experiments, the water inlet temperature was allowed to vary between 24.1 and 23.9 °C.The experiments were started by switching on the DC power supply and the data logger.Then, the experiment was monitored for reaching the temperature steady state condition.The steady state was achieved when the temperature variations in all thermocouple measurement were less than 0.1 °C.The pressure drop measurements were recorded manually.
A list of controlling parameters are shown in Table 3.It shows that the number of the stacked layers in the heat sink is varied from one to four layers.For a fixed layer, the heating power is varied in a range between 5 and 30 W and the HTF flow rate varied from 30 to 60 mL/min.

Test Procedures and Uncertainty Analysis
A typical test run was started when the test module was assembled and all instruments were connected.The modified CPV assembly was attached to the top layer of the heat sink and a thin layer of silicone thermal transfer paste was applied for better thermal conductivity between the two plates.Water was used here as HTF throughout all experiments.Before starting any test, water was circulated from the right to left direction at a temperature of 24 • C and to the desired flow rate for 1 h to ensure the temperature of the heat sink was uniform and at its initial conditions, and also to remove the trapped air bubbles from the loop.Throughout the experiments, the water inlet temperature was allowed to vary between 24.1 and 23.9 • C. The experiments were started by switching on the DC power supply and the data logger.Then, the experiment was monitored for reaching the temperature steady state condition.The steady state was achieved when the temperature variations in all thermocouple measurement were less than 0.1 • C. The pressure drop measurements were recorded manually.3. It shows that the number of the stacked layers in the heat sink is varied from one to four layers.For a fixed layer, the heating power is varied in a range between 5 and 30 W and the HTF flow rate varied from 30 to 60 mL/min.The volume flow measurements were verified manually using a conventional glass water column and stop watch method.Table 4 shows the measurements variations between the two methods.In addition, the mass flow rate was measured using a weight scale.The table shows the maximum variation between the flow meter and the water column of 3.33%.Two structures of thermocouples were used; a wired and a prop type thermocouple base on the measuring location.The prop thermocouples type was used to measure the temperature of fluid flowing inside the tubes, the inlet and the outlet temperature at the test section.The tip of the thermocouple was the place of taking measurements and this was ensured to be in full contact with the fluid.The calibrated props type thermocouples for the inline fluid flow temperature measurements in the HTF inlet and outlet temperature in the heat sink and the middle microchannel temperature was allowed to be permanently fixed in its locations for all experiments.
Wire thermocouples were produced from a roll of thermocouple to be cut according to the desired length.The tips of the thermocouples were made on our own using a thermocouple welder using the L60 thermocouple welder.Three important factors affected the measurements accuracy while making the tips of the thermocouples; the tip length, perfect connection and free of any oxide particles.Therefore, the tips were made as short as possible.Visual inspection using an eye magnifier was used to check the tips' perfection and if it was free from any oxide particles.A voltmeter was used to check the connection of the thermocouples.The thermocouples and the data acquisition system were calibrated using two points of temperatures; boiling water and liquid-ice water mixture against a precision mercury thermometer to an uncertainty of ±0.1 • C.
The power dissipation by the ERH was determined from the product of the voltage and current measured at the electrical terminals.The measured voltage and current ranged between 8 and 21.5 volts and current of ranged between 570 and 1450 mA, respectively.Uncertainties for direct voltage measurement were negligible due to the low voltage drop between the measured terminals.In this experiment, voltage input for the heaters were measured across the ERH and the current meter with 0.1 Ω shunt resistance.For the current settings we had, this caused 0.1 volts loss in the actual ERH.For a 21.5 volts case, this results in ±0.46% uncertainty in voltage measurement.As per the product manual, the current measurement has an uncertainty of ±0.1%.These uncertainties caused a ±0.56% in power input measurement.

Numerical Study
The MLM heat sink was taken further to be analyzed numerically for better understanding of the heat sink thermal performance.COMSOL Multiphysics (version 5.2 COMSOL Inc., Stockholm, Energies 2019, 12, 122 10 of 25 Sweden) was used to simulate the model using the equations of conservation of mass, momentum and energy for single-phase laminar flow [49].

Model Description and Mathematical Formulation
The modeled system is illustrated in Figure 5.The CPV solar cell was replaced by the ERH.Table 5 shows the CPV assembly dimensions, thickness and materials [10].The two by-pass diodes and two electrical terminals were not considered in the model for simplicity.

Model Description and Mathematical Formulation
The modeled system is illustrated in Figure 5.The CPV solar cell was replaced by the ERH.Table 5 shows the CPV assembly dimensions, thickness and materials [10].The two by-pass diodes and two electrical terminals were not considered in the model for simplicity.The thermal boundary conditions of the studied model are shown in Figure 6 and listed in Table 6.As can be seen that the top surfaces of the assembly are exposed to the ambient temperature at a heat transfer coefficient of 15 W/m 2 K.The ERH is assumed to be a heat source object.The surfaces at the bottom and sides are assumed to be thermally insulated.Additionally, the following assumptions were considered:


The heat source is uniform.


The ambient temperature is 25 °C.


Water is used as the HTF and its properties varies with temperature [50,51].The water inlet temperature is 24 °C. HTF varies between 30 and 60 mL/min which corresponds to a low Reynolds numbers.The flow is assumed to be steady, laminar and fully developed in each channel.The thermal boundary conditions of the studied model are shown in Figure 6 and listed in Table 6.As can be seen that the top surfaces of the assembly are exposed to the ambient temperature at a heat transfer coefficient of 15 W/m 2 K.The ERH is assumed to be a heat source object.The surfaces at the bottom and sides are assumed to be thermally insulated.Additionally, the following assumptions were considered:

•
The heat source is uniform.

•
The ambient temperature is 25 • C.

•
Water is used as the HTF and its properties varies with temperature [50,51].The water inlet temperature is 24 • C.

•
HTF varies between 30 and 60 mL/min which corresponds to a low Reynolds numbers.The flow is assumed to be steady, laminar and fully developed in each channel.

Model Description and Mathematical Formulation
The modeled system is illustrated in Figure 5.The CPV solar cell was replaced by the ERH.Table 5 shows the CPV assembly dimensions, thickness and materials [10].The two by-pass diodes and two electrical terminals were not considered in the model for simplicity.The thermal boundary conditions of the studied model are shown in Figure 6 and listed in Table 6.As can be seen that the top surfaces of the assembly are exposed to the ambient temperature at a heat transfer coefficient of 15 W/m 2 K.The ERH is assumed to be a heat source object.The surfaces at the bottom and sides are assumed to be thermally insulated.Additionally, the following assumptions were considered:


The heat source is uniform.


The ambient temperature is 25 °C.


Water is used as the HTF and its properties varies with temperature [50,51].The water inlet temperature is 24 °C. HTF varies between 30 and 60 mL/min which corresponds to a low Reynolds numbers.The flow is assumed to be steady, laminar and fully developed in each channel.The simulations were conducted in steady state study conditions.The electrical power (Q ele ) supplied to the ERH was assumed to be converted to heat and was calculated by [10]: All the top surfaces of the copper plate and the heat source surfaces release heat to the outer environment using natural convection mode and its heat transfer rate (q conv ) is represented by: where h is the convection heat transfer coefficient (W/m 2 K), A is the surface area exposed to the outer environment (m 2 ) and ∆T is the temperature difference between the object surface and the outer environment (K).In addition, the heat is transferred to the environment using radiation mode and its transfer rate (W) is expressed by [52]: where ε is the emissivity of the materials, σ is Stefan-Boltzmann constant, A is the radiated surface area, T surf is the surface area and T amb is the ambient temperature.The coolant fluid flows through the channels at mass flow rate (kg/s).
The governing equations for this 3-D conjugated heat transfer are as following [33,34]: Continuity equation for the HTF: where u, v and w are the velocity components in the x, y and z directions, respectively.Momentum equation for the HTF: where p is the HTF pressure, ρ f and µ f are the density and dynamic viscosity of the HTF, respectively, and velocity components in the x, y and z directions, respectively.The energy equation for the HTF: where T f , C p, f and k f are the temperature, specific heat and thermal conductivity of the HTF, respectively.The energy equation for the solid domain: where T s and k s are the temperature and thermal conductivity of the solid domain.The associated boundary conditions for the governing equations for the channels in all surfaces as follows: HTF inlet: HTF outlet: HTF-solid interface: The correlations of the HTF density (ρ f ), specific heat C p, f , thermal conductivity k f , and dynamic viscosity (u f ), as functions of HTF temperature are described as follows [43,51]: where, T f is the temperature in Kelvin (K).
As mentioned earlier, the electrical performance of solar cell is dependent on its temperature.Therefore, the temperature measurement of heat source (ERH) represents an important parameter.The total amount of heat absorbed by water (Q out ) in the MLM heat sink is found by .
m is the HTF mass flow rate which is found by V, C p is the specific heat of the HTF, T in and T out are water inlet and outlet temperature, respectively.Therefore, the total heat lost to the environment Q loss is found by Thermal resistance of the heat sink is another important parameter to quantify the thermal performance of the heat sinks.The total thermal resistance (R tot ) of the heat sink is the ratio between the maximum temperature difference and the heat source flux and calculated by: where ∆T max is the temperature difference between the maximum surface temperature of the ERH (T s,max ) and the temperature of water at the inlet T f ,in q is the heat flux per unit area, A s is the heat sink surface area.The temperature uniformity within the heat source is measured using the on-chip thermal resistance R chip and is calculated by [13]: where T j,max and T j,min are the maximum and the minimum temperatures within the heat source, respectively.Water pressure drop across the heat sink is another criteria to evaluate the heat sink in terms of the fluid flow as this is used to determine the required power to pump the coolant [41].
The pumping power (P) is calculated by V is the volume flow rate, ∆p is the difference between the inlet pressure p in and the outlet pressure p out and η p is the pump efficiency.The Reynolds number is used to characterize the flow inside the microchannel and is calculated by where ρ and µ are the density and viscosity of the fluid respectively and are evaluated at the mean water temperature [17].V avg is the average velocity of the fluid and is calculated by where, N c is the total number of channels in the heat sink, and A c is the microchannel cross-sectional area.The average heat transfer coefficient by convection is calculated by [19]: The Nusselt number (Nu) is calculated by [53]: where D h is the hydraulic diameter and k is the fluid thermal conductivity.

Grid Independence Study
The computational domain was meshed using a free tetrahedral grid system.Different mesh sizes were used to the model in order to confirm that the solution is independent of the size of the mesh for each case.Table 7 shows the variation of the maximum heat source temperature and HTF outlet temperature with the number of elements (NoE).The number of elements was varied between 560,765 and 1,669,580.The maximum heat source temperature remained unchanged for the number of elements between 1,510,497 and 1,669,580.On the other hand, it no change of the HTF outlet temperature was noticed when using the NoE between 879,992 and 1,669,580.Therefore, the mesh of NoE with 1,510,497 was selected for numerical analysis.

Model Validation
The CFD model was validated using experimental results of the same parameters.The HTF outlet temperature and the resistor temperature for the heat source was used for the validation for all the layers arrangements and of heat powers ranging from 5 to 30 W. Figure 7 compares the experimental and simulation of the heat source and HTF outlet temperatures for a heat sink layers ranging between one and four-layers and of a HTF volume flow rate of 30 mL/min and HTF inlet temperature of 24 • C. The average difference between the numerical and experimental of the HTF outlet temperature is 1.25 • C and a maximum difference of 3.30 • C occurred at 25 W of the three-layers MLM heat sink.However, a more average difference of the heat source temperature between the numerical and experimental is found of 2.58 • C.This shows a good estimation of the numerical model and will be used to explore the thermal characteristics of the system.

Effects of Number of Layers
The effect of number of layers on the performance of the MLM heat sink is investigated in this section for a variable number of layers ranging between one and four.Figure 8 shows the thermal resistance of the MLM heat sink and heat source surface temperature for different numbers of layers

Effects of Number of Layers
The effect of number of layers on the performance of the MLM heat sink is investigated in this section for a variable number of layers ranging between one and four.Figure 8 shows the thermal resistance of the MLM heat sink and heat source surface temperature for different numbers of layers of a HTF flow rate of 30 mL/min and 15 W power rate.The thermal resistance decreases by 17% as the number of layers increases from one-layer to two-layers.However, a slight decrease of the thermal resistance of 2% is noticed when the number of layers increased to three-layers and four-layers.In addition, the figure shows the heat source surface temperature for the various number of layers arrangement.The temperature decreases as the number of layers increases from one-layer to two-layers by 5.21 • C and less temperature reduction is noticed when increasing the number of layers from three to four.Further investigation on the temperature uniformity on the heat source top surface temperature was conducted.Figure 9 shows the measured temperature distributions of the surface of the heat source in nine sections for the different number of layers heat sinks.The results show that the average and the standard deviation of the measured temperatures decreases as the number of layers decreases.It can be clearly noticed that the center of the heat source (section 1) suffers the maximum temperature of all the other sections and a wide difference between this temperature and the other locations within the heat source of a maximum difference of 12.01 °C in the two-layers MLM heat sink.This temperature difference causes an increase of temperature non-uniformity across the heat source surface and hence the standard deviation.Despite the maximum temperature value (section 1), the single layer MLM heat sink is considered to be the most temperature variable around the heat source (sections 2-9) with a standard deviation of 4.30 °C.Interestingly, two-layers MLM heat sinks experiences less temperature non-uniformity compared to the three-layers MLM heat sink with a standard deviation of 3.86 and 4.05 °C, respectively.investigation on the temperature uniformity on the heat source top surface temperature was conducted.Figure 9 shows the measured temperature distributions of the surface of the heat source in nine sections for the different number of layers heat sinks.The results show that the average and the standard deviation of the measured temperatures decreases as the number of layers decreases.It can be clearly noticed that the center of the heat source (section 1) suffers the maximum temperature of all the other sections and a wide difference between this temperature and the other locations within the heat source of a maximum difference of 12.01 • C in the two-layers MLM heat sink.This temperature difference causes an increase of temperature non-uniformity across the heat source surface and hence the standard deviation.Despite the maximum temperature value (section 1), the single layer MLM heat sink is considered to be the most temperature variable around the heat source (sections 2-9) with a standard deviation of 4.30 • C. Interestingly, two-layers MLM heat sinks experiences less temperature non-uniformity compared to the three-layers MLM heat sink with a standard deviation of 3.86 and 4.05 • C, respectively.
A close look into the temperature distribution across the heat source surface was obtained using the contour representation from the numerical solution and is shown in Figure 10.It can be noticed that the hot spot shifted to the center of the heat source as the number of layers increases.In addition, the temperature distribution of the heat source edge near to the downstream becomes more uniform as the number of layers increased.There is little visible difference between the three-layers and four-layers heat sinks in terms of the hot spot distribution and the temperature uniformity.

Inlet Outlet
sink.This temperature difference causes an increase of temperature non-uniformity across the heat source surface and hence the standard deviation.Despite the maximum temperature value (section 1), the single layer MLM heat sink is considered to be the most temperature variable around the heat source (sections 2-9) with a standard deviation of 4.30 °C.Interestingly, two-layers MLM heat sinks experiences less temperature non-uniformity compared to the three-layers MLM heat sink with a standard deviation of 3.86 and 4.05 °C, respectively.A close look into the temperature distribution across the heat source surface was obtained using the contour representation from the numerical solution and is shown in Figure 10.It can be noticed that the hot spot shifted to the center of the heat source as the number of layers increases.In addition, the temperature distribution of the heat source edge near to the downstream becomes more uniform as the number of layers increased.There is little visible difference between the three-layers and fourlayers heat sinks in terms of the hot spot distribution and the temperature uniformity.An isotherm of the full heat sink module for the different layer heat sinks is shown in Figure 11.It can be noticed that less temperature is experienced by the microchannel plates especially at the edges.This implies an opportunity of using the MLM heat sink in the CPV applications by decreasing the CPV receiver assembly to a smaller size so that fewer heat sinks materials are used and less fluid flow.An isotherm of the full heat sink module for the different layer heat sinks is shown in Figure 11.It can be noticed that less temperature is experienced by the microchannel plates especially at the edges.This implies an opportunity of using the MLM heat sink in the CPV applications by decreasing the CPV receiver assembly to a smaller size so that fewer heat sinks materials are used and less fluid flow.The thermal efficiency of the MLM heat sink and the HTF outlet temperature are shown in Figure 12.As seen from the figure, the thermal efficiency has been improved significantly by 20% when the number of layers increases from one-to three-layers.This is due to the increase of the HTF outlet temperature from 29.7 to 31.1 °C for the one-layer and three-layers MLM heat sinks, respectively.However, no effect on the thermal efficiency and the HTF outlet temperature was noticed when the number of layers increased from three-to four-layers.The hydrodynamic characteristics are as important as the thermal performance in evaluating the effect of the number of layers in the MLM heat sink.Pressure drop across the heat sink is one of the characteristics that effects the pumping power.Figure 13 shows the effect of layers in the pressure drop.It can be noticed that the pressure drop decreased as the number of layers increased.An approximate 39.22 Pa pressure drop reduction occurred with the increase of each layer, where the pressure for the MLM heat sink of one-layer and three-layers are 162.79 and 81.30Pa, respectively.The thermal efficiency of the MLM heat sink and the HTF outlet temperature are shown in Figure 12.As seen from the figure, the thermal efficiency has been improved significantly by 20% when the number of layers increases from one-to three-layers.This is due to the increase of the HTF outlet temperature from 29.7 to 31.1 • C for the one-layer and three-layers MLM heat sinks, respectively.However, no effect on the thermal efficiency and the HTF outlet temperature was noticed when the number of layers increased from three-to four-layers.The thermal efficiency of the MLM heat sink and the HTF outlet temperature are shown in Figure 12.As seen from the figure, the thermal efficiency has been improved significantly by 20% when the number of layers increases from one-to three-layers.This is due to the increase of the HTF outlet temperature from 29.7 to 31.1 °C for the one-layer and three-layers MLM heat sinks, respectively.However, no effect on the thermal efficiency and the HTF outlet temperature was noticed when the number of layers increased from three-to four-layers.The hydrodynamic characteristics are as important as the thermal performance in evaluating the effect of the number of layers in the MLM heat sink.Pressure drop across the heat sink is one of the characteristics that effects the pumping power.Figure 13 shows the effect of layers in the pressure drop.It can be noticed that the pressure drop decreased as the number of layers increased.An approximate 39.22 Pa pressure drop reduction occurred with the increase of each layer, where the pressure for the MLM heat sink of one-layer and three-layers are 162.79 and 81.30Pa, respectively.The hydrodynamic characteristics are as important as the thermal performance in evaluating the effect of the number of layers in the MLM heat sink.Pressure drop across the heat sink is one of the characteristics that effects the pumping power.Figure 13 shows the effect of layers in the pressure drop.It can be noticed that the pressure drop decreased as the number of layers increased.An approximate 39.22 Pa pressure drop reduction occurred with the increase of each layer, where the pressure for the MLM heat sink of one-layer and three-layers are 162.79 and 81.30Pa, respectively.

Effects of Power Rates
Total thermal resistances for different power input rate for the three-layers MLM heat sink with 30 mL/min flow rate are shown in Figure 14.As can see from the figure, the total thermal resistance of the heat sink increases slightly with increasing the power rate from 5 to 30 W for the same HTF flow rate of a difference of 0.08 K/W.This shows the ability of the heat sink to undertake a wide range of power rates with a slight change in the thermal resistance.
As expected, the heat source temperature increases significantly as heat source power rate increases from 5 to 30 W for the same HTF flow rate.The temperature of the heat source is doubled when the power rate increases from 5 to 30 W, with a difference of 45 °C.However, the heat source temperature on the worst case at 30 W is 77.48 °C , which is within the recommended temperature of the solar cell manufacturer (80 °C). Figure 15 shows the temperature distributions on the surface of the heat source, average temperature and standard deviation for the three-layers MLM heat sink of 30 mL/min HTF flow rate for various power rates.The results show an increase of the temperature non-uniformity within the heat source as the power rate increases and this is the same for the standard deviation.For the 5 W power rate, the temperature varies with a standard deviation of 1.25 °C and with an average temperature of 29.18 °C.A significant increase of the standard deviation in the 30 W power rate of 9.03 °C in the 30 W power rate and of an average temperature of 54.09 °C.

Effects of Power Rates
Total thermal resistances for different power input rate for the three-layers MLM heat sink with 30 mL/min flow rate are shown in Figure 14.As can see from the figure, the total thermal resistance of the heat sink increases slightly with increasing the power rate from 5 to 30 W for the same HTF flow rate of a difference of 0.08 K/W.This shows the ability of the heat sink to undertake a wide range of power rates with a slight change in the thermal resistance.

Effects of Power Rates
Total thermal resistances for different power input rate for the three-layers MLM heat sink with 30 mL/min flow rate are shown in Figure 14.As can see from the figure, the total thermal resistance of the heat sink increases slightly with increasing the power rate from 5 to 30 W for the same HTF flow rate of a difference of 0.08 K/W.This shows the ability of the heat sink to undertake a wide range of power rates with a slight change in the thermal resistance.
As expected, the heat source temperature increases significantly as heat source power rate increases from 5 to 30 W for the same HTF flow rate.The temperature of the heat source is doubled when the power rate increases from 5 to 30 W, with a difference of 45 °C.However, the heat source temperature on the worst case at 30 W is 77.48 °C , which is within the recommended temperature of the solar cell manufacturer (80 °C). Figure 15 shows the temperature distributions on the surface of the heat source, average temperature and standard deviation for the three-layers MLM heat sink of 30 mL/min HTF flow rate for various power rates.The results show an increase of the temperature non-uniformity within the heat source as the power rate increases and this is the same for the standard deviation.For the 5 W power rate, the temperature varies with a standard deviation of 1.25 °C and with an average temperature of 29.18 °C.A significant increase of the standard deviation in the 30 W power rate of 9.03 °C in the 30 W power rate and of an average temperature of 54.09 °C.As expected, the heat source temperature increases significantly as heat source power rate increases from 5 to 30 W for the same HTF flow rate.The temperature of the heat source is doubled when the power rate increases from 5 to 30 W, with a difference of 45 • C.However, the heat source temperature on the worst case at 30 W is 77.48 • C, which is within the recommended temperature of the solar cell manufacturer (80 • C).
Figure 15 shows the temperature distributions on the surface of the heat source, average temperature and standard deviation for the three-layers MLM heat sink of 30 mL/min HTF flow rate for various power rates.The results show an increase of the temperature non-uniformity within the heat source as the power rate increases and this is the same for the standard deviation.For the 5 W power rate, the temperature varies with a standard deviation of 1.25  The thermal efficiency of the heat sink and the HTF outlet temperature for various power rates for the three-layers MLM heat sink of 30 mL/min HTF flow rate are shown in Figure 16.As can be seen from the figure, the thermal efficiency drops significantly by 6% when the power rate increases from 5 to 10 W. However, the thermal efficiency decreases slightly when the power rate increases from 10 to 30 W and reaches 90% for the 30 W power rate.Additionally, the HTF outlet temperature increases as the power rate increases from 5 to 30 W of 26.45 to 37.88 °C.This shows a higher heat dissipation of a higher power rate.

Effects of HTF Flow Rates
In this section the effect HTF flow rates on the thermal performance of the heat source is presented for various number of layers MLM heat sinks.Figure 17 shows effect of the HTF flow rates on the thermal resistance of the MLM heat sink for the 5 W power rate.It can be seen that for a single layer heat sink, the thermal resistance decreases from 2.04 to 1.78 K/W as HTF flow rate increases from 30 to 60 mL/min, respectively.A significant decrease of the thermal resistance is noticed as the HTF flow rate increases from 30 to 40 mL/min of the three-and four-layers heat sink with a difference The thermal efficiency of the heat sink and the HTF outlet temperature for various power rates for the three-layers MLM heat sink of 30 mL/min HTF flow rate are shown in Figure 16.As can be seen from the figure, the thermal efficiency drops significantly by 6% when the power rate increases from 5 to 10 W. However, the thermal efficiency decreases slightly when the power rate increases from 10 to 30 W and reaches 90% for the 30 W power rate.Additionally, the HTF outlet temperature increases as the power rate increases from 5 to 30 W of 26.45 to 37.88 • C.This shows a higher heat dissipation of a higher power rate.The thermal efficiency of the heat sink and the HTF outlet temperature for various power rates for the three-layers MLM heat sink of 30 mL/min HTF flow rate are shown in Figure 16.As can be seen from the figure, the thermal efficiency drops significantly by 6% when the power rate increases from 5 to 10 W. However, the thermal efficiency decreases slightly when the power rate increases from 10 to 30 W and reaches 90% for the 30 W power rate.Additionally, the HTF outlet temperature increases as the power rate increases from 5 to 30 W of 26.45 to 37.88 °C.This shows a higher heat dissipation of a higher power rate.

Effects of HTF Flow Rates
In this section the effect HTF flow rates on the thermal performance of the heat source is presented for various number of layers MLM heat sinks.Figure 17 shows effect of the HTF flow rates on the thermal resistance of the MLM heat sink for the 5 W power rate.It can be seen that for a single layer heat sink, the thermal resistance decreases from 2.04 to 1.78 K/W as HTF flow rate increases from 30 to 60 mL/min, respectively.A significant decrease of the thermal resistance is noticed as the

Effects of HTF Flow Rates
In this section the effect HTF flow rates on the thermal performance of the heat source is presented for various number of layers MLM heat sinks.Figure 17 shows effect of the HTF flow rates on the thermal resistance of the MLM heat sink for the 5 W power rate.It can be seen that for a single layer heat sink, the thermal resistance decreases from 2.04 to 1.78 K/W as HTF flow rate increases from 30 to 60 mL/min, respectively.A significant decrease of the thermal resistance is noticed as the HTF flow rate increases from 30 to 40 mL/min of the three-and four-layers heat sink with a difference of 0.15 and 0.14 K/W, respectively.An interesting note is that increasing the number of layers from three to four has a slight effect on the thermal resistance for the HTF flow rate of 30, 40, and 60 mL/min and this shows that the optimum number of layers in the 5 W power rating is three-layers.
Energies 2019, 12, x FOR PEER REVIEW 20 of 25 of 0.15 and 0.14 K/W, respectively.An interesting note is that increasing the number of layers from three to four has a slight effect on the thermal resistance for the HTF flow rate of 30, 40, and 60 mL/min and this shows that the optimum number of layers in the 5 W power rating is three-layers.Figure 18 shows the effect of the number of layers and HTF flow rates in reducing the heat sink thermal resistance for a higher power rate of 30 W. At this power rating, both the increase of the number of layers and HTF flow rate have a major impact in reducing the thermal resistance of the heat sink HTF.The thermal resistance reduces by 23.26% and 10.25% when the number of layers increases from one to four layers for 30 to 60 mL/min, respectively.This shows that the power rate is critical in determining the optimum number of layers for the lower thermal resistance.A summary of the measurements of the heat source temperature for various HTF flow rates are shown in Table 8.It is noticed that the heat source temperature decreases significantly as the HTF flow rate increases for all the MLM heat sinks arrangements for power rates of more than 15 W. In addition, the heat source experiences a temperature of more than 80 °C when using a single layer heat sink for power rate of 30 W for all the tested HTF flow rates ranging from 30 to 60 mL/min, Figure 18 shows the effect of the number of layers and HTF flow rates in reducing the heat sink thermal resistance for a higher power rate of 30 W. At this power rating, both the increase of the number of layers and HTF flow rate have a major impact in reducing the thermal resistance of the heat sink HTF.The thermal resistance reduces by 23.26% and 10.25% when the number of layers increases from one to four layers for 30 to 60 mL/min, respectively.This shows that the power rate is critical in determining the optimum number of layers for the lower thermal resistance.
Energies 2019, 12, x FOR PEER REVIEW 20 of 25 of 0.15 and 0.14 K/W, respectively.An interesting note is that increasing the number of layers from three to four has a slight effect on the thermal resistance for the HTF flow rate of 30, 40, and 60 mL/min and this shows that the optimum number of layers in the 5 W power rating is three-layers.Figure 18 shows the effect of the number of layers and HTF flow rates in reducing the heat sink thermal resistance for a higher power rate of 30 W. At this power rating, both the increase of the number of layers and HTF flow rate have a major impact in reducing the thermal resistance of the heat sink HTF.The thermal resistance reduces by 23.26% and 10.25% when the number of layers increases from one to four layers for 30 to 60 mL/min, respectively.This shows that the power rate is critical in determining the optimum number of layers for the lower thermal resistance.A summary of the measurements of the heat source temperature for various HTF flow rates are shown in Table 8.It is noticed that the heat source temperature decreases significantly as the HTF flow rate increases for all the MLM heat sinks arrangements for power rates of more than 15 W. In addition, the heat source experiences a temperature of more than 80 °C when using a single layer heat sink for power rate of 30 W for all the tested HTF flow rates ranging from 30 to 60 mL/min, which is more than the recommended maximum temperature specified by the CPV assembly  A summary of the measurements of the heat source temperature for various HTF flow rates are shown in Table 8.It is noticed that the heat source temperature decreases significantly as the HTF flow rate increases for all the MLM heat sinks arrangements for power rates of more than 15 W. In addition, the heat source experiences a temperature of more than 80 • C when using a single layer heat sink for power rate of 30 W for all the tested HTF flow rates ranging from 30 to 60 mL/min, which is more than the recommended maximum temperature specified by the CPV assembly manufacturer for performing a high electrical efficiency.A summary of the HTF outlet temperature i.e., the dissipated heated for the different flow rates, is presented in Table 9. Increasing the HTF flow rate decreases the HTF outlet temperature where the HTF outlet temperature in the 30 mL/min reaches the maximum among all the flow rates.This is more noticeable in the high power rates where the HTF outlet temperature difference is around 3.5 • C when the HTF flow rate increases from 30 to 40 mL/min for all numbers of layers.The maximum HTF outlet temperatures reaches its maximum in the four-layers heat sink of 38.36 • C. The pressure drop of the MLM heat sink for different number of layer is shown in Figure 19.It is observed that the pressure drop for the flow rate of 30 mL/min is the maximum at the single layer MLM heat sink of 162.79 Pa and its minimum in the four-layers heat sink of 32.75 Pa which is a reduction of three times compared to the single layer heat sink.The heat sink of a single layer experiences a significant increase in the pressure drop when the HTF flow rates are increased from 30 to 60 mL/min, with an increase of 44.13 Pa.However, a slight increase in the pressure drop is noticed when the flow rate is increased for the MLM heat sink of the two-, three-and four-layers MLM heat sinks, with increases of only of 14.81, 9.90, and 15.49Pa, respectively.

Effects of Number of Layers on the Pressure Drop
The pressure drop of the MLM heat sink for different number of layer is shown in Figure 19.It is observed that the pressure drop for the flow rate of 30 mL/min is the maximum at the single layer MLM heat sink of 162.79 Pa and its minimum in the four-layers heat sink of 32.75 Pa which is a reduction of three times compared to the single layer heat sink.The heat sink of a single layer experiences a significant increase in the pressure drop when the HTF flow rates are increased from 30 to 60 mL/min, with an increase of 44.13 Pa.However, a slight increase in the pressure drop is noticed when the flow rate is increased for the MLM heat sink of the two-, three-and four-layers MLM heat sinks, with increases of only of 14.81, 9.90, and 15.49Pa, respectively.

Conclusions
An experimental and numerical study of a multi-layered micro channel (MLM) heat sink for CPV thermal regulation application was carried out.The effect of the number of layers and HTF flow rate on the performance of the MLM heat sink at various heat fluxes was investigated.The number of layers and heat fluxes were varied from one to four and 5 to 30 W/cm 2 , respectively.The following is concluded from this work:


Increasing the number of layers has a major impact in reducing the heat sink thermal resistance and the heat source maximum temperature.The thermal resistance decreases by 17% as the number of layers increases from one-layer to two-layers.A slight decrease of the thermal resistance is observed when the number of layers increased to three and four.


The temperature non-uniformity decreases by 0.5 °C when the number of layers increases from one-layer to two-layers.However, the two-layers MLM heat sink experiences less temperature

Conclusions
An experimental and numerical study of a multi-layered micro channel (MLM) heat sink for CPV thermal regulation application was carried out.The effect of the number of layers and HTF flow rate on the performance of the MLM heat sink at various heat fluxes was investigated.The number of layers and heat fluxes were varied from one to four and 5 to 30 W/cm 2 , respectively.The following is concluded from this work:

•
Increasing the number of layers has a major impact in reducing the heat sink thermal resistance and the heat source maximum temperature.The thermal resistance decreases by 17% as the number of layers increases from one-layer to two-layers.A slight decrease of the thermal resistance is observed when the number of layers increased to three and four.

•
The temperature non-uniformity decreases by 0.5 • C when the number of layers increases from one-layer to two-layers.However, the two-layers MLM heat sink experiences less temperature non-uniformity compared to the three-layers.

•
The thermal efficiency is improved significantly by increasing the number of layers from one-layer to three-layers by 20%, due to the increase in the HTF outlet temperature.

•
Increasing the input power from 5 to 30 W slightly increases the heat sink total thermal resistance and the heat source temperature increases.This shows the ability of the heat sink to accommodate a wide range of power rates with a slight change in the thermal resistance.
This study is one work of CPV-MLM investigations.The future work will consider applying the MLM heat sink in single solar cell CPV module.The investigation will include both the electrical and thermal performance using the indoor and outdoor characterizations.

Figure 2 .
Figure 2. (a) Schematic of the test module parts, (b) a photograph of the test module showing the thermocouple locations at; (1) water inlet temperature, (2) water temperature at the middle channel of the first layer, (3) heat source middle surface temperature, (4) water outlet temperature.

Figure 2 .
Figure 2. (a) Schematic of the test module parts, (b) a photograph of the test module showing the thermocouple locations at; (1) water inlet temperature, (2) water temperature at the middle channel of the first layer, (3) heat source middle surface temperature, (4) water outlet temperature.
thermocouple locations at; (1) water inlet temperature, (2) water temperature at the middle channel of the first layer, (3) heat source middle surface temperature, (4) water outlet temperature.

Figure 4 .
Figure 4. (a) Schematic of the three stacked layers microchannel: width of channels (Wc), height of fins (Hc), thickness of fins (Wf), length of plate (L), width of plate (W), base thickness (tb); (b) microscopic images of the side of the plate; (c) surface finish under the microscope.

Figure 4 .Table 2 .
Figure 4. (a) Schematic of the three stacked layers microchannel: width of channels (W c ), height of fins (H c ), thickness of fins (W f ), length of plate (L), width of plate (W), base thickness (t b ); (b) microscopic images of the side of the plate; (c) surface finish under the microscope.Table 2. The dimensions of the microchannel.

25 A
list of controlling parameters are shown in Table

Figure 5 .
Figure 5.The schematic of the multi-layer microchannel heat sink and CPV receiver.

Figure 5 .
Figure 5.The schematic of the multi-layer microchannel heat sink and CPV receiver.

Figure 5 .
Figure 5.The schematic of the multi-layer microchannel heat sink and CPV receiver.

Figure 6 .
Figure 6.Schematic of the computational domain and boundary conditions: 1 copper layer surface; 2 ERH; 3 top and bottom multi-layered outer surfaces.

Energies 2019 ,Figure 7 .
Figure 7. Validation of simulation results compared with experimental data for 30 mL/min water flow rates for layers ranged from one to four for (a) the temperature of the heat source in the middle top surface; (b) HTF outlet temperature versus the power.

Figure 7 .
Figure 7. Validation of simulation results compared with experimental data for 30 mL/min water flow rates for layers ranged from one to four for (a) the temperature of the heat source in the middle top surface; (b) HTF outlet temperature versus the power.

Energies 2019 ,
12,  x FOR PEER REVIEW 15 of 25 thermal resistance of 2% is noticed when the number of layers increased to three-layers and fourlayers.In addition, the figure shows the heat source surface temperature for the various number of layers arrangement.The temperature decreases as the number of layers increases from one-layer to two-layers by 5.21 °C and less temperature reduction is noticed when increasing the number of layers from three to four.

Figure 8 .
Figure 8. Effects of number of layers on the thermal resistance and heat source temperature of 30 mL/min HTF flow rate and 15 W power rate.

Figure 8 .
Figure 8. Effects of number of layers on the thermal resistance and heat source temperature of 30 mL/min HTF flow rate and 15 W power rate.

Figure 9 . 25 Figure 9 .
Figure 9.The heat source temperature uniformity of 30 mL/min HTF flow rate and 15 W power rate for (a) one-layer; (b) two-layers; (c) three-layers MLM heat sink; (d) Schematic of the thermocouple distribution on heat source surface.

Figure 10 .
Figure 10.The heat source surface temperature contours MLM heat sink of 30 mL/min HTF flow rate and 15 W power rate for (a) one-layer; (b) two-layers; (c) three-layers; (d) 4-layers.

Figure 10 .
Figure 10.The heat source surface temperature contours MLM heat sink of 30 mL/min HTF flow rate and 15 W power rate for (a) one-layer; (b) two-layers; (c) three-layers; (d) 4-layers.

Figure 12 .
Figure 12.Effects of number of layers on the HTF outlet temperature and thermal performance of 30 mL/min HTF flow rate and 15 W power rate.

Figure 12 .
Figure 12.Effects of number of layers on the HTF outlet temperature and thermal performance of 30 mL/min HTF flow rate and 15 W power rate.

Figure 12 .
Figure 12.Effects of number of layers on the HTF outlet temperature and thermal performance of 30 mL/min HTF flow rate and 15 W power rate.

Figure 13 .
Figure 13.Effects of number of layers on the pressure drop across the heat sink of 30 mL/min HTF flow rate and 15 W power rate.

Figure 14 .
Figure 14.Effects of number of layers on the thermal resistance and heat source temperature of 30 mL/min HTF flow rate and 15 W power rate.

Figure 13 .
Figure 13.Effects of number of layers on the pressure drop across the heat sink of 30 mL/min HTF flow rate and 15 W power rate.

Energies 2019 , 25 Figure 13 .
Figure 13.Effects of number of layers on the pressure drop across the heat sink of 30 mL/min HTF flow rate and 15 W power rate.

Figure 14 .
Figure 14.Effects of number of layers on the thermal resistance and heat source temperature of 30 mL/min HTF flow rate and 15 W power rate.

Figure 14 .
Figure 14.Effects of number of layers on the thermal resistance and heat source temperature of 30 mL/min HTF flow rate and 15 W power rate.

Figure 15 .
Figure 15.The heat source temperature uniformity for 3-layers MLM heat sink of 30 mL/min HTF flow rate for a power rate of (a) 5 W, (b) 15 W and (c) 30 W; (d) Schematic of the thermocouple distribution on heat source surface.

Figure 16 .
Figure 16.Effects of power rating on the HTF outlet temperature and thermal performance of threelayers MLM heat sink of 30 mL/min HTF flow rate.

Figure 15 .
Figure 15.The heat source temperature uniformity for 3-layers MLM heat sink of 30 mL/min HTF flow rate for a power rate of (a) 5 W, (b) 15 W and (c) 30 W; (d) Schematic of the thermocouple distribution on heat source surface.

Figure 15 .
Figure 15.The heat source temperature uniformity for 3-layers MLM heat sink of 30 mL/min HTF flow rate for a power rate of (a) 5 W, (b) 15 W and (c) 30 W; (d) Schematic of the thermocouple distribution on heat source surface.

Figure 16 .
Figure 16.Effects of power rating on the HTF outlet temperature and thermal performance of threelayers MLM heat sink of 30 mL/min HTF flow rate.

Figure 16 .
Figure 16.Effects of power rating on the HTF outlet temperature and thermal performance of three-layers MLM heat sink of 30 mL/min HTF flow rate.

Figure 17 .
Figure 17.Heat sink thermal resistance versus number of layers for different HTF flow rates for 5 W.

Figure 18 .
Figure 18.Heat sink thermal resistance versus number of layers for different HTF flow rates for 30 W.

Figure 17 .
Figure 17.Heat sink thermal resistance versus number of layers for different HTF flow rates for 5 W.

Figure 17 .
Figure 17.Heat sink thermal resistance versus number of layers for different HTF flow rates for 5 W.

Figure 18 .
Figure 18.Heat sink thermal resistance versus number of layers for different HTF flow rates for 30 W.

Figure 18 .
Figure 18.Heat sink thermal resistance versus number of layers for different HTF flow rates for 30 W.

Figure 19 .
Figure 19.The pressure drop versus number of layers for different HTF flow rates.

Figure 19 .
Figure 19.The pressure drop versus number of layers for different HTF flow rates.

Table 1 .
Comparative study of different multi-layered microchannel (MLM) heat sink literatures.

Table 3 .
The volume and mass flow rates of water.

Table 4 .
The volume and mass flow rates of water.

Table 6 .
The model thermal boundary conditions.

Table 7 .
Calculated temperature variation with number of elements for three-layers, 30 mL/min and 15 W heating power.
• C and with an average temperature of 29.18 • C. A significant increase of the standard deviation in the 30 W power rate of 9.03 • C in the 30 W power rate and of an average temperature of 54.09 • C.

Table 8 .
The heat source maximum temperature for various HTF flow rates and power rates for the different layers MLM heat sinks.

Table 9 .
The HTF outlet temperature for various HTF flow rates and power rates for the different layers MLM heat sinks.

Table 9 .
Cont.Effects of Number of Layers on the Pressure Drop