Next Article in Journal
A Spectral Method Algorithm for Modeling the Dispersion of Non-Axisymmetric Modes in Fluid-Filled Elastic Tubes
Previous Article in Journal
Multi-View Masked Autoencoder for General Image Representation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

One Hundred and Fifty Years of Skin Effect

Department of Electrical Power Engineering, Brno University of Technology, Technicka 12, 61600 Brno, Czech Republic
Appl. Sci. 2023, 13(22), 12416; https://doi.org/10.3390/app132212416
Submission received: 9 October 2023 / Revised: 8 November 2023 / Accepted: 10 November 2023 / Published: 16 November 2023

Abstract

:
In 1873, J.C. Maxwell derived relations for current density and current in a solitary long cylindrical conductor supplied with variable current. According to Maxwell, the current density in a conductor increases towards the conductor surface. This phenomenon is called the skin effect. The skin effect affects in particular the inductance of a line formed by one or several long parallel conductors. A number of papers have been published regarding the skin effect and its effect on inductance. Another phenomenon, closely related to the skin effect, is the proximity effect, which refers to the dependence of the current density in a conductor on the proximity of other conductors through which a time-variable current is flowing. Many published papers deal with the calculation of current density in two conductors, using the method for calculating current density in a solitary conductor. All the phenomena given above can be analysed and quantitatively described based on the knowledge of current density in conductors, and therefore the method for calculating current density in a group of conductors is of fundamental significance. It follows from the analysis performed that the skin effect is not a general characteristic of current density in long conductors, except for in the solitary conductor. This conclusion affects the knowledge of the phenomena associated with the skin effect.

1. Introduction

Current density in a conductor through which variable current passes increases towards the conductor surface. This phenomenon is called the skin effect. The skin effect is inextricably linked with current density. It is therefore strange that the existence of this phenomenon is generally assumed, even though an exact method for the calculation of current density was until recently known only for the infinitely long solitary conductor of circular cross section with constant resistivity. The aim of the present paper is an analysis of current density in long parallel conductors of arbitrary cross section with respect to the skin effect and other related phenomena. The starting point of the analysis has been the recently proposed accurate method for the calculation of current density in the conductors under investigation.
In 1873, J.C. Maxwell published the book “A Treatise on Electricity and Magnetism” [1]. The book is in two volumes and divided into articles; it was published in three editions [1,2,3]. The next editions, for example [4], are unabridged republications of the third edition [3]. Articles 689 and 690 in [1,2,3] deal with the calculation of current density and current in a solitary straight infinitely long conductor of circular cross section supplied with variable current. According to Maxwell, the amplitude of the current density of alternating current (of non-zero frequency) passing through a conductor increases towards the surface of the conductor. This phenomenon is called the skin effect.
The term skin effect is not used in [1,2,3]; it was probably introduced later. For example, according to [5], Horace Lamb published a paper in 1883 [6] which dealt with the motions of electricity produced in a spherical conductor by any electric or magnetic operations outside it, which was later to be known as the skin effect [5]. J.J. Thomson, who prepared the third edition of Maxwell’s book, added in the form of a footnote a derivation of the formula for calculating the inductance of a cylindrical conductor supplied with sinusoidal current. As part of this derivation, he also obtained a relation for the dependence of conductor resistance on the frequency of the supply current.
Another phenomenon closely related to the skin effect is the proximity effect, which refers to the dependence of current density in the conductor on the proximity of further conductors through which a time-varying current flows. Many published papers deal with the approximate calculation of current density in two conductors using the method for calculating current density in a solitary conductor. What is addressed in papers on the calculation of current density in two conductors is usually the skin effect, the proximity effect and the dependence of resistance on frequency, and their effect on self-inductance. Since the publication of the book [1], a plethora of articles and parts of books have been published that deal with the subject matter of the present review. In [7], published in 1959, the list of references contains 454 items, and one recent item is [8]. Thus, it is not possible to evaluate all the relevant papers.
All the phenomena given above can be analysed and quantitatively described based on the knowledge of current density in conductors, and therefore the method for calculating current density in a group of conductors is of fundamental significance. An accurate method for calculating current density in conductors forming a group has only recently been published [9,10,11,12]. The subject of investigation in [9,10,11,12] is a physical and mathematical model of a task referred to as the basic task.
The basic task is the calculation of the current density J bas ( x , y , z ) , in a group of conductors that do not move and have the following properties;
  • Each of them has a long part parallel to the axis z in the coordinate system x y z ;
  • their cross section is arbitrary but does not change along the long part;
  • some are connected at the end to a voltage or current source;
  • some are connected at the end to other conductors in the group;
  • the permeability of the conductors and their surroundings equals the vacuum permeability, μ 0 ;
  • they are placed in a magnetic field produced by a number of conductors that belong to the group or by an external source.
Solving the basic task is trivial if the magnetic field in which the conductors are placed does not change with time. In the case of a magnetic field varying with time, the solving of the basic task is problematic if some parts of the conductors are not parallel to z. In the physical model of the basic task it is assumed that
  • the conductors in the group have only a long part, which is infinitely long, and each conductor is therefore determined by its cross section;
  • the conductor resistivity, ϱ , is given and only depends on x and y;
  • the phenomena are quasi-stationary, i.e., the propagation speed of current and magnetic field is infinitely high;
  • the displacement current and leakage current between conductors are neglected.
