Optimal Electric Arc Furnace Model’s Characteristics Using Genetic Algorithm and Particle Swarm Optimization and Comparison of Various Optimal Characteristics in DIgSILENT and EMTP-RV

-e stochastic and nonlinear characteristics of electric arc furnaces (EAFs) lead to power quality challenges in the power system. In studying EAF behaviors, having optimized characteristics/models, selecting a suitable and optimum model that adapts to the actual characteristics of EAFs, and investigating simulation software’s capability for implementing EAF models are essential. However, the literature shows a research gap in investigating EAF simulations in various software products based on different models. -is paper studies several time-domain models, such as piece-wise linear, modified piece-wise linear, hyperbolic, exponential, and exponential-hyperbolic models, for EAF modeling and simulation. -e optimal estimation of parameters for the introduced models is necessary to adapt actual EAF characteristics. -us, one of the studies taken in this paper is optimizing the EAF model’s characteristics. -e proposed optimization problem is solved using the genetic algorithm (GA) and particle swarm optimization (PSO). Moreover, the optimized models are simulated in DIgSILENT and EMTP-RV to investigate different EAF models from the viewpoint of accuracy and efficiency. -e optimization of different EAF models’ characteristics and comparison of EMTP-RV and DIgSILENT in simulating EAF behavior are the contributions of this paper. -e proposed method is validated based on the actual data of a realistic EAF-based steel company in Iran. -e obtained results show that the modified piece-wise linear model has the most accuracy in identifying the EAF behavior. -e test results based on DIgSILENT and EMTP-RV simulations imply that the EAF could be simulated with high accuracy using modified piece-wise linear and piece-wise linear models. In general, EMTP-RV has expressed more accuracy in simulating different EAF models, and the simulation execution speed of EMTP-RV is around 2.5 times faster than DIgSILENT. In contrast, DIgSILENT is more suitable to facilitate the power system studies of EAF according to its extensive study tools and library.


Motivation and Incitement.
e deployment of electric arc furnaces (EAFs) is one of the well-known and efficient approaches in steel industries [1,2]. e EAF-based technologies have ranked as the second steelmaking process [3]. Due to the stochastic changes in the arc length and other behaviors of electric arcs during the fusion and refining process, EAF's characteristics would be dynamic, nonlinear, and time-variant [4,5]. Accordingly, EAF causes power quality problems in the supply network, such as voltage and current imbalances [6], harmonics [7,8], and voltage flicker [9][10][11]. To investigate the EAF's adverse effects in the point of common coupling (PCC) and provide appropriate solutions to deal with the disturbances, using precise EAF mathematical models and implementing the EAF in specialized simulation software is essential [12][13][14]. Although several studies have been done in the field of EAFs, there is a research gap in selecting an appropriate model incorporating a suitable simulation environment. is article tries to fill such a research gap by studying various EAF models, using the optimized estimated parameters in EMTP-RV and DIgSILENT.

