3.1. The Method of Moments: Attributes and Limitations
One approach to modeling a substation, as suggested by CIGRE, is to use a large number of dipoles [
1]. However, surfaces, such as the metallic shielding surfaces that surround a valve hall, are not readily modeled merely by dipoles. Given this, several authors have modeled substations using the method of moments (MoM), which provides a method for describing the electromagnetic fields of a collection of conductors with specified sources [
12,
14,
15,
16,
17,
18,
28,
29,
30]. The conductors are divided into small segments (effectively acting like dipoles) or surface patches (thus permitting the modeling of conducting enclosures such as those associated with a valve hall). The MoM is inherently a frequency-domain technique and thus the transient nature of the sources in an HVDC station must be accounted for by specifying the correct magnitude and phase of the spectral components of the excitation and transforming them to the time domain. The MoM does have several compelling attributes in terms of its modeling abilities. First, it permits the simple and accurate modeling of thin conductors such as those one would find throughout an HVDC converter station. Additionally, one can easily incorporate lumped elements at a point on a line or distribute a load over some portion of a line. However, caution must be exercised when modeling using a lumped circuit, since the modeling may require that the physical elements of the structure responsible for its lumped circuit behavior be eliminated. Finally, one need only determine the currents that flow on wires and surfaces. Once these are obtained, the fields can be obtained throughout the space.
MoM was used by Dallaire and Maruvada [
11] to show that it is critically important to create a full-wave model of the enclosure with the associated conductors passing through the shielding. As mentioned, they showed that there are situations in which the radiated field in a particular direction may be
enhanced because of the shielding enclosure [
11]. As also mentioned, this contrasts, for example, with the statements in [
5,
12] that indicate a valve hall provides at least 40 dB of reduced radiation. Again, the statements in [
5,
12] are based on measurements and the authors do not provide a full-wave analysis of an electromagnetic environment representative of an actual valve hall.
Because MoM has frequently been used to model emissions from converter stations, the authors decided to explore the degree to which the most recent version of the MoM-based NEC code, i.e., NEC-5, could be expected to yield meaningful results for scenarios related to HVDC converter stations [
24,
25]. More specifically, a simple test was created where a dipole source was enclosed in a simplified valve hall represented by a hollow metallic cube with 10 m sides. The dipole was aligned with the axes of the cube and centered in two of the dimensions but offset from the center by 3 m in one of the dimensions. The valve hall was “perfect”, in that there were no holes in its walls nor conductors passing from the inside to the outside. The magnitude of the electric field was recorded over a circle that was 200 m in radius and centered about the
z axis. If the numeric code was performing perfectly, the exterior fields would be identically zero at all points throughout the space exterior to the cube. However, that was not the result that was obtained.
Since the MoM is a numerical method and the currents on the patches that represent the walls of the enclosure are discretized, the currents are always approximations and one would not expect perfect cancellation beyond the enclosure. Realistically, the best that one could hope for is that the exterior fields would be “small”, but what constitutes “small”, and the degree to which MoM can achieve that, appears open to question. In many applications, when comparing numerical results to either theoretical or measured results, an error of one percent is quite acceptable. Another way of thinking of one percent is as two orders of magnitude or 40 dB. But in the scenario described above, with a source in a perfect valve hall, where the exterior field should be zero, the difference between any non-zero field and the correct value is infinite on a dB scale. Alternatively, one might compare the “numerically leaked” field at the exterior to the field when no valve hall is present. This was carried out for excitation frequencies of 200 kHz, 1 MHz, 5 MHz, and 25 MHz.
Table 1 shows the reduction of the measured field caused by the valve hall compared to that when the valve hall is not present (a finitely conducting ground beneath a perfect valve hall is present for all calculations with
and
mS/m).
Because the hall has no apertures, the correct change is
on a dB scale, i.e., the lower the value the better. At 200 kHz, we note that the reduction is in the order of 40 dB, which is the value often given as the effect of the valve hall. Recall that a 40 dB reduction appeared in [
5,
12] and was based on measurement. However, the approximate 40 dB reduction observed at 200 kHz with this NEC-5 simulation is purely a numeric artifact. The fact that this may agree with the measurements should be considered coincidental and not indicative of this software capturing the full physics of the effect of the valve hall. In [
25], the authors also report a “shielding effectiveness” at 15 MHz of
dB. That matches the value of 5 MHz in
Table 1 remarkably closely. But, again, this is mere coincidence. If one continues to higher frequencies (up to 100 MHz), ref. [
12] reports that the shielding effectiveness further increases to
dB, whereas the MoM results mysteriously drop to
dB at 25 MHz (the highest frequency that was considered). In contrast to the behavior at these frequencies, the results obtained at 25 MHz are unacceptably high (predicting only about
dB of shielding from the perfect valve hall) and apparently nearly immune to an increase in the discretization of the valve hall walls. This is a cause for concern and at present we have no suitable explanation for this behavior.
The results in
Table 1 were obtained using a discretization that adheres to the recommendations in the NEC-5 user’s manual [
25]. Nevertheless, keeping in mind that any non-zero field represents an error and observing that the results ran into difficulties at 25 MHz, the test was run using a finer discretization. The results in
Table 1 were obtained using 400 patches per face of the cube (i.e., 20 patches along a side of the cube for
patches). There was no significant difference in the results.
3.2. The Finite-Difference Time-Domain Method: Initial Insights
from Simple Geometries
The finite-difference time-domain (FDTD) method [
31,
32] employs a fundamentally different approach. One must discretize not only the sources themselves, but also the space in the vicinity of the sources and the space that surrounds any physical features that might affect the field of interest. In addition to discretizing space (where the discretization assigns the appropriate material properties to each point in the discretization), the FDTD method also discretizes time. In this way, all the derivatives in Maxwell’s equations (specifically the curl equations) are replaced with central differences. This yields a set of equations where future (unknown) fields are expressed in terms of past (known) fields. One then advances by incremental steps in time, revealing the time-domain behavior of the fields. While the MoM is a frequency-domain technique that allows one to obtain fields throughout space at a single frequency, the FDTD method is a time-domain technique that can potentially yield results over a broad spectrum with a single simulation. However, the FDTD method will only directly yield the fields within the space that has been discretized. If one is interested in fields outside of the discretized space, one must carry out a transformation of the known fields to the point of interest. The FDTD method is generally considered to be quite computationally demanding. Furthermore, special consideration is needed to handle any number of features of a particular scenario, such as a relevant physical structure that is small compared to the discretization that is used.
Several FDTD simulations were run using code developed in-house that explored various aspects of HVDC converter station modeling with an emphasis on the effect of the valve hall. The first considered was a “perfect” valve hall. As with the initial MoM testing, this hall had no openings or any conductors running from the interior to the exterior. However, instead of using a cube with 10 m sides, the dimensions were based on those of an actual converter station. The valve hall dimensions were (x) 42 m by (y) 68 m by (z) 19 m. The hall and surrounding space were discretized with a step size of m. Initial simulations, which incorporated part of the AC yard and the DC hall (reactor hall), used a computational domain of cells in the x, y, and z directions, respectively.
The computational domain is depicted in
Figure 1, where the scale is in terms of FDTD cells (one can divide by two to convert to meters). Red lines correspond to conductors and the rectangular structure near the center of the figures corresponds to the valve hall. Conductors are modeled as a line of electric field nodes set to zero. Several researchers have reported that the “effective” radius
of such a line of cells is approximately
[
33,
34,
35,
36,
37], ranging from
[
33] to
[
35]. In [
33,
35,
36,
37], the researchers proposed ways to modify the FDTD update equations to model wires that had radii that were either smaller or larger than this effective radius. In [
38], Taniguchi et al. reexamined the effective radius in the context of surge impedance and compared FDTD results to the theoretical work of Chen [
39]. They obtained an effective radius of
. This is taken to be the most accurate value of the effective radius. Applying that to the modeling performed here yields an effective conductor radius of roughly 10 cm.
The green plane on which the valve hall sits is a perfectly conducting ground plane.
Figure 1a shows its projection onto the
plane (the green ground plane has been removed from this view). The AC yard is to the left (west) of the valve hall. The three lines in the AC yard ultimately are shorted together at what would be, in practice, the location of a transformer. Although this is a crude approximation of the electromagnetic behavior of the transformer at these frequencies, it is not an unreasonable one given the lack of a validated model of such a transformer in the open literature. Nevertheless, because of the shorting of the lines at this location, the model presented here could not predict the high-frequency currents on the AC lines exiting the AC yard.
The conductors to the right (east) of the valve hall correspond to several of the lines that exist in the DC hall. The walls of the DC hall are not modeled. The 18 rectangular conductors seen within the valve hall serve as coarse approximations of the valves themselves. No attempt was made to model the detailed physical structure of the valves, but rather we modeled their overall size and shape to approximate how electromagnetic energy might propagate in their presence. A horizontal line passes through the center of each of the six sets of three “valves”, moving from the AC to the DC side of the hall.
Figure 1b shows a projection of the computational domain onto the
plane while
Figure 1c shows an oblique perspective projection of the computational domain.
Energy is introduced into the computational domain by specifying the electric field at a particular FDTD node along one of the conductors. Specifically, the value of an
node along a line that transits from the AC to the DC side of the hall is specified by a source function at each time step of the simulation. This is the equivalent of specifying the voltage at that location. The source function was a unit-amplitude Ricker wavelet. (A Ricker wavelet is obtained by taking the second derivative of a Gaussian. In this way, the Ricker wavelet has somewhat similar spectral properties to a Gaussian pulse, but no DC component.) The Ricker wavelet had its peak spectral content at 10 MHz. This source can be thought of as introducing energy over a range of frequencies that can be associated with commutation events. But, the simulations are performed with a single “firing” of the Ricker pulse at a single location. For this analysis, the source location is near the right edge of the center “valve” at the northern-most (greatest
y) side of the valve hall, as shown in
Figure 1a.
Although the computational domain is finite, an infinite space is approximated by surrounding the computational domain with a perfectly matched convolutional layer [
40]. Simulations were typically run for
time steps. Using the computational domain mentioned above (
cells), simulations took approximately one hour and 49 min to complete on Apple Mac Studio with an M1 Ultra processor and 128 GB of memory. For the discretization that was used, where the time step
was
ps,
time steps correspond to a total of 63.06
s. (The Courant number
was the 3D limit of
.)
In the FDTD method, boundary conditions at material interfaces “take care of themselves”, as one merely specifies the material properties at all points throughout the discretized space, i.e., throughout the FDTD grid, and then advances in time via the discretized version of Maxwell’s equations. However, perfect electric conductors (or nearly perfect conductors), such as those that pertain to the metal walls of a valve hall, are realized not by specifying material values per se, but rather by setting the electric fields tangential to such conductors to zero and leaving them at zero. (This is similar to the way in which conductors are realized, where a line of electric field nodes is set to zero.)
Simulations were performed of a perfectly intact valve hall where the conductors (red lines) shown in
Figure 1 were present but terminated by unbroken conducting walls. The FDTD model of such an intact valve hall yields results consistent with our expectations: (1) In contrast to the MoM solution, the fields exterior to the hall are identically zero. (2) Once energy is introduced into the valve hall, either in the form of fields that radiate away from the source or currents induced on the line in which the source is embedded (and the fields associated with such currents), this energy simply oscillates inside the hall until the simulation is terminated. (3) Strong resonances occur owing to the structure of the hall and its associated interior conductors. When the model consists solely of perfectly conducting wires and walls, there is no loss mechanism and the fields resonate indefinitely in this rectangular cuboid (The simulation does allow for “loss” in the form of radiation and subsequent absorption in the perfectly matched layer that surrounds the computational domain. However, such a loss can only take place if fields escape the valve hall. And, with the hall intact, no fields escape.).
Nine holes were then introduced into the valve hall walls corresponding to the locations of each of the lines seen in
Figure 1 that transitioned from the interior to the exterior of the valve hall (three on the AC side of the hall and six on the DC side). The size of the holes was not the minimum possible, but rather the minimum possible that would subsequently let us model a conductor passing through the hole. For the somewhat coarse discretization that was used (
m), the hole could be said to be approximately 1 m in diameter. Given that the peak spectrum of the pulsed excitation was at 10 MHz, i.e., a wavelength of 30 m, the holes were small compared to the wavelength of the bulk of the excitation.
With the holes thus created, multiple simulations were performed with various modifications. First, in addition to introducing the holes, a 2 m section of each conductor passing through the walls was removed which was centered about the valve hall walls, i.e., a 1 m section of each conductor was removed to either side of the wall it passed through. Second, the conductors were left intact and merely passed through the holes. For both simulations, the magnitude of the electric field at two points was recorded, one point inside the hall and one point outside. These points were selected somewhat arbitrarily. Relative to the hole that was introduced to accommodate the northern-most line in
Figure 1a, which passes from the valve hall to the DC hall, the “inside point” was 4 m away in
x (i.e., 4 m further into the hall relative to the wall with the hole), 3 m away from the hole in
y (3 m closer to the center of the hall), and 2 m below the hole in the
z direction. The outside point was 4 m away in
x (i.e., 4 m away from the wall), but, as with the inside point, 3 m away in
y and 2 m away (below) in
z. Thus, each observation point was approximately
m away from the center of the hole. This northern-most line is also the one along which the commutation event is simulated.
Figure 2a shows the field measured (for the remaining discussion of the results obtained via the FDTD method, measured and measurement should be considered synonyms for calculated and calculation) at the inside and outside points when the 2 m segments are removed from each line. The fields are plotted on a log scale over the first 5
s of the simulation. Note that the field inside the hall becomes non-zero prior to the field outside the hall due to its closer proximity to the source. The field inside the hall varies in the vicinity of roughly
db, while the field outside the hall varies in the vicinity of roughly
db, i.e., a reduction in the field of approximately 50 db.
Figure 2b shows the field at the same two points, using the same excitation and the same geometry, but with the conductors that transition through the holes in the valve hall walls now intact. Again, the field arrives at the exterior point after arriving at the interior point. However, in this case, the fields are almost comparable in magnitude. The inside field is generally larger, as would be anticipated, but only slightly so, and there are instances where the field at the point outside the hall is greater than inside the hall. Another way to think about the contrast between
Figure 2a,b is that having intact conductors raises the field exterior to the hall by approximately 50 dB.
The simulation in
Figure 2b merely had the conductors passing through holes surrounded by free space. Another simulation was performed that modeled the bushings associated with each hole. These bushing were assumed to be cylindrical, extend 1 m to either side of the wall, and have a relative permittivity,
, of 3. The radius of the bushings was such that they filled the holes through which the conductors passed. The bushings made essentially no difference to the overall observed behavior seen without the bushings. Another simulation was performed where, in addition to having the bushings present, lumped-element inductors were incorporated into the FDTD grid to represent the smoothing reactors found in the DC hall. The inductance used was
mH (in keeping with the inductance provided in [
21]; its implementation followed the work described in [
41]). Again, this made effectively no difference to the roughly similar field observed inside and outside the valve hall.
In
Figure 2a, note that after about 1
s, the fields decrease very little. Contrast this to the slight, but clearly evident, decrease in the field as time progresses in
Figure 2b. In both simulations there are nine holes in the valve hall walls, but only in the simulation for
Figure 2b do conductors pass through the holes. In the case of
Figure 2a, given the relatively small size of the holes compared to the wavelength, the fields escaping the hall are primarily evanescent. These fields can couple onto the conductors that are present outside the hall, but this represents a nearly vanishingly small radiation of the fields, which explains why there is seemingly no decay (with time) of the fields in
Figure 2a. (
Figure 2 only displays the first 5
s of the simulation, but this persistence of the field is evident for the entire 66.03
s simulation.) On the other hand, when the conductors are intact, they serve to create what is effectively a transverse electromagnetic (TEM) transmission line for the fields to pass from the interior to the exterior of the valve hall. There is no cut-off frequency for this transmission line (which needs to be the case to move DC and low-frequency energy in and out of the hall).
A conductor above the ground plane serves as a TEM transmission line. The lines in the AC yard can be thought of as segments of such lines. As a line passes through the valve hall wall, one might envision this transit as a very short co-axial line, where the conductor is the inner core of the co-axial line and the wall is the outer conductor. There is an impedance mismatch going from one form of transmission line to another, but given the short extent of the intervening “co-axial” line, this mismatch does not cause significant reflections, i.e., fields directed out of the hall largely continue to exit the hall, a point which will be returned to when the currents in the conductors are considered.
The
component of the electric field at these observation points when the conductors are intact will now be considered.
Figure 3a shows the log magnitude (which is twenty times the log base 10 of the magnitude) of the Fourier transform of the
component of the electric field, measured at the points inside and outside the hall over a range of frequencies from 0 to 35 MHz. For reference, the Fourier transform of the source function is also provided. One can see that peak energy occurs at 10 MHz. Note that below approximately 14 MHz, there is very little difference between the fields inside and outside the hall and, at several frequencies, the field outside the hall is greater than that inside the hall. For frequencies greater than 14 MHz, the field inside the hall is generally larger than that outside the hall. This appears to indicate that the valve hall is serving as a filter. However, such an assertion can be called into question, as discussed below. Note that the lowest-order resonant frequency for a valve hall of this size is
MHz. The strongest resonances generally occur below this frequency. Thus, these resonances must be associated with the overall structure of the station, including the conductors in the AC yard and DC hall, and not merely with the valve hall itself. As will be seen, these “low-frequency” resonances remain rather persistent regardless of what is done to the valve hall.
Consider the fields normalized by the spectral content of the source function, as shown in
Figure 3b. Here, the similarity of the fields below 14 MHz and their apparent filtering above 14 MHz are somewhat easier to discern. However, also consider the results shown in
Figure 3c, which are identical to those of
Figure 3b, except that now the valve hall has been removed. (These curves will be labeled “Inside” and “Outside”, although in
Figure 3c there is no valve hall present to distinguish inside from outside while the location of the observation points have not changed.) Several of the peaks, especially those below about 7 MHz, remain nearly unchanged. The persistence of these resonances, independent of the existence of the valve hall, indicates that they are associated with the conductors throughout the model. For frequencies above this, the results in
Figure 3c are much “quieter” than they are when the valve hall is present. Consider the “outside” curves in
Figure 3b,c; although the curve in
Figure 3b is much noisier than that of
Figure 3c, the fluctuations in
Figure 3b appear to be roughly about the “baseline” shown in
Figure 3c. The same cannot be said of the results for the “inside” curves. It is certainly true that the results for the inside point in
Figure 3c are quieter than those for
Figure 3b, but one would not say that the fluctuations in
Figure 3b are roughly about the curve present in
Figure 3c. Instead, when the valve hall is present, rather than filtering the outside field, its presence tends to elevate the interior field.
Putting these observations together, we can say that the valve hall does not diminish the exterior field but rather enhances the interior field. Note, however, that this claim pertains to these two particular observation points, and both are in relatively close proximity to the line on which the source pulse exists. Were one to consider an exterior point located close to one of the two walls with no conductors passing through them, the field measured with the presence of the valve hall would certainly be reduced relative to having no hall present. Nevertheless, the interior field is enhanced by the presence of the valve hall (independent of where the field is observed in the interior). This enhancement may not be evident at a particular time or at a particular frequency. Rather, if one integrates the square of the field over time, the overall value will be higher when the hall is present due to the way in which it hinders the radiation of energy in all directions except via the transmission lines that pass out of the hall. Because there is no loss present within the valve hall, the only place energy can go after it is introduced into the valve hall is out along the conductors that pass through the valve hall walls, i.e., along the transmission lines out of the hall. At these frequencies, these lines serve as antennas which then radiate fields into the surrounding environment (as will be considered below).
Instead of a component of the electric field, consider the current on a conductor measured both inside and outside the valve hall. The current is obtained by taking the line integral of the magnetic field around the conductor at the desired measurement point. The current was calculated on the same line on which the source exists (the northern-most line exiting the valve hall to the right in
Figure 1a). The inside and outside measurement points are each 2 m from the wall (with one being inside the valve hall and the other being in the DC hall). Temporal plots of the currents at these points reveal a difference in the arrival time of fluctuations, which is consistent with the different locations of the observation points. There are some other differences in the currents, but these are relatively minor. Rather than showing a temporal plot,
Figure 4a shows the log magnitude of the normalized spectrum of the currents when the valve hall is present.
Figure 4b is the same scenario as
Figure 4a, except the valve hall has been removed. Although there are some slight differences in the features observed in
Figure 3 and
Figure 4, overall, the results in
Figure 4 serve to confirm the statement made in connection to the measurement of the electric field: the valve hall does practically nothing in terms of confining high-frequency energy to the interior of the valve hall. One could say that the hall serves to redirect energy while not eliminating it (i.e., delaying energy’s eventual escape from the hall). Once a current has been established on a line within the valve hall, if that line exits the valve hall, the current found outside the hall will be nearly the same as that inside.
Having established that the valve hall does little in terms of reducing the ultimate escape of high-frequency currents and fields, attention is turned to measurements more aligned with those related to compliance. CIGRE [
1] is frequently concerned with the fields that exist “200 m from the closest active part of the station”. The computational domain was extended 200 m (400 cells) in the
y direction (simulations of
time steps now required approximately five and half hours to complete on the Apple Mac Studio mentioned before) and the field was measured at several observation points that were 200 m “south” of the southern-most valve hall wall depicted in
Figure 1a (i.e., 400 cells away from the valve hall wall in the
y direction). Recall that the entire computational domain sits atop a perfectly conducting ground plane. All components of the electric and magnetic fields were measured at 15 different locations with varying
x and
z values. The observation points were separated 15 m in
x and had vertical displacements (
z values) either 2, 3, or 4 m above the ground plane. In the discussion that follows, any, or all, of these points could have been selected and would lead to the same observations. For the sake of concreteness, however, an observation point was selected that is
m to the “west” of the western-most valve hall wall (i.e., the one leading to the AC yard), 200 m south of the southern-most wall, and 2 m above the ground plane.
Figure 5a shows the electric field magnitude over the initial 5
s of two simulations, one when the valve hall is present and the other where all the conductors are unchanged but the valve hall has been removed. Unsurprisingly, when the valve hall is not present, the initial fields are stronger than when the hall is present (please refer to the interval between 1 and 2
s). This is a consequence of the fields being free to directly radiate in all directions (other than “down” into the perfectly conducting ground) from the commutation event. However, after this initial burst, the fields without the valve hall drop below those observed when the hall is present. This is because without the hall present much of the energy is quickly lost to radiation throughout the surrounding environment. On the other hand, with the hall present (and there being no loss mechanisms contained within the hall) and the energy unable to radiate in all directions, the energy confined to the hall can only escape slowly via the conductors which, in turn, serve as antennas to radiate the field to the surrounding environment. As previously mentioned, including inductances (which model the smoothing reactors) and bushings in the simulation has little effect on the radiated field (affecting mostly its phase but not its magnitude).
Two other scenarios were considered. In one, the walls of the valve hall remained in place, but the roof was removed. Potentially, this could have directed significant amounts of energy “up” so that this energy would not be present for terrestrial measurements, such as at the observation point being considered. However, at the frequencies being considered, it appears that the height of the valve hall walls was not sufficient to direct the energy up. Instead, the fields largely diffract over the walls so that the observed fields, in this case, are similar to having no valve hall present at all, and thus these results are not shown here.
The other scenario we considered introduced absorption loss into the valve hall. This took the form of a 2 m thick (in the
z direction) dielectric slab adjacent to the ceiling of the intact valve hall. The slab nearly spanned the
x and
y extents of the hall and was displaced from the walls and roof by
m. The slab had a relative permittivity,
, of
and a conductivity,
, of
mS/m.
Figure 5b again shows the electric field magnitude at the observation point over the initial 5
s of two simulations. The results are again shown for an intact valve hall (the same as in
Figure 5a) but also for the case in which the lossy dielectric slab is in the valve hall. In this case, the initial fields in the interval from 1 to 1.5
s are nearly identical, but with the lossless case having slightly higher peaks. However, beyond that, the field for the lossy hall decreases significantly relative to the lossless case.
Figure 6 shows the log magnitude of the normalized spectrum of the
x component of the magnetic field for the three scenarios considered in
Figure 5. (If
is scaled by the impedance of free space, nearly identical results are obtained compared to plotting the spectrum of
. In practice, the horizontal component of the magnetic field is measured using a loop antenna, which is what motivated the plotting of
in
Figure 6). Across much of the spectrum, when the lossy slab is present, the field is smaller than in the other two scenarios and, over certain ranges of frequencies, significantly so. Above approximately 14 MHz, the lossy results are below those of the intact hall by roughly 10 dB. Between 7 MHz and 14 MHz, the difference between these results varies, but the lossy hall consistently out-performs the intact hall in terms of having lower fields. Below 5 MHz, the results are all quite similar. Again, this range of frequencies is largely below the lowest-order resonance for a cuboid the size of the valve hall.