Abstract

There are many alternative gas-fired combi boilers that can be used to heat residential homes. Evaluation and selection of gas-fired combi boilers for buildings is an intricate multi-criteria decision-making (MCDM) problem involving perhaps contradictory quantifiable and qualitative criteria. In this research, as the MCDM approach, hesitant fuzzy linguistic analytic hierarchy process (HF-AHP) and hesitant fuzzy linguistic “multiple objective optimization based on ratio analysis plus full multiplicative form (MULTIMOORA)” (HF-MULTIMOORA) are integrated to assess and rank combi boiler alternatives for buildings. First, with HF-AHP, fuzzy criteria weights are determined and then with HF-MULTIMOORA, boiler alternatives are ranked from best to worst. In this integrated HF-AHP-MULTIMOORA method, evaluations of decision-makers are combined with fuzzy envelope approach and then triangular fuzzy numbers are utilized. For comparison analysis, HF-AHP-TOPSIS method is also applied to the same problem. A case study in Turkey is presented where ten combi boiler alternatives are assessed based on fifteen criteria by five decision-makers. We have used various selection criteria for boilers ranging from maximum temperature, heating capacity up to environmental effects and decided on the best combi boiler for heating residential buildings in Turkey.

1. Introduction

Buildings and activities in buildings play a very critical role in energy consumption around the world. The amount and quality of energy used in the construction industry have a great impact on environmental and global activities, consumption of natural resources, emissions to the environment, and even human health. Building energy consumption continues to increase in recent years due to rapid population growth, global climate changes, and increasing demand for a healthy, comfortable, and efficient indoor environment. Nowadays, buildings are responsible for over one-third of global final energy consumption and around 40% of CO2 emissions.

Most of the energy use in buildings is to provide heating, ventilation and air conditioning (HVAC). Natural gas is one of the methods used in air conditioning and has become more popular recent years. Also natural gas consumption is one of the most critical sources of energy saving in building energy consumption.

After the initial use of natural gas, the heating sector met with combi boilers and the use of combi boilers became widespread. The combi boilers that have been produced since 1992 have been diversified according to the developing technology, today’s needs, and customers’ requests. Mistakes made in the selection of combi boiler devices, which have many differences in terms of their features, bring about heating problems; moreover, they may negatively affect fuel consumption.

Combi boilers are used to supply heat and hot water to households. These types of boilers are very popular in Turkish residential buildings because they are effective and take up little room. These boilers are usually placed on kitchen or bathroom walls. Combi boilers are energy efficient because they do not store water, and they heat water for baking and bathing as an open system and space heating as a closed system.

Another characteristic of the recent combi boilers that makes it more effective than conventional boilers is that the flue gas produced during the heating process is not completely lost with the combi boilers, as the gas is extracted and reused (condensation technique). Combi boilers do not require complex installation, and they are convenient to be installed. There is no dedicated tank for hot water as water is heated directly and used instantly. Combi boilers are directly connected to the main water source to supply hot water.

Condensation technique for combi heating was a major step as it reduces gas consumption in apartment houses and independent houses. Moreover, sensitive heat losses are significantly reduced due to the resulting low temperature of the flue gas. Condensing boilers are aimed to draw heat out of a given quantity of fuel as the temperature of the flue gases is reduced to a point where water vapor produced during combustion is condensed, releasing latent heat that would otherwise escape with the flue gases. Since the 1970s, condensing boilers have gone through many design and installation changes to improve their performance.

Natural gas has been used in Turkey since 1988. Around 25,000 m3 natural gas is consumed per capita for individual heating in addition to central heating of the residential buildings. The individual heating in the form of combi boilers is commonly used for space and water heating. The life of the combi boiler unit is around 10 years. The two main factors are important for combi boiler selection for houses. The first one is the economic one that directly depends on the fuel consumption of combi boilers. The second one is related to the heating capacity of boilers that usually depends on the size of the house that needs to be heated. Combi boiler systems are installed to produce both hot water and space heating. Therefore, the maximum temperature of the hot water is an important criterion.

Condensing boilers utilize the latent heat contained in the flue gas. Therefore, the boiler efficiency is another important criterion. In the condensing combi boiler, the temperature of the heating medium must be as low as to create some temperature difference between the flow pipe and the return temperature. Therefore, the temperature of the heating medium and outdoor temperature are other important criteria. The supply and return temperatures vary according to the ambient temperature.

