Experimental and Theoretical Analysis of Energy Efﬁciency in a Flat Plate Solar Collector Using Monolayer Graphene Nanoﬂuids

: Flat-plate solar collectors are one of the cleanest and most efﬁcient heating systems available. Studies on the presence of covalently functionalized graphene (Gr) suspended in distilled water as operating ﬂuids inside an indoor ﬂat-plate solar collector (FPSC) were experimentally and theoretically performed. These examinations were conducted under different testing conditions namely 0.025 wt.%, 0.05 wt.%, 0.075 wt.%, and 0.1 wt.%, 0.5, 1, and 1.5 kg/min, 30, 40, and 50 ◦ C, and 500, 750, and 1000 W/m 2 . Various techniques were used to characterize the functionalized nanoﬂuids’ stability and morphological properties namely UV/Vis spectrophotometry, EDX analysis with a Scanning Electron Microscope (SEM), zeta potential, and nanoparticle size. The results showed that the collected heat improved as the percentage of GrNPs and the ﬂuid mass ﬂow rates increased, although it decreased as the reduced temperature coefﬁcient increased, whereas the maximum increase in collector efﬁciency at higher concentration was 13% and 12.5% compared with distilled water at 0.025 kg/s. Finally, a new correlation was developed for the base ﬂuid and nanoﬂuids’ thermal efﬁciency as a function of dropped temperature parameter and weight concentration with 2.758% and 4.232% maximum deviations.


Introduction
Flat-plate solar collectors (FPSCs) are the most basic and studied systems for solarpowered hot water technologies in domestic water. The usage of FPSC is widespread for low-and medium-heat applications. Hence, many advancements and continuous improvements have been undertaken to reduce the size of FPSCs and increase their performance. A black surface absorber, riser, and header pipes with coolant, a glazing cover, coolant, and thermal insulation are the essential components of flat plate collectors. The investigations of Hottel and Woertz in 1942 [1] and Hottel and Whiller in 1958 [2], were among the first studies on FPSC. In 1955, Tabor [3] used selected black surfaces to enhance the heat collected performance. His optical absorption studies showed an optical concentration capability to generate high-pressure steam. As a consequence, numerous experiments conditions, study [23] was conducted to address pentaethylene glycol (PEG) behavior decorated GNPs-H 2 O nanofluids inside an indoor flat plate solar collector. Improvements in collector performance were recorded at around 10.7%, 11.1%, and 13.3%, respectively, at 0.00833, 0.01667, and 0.025 kg/s. The most thermal-efficient solar collector improved up to 85% with (CF-MWCNTs + CF-GNPs + h-BN) hybrid nanofluid as the absorption medium at 4 L/min flow rate [24]. A comparison of 0.1 wt.% water-based nanofluids can be sequenced f-GNPs > ZnO > SiO 2 because of a percentage improvement of thermal efficiency of the flat-plate solar collector obtained at a mass flow rate of 1.6 kg/min with values of 17.45 > 13.05 > 12.36%, respectively in comparison to water [25]. The maximum enhancement in LFPSC efficiency was 24.09% for 0.1 wt.% GGNPs and flow rate of 1.5 L/min more than distilled water. Meanwhile, exergy efficiency was maximum for 0.1% GGNP concentration and flow rate of 0.8 L/min [26].
It can be seen from the above literature review that the enhancement of FPSC research using graphene nanofluids is still in progress, and needs to be thoroughly investigated. Monolayer graphene is a carbon-based material with thermal conductivity in the range of 2000-4000 W/m K. In comparison to other nanoparticles, this value is very high. As a result, the current research focuses on evaluating collector performance both experimentally and theoretically. Graphene nanomaterials were synthesized and characterized with different techniques, namely spectrophotometry, Electron Microscope, zeta potential, and nanoparticle sizing. Four samples of based-DW were tested with different testing parameters, namely mass flow rates, input temperatures, and heating fluxes. Besides this, new empirical correlations were developed for predicting the thermal efficiency of DW and Gr-nanofluids as a function of dropped temperature parameter and weight concentration with three coefficients at different flow rates.

Experimental Approach
In this section, the nanomaterials and chemicals for covalent functionalization are introduced. Moreover, the procedure used to prepare the graphene-based water nanofluid is discussed. Then, the experimental indoor setup is presented. In addition, characterization and thermophysical properties are specified.

