Transient Thermal Performance of Power Cable Ascertained Using Finite Element Analysis

: This paper presents the computation of the cable ampacity and the temperature distribution through long duration based on the equivalent thermal circuit based on IEC 60287 standard and the Finite element method using COMSOL (Multiphysics environment, version 5.5). This study investigated the cable ampacity and the temperature rise of the cable core and sheath at steady state and emergency conditions. The cable ampacity was investigated at different conditions such as the variation of cable depth, soil properties, and soil temperature. The results conﬁrmed the adaptation between the thermal circuit results and the COMSOL results as well as the effectiveness of using the numerical method to compute the cable ampacity. Using the COMSOL-based thermal properties evaluations, the transient performance of the cable is ascertained. The transient study is performed for different cable loading currents and dry zone sizes.


Introduction
The installation and maintenance of power cables are more expensive than overhead lines. The huge installation cost makes it imperative to use its full capacity. The ampacity of the cable is limited based on the cable temperature, which does not exceed 90 • C for XLPE cables and controls the cable lifetime. The cable insulation and its jacket subjected to damage with the excessive cable temperature [1]. The model of the power cable ampacity computation was proposed by Neher-McGrath in 1957 [2,3]. Utilities and cable manufactures depended on the international electrotechnical commission IEC 60287 standard, which used several simplifications and had its limitations. In addition, rating computations of the insulated cables with considering the non-adiabatic heating effects are explained in detail in [4,5].
Recently, several attempts were accomplished to modify the limitations in IEC 60287 standard. Aras et al. [1], compared the results of applying finite element method (FEM) to calculate the cable ampacity with the results of IEC 60287 standard and the experimental work results. Freitas, et al. [6], presented a numerical simulation for cable ampacity computation taking into consideration the effect of soil heat and moisture migration on the cable ampacity. Hiranandani [7], solve the problem of cable rating using Finite difference method (FDM) concerning the cables touching and separation of the cables using insulation materials. Hwang [8], used FEM to compute the cable rating in ducts.
The parameters of the thermal circuit representing the heat exchange of the cable under load condition were addressed in many works. Zaytsev [9] developed a technique to compute the thermal equivalent circuit of the high voltage cable buried into ground. The thermal circuit parameters are discussed simulating the environment of the cable in trefoil and plane configuration. In addition, the prediction of the heating of the cable is developed using the proposed method to protect the cable from insulation degradation. Yang, et al. [10], calculated the conductor temperature based on the thermal sensor position for high voltage cables. Four thermal circuit parameters were considered the temperature of insulation shield, the aluminum sheath, center of the waterproof compound, and the jacket surface to compute the conductor temperature. Simulation was presented to investigate the accuracy of the temperature computation of the cable. Diaz-Aguiló and de León, [11], modeled the transient thermal behavior of mutual heating effects between the power cables. They referred that although the proposed accurate model of the thermal circuit in IEC 60287, it is difficult to implement in electrical software and not consistent for layered model of the thermal resistance. In [11], an alternative model was developed consisting of injecting the correct current at the right position of the RC circuit representing the soil. The proposed model is validated using the FEM simulation of realistic cable establishment.
Nelson, et al. [12], the thermal behavior considering the cable center as a hottest spot was analyzed based on the ladder thermal circuit method. A suitable cross section area of the high voltage submarine cable was the main objective of this work, as well as a method of dynamic rating of the cable was proposed and accomplished via iterative process. The FEM results was compared with that of the proposed method to validate the proposed technique. Xiong, et al. [13], indicated that the studied the temperature field distribution and the cable ampacity improving the guidelines of cable manufactures. Thermal circuit in IEC 60287 and numerical simulation were used to compute the cable temperature. The numerical simulation results were close to that from the actual measurement. Žiga, et al. [14], investigated the influence of environmental conditions on the core heating of reinforced aluminum conductor steel. The computations of temperature distribution in steady state were performed based on the balance of powers not energies. The results showed the reliability and the ability of the proposed model to give satisfactory results according to the actual measured cable temperature.
Due to several simplifications and limitations of the IEC 60287 standard, this paper modelled the cable, based on the equivalent thermal circuit, using COMSOL environment. First, paper used the analytical model to compute the cable ampacity under different conditions and the temperature through a cable, using the two-loop thermal circuit of the cable, based on IEC 60287 standard. Since the analytical method considers some ideal assumptions, COMSOL Multiphysics was used to model the cable, to evaluate the temperature performance of the cable and sheath, concerning thermal model at different loading current.
In this work, the ladder thermal circuit is used to compute the cable ampacity and the distribution of temperature through the cable. The effect of the soil resistivity and cable depth on the cable ampacity is investigated. In addition, the transient rating calculations can be accomplished. The results from the thermal circuit were compared with that using numerical simulation using COMSOL software, which is build based on FEM. Transient thermal performance is evaluated at different load performances covering the conditions from maximum continuous load cycle to typical one. Concerning the maximum continuous load for long period (1000 h), impact of dry zone formation around the cable on the thermal performance is investigated at different dry zone thickness. The results indicated a good correlation between the analytical and numerical results.

