Optimization of Lumped Parameter Models to Mitigate Numerical Oscillations in the Transient Responses of Short Transmission Lines

: The Lumped Parameter Model (LPM) is a known approach to represent overhead transmission lines (TLs), especially when these elements comprehend a few tens of kilometers. LPMs employ a large number of cascaded π -circuits to compute accurately the transient responses. These responses contain numerical spurious oscillations (NSO) characterized by erroneous peaks which distort the transient responses, mainly their peak values. Two modiﬁed LPM topologies composed of damping resistances inserted along the longitudinal or transversal branches of the cascaded π -circuits offer signiﬁcant mitigations in the NSO. In this paper, in an effort to have the maximum mitigation of the NSO and low distortion in the transient responses, two modiﬁed topologies with optimized damping resistances are proposed to represent short TLs. Results demonstrate expressive attenuation in the peaks of NSO which reﬂect good agreement in comparison with the responses computed by the Bergeron’s line model. The mitigation of the NSO is carried out directly in the time domain and it does not require either analog or digital ﬁlters.Furthermore, no frequency-to-time transformations are necessary in this procedure. These alternative topologies can be incorporated into any electromagnetic transient program to study switching operations in power systems.


Introduction
Electromagnetic transients can be generated by switching operations, short-circuit faults, and lightning strikes on overhead transmission lines (TLs). As a consequence, the overvoltages and overcurrents produced may affect the electrical supportability of insulator strings and surge arresters, may trigger undesirable operation or malfunction on protection devices, or establish high current faults in the power system [1].
LPMs can be applied to represent short TLs when submitted to switching maneuvers, characterized by a low-frequency content [2]. In this frequency content, the TL electrical parameters can be considered constants [3]. LPMs are employed to simulate electromagnetic transients in many conditions, such as (i) for switching maneuvers (energization) for which frequency content varies from 50/60 Hz up to 20 kHz [2,4]; (ii) when nonlinear loads are connected (e.g., surge arresters, circuit breakers, preinsertion resistors) in the power system; (iii) when fault currents are in different locations along the line; (iv) when Corona effect must be included in the analysis [5]. Additionally, LPMs are inserted in most Electromagnetic Transient (EMT)-type programs, such as ATP ® (ATPDraw, Trondheim, Norway), PSCAD ® (Manitoba Hydro International Ltd., Manitoba, Canada) or EMTP-RV ® (EMTP Alliance, Montréal, Québec, Canada).
On the other hand, numerical spurious oscillations (NSO) are present in the transient responses computed with LPM [6]. The NSO occurs due to the lumped representation of the transmission line parameters, which by nature are distributed associated with numerical integration methods to compute the transient responses [5,7]. These NSO do not correctly represent the transient responses which may provide incorrect analysis and affect the adequate operation of protection devices. These erroneous peaks may lead to an overestimation in the insulation level of equipment, such as insulation strings, preinsertion resistors and circuit breakers [8,9]. The NSO depends on the number of π-circuits in the cascade, line parameters and the step size of simulations. Many methods have been used to mitigate these peaks, such as: (i) digital filters [6,10]; (ii) analog filters [9]; (iii) numerical methods [11] and (iv) modified topologies of the LPM [6,12].
Two modified LPM topologies are available in the literature to mitigate the NSO [6,12]. In these topologies, an additional damping resistance is inserted in the cascaded π-circuit which acts as a low-pass filter of the NSO. One topology is composed of a damping resistance in parallel to the RL branch of the traditional LPM. The other is composed of a damping resistance inserted in series to the shunt CG branch of the traditional LPM, as further detailed. Both these damping resistances are dependent on adjustable factors, step size, and distributed line parameters [6]. The adjustable factor is chosen randomly by the user, which can be a cumbersome task with no maximum attenuation and large distortion in the transient responses [6,13]. Besides, these modified cascaded structures keep the same advantages as the traditional LPM. In the traditional and modified LPMs, the transient voltages and currents can be expressed in terms of state-space systems, the solutions of which are computed by any numerical integration method. The main advantage of modified topologies is the direct reduction of the NSO in the time domain, without any frequency-to-time conversion tools, Laplace or Fourier transforms. These transformations are carried out by numerical methods which require a high number of samples in the frequency domain and window functions which can affect the precision of the transient responses.
In order to investigate the performance of the modified topologies in the transient responses, concerning the maximum attenuation of NSO with the lowest distortion, this paper presents a methodology to compute the best adjustable factors for each modified topology based on the LPM. These adjustable factors depend on the line length and number of π-circuits/km, where fitted 3-D surface equations are obtained using Curve Fitting in Matlab ® . The NSO are significantly mitigated in the transient responses when single-and three-phase TLs, with lengths varying from 20 km up to 80 km subjected under energization and fault conditions, are represented by these modified topologies. Simulation results demonstrate that the transient responses obtained with the optimized adjustable factors are in good agreement with those obtained with Bergeron's line model (BM) which is a traditional approach in power system analysis. The modified topologies with the best adjustable factors can be incorporated into EMT-type programs in which neither digital nor analog filters are required since the filtering process is carried out directly in the time domain. This procedure is very advantageous due frequency-to-time domain transformations such that Fourier transforms are not necessary. This paper is organized as follows: in Section 2, the traditional LPM is presented; in Section 3, two modified topologies are described, as well as the state-equations to compute the transient voltages and currents in the TL. In Section 4, the optimized methodology to compute the best adjustable factors, 3-D graphs, and fitted equations are provided. The numerical results for the mitigated NSO in the transient responses computed for singleand three-phase lines and discussions are presented in Section 5. Finally, the conclusions are summarized in Section 6.