The physical model is based, in the first place, on Faraday’s law of electromagnetic induction and is an application of the Biot and Savart law, loop current method, Kirchhoff’s voltage law, and Ohm’s law. Applying these laws to a finite segment of conductors yields a system of equations whose solution is a current density whose nonzero component is only the z-component:
J ( x , y ) = 0 , 0 , J ( x , y ) .
In the mathematical model, the conductor cross sections are approximated by the a ( N , d ) aggregate of disjunct rectangles A i ,
a ( N , d ) = i = 1 N A i ,
in which only the resistivity, ϱ i , and the current density, J i , have constant values, and
d = max i = 1 , 2 , , N d i ,
where d i is the length of the diagonal of rectangle A i . The mathematical model is a system of N ordinary differential equations of the first order [9,10,11,12]. A detailed description of the physical and mathematical models is given in [9,10,11,12], together with the solution to specific cases. An exact solution to the mathematical model is the set
J mat = { J 1 , J 2 , , J N } .
J mat can also be considered a piecewise constant function, J mat = J mat ( x , y , d ) , defined in the plane x y :
J mat ( x , y , d ) = 0 f o r ( x , y ) a ( N , d ) J i f o r ( x , y ) A i .
It holds that
lim d 0 J mat ( x , y , d ) = J ( x , y ) .
This means that J mat can approximate the solution J ( x , y ) of the physical model, theoretically with an arbitrary accuracy. In practice, the accuracy is limited by the level of the calculation means used. For a sufficiently small d, J mat ( x , y , d ) is an accurate solution to the physical model and can be considered a solution to the basic task. The method for the calculation of the current density is numerical. In the calculation, the number of rectangles that approximate the cross sections of conductors is approximately 10,000. The accuracy of the calculation is assessed via comparing the results for different values of d.
Part of the calculation of current density in conductors is the calculation of independent magnetic fluxes. The determination of independent fluxes is described in [9,10,11,12]. If in the calculation of current density in a solitary conductor or in a loop formed by two conductors the conductor cross sections are symmetrical and the values of conductor resistivity are also symmetrical, the determination of independent fluxes is simple if symmetry is used [13,14,15,16,17,18,19,20,21].

2. Solitary Conductor Concept