Heating requirements of residential homes, in particular space heating of residential homes, are generally a significant component of energy consumption in Turkey. Several heating appliance manufacturers are concerned about the search for fuel efficient boilers. Decisions on how to heat residential homes include many technological alternatives that involve many conflicting criteria. First of all, high life and capital costs are important economic indicators. Second, the solutions that are most beneficial for both capital and lifecycle costs cannot be implemented for noneconomic reasons. Another issue is decision-making: costs represent quantitative information, but there is also qualitative information that influences final choices. In order to make decisions where there are multiple criteria, a number of formal methods, known as MCDM methods, are available. A recent review of these methods and their application in the energy and building field can be found in [1]. Our study aims to show, through a case study, how MCDM can enable a conscious selection from among the combi boiler options for residential heating. In our paper, we will point out some distinctive features of choice of heating systems in a residential context.

In this study, as the MCDM method, HF-AHP is integrated with HF-MULTIMOORA. Here, first, HF-AHP is utilized to calculate the fuzzy criteria (importance) weights since HF-MULTIMOORA does not include a step to determine (fuzzy) weights and assumes readily available (fuzzy) weights. With HF-AHP, as a result of pairwise comparisons of criteria, consistent, therefore reliable fuzzy weights are obtained. Afterward, HF-MULTIMOORA is implemented to rank combi boiler alternatives since MULTIMOORA and its extensions are appraised as the most robust MCDM method with a wide range of real-life applications [24]. Here, in both HF-AHP and HF-MULTIMOORA methods, “hesitant fuzzy linguistic term sets” and “hesitant fuzzy set concepts” are employed with “fuzzy envelope approach” to “reflect the uncertainty, ambiguity, and hesitations DMs might have in their assessments” [5, 6] of criteria and combi boiler alternatives. Thus, the decision-making process becomes more flexible and “close to real” [5, 6]. As the novelty of the research, to the best of our knowledge, integration of HF-AHP with HF-MULTIMOORA (HF-AHP-MULTIMOORA) has never been studied, especially for the selection of combi boilers. Therefore, this study has the aim of attaining a consensus by the DMs about the best combi boiler among the alternatives considered, applying HF-AHP-MULTIMOORA. In the following section, associated literature review is presented. In section 3, details of HF-AHP-MULTIMOORA are given, along with the case study in Section 4 and conclusions in Section 5.

2. Literature Review

The selection of heating systems for the various buildings is a new area of application of MCDM methods. MCDM has been used extensively in the energy strategy [1], but there are few applications of MCDM methods to plan residential energy systems [1]. Multi-objective optimization tools for building design [7] or heating and cooling devices [8] have only been quantitatively investigated. Multi-objective optimization models are adapted to the solution of continuous problems [9]; on the other hand, multi-characteristic methods require a limited number of choices. Solving a multi-characteristic decision-making problem can be assumed to be a selection problem [9]. So far, only a few examples of heating system decision-making problems with multiple-characteristic approaches can be found in the literature. A widely used multi-characteristic approach, called ELECTRE method [10], was applied for the selection of the heating devices by Thiel and Mroz [11]: three alternative heating devices for a building based on quantitative primary energy consumption, total cost criteria, and qualitative criteria. Only conventional devices such as fossil fuels are considered in the study. Emphasis was placed on heat distribution systems and their interaction with the building structure. Kaklauskas et al. [12] applied their multi-characteristic method (COPRAS) to residential homes, focusing mainly on building design and renovation problems (e.g., Kaklauskas et al. [13]) to guide the selection of contractors to renovate buildings. Jaber et al. [14] and Alanne et al. [15] focused on policy designs: Jaber et al. [14] used the AHP method [16] and fuzzy sets [17, 18] to evaluate conventional and alternative energy sources for residential heating systems, and Alanne et al. [15] focused on combined heat and power systems and considered uncertainties in the decisions of decision-makers by means of the MCDM method [15]. Until now, the AHP method has been widely applied in the energy strategy field, especially for conventional and renewable power generation (see [19] and [2022] for some more examples). However, the selection of combi boilers for residential heating has never been investigated in the literature from the point of view of the MCDM.