Computation of the Steady State and Transient Cable Rating
The cable system as in Figure 1 can be presented using the ladder thermal circuit to compute both steady state and transient rating equations of the cable [15]. These equations are suitable for a single cable including one single core or multiple cores. Figure 2 illustrates the ladder thermal circuit representing the cable system and the parameters T and Q that refer to the thermal resistance and capacitance of the different cable components, respectively. The IEC standard [16,17] computes the cable ampacity based on reducing the ladder thermal circuit into a two-loop network to simplify the ampacity computation and standardize the procedures for basic cable types. Therefore, the transient response computation of a cable under load variations can be performed using a simplified two-loop ladder thermal circuit. If the cable system is presented using the ladder thermal circuit as in Figure 3, then the equivalent two-loop thermal circuit is illustrated in Figure 4 [15]. The parameters T A is used to express T α and Q A representing Q α to maintain the correct response of short duration. The parameters T B and Q B can be expressed as follows [18,19], Processes 2021, 9,438 3 of 16 loop ladder thermal circuit. If the cable system is presented using the ladder thermal circuit as in Figure. 3, then the equivalent two-loop thermal circuit is illustrated in Figure. 4 [15]. The parameters TA is used to express Tα and QA representing Qα to maintain the correct response of short duration. The parameters TB and QB can be expressed as follows [18,19],   loop ladder thermal circuit. If the cable system is presented using the ladder thermal circuit as in Figure. 3, then the equivalent two-loop thermal circuit is illustrated in Figure. 4 [15]. The parameters TA is used to express Tα and QA representing Qα to maintain the correct response of short duration. The parameters TB and QB can be expressed as follows [18,19],     The heat flow in the ladder thermal circuit increases in steps, then, the total joule losses can be expressed as follows, where, Wc refers to the conductor losses, Ws refers to the sheath losses, Wa refers to the armor losses, λ1 refers to the sheath loss factor, and λ2 expresses the armor loss factor. The cable temperature rise (Δθ) above the ambient temperature for the ladder thermal circuit can be expressed as Δ = + 1 2 + 1 + + where, Wd is the dielectric losses, T1, T2, T3, and T4 are the thermal resistances per unit length between one conductor and the sheath, between sheath and armor, the external serving of the cable, and between the cable surface and the surrounding medium, respectively. The factor n refers to the number cable conductors. Based on the fact that Wc=I 2 R, the permissible current rating of the cable without moisture migration is developed from (3) and can be expressed as follows, where, R refers to the ac resistance of the cable. The equations of the thermal network of Fig. 5 can be solved to compute the transient rating of the power cables [15]. These equations are solved to determine the variation of the conductor temperature rise with time (θ(t)), where this temperature is not a function of I(t) as in steady state case.   The heat flow in the ladder thermal circuit increases in steps, then, the total joule losses can be expressed as follows, where, Wc refers to the conductor losses, Ws refers to the sheath losses, Wa refers to the armor losses, λ1 refers to the sheath loss factor, and λ2 expresses the armor loss factor. The cable temperature rise (Δθ) above the ambient temperature for the ladder thermal circuit can be expressed as Δ = + 1 2 + 1 + + where, Wd is the dielectric losses, T1, T2, T3, and T4 are the thermal resistances per unit length between one conductor and the sheath, between sheath and armor, the external serving of the cable, and between the cable surface and the surrounding medium, respectively. The factor n refers to the number cable conductors. Based on the fact that Wc=I 2 R, the permissible current rating of the cable without moisture migration is developed from (3) and can be expressed as follows, where, R refers to the ac resistance of the cable. The equations of the thermal network of Fig. 5 can be solved to compute the transient rating of the power cables [15]. These equations are solved to determine the variation of the conductor temperature rise with time (θ(t)), where this temperature is not a function of I(t) as in steady state case. The heat flow in the ladder thermal circuit increases in steps, then, the total joule losses can be expressed as follows, where, W c refers to the conductor losses, W s refers to the sheath losses, W a refers to the armor losses, λ 1 refers to the sheath loss factor, and λ 2 expresses the armor loss factor. The cable temperature rise (∆θ) above the ambient temperature for the ladder thermal circuit can be expressed as where, W d is the dielectric losses, T 1 , T 2, T 3 , and T 4 are the thermal resistances per unit length between one conductor and the sheath, between sheath and armor, the external serving of the cable, and between the cable surface and the surrounding medium, respectively. The factor n refers to the number cable conductors. Based on the fact that W c = I 2 R, the permissible current rating of the cable without moisture migration is developed from (3) and can be expressed as follows, where, R refers to the ac resistance of the cable.
The equations of the thermal network of Figure 5 can be solved to compute the transient rating of the power cables [15]. These equations are solved to determine the variation of the conductor temperature rise with time (θ(t)), where this temperature is not a function of I(t) as in steady state case. Van Wormer coefficient p* is expressed as where, Di refers to the insulation diameter and dc refers to the cable diameter.
In addition, the other Van Wormer coefficient p' can be expressed as, where, De refers to the cable diameter.
The internal temperature rise of the conductor due to internal parts of the cable as a function of time can be computed as follows, where, Wc is the power loss per unit length in a conductor based on the maximum conductor temperature attained. Ta and Tb refer to the thermal resistances of two-loop thermal circuit and can be expressed as follows, where, TA and TB refer to the insulation layer and sheath layer thermal resistances. QA refer to the heat capacitance of the insulation. The parameters a and b are the poles of the transfer function (θc/Wc) where a and b can be expressed as follows [20], where, The power loss is assumed to be constant during the step of the transient. Further, the attainment factor of the conductor can be expressed as follows, Based on the Kennelly hypothesis, the temperature rise at any point in the soil at any time is the sum of the temperature cause by heat sources Wt and by the fictitious image Van Wormer coefficient p* is expressed as where, D i refers to the insulation diameter and d c refers to the cable diameter.
In addition, the other Van Wormer coefficient p' can be expressed as, where, D e refers to the cable diameter.
The internal temperature rise of the conductor due to internal parts of the cable as a function of time can be computed as follows, where, W c is the power loss per unit length in a conductor based on the maximum conductor temperature attained. T a and T b refer to the thermal resistances of two-loop thermal circuit and can be expressed as follows, where, T A and T B refer to the insulation layer and sheath layer thermal resistances. Q A refer to the heat capacitance of the insulation. The parameters a and b are the poles of the transfer function (θ c /W c ) where a and b can be expressed as follows [20], where, The power loss is assumed to be constant during the step of the transient. Further, the attainment factor of the conductor can be expressed as follows, Based on the Kennelly hypothesis, the temperature rise at any point in the soil at any time is the sum of the temperature cause by heat sources W t and by the fictitious image charge placed symmetrically with the earth surface. Therefore, the temperature rise at the outer surface of the cable is then given by, where, D* = external surface diameter of cable (m), L* = axial burial depth of the cable (m), δ = soil diffusivity (m 2 /s). The Ei refers to the exponential integral which can be expressed as, For And The total temperature rise of the conductor (θ r (t)) above ambient temperature is computed as follows, where, θ c (t) refers the internal temperature rise of the conductor due to internal parts of the cable (Equation (5)), θ e (t) refers the temperature rise at the outer surface of the cable (Equation (13)), and α(t) refers to the attainment factor of the conductor (Equation (12)). Then, the total conductor temperature (θ t (t)) can be expressed as follows, Based on the cable internal parts, the internal temperature rise of the sheath can be computed as follows, where, The total sheath temperature rise above the ambient temperature can be calculated as follows, And then the total temperature of the sheath To compute the transient analysis of temperature of the cable, the conductor heat loss must be updated based on the calculated new conductor temperature then, The reference temperature (θ ref ) is considered as 20 • C. The α cu is the temperature coefficient of copper. The parameter (I 0 ) refers to the loading current or the cable ampacity.

