Mohamed Y.M.Mohsen?Mohamed A.E.Abdel-Rahman?A.Abdelghafar Galahom

Abstract The possibility of utilizing thorium as a fuel in a pressurized water reactor(PWR)has been proven from the neutronic perspective in our previously published work without assessing the thermal hydraulic(TH)and solid structure performances.Therefore,the TH and solid structure performances must be studied to confirm these results and ensure the possibility of using a thorium-based fuel as an excellent accident-tolerant fuel.The TH and solid structure performances of thorium-based fuels were investigated and compared with those of UO2.The radial and axial power peaking factors(PPFs)for UO2,(232Th,235U)O2,and(232Th,233U)O2 were examined with a PWR assembly to determine the total PPF of each one.Both Gd2O3 and Er2O3 were tested as burnable absorbers(BAs)to manage the excess reactivity at the beginning of the fuel cycle(BOC)and reduce the total PPF.Er2O3 resulted in a more significant reduction to the total PPF and,therefore,a greater reduction to the temperature distribution compared to Gd2O3.Given these results,we analyzed the effects of adding Er2O3 to thorium-based fuels on their TH and solid structure performances.

Keywords Thermal hydraulic(TH)?Solid structure?Thorium-based fuel?Gd2O3?Er2O3

1 Introduction

Since the nuclear accident of Fukushima Daiichi,many efforts have been made to find an accident-tolerant fuel(ATF)that can safely operate under harsh conditions inside the reactor.This type of fuel would prevent the reactor core from melting down and causing a disaster[1,2].ATFs must have better properties than traditional fuels to perform essential tasks during accident conditions,for example,decreasing oxidation kinetics and hydrogen generation rates.

Thorium-based fuels have advantages that invite researchers to investigate the possibility of using them as AFTs,such as their high resistance to nuclear proliferation and the low quantity of long-lived minor actinides(MAs)and non-actinides in the discharged fuel from the reactor[3–8].Various forms of thorium-based fuels have been certified,such as oxides,metals,nitrides,carbides,and molten salts.These fuels have also been nominated as nuclear fuels for different types of reactors,including heavy-water reactors,light-water reactors,accelerator-driven systems,and high-temperature gas-cooled reactors.The breeding capability of Th-U fuel material in accelerator-driven subcritical molten salt reactors has already been studied,where researchers determined the conversion ratio and net production of233U and the inventories of the various MAs[9].

Thorium dioxide carries significant research interest as an alternative to the traditional fuel material of uranium dioxide in Generation IV reactors due to many reasons,including its high abundance in the earth’s crust,its effectiveness in power flattening for nuclear reactors,its usefulness as a breeder blanket to generate233U,and its capability to suppress neutron leakage from the reactor core[5,10,11].This interest also comes from its physical advantages,such as its high melting point,high corrosion resistivity,low thermal expansion coefficients,and high thermal conductivity,compared to traditional fuel materials[4,12,13].Moreover,the fission gas release is lower than in uranium dioxide,which keeps the fuel pellet away from any defects that could have side effects on the TH and solid structure performances.These advantages provide the fuel material with greater stability during severe operational conditions from the neutronic,TH,and solid structure viewpoints.Spark plasma sintering(SPS)is an example of the technologies used to produce high-density(232Th,U)O2pellets,in which electric current flows through the sample and a conductive die provides rapid sintering and low temperatures[2].

Thorium dioxide has a melting point of 3651 K and is used as a ceramic fuel pellet in nuclear reactors.Its melting point is the highest known melting point for an oxide.Thorium dioxide also has an identical fluorite structure to cerium,plutonium,and uranium dioxide.However,thorium dioxide has a higher oxidation state than uranium and plutonium dioxides,making it exceptionally stable in the presence of oxygen or water and boosting its corrosion resistance[14].The existence of uranium in thorium dioxide can occur in two ways:either,through the doping process,in which,uranium is added to the thoria in a mixed oxide fuel(MOX);or,through the transmutation of the232Th isotope into the fissile233U isotope by bombarding it with a neutron,resulting in two successive beta emissions[13].