To capture the ambiguity and fuzziness in the decision-making process, fuzzy set theory concepts [18, 23] were integrated where AHP and fuzzy AHP were utilized in many MCDM applications. In this study, HF-AHP is favored instead of AHP or fuzzy AHP to determine fuzzy criteria weights since with HF-AHP, the ambiguity and vagueness on DMs’ decisions are captured, and moreover, the degree of hesitancy DMs may have reflected with the usage of hesitant fuzzy set concepts and hesitant fuzzy linguistic term sets, which are introduced in [24, 25]. Overall, all these features provide flexibility in the decision-making. HF-AHP has been applied to a varied series of MCDM problems in the literature in order to determine (fuzzy) criteria weights such as evaluation of COVID-19 intervention strategies [6], CNC routers [26], innovation projects [27], power generation initiatives [28], summer sport schools [29], transformers in a power distribution project [30], performances of bank regions [31], and a cargo company [32].

After determination of fuzzy criteria weights with HF-AHP, in this study HF-MULTIMOORA is utilized to rank alternatives. Fuzzy MULTIMOORA (F-MULTIMOORA) is acclaimed as the most robust MCDM method and it is widely applied in real-life problems [24]. The parts of MULTIMOORA, the “ratio system,” “reference point approach,” and the “full multiplicative form” were first described by Brauers in 2004 [33]. Afterward, this method was named MOORA [34] and then MULTIMOORA was developed [35], which is combination of MOORA and the “full multiplicative form.” The fuzzy extension, F-MULTIMOORA, was launched by Brauers et al. in 2011 [36]. MULTIMOORA and F-MULTIMOORA have been implemented in many MCDM problems such as assessment of efficient farming types (F-MULTIMOORA) [37], decisions about bank loans to buy property (MOORA) [36], evaluation of performance of smart-bike sharing programs (F-MULTIMOORA) [4], comparing EU member states in reaching Lisbon Strategy 2000–2008 goals [36] (F-MULTIMOORA), evaluation of personnel (F-MULTIMOORA with group decision-making) [38], selection of sustainable energy crop [39], and evaluation of innovative ability of universities (MULTIMOORA) [40]. An extensive overview of MULTIMOORA with theory and applications is given in Hafezalkotob et al.’s study [41].

Extension of MULTIMOORA with hesitant fuzzy sets was proposed by Li in 2014 [42]. Afterward, Liang et al. [43] developed the “dual hesitant fuzzy extended Bonferroni mean”-based MULTIMOORA and applied it to the renewable energy technology selection problem. Gou et al. [44] utilized the double hierarchy hesitant fuzzy linguistic MULTIMOORA to select the best city in China in terms of implementation of haze control measures. Liao et al. [45] applied unbalanced hesitant fuzzy linguistic MULTIMOORA to an investment problem regarding the shared bicycles. Later, Lia et al. [46] implemented the hesitant fuzzy linguistic Choquet integral-based MULTIMOORA to develop human resources and select talents. Zolfaghari and Mousavi [47] applied interval-valued hesitant fuzzy linguistic MULTIMOORA to rank healthcare system failures. Saraji et al. [48] worked on “stepwise weight assessment ratio analysis (SWARA)-MULTIMOORA” with hesitant fuzzy sets to assess defies of online education during COVID-19 and rank higher education institutions.

In the literature, there is no study focusing on integration of HF-AHP with HF-MULTIMOORA (HF-AHP-MULTIMOORA). Moreover, a MCDM method has never been implemented to combi boiler selection for residential heating. With HF-AHP-MULTIMOORA, reliable fuzzy criteria weights and a robust ranking are obtained. Here, first, HF-AHP is utilized to calculate fuzzy criteria weights, and afterward, HF-MULTIMOORA is used to rank combi boiler alternatives, applying acquired fuzzy weights. In the next section, details of the HF-AHP-MULTIMOORA are given.

3. HF-AHP-MULTIMOORA

3.1. Definitions

In the proposed HF-AHP-MULTIMOORA and HF-AHP-TOPSIS (applied for comparison analysis), due to its practicality, triangular fuzzy numbers (TFNs) are employed. A fuzzy number is a special fuzzy set . Here, is a continuous mapping from R to [0, 1]. A TFN , where , has the triangular-type membership function.

Basic operations between two positive TFNs , are given as [6, 27, 30, 49, 50]

and are TFNs (not essentially positive).