Literature Review.
Several studies have proposed various models to identify the EAF behavior and characteristics [15,16]. Models used to describe EAF have been classified into two categories, e.g., frequency-domain and time-domain models. In the frequency-domain models, the EAF is modeled using a harmonic current source or a harmonic voltage source [17]. On the other hand, in time-domain models, the EAF models are presented based on the nonlinear differential equations of the furnace [18,19] or the nonlinear approximation of the EAF's voltage-current characteristic [4,20]. Also, in [21], the EAF dynamic models have been proposed, considering the stochastic variations in arc length.
In [18], the EAF has been modeled using a cascade combination of variable resistance and inductance in parallel with a current source. e variable resistance and inductance have modeled the EAF dynamic and time-varying behavior. Also, the current source has modeled EAF current harmonics. Reference [19] reported a current-controlled voltage source to model the EAF and calculate the voltage flicker.
Based on the voltage-current characteristic of an actual EAF presented in [22], references such as [23,24] have presented mathematical models approximating the EAF voltage-current characteristic. e approximated EAF models' characteristics have been used in the subsequent investigations for EAF studies. In [25], EAF voltage-current characteristics have been used for flicker studies. EAF power quality has been studied in [26] using voltage-current characteristics. e EAF voltage-current characteristics in [23] have been followed from the flicker studies viewpoint, considering the EAF transformers and cable system. In other articles, such as [27], the EAF voltage-current characteristics have been used for power quality investigations. In [28], the exponential-hyperbolic model has been proposed, and the model on the power system was studied using the PSCAD. Voltage harmonics and other power quality issues have been studied in [29,30], using EAF voltage-current characteristics.
Practical approaches are associated with the conformity of EAF models to the actual characteristics. is accuracy and conformity can be provided by estimating parameters for EAF models. Reference [31] has studied parameter identification of EAF models, which adapted arc measured values to estimate parameters of the EAF models. Despite the importance of representing EAF with an accurate model, parameter estimation and optimization of the EAF model's characteristics have not received much attention in the literature. e frequency-domain models have been mainly used for harmonic analysis in the EAF upstream networks [32]. On the other hand, time-domain models are suitable for studying EAF in power systems. us, several time-domain models derived by approximating the EAF voltage-current characteristic, such as piece-wise linear, modified piece-wise linear, hyperbolic, exponential, and exponential-hyperbolic models, have been introduced in the literature for studying the EAF behaviors [22,33]. e literature review in the field of EAF modeling and simulation is summarized in Table 1. As seen, less attention has been paid to optimizing EAF's various models' characteristics in the literature. Moreover, according to Table 1, in most recent studies, such as [19,32,34], MATLAB has been used for EAF modeling. Despite the many capabilities of MATLAB, this software does not have the conditions and requirements for specialized studies in the power system fields. Several references have studied the EAF models in PSCAD and EMTP-RV. PSCAD is a power system simulation software that has been utilized for EAF simulation in some research works, such as [4,35]. Also, EMTP-RV is a technically advanced software for transient analysis in power systems. In more recent papers, like [4,36], EMTP-RV has been used for EAF simulation studies. Moreover, DIgSI-LENT is a leading power system software that can be used for simulating EAF models. However, less attention has been paid to DIgSILENT to simulate the EAF and corresponding studies. e literature shows a research gap in investigating EAF simulations in various software programs based on different models with optimized estimated model parameters.

Contributions and Paper Organization.
As discussed, different time-domain static models have been developed for identifying the EAF behaviors in the power system. e optimal estimation of parameters for the EAF models is necessary to adapt actual EAF characteristics. One of the studies taken in this paper is optimizing the EAF model's characteristics by estimating the models' parameters. e proposed optimization problem is solved using the GA and PSO. e appropriate and effective simulation environment/ software for implementing the EAF models is another essential aspect of EAF's accurate simulation. In recent works, EAF models have been studied in general simulation environments like MATLAB rather than specialized software for power system simulations. To fill this knowledge gap, this research implements various optimized models' characteristics of EAF in leading power system software, such as EMTP-RV and DIgSILENT. Moreover, EMTP-RV and DIgSILENT performance in simulating EAF models have been compared in terms of accuracy and speed. e most important contributions of this paper in the research field of EAF modeling and simulation are listed as follows: Parameter estimation of EAF different models to adapt measured arc characteristics through optimization Optimizing the EAF models' characteristics, using the GA and PSO Implementation of different models for the EAF in EMTP-RV and DIgSILENT Comparing EMTP-RV and DIgSILENTperformance in simulating EAF models, considering accuracy and speed aspects Applying the proposed studies to an actual EAF test system located in Iran e rest of this paper is organized as follows. In Section 2, the problem statement is drawn. In Section 3, various EAF static models in the time-domain are introduced. In Section 4, optimization of EAF models' characteristics using GA and PSO is discussed. Section 5 presents the implementation and simulation of EAF models. In Section 6, the optimization and simulation results are proposed. Finally, the conclusion is drawn in Section 7.

Problem Statement and Research Overview
e nonlinear and stochastic nature of EAF during its operation causes power quality issues, such as voltage and current imbalances, harmonics, flicker, and other disturbances. To investigate the EAF's abnormal effects and provide appropriate solutions to deal with them, modeling and simulation of EAF characteristics are essential. us, precise mathematical models and specialized simulation software that can simulate the stochastic, time-variant, and nonlinear behavior of EAF are primarily required in EAF studies.
Selecting an appropriate model that adapts to the actual characteristics of EAFs is an essential issue. In addition, optimizing the estimated parameters of each model based on the actual measured values of the voltage and current of the EAF is inevitable. Moreover, investigating simulation software's capability for implementing EAF models is another concern that received less attention in the literature. In this paper, the EMTP-RV and DIgSILENT simulations of EAF based on the optimized estimation of various models are examined in the following steps: Step 1: Estimation of EAF models' parameters, e.g., piece-wise linear, modified piece-wise linear, hyperbolic, exponential, and exponential-hyperbolic models, by the GA and PSO.
Step 2: Simulating the EAF based on various models with the optimized estimation of their parameters in EMTP-RV and DIgSILENT.
Step 3: Analyzing the simulation results of EMTP-RV and DIgSILENT based on various models and different aspects, like the conformity index.
Step 4: Selecting suitable models for EAF simulation in EMTP-RV and DIgSILENT.