In ref. [1], Art. 689 and 690, Maxwell proposed a method for the calculation of current density, J sol , in a solitary long cylindrical conductor. In the following analysis of Maxwell’s result, the quantities are denoted by currently used symbols. In ref. [1], J sol is the solution to a differential equation; it depends on the distance, r, from the conductor axis and on t,
J sol ( r , t ) = μ 0 4 π ρ k = 0 μ 0 r 2 4 ρ k i ( k + 1 ) ( t ) ( k ! ) 2 ,
where
i ( t ) = i ( 0 ) ( t ) , i ( k + 1 ) ( t ) = d k + 1 d t k + 1 i ( t ) .
i ( t ) is the given function and, according to (19) in [22], it holds that
J sol ( 0 , t ) = μ 0 4 π ϱ i ( 1 ) ( t ) .
If the function i ( t ) is sinusoidal, then the function J sol is also sinusoidal,
J sol ( r , t ) = J ^ sol ( r ) sin [ ω t + ε ( r ) ] , ω = 2 π f ,
and its phasor, J ̲ sol , is a solution to the equation
d 2 J ̲ sol d r 2 + 1 r d J ̲ sol d r j κ 2 J ̲ sol = 0 , κ 2 = ω μ 0 ϱ ,
which has been derived, for example, in ref. [16].
Using Maxwell’s equations in differential form, an equation can be derived for the solitary conductor with an arbitrary cross section, which is satisfied by current density over the cross section of a long conductor,
2 J μ 0 γ J t = 0 ,
where γ is the electric conductivity of the conductor. Expressing current density in (7) via the vector potential, A , will yield an equation that is similar to (7).
The solution to Equation (6) can be expressed using the Bessel or the Kelvin functions [23], Art. 16.4. Since the publication of [1], a number of methods have been published based on the solitary conductor concept; some of the recent ones are, for example, [24,25,26,27,28].
Example 1.
Cylindrical Al conductor ( ϱ = 2.650 × 10 8 Ω · m at 20 C [29]), whose axis is the axis z and the cross section radius r 1 = 6 mm. The function i ( t ) has been chosen such that it holds that J sol ( 0 , t ) = i ^ sin ω t , where i ^ = 1 A·m 2 .
The solution to Equation (6) in Example 1 is given in Figure 1.
The curve ε ( r ) for f = 10 4 Hz in Figure 1 has two parts and is discontinuous. It is only a seeming discontinuity since ε is modified such that ε ( π , π ] , which is possible because the period of the function sin is 2 π ; ε ( r ) is continuous for f = 10 3 Hz. The number of seeming discontinuities of ε ( r ) increases with increasing f. The function J ^ sol ( r ) is increasing (monotonous) in contrast to the function J sol ( r ) , which is not monotonous, and for the discontinuous function ε ( r ) it is oscillating. This property of current density is not limited to the solitary conductor of circular cross section; it also occurs, for example, with the loop formed by two conductors of rectangular cross section. If the surface of the conductors is parallel with the planes x y and x z , then the current density, J ( x , y ) , may oscillate for y = c o n s t , as can be seen in Figure 3 in ref. [21].
The amplitude, J ^ sol , is the largest on the conductor surface, r = r 1 , and it decreases with deceasing r; for r = r 1 δ , its value for high frequencies is e-times smaller than on the conductor surface, where δ is the so-called penetration depth [30]:
δ = 2 ϱ ω μ 0 .
Related to sinusoidal current density in a solitary conductor is the technical term skin effect, which should characterize current density in conductors. According to ref. [31], item 121-13-18, “Skin effect is a phenomenon in which the current density is greater near the surface than in the interior of the conductor for an alternating electric current in the conductor”. In the literature are further definitions, which are either similar to the definition in ref. [31] or they define the skin effect as a phenomenon in which most of the current in the conductor flows through a thin layer near the conductor surface. In all the definitions, only the word ‘conductor’ is given, although all the definitions were proposed in analyses of the current density in a ‘solitary cylindrical conductor’. Examples are given below that raise doubts as to the given skin effect definitions.
Figure 2 gives the current density J sol = J sol ( r , t ) in Example 1 at five instants, t = θ T , where θ = 0 , 0.3 , 0.4 , 0.5 , and T = 1 / f , f = 10 4 Hz.
It follows from Figure 2 that for a solitary cylindrical conductor the conductor surface exhibits the maximum J ^ sol but not always the maximum current density. For a solitary conductor of rectangular shape, the situation is similar.
Example 2.
Solitary Al conductor ( ϱ = 2.650 × 10 8 Ω · m at 20 C [29]) with the rectangle cross section [ 0 , 10 ] mm × [ 0 , 16 ] mm. The current through the conductor is
I rec ( t ) = I ^ rec sin ω t , where I ^ rec = 1 A , f = 10 4 Hz .
This example is given, inter alia, because the current density,
J rec ( x , y , t ) = J ^ rec ( x , y ) sin [ ω t + ε rec ( x , y ) ] ,
in a solitary rectangular conductor has not been calculated using the method proposed by Maxwell nor by methods derived from his method. The reason for this is ignorance of the initial and boundary conditions necessary for the solution of the respective differential equation. Figure 3 provides the dependence of the amplitude, J ^ rec , on x for several values of y in Example 2. The values of J ^ rec were obtained by a method published in refs. [9,10,11,12], and in Figure 3 they are only given in one quarter of the conductor cross section because not only the amplitude but also the initial current density phase are symmetrical along the symmetry axes of the cross section x = 5 mm and y = 8 mm.
Figure 4 gives the current density in a quarter of the rectangular cross section of a solitary conductor at the instant t = 0.39 T in Example 2.
For a solitary tubular conductor with internal non-zero radius and constant resistivity, the skin effect obviously does not occur. J ^ sol is maximal and minimal near the outer and the inner surface, respectively; see Figures 2, 4, and 6 in ref. [16] and Figure 3 in ref. [17]. The skin effect is associated in practice not only with solitary conductors but also with two conductors, although in its definition (item 121-13-18 in ref. [31]) there is only ‘conductor’. It does not occur in the loop formed by two coaxial tubular conductors; see Figures 3, 5, and 7 in ref. [13] and Figures 4, 5, and 6 in ref. [17]. In some cases, it does not occur even in the loop formed by two parallel conductors; see Figures 10 and 11 in ref. [18].

3. Solitary Conductor versus Reality