Synthesis of Gr-Nanofluids
Monolayer Graphene (Gr) was purchased abroad from VCN-Materials Co., Ltd., Bushehr, Iran. While the chemical reaction approach was carried out using different reagents and materials were locally obtained from Sigma-Aldrich: Pentaethylene glycol (PEG), strong medium (AlCl 3 :HCl), tetrahydrofuran (THF), and N, N-dimethylformamide (DMF). The adding of (AlCl 3 :HCl) to the GrNPs introduced specific groups like (-OH) and (-COOH) on the surface affinity; the details can be found in [27] (see Figure 1). To prepare the experimental samples of GrNPs-DW, Pioneer Precision Electronic Balance (OHAUS PX163, Parsippany, NJ, USA) with Std. Dev. of 0.1 mg and a probe (Vibra-Cell, Sonics, VC 750, Newtown, CT, USA) were used.

Experimental Setup
The setup configuration of an indoor-FPSC was used here to study the behavior of using graphene nanocoolants rather than conventional fluids (see Figure 2). The testing system consists of an absorber plate, controlling and measuring tools, an electrical centrifugal pump, a water-cooling chiller bath (WiseCircu WCR-P22, Seoul, Korea), a flow piping loop, and a digital data recorder (EC18, paperless Ecolog 18-channel recorder, Japan). The absorbing copper sheet was completely soldered over the contact length of the copper riser pipes. An isowool ceramic fiber blanket was installed with a thickness of 50 mm (back), 30 mm (sides), and thermal conductivity of 0.07 W/m-K at 400 • C to minimize the heat loss surroundings. A configuration of the flexible adhesive heater/adjustable voltage transformer (KH-1012/(2.5)-P, from Omega, Norwalk, CT, USA) was used to control the input power (heat flux intensity). Digital-sense T-thermocouple self-adhesive surface probes (model; SA1XL-T-72, Omega, USA) were installed at twelve locations to estimate the riser's surface temperatures of pipes and flat plate. Furthermore, the bulk temperature was measure by inserting two RTDs (PT100-type, Omega, USA) at the fluid in/out.

Characterization and Thermophysical Properties
The measurements of ultraviolet-visible spectroscopy (UV/Vis) and nanoparticle sizing were carried out to verify the stability of the graphene dispersions. UV/Vis spectroscopy analysis of light absorbance for a suspension can be used to provide a quantifiable stability characterization. Perkin Elmer (Lambda 750, Waltham, MA, USA) was used for measuring the UV/Vis within the range of (190-3300) nm. The Anton Paar instrument (Litesizer 500, Graz, Austria) was used with the Electrophoretic Light Scattering (ELS) concept for calculating zeta potential and nanoparticle size measurements. For zeta-potential analysis [28], the degree of aversion between near-nanoparticles shows the same nanofluid dispersal effect. The Tescan VEGA3 (Brno, Czech Republic) tool was used to examine the morphology and elemental properties of GrNPs.
In this work, to prepare the GrNPs-nanofluids with the desired concentration, nanoparticles were equivalently dispersed into DW [29]. The viscosity, thermal conductivity, specific

Data Collection and Error Analysis
The useful heat gain of the FPSC-fluid (Q u ) was calculated using the following Equations (1) and (2) [30][31][32]: In this case, .
(m) is the mass flow rate, (C p ) indicates the fluid-specific heat. (T out and T in ) are the output and input fluid temperatures. (F R ) refers to the heat removal coefficient, (A c ) indicates the gross area of the collector, (G T ) states the heat flow intensity, (τα) specifies the effective absorptance-transmittance product, (U L ) shows the overall heat loss coefficient. (T amb ) is the ambient temperature. (η) denotes the FPSC-efficiency as the ratio of the usable energy gain to the solar energy incident over a given time [30][31][32]: Experimental uncertainty assessment deals with assessing the uncertainty in a measurement. The uncertainty for FPSC-efficiency is an equation of five factors, namely mass flow rate, specific heat capacity, the intensity of solar radiation, collector area, and the difference between output/input fluid temperatures. Equation (5) shows the combination formula of the uncertainty for FPSC-efficiency [20]: Table 1 demonstrates the ranges and accuracies of the instruments, equipment, and working fluid properties; after the measurement procedure, the overall uncertainty of the FPSC-efficiency is around 2.92%.