Model 1: Piece-Wise Linear Model.
e piece-wise linear model is one of the simplest EAF models in the time-domain based on the linear approximation of the EAF's voltagecurrent characteristics. In this model, the relation between arc voltage and arc current is expressed by (1) and (2) [29,37].
In (1) and (2), v a and i a denote the arc voltage and arc current. In addition, v ig and v ex represent the arc ignition voltage and arc extinction voltage, respectively. Moreover, i 1 and i 2 are used to express the arc ignition current and the arc extinguishing current. As shown in Figure 1, R 1 and R 2 would be the slope of OA and AB segments in the voltagecurrent characteristic of the piece-wise linear model.

Model 2: Modified Piece-Wise Linear Model.
e modified piece-wise linear model is developed based on more precise approximations of the EAF's characteristics. e mathematical expression of the modified piece-wise linear model has been given in [29] v a � In (3), the expressions inc and dec indicate increasing and decreasing parts in the arc current, respectively. In addition, v st is a voltage between the arc ignition voltage (v ig ) and the arc extinction voltage (v ex ). Also, v st happens due to a sudden voltage drop across the EAF electrodes during the melting stage. is voltage causes the arc voltage to decrease exponentially from v ig to v st . Moreover, R 1 , R 2 , and R 3 are the slope of the OA, CB, and BD segments in the voltagecurrent characteristic of the modi ed piece-wise linear model, respectively. In (4), i 1 , i 2 , and i 3 , which are constants, are de ned. Figure 2 shows the voltage-current characteristic of the modi ed piece-wise linear model.

Model 3: Hyperbolic Model.
e hyperbolic model is one of the most widely used time-domain EAF models. In this model, the voltage rising time is neglected and causes sudden arc voltage changes during the zero-crossing arc current. e voltage arc in the hyperbolic model is formulated using [30] v a v at + C D + i a signum i a .
In (5), C andD are constant parameters related to the arc power and arc current. Also, v at is the voltage that depends on the arc length. e voltage-current characteristic of the hyperbolic model is shown in Figure 3.

Model 4: Exponential Model.
In the EAF's exponential model, as shown in (6), an exponential function is used to approximate the EAF's voltage-current characteristic [28].
In (6), the constant I 0 denotes the slope of positive and negative currents. Moreover, the voltage-current characteristic of EAF's exponential model is depicted in Figure 4.

Model 5: Hyperbolic-Exponential Model.
e hyperbolicexponential model is a combination of hyperbolic and exponential models for EAFs. e mathematical expressions and formulas of the hyperbolic-exponential model are given in [28] v a v at According to (7), the arc voltage during current rise is modeled with a hyperbolic function, and during current descent it is modeled with an exponential function. e voltage-current characteristic of the hyperbolic-exponential model is shown in Figure 5. denote the modeled arc voltage, the measured arc voltage, the modeled arc current, and the measured arc current, respectively. e proposed optimization problem aims to minimize the error between the voltage and current waveforms of the model and the measured voltage and current waveforms.

Optimization of EAF Models' Characteristics
e optimization problem of the EAF model's characteristics is a nonconstrained minimization problem. However, specifying the upper and lower boundaries for each model input parameter reduces the problem search space and leads to better convergence to the optimal global value. e parameter boundaries can be assigned using the de nition of parameters in the EAF model and based on the results of previous studies.