Traditional Lumped Parameter Model
The traditional LPM has been employed to represent short TLs and cables in power systems in which the transient responses are computed directly in the time domain [5,7]. LPMs are employed when nonlinear loads, faults, or Corona effect must be considered, or when transmission line parameters must be estimated (used as an equivalent π-circuit) in the time-domain transient analysis [5,14]. However, the accuracy from LPMs is restricted to TLs within a few tens of kilometers submitted to low-frequency disturbances [15]. When switching maneuvers (energization) are involved, the TL can be considered as composed of frequency-independent parameters and the LPM constitutes a reasonable approach [2,4,8].
In practical applications, the resistance and inductance can be calculated at a standard frequency (50/60 Hz), the capacitance is assumed constant and G can be neglected in short TLs. Thus, the LPM can be employed for transient analysis [16] or in the steadystate analysis of short lines (distribution network) for power analysis [3]. The TLM based on traveling waves such as, BM, JMarti's line model, and the Universal Line Model, present limitations concerning short TLs due to the time step ∆t that must be less than the traveling time τ (τ = d/v), where d and v are the line length and the propagation velocity, respectively [17]. Furthermore, in EMT-type programs, when several TLs are involved in the transient simulation, the minimum ∆t is adopted based on the shortest line length [17]. The LPMs present the ∆t chosen independently of the τ, as seen in [5,7,18]. This fact can be used in transient simulations where ∆t is greater than τ. As an example, where short TLs present a τ smaller than ∆t, the LPM can be a reasonable approach for energization studies without excessive loss of accuracy [2]. In order to show the influence of the number of π-sections in the transient responses, a coefficient k named the discretization factor is defined as follows: where n is the number of the π-circuits to represent the TL and d is the line length (km). A generic single-phase TL is used to illustrate the LPM as depicted in Figure 1 of [9]. It is assumed that the height of the line conductors is kept constant, the conductors are parallel to the ground and the line is located above a soil represented by its frequencyconstant parameters (resistivity and permittivity). Each section is characterized by a series impedance, consisting of a lumped resistance R with an inductance L, and by a transversal (shunt) admittance, consisting of shunt lumped capacitance C and conductance G. These parameters are given as follows: where the elements R , L , C and G are the per-unit-length (p.u.l) resistance, inductance, capacitance and conductance, respectively. In this representation, a large n is required to satisfactorily represent the distributed nature of a transmission line and, as a consequence, it also provides a better transient response [16,19]. When LPMs are used to represent TLs, the transversal voltages and longitudinal currents along the cascaded π-circuit sections can be expressed in terms of a state-space system given as follows [7,18]: where in (3), x 0 (t) is the state variable vector, u(t) is the input vector, where it includes ideal current or voltage sources connected at the sending end, A 0 is the state (or system) matrix (2n × 2n) and B 0 is the input matrix (2n × 1). The matrices A 0 and B 0 and vector x 0 (t) are generically represented by: The voltages and currents along the line can be obtained by any numerical integration method to solve (3). In order to illustrate the influence of k in the transient responses, as well as the NSO peaks, a 20-km single-phase TL is employed, the configuration of which is shown in Figure 1a. In this case, the TL is energized by a 1-p.u. DC voltage source (energization maneuver) with an open-circuit receiving end. The computed line parameters are (at 60 Hz): R = 0.075 Ω/km, L = 1.722 mH/km, C = 9.638 nF/km and G is neglected. The TL is represented by n = 10, 20, 40 and 80 π-circuits which results in k = 0.5, 1, 2 and 4 π/km, respectively. It is adopted that ∆t = 1 µs for these simulations. The transient voltages are calculated by LPM and compared with the responses obtained by BM, as depicted in Figure 1b. It can be noted as k increases, the better is the voltage waveform despite the NSO peaks being more pronounced. As seen, the NSO may exceed by 25% the transient voltage peak for the highest k, which may lead to an erroneous analysis by users. The rise time is better represented as k increases; however, the computational time required to compute the transient responses is higher in comparison with those obtained at lower k. In the next section, modified topologies to mitigate the NSO are presented.