Significant for practice are the theoretical models that are the most faithful image of reality. An infinitely long cylindrical conductor through which steady sinusoidal current flows is the basis of the solitary conductor concept. A conductor connected to the ideal current source can be considered solitary [32]. The ideal current source cannot be realized because, on a finite segment of the conductor, it produces an infinitely large voltage, but in spite of this a finite voltage is assumed in ref. [1], Art. 690. J.J. Thomson, who prepared the third edition [3] of [1], added in the form of a footnote a derivation of the formula for the calculation of inductance of a solitary conductor and the dependence of effective resistance on the frequency of the sinusoidal source the examined conductor is connected to. The skin effect is defined for the solitary conductor. Questioning the existence of solitary conductor also means questioning the skin effect.
Closest to the cylindrical conductor is in fact a pair of coaxial conductors connected to a source. The relation between current density in a solitary cylindrical conductor and in a pair of coaxial conductors is demonstrated in the following example.
Example 3.
A pair of coaxial conductors at a temperature of 25 C, whose axis is the axis z, is formed by an inner cylindrical Cu conductor: cross section radius r 1 = 10 mm, resistivity ϱ 1 = 1.712 × 10 8 Ω · m [29], and by an outer tubular Al conductor: inner radius r 2 = 11 mm, resistivity ϱ 1 = 2.709 × 10 8 Ω · m [29], wall thickness Δ r . A source of steady sinusoidal voltage with f = 10 3 Hz causes the drop V ^ sin ω t , where V ^ = 1 V·m 1 .
The steady current density in coaxial conductors,
J coa ( r , t ) = J ^ coa ( r ) sin [ ω t + ε coa ( r ) ] ,
is symmetrical with respect to the axis z; for the given parameters r 1 , r 2 , Δ r , ϱ 1 , ϱ 2 , f , it only depends on the distance, r, from the axis z and on t, and is directly proportional to the voltage drop on a segment of both conductors. According to refs. [13,18,22], it holds for a given voltage drop:
lim r 2 J ^ coa ( r ) = 0 , for all r .
It should be noted that (9) holds for the quasi-stationary state, which is assumed in the present review ( see Section 1), but assuming the quasi-stationary state is not adequate for a large r 2 . According to ref. [33], the assumption of the quasi-stationary state is only adequate when the maximum dimension max of the system in the plane x y is much smaller than the free space wavelength, i.e.,
max c / f , where c is the speed of light .
If the maximum magnitudes of r 1 and Δ r are several centimetres, then the inequality, (10) corresponds to the inequality r 2 c / f .
In Figure 5, the function J ^ coa ( r ) in Example 3 is illustrated for three values of Δ r . The method for calculating current density in two coaxial conductors was described in refs. [9,12,13,14,15]. In Example 3, r 2 = 0.02 c / f was chosen, i.e., r 2 = 6000 m, in agreement with (10). It was established by calculation that for Δ r = 5 mm and for all r, J ^ coa is smaller than 6.73 × 10 5 A·m 2 and 1.13 × 10 3 A·m 2 in the inner and outer conductor, respectively. The maximum value of J ^ coa for r 2 = 6000 m is only 3% of the maximum value of J ^ coa for r 2 = 11 mm.
From the above, it follows that a solitary conductor connected to the ideal current source cannot be implemented. The inner conductor in a pair of coaxial conductors cannot be considered a solitary conductor either if the outer conductor is far from the inner conductor. The author of ref. [34] is of the opposite opinion when he says:“… the solitary conductor is a realizable abstraction. A solitary conductor is one where the radially varying inner electric and magnetic fields do not suffer any azimuthal variations, since no other neighbouring conductor exists. However, in a centred coaxial cable of long length, the shield conductor has no influence on the inner cylindrical conductor…”. The author of ref. [34] does not say what effect on the inner conductor is concerned. If a magnetic field is concerned that is produced by coaxial conductors supplied by a dc source, then the author of ref. [34] is absolutely right. However, the subject of this review is the effect of the magnetic field produced by a source varying with time. Current density in conductors is a superposition of a source-driven current density and the density of eddy currents [12,14,15]. The eddy currents are determined by magnetic fluxes produced by the two conductors. The effect of the outer conductor on J ^ coa in the inner conductor in Example 3 can be seen in Figure 5. It is difficult to imagine calculating current density in conductors that are not of circular cross section using the current density in a cylindrical conductor.
Current density in a conductor is connected with the vague term proximity effect. According to ref. [31], item 121-13-19, “The proximity effect is a non-uniform distribution of electric current density in a conductor, attributable to electric currents in neighbouring conductors”. It is necessary to add that the cause of the non-uniform density of current in the conductor can also be the resistivity, which is not constant over the conductor cross section. If the cause of the non-uniform distribution of current density in the conductor is the current induced according to Faraday’s law, then the induced current can also be produced by the conductor itself or by distant conductors, and also by, for example, a moving permanent magnet. The neighbouring conductor may be solitary, unconnected to the source, or it may form with the respective conductor a circuit in which there is a source, or it may form a part of a circuit with a source, or the neighbouring conductor is connected to the same pole of the source as the respective conductor. According to [9,10,11,12], the neighbouring conductor can be included in the calculation of current density in the respective conductor, and thus it is redundant to deal with the proximity effect.

4. Self-Inductance of a Pair of Conductors