Proposed Optimization Problem
Solving. According to the nonlinear nature of the proposed optimization problem, evolutionary algorithms are more appropriate to solve it. PSO and GA are e cient optimization techniques for solving nonlinear problems in terms of speed, simpli ed implementation, robustness, and capability of converging to global optima. In this paper, the GA and PSO algorithms are utilized. e statistical analysis will be performed to evaluate the performance of both PSO and GA in solving the proposed optimization problem. e statistical analysis would be helpful for examining whether the optimization algorithms have overcome the local optimization traps around the feasible region or not. Moreover, the convergence and the minor deviation in reaching the global optima in various runs and di erent cases should be used to investigate the suitability of the algorithms for solving the proposed optimization problem. It should be noted that the proposed optimization problem can be solved by other optimization algorithms. However, the examination of the capability of the selected optimization algorithms to nd the global optima and other performance criteria is necessary. e owcharts of solving the proposed optimization problem to estimate the parameters of the EAF's models by the GA and PSO algorithms are shown in Figures 6 and 7, respectively. As illustrated in Figures 6 and 7, the network data, problem constraints, and EAF measured values should be imported to de ne the optimization problem. Network data include the resistance and inductance values of exible cables and the secondary voltage level of the EAF transformer. According to the EAF model, the allowed range for decision variables (model parameters) and the problem constraints are determined. Also, measured arc voltage and arc current in a cycle are imported to the proposed optimization problem as the base waveforms. Finally, the optimization process starts using one of the heuristic algorithms.

GA.
e GA is a numerical problem-solving technique that is inspired by Darwin's evolution theory. In GA, each chromosome is a potential solution to the problem, and the genes are related to the problem variables [38]. e GA initializes with a randomly generated population of chromosomes. In each iteration, OF value of each chromosome is evaluated. en, the selection operator selects some of the best chromosomes using a stochastic process [39]. e crossover and mutation operators are applied to the selected chromosomes, producing a new generation of chromosomes [40,41]. is process continues until a certain number of iterations or convergence criteria is reached. Eventually, the algorithm converges to the ttest individuals over the iterations.

PSO.
e PSO is a metaheuristic optimization algorithm developed based on the social behavior of birds ocking in search of food sources. In the PSO, particles are

Arc Voltage
Arc Current Figure 4: Voltage-current characteristic of EAF's exponential model [28]. 6 International Transactions on Electrical Energy Systems search agents in the feasible space, and the particles' position represents valid solutions to the problem. Particles move to new positions to discover better positions by tracking optimal ones [42]. In the beginning, particles are scattered randomly in the feasible region. As shown in (9), x i (k) and v i (k) are respectively the position and velocity of the ith particle in the kth iteration. In each iteration, the OF value corresponding to each particle is evaluated, and the global best position denoted as g best and the individual best position denoted as P i best are updated. e velocity of particles is updated using (9), and they are transferred to the new positions using (10).
Measuring arc voltage and current waveforms in an actual steel company using Fluke power quality analyzer .

Start
The under study network data including upstream network's specifications, transformers, lines, etc.

Determine the EAF model
Determine the allowable range of decision variables in the proposed optimization problem Input network data Import the voltage and current waveforms Selecting one of the time-domain models: Model (1) to (5).
Using the definition of parameters in the EAF model and based on the results of previous studies.
Random production of initial population based on allowable range determined for optimization problem decision variables.

Generate initial population
Evaluate the objective function for each particle  Measuring arc voltage and current waveforms from the understudy test system using Fluke power quality analyzer.

Start
The understudy network data, including upstream network's specifications, transformers, lines, etc.

Determine the EAF model
Determine the allowable range of decision variables in the proposed optimization problem Input network data Import the voltage and current waveforms Selecting one of the time-domain models: Model (1) to (5).
Using the definition of parameters in the EAF model and based on the literature.
Generating random initial population of answers as chromosomes based on problem constraints.

Generate initial population
Evaluate the objective function for each chromosome Using the Roulette wheel selection operator Select superior chromosomes using selection operator Using bitwise mutation, uniform crossover, single point crossover and double point crossover.
Apply mutation and crossover operators on selected chromosomes Receiving the required data from power quality measurements, and defining the optimization problem in MATLAB Evaluate the software performance in simulating EAF models Evaluation criteria: conformity index and instantaneous error Implement the EAF model using coding in DSL environment Implement EAF element in the simulated supply system Simulate EAF using optimal model parameters and perform harmonic-based analysis EMTP-RV Implement the EAF model using block diagrams Implement EAF element in the simulated supply system Simulate EAF using optimal model parameters and perform harmonic-based analysis According to (9), r 1 and r 2 are random vectors in the range of [0 1]. Also, c 1 and c 2 indicate the personal learning coefficient and global learning coefficient, respectively. In (9), ω represents the inertia weight which decreases during iterations using a Damp coefficient through is process continues until a certain number of iterations or convergence criteria is reached. Finally, particles would converge to the optimal position over successive movements [43,44].

