Fundamentals of the Thermal Analysis of Complex Arrangements of Underground Heat Sources

: This paper introduces a new analytical method for the 3-dimensional analysis of heat sources installed underground. Such sources include primarily electric power cables and steam pipes in urban areas. For complex arrangements of the heat sources, the heat transfer and cable rating calculations require the application of numerical methods, which call for high level of expertise and are generally difﬁcult to use. The computational algorithm presented in this paper uses a point sources approach and is an extension of the work published by the authors based on line source modeling, with the goal of simplifying the numerical calculations. The proposed approach is applicable for all heat sources, which are directly buried in a uniform or a partially dried out soil. The method is illustrated by several numerical examples. These examples are used for comparison with the existing standard approach.


Introduction
Knowledge of the current carrying capacity of electric power cables, or ampacity (or rating) for short, is needed when the cable system is designed and it is an important parameter during its operation. Cable ampacity is dependent on its construction, laying conditions and the duty cycle. The heat generated by the loaded cables is dispersed in the surrounding environment. Thus, the analysis involves two types of calculations: (1) evaluation of the heat generated inside the cables, and (2) heat transfer through the cable components and the surrounding soil.
The heat is generated by the load current in the conductor and, for alternating currents, by the induced currents in other metallic parts of the cable. Additionally, the alternating voltage generates dielectric losses in the insulation. The heat transfer inside and outside the cable is described by partial differential equations of heat conduction, which are in principle nonlinear, because of the temperature-dependency of the losses. For simple cable arrangements, these equations can be solved analytically with some simplifications. For such cases, the generation of the heat and the heat transfer equations are described in the national and international standards. Perhaps the most widely used standards are those published by the International Electrotechnical Commission (IEC) [1][2][3]. The derivation of the equations found in the IEC standards, as well as the long history of their development, can be found in [4].
However, in many cities around the world, cable infrastructure can be very complex, especially in the vicinity of substations, and in those situations, only numerical approaches are available for performing ampacity studies. The finite element method (FEM) or the finite difference approach (FDA) are the solutions of choice, with the former being the most widely used.
In recent years FEM has been used to analyze cable losses, with several different software packages, in 2D [5] and in 3D [6]. The method has also been used to combine electromagnetic models with thermal conduction, in 2D [7] and 3D [8]. FEM has been used to study convection and radiation heat transfer in risers (J/I-tubes) [9] and horizontal pipes (HDD) [10] with computational fluid dynamics and with heat transfer coefficients [11].
However, computer programs based on applying the finite element method to solve partial differential equations describing heat generation and heat transfer of the power cables require great skill in their application and are generally difficult to use. The first analytical description was developed in [12], in which straight sections of the cables are represented as line sources with finite lengths. Attempts to address cables crossing other heat sources for rectilinear cables were described in [13][14][15]. The concept of finite line sources in [12] was recently used to calculate the thermal fields of such complicated cable structures, taking the longitudinal heat flows along the cable conductors into account [16].
However, during some numerical applications, several important disadvantages of the method of finite line sources became evident. First, the numerical handling of the line sources, which result from subdividing straight cable sections into subsections, is elaborate and somewhat error-prone. The same holds true for the analysis of the finite line sources (the number of which may be more than 10,000 per cable), where for each of them, the coordinate system must be defined to consider each field point by means of a transformation of the coordinates. Moreover, for sloped cable sections (where the laying depth is changing), no representation by line sources is known. Furthermore, a severe disadvantage is that no formulations of the transient behavior of finite line sources are known, so that an analysis of such arrangements for time-dependent loads is not possible.
To overcome these problems, a new model is presented in this paper. This model uses point sources instead of line sources. The method can be abbreviated as PSM, referring to the point sources method. As shown below, it is easy to use and is applicable to all possible configurations of buried heat sources. The method not only enables steady-state analysis, taking into account the drying-out of the soil and/or longitudinal heat flows, but also allows the calculation of the transient thermal behavior with time-varying heat inputs.