The Solar Collector Geometry
The essential elements of the FPSC are the glazing frame with one/more glass sheets, heat-absorbing backplate, heat-absorbing riser tubes with lower manifold collections tubes, sides, and bottom insulation. For the best design, various materials for all these components need to be measured for their properties and configurations. It is also necessary to study the processes and techniques to increase its efficiency further. Some assumptions for modelling solar collectors on the flat plate are as follows [33,34]: (1) The FPSC works under stable state conditions. (2) The covered area of header pipes can be ignored.

The Collector Heat Loss
As the collector absorbs heat and heat is lost to the atmosphere by convection and radiation, the FPSC-surface temperature becomes higher than that of the ambient. The heat loss rate (Q loss ) depends on the overall heat transfer coefficient (U L ) and the flat plate mean temperature (T pm ) [33,34]: where, (t, b, and e) denote the subscripts for top, back, and edge, respectively. An empirical equation for estimating the maximum collector-to-environment loss coefficient (U t ) for both manual and computer measurements is given by [34]: In this case, (ff ) = (1− 0.089 h wind + 0.1166 h 2 wind ε ap ) (1 + 0.078661N), (CC) = 520(1 − 0.000051 ϕ 2 ), (ee) = 0.430(1 − 100/T pm ).
The formulas below used to estimate the back and edge heat loss (Q b and Q e ) [33,34]: where (L b and k b ) refer to the thickness and thermal conductivity of back-insulation, respectively. (k e , L e, and A e ) indicate the thermal conductivity, thickness, and area of the collector edge-insulation.

Collector Efficiency Factor
Collector Efficiency Factor (F) signifies the temperature distributions in an absorbing sheet at a right angle to the fluid flow. The configuration of the plate-tube collector absorber is exhibited in Figure 3. Through solving the problem of the standard fin, (F) can be obtained by assuming an insignificant temperature tendency in the flow direction. The plate's temperature is believed to be (T b ), and the centerline between two riser tubes is (W-D)/2. The fin formula can be solved by adding the energy balance in the direction of flow into an elementary area of (dx) width and unit length [33,34]: Here, (Ú L ) is the total loss coefficient by the aperture size. To simplify this 2nd-order differential equation, the following necessary boundary conditions must be symmetry at the centerline as well as the plate's temperature.
The energy transmitted by conduction through the flow-direction of the fin-base can be evaluated by:q where (F) is the rectangular fin efficiency and given by Equation (14): Here, (m) is the variable arrangement of fin-air and defined as m = Ú L kδ . The collector, useful energy gain, covers the heat collected above the pipe area: The useful gain in the flow-direction for the pipe and fin-per-unit is the summation of (q f in andq tube ):q The useful gain is transported to the operating liquid. The heat-flow resistance to the absorbent medium combines the absorbing medium's bond and tube resistance. The useful energy gain per unit of length can be seen in the flow-direction: where, (T b and T f ) are the base and fluid temperatures, (h fi ) is the internal heat transfer coefficient, (D i ) is the inner diameter of the riser pipe, (C b ) is the bond thermal conductance. The useful gain in the flow direction for the tube-fin per unit of length is the summation of Equations (16) and (17) where (F) refers to the collector efficiency factor: The consequence of Equation (19) for the experiment is a physical description of the collector performance factor. The local fluid temperature (F ) represents the relation of the actual energy gain to the useful energy gain in a specific area.

Heat Removal and Flow Factors
The FPSC-heat removal factor (F R ) is a quantity that relates to the actual beneficial energy gain of a collector to the useful gain if the whole collector surfaces were at the fluid inlet temperature [33,34]: Meanwhile, the flow factor (F ) refers to the relation between the heat removal factor and the FPSC-efficiency factor:

Mean Fluid and Flat Plate Temperatures
Knowing that the energy loss coefficient (U L ) and the internal heat transfer coefficient (h fi ) are essential for calculating the heat collecting efficiency, the mean fluid temperature can be described as [33,34]: One of the critical properties of the absorbing medium is the evaluation of an appropriate temperature. Owing to the heat transfer resistance between the collector and the working medium, the mean fluid temperature (T fm ) is lower than the mean plate temperature (T pm ). (T pm ) and can be measured with Equation (24):