EAF Simulation in DIgSILENT and EMTP-RV
According to the mathematical models of the EAF, the arc voltage at any time step is calculated based on the arc current at that moment. us, in EAF models, a new relationship is required to calculate the arc current at any time. According to the EAF power supply circuit illustrated in Figure 8, using the KVL relation from the secondary side of the EAF transformer to the electrode ends, (12) is extracted. In (12), L eq and R eq are the equivalent resistance and inductance of the network.
Since the EAF models in DIgSILENT and EMTP-RV are implemented discretely, converting (12) to the discrete relation in (13) is necessary.
In (13), Δt and T denote the duration of simulation steps and the simulation step number. By arranging (13), the arc current is calculated at any moment based on the arc current and arc voltage at the previous step as given in Δt i a (T) + R eq i a (T) .
Using (14) and assuming an initial value for arc current, the EAF models can be implemented in different software. e EAF mathematical models are simulated in DIgSILENT using the DSL environment and in EMTP-RV using control block diagrams.
In DIgSILENT and EMTP-RV libraries, no power element simultaneously controls its voltage and current. us, in the EAF element structure, a current source in series with a variable resistor is used, as shown in Figure 9.
As seen in Figure 9, if the arc current waveform (i a ) is applied to the current source and the arc resistance waveform (v a /i a ) is applied to the variable resistor, the EAF voltagecurrent characteristic is created at the PCC (Bus (A)).
Indeed, the voltage of Bus (A) and the current passing through this bus are precisely equal to the EAF voltage and current. e current of Bus (A), which is the same as the arc's current, causes the arc's voltage in Bus (A) to pass through the arc resistance. e EAF implementation using the mentioned elements in DIgSILENT and EMTP-RV is shown in Figure 10.

Test System Descriptions and Measurement Results.
e proposed method is tested on an existing EAF plant in Iran. e single-line diagram of the understudy EAF's supply system is shown in Figure 11. As depicted, the studied EAF has been connected to the external grid through two sets of transformers, including the 400/63 kV power transformer and the EAF transformer (63/0.9 kV) [15]. A 63 kV cable system with a nominal current of 1100 A with 1320 m length has been installed for the power transmission section. e EAF electrodes are connected to the supply system via flexible cables. e electrodes and flexible cables are modeled through l eq and r eq impedances.
e EAF voltage, current waveforms, and harmonic distortion have been measured/recorded to validate the proposed study. e Fluke 437-II Power quality monitor and energy analyzer have been used in this study. In Figure 12, the EAF voltage and current measurements are shown. Also, Figure 13 presents the recorded harmonic distortion of the studied EAF [8].
e measured waveforms are used as the reference waveforms of voltage and current in the EAF optimization problem. In Figure 14, the measured values and characteristics of understudy EAF are shown.

Optimization Results.
e measurement values are imported into the optimization problem. Also, the GA and PSO setting parameters are selected, as given in Tables 2 and  3, for solving the proposed optimization problem. e convergence diagrams of GA and PSO in solving the optimization problem for different models are shown in Figure 15.
As revealed by test results shown in Figure 15, the convergence of GA is much higher than PSO. In addition, it can be concluded that as the number of decision variables in the optimization problem of the EAF model characteristics increases, the PSO convergence speed decreases. e statistical analysis of solving the proposed optimization problem by the GA and PSO in optimizing the EAF model characteristics is presented in Figure 16 and Table 4. e GA would be suitable for solving discrete problems, while the PSO algorithm has a better performance in solving continuous problems. Because the optimization problem of the EAF model's characteristics is continuous, it is expected that the PSO results in better performance in finding the optimal points. Statistical analysis shows that although the PSO has a relatively slow convergence rate than the GA, it is more robust in finding the global optimal point. For instance, in the statistical analysis diagram of the hyperbolic model, although the obtained optimal values are equal to each other with one-tenth of accuracy, the GA has not been able to find the global optimal point with one hundred-thousandths of accuracy. 8 International