Results of the Simulation Model
MATLAB software (R2020b/17 September 2020) was used to compute the cable ampacity under different conditions and the transient temperature of the cable for long duration based on the two-loop thermal circuit of the cable. Figure 6 illustrates the flowchart to compute the cable ampacity. The required cable data, soil data, and temperature data are used as the input to compute the R, W d , the sheath and armor factors, T, and Q. Furthermore, the type of cable configurations is identified (Flat or Trefoil), then the ampacity is computed for single and double bond (steady state and emergency). Several parameters that influence on the cable ampacity are investigated such as the cable depth, soil resistivity, and soil temperature.
pacity under different conditions and the transient temperature of the cable for ration based on the two-loop thermal circuit of the cable. Figure 6 illustrates the to compute the cable ampacity. The required cable data, soil data, and tempera are used as the input to compute the R, Wd, the sheath and armor factors, T, an thermore, the type of cable configurations is identified (Flat or Trefoil), then the is computed for single and double bond (steady state and emergency). Several p that influence on the cable ampacity are investigated such as the cable depth, so ity, and soil temperature.
On the other hand, Figure 7. illustrates the flowchart explaining the temper tribution through the cable element, which is determined based on the cable loa pacity), the cable data, and burial depth of the cable. The parameters of the twomal circuit are identified as well as the conductor losses (Wc). The temperatu cable and the sheath are now computed based on the variation of the conduc with the time (transient state). The specifications of the cable and soil are illu Table 1.  On the other hand, Figure 7. illustrates the flowchart explaining the temperature distribution through the cable element, which is determined based on the cable loading (ampacity), the cable data, and burial depth of the cable. The parameters of the two-loop thermal circuit are identified as well as the conductor losses (W c ). The temperatures of the cable and the sheath are now computed based on the variation of the conductor losses with the time (transient state). The specifications of the cable and soil are illustrated in Table 1.