This study has three primary sections.In the first section,we investigate and compare the TH and solid structure performances of the suggested thorium-based fuel types with those of traditional fuels;then,in the second section,we determine which of the two different burnable absorbers(BAs)(Gd2O3and Er2O3)added to the standard fuel was most suitable;and lastly,in the third section,we investigate the optimum BA for thorium-based fuel.

2 Methodology

Version 2.7 of the Monte Carlo N-Particle transport code system(MCNPX)was used with a cross-section library called ENDF/B-VII to design a three-dimensional model of the PWR fuel assembly(FA).The MCNPX code,based on the Monte Carlo method,is a multi-purpose code that can simulate almost all particles in the PWR environment.It was used to investigate the radial power distribution and specify the hot channel where the TH and solid structure analyses would be performed[15,16].Figure 1 illustrates the horizontal cross-section assembly model of the PWR,as generated by MCNPX.In total,70 000 histories per cycle,20 ineffective cycles,and 100 effective cycles were used in the neutronic computations.The boundary surfaces of the fuel assembly are reflection surfaces that reflect the neutrons inside the fuel assembly.The calculations were carried out under hot full-power conditions(HFP).

The COMSOL Multiphysics computer software,which is dependent on the computational fluid dynamics(CFD)technique,was applied to investigate the TH and solid structure performances.The rod-centered subchannel model was utilized and solved numerically in the TH study.The temperature distribution for the fuel,cladding,and coolant materials,the coolant pressure drop,and the departure from nucleate boiling ratio(DNBR)was investigated using the rod-centered subchannel model.The most prominent advantage of the CFD approach is the coupling of various types of physics,which makes the solution more realistic than the analytical approach.The coupling of heat transfer and turbulence(K-?)physics mimics the effects of coolant velocity fluctuation in the axial and radial directions,and therefore,the variation of the coolant heat transfer convection coefficient,on the temperature distribution of the fuel and cladding materials.For these reasons,we used this coupling to perform our TH analysis.The coupling of solid structure and heat transfer physics helps to investigate the effects of the thermal expansion of solid materials(fuel and cladding)on von Mises stress and total displacement.To simulate the effect of the coolant pressure on the solid materials,the turbulence and solid structure physics coupling was considered.We used the 2D axisymmetric space dimension,rather than the 3D space dimension,to simulate the TH and solid structure for the hot channel in each FA to simplify the simulation and reduce the calculation time and cost.To calculate the distribution of the DNBR,the EPRI-1 correlation for PWR was used[18];this process has been illustrated in detail in our previously published article[17].Table 1 lists the properties of the examined fuel types.In this work,the TH and solid structure performances of UO2were investigated and compared with the thorium-based fuel types.Table 2 shows the main design features and input data used for simulating the TH and solid structure performances using COMSOL Multiphysics.Table 3 presents the correlations that describe the thermophysical and solid structure properties of the applied materials.

Table 1 Characteristics of the examined ATFs

Table 2 TH and solid structure main design features

3 Results and discussion

The TH and solid structure performances were investigated to provide an integrated study on the effectiveness of using(232Th,235U)O2and(232Th,233U)O2as nuclear fuels in the PWR.Performances in TH and solid structures have been investigated at the HFP near the BOC.This state is characterized by a higher PPF relative to the PPF value during the rest of the fuel cycle.This is true because the PPF decreases as the fuel burn-up increases.The TH and solid structure performances of(232Th,235U)O2and(232Th,233U)O2,without burnable absorbers,will be discussed and compared with UO2in the following section.

3.1 Comparing thorium-based fuel with UO2 without a burnable absorber

The TH and solid structure performances of(232Th,235U)O2and(232Th,233U)O2were investigated and compared with those of UO2.To study the TH and solid structure of the proposed fuel types,it is important to analyze the radial and axial power distributions through the FA.The power distribution was determined using the MCNPX code.Figure 2 presents the radial and axial PPF distributions of the investigated fuel types,with the hot channels where the power is peaking marked by orange cells.The total PPF is calculated from the axial and radial PPFs using Eq.(1).

Table 4 lists the radial,axial,and total PPF for the investigated fuel types.It is observed that the total PPF in the case of(232Th,233U)O2is smaller than that of the other fuel materials studied.This is because the infinite multiplication factor(Kinf)at the HFP and the beginning of the fuel cycle length for(232Th,233U)O2is greater than that of the other investigated fuel assemblies[3].The normalization factor used in the multiplier card in the tallies(FM)is inversely proportional to theKinffor the same thermal power.Therefore,to produce almost the same thermal power from the three investigated fuel assemblies,the normalization factor for both the flux and power for(232Th,233U)O2is smaller than that of the other fuel materials studied.