LPM with Damping Resistances in the RL Branch
A modified topology (T 1 ) in the literature to mitigate the NSO is presented in Figure 2a [6,20]. To reduce the NSO, a damping resistance R d1 is connected in parallel to the RL branch of the traditional π-circuit section, as depicted in Figure 2a. In this topology, the artificial damping resistance R d1 and its conductance G d1 , are given by: where K d is an adjustable factor that varies between 2 and 10 [6], L is the lumped inductance and ∆t is the step size adopted in the simulations. Applying Kirchhoff's circuit laws in this modified π-circuit section, the voltage in the shunt CG branch (v k ) and the current in the series RL branch (i k ) are written as follows in (6). Expanding to a TL represented by n-π-circuit sections, the generic v k (t) and i k (t) are expressed as a linear system of space-state variables, as follows: The matrices A 1 and B 1 in (7), and vector x 1 (t) are expressed by:

LPM with Damping Resistances in the GC Branch
The modified topology (T 2 ) is depicted in Figure 2b, in which a damping resistance R s1 is inserted in series to the GC branch of the π-circuit section. The NSO are mitigated with the insertion of a damping resistance R s1 given by: where C is the lumped shunt capacitance of the π-circuit, ∆t is the step size of the simulations and the K s is an adjustable factor. Applying the same procedure in Figure 2b, the current (i k ) and voltage (v k ) are computed as follows: where α 2 and β 2 are defined by (11).
The voltages v k (t) and currents i k (t) are written in a state-equation form, given by: where the matrices A 2 and B 2 and vector x 2 (t) are given by (13).
Although these two modified topologies are known, just a few papers dealing with the influence of these adjustable factors on the reduction of the NSO are presented in the literature [6]. It is worth mentioning that the adjustable factors K s and K d are selected randomly by the users which may affect significantly the NSO peaks and the rise/fall times in the transient responses. To find the best benefit concerning NSO peaks and distortion at the rise/fall times, a proposed methodology to select the best adjustable factors K * s and K * d is presented in this work.