Results of Cable Ampacity with Soil Resistivity Variation
The thermal soil resistivity measurement depends on the thermal constants of the medium and can be used to install and design the underground cables. The variation of the thermal soil resistivity is attributed to the variation of the moisture content of the soil. An increase of the soil moisture content results in a reduced in the thermal soil resistivity [21]. Figure 8 illustrates the influence of the thermal soil resistivity in (m· • C/W) on the cable ampacity for different cable configurations. The trefoil configuration provides cable ampacity higher than that for flat configuration at a constant soil resistivity.

Results of cable ampacity with soil resistivity variation
The thermal soil resistivity measurement depends on the thermal constants of the medium and can be used to install and design the underground cables. The variation of the thermal soil resistivity is attributed to the variation of the moisture content of the soil. An increase of the soil moisture content results in a reduced in the thermal soil resistivity [21]. Figure 8 illustrates the influence of the thermal soil resistivity in (m. o C/W) on the cable ampacity for different cable configurations. The trefoil configuration provides cable ampacity higher than that for flat configuration at a constant soil resistivity.  Figure 9 explains the variation of the cable ampacity with the variation of the soil temperature for different cable configurations. At the same soil temperature, the cable ampacity for trefoil configuration is higher than that for the flat configuration. An increase in the cable temperature corresponds to the decrease of the cable current (cable loading).  Figure 9 explains the variation of the cable ampacity with the variation of the soil temperature for different cable configurations. At the same soil temperature, the cable ampacity for trefoil configuration is higher than that for the flat configuration. An increase in the cable temperature corresponds to the decrease of the cable current (cable loading). At the same cable loading, the soil temperature for trefoil configuration of the cable is higher than that for flat configuration.