(graded mean integration approach [51])

Hesitant fuzzy set (HFS) concept determines the membership degree of an element when a DM has hesitation in evaluations, since there might be more than one potential value [25, 52]. The “membership degree” of an element to a given set is articulated by numerous likely TFNs in triangular fuzzy HFS.

If is a fixed set, the HFS on returns a subset of . Here, is the possible membership degrees of element to set by taking values in .

The lower and upper bounds are determined as

HFS basic operations with h1 and h2 are presented as

An ordered weighting averaging (OWA) operator can be applied as below:

, is the jth largest of [53],

In the presented HF-AHP-MULTIMOORA, DM assessments are combined with fuzzy envelope approach [6, 27, 30, 54]. Here, and are the lowest and highest scales for the assessments. If the assessments are between and , then . Based on the hesitant fuzzy linguistic term sets, linguistic expressions can be represented by a , where d, e, and f are

Weight vector in OWA operator is defined as [55]

Here, b is the number of terms in DM’s assessment scale (in Tables 1 or 2), j is the rank of the highest, and i is the rank of the lowest assessment value. i and j can be ranks from 0 to b and n=j-i. [53].

3.2. HF-AHP

In HF-AHP-MULTIMOORA, first fuzzy criteria weights are calculated with HF-AHP. DMs compare criteria (pairwise) with the linguistic terms in Table 1 [6, 30, 56].

Steps of HF-AHP [6, 27, 30] are presented as(1)For DMs’ assessments with hesitant linguistic terms, fuzzy envelope approach is applied with equations (11)–(13) and corresponding TFNs are obtained. Since there are K DMs, are calculated. Here, is the corresponding TFN of the Kth DM. There are n criteria, so nxn matrix with elements is normalized and with elements is determined. Fuzzy criteria weight vector is calculated by taking the average of the entries at each row of .(2) is defuzzified with equation (5), and is the normalized principal eigenvector of X. Principal eigenvalue is obtained with .(3)“Consistency index” (CI) is calculated with , and “consistency ratio” (CR) is CR = CI/RI, where RI is “random index.” If CR < 0.1, the comparisons are “consistent” [57].

obtained with HF-AHP is then utilized in HF-MULTIMOORA.

3.3. Hesitant Fuzzy Linguistic MULTIMOORA (HF-MULTIMOORA)

In HF-MULTIMOORA, DMs assess alternatives with respect to each criterion with the linguistic terms in Table 2 [6, 30, 56].

Steps of HF-MULTIMOORA [34, 36, 41, 42, 58, 59] are given as follows:(1)For DMs’ assessments with hesitant linguistic terms, fuzzy envelope approach is applied with equations (11)–(13) and corresponding TFNs are calculated. Determine where is the corresponding TFN of the Kth DM. Since there are m alternatives and n criteria, mxn fuzzy evaluation matrix with positive TFN elements is determined.(2) is normalized with the following equations [60, 61]:where is the normalized version of the . Here, .(3)Ranking based on fuzzy ratio system (RS): To compute the utility of the fuzzy RS, fuzzy weighted normalized ratings are added for benefit (maximization) criteria and deducted for cost (minimization) criteria for each alternative i. is defuzzified with equation (5), and is determined for each alternative i. Based on the fuzzy RS, the best alternative is the one with the maximum utility . The ranking is obtained for alternatives A1, A1, …Am in descending order as(4)Ranking based on fuzzy reference point approach (RP) :RP is based on Tchebycheff min-max metric [36, 37, 41]. In fuzzy RP, first, the maximal objective reference point vector is determined asThen, ranking of alternatives is done based on deviation from the reference point vector and this is determined with the Tchebycheff min-max metric asIn equation (20), distance (deviation) is calculated with the vertex method (equation (7)). Based on the fuzzy RP, the best alternative has the minimum . The ranking is obtained for alternatives A1, A1,…,Am in ascending order as(5)Ranking based on fuzzy full multiplicative form (FMF): Utility of fuzzy FMF is determined by dividing the product of fuzzy weighted normalized alternatives ratings on maximization criteria by the product of fuzzy weighted normalized alternative ratings on minimization criteria as follows:is defuzzified with equation (5), and is obtained. Based on the fuzzy FMF, the best alternative is the one with the maximum . The ranking is obtained for alternatives A1, A1, …Am in descending order as(6)After obtaining ranking of alternatives with these methods (), these rankings are fused to determine the final ranking of alternatives. Here, improved Borda rule (IMB) is used for this since it is superior to dominance theory in the sense that it integrates both utilities and rankings [41, 59].