The Solar Collector Efficiency
Instantaneous efficiency is a useful metric to characterize the heat collected performance. Instantaneous solar collector efficiency η i represents the ratio of solar energy collected divided by the solar energy available [33,34]: where I T is the incident solar radiation intensity, the instantaneous efficiency becomes: where the absorbed energy S c is changed to S c = I T (τα) according to the gross absorbed area. MATLAB package was applied to develop the mathematical model algorithm for the FPSC. The algorithm flowchart is shown in Figure 4.  Figure 5 displays the Ultraviolet-Visible spectroscopy for covalently functionalized GrNPs suspended in a base fluid (distilled water). It can be shown that in the wavelength range of 260-270 nm, the peak absorption is produced by the existence of GrNPs in all of the synthesized nanoparticles. The maximum absorption peak indicated at~264 nm is possibly due to the π→π* transition of C=C bonding for the graphene [35,36]. Following this peak and within the spectrum of wavelengths revealed in Figure 5, a reduction in absorbance was noticed in all tests at a minimum peak of~974 nm [37]. Zeta potential measurement (surface charge analysis) is a technique used for determining the surface charge of GrNPs in a colloidal solution. Charged particles may exhibit various electrokinetic effects under the influence of an electrical field, including electrophoresis, electroosmosis, streaming potential, or settlement potential [38,39]. The electrophoretic performance and additional details were performed by zeta potential evaluation to consider the water based GrNPs nanofluids' dispersion behavior. Zeta-potential is one of the standard approaches for determining the long-term stability of dispersions; the solution is considered stable if the zeta-potential value exceeds (−30 mV) or exceeds (+30 mV) [40]. Figures 6 and 7 demonstrate the examination of zeta potential and nanoparticle size distributions for the GrNPs-H 2 O nanofluids. At pH = 7, Figure 6 shows the zeta-potential and the polydispersity index (PDI) for GrNPs-H 2 O nanofluid. To create a natural repulsive force between the nanoparticles, the zeta potential must be as high as possible (+/−) to form a traditional repulsive force between the GrNPs. It can be seen that GrNPs show a negative charge of about (−39.6 mV) in the test condition of 25 • C and one hour sonication. To calculate the average size of the synthesized GrNPs, a nanoparticle size-distribution analysis was performed via the principle of dynamic light scattering (DLS) in aqueous colloidal dispersions. The average particle size of GrNPs was 289.6 nm, as illustrated in Figure 7. DLS results showed that the scale ranged from 71.4-953.7 nm with a low polydispersity index of 0.248. This implies that the solution had the most consistent and uniform distribution of particle size. Figure 8 depicts the representative scanning electron microscope (SEM) micrographs of GrNPs. SEM is a valuable method for distinguishing the size-scale, distribution, and shape of nanomaterials. SEM micrographs visualize the dispersion and stabilization effect of GrNPs-H 2 O nanofluid. Ultra-Flat Silicon Wafers is an excellent substrate for imaging, experiments, and microfabrication applications to obtain SEM characterization. Figure 8 shows that nanoparticles do not accumulate, so they appear in an aggregation mechanism and are well distributed. SEM micrographs also reveal that the covalent synthesizing of Gr actively allows wrinkled structures. Figure 8 shows the spectral study of GrNPs with energy-dispersive X-ray (EDX). The EDX measurements show a carbon content of 94.03% and an atomic oxygen content of 5.86%. While the atomic content of Si and S is 0.10% and 0.01%, respectively. The findings presented demonstrate the high quality of the shaped nanomaterial and agree well with other reported data [41][42][43].   [27,29]. The dynamic viscosity, thermal conductivity, and density were directly proportional to the percentage increase in GrNPs, while the specific heat capacity demonstrated different characteristics. Relative to the base fluid (DW) data, the average thermal conductivity changes were 20.70%, 25.87%, and 30.40% for 30 • C, 40 • C, and 50 • C, respectively, at 0.1 wt.% The corresponding highest dynamic viscosity rises were 38.30%, 36.01%, and 42.44% for 30 • C, 40 • C, and 50 • C, respectively, at 0.1 wt.% Slight changes in the density were observed at 0.14%, 0.10% and 0.11% for 30 • C, 40 • C, and 50 • C, respectively, at 0.1 wt.% While the specific heat capacity results show decreases of 1.47%, 1.50% and 1.56% for 30 • C, 40 • C, and 50 • C, respectively, at 0.1 wt.% Table 2. Thermophysical properties of Gr-nanofluids at various mass fractions and temperatures [23,27].