Cable Ampacity under Soil Temperature Variation
Processes 2021, 9,438 10 of 16 Figure 9. Influence of the soil temperature on the cable ampacity. Figure 10 illustrates the effect of the cable depth on the cable ampacity. An increase in the depth of the cable decreases the ampacity of the cable. The trefoil configuration provides cable ampacity higher than that for flat configuration at the same cable depth. An increase of the burial depth of the cable leads to the increase of the thermal resistance of the surrounding soil, which reduces the cable ampacity. Furthermore, at a great depth the laying cables has a significant effect on the cable rating. This effect relies on the ratio of the depth of lying to the tunnel/duct/cable diameter [19].   Figure 10 illustrates the effect of the cable depth on the cable ampacity. An increase in the depth of the cable decreases the ampacity of the cable. The trefoil configuration provides cable ampacity higher than that for flat configuration at the same cable depth. An increase of the burial depth of the cable leads to the increase of the thermal resistance of the surrounding soil, which reduces the cable ampacity. Furthermore, at a great depth the laying cables has a significant effect on the cable rating. This effect relies on the ratio of the depth of lying to the tunnel/duct/cable diameter [19].

Cable Ampacity under Varying of Cable Depth
Processes 2021, 9,438 10 of 16 Figure 9. Influence of the soil temperature on the cable ampacity. Figure 10 illustrates the effect of the cable depth on the cable ampacity. An increase in the depth of the cable decreases the ampacity of the cable. The trefoil configuration provides cable ampacity higher than that for flat configuration at the same cable depth. An increase of the burial depth of the cable leads to the increase of the thermal resistance of the surrounding soil, which reduces the cable ampacity. Furthermore, at a great depth the laying cables has a significant effect on the cable rating. This effect relies on the ratio of the depth of lying to the tunnel/duct/cable diameter [19].

Effect of the Cable Bond Configuration on the Cable Ampacity at Different Cable Temperature
There are two category ratings based on the load rating of the cable (cyclic ratings and emergency ratings). Cyclic rating is used the load varies daily over 24 h and it is the same for each day. The emergency ratings are considered when the load varies between ten minutes to one hour. Therefore, the emergency refers to the transient period when the load varies in the period less than one hour. The cable must operate at a load (temperature) higher than its normal load (temperature) for a short time. Therefore, the cable is designed to withstand the emergency overload rating that increases the cable temperature more than the operating cable temperature. In short circuit case occurs for half a second the 90 • C cable must withstand the highest temperature up to 250 • C.
The effect of cable temperature on the cable ampacity is investigated at different cable configurations and different mode operation (continuous and emergency operations). The operation modes refer to the maximum conductor operating temperatures. Figure 11 shows that the single bond for flat configuration represents higher ampacity rather than the other bond and cable configuration at the same cable temperature. For the same loading current, the temperature of the cable for single bond flat configuration is less than the other cable configuration.