Optimization of the Adjustable Factors
To minimize the peaks of the NSO in the transient responses, the best (optimum) adjustable factors must be computed based on the TL parameters and step size. The optimum adjustable factors are indicated by (*), as K * d and K * s in this work. For this purpose, two algorithms are presented to compute the K * d and K * s , based on the voltage peaks calculated in each topology. To illustrate these optimization methods, the 20-km single-phase TL from Figure 1a is represented by the traditional LPM and by the alternative topologies, which can be the modified T 1 or T 2 . This TL is energized by a 1-p.u. DC voltage source and the receiving end is left open. The computed voltages are illustrated in Figure 3. As noted in Figure 3, the NSO in the transient responses obtained with the alternative topology are significantly mitigated in comparison with those computed with the LPM and also closer to those obtained with BM. Aiming to obtain the optimum adjustable factors K * d and K * s , the amplitudes of the voltage peaks must satisfy the condition in (14), as follows: where A 1 and A 2 are the peaks of the transient voltages obtained for the LPM and modified topology in relation to BM, as detailed in Figure 3. The steps of the optimization algorithm are shown in the flow-chart depicted in Figure 4.  In these algorithms, the R , L , C and G are the TL per-unit-length parameters, d is the line length, and ∆t is the step time. As a numerical example, the single-phase TL from Figure 1a is represented by several discretization factors k varying from 0.5 up to 3.50 π-circuits/km and line lengths d varying from 20 up to 80 km where the 3-D graphs (abacuses) are shown in Figure 5. Based on the flow-chart in Figure 4a, the K * d in a 3-D graph is depicted in Figure 5a. The red curve represents the K * d s for a line of 50-km and varies from 0.3 up to 1.5. For low values of k, the K * d s are high due to the NSO peaks that increase as k increases in the traditional LPM (see Figure 1b). As a consequence, the K * d tends to decrease, resulting in a lower damping resistance R d1 and in a higher reduction in the NSO peaks (as seen in Figure 7 in [6]). However, based on Figure 4b, the K * s in a 3-D graph is illustrated in Figure 5b. The points along the purple line represent a TL of 50-km obtained for different k (π/km), in which the K * s varies from 1.20 up to 2.50. As k increases, the optimum K * s also increases which reflects on the damping resistances R s1 and NSO peaks which are more reduced. In an effort to find an analytical representation of these abacuses, using the Toolbox Curve Fitting in MatLab™, fitted equations for each 3-D surface K * s and K * d in Figure 5 are generated, given by a general polynomial form as follows: where P ij are coefficients for the optimum adjustable factors K * s and K * d . The P ij for K * d and K * s are shown in Tables 1 and 2, respectively. In both fitted polynomial equations of this work, N = 4 and M = 3.

Interface with ATP Software
As depicted in the Figure 6, an alternative to simulate the modified topologies of the LPM consist of employing the programing language MatLab combined with ATPDraw. In this case, using a programing code in MatLab, the per-unit-length line parameters (R', L', C'and G') are calculated for a single-phase line. Then, using the steps described in the flow-chart from Figure 5, the optimized topology using the damping adjustable factor (Kd* or Ks*), the number of π-circuits (n) and the discretization factor (k) are selected. An external code is applied to save these parameters (R', L', C', G', n, k, Kd* and Ks*), that will represent a line with a file in the format *.lib, as described in [21]. Finally, this *.lib file is loaded in the LIB icon in ATPDraw which will generate the modified topologies

Numerical Results and Discussions
The numerical results are divided into three sections organized as follows: in Section 5.1, the frequency responses of the modified topologies are compared to those obtained by the traditional LPM. In Section 5.2, the transient responses are computed for three scenarios in single-phase TLs. In Section 5.3, transient responses are computed for two scenarios using three-phase TLs. Due to the low-frequency range of the switching maneuver, the line parameters are considered constants for this study. The numerical method of Heun is adopted to solve all the state equations of each topology.