International Transactions on Electrical Energy Systems
Since all models include a few decision variables, both algorithms reach the global optima in almost all cases. erefore, the standard deviations in optimum OFs and solutions would not be signi cant in such small-sized problems. However, the smaller the standard deviation for an algorithm in solving an optimization problem, the algorithm is more reliable in nding the global optima.
According to the statistical test results and the convergence diagrams in di erent cases, it is possible to ensure that the GA and PSO techniques have suitable performance in solving the proposed problem. It can be concluded that the converged solutions are global optimums. Since the presented research focuses on optimizing di erent EAF models' characteristics rather than the optimization method, it has been con ned to introducing GA and PSO to solve the optimization problem. However, solving the optimization problem of EAF models' characteristics is not restricted to the presented techniques, and other optimization algorithms can be used to solve it. e optimized parameters for EAF models' characteristics using the GA and PSO are given in Table 5. Moreover,  the optimum voltage-current characteristics of di erent EAF models are shown in Figure 17. e harmonics content of voltage and current in EAF models has been compared with measured values in Tables 6  and 7. e OF value for di erent models of the EAF comparing the harmonics contents with the measured values shows that model 2 (modi ed piece-wise linear) has the atmost accuracy in modeling the EAF behavior.
According to the problem conditions and the intended studies, other models can be used for EAF simulation. For

EMTP-RV and DIgSILENT Simulation
Results. Models 1 to 5 are implemented in DIgSILENT and EMTP-RV using DSL and control circuits in this section. It is worth noting that simulated characteristics under di erent models in both software are accurately similar to the optimized models' characteristics. e EAF element should be simulated with the same accuracy as the mathematical model in the understudy system using e cient specialized simulation software. e critical point is that the simulated EAF must implement the voltage and current signals generated by the model with the same accuracy and quality in the power grid.
In simulation studies, the EAF element is implemented in the power grid indicated in Figure 11. e power system parameters and speci cations of the actual test system have been presented in Table 8. e EAF element is simulated utilizing various models in DIgSILENT and EMTP-RV. Figure 18 shows the voltage      us, the EAF current waveforms in both software are the same as the current waveforms in the respective model. erefore, the arc voltage waveform of the EAF element, created by passing the arc current through the variable resistance at the furnace bus, is compared with the arc voltage waveform of the respective model. However, a large number of data samples contained in the measured, modeled, and simulated waveforms of the arc voltage complicate this comparison. erefore, the conformity index and instantaneous error have been used to quantitatively compare the performance of DIgSILENT and EMTP-RV in EAF simulations based on different models. e conformity index for two parameters of the same material is defined according to the least square method. e OF of the optimization problem of the EAF model's characteristics is the sum of the conformance indices between the measured values and the model of arc voltage and arc current in one cycle. For comparing the performance of DIgSILENT and EMTP-RV in EAF simulations based on different models with the optimized estimation of their parameters, the conformity index between the measured values, the mathematical model, and the simulated arc voltage waveform per one cycle is given in Table 9.
As the conformity index considers the conformance of the total cycle in the evaluation, the instantaneous error represents another criterion that observes the conformance of the arc voltage waveform sample by sample. e instantaneous error is helpful in detecting the significant deviations from the reference waveforms. For comparing the performance of DIgSILENT and EMTP-RV, the maximum and minimum instantaneous error between the measured values, the mathematical model, and the simulated arc voltage waveform is given in Table 10.
Comparing the conformity index between the model values and the simulated values by DIgSILENT, the DIgSILENT has simulated only models 1 and 2 with high accuracy. e significant increase in conformity index for model 5 is due to the instantaneous peaks generated in DIgSILENT, which can be explained using the instantaneous error. According to Table 10, the maximum instantaneous error between the model values and the simulated values by DIgSILENT under model 5 is greater than in other models. e conformity index and instantaneous error between the model values and the simulated     Due to the large number of data samples in the understudy waveforms, the instantaneous error could be better when considering the instantaneous peaks generated in both software programs in comparison. e instantaneous error between the model values and the simulated values demonstrates that the instantaneous peaks under simulated models 1, 2, and 4 are less than simulated models 3 and 5.
According to the nature of models 3 and 5, in simulation with simulation intervals (0.04 μs), these models undergo rapid changes at the moment of passing zero. e rapid changes cause high-frequency oscillations in the DIgSI-LENT; as in model 5, the oscillation peak increases to 12 kA. In comparison, EMTP-RV is more accurate and can control these rapid changes in the power system. e comparison of simulation execution time shows that the speed of EMTP-RV is 2.5 times faster than DIgSILENT. Also, unlike EMTP-RV, which has been often used for transient studies, DIgSILENT is a suitable platform for implementing power systems according to its powerful user interface, extensive study tools, and comprehensive industrial library. erefore, to facilitate the power system studies of EAF test systems and steel company electrical grids, it is preferred to implement the EAF load in the DIgSILENT. Since model 2 precisely expresses the EAF behavior and has been implemented with equal accuracy in DIgSILENT and EMTP-RV, it is concluded that model 2 can be suggested as the best alternative for implementing the EAF-based test systems in DIgSILENT.