Table 3 Correlations used for describing the thermophysical and solid structure properties for the applied materials

3.1.1 Thermal–hydraulic performance

The subchannel model based on the CFD technique was used to solve the equations of the coupled system(partial differential heat generation equation for heat transfer and Navier–Stokes equation for fluid dynamics).The rod-centered subchannel model is discussed in detail in our previously published work[18].The primary goal of the TH analysis was to obtain the DNBR distribution and the temperature distribution of the coolant,clad,and fuel.Figure 3 presents the temperature distribution of the coolant,cladding,and fuel materials,where the maximum fuel temperatures for the examined FAs of UO2,(232Th,235U)O2,and(232Th,233U)O2are 1565.21 K,1497.56 K,and 1403.49 K,respectively.The safety limits for the fuel temperatures are less than 2873 K for UO2and 3500 K for(U,Th)O2,which ensures safe operation and the impossibility of fuel meltdown at the BOC.For the maximum clad inner surface temperature,the UO2FA had a maximum temperature of 611.859 K,while for the thoriumbased fuels,the maximum temperatures were 610.29 K for(232Th,235U)O2and 606.18 K for(232Th,233U)O2.Additionally,the maximum clad temperatures indicate that the clad material meltdown is impossible because the maximum clad temperatures are less than 2123 K,the melting point of zircaloy-4.Furthermore,these temperatures ensure the impossibility of the pellet cladding chemical interactions(PCCI)because they are less than 1135 K.The EPRI-1 correlation was utilized to determine the critical heat flux(CHF),which was divided by the actual heat flux to calculate the DNBR.The DNBR was then studied to explain how far the coolant is from the nucleate boiling hazards by comparing the minimum DNBR with the PWR’s safety margin of 1.75.The minimum DNBR values were 1.767,1.765,and 1.8615 for UO2,(232Th,235U)O2,and(232Th,233U)O2,respectively,which emphasizes that the formulation of nucleate boiling is impossible.It can be noticed that the temperature distributions of the investigated fuel materials are similar.This is primarily because the values of thermal conductivity for the UO2and the two suggested thorium-based fuels are almost equal,as shown in Table 3.However,the thorium-based fuels show a safe performance relative to that of UO2because of their higher melting points compared to those of the traditional fuel.There are various types of temperature changes with radius within the radial temperature distribution,where it is parabolic(1/r2)across the fuel materials.Furthermore,in both the He gap and clad materials,the temperature changes exponentially(ln(1/r)).Finally,in the coolant material,the temperature changes with 1/r.Figure 4 presents the distributions of the temperature along the channels that have the maximum PPF inside the fuel types under examination.

Table 3 continued

Table 4 Radial PPF,axial PPF,and the total PPF for the investigated FAs

3.1.2 Solid structure performance

In this section,the primary solid structure parameters(von Mises stress and the displacement of the outer fuel surface)were simulated.The von Mises stress is the total stress that affects the materials and is caused by either thermal or pressure loads.The thermal loads are due to the fission energy that is transferred to the coolant through the fuel and clad materials,which causes thermal expansion.The pressure loads were caused by the coolant pressure.If the von Mises stress is larger than the yield stress,the fuel material transfers from the elastic zone to the plastic or ductile zone.Therefore,for safe operation,the von Mises stress should not exceed the yield stress.The possibility of PCCI can be determined from the displacement of the fuel outer surface.The main causes of PCCI are the temperature of the clad materials reaching 1135,and the surfaces of both fuel and clad materials coming into contact.Under these conditions,the uranium atoms interact chemically with the zirconium atoms to form PCCI.Another type of interaction between the fuel and clad materials is called pellet clad mechanical interaction(PCMI).This type of interaction occurs when the outer surface of the fuel pellets expands thermally to reach the surface of the clad materials,resulting in excessive compression stress on the clad material,or where the clad exerts excess compression stress on the fuel material owing to its thermal expansion.Figure 5 shows that the highest von Mises pressures occurring on the investigated FAs of UO2,(232Th,235U)O2,and(232Th,233U)O2were 144.42,111.06,and 103.08 MPa,respectively.The displacements of the fuel’s outer surface for the examined FAs of UO2,(232Th,235U)O2,and(232Th,233U)O2were 0.145,0.098,and 0.089 mm,respectively.Although the displacements of the fuel’s outer surface for the examined fuel materials were greater than the thickness of the He gap,the occurrence of the PCCI is impossible because the maximum clad temperatures in the aforementioned fuel materials were less than 1135 K.The obtained solid structure results show a remarkable superiority of the two thorium-based fuel over the traditional UO2fuel material:The thermal expansion coefficient for the thorium-based fuel is lower than that for the traditional fuel material UO2,which reduces the thermal stress occurring on the fuel materials and decreases the displacement of the fuel’s outer surface.Because of the excess reactivity at the BOC that occurred in the investigated FAs of UO2,(232Th,235U)O2,and(232Th,233U)O2,it is important to use BAs to reduce this excess reactivity.Therefore,two BAs are investigated in the following section to distinguish which is more suitable from the TH and solid structure viewpoints.