Frequency Responses for the TL Models
The single-phase TL is employed to calculate the frequency response of the transient voltages V m (s) at the open receiving end. For these simulations, a 1-p.u. DC voltage source is applied at t = 0 s at the sending end and the transient voltages are obtained with the traditional LPM and modified topologies T 1 and T 2 . Figure 7a shows a black-box model for transfer function H(s) represented for each previous topology. This transfer function is given by: where V k (s) and V m (s) are voltages at the sending and receiving ends in the frequency domain and s = jω is the complex angular frequency (rad/s). To clarify the results, the abscissa is normalized by the fundamental f 0 and harmonic f h frequencies expressed by [6,22]: where τ is the traveling time, L' and C' are the p.u.l. inductance and capacitance and d is the line length. The frequency response is calculated for the TL length of 20-km depicted in Figure 7a. This line is represented by the traditional LPM with k = 3 π/km and p.u.l. parameters of R = 0.075 Ω/km, L = 1.722 mH/km, C = 9.638 nF/km and G neglected. The TL represented by the modified topologies T 1 and T 2 , using fitted Equations (15) and 3-D abacus in Figure 5, present optimum adjustable factors of K * d = 0.526 and K * s = 1.943. The computed fundamental frequency f 0 is approximately 3068 Hz. The receiving-end voltage V m (s) is computed by a step function V k (s) = 1/s applied at the sending end. The computed V m (s) as a function of the normalized h = f h / f 0 is illustrated in Figure 7b where f h is the harmonic frequency given in (17), for m integer. The transient peaks occur at the odd harmonic frequencies (h) such as: 1 f 0 ≈ 3068 Hz, f 3 ≈ 9204 Hz, f 5 ≈ 15,304 Hz and so on. It can be observed that the voltage peaks for T 1 and T 2 topologies present lower values in comparison with those obtained by the LPM as the harmonic frequency increases. As a consequence, the peaks of the transient responses will have a pronounced reduction in the NSO when these modified topologies are used (further presented). The harmonic peaks of the H(s) are lower than those computed with the LPM and BM, for a given harmonic frequency, which represents an expressive attenuation in the time-domain responses.

Performance of the Modified Topologies
To investigate the performance in the transient responses with optimum adjustable factors in the modified topologies, single-phase TLs are represented by the traditional LPM and modified topologies T 1 and T 2 in three distinct scenarios, as depicted in Figure 8. The transient responses are compared with those obtained with the classical BM [18,23]. The TL parameters are: R = 0.075 Ω/km, L = 1.722 mH/km, C = 9.638 nF/km and G is neglected. The three scenarios studied are: In each scenario, a given line length d and discretization factor k are chosen. Based on Equation (15) and on Figure 5, the optimum adjustable factors are calculated and these values are organized in Table 3. Other values of K d and K s are employed to compute the transient responses which are described in each correspondent figure. The simulated responses obtained in the scenarios 1, 2 and 3 are illustrated in Figures 9-11, respectively. As seen in these figures, the NSO are significantly mitigated when the modified topologies are employed; however, high distortions in the rise and fall times are associated for each topology. In comparison with the responses for random factors K d and K s , the responses with the maximum attenuation and minimum distortion are obtained for optimum adjustable factors K * d and K * s in each of the modified topologies (T 1 ) and (T 2 ). The transient responses obtained with the K * d and K * s are in good agreement with those responses obtained by BM, assumed here as a reference. To verify the performance of the modified topologies, the percentage error (%) obtained with the LPM and with the modified topologies T 1 and T 2 are computed for each scenario. The percentage error (%) is given by where R represents the transient response (voltage or current) and the index mod − top is the transient response for each modified topology. The peaks P (p.u.) of each transient response and (%) obtained in the scenarios 1, 2 and 3, are organized in Table 4. It can be seen in this table, the response peaks computed with topologies T 1 and T 2 using the optimum K * d and K * s have provided the highest reduction in the NSO, resulting in the lowest percentage errors, as highlighted in light-green. These errors are lower than 1% in comparison with those calculated with the LPM, which gave errors of 3.20%, 11% and 26% as noted in Table 4. These numerical results confirm that the best agreement is obtained with the optimum adjustable factors K * d and K * s for the topologies T 1 and T 2 , respectively. In the next section, transient voltages are analyzed for the three-phase TLs under two different scenarios.