First, vector normalization [34, 41] is used to normalize values as shown below:

Then, assessment value of IMB () for each alternative i is calculated aswhere are the rankings of fuzzy RS, RP, and FMF [41, 59].

As the final ranking of alternatives (), the best alternative is the one with the maximum . The ranking is obtained for alternatives A1, A1, …Am in descending order as

Summary of the steps of HF-AHP-MULTIMOORA is given in Figure 1.

4. HF-AHP-TOPSIS

Steps of HF-AHP-TOPSIS are as follows [30, 62]:(1)After fuzzy evaluation matrix is determined (in Step 1 of HF-MULTIMOORA), elements of are defuzzified with equation (5) and with elements is obtained. “Normalized decision matrix” D with elements is calculated. Given (calculated in Step 2 of HF-AHP), “weighted normalized decision matrix” E is obtained with elements .(2)The “positive ideal” and “negative ideal” solutions are determined. where . , .(3)The “separation measures” of ach alternative are calculated. “Separation from ideal solution is” , and “separation from negative ideal solution” is .(4)“Closeness index values” for each alternative are calculated as . Alternatives are ranked based on (highest is best).

5. Case Study

In this research, 10 combi boiler alternatives that are frequently used in Turkey (A1, A2,…, A10) are evaluated with respect to 15 benefit (maximization) criteria (C1, C2, …, C15) by 5 decision-makers (DM1, DM2, …, DM5). Evaluation criteria are listed in Table 3 with explanations. The combi boiler alternatives that are ranked in this study are Vaillant ecoTEC Pro VUW 286/5-3 (A1), Baymak Brötje Novadens 42 Fi (A2), E.C.A. Proteus Premix 28 HM (A3), Baymak Brötje Novadens 24 Fi (A4), Demirdöküm Nitromix P 28 (A5), Viessmann Vitodens 050-W 24 (A6), Buderus Logamax Plus GB012-25K V2 (A7), Buderus Logamax Plus GB122i 24KD (A8), E.C.A. Calora Premix 24 (A9), and Bosch Condens 2200i W 24 (A10). DMs are thermal engineering graduate students from Istanbul Technical University (Turkey) and Middle East Technical University (Turkey).

In HF-AHP-MULTIMOORA, fuzzy criteria weight vector is determined with HF-AHP. Here, DMs’ pairwise comparisons of criteria with the linguistic terms in Table 1 are presented in Table 4.

For DMs’ assessments with hesitant linguistic terms, fuzzy envelope approach is applied with equations (11)–(13) and corresponding TFNs are obtained. Afterward, average of 5 DMs TFNs is taken and fuzzy evaluation matrix in Table 5 is determined. is obtained by normalizing . Then, in Table 6 is calculated by taking the average of each row of .

is defuzzified with equation (5), and is calculated. , , random index (RI) = 1.59 for n = 15, and CR = CI/RI = 0.065/1.59 = 0.041 < 0.1, so the comparisons are consistent.

obtained with HF-AHP is then utilized in HF-MULTIMOORA to rank combi boilers. DMs’ assessments of alternatives with respect to each criterion using the linguistic terms in Table 2 are presented in Table 7.

For DMs’ assessments with hesitant linguistic terms, fuzzy envelope approach is applied with equations (11)–(13) and TFNs are obtained. Then, 5 DMs TFNs are aggregated by taking average and fuzzy evaluation matrix in Table 8 is determined. Afterward, is normalized with equations (14)–(16) and Table 9 is obtained.

In steps 3, 4, and 5 of HF-MULTIMOORA, alternatives’ rankings based on the fuzzy RS, fuzzy RP, and fuzzy FMF are determined, based on , , and values, respectively. For and , bigger values are better, and for smaller values are better in terms of rankings. In the last step, final ranking of alternatives is obtained applying the improved Borda rule. Alternatives are ranked based on (higher values are better). Consequently, based on HF-AHP-MULTIMOORA, alternatives are ranked from best to worst as Bosch Condens 2200i W 24 (A10), Baymak Brötje Novadens 24 Fi (A4), Buderus Logamax Plus GB012-25K V2 (A7), Buderus Logamax Plus GB122i 24KD (A8), Baymak Brötje Novadens 42 Fi (A2), Viessmann Vitodens 050-W 24 (A6), Vaillant ecoTEC Pro VUW 286/5-3 (A1), E.C.A. Calora Premix 24 (A9), E.C.A. Proteus Premix 28 HM (A3), and Demirdöküm Nitromix P 28 (A5). These results are presented in Table 10.