Fig.2 (Color online)Radial and axial PPF distributions for the examined fuel assemblies(FAs)

3.2 Determination of the most suitable IFBA

The main goal of utilizing the BA is to suppress the excess reactivity at the BOC.This section aims to determine the most suitable integral fuel burnable absorber(IFBA)material to be applied to the investigated fuel types in an attempt to enhance their performance in both TH and solid structures,in addition to achieving suitable reactivity management.Both Er2O3and Gd2O3at a concentration of 2%were mixed with UO2.The IFBA rods are marked in this work by the green cells in the FA.As shown in Fig.6,the IFBAs were distributed near the guide tubes(coolant channels),where they occupied positions near the regions with the highest probability of neutron moderation,the highest probability of fission reactions,and hence the highest excess reactivity.The main goal of this distribution is to reduce the excess reactivity in this region and achieve fine-tuning of the power distribution across the FA.This is what happened after the use of(UO2+Er2O3),where the total PPF decreased from 1.692 for UO2(without BAs)to 1.571 and the power distribution became flatter.In contrast,after the use of(UO2+Gd2O3),the total PPF increased from 1.692 to 2.116.This is due to the large thermal neutron absorption cross section(49,000 barn)[24].The three examined FAs were designed to produce nearly the same average thermal power.Therefore,the significant decrease in the power values of the IFBA rods in the case of(UO2+Er2O3)was compensated by increasing the power values at the periphery of the FA.

Fig.3 (Color online)Axial and radial temperature distributions along the channels that have the maximum PPF inside the fuel types under examination

Fig.6 (Color online)Radial and axial PPF distributions for UO2,(UO2+Gd2O3),and(UO2+Er2O3)

Table 5 presents the maximum PPF in both the axial and radial directions and the total PPF for UO2with and without IFBAs.The total PPF in the case of UO2without IFBAs was 1.692.This value increases to 2.116 when Gd2O3is used.The opposite occurred after applying Er2O3,where the total PPF was reduced to 1.571.Therefore,the use of Er2O3not only reduces excess reactivity at the BOC,as proved in the first annex of this paper[3],but also decreases the total PPF and fine-tunes the power density distribution,which enhances the TH and solid structure performance.

Table 5 Maximum radial PPF,axial PPF,and the total PPF for the investigated FAs((UO2+Gd2O3)and(UO2+Er2O3))

3.2.1 Thermal hydraulic performance

Figure 7 presents the results of the main TH parameters,where the maximum fuel temperatures for assemblies fueled with UO2,(UO2+Er2O3),and(UO2+Gd2O3)were 1646.67 K,1585.64 K,and 1760.87 K,respectively.The maximum clad temperatures for the three investigated FAs of UO2,(UO2+Er2O3),and(UO2+Gd2O3)were 629.46 K,623.77 K,and 643.25 K,respectively.The MDNBR values for the FAs of UO2,(UO2+Er2O3),and(UO2+Gd2O3)were 1.767,1.847,and 1.530,respectively.In the case of using Gd2O3,the MDNBR value is lower than 1.75,which is the safety limit for the MDNBR according to the EPRI-1 correlation.Figure 8 displays the temperature distributions in 3D,where the temperature is maximum at the center of the fuel rod and decreases in the direction of the coolant.The results of the TH parameters indicate that the use of Er2O3decreases the total PPF,which reduces the maximum fuel and clad temperatures and increases the DNBR safety margin.Therefore,the use of Er2O3led to a more significant improvement in the thermal hydraulic performance than using Gd2O3.