Effect of the cable bond configuration on the cable ampacity at different cable temperature
There are two category ratings based on the load rating of the cable (cyclic ratings and emergency ratings). Cyclic rating is used the load varies daily over 24 hours and it is the same for each day. The emergency ratings are considered when the load varies between ten minutes to one hour. Therefore, the emergency refers to the transient period when the load varies in the period less than one hour. The cable must operate at a load (temperature) higher than its normal load (temperature) for a short time. Therefore, the cable is designed to withstand the emergency overload rating that increases the cable temperature more than the operating cable temperature. In short circuit case occurs for half a second the 90 o C cable must withstand the highest temperature up to 250 o C.
The effect of cable temperature on the cable ampacity is investigated at different cable configurations and different mode operation (continuous and emergency operations). The operation modes refer to the maximum conductor operating temperatures. Figure 11 shows that the single bond for flat configuration represents higher ampacity rather than the other bond and cable configuration at the same cable temperature. For the same loading current, the temperature of the cable for single bond flat configuration is less than the other cable configuration. Figure 11. Effect of the bond category on the cable ampacity with the variation of cable temperature.

Discussions
In IEC-60287 standard, the analytical method is utilized to compute the ampacity of the power cables. The analytical method is suitable when the cable is buried directly into ground and it considered the ground surface as a constant temperature boundary. In addition, it assumed constant temperature for each cable layer. The previous assumptions are ideal, and it did not meet the practical operational environment [13]. Moreover, these ideal assumptions that were introduced in various structural details and material models reflected to some extent for the obtained results [22]. Therefore, many advantages can be obtained when using the numerical methods such as boundary element method, finite difference method, and finite element method for computing the ampacity and temperature distribution of the power cables. The advantages include the less error of the computation results compared with the results from IEC 60287 standard and using the Multiphysics filed coupling analysis to enhance the computation operation process. Therefore, Figure 11. Effect of the bond category on the cable ampacity with the variation of cable temperature.

Discussions
In IEC-60287 standard, the analytical method is utilized to compute the ampacity of the power cables. The analytical method is suitable when the cable is buried directly into ground and it considered the ground surface as a constant temperature boundary. In addition, it assumed constant temperature for each cable layer. The previous assumptions are ideal, and it did not meet the practical operational environment [13]. Moreover, these ideal assumptions that were introduced in various structural details and material models reflected to some extent for the obtained results [22]. Therefore, many advantages can be obtained when using the numerical methods such as boundary element method, finite difference method, and finite element method for computing the ampacity and temperature distribution of the power cables. The advantages include the less error of the computation results compared with the results from IEC 60287 standard and using the Multiphysics filed coupling analysis to enhance the computation operation process. Therefore, in the next section the temperature distribution at specified loading current is computing using the COMSOL Multiphysics software and comparing its results with that obtained using the IEC 60287 standard.

Temperature Rise Evaluation Using COMSOL Multiphysics
The thermal performance of the cable is evaluated using COMSOL Multiphysics where the prediction of the cable temperatures depends on the generated heat losses due to the conductor current, soil thermal resistivity and the ambient temperature of the cable. Therefore, Heat transfer in solids is used as a physics interface in COMSOL Multiphysics to determine the temperature performance through the underground cable. The cable can be modeled as a 2D using COMSOL Multiphysics where it is very long. This section covers temperature rise of the cable considering different conditions such as steady state load (maximum continuous load), transient loads, and existing dry zone around the cable.
For interface of heat transfer in solid in COMSOL Multiphysics, the following equation is solved [23]: where, α, S, ρ, C p , T, u trans , q, q r , and Q are thermal expansion coefficient, second Piola-Kirchhoff stress tensor, density, specific heat capacity, temperature, velocity vector, conduc tion-based heat flux, heat flux by radiation, and additional heat sources, respectively.