Conclusion
Much attention has been paid to EAF modeling in the literature. Different mathematical time-domain models have been developed to describe the EAF's time-variant and nonlinear behavior. In this paper, the parameters and characteristics of various time-domain models, such as piece-wise linear, modified piece-wise linear, hyperbolic, exponential, and hybrid models, have been optimized using GA and PSO. e actual measurements of a realistic EAF in Iran have been used to validate the obtained optimum characteristics and models. It has been concluded that the modified piece-wise model is adequately efficient in providing suitable accuracy with the actual EAF characteristics compared to other models. is paper has also addressed the simulation and implementation of different mathematical models in specialized power system software. Investigating the simulation results of DIgSILENT and EMT-RV for various EAF models is another contribution of this research. e test results extracted from DIgSILENT and EMTP-RV indicated that the EAF could be simulated with high accuracy using modified piece-wise linear and piece-wise linear models. However, the simulation of other EAF models in DIgSI-LENT might not be adequately precise. In contrast, the simulation results in EMTP-RV by different EAF models have been satisfying. Moreover, the simulation execution speed of EMTP-RV is around 2.5 times faster than DIgSI-LENT. Test results and comparative analysis illustrate this research's importance in optimizing the EAF models and investigating DIgSILENT and EMTP-RV simulations by various models.

Nomenclature
Parameters i a : Electric arc furnace (EAF)/arc current v a : EAF/arc voltage v ig : Arc ignition voltage v ex : Arc extinguishing voltage i 1 : Arc ignition current i 2 : Arc extinguishing current R 1 , R 2 , R 3 : Slops in voltage-current characteristics of EAF in various models v st : Arc voltage between the arc ignition voltage and the arc extinction voltage in the modified piecewise linear model i 3 : Arc current at v st in the modified piece-wise linear model C: Constant parameter related to the arc power in EAF's hyperbolic and hyperbolic-exponential models D: Constant parameter related to the arc current in EAF's hyperbolic and hyperbolic-exponential models e voltage, depending on the arc length in EAF's hyperbolic, exponential, and hyperbolicexponential models I 0 : e slope of positive and negative currents in EAF's exponential and hyperbolic-exponential models v model a : Modeled arc voltage v actual a : Measured arc voltage i model a : Modeled arc current i actual a : Measured arc current S sc : e short circuit power at any bus of the power system OF: Objective function V: e evenin voltage of the understudy point (EAF's connection point) X th : e evenin reactance of the understudy point (EAF's connection point) V nom : Nominal voltage S nom : Nominal apparent power %UK: Rated short circuit voltage/series reactance of the EAF's transformer R: Resistance X: Reactance L: Inductance THD V : Total voltage harmonic distortion THD I : Total current harmonic distortion v i (k): e velocity of the ith particle in the kth iteration of the PSO algorithm x i (k): e position of the ith particle in the kth iteration of the PSO algorithm g best : e global best position discovered in the PSO algorithm P i best : e best-experienced position of the ith particle in the PSO algorithm r 1 and r 2 : e random vectors in the range of [0 1] in the PSO algorithm ω: Inertia weight in the PSO algorithm Damp: Inertia weight damping rate in the PSO algorithm c 1 : Personal learning coefficient of particles in the PSO algorithm c 2 : Global learning coefficient of particles in the PSO algorithm Δt: Duration of simulation steps N: Total number of measured samples n: Index of measured sample T: Simulation step number L eq : e equivalent inductance of the upstream network R eq : e equivalent resistance of the upstream network l eq : e equivalent inductance of the flexible cables and EAF electrodes r eq : e equivalent resistance of the flexible cables and EAF electrodes R a : Arc resistance. Data Availability e data used to support the findings of the study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.