Fig.4(Color online)3D temperature distribution along the channels that have the maximum PPF inside the fuel types under examination

Fig.5 (Color online)Results of the primary solid structure parameters for hot channels in the investigated FAs of UO2,(232Th,235U)O2,and(232Th,233U)O2

3.2.2 Solid structure performance

Figure 9 illustrates the variation in the main solid structure parameters(the von Mises stress operating on the fuel material and the displacement of the fuel outer surface)along the axial direction.The greatest von Mises stresses operating on the fuel materials of UO2,(UO2-+Er2O3),and(UO2+Gd2O3)are 144.42,135.70,and 158.76 MPa,respectively.Additionally,the maximum displacement of the fuel’s outer surface for the FAs of UO2,(UO2+Er2O3),and(UO2+Gd2O3)was 0.145,0.130,and 0.172 mm.The solid structure results indicate the safe operation of the three investigated FAs because the maximum von Mises stress acting on the fuel materials was less than the yield stress of UO2.Furthermore,although the displacement of the fuel outer surface is greater than the thickness of the He gap,PCCI will not occur because the maximum clad temperature in the three investigated FAs is smaller than the safety limit of 1135 K.Er2O3provided the best solid structure performance compared to UO2and the suggested FA(UO2+Gd2O3).

The results of the TH and solid structure parameters confirmed our findings from neutronic calculations,in which Er2O3gave the best performance and now provided the best and safest TH and solid structure performances relative to Gd2O3in the case of UO2.Therefore,Er2O3is applied to the investigated thorium-based fuel types in the following section.

Fig.7 (Color online)Radial and axial temperature distributions for the investigated FAs of UO2,UO2+Er2O3,and UO2+Gd2O3

3.3 Effect of Er2O3 on thorium-based fuel types

After investigating the TH and solid structure performances of the traditional FA with Er2O3and Gd2O3,the former proved to be the best IFBA for application to thorium-based fuels.Figure 10 presents the distribution of the radial and axial PPFs for the examined FAs(UO2+Er2O3,(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3).As shown in Table 6,the use of Er2O3reduces the total PPF for(232Th,235U)O2from 1.695 to 1.663 and reduces the total PPF for(232Th,233U)O2from 1.552 to 1.469.This reduction in the total PPF reflects the effective role of Er2O3in enhancing the neutronic performance,as discussed in the first annex of this work,and the performance of the TH and solid structure.The main objective of this section is to reduce the maximum fuel and cladding temperatures and increase the value of MDNBR by applying Er2O3to the two thorium-based fuel types and the standard fuel material.

Table 6 Maximum radial PPF,maximum axial PPF,and the total PPF for the investigated thorium-based FAs of(232Th,235U)O2-+Er2O3 and(232Th,233U)O2+Er2O3)

3.3.1 Thermal hydraulic performance

Fig.8 (Color online)3D temperature distributions for UO2,(UO2+Er2O3),and(UO2+Gd2O3)

Fig.9 (Color online)Results of the main solid structure parameters for the investigated FAs of UO2,(UO2+Er2O3),and(UO2+Gd2O3)