Thermal Behavior of FPSC-Based Water
Distilled water was tested first inside the indoor-FPSC experimentations to validate the collected data. The setup and MATLAB package data were compared throughout the current study. The maximum and average errors were 4.56% and 1.95%, respectively, between the MATLAB and experimental data. The heat collected efficiency is experimentally and theoretically displayed in Figure 9 against the dropped temperature coefficient ((T in −T amb )/G T ) at multiple mass flow rates. As the water flow rate increased from 1 to 1.5 kg/min experimentally and theoretically, the heat collected energy increased between (1.81% to 2.42%) and (1.8% to 2.41%). Therefore, it was also apparent that the rise in water mass flow rates resulted in the FPSC-efficiency improvement and decreased the absorber's surface temperature. The overarching explanation was impacted by the collector's reduced heat loss, thus increasing efficiency. F R (τα) and F R U L are also calculated by drawing the collected data with curve fittings to present the thermal behavior of the FPSC. Table 3 gives the values of F R (τα) and F R U L for numerous water mass flow rates. Once the water flow rate had risen from 1 to 1.5 kg/min, the value of F R (τα) increased by about 2.4% and 2.1% for the experimental and theoretical approaches. Simultaneously, the experimental and theoretical decrements for the value of F R U L were 21.44% and 15.14%, respectively.

Thermal Behavior of FPSC-Based Nanofluids
Owing to the high ratio of surface area to volume, suspending solid nanoparticles in distilled water will not result in a simple mixture. The nanoparticles will tend to accumulate over time due to high surface tension between them. The nanoparticles will settle and block the flow channels, and the thermal conductivity will be reduced as a result of this agglomeration. To synthesize a nanofluid with high thermal conductivity and a high heat transfer coefficient, the lowest possible concentration of nanoparticles with relatively high thermal conductivity should be used. This method is critical because nanofluids with lower concentrations of dispersed nanoparticles are known to be more stable. Viscosity would be higher in nanofluids with higher nanoparticle concentrations. The higher the viscosity of the nanofluid compared to the base fluid, the higher the pressure drop, which is directly proportional to the pumping power needed. Figure 10a-c shows the experimental and theoretical data for base fluid and nanofluids at different percentages to describe the heat collected efficiency against ((T in −T amb )/G T ) for multiple GrNPs-H 2 O flow rates. The energy performance increased when the GrNPs nanomaterials were suspended in the distilled water at different nanoparticle percentages. Relative to DW, the experimental data of the FPSC-efficiency increased by 10.4%, 10.9%, and 13% for 0.5, 1, and 1.5 kg/min, respectively. Furthermore, the theoretical efficiency improved by 9.9%, 10.4%, and 12.5% for the case of 0.5, 1, and 1.5 kg/min. The maximum and average errors between the theoretical and experimental results were 4.03% and 1.74%, respectively.  The current findings are consistent with the previous reports of Karami et al. [44][45][46]. The rise in the weight percentages of Gr-H 2 O increased the input energy absorption, resulting in a higher efficiency enhancement. The heat supplied was absorbed uniformly over the nanofluid layers for the lower mass fraction of Gr-H 2 O. Consequently, the amount of thermal loss at the boundaries was smaller than the maximum heat absorption in the top layer of the nanofluid for the high mass fraction of GrNPs. This resulted in a region of higher temperatures near the top wall, which reduced the heat gained efficiency due to increased heat losses [47]. Table 4 provides a comparative summary of the earlier investigations of FPSC performance using nanofluids with carbon nanostructures and the outcomes of the present results. Table 3 displays (F R (τα)) and (F R U L ) for Gr-H 2 O nanofluids. At the test conditions of 0.5 kg/min and 0.1 wt.%, (F R (τα)) were experimentally and theoretically enhanced by 10.89% and 10.38%, respectively. The heat loss coefficient (FRUL) values experimentally and theoretically increased by 10.97% and 6.09%, respectively. In the same order for 1 kg/min, the experimental increases in the (F R (τα)) and (F R U L ) were 11.39% and 8.06%, respectively. Besides, the theoretical increments in the (F R (τα)) and (F R U L ) were 10.86% and 2.77%, respectively. Further, Table 3 indicates the experimental increases in the values of (F R (τα)) and (F R U L ) which were 12.2% and 7.28%, respectively, at 1.5 kg/min. The MATLAB data enhancements for the heat gain and the heat loss coefficients were 11.67% and 7.58% for the same condition.