For comparison analysis, HF-AHP-TOPSIS is implemented to the same problem. First, fuzzy evaluation matrix in Table 8 is defuzzified with equation (5) and normalized. Then, after multiplication with the weights (), “weighted normalized decision matrix” E is determined as in Table 11. “Positive and negative ideal solutions” () are also presented in Table 11. “Separation measures” and “closeness index of each alternative” and ranking of alternatives based on values are given in Table 12. As shown in Tables 10 and 12, HF-AHP-TOPSIS ranking is the same as the ranking obtained in fuzzy full multiplicative form (FMF) of HF-AHP-MULTIMOORA and the best alternative is determined by both HF-AHP-MULTIMOORA and HF-AHP-TOPSIS as Bosch Condens 2200i W 24 (A10).

6. Conclusions

As the proposed combi boiler alternatives for building heating were so diverse, our objective was to help select the most suitable option for residential customers; therefore, our aim was not to design the system or suggest further alternatives, but to evaluate the available ones. Basically, we can group heat distribution systems into two classes, depending on the temperature of the flowing heat carrier, as high and low temperature systems. These systems are both radiant heating systems which fundamentally consist of a burner unit coupled with a heat exchanger tube where combustion flue gases flow; however, their heat capacities differ. In this study, we targeted low temperature systems.

On the technical side, yearly energy requirements can be calculated before deciding on the combi boiler by estimating the local degree days, the indoor design temperature, and an operating time. The difference between the indoor design temperature and the outdoor design temperature is also an important criterion before selecting the combi boiler. However, one still needs to consider nontechnical criteria for selection of a combi boiler for residential heating.

The goal in this research was to choose the best combi boiler for residential heating according to the technical and nontechnical criteria listed in Table 3. These criteria were elicited from the expert decision-makers through a number of interviews, firstly defining a wider set of criteria, then identifying the most important through a tentative evaluation of criteria weights by pairwise comparisons, and removing the more minor with the aim of obtaining no more than fifteen final criteria.

In the end, the following fifteen criteria were retained: 1. space and water heating capacity is the spatial heating capacity of the combi boiler; 2. maximum space and water heating temperature is the maximum space and water temperature that the combi boiler provides; 3. natural gas and electricity consumption reduction is the natural gas and electricity consumption level of the product; 4. noise level reduction is the noise level of the product; 5. NOx emission reduction is the NOx emission level of the product; 6. condensing boiler is whether the combi boiler is condensing one or not; 7. price efficiency is the price of the product compared to similar products in the market; 8. product design is the product design and its attractiveness; 9. digital level of the product is how advanced the digital level of the product; 10. security level is the gas and water security level of the product; 11. warranty period is the warranty period that manufacturing company provides; 12. reliable brand is the brand perception among customers; 13. end-user satisfaction is the customer satisfaction level; 14. customer service is the customer service provided by the manufacturing company; and 15. after-sales service is the service provided by the company after sale. With HF-AHP, the relative importance of these criteria was assessed by pairwise comparisons elicited from expert decision-makers and fuzzy criteria weights were obtained.

As a result of HF-AHP-MULTIMOORA method and HF-AHP-TOPSIS method, the best alternative is found as Bosch Condens 2200i W 24 (A10). This combi boiler has a good customer service and reputation together with a good price compared to its heating capacity. This brand is well known in Turkey, and its other household appliances have good end-user satisfaction. Its natural gas consumption is low compared to other alternatives, and its space and water heating capacity is higher compared with other alternatives.

In this paper, we showed that HF-AHP-MULTIMOORA method can be successfully used to select building energy systems by evaluating technical and nontechnical selection criteria.

6.1. Statements and Declarations

The authors have no competing interests to declare that are relevant to the content of this article. All the data needed are included in the manuscript. No funding was received for conducting this study.

Data Availability

All the data needed are included in the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.