Self-induction of a pair of conductors is characterized by the quantity self-inductance, L. L determines the voltage induced in the conductors and is part of conductor impedance. There is no universal method for the calculation of L, because its value, among other things, depends on the current density in the conductors. The simple formulae given in the literature have been derived on the assumption of constant current density. This is a rather simplifying assumption because, in such a case, the induction phenomenon does not occur.
In the calculation of L, it is necessary to respect Formula (10), according to which the maximum distance of the conductors and the frequency of the current flowing through the conductors are limited. With a conductor distance of several metres, the frequency should be less than 10 7 Hz. In ref. [13], the values are given for the energy in the magnetic and electric fields of coaxial conductors. From a comparison of these values, it results that neglecting the displacement current in the insulation between the conductors is acceptable for supply voltage frequencies not exceeding 10 7 Hz if there is vacuum between the conductors. For a different insulation, the upper frequency limit would be lower. It is therefore necessary to assume the fulfilment of the condition f < 10 7 Hz; similar to ref. [35], p. 132, frequencies of up to 1 MHz are assumed so that quasi-stationary approximation may be assumed.
Formulae for the calculation of the self-inductance of a pair of long conductors are given in the literature, sometimes also with a derivation, e.g., [33,36,37,38,39,40,41,42,43,44,45,46,47]. Papers regarding the calculation of current density in a pair of long conductors often contain the method for calculating the self-inductance of a segment of the two conductors or the equivalent impedance of this segment [1,2,4,14,15,18,20,22,24,34].
Figure 6 gives the equivalent schematic of the segment of conductors between the planes z = z 1 and z = z 2 , where z 2 > z 1 .
The inductor is characterized by the self-inductance coefficient L, which is used to calculate the voltage induced on the inductor,
U L ( t ) = L d I ( t ) d t ,
where I ( t ) is the current in the loop segment or in the circuit in Figure 6. The result of solving the basic task is the vector of current density in the loop conductors, and I ( t ) is the flux of this vector through the conductor cross section. Solving the basic task thus also yields I ( t ) . If the voltage U ( t ) = V ( z 1 , t ) V ( z 2 , t ) is given, then Kirchhoff’s voltage law holds for the loop and for the circuit in Figure 6,
U L ( t ) = U ( t ) U R ( t ) , where U R ( t ) = R I ( t ) ,
with R being the given resistance of the loop segment. This means that determining L is redundant, and in the equivalent diagram (Figure 6) the inductor characterized by L is replaced by the voltage source U L ( t ) . Unlike the calculation using L, the calculation of U L ( t ) using (12) is accurate, in particular in cases when the proximity effect manifests itself or when U L ( t ) is not sinusoidal.
In the literature, the calculation of L is, in most cases, part of the calculation of the impedance of the circuit in Figure 6,
Z ̲ = ( Z ̲ ) + j ω L ,
where ω is the angular frequency of all the voltages and currents in the circuit. The application of (13) is conditional on the steady state of the circuit and on the distance of the conductors at which the proximity effect can be neglected. Another possible calculation of L consists of replacing each of the conductors with a solitary conductor and replacing L by the sum of internal and external inductance [43,48,49,50,51,52]. In the calculation of internal inductance, what is taken into account is only the magnetic flux inside the conductor and produced by this conductor. In the calculation of external inductance, only the magnetic flux between conductors is taken into consideration. This calculation procedure is mainly applied when the conductors are cylindrical because current density in solitary cylindrical conductors is easy to calculate.
The real component ( Z ̲ ) increases with increasing ω = 2 π f , and this is the reason for the incorrect statement “resistance increases with frequency”, which is frequently given in the literature and applied in technical practice. Unfortunately, even in the definition of skin effect, in item 121-13-18 [31], there is a note 1: “The skin effect increases the resistance and decreases the inductance of a conductor with the frequency of the electric current”. The standard [53] defines resistance (1) and resistance (2). Resistance (1) is determined by Ohm’s law, which holds for metallic conductors over a large range of parameters. The dependence of resistance (1) on f would above all mean the dependence of resistivity on f. According to ref. [40], Ohm’s law holds for metallic conductors if f 10 14 Hz. The value of resistance (2) is equal to ( Z ̲ ) according to [53], and therefore for f > 0 it is necessary to distinguish between R and ( Z ̲ ) .

5. Conclusions

The subject of the present paper is current density and current in a group of long conductors. The group is described in Introduction. The publication of Maxwell’s book [1] marked the beginning of a one hundred and fifty year history of attempts at proposing a method for the calculation of current density in long parallel solid conductors.
Maxwell’s method can only be applied to a solitary conductor connected to an ideal current source, which cannot be implemented. If it is assumed that the current density in the conductors of a loop formed by two conductors is the current density of two solitary conductors, then the accuracy of such an approximation increases with increasing distance of the conductors. In that case, the proximity effect cannot be analysed, and the effect of current density on the result of the calculation of inductance is small because the decisive influence is that of the magnetic flux through the area between the conductors.
All the papers published up to now, with the exception of those by the author of this review and his fellow workers, deal with the solitary conductor, or they try to transfer the result of the solitary conductor to a group of conductors or they solve Equation (6) or (7). The current density calculated by these methods is approximate, and its application in the analysis of the proximity effect and in the calculation of inductance is therefore also approximate. The knowledge of exact current density obtained via solving the basic task makes it possible to quantify the proximity effect and thus not only to state that the current density in the conductor is affected by current in the neighbouring conductors. The proximity effect affects the voltage on the inductor, and the calculation of this voltage in terms of self-inductance is highly inaccurate and therefore it is suitable to characterize the inductor by a time-dependent induced voltage instead of self-inductance.
The definition in ref. [31] is special in that the skin effect is not defined for conductors but only for the conductor. On the other hand, it is very general because it does not contain any assumption as to the conductor length and shape. It is assumed in the definition that the current in the conductor is alternating, but in ref. [31], item 121-13-18, there is a note 2: “The skin effect occurs also in the more general case of any time-varying current”. It follows from the specific cases quoted above that the skin effect is not an entirely general phenomenon, the more so because no analysis has so far been made of current density in conductors during transient phenomena and with conductor resistivity that is not constant over the conductor cross section. When using the term skin effect, it would therefore be advisable to state the relevant conditions: the form of the dependence of the voltage or current source on time, the shape of the conductor cross section, the position of the conductors in space, and in the case of transient phenomena, the initial current density. Playing with words in connection with the skin effect and the proximity effect can be replaced by an indisputable statement: “Current density is not constant in conductors that are in a time-varying magnetic field”.
The skin effect is connected with the so-called resistance dependence on the frequency of the passing current. A strict distinction needs to be made between the resistance according to Ohm’s law and the real component of impedance. For alternating current with non-zero frequency, the resistance according to Ohm’s law does not depend on frequency and its value is not equal to the value of the real component of impedance, which depends on frequency.
Knowing the current density obtained by the solution of the basic task enables analysis of the phenomena connected with the skin effect, in particular analysing an accurate calculation of voltage on the conductor. This does not mean a refutation of the papers published so far. The methods for inductance calculation in these papers need to be taken as approximate but sufficiently accurate, in particular for larger conductor distances. Nevertheless, the accuracy of the results obtained by these methods can be obtained only via a comparison with the results obtained by the solution of the basic task.