Basic Principles
The main idea is to subdivide all heat sources in a given installation into short sections (normally of equal length of, e.g., 0.01 m), which will be represented as point sources with the emission of the local losses W.

Steady-State Temperatures
The steady-state temperature rise ∆θ caused by a point source is given by: where r-distance between a point source (with losses W) and the considered point, λ-thermal conductivity of the soil, and c-constant.
Similar to the Kennelly formula, which is normally used by cable engineers for long, straight cables and is based on an infinite line source, we consider the influence of the isothermal soil surface by reflecting the point source and weighting it with the losses -W, as illustrated in in Figure 1. Thus, we can obtain the steady-state temperature rise Δθ caused by a point source in the soil: where , _ are the distances between the source (with losses W) and the reflected source (with losses -W) in relation to the considered point. Note that summarizing an infinite number of point sources by means of (2), along an infinite line (which is also reflected at the soil surface) leads to the model shown in Figure  2 and (for constant losses) to the well-known Kennelly equation for the temperature rise caused by an infinite line source in the soil: with the distances between the considered point P (xP, yP, zP) and the point source (x1, y1, z1) expressed as: , In this paper, the losses per unit length (p.u.l.) will be denoted by W' and the losses of a point source by W. The notation z' pertains to the local coordinate system associated with the considered heat source. Thus, we can obtain the steady-state temperature rise ∆θ caused by a point source in the soil: where r + , r _ are the distances between the source (with losses W) and the reflected source (with losses -W) in relation to the considered point. Note that summarizing an infinite number of point sources by means of (2), along an infinite line (which is also reflected at the soil surface) leads to the model shown in Figure 2 and (for constant losses) to the well-known Kennelly equation for the temperature rise caused by an infinite line source in the soil: with the distances between the considered point P (x P , y P , z P ) and the point source (x 1 , y 1 , z 1 ) expressed as: From the above developments, we have the possibility of describing any section of a heat source-straight or curved-through the superposition of point sources with constant or varying losses. For this purpose, the considered heat source section must be subdivided into small sections of equal lengths. The next question is what the maximum length of these sections should be so that an approximate representation by the point sources works with a sufficient accuracy. To examine this question, consider Figure 3,  In this paper, the losses per unit length (p.u.l.) will be denoted by W and the losses of a point source by W. The notation z pertains to the local coordinate system associated with the considered heat source.
From the above developments, we have the possibility of describing any section of a heat source-straight or curved-through the superposition of point sources with constant or varying losses. For this purpose, the considered heat source section must be subdivided into small sections of equal lengths. The next question is what the maximum length of these sections should be so that an approximate representation by the point sources works with a sufficient accuracy. To examine this question, consider Figure 3, which shows the temperature rise caused by a line source extending from z = −10 m to z = +10 m and the losses (p.u.l.) of W = 100 W/m at the distance x P from the axis of the line source (at z = 0). The line source is represented by the point sources of length ∆l and losses of W ∆l The depth of laying is y = 2.0 m. From the above developments, we have the possibility of describing any section of a heat source-straight or curved-through the superposition of point sources with constant or varying losses. For this purpose, the considered heat source section must be subdivided into small sections of equal lengths. The next question is what the maximum length of these sections should be so that an approximate representation by the point sources works with a sufficient accuracy. To examine this question, consider Figure 3, which shows the temperature rise caused by a line source extending from z′ = −10 m to z′ = +10 m and the losses (p.u.l.) of W′ = 100 W/m at the distance xP from the axis of the line source (at z′ = 0). The line source is represented by the point sources of length Δl and losses of ∆ The depth of laying is y = 2.0 m. The dotted red line shows the exact solutions, resulting from (3), for different points (xP). We can observe that even for the very small lengths Δl, there is a deviation from the exact solution. This deviation is a result of the limited length of the examined line source (i.e., from the "missing ends"). This topic will be discussed later in the text.
If the value of xP is interpreted as half of the cable or heat pipe diameter (assuming that the conductor's center is on the z-axis), Figure 3a,b describe cables/pipes with diameters between 20 and 100 mm. It can be concluded that even for the very small cables, the length of Δl = 0.01 m = 1 cm will be sufficient.  The remaining deviation (as compared with the infinite line source) can be explained by examining Figure 4. In the right-hand part of this figure, cable/pipe is represented by a semi-infinite line source, which again is represented by "point sources" of length Δl. Obviously, for calculating the temperature rise ΔθP at the point at a distance xP from the heat source axis and in a distance Δz to the end of the section, the deviation of ΔθP from the exact solution for an infinite line source will increase with decreasing distance Δz, caused by the left-hand "missing end". The dotted red line shows the exact solutions, resulting from (3), for different points (x P ). We can observe that even for the very small lengths ∆l, there is a deviation from the exact solution. This deviation is a result of the limited length of the examined line source (i.e., from the "missing ends"). This topic will be discussed later in the text.
If the value of x P is interpreted as half of the cable or heat pipe diameter (assuming that the conductor's center is on the z-axis), Figure 3a,b describe cables/pipes with diameters between 20 and 100 mm. It can be concluded that even for the very small cables, the length of ∆l = 0.01 m = 1 cm will be sufficient.
The remaining deviation (as compared with the infinite line source) can be explained by examining Figure 4. In the right-hand part of this figure, cable/pipe is represented by a semi-infinite line source, which again is represented by "point sources" of length ∆l. Obviously, for calculating the temperature rise ∆θ P at the point at a distance x P from the heat source axis and in a distance ∆z to the end of the section, the deviation of ∆θ P from the exact solution for an infinite line source will increase with decreasing distance ∆z, caused by the left-hand "missing end". The remaining deviation (as compared with the infinite line source) can be explained by examining Figure 4. In the right-hand part of this figure, cable/pipe is represented by a semi-infinite line source, which again is represented by "point sources" of length Δl. Obviously, for calculating the temperature rise ΔθP at the point at a distance xP from the heat source axis and in a distance Δz to the end of the section, the deviation of ΔθP from the exact solution for an infinite line source will increase with decreasing distance Δz, caused by the left-hand "missing end". To examine this effect further, let us consider a heat source that is 130 m long. The heat source is represented by point sources (with Δl = 0.01 m) and we are interested in the temperature at the same distances xp from the z-axis, as in Figure 4. Computed temperature rises are shown in Figure 5 as a function of Δz (keep in mind that in Figure 3 we have To examine this effect further, let us consider a heat source that is 130 m long. The heat source is represented by point sources (with ∆l = 0.01 m) and we are interested in the temperature at the same distances x p from the z-axis, as in Figure 4. Computed temperature rises are shown in Figure 5 as a function of ∆z (keep in mind that in Figure 3 we have a distance of ±10 m to both ends, which will double the deviation). From Figure 5, it can be concluded that if a point of interest is in the range of 10 . . . 20 m remote from the line source end, the calculation errors caused by the missing cable ends will become negligible. a distance of ± 10 m to both ends, which will double the deviation). From Figure 5, it can be concluded that if a point of interest is in the range of 10…20 m remote from the line source end, the calculation errors caused by the missing cable ends will become negligible. Thus, we have found an algorithm which can easily be applied to all geometries of a heat source run. A precondition of the calculation is a subdivision of the run into small sections of equal length (of, e.g., 0.01 m) and to define the coordinates of the center points of these subsections, which are treated as point sources. Moreover, the emitted losses (p.u.l.) W′ are to be defined first, which can consequently be locally handled in a function of temperature. This will be discussed in detail in the following chapters. Thus, we have found an algorithm which can easily be applied to all geometries of a heat source run. A precondition of the calculation is a subdivision of the run into small sections of equal length (of, e.g., 0.01 m) and to define the coordinates of the center points of these subsections, which are treated as point sources. Moreover, the emitted losses (p.u.l.) W are to be defined first, which can consequently be locally handled in a function of temperature. This will be discussed in detail in the following chapters.

Transient Temperatures
A point source with a loss step ∆W at the time t = 0 will cause the following temperature rise in a considered point at the distance r [17]: with δ thermal diffusivity of the soil. As shown in [3], this variable can be approximated by means of the following equation: The error function erf(x) and its complementary function erfc(x) are defined by: Following [18], a numerical solution can be achieved by: with the maximum absolute value of the error of | | < 3.10 −7 . The coefficients a i are equal to: [0.0705230784, 0.0422820123, 0.0092705272, 0.0001520143, 0.0002765672, 0.0000430638] Since erfc(0) = 1, we see from (5) that for large values of time t and/or for large distances r, the solution approaches the steady-state description of (1).
Reflecting the point source at the soil surface again, as was done for the steady-state analysis, and superimposing all N point sources along a line section, gives a temperature rise of: where the point sources i at places z = i ∆z have loss steps ∆W i at times t i . A summation in (9) over an infinite number of point sources along an infinite line (which is also reflected at the soil surface) leads to the model of Figure 2 and (for constant losses) to the well-known equation for the transient temperature rise caused by an infinite line source with the loss step (p.u.l.) ∆W : with the exponential-integral-function E1: Figure 6 shows the correspondence of the solutions (10) for the line source (red lines) and (9) for the point sources (green lines). It shows the temperature rise at a point at the distance ∆y below the line caused by a line of point sources at the depth of y = 2.0 m and with the end coordinates of ±30 m. We consider a loss step at t = 0 in all point sources with ∆W = 100 W/m. It can be seen that both solutions are practically identical.
where the point sources i at places z′ = i Δz have loss steps ∆ at times ti. A summation in (9) over an infinite number of point sources along an infinite line (which is also reflected at the soil surface) leads to the model of Figure 2 and (for constant losses) to the well-known equation for the transient temperature rise caused by an infinite line source with the loss step (p.u.l.) ΔW': with the exponential-integral-function E1: Figure 6 shows the correspondence of the solutions (10) for the line source (red lines) and (9) for the point sources (green lines). It shows the temperature rise at a point at the distance Δy below the line caused by a line of point sources at the depth of y = 2.0 m and with the end coordinates of ±30 m. We consider a loss step at t = 0 in all point sources with ΔW′ = 100 W/m. It can be seen that both solutions are practically identical.

Some Basic Examples
In this chapter, some basic examples of the application of the calculation model will be presented.

Line Source with Locally Varying Losses
Consider a line source shown in Figure 7 with z′ = −30 m… + 30 m and losses (p.u.l.) of 100 W/m for the section z′ = −30 m…0 m and of 150 W/m for z′ = 0 m…30 m (the technical background may be a cable with a changing cross section or a steam pipe with a changing temperature at z′ = 0 m). These losses arise as power steps at t = 0.

Line Source with a 90° Bend
In the following, the thermal behavior of a line source with a 90° bend will be discussed, as shown in Figure 9. Such situations often arise in practice but there are no analytical tools to handle them. Figure 9 shows the geometry of the bend with the radius R. The line has a depth of y = 2.0 m and is represented along its axis (z′ = 0…100 m) by point sources, again defined by a length of Δl = 0.01 m, a depth of y = 2.0 m and a thermal conductivity of the soil of λ = 1.0 W/(K m).

Line Source with a 90 • Bend
In the following, the thermal behavior of a line source with a 90 • bend will be discussed, as shown in Figure 9. Such situations often arise in practice but there are no analytical tools to handle them. Figure 9 shows the geometry of the bend with the radius R. The line has a depth of y = 2.0 m and is represented along its axis (z = 0 . . . 100 m) by point sources, again defined by a length of ∆l = 0.01 m, a depth of y = 2.0 m and a thermal conductivity of the soil of λ = 1.0 W/(K m).

Line Source with a 90° Bend
In the following, the thermal behavior of a line source with a 90° bend will be discussed, as shown in Figure 9. Such situations often arise in practice but there are no analytical tools to handle them. Figure 9 shows the geometry of the bend with the radius R. The line has a depth of y = 2.0 m and is represented along its axis (z′ = 0…100 m) by point sources, again defined by a length of Δl = 0.01 m, a depth of y = 2.0 m and a thermal conductivity of the soil of λ = 1.0 W/(K m).  Figure 10 shows the temperature rises along the bend, 0.05 m below (y = 2.05 m) the line source, as a function of the axial coordinate z′ with the bending radius R as a parameter. With increasing z′, the bend begins at z′ = 50.0 m (see Figure 9).  Figure 10 shows the temperature rises along the bend, 0.05 m below (y = 2.05 m) the line source, as a function of the axial coordinate z with the bending radius R as a parameter. With increasing z , the bend begins at z = 50.0 m (see Figure 9).  Figure 10 illustrates that at both sides of the bend we have the maximum values of the temperature rise, which are dependent on the bending radius. A minimum value of the temperature rise appears in the middle of the bend, but this remains higher than the compared temperature rise of a straight line without a bend (69.9 K). This value is again approached by all the characteristic lines for small and large values of z′.
For the points with temperature maxima, which are marked in Figure 9 by small circles, the transient temperature rises are calculated, using (7), for a loss step (p.u.l.) of W' = 100 W/m at t = 0. The results are shown in Figure 11, in which the temperature runs are well approaching the stationary values of Figure 10.
Typical minimal bending radii Rmin of power cables are 15-25 times the cable radius. This means that for a small medium voltage cable with a diameter, e.g., of 30 mm: Rmin = 0.45…0.75 m, whereas an extra-high-voltage cable with a diameter of, e.g., 140 mm, has Rmin = 2.1…3.5 m. From Figure 10, we can conclude that in the first example the temperature rise caused by the bend is appreciable (>5%), whereas in the second example, it is nearly negligible.  Figure 10 illustrates that at both sides of the bend we have the maximum values of the temperature rise, which are dependent on the bending radius. A minimum value of the temperature rise appears in the middle of the bend, but this remains higher than the compared temperature rise of a straight line without a bend (69.9 K). This value is again approached by all the characteristic lines for small and large values of z .
For the points with temperature maxima, which are marked in Figure 9 by small circles, the transient temperature rises are calculated, using (7), for a loss step (p.u.l.) of W = 100 W/m at t = 0. The results are shown in Figure 11, in which the temperature runs are well approaching the stationary values of Figure 10.  Figure 11. Transient temperature rise in the considered points of Figure 9 with maximum temperature rise for a loss step (p.u.l.) of W′ = 100 W/m at t = 0; parameter: bending radius R.

Temperature and Rating Analysis of a Complicated Cable Route
The foregoing basic explanations of the new method for determining the temperature rise in the soil surrounding heat sources can be applied for the temperature and rating analysis of complicated cable arrangements.

Steady State, without Longitudinal Heat Fluxes
For determining the steady-state conductor temperature of a cable, the normal relations of the IEC Standard 60287 [1,2] can be used, but applied to a predefined section of cable i. For a single-core cable i, installed in a group of N cables, we obtain the conductor temperature in section by superposing the temperature rises caused by all sections 1, … in all cables 1…N: , , , , , , , , here: , , ∆ , = temperature rise of conductor section because of the dielectric losses of cable, (K) T1, T2, T3 = thermal resistance of the insulation, armor bedding and the jacket, respectively, (K.m/W). λ1, λ2 = sheath and armor loss factors, respectively = total losses of section js in cable j, (W). δi,j is the Kronecker-symbol with δi,j = 1 for i = j and δi,j = 0 otherwise. , , is the thermal resistance between section js of cable j and section is of cable, (K.m/W). Following (3), with Ns,j sectors (i.e., point sources) of the influencing cable j, we have: Figure 11. Transient temperature rise in the considered points of Figure 9 with maximum temperature rise for a loss step (p.u.l.) of W = 100 W/m at t = 0; parameter: bending radius R.
Typical minimal bending radii R min of power cables are 15-25 times the cable radius. This means that for a small medium voltage cable with a diameter, e.g., of 30 mm: R min = 0.45 . . . 0.75 m, whereas an extra-high-voltage cable with a diameter of, e.g., 140 mm, has R min = 2.1 . . . 3.5 m. From Figure 10, we can conclude that in the first example the temperature rise caused by the bend is appreciable (>5%), whereas in the second example, it is nearly negligible.

Temperature and Rating Analysis of a Complicated Cable Route
The foregoing basic explanations of the new method for determining the temperature rise in the soil surrounding heat sources can be applied for the temperature and rating analysis of complicated cable arrangements.

Steady State, without Longitudinal Heat Fluxes
For determining the steady-state conductor temperature of a cable, the normal relations of the IEC Standard 60287 [1,2] can be used, but applied to a predefined section i s of cable i. For a single-core cable i, installed in a group of N cables, we obtain the conductor temperature in section i s by superposing the temperature rises caused by all sections j s = 1, . . . N sj in all cables 1 . . . N: here: ∆θ d,i s = temperature rise of conductor section i s because of the dielectric losses of cable, (K) T 1 , T 2 , T 3 = thermal resistance of the insulation, armor bedding and the jacket, respectively, (K.m/W). λ 1 , λ 2 = sheath and armor loss factors, respectively W j s = total losses of section j s in cable j, (W). δ i,j is the Kronecker-symbol with δ i,j = 1 for i = j and δ i,j = 0 otherwise. T 4,j s, i s is the thermal resistance between section j s of cable j and section i s of cable, (K·m/W). Following (3), with N s,j sectors (i.e., point sources) of the influencing cable j, we have: where the radii r + and r --are the distances between the considered section i s and the influencing point source j s (and its image source, respectively).
We will now discuss a practical example, shown in Figure 12, which illustrates a crossing area with nine 220-kV-XLPE-cables (3 × 3 × 1 × 2500 mm 2 copper; 60 Hz, the cable construction is described in detail in [16]). The laying depth of all cables is y The cable losses in (10) are temperature-dependent. Therefore, they start with any estimated losses in all sections, e.g., for conductor temperatures of 90 °C or 20 °C, and they are iterated (3…6 times) after executing the calculation of (12)- (14), until a maximum sum of the temperature deviations in all 45,000 points of smaller than 0.5 K. For this, an adapted temperature coefficient α20* of the conductor losses is defined, which is able to take the temperature dependences of the different loss factors in the cable (ys, ys, λ1, λ2) into account. This new coefficient was defined in [19] as: In Figure 12, some special points along the axes of the cables are defined, which appear later in the temperature diagrams and which are explicitly explained in Table 1. Table 1. Explanations of the special points in Figure 12. The cable losses W j s in (10) are temperature-dependent. Therefore, they start with any estimated losses in all sections, e.g., for conductor temperatures of 90 • C or 20 • C, and they are iterated (3 . . . 6 times) after executing the calculation of (12)- (14), until a maximum sum of the temperature deviations in all 45,000 points of smaller than 0.5 K. For this, an adapted temperature coefficient α 20 * of the conductor losses is defined, which is able to take the temperature dependences of the different loss factors in the cable (y s , y s , λ 1 , λ 2 ) into account. This new coefficient was defined in [19] as: with F c,θ = 1 + k y s,θ + y p,θ ; F sa,θ = (1 + λ 1,θ + λ 2,θ )F θ = F c,θ F sa,θ The values of F are obtained from the steady state analysis for θ = 20 and θ = 90 • In our example we have: Thus, α * 20 = 0.00247/K. The dotted lines in Figure 13 show the distributions of the conductor temperatures of the middle cables 2, 5 and 8 of Figure 12 along the particular cable axes (z ). For explanations of the special points, see Figure 12 and Table 1.  Figure 12 along the particular cable axes (z′); for explanations of the special points, see Figure 12 and Table 1. Abbreviations: transv. refers to the radial heat transfer only, whereas longit represents both radial and longitudinal heat flows.

Steady State, with Consideration of the Longitudinal Heat Fluxes
Since in complex cable arrangements the conductor temperature will vary along its length because of the influence of cable crossings and changes in the surrounding environment, the longitudinal heat flow should be taken into account. This is discussed briefly below, as this subject was addressed in [16].
As before, we will assume that each cable section is (with a typical length of, e.g., 0.01 m) has constant conductor losses Wi. The temperature along the section was obtained in [16] as: Where , denotes the uninfluenced temperature, that is, without the longitudinal heat flow, and γ Let us now consider section ν with the right edge having the coordinate z′ = ν Δz. The constants A and B in (17) are obtained iteratively as [16]: The calculation starts with z′ = 0 at the left border of Subsection 1 in Section 1, that is, at a point where a negligible longitudinal heat flow can be assumed (typically at a distance of 10…20 m from the critical area). This means that: Figure 13. Distributions of the conductor temperatures of the middle cables 2, 5 and 8 in Figure 12 along the particular cable axes (z ); for explanations of the special points, see Figure 12 and Table 1. Abbreviations: transv. refers to the radial heat transfer only, whereas longit represents both radial and longitudinal heat flows.

Steady State, with Consideration of the Longitudinal Heat Fluxes
Since in complex cable arrangements the conductor temperature will vary along its length because of the influence of cable crossings and changes in the surrounding environment, the longitudinal heat flow should be taken into account. This is discussed briefly below, as this subject was addressed in [16].
As before, we will assume that each cable section i s (with a typical length of, e.g., 0.01 m) has constant conductor losses W i . The temperature along the section i s was obtained in [16] as: where θ u,i s denotes the uninfluenced temperature, that is, without the longitudinal heat flow, and Let us now consider section ν with the right edge having the coordinate z = ν ∆z. The constants A and B in (17) are obtained iteratively as [16]: The calculation starts with z = 0 at the left border of Subsection 1 in Section 1, that is, at a point where a negligible longitudinal heat flow can be assumed (typically at a distance of 10 . . . 20 m from the critical area). This means that: and: The computational algorithm starts with the assumed temperatures for all cable sections. Next, through an iterative process, new temperatures are obtained and the corresponding losses are computed. Usually, four to five iterations are required to achieve convergence with 0.1 K accuracy.
The full lines in Figure 13 show the resulting temperature distributions, compared with the foregoing distributions (thin, dotted). Comparing cable 2 with the other cables, 5 and 8, leads to the conclusion that the broad saddle of high temperatures in cable 2 is hindering longitudinal heat fluxes to the colder ends, so that the resulting temperatures are only slightly reduced.
In order to test the PSM model, finite element (FEM) studies have been performed using the same conditions as in Figure 12. The mesh generated by the COMSOL Multiphysics program is shown in Figure 14 with the comparative results shown in Figure 15.
We can observe that the results are very close, within 3 K difference. Keeping in mind that through FEM an additional longitudinal heat flow along the metal screens of the cables will be taken into account-which, e.g., in [15] caused a rating difference of nearly 1%, which again may cause a temperature difference of 1 . . . 2 K-these temperature differences are further relativized. and: The computational algorithm starts with the assumed temperatures for all cable sections. Next, through an iterative process, new temperatures are obtained and the corresponding losses are computed. Usually, four to five iterations are required to achieve convergence with 0.1 K accuracy.
The full lines in Figure 13 show the resulting temperature distributions, compared with the foregoing distributions (thin, dotted). Comparing cable 2 with the other cables, 5 and 8, leads to the conclusion that the broad saddle of high temperatures in cable 2 is hindering longitudinal heat fluxes to the colder ends, so that the resulting temperatures are only slightly reduced.
In order to test the PSM model, finite element (FEM) studies have been performed using the same conditions as in Figure 12. The mesh generated by the COMSOL Multiphysics program is shown in Figure 14 with the comparative results shown in Figure 15.   with its own thermal capacitance. In order to account for the thermal capacitances of the internal cable components, the total temperature rise of the cable is obtained as [3]: where the attainment factor, defined as, reflects the fact that the step response is initially (in a range of some hours) dominated by the inner temperature rise of the cable. We can observe that the results are very close, within 3 K difference. Keeping in mind that through FEM an additional longitudinal heat flow along the metal screens of the cables will be taken into account-which, e.g., in [15] caused a rating difference of nearly 1%, which again may cause a temperature difference of 1…2 K-these temperature differences are further relativized.

Transient Behaviour (without Longitudinal Heat Fluxes)
Applying (9) offers the possibility of analyzing the transient behavior of complex cable systems. The theory presented in Section 2.2 dealt with heat generated by a line source. A typical power cable is composed of several conducting and nonconducting layers, each with its own thermal capacitance. In order to account for the thermal capacitances of the internal cable components, the total temperature rise of the cable is obtained as [3]: where the attainment factor, defined as, reflects the fact that the step response is initially (in a range of some hours) dominated by the inner temperature rise of the cable. For the arrangement of Figure 12, this is done under the assumption that the cable system 1 (the cables 1…3) are steadily loaded with a current of 1600 A, whereas the other six cables are switched on with 1600 A each at the time t = 0. The losses during the considered times at t = 100 h, 1000 h, 10,000 h and 100,000 h (t  ∞) are taken from the final situation, which is shown in Figures 13 and 15. For the arrangement of Figure 12, this is done under the assumption that the cable system 1 (the cables 1 . . . 3) are steadily loaded with a current of 1600 A, whereas the other six cables are switched on with 1600 A each at the time t = 0. The losses during the considered times at t = 100 h, 1000 h, 10,000 h and 100,000 h (t ⊂ ∞) are taken from the final situation, which is shown in Figures 13 and 15.

Transient Behaviour with Consideration of the Longitudinal Heat Fluxes
As an example, the installation of the nine 220-kV cables in Figure 12 is chosen and the transient case of Section 4.3. Figure 16 shows the temperature distributions for the times named above with a steady load of 1600 A in system I, whereas cables 4 . . . 9 in systems II (here: cable 5) and III (here: cable 8) are switched on at t = 0 with a load step of 1600 A. It can be seen that for t = 10,000 h the final distributions are close to the steady-state ones and for t = 100,000 h (t ⊂ ∞) are almost exactly equal to the steady-state values. 1600 A. It can be seen that for t = 10,000 h the final distributions are close to the steadystate ones and for t = 100,000 h (t  ∞) are almost exactly equal to the steady-state values. Figure 16. Conductor temperatures of the middle cables 2, 5 and 8 Figure 12: in system I (cable 2) with a steady load of 1600 A, whereas in systems II (cable 5) and III (cable 8) with load steps of 1600 A at t = 0 in cables 4…9; temperature distributions for (a) t = 100 h, (b) 1000 h, (c) 10,000 h, (d) 100,000 h; for explanations of the special points: see Figure 12 and Table 1.

Conclusions
A new method for temperature and rating calculations in arbitrarily complicated cable arrangements was introduced in this paper. The PSM method results in much simpler equations than the traditional approaches. For example, for the steady-state solution, 1/r replaces the logarithm of r and in the transient case, (1/r) erfc(r, t) replaces the exponential integral EI (r, t).
The most important conclusions can be summarized as follows:  The most simple geometric definition is achieved by using equidistant points along the cable axis.  This method is much easier to handle compared with the finite line sources with their special systems of coordinates. with a steady load of 1600 A, whereas in systems II (cable 5) and III (cable 8) with load steps of 1600 A at t = 0 in cables 4 . . . 9; temperature distributions for (a) t = 100 h, (b) 1000 h, (c) 10,000 h, (d) 100,000 h; for explanations of the special points: see Figure 12 and Table 1.

Conclusions
A new method for temperature and rating calculations in arbitrarily complicated cable arrangements was introduced in this paper. The PSM method results in much simpler equations than the traditional approaches. For example, for the steady-state solution, 1/r replaces the logarithm of r and in the transient case, (1/r) erfc(r, t) replaces the exponential integral EI (r, t).
The most important conclusions can be summarized as follows: • The most simple geometric definition is achieved by using equidistant points along the cable axis. • This method is much easier to handle compared with the finite line sources with their special systems of coordinates. • It provides a very simple formulae for expressing temperature rises, relating only to the distance r between the considered point and the point source.  Funding: This research received no external funding.

Institutional Review Board Statement:
There were no studies on human or animals in this work.