New Correlations of Thermal Efficiency
A new thermal efficiency correlation was developed as a function of the dropped temperature coefficient ((T in −T amb )/G T ) and various mass percentages at multiple flow rates, as shown in Equations (27) and (28). A linear model with 95% confidence bounds using MATLAB package was generated to predict the coefficients of FPSC-based H 2 O and Gr-H 2 O nanofluids. As shown in Figure 11a,b, the average and maximum deviations between the predicted and experimental data were (±2.653% and ±3.353%) and (±2.758% and ±4.232%) for H 2 O and Gr-H 2 O nanofluids, respectively. The coefficients of the new correlation and regression statistics are presented in Table 5. Table 3 displays (FR(τα)) and (FRUL) for Gr-H2O nanofluids. At the test conditions of 0.5 kg/min and 0.1 wt.%, (FR(τα)) were experimentally and theoretically enhanced by 10.89% and 10.38%, respectively. The heat loss coefficient (FRUL) values experimentally and theoretically increased by 10.97% and 6.09%, respectively. In the same order for 1 kg/min, the experimental increases in the (FR(τα)) and (FRUL) were 11.39% and 8.06%, respectively. Besides, the theoretical increments in the (FR(τα)) and (FRUL) were 10.86% and 2.77%, respectively. Further, Table 3 indicates the experimental increases in the values of (FR(τα)) and (FRUL) which were 12.2% and 7.28%, respectively, at 1.5 kg/min. The MATLAB data enhancements for the heat gain and the heat loss coefficients were 11.67% and 7.58% for the same condition.

New Correlations of Thermal Efficiency
A new thermal efficiency correlation was developed as a function of the dropped temperature coefficient ((Tin−Tamb)/GT) and various mass percentages at multiple flow rates, as shown in Equations (27) and (28). A linear model with 95% confidence bounds using MATLAB package was generated to predict the coefficients of FPSC-based H2O and Gr-H2O nanofluids. As shown in Figure 11a,b, the average and maximum deviations between the predicted and experimental data were (±2.653% and ±3.353%) and (±2.758% and ±4.232%) for H2O and Gr-H2O nanofluids, respectively. The coefficients of the new correlation and regression statistics are presented in Table 5.

Conclusions
In conclusion, experimental and theoretical methods were performed to report the heat collected performance of the indoor configuration of the FPSC-system using GrNPs-H 2 O nanofluid. This work was carried out under different testing conditions, namely different weight concentrations, various inlet fluid temperatures, multiple flow rates, and changing heat flux intensities. Ultraviolet-Visible spectroscopy, zeta potential, nanoparticle size distribution, and SEM-EDX were used as characterization instruments to analyze covalentfunctionalized graphene nanomaterials. The analysis led to the following conclusions: (1) The maximum absorption peak was indicated at~264 nm, a reduction in the absorbance was noticed in all tests at a minimum peak of~974 nm. GrNPs showed a negative charge of about (−39.6 mV) in the test condition of 25 • C and one hour sonication. (2) The EDX measurements showed a carbon content of 94.03% and an atomic oxygen content of 5.86%, while the atomic contents of Si and S were 0.10% and 0.01%, respectively. (3) The maximum experimental and theoretical heat collected efficiency was achieved at 76.13% and 79.2% at 0.1 wt.% and 1.5 kg/min, which was 13% and 12.5% more than DW at the same conditions. (4) The overall and average errors between the theoretical and experimental values were 4.03% and 1.74%, respectively. (5) The maximum experimental and theoretical rise in the (F R (τα)) was 12.2% and 11.67%, respectively, for 1.5 kg/min and 0.1 wt.%, while the corresponding experimental and theoretical values of (F R U L ) increased by 10.97% and 7.58% at 0.5 and 1.5 kg/min, respectively, at 0.1 wt.%.  Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Data available on request due to restrictions e.g., privacy or ethical.
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the need for further research.

Conflicts of Interest:
The authors declare no conflict of interest.