Modified Topologies in Three-Phase Lines
To extend the validation of the modified topologies for transients, the transient responses for a 50-km and 80-km three-phase TL are carried out, where modal decomposition technique is applied in each TL [24]. In this technique, a transformation matrix decouples a three-phase transmission line into three distinct single-phase lines (or propagation modes) where each one is represented by the traditional LPM and modified T 1 and T 2 topologies [25]. In this work, the Clarke's matrix is employed to decouple the three-phase TLs. Once the transient responses are obtained for each independent mode, they are converted to a phase domain employing inverse Clarke's matrix transformation. The ideally transposed three-phase TL is composed of bundles of 0.6-m quad-conductor per phase (phases 1, 2 and 3) plus two ground-wire conductors are installed. The spacing, including average conductor sag relative to the center of the tower and tower dimensions, are shown in Figure 12. The soil resistivity is ρ e = 100 Ωm, relative permittivity ε r = 1, and the ground-return impedances are computed by the Carson's approach. The longitudinal impedance and transversal admittance matrices ([Z] (Ω/km) and [Y], (µS/km), respectively), computed at 60 Hz, are given as follows: The three-phase TL configuration is provided in Figure 13, where two different scenarios are studied: (a) Scenario 4: The 50-km three-phase TL is energized by a 1-p.u./60-Hz voltage source only in phase 1 and the other phases are short-circuited. The receiving end is left open, as depicted in Figure 13a; (b) Scenario 5: The 80-km three-phase line is energized by a 1-p.u./60-Hz generator at t = 0 s and a balanced load is connected at the receiving end (Z load = 3 kΩ). At t = 33.33 ms, a fault occurs at phase 1, as depicted in Figure 13b.
In both scenarios, the transient voltages V m (t) are computed by the traditional LPM and modified topologies T 1 , T 2 with the best adjustable factors K * d and K * s and compared with BM. In Scenario 4, the 50-km three-phase TL is decoupled in three single-phase lines represented by k = 2π/km where the adjustable factors are K * d = 0.929 and K * s = 1.375, based on (15). In Scenario 5, the 80-km three-phase TL is decoupled in three singlephase lines represented by k = 1π/km where the adjustable factors are K * d = 1.51 and K * s = 1.236. To verify the mitigation of the NSO, simulation results are compared with transient responses computed by BM. The transient voltages V m (t) for Scenario 4 obtained with T 1 are shown in Figure 14 and those obtained with T 2 are shown in Figure 15.    Table 5. In Figures 14 and 15, the energization generates surge waves that travel along the line, resulting in overvoltages at the receiving end. Due to the inductive and capacitive couplings between the phases, induced voltages are produced at phases 2 and 3. It can be seen that when a positive variation in the voltage at phase 1 induces negative voltages at phases 2 and 3, as predicted by Faraday's law. In the steady state, the voltage at phase 1 will reach 1 p.u. while at phases 2 and 3 it is zero, due to the open-circuit condition. In these figures, one observes that the amplitude of NSO is significantly reduced when the three-phase line is represented by the modified topologies T 1 and T 2 . The voltage peaks and errors obtained are shown in Table 5. As noted, for the voltages in phase 1, the errors for the topologies T 1 and T 2 are equal to 0.82%, whereas the LPM provides an error of 12.64%. In relation to phases 2 and 3, the voltages with the LPM have presented an error of 26.82%, whereas the modified topologies T 1 and T 2 have provided errors smaller than 1.50%. These results indicate the good performance in the attenuation of NSO, especially, for the induced voltages at phases 2 and 3. In scenario 5, the transient voltages are generated by the three-phase generator switched on at t = 0 s, and at t = 33 ms an earth-fault occurs at phase 1 at the receiving end (the fault impedance is 5 Ω). The transient voltages at the receiving end V m (t) computed with BM, traditional LPM and modified topologies T 1 and T 2 are illustrated in Figure 16. In order to highlight the performance of the optimum adjustable factors K d and K s , the simulations are detailed as depicted in Figure 17 for T 1 and Figure 18 for T 2 , respectively. The transient voltage peaks in Figures 17 and 18, after the earth-fault at t = 33 ms, and the (%) obtained with the LPM, with the modified topologies T 1 and T 2 for best adjustable factors, and with BM for the scenario 5 are shown in Table 6 (only for phases 2 and 3). Due to the mutual coupling inherent in the three-phase line, the overcurrent in phase 1 produces the transients at the load, where the NSO are presented in the LPM responses. The voltage peaks are shown in Table 6. It can be seen in detailed parts of Figures 17 and 18 that voltage peaks are significantly mitigated in the transient responses when the modified topologies are employed. The errors in phase 2 for the modified topologies T 1 and T 2 are smaller than 1.50%, whereas the LPM has provided 12.78%. For phase 3, the errors are smaller than 1% and an error of 10.45% is seen for the LPM.   In general, the modified topologies behave as low-pass filter circuits that mitigate the NSO in the transient responses directly in the time domain and are in good agreement with those transient responses obtained with BM. Furthermore, no digital or analog filters are needed to mitigate the NSO in the computed transient responses when the simulations are carried out with modified topologies. The proposed algorithms are simple tools that can be implemented, and 3-D surfaces with fitted equations are a fast way to calculate the best adjustable factors of each topology. Based on these best adjustable factors, the maximum attenuation of the NSO and minimum distortion in the time-domain transient responses are obtained when low-frequency disturbances occur in TLs. Additionally, these modified topologies are easily incorporated into EMP-type programs.