Steady State Cable Temperature Response
For maximum continuous loading of the cable, 1000 A, the heat losses are computed using Equation (23) and include as a heat source of the conductor in the thermal model in COMSOL Multiphysics. Figure 12 shows the conductor and sheath temperatures for long time duration such as 1000 h where the evaluations using COMSOL and equations are declared. Equations from (5) to (23) are based on the equivalent two-loop thermal circuit of the cable. A very close results of equations and COMSOL are achieved which confirms the validation of the heat transfer model. The evaluated temperature shows the safety performance of the cable under thermal condition as its temperature is far from the maximum permissible value which is 90 • C. The temperature difference between the conductor and sheath is around 6 • C. The temperature rise percentage is 4.61% when the time changes from 100 h to 1000 h.
Processes 2021, 9, 438 12 of 16 in the next section the temperature distribution at specified loading current is computing using the COMSOL Multiphysics software and comparing its results with that obtained using the IEC 60287 standard.

Temperature Rise Evaluation Using COMSOL Multiphysics
The thermal performance of the cable is evaluated using COMSOL Multiphysics where the prediction of the cable temperatures depends on the generated heat losses due to the conductor current, soil thermal resistivity and the ambient temperature of the cable. Therefore, Heat transfer in solids is used as a physics interface in COMSOL Multiphysics to determine the temperature performance through the underground cable. The cable can be modeled as a 2D using COMSOL Multiphysics where it is very long. This section covers temperature rise of the cable considering different conditions such as steady state load (maximum continuous load), transient loads, and existing dry zone around the cable.
For interface of heat transfer in solid in COMSOL Multiphysics, the following equation is solved [23]: where, α, S, , Cp, T, utrans, q, qr, and Q are thermal expansion coefficient, second Piola-Kirchhoff stress tensor, density, specific heat capacity, temperature, velocity vector, conduction-based heat flux, heat flux by radiation, and additional heat sources, respectively.

Steady state cable temperature response
For maximum continuous loading of the cable, 1000 A, the heat losses are computed using equation (23) and include as a heat source of the conductor in the thermal model in COMSOL Multiphysics. Figure 12 shows the conductor and sheath temperatures for long time duration such as 1000 h where the evaluations using COMSOL and equations are declared. Equations from (5) to (23) are based on the equivalent two-loop thermal circuit of the cable. A very close results of equations and COMSOL are achieved which confirms the validation of the heat transfer model. The evaluated temperature shows the safety performance of the cable under thermal condition as its temperature is far from the maximum permissible value which is 90 o C. The temperature difference between the conductor and sheath is around 6 o C. The temperature rise percentage is 4.61% when the time changes from 100 hr to 1000 hr. Figure 12. The cable and sheath temperature variation through long duration. Figure 12. The cable and sheath temperature variation through long duration.

Transient Cable Temperature
Based on the daily load cycle published in [24], the corresponding load currents similar to the published performance is considered and mentioned as "Load performance 1" as shown in Figure 13. Further performances are considered such as the step changing from 1000 A to 900 A repeatedly and the maximum continuous load 1000 A as shown in Figure 13 and indicated by "Load performance 2" and "Load performance 3", respectively. The corresponding transient temperatures of the conductor and sheath based on these load performances under study are shown in Figure 14a,b, respectively. From the corresponding temperature behavior concerning Load performance 2, there are the accumulations of the conductor and sheath temperatures, where the current for the second maximum current 1000 A has a higher temperature comparing to the period during the first maximum current (temperature period (8-12 h)) although the same loading current of 1000 A. From the corresponding temperature behavior concerning Load performance 1, comparing the temperature performance through the periods (8-12 h (1000 A) and 16-20 h (960 A)) results temperature during 960 A higher than 1000 A. Therefore, changing loads have to be considered especially from lower value to highest one. as shown in Figure 13. Further performances are considered such as the step changing from 1000 A to 900 A repeatedly and the maximum continuous load 1000 A as shown in Figure 13 and indicated by "Load performance 2" and "Load performance 3", respectively. The corresponding transient temperatures of the conductor and sheath based on these load performances under study are shown in Figures 14.a and 14.b, respectively. From the corresponding temperature behavior concerning Load performance 2, there are the accumulations of the conductor and sheath temperatures, where the current for the second maximum current 1000 A has a higher temperature comparing to the period during the first maximum current (temperature period (8-12 hours)) although the same loading current of 1000 A. From the corresponding temperature behavior concerning Load performance 1, comparing the temperature performance through the periods (8hr-12hr (1000 A) and 16hr-20hr (960 A)) results temperature during 960 A higher than 1000 A. Therefore, changing loads have to be considered especially from lower value to highest one.