Funding

This research work has been carried out in the Centre for Research and Utilization of Renewable Energy (CVVOZE). The author gratefully acknowledges financial support from the Ministry of Education, Youth and Sports of the Czech Republic under Brno University of Technology specific research programme (project No. FEKT-S-23-8403).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Maxwell, J.C. A Treatise on Electricity and Magnetism; Clarendon Press: Oxford, UK, 1873. [Google Scholar]
  2. Maxwell, J.C. A Treatise on Electricity and Magnetism; Clarendon Press: Oxford, UK, 1881. [Google Scholar]
  3. Maxwell, J.C. A Treatise on Electricity and Magnetism; Clarendon Press: Oxford, UK, 1891. [Google Scholar]
  4. Maxwell, J.C. A Treatise on Electricity and Magnetism; Dover Publication, Inc.: New York, NY, USA, 1954. [Google Scholar]
  5. Horace Lamb. Available online: https://en.wikipedia.org/wiki/Horace_Lamb (accessed on 7 June 2023).
  6. Lamb, H. On electrical motions in a spherical conductor. Philos. Trans. R. Soc. Lond. 1883, 174, 519–549. [Google Scholar] [CrossRef]
  7. Mukoseev, Y.L. Raspredelenie Peremennogo Toka v Tokoprovodakh (Distribution of Time-Variable Current in Conductors); Gos. Energet. Izdat: Moscow, Russia, 1959. [Google Scholar]
  8. Brandão Faria, J.A.; Araneo, R.; Stracqualursi, E. Surface Wave, Skin Effect, and Per Unit Length Parameters of the Single-Wire Transmission Line at Low Frequency, for Nonmagnetic and Magnetic Wires. IEEE Access 2023, 11, 59621–59635. [Google Scholar] [CrossRef]
  9. Coufal, O. Current density in a group of long parallel conductors. Phys. Scr. 2019, 94, 125504. [Google Scholar] [CrossRef]
  10. Coufal, O. Transient Current Density in a Pair of Long Parallel Conductors. Appl. Sci. 2021, 11, 6920. [Google Scholar] [CrossRef]
  11. Coufal, O. Transient Current Density in a Pair of Long Parallel Conductors. In Prime Archives in Applied Sciences; Henninger, H., Suzdaltsev, A., Eds.; Vide Leaf: Hyderabad, India, 2022; ISBN 978-93-92117-54-1. Available online: https://videleaf.com/product/prime-archives-in-applied-sciences/ (accessed on 2 June 2023).
  12. Coufal, O. Eddy currents in a group of long solid parallel conductors. In Horizons in World Physics; Reimer, A., Ed.; Nova Science Publisher: New York, NY, USA, 2021; Volume 305, pp. 195–222. [Google Scholar] [CrossRef]
  13. Coufal, O. Current density in a pair of solid coaxial conductors. Electromagnetics 2007, 27, 299–320. [Google Scholar] [CrossRef]
  14. Coufal, O. On Resistance and Inductance of Solid Conductors. J. Eng. 2013, 2013, 526072. [Google Scholar] [CrossRef]
  15. Coufal, O. Study on Resistance and Inductance of Solid Conductors. In Technological Innovation in Engineering Research 2; Mahul, N., Ed.; Book Publisher International: Hooghly, India, 2022. [Google Scholar] [CrossRef]
  16. Coufal, O. Current density in a long solitary tubular conductor. J. Phys. A Math. Theor. 2008, 41, 145401. [Google Scholar] [CrossRef]
  17. Coufal, O. Comments on skin effect in solitary solid tubular conductor. Adv. Math. Phys. 2011, 2011, 983678. [Google Scholar] [CrossRef]
  18. Coufal, O. Current density in two solid parallel conductors and their impedance. Electr. Eng. 2014, 96, 287–297. [Google Scholar] [CrossRef]
  19. Freitas, D.; Neves, M.G.; Almeida, M.E.; Machado, V.M. Evaluation of the longitudinal parameters of an overhead transmission line with non-homogeneous cross section. Elect. Power Syst. Res. 2015, 119, 478–484. [Google Scholar] [CrossRef]
  20. Coufal, O. Faraday’s law of electromagnetic induction in two parallel conductors. Int. J. Appl. Electromagn. Mech. 2017, 54, 263–280. [Google Scholar] [CrossRef]
  21. Coufal, O.; Radil, L.; Toman, P. Magnetic field and forces in a pair of parallel conductors. Int. J. Appl. Electromagn. Mech. 2018, 56, 243–261. [Google Scholar] [CrossRef]
  22. Coufal, O. On inductance and resistance of solitary long solid conductor. Acta Tech. 2012, 57, 75–89. [Google Scholar]
  23. Rektorys, K.; Vitasek, E. (Eds.) Survey of Applicable Mathematics; Kluwer: Dordrecht, The Netherland, 1994. [Google Scholar]
  24. Mohos, A.; Ladányi, J.; Divényi, D. Methods to ascertain the resistance of stranded conductors in the frequency range of 40 Hz–150 kHz. Electr. Power Syst. Res. 2019, 174, 105862. [Google Scholar] [CrossRef]
  25. Szczegielniak, T.; Jabłoński, P.; Kusiak, D.; Pia̧tek, Z. Eddy Currents Induced in Two Parallel Round Conductors. Appl. Comput. Electromagn. Soc. J. 2019, 34, 1922–1930. [Google Scholar]
  26. Jabłoński, P.; Szczegielniak, T.; Kusiak, D.; Pia̧tek, Z. Analytical–Numerical Solution for the Skin and Proximity Effects in Two Parallel Round Conductors. Energies 2019, 12, 3584. [Google Scholar] [CrossRef]
  27. Jabłoński, P.; Kusiak, D.; Szczegielniak, T.; Pia̧tek, Z. The Proximity Effect in Twin Line with Round Conductors Placed in Conductive Medium. Energies 2020, 13, 6087. [Google Scholar] [CrossRef]
  28. Jabłoński, P.; Kusiak, D.; Szczegielniak, T. Analytical-Numerical Approach to the Skin and Proximity Effect in Lines with Round Parallel Wires. Energies 2020, 13, 6716. [Google Scholar] [CrossRef]
  29. Lide, D.R. (Ed.) CRC Handbook of Chemistry and Physics, 88th; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  30. Brandão Faria, J.A. Electromagnetic Foundations of Electrical Engineering; Wiley: Chichester, UK, 2008. [Google Scholar]
  31. International Electrotechnical Vocabulary—Part 121: Electromagnetism, IEC 60050-121: 1998. Available online: https://webstore.iec.ch/publication/66847 (accessed on 29 October 2010).
  32. Miranda, E.N. A simple model for understanding the skin effect. Int. J. Electr. Eng. Educ. 1999, 36, 31–36. [Google Scholar] [CrossRef]
  33. Reitz, J.R.; Milford, F.J.; Christy, R.W. Foundations of Electromagnetic Theory; Addison-Wesley Publishing Company, Inc.: San Francisco, CA, USA, 1993. [Google Scholar]
  34. Brandão Faria, J.A. Commemts on ’Current density in two parallel cylindrical conductors and their inductance’. Electr. Eng. 2018, 100, 1535–1536. [Google Scholar] [CrossRef]
  35. Brito, A.I.; Machado, V.M.; Almeida, M.E.; Neves, M.G. Skin and proximity effect in the series-impedance of three-phase underground cables. Elect. Power Syst. Res. 2016, 130, 132–138. [Google Scholar] [CrossRef]
  36. Stratton, J.A. Electromagnetic Theory; McGraw-Hill: New York, NY, USA; London, UK, 1941. [Google Scholar]
  37. Jackson, J.D. Classical Electrodynamics; John Wiley & Sons: New York, NY, USA, 1975. [Google Scholar]
  38. Haňka, L. Electromagnetic Theory; SNTL, ALFA: Prague, Czech Republic, 1975. (In Czech) [Google Scholar]
  39. Tamm, I.E. Fundamentals of the Theory of Electricity; Mir Publisher: Moscow, Russia, 1979. [Google Scholar]
  40. Purcell, E.M. Electricity and Magnetism; Berkeley Physics Course-Volume 2; Mc-Graw-Hill, Inc.: New York, NY, USA, 1985. [Google Scholar]
  41. Heaviside, O. Electromagnetic Theory; AMS Chelsea Publishing: Providence, RI, USA, 2003; Volume 1. [Google Scholar]
  42. Feynman, R.P.; Leighton, R.B.; Sands, M. The Feynman Lectures on Physics; Basic Books: New York, NY, USA, 2010; Volume II. [Google Scholar]
  43. Mayer, D. Applied Electromagnetism; KOPP: České Budějovice, Czech Republic, 2012. (In Czech) [Google Scholar]
  44. Sarkar, T.K.; Djordjević, A.R. Closed-Form Formulas for Frequency-Dependent Resistance and Inductance per Unit Length of Microstrip and Strip Transmission Lines. IEEE Trans. Microw. Theory Tech. 1994, 42, 241–248. [Google Scholar] [CrossRef]
  45. Sarkar, T.K.; Djordjević, A.R. Wideband Electromagnetic Analysis of Finite Conductivity-Cylinders. Prog. Electromagn. Res. 1997, 16, 153–173. [Google Scholar] [CrossRef]
  46. Djordjević, A.R.; Stojilović, M.; Sarkar, T.K. Closed-Form Formulas for Frequency-Dependent Per-Unit-Length Inductance and Resistance of Microstrip Transmission Lines That Provide Causal Response. IEEE Trans. Electromagn. Compat. 2014, 56, 1604–1612. [Google Scholar] [CrossRef]
  47. Higgins, T.J. Formulas for the inductance of rectangular tubular conductors. Electr. Eng. 1941, 60, 1046–1050. [Google Scholar] [CrossRef]
  48. Matick, R.E. Transmission Lines for Digital and Communication Networks; IEEE Press: New York, NY, USA, 1995. [Google Scholar]
  49. Simonyi, K. Foundations of Electrical Engineering; Part III; VEB Verlag: Berlin, Germany, 1956; (In Germany). Pergamon: Oxford, UK, 1963; (In English). Mir Publishers: Moscow, Russia, 1964. (In Russian) [Google Scholar]
  50. Vujević, S.; Boras, V.; Sarajcev, P. A novel algorithm for internal impedance computation of solid and tubular cylindrical conductors. Int. Rev. Electr. Eng. 2009, 4, 1418–1425. [Google Scholar]
  51. Lovrić, D.; Boras, V.; Vujević, S. Accuracy of approximate formulas for internal impedance of tubular cylindrical conductors for large parameters. Prog. Electromagn. Res. M 2011, 16, 171–184. [Google Scholar] [CrossRef]
  52. Brandao Faria, J.A.M. A matrix approach for the evaluation of the internal impedance of multilayered cylindrical structures. Prog. Electromagn. Res. B 2011, 28, 351–367. [Google Scholar] [CrossRef]
  53. International Electrotechnical Vocabulary—Part 131: Circuit theory, IEC 60050-131: 2002. Available online: https://webstore.iec.ch/publication/172 (accessed on 22 December 2011).