Conclusions
This paper has presented two modified LPM topologies based on the cascaded π-circuit for the representation of short transmission lines subjected to energization maneuvers and faults, both being low-frequency disturbances common in power systems. Furthermore, to determine the best adjustable factors for the modified topologies, a simple programing code was described in this paper.
The modified topologies have presented interesting behavior regarding the mitigation of spurious oscillations in the simulated transient responses. It was demonstrated that the optimized adjustable factors provide the maximum mitigation in the NSO, especially the peaks, of the obtained transient responses. The adjustable factors K d and K s of these modified topologies were randomly selected by the authors in the literature. The optimization method proposed here selects the best values of the adjustable factors K * d and K * s for each modified topology. Each adjustable factor is a function of the line length and the number of π-circuit/km employed in the line representation. Using the Toolbox Curve Fitting in MatLab™(MathWorks, Natick, United States), the Matlab fitting tool, the optimum K * d and K * s were plotted in 3-D graphs where fitted equations are provided. Results have indicated expressive attenuation in the erroneous peaks of the NSO in the transient responses developed for single-and three-phase transmission lines when compared to those responses computed by the traditional LPM. Furthermore, the best K * d and K * s provide low distortion in the rise and fall times of the transient responses. The traditional LPM has produced errors around 25% in all simulations concerning singleand three-phase transmission lines. However, the percentage errors obtained with the modified topologies are very small in comparison to the reference responses obtained with Bergeron's line model. Results also demonstrated that the modified topologies produced transient responses which are in good agreement with Bergeron's line model for all transient simulations. As advantages, these modified topologies are expressed by space-state equations in which the solutions can be calculated by any numerical integration method. Further, these structures can be inserted into EMT-type programs. These modified topologies with the best adjustable factors can be applied in short transmission lines for low-frequency disturbances (energization maneuver and faults). The mitigation procedure is performed directly in the time domain and neither digital nor analog filters are required to reduce the NSO in the transient responses. For future works, other sophisticated optimization methods can be applied to this problem. This paper has highlighted the importance of proper transmission line representation where the erroneous peaks may affect the project at the insulation level of equipment, such as, insulation strings, preinsertion resistors, circuit breakers, and impact on the operation of the power system.