Impact of Dry zone formation on temperature performance
Dry zone can be formed around the underground cable due to the heating resulting from high temperature. In this study, dry zone thickness varies from 40 mm to 500 mm around the cable. The resistivity of the formation dry zone is assumed as 2.5 K.m/W. Considering maximum continuous load (1000 A), Figure 15a shows the conductor temperature rise under several thicknesses of dry zone. It shows that the temperature is raised by 15% when the dry zone thickness reaches to 500 mm. The same performance is declared from Figure 15b for sheath temperature.

Impact of Dry Zone Formation on Temperature Performance
Dry zone can be formed around the underground cable due to the heating resulting from high temperature. In this study, dry zone thickness varies from 40 mm to 500 mm around the cable. The resistivity of the formation dry zone is assumed as 2.5 K·m/W. Considering maximum continuous load (1000 A), Figure 15a shows the conductor temperature rise under several thicknesses of dry zone. It shows that the temperature is raised by 15% when the dry zone thickness reaches to 500 mm. The same performance is declared from Figure 15b for sheath temperature.

Conclusions
The cable ampacity and temperature distribution through the cable layer can be computed based on the IEC 60287 standard using an analytical method that was applied on equivalent two-loop thermal circuit. The analytical method considered some ideal assumptions such as using the ground as a constant temperature boundary and considered the temperature of the cable layer is constant. The results of the computed cable ampacity based on the equivalent two-loop thermal circuit indicated that cable configuration, the soil resistivity, the soil temperature, and the cable depth have a great effect on the cable ampacity. Due to the shortage and ideal assumption in the analytical method in IEC 60287 standard, a numerical method as finite element method was used. COMSOL Multiphysics has been used as an accurate software to evaluate the temperature performance of the cable and sheath concerning thermal model. These performances have been studied considering different loads (typical, step, and maximum load cycles). The magnitudes impact of previous loading periods on the temperature performance have been investigated as accumulated temperature attitude. Effect of dry zone formation around the cable on the temperature behavior has been evaluated, as the temperature rise increased with increasing the thickness of the dry zone taking the same pattern along the studied period.

Conclusions
The cable ampacity and temperature distribution through the cable layer can be computed based on the IEC 60287 standard using an analytical method that was applied on equivalent two-loop thermal circuit. The analytical method considered some ideal assumptions such as using the ground as a constant temperature boundary and considered the temperature of the cable layer is constant. The results of the computed cable ampacity based on the equivalent two-loop thermal circuit indicated that cable configuration, the soil resistivity, the soil temperature, and the cable depth have a great effect on the cable ampacity. Due to the shortage and ideal assumption in the analytical method in IEC 60287 standard, a numerical method as finite element method was used. COMSOL Multiphysics has been used as an accurate software to evaluate the temperature performance of the cable and sheath concerning thermal model. These performances have been studied considering different loads (typical, step, and maximum load cycles). The magnitudes impact of previous loading periods on the temperature performance have been investigated as accumulated temperature attitude. Effect of dry zone formation around the cable on the temperature behavior has been evaluated, as the temperature rise increased with increasing the thickness of the dry zone taking the same pattern along the studied period. Funding: This research was funded by TAIF UNIVERSITY RESEARCHERS SUPPORTING PROJECT, grant number TURSP-2020/146" and "The APC was funded by Mahrous Ahmed".