In this section,the gain in TH performance after applying Er2O3to the investigated thorium-based FAs of UO2, (232Th,235U)O2+Er2O3, and (232Th,233U)O2+Er2O3is illustrated.As shown in Fig.11,for UO2,the maximum fuel temperature decreased from 1565.21 K to 1481.99 K,the maximum clad temperature decreased from 622.07 K to 615.95 K,and the MDNBR increased from 1.767 to 1.847.For(232Th,235U)O2,the maximum fuel,clad temperatures,and the MDNBR are 1497.56 K,620.20 K,and 1.765,respectively.These values are reduced to 1474.70 K,619.01 K,and 1.787,respectively,after applying Er2O3.For(232Th,233U)O2,the maximum fuel temperature decreased from 1403.49 K to 1352.20 K,the maximum clad temperature decreased from 615.26 K to 612.51 K,and the MDNBR increased from 1.8615 to 1.919.Figure 12 presents the 3D temperature distributions for the three examined FAs(UO2,(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3).

Fig.1 (Color online)Horizontal cross section of the PWR FA using visual editor version X_24E of MCNPX

Fig.10 (Color online)Radial and axial PPF distributions for the investigated FAs(UO2+Er2O3),(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3

Fig.11 (Color online)Radial and axial temperature distributions in the investigated FAs of UO2+Er2O3,(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3

3.3.2 Solid structure performance

The gain in solid structure performance after applying Er2O3to the investigated FAs(UO2+Er2O3,(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3)is illustrated in Fig.13.For UO2,the von Mises stress and the displacement of the fuel’s outer surface decreased after applying Er2O3from 144.42 MPa and 0.145 mm to 135.70 MPa and 0.130 mm,respectively.For(232Th,235U)O2,the von Mises stress operating on the fuel material decreased from 111.06 to 109.12 MPa and the fuel outer surface displacement decreased from 0.098 to 0.095 mm.Finally,the greatest von Mises stress,acting on the(232Th,233U)O2,decreased from 103.08 to 98.74 MPa and the displacement of its fuel’s outer surface decreased from 0.089 to 0.084 mm after applying Er2O3.

Fig.13 (Color online)Results of the main solid structure parameters for the investigated FAs of UO2+Er2O3,(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3

Table 7 presents a summary of the results obtained from neutronic,TH,and solid structure analyses.This illustrates the effectiveness of each type of fuel and proves that the thorium-based fuels are a very good alternative to UO2.There are several reasons for the difference between the results of the TH and solid structure parameters of the investigated FAs,firstly,the difference in the values of the various thermophysical properties that are used in the solution of the heat generation partial differential equation,such as thermal conductivity(K),heat capacity at constant pressure(Cp),density(ρ),and heat source(total PPF),as illustrated in Table 3.Additionally,the difference in the solid structure properties of the investigated FAs is based on the variation in Young’s modulus(E),Poisson’s ratio(ν),and thermal expansion coefficient(α).The dominant factors that strongly affect the TH and solid structure results are the heat source value(total PPF)and the thermal expansion coefficient,assuming that the remaining factors are almost the same for all the investigated fuel materials.The main turbulence properties are assumed to be the same for all investigated fuel assemblies,such as the coolant mass flow rate of 16,720 kg/s,coolant inlet pressure of 15.41 MPa,turbulence model(K–ε),and wall condition(no-slip).

Fig.12 (Color online)3D temperature distribution for UO2+Er2O3,(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3

Table 7 Results of the main neutronic,TH,and solid structure parameters for the investigated FAs of UO2+Er2O3,(232Th,235U)O2+Er2O3,and(232Th,233U)O2+Er2O3

4 Conclusion

From the neutronic perspective,the results show that the use of(232Th,235U)O2and(232Th,233U)O2reduced the actinides and non-actinides nuclide inventory,which decreased the radioactivity of the spent fuel during the storage step.

From the TH perspective,the high melting point of the thorium(232Th,235U)O2and(232Th,233U)O2is one of its essential thermophysical properties.Additionally,the maximum fuel and cladding temperatures for the(Th,235U)O2and(232Th,233U)O2are lower than those for the standard fuel(UO2),while the MDNBR values of(232Th,235U)O2and(232Th,233U)O2were greater than those of UO2.

From the solid structure perspective,owing to the superiority of the solid structure properties of(232Th,235U)O2and(232Th,233U)O2over that of UO2,the maximum von Mises stress operating on the fuel and cladding materials and the displacement of the fuel’s outer surface are less than that of UO2.Thus,it can be concluded from the neutronic,TH,and solid structure perspectives that(232Th,235U)O2and(232Th,233U)O2are good alternatives for the traditional fuel material UO2.

Authors’contributionsAll authors contributed to the study conception and design.Material preparation,data collection,and analysis were performed by Mohamed Y.M.Mohsen,Mohamed A.E.Abdel-Rahman,and A.Abdelghafar Galahom.The first draft of the manuscript was written by Mohamed Y.M.Mohsen,and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.