Figure 1. Solution to Equation (6) in Example 1. Dependence of amplitude, J ^ sol , on the distance, r, from the conductor axis for f = 10 3 , 10 4 , and 10 5 Hz; dependence of initial phase, ε , on r for f = 10 3 Hz (red line) and for f = 10 4 Hz (blue lines).
Figure 1. Solution to Equation (6) in Example 1. Dependence of amplitude, J ^ sol , on the distance, r, from the conductor axis for f = 10 3 , 10 4 , and 10 5 Hz; dependence of initial phase, ε , on r for f = 10 3 Hz (red line) and for f = 10 4 Hz (blue lines).
Applsci 13 12416 g001
Figure 2. Dependence of instantaneous current density on r for t = θ T ; θ = 0 , 0.3 , 0.4 , and 0.5, in Example 1. The amplitude, J ^ sol , is indicated by the red line and the quantity, δ , is determined by Formula (8). For r [ 0 , 3 ) mm, the current density magnitude is smaller than 6.943 A·m 2 .
Figure 2. Dependence of instantaneous current density on r for t = θ T ; θ = 0 , 0.3 , 0.4 , and 0.5, in Example 1. The amplitude, J ^ sol , is indicated by the red line and the quantity, δ , is determined by Formula (8). For r [ 0 , 3 ) mm, the current density magnitude is smaller than 6.943 A·m 2 .
Applsci 13 12416 g002
Figure 3. Dependence of the current density amplitude on x and y in Example 2 in a quarter of the conductor cross section; δ is defined by (8).
Figure 3. Dependence of the current density amplitude on x and y in Example 2 in a quarter of the conductor cross section; δ is defined by (8).
Applsci 13 12416 g003
Figure 4. Instantaneous current density, J rec ( x , y , t ) , for t = 0.39 T in Example 2.
Figure 4. Instantaneous current density, J rec ( x , y , t ) , for t = 0.39 T in Example 2.
Applsci 13 12416 g004
Figure 5. Dependence of amplitude, J ^ coa , on r in Example 3, for three values of Δ r .
Figure 5. Dependence of amplitude, J ^ coa , on r in Example 3, for three values of Δ r .
Applsci 13 12416 g005
Figure 6. Equivalent schematic of a segment of two conductors between the planes z = z 1 and z = z 2 . U L ( t ) is the induced voltage, U R ( t ) is the voltage across the segment of the two conductors. V ( z 1 , t ) and V ( z 2 , t ) are the voltages across the conductors for z = z 1 and z = z 2 , respectively.
Figure 6. Equivalent schematic of a segment of two conductors between the planes z = z 1 and z = z 2 . U L ( t ) is the induced voltage, U R ( t ) is the voltage across the segment of the two conductors. V ( z 1 , t ) and V ( z 2 , t ) are the voltages across the conductors for z = z 1 and z = z 2 , respectively.
Applsci 13 12416 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Coufal, O. One Hundred and Fifty Years of Skin Effect. Appl. Sci. 2023, 13, 12416. https://doi.org/10.3390/app132212416

AMA Style

Coufal O. One Hundred and Fifty Years of Skin Effect. Applied Sciences. 2023; 13(22):12416. https://doi.org/10.3390/app132212416

Chicago/Turabian Style

Coufal, Oldřich. 2023. "One Hundred and Fifty Years of Skin Effect" Applied Sciences 13, no. 22: 12416. https://doi.org/10.3390/app132212416

APA Style

Coufal, O. (2023). One Hundred and Fifty Years of Skin Effect. Applied Sciences, 13(22), 12416. https://doi.org/10.3390/app132212416

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop