Next Article in Journal
TranStutter: A Convolution-Free Transformer-Based Deep Learning Method to Classify Stuttered Speech Using 2D Mel-Spectrogram Visualization and Attention-Based Feature Representation
Next Article in Special Issue
Smart Grid Outlier Detection Based on the Minorization–Maximization Algorithm
Previous Article in Journal
Integrating Target and Shadow Features for SAR Target Recognition
Previous Article in Special Issue
Decentralized Sensor Fault-Tolerant Control of DC Microgrids Using the Attracting Ellipsoid Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplifying Rogowski Coil Modeling: Simulation and Experimental Verification

Department of Electrical, Electronic and Information Engineering, Guglielmo Marconi, Alma Mater Studiorum, University of Bologna, Viale del Risorgimento 2, 40136 Bologna, Italy
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(19), 8032; https://doi.org/10.3390/s23198032
Submission received: 28 August 2023 / Revised: 21 September 2023 / Accepted: 21 September 2023 / Published: 22 September 2023
(This article belongs to the Special Issue Sensors and Fault Diagnostics in Power System)

Abstract

:
The integration of renewable energy sources, electric vehicles, and other electrical assets has introduced complexities in monitoring and controlling power networks. Consequently, numerous grid nodes have been equipped with sensors and complex measurement systems to enhance network observability. Additionally, real-time power network simulators have become crucial tools for predicting and estimating the behavior of electrical quantities at different network components, such as nodes, branches, and assets. In this paper, a new user-friendly model for Rogowski coils is presented and validated. The model’s simplicity stems from utilizing information solely from the Rogowski coil datasheet. By establishing the input/output relationship, the output of the Rogowski coil is obtained. The effectiveness and accuracy of the proposed model are tested using both simulations and commercially available Rogowski coils. The results confirm that the model is simple, accurate, and easily implementable in various simulation environments for a wide range of applications and purposes.

1. Introduction

In recent years, the network has undergone a significant revolution. On one hand, there have been notable changes in electric assets. For instance, photovoltaic plants are becoming increasingly common on the roofs of industrial, commercial, and private buildings [1,2]. The number of electric vehicles has rocketed in the last few years thanks to commercial maneuvers and country-based incentives [3,4]. Electric scooters have also become a popular means of private transportation, requiring a power supply for charging [5,6]. On the other hand, monitoring techniques and associated software have been evolving. The advancement and widespread adoption of artificial intelligence (AI) and machine learning (ML) techniques have brought about significant changes in power system monitoring techniques [7,8]. The common thread among all these developments is instrumentation. The new assets can be monitored through sensors and distributed measurement systems deployed throughout the power network. These sensors provide measurements for the monitoring system and the learning algorithms used in ML. Therefore, it can be concluded that sensors and instrumentation play a significant role.
Voltage and current measurements are the primary activities conducted at network nodes and branches. However, numerous other quantities are also measured, such as humidity, temperature, power, energy, pressure, etc. [9,10,11]. All these electrical and environmental parameters are measured to obtain the most accurate health status of the grid. As a result, studies and standards dedicated to sensors and instrumentation are published on a daily basis. Regarding standards, the IEC 61869 series comprises fifteen main documents that pertain to instrument transformers (ITs). Readers can refer to documents IEC 61869-1 and -6 [12,13] or the generic specifications applicable to ITs and low-power instrument transformers (LPITs), respectively. In recent literature, a calibration technique associated with partial discharge measurement using a high-frequency current transformer (HFCT) is described in [14]. ITs with digital output are calibrated in [15]. The authors in [16] characterize a measurement chain that includes LPIT and the phasor measurement unit (PMU). The aim in [17] is to accurately measure distorted signals, while [18] deals with design solutions for ITs, and [19] focuses on failure investigation techniques to be applied to ITs.
In this paper, the IT considered is the Rogowski coil. Due to its features, the Rogowski coil is widely adopted in various applications beyond power systems. Consequently, Rogowski coil modeling has been a research topic for the past few decades, and numerous works can be found in the literature. For example, an electro-thermal model written in VHDL-AMS language and based on finite element analysis (FEA) is described in [20]. In [21], the authors design a Rogowski coil from scratch for three-phase current measurement in electric motors. A highly detailed model is presented in [22], relying on equations that require the complete geometry of the device and the space in which it is installed. A four-layer PCB-based Rogowski coil is designed in [23]. A simple model, although based on several preliminary measurements, is detailed in [24]. The same authors introduced a new testing signal, known as the sinc response, to be measured by the Rogowski coil to obtain its input-output response [25]. The high-frequency behavior of the Rogowski coil is studied in [26] using the classical lumped model, while [27] explores another high-frequency modeling approach based on a black-box conceptualization of the Rogowski derived from preliminary measurements on the device.
From the literature review, it becomes apparent that there is a lack of very simple models that do not require extensive preliminary measurements. Furthermore, obtaining information on the internal structure and design of the device is not always feasible. This makes the modeling even more complicated and not generalizable for every commercial Rogowski coil. To address this issue, this paper introduces a new and straightforward method to obtain the frequency model of a commercial Rogowski coil. In practice, it is common to purchase a device for various applications rather than designing a custom Rogowski coil from scratch. The proposed model is based on a few inputs obtained from the device’s datasheets, and a set of equations is developed to derive a transfer function dependent on the cross-sectional shape. The lack of knowledge about commercial devices’ internal structure is an obstacle to many existing models. Subsequently, the obtained model is implemented in the Matlab 2023a environment for preliminary efficacy testing. Finally, the model is experimentally validated through measurements performed on commercial Rogowski coils. The results clearly demonstrate the validity and accuracy of the frequency model, confirming its potential for implementation in simulation environments such as real-time power network simulators.
The remainder of this paper is structured as follows: Section 2 includes the theoretical concepts and motivation behind the work. Section 3 is dedicated entirely to the model description and contextualization. Model validation is presented in Section 4. Finally, Section 5 concludes this paper and provides suggestions for future research.

2. Motivation and Background

2.1. The Rogowski Coil

The Rogowski coil (RC) is an LPIT for the measurement of AC and, with some expedients, also DC current. Manufacturers can refer to the standard IEC 61869-10 [28] for designing details, testing procedures, and accuracy evaluation. RC is widely adopted for power system applications due to its features. They are small, lightweight, and free from iron cores, and often they can be opened to facilitate their installation. Considering other current measurement solutions, e.g., inductive current transformers (CT), the RC features become significant.
The RC working principle is quite simple, and it can be described looking at Figure 1.
A copper wire is wound around an iron-free insulating support. This prevents any saturation behavior, typical of iron-core-based devices. The conductor, whose current I P must be measured, is inserted inside the RC. Its output, a voltage, u S t , is proportional to the derivative of I P :
u S t = M d I P d t
where M is the mutual inductance between the primary and secondary windings. The input-output relationship of the RC raises some comments. On the one hand, being the output voltage, it becomes easier to connect the RC to common analog-to-digital converters (ADC) without any other element in the measurement chain. On the other hand, u S is proportional to the derivative of the current and not to the current itself. Therefore, to obtain the correct value of I P , an analog or digital integration step is required. As for this paper, the integration system is not considered. The aim is to model the RC behavior independently of the integration solution adopted (it might introduce many aspects to be considered to obtain the input-output relationship of the RC).

2.2. Motivation

The need for RC models is supported by their massive adoption in several applications in different fields. As previously mentioned, many models are already available; however, their implementation is typically not straightforward. Therefore, the main goal of this paper is to provide an easy-to-use modeling procedure to obtain the input-output relationship of whatever commercial (or not) RC. The modeling procedure does not involve any preliminary electrical measurements on the RC but only external geometrical ones. Of course, the availability of preliminary information is an added value for the modeling procedure.
The benefits of a simplified modeling procedure are several. For example, an a priori estimate of the RC response would be possible without performing a time-consuming laboratory measurement. Second, the use of digital simulators is spreading; hence, having an accurate RC model would allow for the integration of realistic current measurements inside the simulations. Finally, the optimum exploitation of a model is the creation of a RC digital twin, which would help from the design to the implementation of the RC inside each application.

3. The Modelling Procedure

The proposed modeling procedure starts with the geometrical measurements on the RC and ends with the transfer function (TF) that correlates the frequency relation between the output voltage and the input current. The procedure is summarized in the flowchart depicted in Figure 2.
The first step is to measure the inner diameter 2a and the diameter of the cross-section d. The frequency modeling procedure is developed for three cross-section geometries: circular, square/rectangular, and oval. In the latter two geometries, an extra measurement is needed, as explained in Figure 3. As for the oval case, commercial devices do not present ideal oval shapes. Therefore, the geometry is often schematized, as shown in Figure 3.
The second step consists of extracting some information from the datasheet. The two key parameters needed for the modeling are the accuracy class and the transformation ratio (TR). Other typical parameters, given in the RC datasheet, are listed in Table 1. Note that the TF also exploits the information about the rated burden. However, if not provided by the manufacturer, the standard value of 2 MΩ given in [13] can be used.
The third step allows us to obtain the mutual induction, M starting from (1). Substituting the derivative over time with the angular speed ω = 2 π f :
M = T R 50 ω
where f is the generic frequency and T R 50 is the rated transformation ratio (at 50 Hz). Once M is found, it is possible to obtain the TR at each frequency of interest. For this assumption, the parasitic parameters will be neglected.
The fourth step regards geometry. Using the quantities obtained in the first step, the perimeter and the cross-section area of the RC can be calculated.
At this point, using consolidated expressions of the mutual inductance [29,30], the number of turns of the RC can be obtained from the reversed formula. Note that, from a practical perspective, the lack of knowledge of the number of turns is one of the limiting factors of RC modeling. Hence:
N r e c = 2 π M d μ l o g 2 a + 2 d 2 a
N c i r = 2 M μ 2 a + d 2 a + b a
N o v a = 2 π M μ 2 a + d 2 a + b a + h π log a + b a
where μ is the overall permeability, and N r e c , N c i r , and N o v a are the number of turns of the rectangular, circular, and oval cross-section, respectively. They can be used to obtain, in addition to the geometrical parameters, the length l w of the single coil, section A w , and radius r w of the wire used to wound the RC.
The last step is the calculation of the parasitic parameters of the RC:
C i = 2 ε π 2 2 a + d log A i π r w
L r e c = μ d l o g a + d a N r e c 2 2 π
L c i r = μ N c i r 2 2 a + d 2 a + d a 2
L o v a = μ N o v a 2 2 a + d 2 a + d a + h π log a + d a 2 π
where ε is the overall permittivity (to be fixed depending on the used material), C and L are the parasitic capacitance and inductance, respectively. The suffix i is used to switch among the cross-section parameters ( A i ). As for the winding resistance, R w , it can be easily obtained with the second Ohm’s law, measured, or extracted from the datasheet. Of course, there will be datasheets with extra information that can be used instead of some parameter expressions. However, the proposed method has been generalized to avoid any exceptional cases of missing (extra) inputs. A small note on permeability and permittivity is necessary. The general notation has been used to highlight that, depending on the RC adopted, the user may need to insert or not insert the relative permeability/permittivity. Further tests on these specific parameters demonstrated an almost negligible effect of the permittivity/permeability value on the accuracy parameters.
With the results from (2)–(9), the RC transfer function in the “s” domain can be obtained as:
T F = Z M s Z C L s 2 + s L + R w C Z + R w + Z
where Z is the rated burden, or better, the impedance of the measuring device connected in cascade to the RC. The TF in (10) must be customized with the parameters associated with the studied cross-section. To better contextualize the TF in (10), a schematic with the equivalent circuit is depicted in Figure 4. All the symbols in the picture were previously described.

4. Validation

In this section, two types of validations are presented. The first validation is obtained by simulation, and the second by experimental measurements on commercial RCs.

4.1. Validation by Simulation

To verify the performance of the proposed method, two sets of tests are performed. First, a frequency sweep test is simulated to assess the accuracy vs. frequency. Second, the MonteCarlo method (MCM) is applied to quantify the effect of the sources of uncertainty on the accuracy of the method.

4.1.1. Frequency Sweep

The simulation of the frequency sweep was run in the Matlab 2023 environment. A sinusoidal signal with an amplitude of 1000 A and a 0° initial phase was generated for all the frequencies from 50 Hz to 2500 Hz at 50 Hz steps. The frequency range was selected to cover the power-quality frequency range. As for the acquisition details, the input currents were simulated using a 200 ms window and a 50 kSa/s sampling frequency. Afterwards, each signal was multiplied by the TF in (10) to obtain the RC response. By means of the discrete Fourier transform (DFT), the goodness of the RC response was evaluated with the typical parameters: ratio error ε and phase displacement φ :
ε = T R f I S f I P f I P f 100
φ = I S f ^ I P f ^
where T R f , I S f , and I P f are the transformation ratio, the rms of the secondary current, and the rms of the primary current, respectively, evaluated at the generic frequency f . Analogously, I S f ^ and I P f ^ are the phase angles of the secondary and primary currents, respectively.
The results of this first set of tests are listed in Table 2 and depicted in Figure 5. As for the RC simulated, it is a rectangular cross-section, 0.5 accuracy class device.
From the results, it can be concluded that the estimation of the RC output is accurate in all the PQ frequency ranges. In terms of absolute value, the results are aligned with the accuracy class of the RC simulated. For the limits of ε and φ , the reader can refer to [13], which provides the maximum values allowed for each harmonic up to the 13th.

4.1.2. Sources of Uncertainty

To further validate a method, it is good practice to locate and assess the sources of uncertainty. For the proposed approach, the geometrical measures are those affected by the measurement uncertainty. It is worth mentioning that many influence quantities and many sources of uncertainty might affect the performance of RCs during normal operations. However, the uncertainty analysis run in this paper is focused on the sources that directly affect the proposed procedure. Therefore, the quantities d , h (Figure 3), and 2 a are those corrupted by uncertainty in the MCM implementation. The MCM was implemented with ten thousand iterations, assuming a uniform distribution for the probability density functions associated with the quantities. Considering realistic measurement devices used to measure lengths, the limits of the uniform distribution were fixed at 1%, 10%, and 20% of the measured quantity. To consider the worst case, the measurement of the RC windings’ resistance was also considered. Therefore, the same percentages and distributions apply to the resistance value R w . The results of the MCM implementation are listed in Table 3. For the sake of brevity, the case of the rectangular cross-section is considered.
For each quantity affected by uncertainty and for each percentage of uncertainty, the table contains the mean value μ , the lower limit l 95 % and the upper limit u 95 % of the 95% confidence interval. As expected, the spread of the results increases with the increase in uncertainty associated with the length measurements. However, the parasitic parameters are slightly affected by such a contribution. The most affected parameter is the number of turns, which is already well estimated by the method. For the sake of clarity, in Figure 6, the pdf of the parasitic capacitance C obtained running 10,000 Monte Carlo trials is plotted. Considering that C is obtained from the combination of more than one random variable affected by uncertainty, the resulting pdf is prone to a trapezoidal one.

4.2. Validation by Measurements

To assess the usefulness and applicability of the proposed approach, it is fundamental to perform validation with commercial devices. To this purpose, three RCs, namely R1, R2, and R3, were tested. Table 4 lists the characteristics of the RCs.
The RC were tested with the measurement setup depicted in Figure 7. It includes a Fluke calibrator 6105A and its transconductance Fluke 52120A. They were used together to generate the desired current a reference shunt to measure the input current injected through the RCs. The shunt resistance is, after an accurate characterization, 1.02129 mΩ ± 0.00008 mΩ, up to 2500 Hz. A data acquisition board, NI 9238, whose characteristics are given in Table 5, it was used to collect the three voltage signals from the RCs under test and the reference signal coming from the reference shunt. In Figure 7, a color code is used to distinguish between the current and the voltage signals.
The tests performed with the setup in Figure 7 are here described. The signals described in Table 2 were generated with the Fluke plus transconductance and injected through the reference shunt and the RCs under test. All tests were performed using a primary current of 100 A. The load impedance is given in Table 5, and it is greater than 1 GΩ. This preliminary measurement allows the collection of the actual RC response. In fact, the second part of the test consists of the combination of the injected signal (acquired from the shunt) and the TFs obtained from the modeling approach. The result is an estimation of the RC behavior at each frequency. To validate the estimation process, hence the proposed modeling approach, the phase error (PE) and the percentage voltage difference (VD) indicators are used. The PE is the difference in mrad between the phase displacement obtained in simulation and the one obtained during the experimental measurements (it is not the difference between absolute values; it would be meaningless). The VD, instead, is merely the percentage difference between the voltage measured by the RC and the one estimated by the method (using the rated current as the base of the ratio). The results, for the sake of brevity, are listed in Table 6 for a comprehensive selection of frequencies.
According to the results, it is possible to appreciate how accurately the model can estimate RC behavior. As expected, the accuracy drops when the frequency increases. This is due to the higher uncertainty of the generation system and the lower number of samples per period at those frequencies. However, for all RCs and at all frequencies, the obtained errors are far below the limits given in [13]. The only exception is the VD at 2500 Hz. The limit given in [13] is 20%. However, at this stage, and considering the amount of approximation conducted in the simplified procedure, the results can be considered promising. It is worth mentioning that, even if the results are all above the uncertainty limits (except one), the goodness of the results is increased by the final application in which they are going to be used. For example, applications with limited target uncertainty will benefit the most from the proposed method.

5. Discussion

The obtained results trigger discussion on the proposed method. As a matter of fact, the simulations performed and the experimental measurements provided encouraging results. They were obtained by exploiting almost any input but those included in the datasheet. The direct consequence is obvious. The model would provide far better results if one or more measurements could be performed on the RCs. Furthermore, the modeling approach can be improved if extra information is added to it. For example, imagine a RC coil being characterized vs. temperature. Its results can be included in the model, which will differentiate its output depending on the working temperature. The same consideration can be extended to any other influence quantity that may affect the RC operation. Therefore, the chances of exploitation of the modeling approach are several, such as the simple laboratory, the company that develops power system simulators, the metrological institute, etc.

6. Conclusions

This paper aims to provide a simplified modeling approach to be applied to Rogowski coils. It is designed to be applied at all levels, with the flexibility to be improved depending on the available information. The approach is then tested with both simulations and experimental measurements. The results clearly confirm the validity of the approach and pave the way for further studies needed for its improvement.

Author Contributions

Conceptualization, A.M.; methodology, R.T.; software, C.B.; validation, L.P.; formal analysis, A.M.; investigation, C.B.; resources, L.P.; data curation, C.B.; writing—original draft preparation, A.M.; writing—review and editing, R.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmed, Y.A.; Rashid, A.; Khurshid, M.M. Investigating the determinants of the adoption of solar photovoltaic Systems—Citizen’s perspectives of two developing countries. Sustainability 2022, 14, 1764. [Google Scholar] [CrossRef]
  2. Linhart, M.; Rodin, V.; Moser, S.; Kollmann, A. Citizen participation to finance pv power plants focused on self-consumption on company roofs-findings from an austrian case study. Energies 2021, 14, 738. [Google Scholar] [CrossRef]
  3. Breschi, V.; Ravazzi, C.; Strada, S.; Dabbene, F.; Tanelli, M. Fostering the mass adoption of electric vehicles: A network-based approach. IEEE Trans. Control Netw. Syst. 2022, 9, 1666–1678. [Google Scholar] [CrossRef]
  4. Lv, Z.; Qiao, L.; Cai, K.; Wang, Q. Big data analysis technology for electric vehicle networks in smart cities. IEEE Trans. Intell. Transp. Syst. 2021, 22, 1807–1816. [Google Scholar] [CrossRef]
  5. Hieu, L.T.; Lim, O.T. Effects of the structure and operating parameters on the performance of an electric scooter. Sustainability 2023, 15, 8976. [Google Scholar] [CrossRef]
  6. Kubik, A. The energy consumption of electric scooters used in the polish shared mobility market. Energies 2022, 15, 8193. [Google Scholar] [CrossRef]
  7. Razmi, P.; Ghaemi Asl, M. Machine learning approaches in a real power system and power markets. In Application of Machine Learning and Deep Learning Methods to Power System Problems; Power Systems; Springer: Cham, Switzerland, 2021. [Google Scholar] [CrossRef]
  8. Zhu, Y.; Singh, C. End-to-end topology-aware machine learning for power system reliability assessment. In Proceedings of the 2022 17th International Conference on Probabilistic Methods Applied to Power Systems, PMAPS 2022, Manchester, UK, 12–15 June 2022. [Google Scholar] [CrossRef]
  9. Chen, X.; Xu, J.; Fang, Y.; Li, W.; Ding, Y.; Wan, Z.; Tu, Z. Temperature and humidity management of PEM fuel cell power system using multi-input and multi-output fuzzy method. Appl. Therm. Eng. 2022, 203, 117865. [Google Scholar] [CrossRef]
  10. Wu, W.; Yue, S.; Zhou, X.; Guo, M.; Wang, J.; Ren, L.; Yuan, B. Observational study on the impact of large-scale photovoltaic development in deserts on local air temperature and humidity. Sustainability 2020, 12, 3403. [Google Scholar] [CrossRef]
  11. Peretto, L.; Tinarelli, R.; Ghaderi, A.; Mingotti, A.; Mazzanti, G.; Valtorta, G.; Danesi, S. Monitoring cable current and laying environment parameters for assessing the aging rate of MV cable joint insulation. In Proceedings of the Annual Report—Conference on Electrical Insulation and Dielectric Phenomena, CEIDP, Shenzhen, China, 20–23 October 2018; pp. 390–393. [Google Scholar] [CrossRef]
  12. IEC 61869-1:2010; Instrument Transformers—Part 1: General Requirements. International Standardization Organization: Geneva, Switzerland, 2010.
  13. IEC 61869-6; Instrument Transformers—Part 6: Additional General Requirements for Low-Power Instrument Transformers. International Standardization Organization: Geneva, Switzerland, 2016.
  14. Besharatifard, H.; Hasanzadeh, S.; Muyeen, S.M.; Kamwa, I. Evaluation of a calibration technique in measuring partial discharges inside mineral oils with a high-frequency current transformer (HFCT) sensor: A case study. IET Gener. Transm. Distrib. 2023, 17, 706–715. [Google Scholar] [CrossRef]
  15. Chen, Y.; Crotti, G.; Dubowik, A.; Letizia, P.S.; Mohns, E.; Luiso, M.; Bruna, J. Novel calibration systems for the dynamic and steady-state testing of digital instrument transformers. In Proceedings of the AMPS 2021—2021 11th IEEE International Workshop on Applied Measurements for Power Systems, Virtual, 29 September–1 October 2021. [Google Scholar] [CrossRef]
  16. Mingotti, A.; Peretto, L.; Tinarelli, R.; Angioni, A.; Monti, A.; Ponci, F. A simple calibration procedure for an LPIT plus PMU system under off-nominal conditions. Energies 2019, 12, 4645. [Google Scholar] [CrossRef]
  17. Kaczmarek, M.; Stano, E. Challenges of accurate measurement of distorted current and voltage in the power grid by conventional instrument transformers. Energies 2023, 16, 2648. [Google Scholar] [CrossRef]
  18. Lesniewska, E. Modern methods of construction problem solving in designing various types of instrument transformers. Energies 2022, 15, 8199. [Google Scholar] [CrossRef]
  19. Tee, S.; Liu, Q.; Wang, Z.; Hafid, F.; Tournet, P. Failure investigation and asset management of combined measuring instrument transformers. High Volt. 2021, 6, 61–70. [Google Scholar] [CrossRef]
  20. Estupiñán JS, R.; Vachoux, A.; Pascal, J. Electro-thermal modeling of a rogowski coil sensor system. In Proceedings of the 2015 International Symposium on VLSI Design, Automation and Test, VLSI-DAT, Hsinchu, Taiwan, 27–29 April 2015. [Google Scholar] [CrossRef]
  21. Ibrahim, M.E.; Abd-Elhady, A.M. Design and modeling of a two-winding rogowski coil sensor for measuring three-phase currents of a motor fed through a three-core cable. IEEE Sens. J. 2021, 21, 8289–8296. [Google Scholar] [CrossRef]
  22. Liu, X.; Huang, H.; Jiao, C. Modeling and analyzing the mutual inductance of rogowski coils of arbitrary skeleton. Sensors 2019, 19, 3397. [Google Scholar] [CrossRef] [PubMed]
  23. Ming, L.; Xin, Z.; Liu, W.; Loh, P.C. Structure and modelling of four-layer screenreturned PCB rogowski coil with very few turns for high-bandwidth SiC current measurement. IET Power Electron. 2020, 13, 765–775. [Google Scholar] [CrossRef]
  24. Mingotti, A.; Peretto, L.; Tinarelli, R. A simple modelling procedure of rogowski coil for power systems applications. In Proceedings of the AMPS 2019—2019 10th IEEE International Workshop on Applied Measurements for Power Systems, Aachen, Germany, 25–27 September 2019. [Google Scholar] [CrossRef]
  25. Mingotti, A.; Peretto, L.; Tinarelli, R. A smart frequency domain-based modeling procedure of rogowski coil for power systems applications. IEEE Trans. Instrum. Meas. 2020, 69, 6748–6755. [Google Scholar] [CrossRef]
  26. Orosz, T.; Tamus, Z.A.; Vajda, I. Modeling the high frequency behavior of the rogowski-coil passive L/r integrator current transducer with analytical and finite element method. In Proceedings of the Universities Power Engineering Conference, Cluj-Napoca, Romania, 2–5 September 2014. [Google Scholar] [CrossRef]
  27. Sabiha, N.A.; Alkhammash, H.I. Modeling rogowski coils for monitoring surge arrester discharge currents. Comput. Syst. Sci. Eng. 2021, 42, 439–449. [Google Scholar] [CrossRef]
  28. IEC 61869-10; Instrument Transformers—Part 10: Additional Requirements for Low-Power Passive Current Transformers. International Standardization Organization: Geneva, Switzerland, 2017.
  29. Styblìkovà, R.; Draxler, K. Design and calibration of Rogowski coils. In Proceedings of the Power and Energy Workshop, Noordwijk, The Netherlands, 22–23 March 2011. [Google Scholar]
  30. Moreno, M.V.R.; Robles, G.; Albarracín, R.; Rey, J.A.; Tarifa, J.M.M. Study on the self-integration of a Rogowski coil used in the measurement of partial discharges pulses. Electr. Eng. 2017, 99, 817–826. [Google Scholar] [CrossRef]
Figure 1. Picture of the Rogowski coil and its components.
Figure 1. Picture of the Rogowski coil and its components.
Sensors 23 08032 g001
Figure 2. Flowchart of the modeling procedure.
Figure 2. Flowchart of the modeling procedure.
Sensors 23 08032 g002
Figure 3. Cross-section geometries.
Figure 3. Cross-section geometries.
Sensors 23 08032 g003
Figure 4. Schematic of the equivalent circuit of the RC.
Figure 4. Schematic of the equivalent circuit of the RC.
Sensors 23 08032 g004
Figure 5. ε and φ for each test signal described in Table 2.
Figure 5. ε and φ for each test signal described in Table 2.
Sensors 23 08032 g005
Figure 6. Probability density function obtained for C after 100,000 Monte Carlo trials.
Figure 6. Probability density function obtained for C after 100,000 Monte Carlo trials.
Sensors 23 08032 g006
Figure 7. The measurement setup used for the experimental measurements on the commercial RCs.
Figure 7. The measurement setup used for the experimental measurements on the commercial RCs.
Sensors 23 08032 g007
Table 1. Typical parameters available in the RC datasheet.
Table 1. Typical parameters available in the RC datasheet.
Accuracy ClassLinearity ErrorResistance of the WindingsTransformation Ratio
Max CurrentBandwidthTemperature ClassRated Burden
Table 2. ε and φ for each test signal.
Table 2. ε and φ for each test signal.
Frequency (Hz) ε (%) φ (mrad)
50−0.013.2
100−0.016.3
150−0.019.5
200−0.0112.6
250−0.0115.8
300−0.0218.9
350−0.0222.1
400−0.0225.2
450−0.0228.4
500−0.0331.5
550−0.0334.7
600−0.0337.8
650−0.0441.0
700−0.0444.1
750−0.0547.3
800−0.0550.4
850−0.0653.6
900−0.0656.7
950−0.0759.9
1000−0.0863.0
1050−0.0866.2
1100−0.0969.3
1150−0.172.5
1200−0.1175.6
1250−0.1178.8
1300−0.1281.9
1350−0.1385.1
1400−0.1488.2
1450−0.1591.4
1500−0.1694.5
1550−0.1797.7
1600−0.18100.9
1650−0.19104.0
1700−0.2107.2
1750−0.21110.3
1800−0.23113.5
1850−0.24116.6
1900−0.25119.8
1950−0.26122.9
2000−0.28126.1
2050−0.29129.2
2100−0.3132.4
2150−0.32135.5
2200−0.33138.7
2250−0.35141.8
2300−0.36145.0
2350−0.38148.1
2400−0.39151.3
2450−0.41154.4
2500−0.43157.6
Table 3. Results of the uncertainty analysis.
Table 3. Results of the uncertainty analysis.
Uncertainty Considered (%)Quantity μ l 95 % u 95 %
1 N (-)134013091370
10 N (-)135010541665
20 N (-)14028372092
1 C (F) 3.8 × 10 12 3.8 × 10 12 3.9 × 10 12
10 C (F) 3.9 × 10 12 3.4 × 10 12 4.3 × 10 12
20 C (F) 3.9 × 10 12 3.0 × 10 12 4.6 × 10 12
1 L (H) 4.3 × 10 4 4.2 × 10 4 4.4 × 10 4
10 L (H) 4.3 × 10 4 3.3 × 10 4 5.3 × 10 4
20 L (H) 4.5 × 10 4 2.7 × 10 4 6.6 × 10 4
Table 4. Characteristics of the commercial RCs under test.
Table 4. Characteristics of the commercial RCs under test.
FeatureR1R2R3
T R 50 (mV/kA)100100100
d (mm)8812
a (mm)505721.5
R w (Ω)25638122
Cross-SectionCircularCircularOval
Accuracy Class0.511
Bandwidth1 Hz to 100 kHzNA20 Hz to 5 kHz
Operating Temperature−30 °C to 80 °C−20 °C to 85 °C−20 °C to 70 °C
Rated Current (A)100010,0001000
Table 5. Main characteristics of the data acquisition board NI9238.
Table 5. Main characteristics of the data acquisition board NI9238.
Converter24-bitVoltage Range±500 mV
Sampling Frequency50 kSa/s/ChInput Impedance>1 GΩ
Simultaneous ChannelsYesGain Error±0.07%
Offset Error±0.005%Input Noise 3.9   μ V
Table 6. Results from the comparison of the actual and estimated RC behaviors.
Table 6. Results from the comparison of the actual and estimated RC behaviors.
Frequency (Hz)R1R2R3
PE (mrad)VD (%)PE (mrad)VD (%)PE (mrad)VD (%)
502.920.4−2.800.5−4.411.8
25014.361.2−12.831.0−15.884.2
55031.392.7−26.911.5−33.668.9
85048.344.6−40.711.7−51.7313.5
125070.935.3−59.100.8−75.9113.7
2500141.7715.8116.92−5.1−151.3122.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mingotti, A.; Betti, C.; Tinarelli, R.; Peretto, L. Simplifying Rogowski Coil Modeling: Simulation and Experimental Verification. Sensors 2023, 23, 8032. https://doi.org/10.3390/s23198032

AMA Style

Mingotti A, Betti C, Tinarelli R, Peretto L. Simplifying Rogowski Coil Modeling: Simulation and Experimental Verification. Sensors. 2023; 23(19):8032. https://doi.org/10.3390/s23198032

Chicago/Turabian Style

Mingotti, Alessandro, Christian Betti, Roberto Tinarelli, and Lorenzo Peretto. 2023. "Simplifying Rogowski Coil Modeling: Simulation and Experimental Verification" Sensors 23, no. 19: 8032. https://doi.org/10.3390/s23198032

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

Article Metrics

Back to TopTop