<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1366283</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2024.1366283</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>A DC fault current fast-computing method of MMC-HVDC grid with short circuit protection equipment</article-title>
<alt-title alt-title-type="left-running-head">Zhang et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2024.1366283">10.3389/fenrg.2024.1366283</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhang</surname>
<given-names>Xiong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2619697/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Xu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhuo</surname>
<given-names>Chaoran</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2071670/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>State Key Laboratory of Electrical Insulation and Power Equipment</institution>, <institution>School of Electrical Engineering</institution>, <institution>Xi&#x2019;an Jiaotong University</institution>, <addr-line>Xi&#x2019;an</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Electrical Engineering</institution>, <institution>Xi&#x2019;an University of Technology</institution>, <addr-line>Xi&#x2019;an</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1259738/overview">Jian Zhao</ext-link>, Shanghai University of Electric Power, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1699858/overview">Ma Jianjun</ext-link>, Shanghai Jiao Tong University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2650878/overview">Hua Geng</ext-link>, Tsinghua University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2651218/overview">Jinghua Zhou</ext-link>, North China University of Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Xiong Zhang, <email>zzx81765392@163.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>08</day>
<month>04</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1366283</elocation-id>
<history>
<date date-type="received">
<day>06</day>
<month>01</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>25</day>
<month>03</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Zhang, Yang and Zhuo.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Zhang, Yang and Zhuo</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The multi-terminal modular multi-level converter-based high voltage direct current (MMC-HVDC) grid with short circuit protection equipment (SCPE) is so complex that it is difficult to estimate its fault current and analyze the performance of SCPE by conventional time-domain numerical calculation method, it meets three big obstacles. This paper has made significant progress in overcoming these obstacles. 1). By applying the modern electrical circuit theory, a systematic formulation of the differential equation set for fault current calculation is developed to avoid a lot of complex and cumbersome matrix manual calculations. 2). A novel Y-Delta transformation in the s-domain is proposed to develop an eliminating virtual node approach for a complex MMC-HVDC grid, including the ring, radial, and hybrid topologies. 3). It is difficult to solve the equivalent circuit of MMC-HVDC grid with SCPE since SCPE is a time-variable-nonlinear circuit. A canonical voltage source model of SCPE is established to transform the time-variable-nonlinear circuit into a piecewise linear circuit. Based on the three significant progresses, a DC fault current fast-computing method of MMC-HVDC grid with SCPE is put forward to deal with all kinds of MMC-HVDC grids with several kinds of SCPEs. Then, the performance of several kinds of SCPE is analyzed and compared by this method. Consequently, the proposed DC fault current fast-computing method is a new powerful tool to estimate the fault current of MMC-HVDC grid and analyze the performance of SCPE.</p>
</abstract>
<kwd-group>
<kwd>HVDC</kwd>
<kwd>DC circuit breaker</kwd>
<kwd>dc fault current calculation</kwd>
<kwd>modular multilevel converter</kwd>
<kwd>fault current limiting circuit breaker</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Process and Energy Systems Engineering</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Compared to point-to-point transmission, the multi-terminal high-voltage direct current (HVDC) grid offers several advantages, such as connecting a larger number of renewable energy sources and provides higher power supply reliability, making it a key area of focus in future grid research (<xref ref-type="bibr" rid="B11">Li et al., 2021</xref>). The voltage-source converter (VSC) technologies, specifically those based on the modular multilevel converter (MMC), are rapidly expanding the applications of high-voltage direct current (<xref ref-type="bibr" rid="B21">Shu et al., 2023</xref>). VSC-HVDC systems, which utilize MMC, are particularly well-suited for multi-terminal DC (MTDC) systems due to their ability to flexibly control current direction and redistribute power between terminals (<xref ref-type="bibr" rid="B3">An et al., 2017</xref>).</p>
<p>However, the issue of DC short circuit fault protection poses a significant challenge and acts as a major obstacle to the development and implementation of multi-terminal HVDC grids with overhead lines. When a DC grid experiences a short-circuit fault, it is necessary to insert short circuit protection equipment into the fault path to limit the rising rate and break current of the fault current, as well as isolate the short circuit fault (<xref ref-type="bibr" rid="B1">Ahmad et al., 2022</xref>; <xref ref-type="bibr" rid="B31">Zhang et al., 2022</xref>). The presence of SCPE within a multi-terminal MMC-HVDC grid adds complexity to the system, making it highly challenging to accurately calculate fault currents and analyze the performance of SCPE within the grid.</p>
<p>The Electromagnetic-transient-type (EMT) simulation tool is capable of accurately simulating the performance of an HVDC grid under normal or fault conditions (H. <xref ref-type="bibr" rid="B19">Saad et al., 2013</xref>). The EMT model of the MMC converter accurately represents thousands of IGBT switching events simultaneously. A multi-terminal HVDC system based on MMC consists of multiple MMC converters (<xref ref-type="bibr" rid="B24">Stepanov et al., 2021</xref>). However, simulating a multi-terminal MMC-HVDC grid using EMT model requires a lot of computing resources and is time-consuming, because the multi-terminal MMC-HVDC grid is a time-variable and nonlinear circuit. To improve computational efficiency, several numerical computing approaches have been developed. The literature (<xref ref-type="bibr" rid="B15">LI et al., 2017a</xref>) provides a detailed analysis of MMC arm fault current before the converter blocks during a pole-to-pole fault in a single-terminal HVDC system. Additionally, an RLC equivalent circuit model of the MMC under short fault conditions has been reported by (<xref ref-type="bibr" rid="B5">Belda et al., 2018</xref>; <xref ref-type="bibr" rid="B29">Yang et al., 2018</xref>).</p>
<p>Moreover, the overhead line is approximated as an R-L circuit in (<xref ref-type="bibr" rid="B13">Li et al., 2016</xref>; <xref ref-type="bibr" rid="B14">Li et al., 2022</xref>). Based on the RLC equivalent circuit model of MMC and the R-L circuit of overhead lines, a linear RLC model of MMC-HVDC grid has been established in literature (<xref ref-type="bibr" rid="B12">Li et al., 2017b</xref>) for calculating short-circuit currents. This model is referred to as the RLC equivalent circuit in this paper and has made significant progress in the calculation of short-circuit current by representing the time-variable-nonlinear circuit of the multi-terminal HVDC grid as an RLC linear model. In addition to the RLC equivalent circuit, the literature (<xref ref-type="bibr" rid="B15">Li et al., 2017a</xref>) proposed a general time-domain numerical calculation method for analyzing the pole-to-pole short-circuit fault of the pseudo-bipolar multi-terminal MMC-HVDC grid. Furthermore, the literature (<xref ref-type="bibr" rid="B18">Ning et al., 2019</xref>; <xref ref-type="bibr" rid="B30">Ye et al., 2021</xref>) introduced how to utilize this method to calculate and analyze pole-to-pole and pole-to-ground short-circuit faults of bipolar multi-terminal MMC-HVDC grids, respectively.</p>
<p>The conventional time-domain numerical calculation method can be divided into two steps. The first step is to establish the state equation set or the differential equation set of the RLC equivalent circuit of the MMC-HVDC grid by using Kirchhoff&#x2019;s Voltage Law (KVL), Kirchhoff&#x2019;s Current Law (KCL), and V&#xa0;A Relation (VAR). The second step is solving the equation to obtain the transient value of the fault current of the multi-terminal MMC-HVDC grid. However, despite being an important progress, the general short-circuit current calculation method of a DC grid presents several problems that need to be addressed. Please note that in the following section, the term &#x201c;DC grid&#x201d; specifically refers to the MMC-HVDC grid.<list list-type="simple">
<list-item>
<p>1) It is challenging to write the state equation set or differential equation set for the RLC equivalent circuit due to the involvement of multiple matrix calculations, making the writing process complex and cumbersome. To address this limitation, this paper presents a systematic formulation of the differential equation set.</p>
</list-item>
<list-item>
<p>2) If the structure of the multi-terminal DC grid is a ring, the number of equations is identical to that of the state variables in the RLC equivalent circuit. In this case, one can easily use MATLAB to solve the equations and predict the short-circuit current curves. However, if the structure is radial or hybrid, there is a virtual node in the RLC equivalent circuit, resulting in fewer equations than state variables. It becomes difficult to obtain the transient solution of the equations. Therefore, this paper proposes a novel Y-Delta transformation in the s-domain to develop an approach that eliminates the virtual node in a complex DC grid, including the ring, radial, and hybrid topologies.</p>
</list-item>
</list>
</p>
<p>In the short-circuit fault current generic calculation method based on short-circuit equivalent RLC model of the DC grid, although the SCPEs have been equipped in the DC grid, these have to be ignored since it is difficult to model the SCPEs. In this case, this method is mainly used to analyze and evaluate the development of fault currents after short-circuiting faults in DC grids without considering the effect of SCPEs, resulting in bigger analysis errors. In order to improve the analyzing accuracy, a short-circuit fault current calculation method of DC grid with SCPE is developed in this paper.</p>
<p>In this paper, the SCPEs are classified into three categories: fault current limiter (FCL) (<xref ref-type="bibr" rid="B20">Safaei et al., 2020</xref>), current limiting circuit breaker (CLCB), and DC circuit breaker (DCCB). The protection scheme for short-circuit faults in the DC grid is divided into three options: FCL &#x2b; DCCB scheme, DCCB scheme, or only CLCB scheme. However, regardless of the chosen protection scheme, establishing a short-circuit fault equivalent model for DC grids with SCPEs is challenging due to the presence of solid-state switches, mechanical ultra-fast disconnectors (UFD), and metal-oxide arresters (MOA) in the SCPEs.</p>
<p>In the literature (<xref ref-type="bibr" rid="B28">Xu et al., 2019</xref>), an improved general calculation method is presented for calculating the DC fault current of a DC grid with hybrid FCL. The CLCB integrates isolating fault and limiting fault currents into a single device, reducing the cost of SCPE. According to the current limiting mechanism, the published short circuit protection equipment for the MTDC grid can be divided into MOA-type current limiting circuit breaker (CLCB) (<xref ref-type="bibr" rid="B9">Hedayati and Jovcic, 2018</xref>; <xref ref-type="bibr" rid="B23">Song et al., 2019</xref>), inductance type CLCB (<xref ref-type="bibr" rid="B10">Heidary et al., 2020</xref>; <xref ref-type="bibr" rid="B25">Wang et al., 2020</xref>) and capacitor type CLCB (<xref ref-type="bibr" rid="B27">Wu et al., 2020</xref>; <xref ref-type="bibr" rid="B26">Wu et al., 2022</xref>; <xref ref-type="bibr" rid="B32">Zhang et al., 2023</xref>).</p>
<p>Although there are several types of CLCBs, a canonical voltage source model is proposed in this paper to describe their external electrical characteristic. Importantly, this model can also be used to describe the external characteristics of other protection schemes employed in the DC grid, such as FCL &#x2b; DCCB schemes or only DCCB schemes. By utilizing the canonical voltage source model of SCPEs, along with the RLC equivalent circuit of MMC and the RL equivalent circuit of the overhead line, the short-circuit fault equivalent model of the DC grids with SCPEs can be transformed into a linear circuit.</p>
<p>Hence, based on three big progresses, a fault current fast-computing method of MMC-HVDC grid with SCPE is proposed. This method is designed to handle various structures of DC grids with SCPEs and offers several advantages, including generality, ease of implementation, and fast computation. So, this method helps design and quickly screen and evaluate the topology, configuration scheme and timing logic of the SCPEs. In this paper, the proposed calculation method analyzes and compares the performance of three typical CLCBs.</p>
<p>This paper is organized as follows: A systematic formulation of the differential equation set of DC grid and eliminating virtual node approach for a radial DC grid are respectively proposed in <xref ref-type="sec" rid="s2">Section 2</xref> and <xref ref-type="sec" rid="s3">Section 3</xref>. <xref ref-type="sec" rid="s4">Section 4</xref> establishes a canonical voltage source model of CLCBs to describe the external electrical characteristics of three typical CLCBs. The short circuit fault current calculation method of the DC grid with CLCB is developed in <xref ref-type="sec" rid="s5">Section 5</xref>, and its computing accuracy is validated in <xref ref-type="sec" rid="s6">Section 6</xref>. <xref ref-type="sec" rid="s7">Section 7</xref> analyses and compares the performances of three typical CLCBs by the proposed method, and the conclusion is introduced in <xref ref-type="sec" rid="s8">Section 8</xref>.</p>
</sec>
<sec id="s2">
<title>2 Formulation of the differential equation set</title>
<sec id="s2-1">
<title>2.1 RLC equivalent circuit of MTDC grid</title>
<p>In this paper, a four-terminal DC grid system is taken as an example, as shown in <xref ref-type="fig" rid="F1">Figure 1A</xref>. The system consists of four VSCs, referred to as VSC<sub>1</sub>, VSC<sub>2</sub>, VSC<sub>3</sub>, and VSC<sub>4</sub>. There are four pairs of overhead lines, and their RLC equivalent circuit model (<xref ref-type="bibr" rid="B12">Li et al., 2017b</xref>) is shown in <xref ref-type="fig" rid="F1">Figure 1B</xref> for short-circuit fault current calculation. According to (<xref ref-type="bibr" rid="B5">Belda et al., 2018</xref>; <xref ref-type="bibr" rid="B29">Yang et al., 2018</xref>), MMC converters can be represented by an RLC equivalent circuit under short-circuit faults. In this model, a series RLC branch is used to represent the VSC<sub>
<italic>i</italic>
</sub>, which includes an internal resistance <italic>R</italic>
<sub>
<italic>i</italic>
</sub>, limiting current inductance <italic>L</italic>
<sub>
<italic>i</italic>
</sub> and discharging capacitance <italic>C</italic>
<sub>
<italic>i</italic>
</sub>, where <italic>i</italic> &#x3d; 1,2,3,4. According to transmission line theory (<xref ref-type="bibr" rid="B13">Li et al., 2016</xref>), a pair of overhead line is modeled as an RL series circuit, with <italic>R</italic>
<sub>
<italic>ij</italic>
</sub> and <italic>L</italic>
<sub>
<italic>ij</italic>
</sub> representing the overhead line between VSC<sub>
<italic>i</italic>
</sub> and VSC<sub>
<italic>j</italic>
</sub>. The nodes <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:msubsup>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>/</mml:mo>
<mml:msubsup>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> representnthe positive/negative pole of VSC<sub>
<italic>i</italic>
</sub>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>
<bold>(A)</bold> four-terminal DC grid; <bold>(B)</bold> RLC equivalent circuit.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g001.tif"/>
</fig>
<p>To establish the RLC equivalent circuit, it is assumed that a short-circuit fault occurs in the middle position of the power transmission line between VSC<sub>1</sub> and VSC<sub>2</sub>. So, the number of pairs of overhead lines increases from four to five. Node <italic>n</italic>
<sub>0</sub> is the fault point in <xref ref-type="fig" rid="F1">Figure 1B</xref>, and <italic>R</italic>
<sub>10</sub>/<sub>20</sub> and <italic>L</italic>
<sub>10</sub>/<sub>20</sub> are the equivalent resistance and inductance of the overhead line between nodes <italic>n</italic>
<sub>1</sub>/<sub>2</sub> and <italic>n</italic>
<sub>0</sub>. The parameter calculation method for this model has been developed in (<xref ref-type="bibr" rid="B15">Li et al., 2017a</xref>).</p>
</sec>
<sec id="s2-2">
<title>2.2 Systematic formulation of the loop equations</title>
<p>In the RLC equivalent circuit, the circular arrows indicate the orientations of the loops chosen for writing the KVL equation. We conceive of fictitious circulating loop currents, with references given by the loop orientations. Examination of the RLC equivalent circuit shows that these loop currents are identical with the branch currents <italic>i</italic>
<sub>1</sub>, <italic>i</italic>
<sub>2</sub>, <italic>i</italic>
<sub>3</sub>, <italic>i</italic>
<sub>4</sub> and <italic>i</italic>
<sub>5</sub>, which are the currents through the equivalent inductances <italic>L</italic>
<sub>10</sub>, <italic>L</italic>
<sub>20</sub>, <italic>L</italic>
<sub>24</sub>, <italic>L</italic>
<sub>34</sub>, and <italic>L</italic>
<sub>13</sub>. So, these loop currents are the state variable of the equivalent inductances of the overhead lines.</p>
<p>On the other hand, the voltages (<italic>u</italic>
<sub>c1</sub>, <italic>u</italic>
<sub>c2</sub>, <italic>u</italic>
<sub>c3</sub> and <italic>u</italic>
<sub>c4</sub>) across the discharging capacitance <italic>C</italic>
<sub>1</sub>, <italic>C</italic>
<sub>2</sub>, <italic>C</italic>
<sub>3</sub>, and <italic>C</italic>
<sub>4</sub> are also another set of state variables. To obtain a compact state equation, the loop currents rector <bold>
<italic>i</italic>
</bold>, the capacitance voltages rector <bold>
<italic>u</italic>
</bold> and its current rector <bold>
<italic>i</italic>
</bold>
<sub>
<italic>c</italic>
</sub> are defined as follows,<disp-formula id="e1">
<mml:math id="m2">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>4</mml:mn>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>4</mml:mn>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>By using KVL and VAR, the loop current equation can be obtained in matrix form,<disp-formula id="e2">
<mml:math id="m3">
<mml:mrow>
<mml:mi mathvariant="bold">A</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">R</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold">L</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>Where <bold>A</bold> is an incidence matrix, <bold>R</bold> and <bold>L</bold> are parameter matrices, referred to respectively as the resistance matrix and the inductance matrix, which will be introduced later.</p>
</sec>
<sec id="s2-3">
<title>2.3 Incidence and parameter matrix</title>
<sec id="s2-3-1">
<title>2.3.1 The incidence matrix A</title>
<p>For an RLC equivalent circuit, if there are <italic>n</italic> capacitors and <italic>b</italic> pairs of overhead lines, an incidence matrix <bold>A</bold> &#x3d; [<italic>a</italic>
<sub>
<italic>ij</italic>
</sub>] is an <italic>n</italic>&#xd7;<italic>b</italic> rectangular matrix. In this paper, the RLC equivalent circuit shown in <xref ref-type="fig" rid="F1">Figure 1A</xref> will be taken as an example to introduce how to form the incidence matrix <bold>A</bold> &#x3d; [<italic>a</italic>
<sub>
<italic>ij</italic>
</sub>].</p>
<p>Firstly, it is necessary to draw an oriented graph of the RLC equivalent circuit. The oriented graph associated with the RLC equivalent circuit in <xref ref-type="fig" rid="F1">Figure 1A</xref> is illustrated in <xref ref-type="fig" rid="F2">Figure 2A</xref>. Each node is associated with a discharging capacitor. For example, the capacitance <italic>C</italic>
<sub>
<italic>i</italic>
</sub> is repressed as node n<sub>
<italic>i</italic>
</sub>, where <italic>i</italic> &#x3d; 1, 2, 3, 4. Each branch is associated with a pair of overhead lines, and its orientation is identical to the assumed current flow. The five pairs of overhead lines are expressed as b<sub>1</sub>, b<sub>2</sub>, b<sub>3</sub>, b<sub>4</sub> and b<sub>5</sub>, respectively, in <xref ref-type="fig" rid="F2">Figure 2A</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Oriented graph and its incidence matrix: <bold>(A)</bold> oriented graph <bold>(B)</bold> Incidence matrix.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g002.tif"/>
</fig>
<p>Based on the oriented graph, the elements of incidence matrix <bold>A</bold> have the following values (<xref ref-type="bibr" rid="B12">Li et al., 2017b</xref>):</p>
<p>
<italic>a</italic>
<sub>
<italic>ij</italic>
</sub> &#x3d; 1, if branch <italic>j</italic> is incident at node <italic>i</italic> and oriented away from it; <italic>a</italic>
<sub>
<italic>ij</italic>
</sub> &#x3d; &#x2212;1, if branch <italic>j</italic> is incident at node <italic>i</italic> and oriented toward it; <italic>a</italic>
<sub>
<italic>ij</italic>
</sub> &#x3d; 0, if branch <italic>j</italic> is not incident at node <italic>i</italic>. For the oriented graph shown in <xref ref-type="fig" rid="F2">Figure 2A</xref>, the incidence matrix <bold>A</bold> is shown in <xref ref-type="fig" rid="F2">Figure 2B</xref>.</p>
</sec>
<sec id="s2-3-2">
<title>2.3.2 Parameter matrices of <bold>R</bold> and <bold>L</bold>
</title>
<p>In this paper, a systematic formulation approach for the parameter matrices of <bold>R</bold> and <bold>L</bold> is also proposed as follows. By relating parameter matrices to the RLC equivalent, the elements of the parameter matrices can be obtained using the following straightforward way.</p>
<p>Each term on the main diagonal is the sum of the resistance/inductance value of the branches on the corresponding loop. Each off-diagonal term is plus or minus the resistance/inductance value of branches common between two loops. The sign is positive if the loop currents traverse the common branch with the same orientation, and negative if they traverse the common branch with opposite orientations. Verify the <bold>R</bold> parameter matrix for the example by using this method.<disp-formula id="e3">
<mml:math id="m4">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold">R</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>10</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>20</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>24</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>34</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>13</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>p</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>p</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>p</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>p</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>p</mml:mi>
<mml:mn>5</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
</sec>
</sec>
<sec id="s2-4">
<title>2.4 Systematic formulation of the node equation</title>
<p>In <xref ref-type="fig" rid="F1">Figure 1B</xref>, four capacitances, such as <italic>C</italic>
<sub>1</sub>, <italic>C</italic>
<sub>2</sub>, <italic>C</italic>
<sub>3</sub> and <italic>C</italic>
<sub>4</sub>, are associated with four independent capacitance-voltage state variables. On the other hand, although there are fourteen inductances, there are only five independent inductance current state variables since each independent loop corresponds to a unique independent inductance (<xref ref-type="bibr" rid="B4">Balabanian and Bickart, 1969</xref>). Therefore, the RLC equivalent circuit consists of nine state variables, but the Eq. <xref ref-type="disp-formula" rid="e2">2</xref> include only five loop current equations. Thus, it is necessary to write the remaining four equations for the state variables.</p>
<p>Application of KCL, the node current equations for nodes n<sub>1</sub>, n<sub>2</sub>, n<sub>3</sub> and n<sub>4</sub> can be written as<disp-formula id="e4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi mathvariant="bold">A</mml:mi>
<mml:mi mathvariant="bold">T</mml:mi>
</mml:msup>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>It is observed from the Eq. <xref ref-type="disp-formula" rid="e4">4</xref> that the capacitance current is not an independent state variable, which the independent inductance current state variables can express.</p>
<p>According to the constraint relationship of the capacitance voltage and current through each VSC equivalent circuit, their constraint equations are written as<disp-formula id="e5">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">C</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">C</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2260;</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>Where C is a diagonal matrix.</p>
<p>Substituting the equation Eq. <xref ref-type="disp-formula" rid="e4">4</xref> into Eq. <xref ref-type="disp-formula" rid="e5">5</xref>, we can get the node equation,<disp-formula id="e6">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">P</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">P</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold">C</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold">A</mml:mi>
<mml:mi mathvariant="bold">T</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>The final differential equation set are given as<disp-formula id="e7">
<mml:math id="m8">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold">A</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">R</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold">L</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">P</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>By using the Eq. <xref ref-type="disp-formula" rid="e7">7</xref>, the DC short fault currents of MTDC grid can be promptly and accurately calculated.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Eliminating virtual node approach for a radial DC grid</title>
<p>There are three general topologies of DC grid, referred to as the ring, radial, and hybrid topologies. According to the definition of node attributes (<xref ref-type="bibr" rid="B15">Li et al., 2017a</xref>), nodes n<sub>1</sub>, n<sub>2</sub>, n<sub>3</sub>, and n<sub>4</sub> are real nodes, but node n<sub>5</sub> is defined as a virtual node in the radial shown in <xref ref-type="fig" rid="F3">Figure 3</xref> because it is not connected to any VSC. This virtual node will result in an important issue. The number of the differential equations of Eq. <xref ref-type="disp-formula" rid="e7">7</xref> will be less than that of state variables in the RLC equivalent circuit, which forms an underdetermined equation set. The reason is that we cannot find any constraint equation of the virtual node voltage and it&#x2019;s current. An eliminating virtual node approach has been investigated using the Y-Delta transformation (<xref ref-type="bibr" rid="B2">Alexander and Sadiku, 2013</xref>) theory in this section for this type of DC grid with a virtual node.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>A radial four-terminal DC grid with virtual node and its equivalent circuit: <bold>(A)</bold> A radial four-terminal DC grid; <bold>(B)</bold> RLC equivalent circuit.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g003.tif"/>
</fig>
<sec id="s3-1">
<title>3.1 Y-Delta Transformation in the s-domain</title>
<p>
<xref ref-type="fig" rid="F4">Figure 4A</xref> shows the equivalent circuit of Y-Delta transformation in the <italic>s</italic>-domain for the virtual node <italic>n</italic>
<sub>5</sub> in <xref ref-type="fig" rid="F3">Figure 3</xref>, which is a time-domain circuit. There are two circuits named interior and exterior circuit, respectively. The interior circuit is a Why-type circuit, and the exterior circuit is its Delta-type equivalent circuit. In the <italic>s</italic>-domain circuit, an inductance <italic>L</italic> is expressed as an inductance <italic>L</italic> series with a voltage source <italic>L</italic>&#xd7; <italic>i</italic> (0), where <italic>i</italic> (0) is the initial current value through the inductance <italic>L</italic>. The resistance <italic>R</italic> shares an identical form in the time domain.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>
<bold>(A)</bold> Equivalent circuit of Y-Delta Transformation for the virtual node n5. <bold>(B)</bold> Natural response circuit <bold>(C)</bold> Force response circuit.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g004.tif"/>
</fig>
<p>In order to determine the parameters of the Delta equivalent circuit, labelled as <italic>L</italic>
<sub>&#x394;</sub>, <italic>R</italic>
<sub>&#x2206;</sub> and <italic>i</italic>
<sub>&#x394;</sub>, a circuit complete response can be broken into the natural response and the forced response (<xref ref-type="bibr" rid="B2">Alexander and Sadiku, 2013</xref>). Therefore, the circuit shown in <xref ref-type="fig" rid="F4">Figure 4A</xref> can also be separated into the natural and forced response circuits, as shown in <xref ref-type="fig" rid="F4">Figures 4B,C</xref>.</p>
<p>For the natural circuit shown in <xref ref-type="fig" rid="F4">Figure 4B</xref>, the resistance <italic>R</italic> and inductance <italic>L</italic> are considered short circuits, because the input currents <italic>i</italic>
<sub>1</sub>, <italic>i</italic>
<sub>2</sub> and <italic>i</italic>
<sub>3</sub> are zero. By using KVL, the inductance current initial value of the Delta equivalent circuit can be obtained as,<disp-formula id="e8">
<mml:math id="m9">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>41</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>54</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>41</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>34</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>54</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>54</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>34</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>13</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mo>&#x2212;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>13</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>For the force circuit shown in <xref ref-type="fig" rid="F4">Figure 4C</xref>, all voltage sources can be considered short circuits since the initial value of inductance current equals zero. By using the conversion rule for Wye to Delta (<xref ref-type="bibr" rid="B2">Alexander and Sadiku, 2013</xref>), the following formula can be obtained,<disp-formula id="e9">
<mml:math id="m10">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>13</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>13</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>34</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>34</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>41</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>41</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>In DC grid, the equivalent resistance of an overhead line is much bigger than the equivalent inductance, Eq. <xref ref-type="disp-formula" rid="e9">9</xref> can be simplified, and the <italic>R</italic> and <italic>L</italic> expressions of the Delta equivalent circuit can be written as.<disp-formula id="e10">
<mml:math id="m11">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>13</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2248;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>13</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2248;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>34</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2248;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>34</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2248;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>41</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2248;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>35</mml:mn>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mn>41</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2248;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>15</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>45</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
</sec>
<sec id="s3-2">
<title>3.2 Equivalent circuit without virtual node</title>
<p>By using the proposed Y-Delta transformation in the <italic>s</italic>-domain, the virtual node <italic>n</italic>
<sub>5</sub> in the radial four-terminal dc grid shown in <xref ref-type="fig" rid="F3">Figure 3</xref> can be eliminated, obtaining its equivalent circuit without virtual node, as shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. In this equivalent circuit, there is not any virtual node. So, one can use the proposed systematic formulation of the differential equation set developed in <xref ref-type="sec" rid="s2">Section 2</xref> to write its differential equation set.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Equivalent circuit without virtual node.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g005.tif"/>
</fig>
</sec>
</sec>
<sec id="s4">
<title>4 Canonical voltage source model of CLCB</title>
<p>According to the current limiting mechanism, the published DC circuit breaker with current limiting function can be divided into three types: MOA type CLCB, inductance type CLCB, and capacitor type CLCB, referred to respectively as MOA-CLCB, inductance-CLCB, and capacitor-CLCB. Although there are several types of CLCB, this paper proposes a canonical voltage source model, as shown in <xref ref-type="fig" rid="F6">Figure 6A</xref>, to describe their external electrical characteristics. In the CLCB topology, MOA is generally employed to provide different voltages counteracting the fault current (<xref ref-type="bibr" rid="B17">Mohammadi et al., 2021</xref>).</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Canonical voltage source model of CLCB, I-V curve of MOA and its piecewise linear equivalent circuit: <bold>(A)</bold> Canonical voltage source model of CLCB. <bold>(B)</bold> I-V curve of MOA. <bold>(C)</bold> Piecewise linearization V-I characteristic of the MOV. <bold>(D)</bold> Piecewise linearization equivalent circuit of the MOV.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g006.tif"/>
</fig>
<sec id="s4-1">
<title>4.1 Piecewise linear model of MOA</title>
<p>
<xref ref-type="fig" rid="F6">Figure 6B</xref> shows a typical I-V characteristic curve of MOA, which is generally separated into four regions: the normal operational region, pre-breakdown region, breakdown region, and high current region (<xref ref-type="bibr" rid="B16">Martinez and Durbak, 2005</xref>). There are two voltage parameters, <italic>U</italic>
<sub>r</sub> is the rated voltage and <italic>U</italic>
<sub>ref</sub> is the reference voltage.</p>
<p>In CLCB, MOA is utilized to clamp the maximum to <italic>U</italic>
<sub>ref</sub> and dissipate the energy of fault current. It is difficult to analyze a CLCB circuit because MOA is a nonlinear circuit element. A piece-wise linear I-V curve of MOA, as shown in <xref ref-type="fig" rid="F6">Figure 6C</xref>, is used for the pre-breakdown region and breakdown region to simplify analysis of CLCB circuit, since it is allowed to operate in the normal operational region, pre-breakdown region and breakdown region.</p>
<p>As shown in <xref ref-type="fig" rid="F6">Figure 6C</xref>, four piece straight lines are used to describe the I-V curve, and its equivalent circuit is depicted in <xref ref-type="fig" rid="F6">Figure 6D</xref>. When the MOA operates in the normal operational region, the current flowing through it is approximately several mA, it can be considered as an open circuit. Consequently, the switch S<sub>0</sub> remains open in the equivalent circuit. If <italic>U</italic>
<sub>r</sub> &#x2266; <italic>u</italic>
<sub>moa</sub> &#x3c; <italic>U</italic>
<sub>a</sub>, both switches of <italic>S</italic>
<sub>
<italic>0</italic>
</sub> and <italic>S</italic>
<sub>1</sub> close, but the switches of <italic>S</italic>
<sub>2</sub>, <italic>S</italic>
<sub>3</sub>, and <italic>S</italic>
<sub>4</sub> remain open. In this scenario, the equivalent circuit is a voltage source <italic>U</italic>
<sub>
<italic>r</italic>
</sub> series with the resistor <italic>R</italic>
<sub>a</sub>. Hence, the piecewise linear model of MOA can be respectively expressed as.<disp-formula id="e11">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x221e;</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>g</mml:mi>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>g</mml:mi>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>g</mml:mi>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m13">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
</sec>
<sec id="s4-2">
<title>4.2 Modelling of inductance-CLCB</title>
<p>
<xref ref-type="fig" rid="F7">Figure 7A</xref> illustrates an inductive-CLCB circuit [21]. The circuit consists of <italic>L</italic>
<sub>
<italic>1</italic>
</sub>, which functions as the limiting current inductor of the VSC and serves as the primary side of the transformer, and <italic>L</italic>
<sub>
<italic>2</italic>
</sub>, which acts as the secondary side. Under normal conditions, IGBT <italic>T</italic>
<sub>1</sub>, <italic>T</italic>
<sub>3</sub> and UFD<sub>1</sub> are turned on to the normal current to flow through <italic>L</italic>
<sub>1</sub>-<italic>T</italic>
<sub>1</sub>-UFD<sub>1</sub>-<italic>T</italic>
<sub>3</sub>, as indicated by the blue dotted line in <xref ref-type="fig" rid="F7">Figure 7A</xref>. In this case, there is no current flowing through <italic>L</italic>
<sub>2</sub>, because it is a high impedance branch compared with the turned on <italic>T</italic>
<sub>1</sub> branch. The sequential control strategy of inductive-CLCB is illustrated in <xref ref-type="fig" rid="F7">Figure 7A</xref>. According to the strategy, the modelling process of inductance-CLCB will be developed as follows.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Three typical CLCBs CLCB and its sequential control strategy: <bold>(A)</bold> Inductance-CLCB circuit. <bold>(B)</bold> MOA-CLCB circuit. <bold>(C)</bold> Capacitance-CLCB.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g007.tif"/>
</fig>
<p>Assuming a short-circuit fault occurs at time <italic>t</italic>
<sub>
<italic>0</italic>
</sub>, the inductance-CLCB does not operation in this interval due to the delay in the detecting signal. As a result, the fault current follows the same path as in the normal condition.</p>
<sec id="s4-2-1">
<title>4.2.1 Current stepping interval (<italic>t</italic>
<sub>
<italic>1</italic>
</sub>
<italic>&#x223c;t</italic>
<sub>
<italic>2</italic>
</sub>)</title>
<p>At <italic>t</italic>
<sub>1</sub>, the CLCB receives the break signal, and lets IGBT <italic>T</italic>
<sub>1</sub> be turned on, resulting in the fault current being forced to flow through the coupling inductor <italic>L</italic>
<sub>2</sub>. At the same time, UFD<sub>2</sub> also receives a close signal and starts to close in preparation for the active short circuit operation.</p>
<p>When the coupling inductor <italic>L</italic>
<sub>2</sub> is inserted into the circuit, the current through <italic>L</italic>
<sub>2</sub> will experience a step from its initial value of zero. This sudden change in current will result in an extremely high overvoltage, which can potentially damage IGBTs <italic>T</italic>
<sub>1</sub>. Therefore, it is necessary to parallel MOA.1 with <italic>L</italic>
<sub>2</sub> to limit the amplitude of the overvoltage and dissipate its energy. The equivalent circuit is depicted in <xref ref-type="fig" rid="F8">Figure 8A</xref>. The voltage across <italic>L</italic>
<sub>2</sub> is equals the voltage <italic>u</italic>
<sub>moa.1</sub> of the MOA.1. Hence, the current through <italic>L</italic>
<sub>2</sub> can be calculated by,<disp-formula id="e13">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:msubsup>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Equivalent circuit of inductance-CLCB: <bold>(A)</bold> equivalent circuit of (t<sub>
<italic>1</italic>
</sub>&#x223C;t<sub>
<italic>2</italic>
</sub>). <bold>(B)</bold> Simplified equivalent circuit of (t<sub>
<italic>1</italic>
</sub>&#x223C;t<sub>
<italic>2</italic>
</sub>). <bold>(C)</bold> equivalent circuit of (t<sub>
<italic>2</italic>
</sub>&#x223C;t<sub>
<italic>3</italic>
</sub>).</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g008.tif"/>
</fig>
<p>In this interval, the equivalent circuit of the CLCB is <italic>L</italic>
<sub>1</sub> in series with a voltage source <italic>u</italic>
<sub>moa.1</sub>, as shown in <xref ref-type="fig" rid="F8">Figure 8B</xref>. The CLCB voltage in this interval is<disp-formula id="e14">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
<p>At <italic>t</italic>
<sub>2</sub>, two inductances <italic>L</italic>
<sub>1</sub> and <italic>L</italic>
<sub>2</sub> share an identical current. By using constant-flux-linkage theorem, <italic>i</italic>
<sub>L2</sub> (t<sub>2</sub>) can be determined<disp-formula id="e15">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
</p>
<p>Where <italic>m</italic> is the mutual coefficient, and <italic>i</italic>
<sub>L1</sub> (t<sub>1</sub>) is the current value through <italic>L</italic>
<sub>1</sub> before the IGBT <italic>T</italic>
<sub>1</sub> is turned on.</p>
</sec>
<sec id="s4-2-2">
<title>4.2.2 Inductive current limiting interval (<italic>t</italic>
<sub>
<italic>2</italic>
</sub>&#x223c;<italic>t</italic>
<sub>
<italic>3</italic>
</sub>)</title>
<p>Since <italic>L</italic>
<sub>1</sub> and <italic>L</italic>
<sub>2</sub> share identical current at <italic>t</italic>
<sub>2</sub>, no current passes through MOA.1, resulting in that MOA.1 opening. Therefore, the fault current only flows through <italic>L</italic>
<sub>1</sub> and <italic>L</italic>
<sub>2</sub> and the equivalent circuit of CLCB is shown in <xref ref-type="fig" rid="F8">Figure 8C</xref>, and <italic>L</italic>
<sub>CL</sub> is the decoupled equivalent inductance of <italic>L</italic>
<sub>1</sub> and <italic>L</italic>
<sub>2</sub>. <italic>L</italic>
<sub>CL</sub> and voltage <italic>u</italic>
<sub>CLCB</sub> (t) across the CLCB can be calculated.<disp-formula id="e16">
<mml:math id="m17">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
</p>
</sec>
<sec id="s4-2-3">
<title>4.2.3 Current shifting by active short-circuit (<italic>t</italic>
<sub>
<italic>3</italic>
</sub>&#x223c;<italic>t</italic>
<sub>
<italic>4</italic>
</sub>)</title>
<p>At <italic>t</italic>
<sub>3</sub>, UFD<sub>2</sub> has been closed completely, <italic>T</italic>
<sub>2</sub> is turned on and UFD<sub>1</sub> opens. As a result, the CLCB circuit forms a low impedance ground branch, consisting of <italic>L</italic>
<sub>2</sub>, UFD<sub>2</sub> and <italic>T</italic>
<sub>2</sub>. This creates a new current flow path, represented by the green dotted line in <xref ref-type="fig" rid="F7">Figure 7A</xref>. Hence, the equivalent circuit and voltage <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub>
<italic>(t)</italic> expression are identical to the inductive current limiting interval.</p>
</sec>
<sec id="s4-2-4">
<title>4.2.4 MOV current limiting interval (<italic>t</italic>
<sub>
<italic>4</italic>
</sub>&#x223c;<italic>t</italic>
<sub>
<italic>5</italic>
</sub>)</title>
<p>At <italic>t</italic>
<sub>4</sub>, <italic>T</italic>
<sub>2</sub> is turned off to force the fault current to pass through MOA.2. So that, <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub>
<italic>(t)</italic> can be calculated by<disp-formula id="e17">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
</p>
<p>Since the reference voltage <italic>U</italic>
<sub>
<italic>ref</italic>
</sub> of MOA.2 is approximately 1.5 times the DC voltage source <italic>U</italic>
<sub>dc</sub> of VSC, MOA.2 is utilized to provide different voltages that counteract the fault current (<xref ref-type="bibr" rid="B17">Mohammadi et al., 2021</xref>). At <italic>t</italic>
<sub>4</sub>
<italic>,</italic> the fault current will be reduced to zero.</p>
<p>Based on the above modelling process, the voltage <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub>(t) across CLCB can be written as a united form as the followings<disp-formula id="e18">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>Therefore, the external electrical characteristic of Inductance-CLCB can be modelled as a time-variable voltage source shown in <xref ref-type="fig" rid="F6">Figure 6A</xref>, and <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub> can be calculated by Eq. <xref ref-type="disp-formula" rid="e18">18</xref>.</p>
</sec>
</sec>
<sec id="s4-3">
<title>4.3 Modelling of sequential MOA-CLCB</title>
<p>
<xref ref-type="fig" rid="F7">Figure 7B</xref> shows the MOA-CLCB circuit, proposed by ABB (<xref ref-type="bibr" rid="B17">Mohammadi et al., 2021</xref>). It consists of two paths referred to as LCS and main breaker. Under normal conditions, the current is allowed to pass through the LCS path. However, while a short circuit fault occurs, all IGBTs in the main breaker path are turned on to force the fault current to pass through the main breaker, comprised of several identical modules paralleled with MOA. Once the main breaker establishes a conducting path, the UFD will be opened (<xref ref-type="bibr" rid="B17">Mohammadi et al., 2021</xref>). When the UFD complete open operation, all IGBTs in the main breaker path are turned off in the conventional switching strategy, and then the fault current is reduced to zero. However, this conventional switching strategy can result in a high voltage.</p>
<p>It is noted that the UFD is a special element with non-linear-time variable characteristics, which has two moving contactors. When the contactors begin to separate, the distance between the two contactors increases linearly with time (<xref ref-type="bibr" rid="B22">Skarby and Steiger, 2013</xref>; <xref ref-type="bibr" rid="B8">Hedayati and Jovcic, 2017</xref>). By utilizing this non-linear-time variable characteristic, a sequential switching strategy of IGBTs in the main breaker was proposed, as shown in <xref ref-type="fig" rid="F7">Figure 7B</xref>. This strategy aims to reduce the peak fault current, overvoltage, and fault clearance time, and is referred to as sequential MOA-CLCB (<xref ref-type="bibr" rid="B9">Hedayati and Jovcic, 2018</xref>; <xref ref-type="bibr" rid="B23">Song et al., 2019</xref>).</p>
<p>At <italic>t</italic>
<sub>0</sub>, a short circuit fault occurs, and the fault current flows through the LCS path. At <italic>t</italic>
<sub>
<italic>d</italic>
</sub>
<italic>,</italic> let all IGBTs (<italic>T</italic>
<sub>m1</sub>&#x223c;<italic>T</italic>
<sub>
<italic>mn</italic>
</sub>) in these sub-modules be turned on to form a new fault current path. During the interval [<italic>t</italic>
<sub>d</sub>
<italic>, t</italic>
<sub>1</sub>], the fault current commutates from the LCS path to the main breaker, and <italic>u</italic>
<sub>CLCB</sub> is about zero.</p>
<p>At <italic>t</italic>
<sub>
<italic>1</italic>
</sub>
<italic>,</italic> let IGBT <italic>T</italic>
<sub>m1</sub> be on to force the fault current to pass through MOA.1. Hence, <italic>u</italic>
<sub>CLCB</sub> is identical to the MOA.1 voltage, expressed as<disp-formula id="e19">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>
</p>
<p>At <italic>t</italic>
<sub>
<italic>2</italic>
</sub>, <italic>T</italic>
<sub>m2</sub> is turned on and <italic>u</italic>
<sub>CLCB</sub> (t) can be written as<disp-formula id="e20">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(20)</label>
</disp-formula>
</p>
<p>The rest can be done in the same manner. Hence, the voltage <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub> across the sequential MOA-CLCB is written as a canonical form, such as<disp-formula id="e21">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(21)</label>
</disp-formula>
</p>
</sec>
<sec id="s4-4">
<title>4.4 Modelling of capacitance-CLCB</title>
<p>
<xref ref-type="fig" rid="F7">Figure 7C</xref> shows the capacitance-CLCB circuit and its switching strategy [24]. At <italic>t</italic>
<sub>0</sub>, the MTDC grid experiences a short circuit fault. At <italic>t</italic>
<sub>
<italic>d</italic>
</sub>, the MDTC grid detects the short circuit fault and sends a break command to capacitance-CLCB. Firstly, all IGBTs <italic>T</italic>m in the main breaker branch is turned on to form the fault current transfer path. After that, the LCS branch receives the turn-off drive signal, and the UFD starts to open to commutate the fault current from the LCS branch to the main breaker. So, during the interval [<italic>t</italic>
<sub>0</sub>
<italic>, t</italic>
<sub>1</sub>], the voltage <italic>u</italic>
<sub>CLCB</sub> across the capacitance-CLCB is approximately equal to zero.</p>
<p>At <italic>t</italic>
<sub>1</sub>, all IGBTs <italic>T</italic>
<sub>m</sub> is turned off, so that the fault current is forced to transfer to the capacitance <italic>C</italic>
<sub>m</sub> through freewheeling diodes <italic>D</italic>
<sub>m</sub>. So that <italic>C</italic>
<sub>m</sub> begins to be charged by the fault current, and the voltage across the capacitances starts to increase from its initial value of zero. Hence, the voltage <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub> is.<disp-formula id="e22">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x222b;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mi>t</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>F</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(22)</label>
</disp-formula>
</p>
<p>Where <italic>N</italic> is the number of sub-modules in the main breaker.</p>
<p>At time <italic>t</italic>
<sub>2</sub>, the voltage of <italic>C</italic>
<sub>m</sub> is charged to be equal to the rated voltage <italic>U</italic>
<sub>r</sub> of MOA, and the part of the fault current will pass through the MOA. Hence, the voltage <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub>
<italic>(t)</italic> is<disp-formula id="e23">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x222b;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mi>t</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>F</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(23)</label>
</disp-formula>
</p>
<p>At <italic>t</italic>
<sub>3</sub>, the voltage of <italic>C</italic>
<sub>m</sub> is equal to the reference voltage <italic>U</italic>
<sub>ref</sub> of MOA, and MOA acts as a voltage source. Hence, the voltage <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub> is<disp-formula id="e24">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mtext>&#x2003;</mml:mtext>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(24)</label>
</disp-formula>
</p>
<p>Based on above modelling process, the voltage <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub>(t) can be written as a united form as the followings<disp-formula id="e25">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x222b;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mi>t</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>F</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x222b;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mi>t</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>F</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(25)</label>
</disp-formula>
</p>
<p>In this section, the models of three typical CLCBs have been established, and the main contributions are the following. A canonical voltage source model, shown in <xref ref-type="fig" rid="F6">Figure 6A</xref>, is proposed to describe the external electrical characteristic of the three typical CLCBs. The voltage source expressions are developed and shown in Eqs <xref ref-type="disp-formula" rid="e18">18</xref>, <xref ref-type="disp-formula" rid="e21">21</xref>, <xref ref-type="disp-formula" rid="e25">25</xref>.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Short circuit fault current calculation method of DC grid with CLCB</title>
<p>In this section, the four-terminal DC grid system shown in <xref ref-type="fig" rid="F1">Figure 1A</xref> will still be taken as an example to investigate the short circuit fault current calculation method of DC grid with CLCB (<xref ref-type="bibr" rid="B9">Hedayati and Jovcic, 2018</xref>; <xref ref-type="bibr" rid="B23">Song et al., 2019</xref>). <xref ref-type="fig" rid="F9">Figure 9A</xref> illustrates the four-terminal DC grid system with eight CLCBs. Under normal condition, the CLCB is considered as a short circuit, so that the DC grid with CLCB shown in <xref ref-type="fig" rid="F9">Figure 9A</xref> is identical with that of the DC grid shown in <xref ref-type="fig" rid="F1">Figure 1A</xref>. However, in the case of a short-circuit fault, the situation is totally different. For example, assume that a short-circuit fault occurs in the middle position of the power transmission line OL<sub>12</sub> between VSC<sub>1</sub> and VSC<sub>2</sub>, the CLCB<sub>a</sub> and CLCB<sub>b</sub> will respond, but the other CLCBs will remain silent.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>
<bold>(A)</bold> Four-terminal DC gird with CLCB, <bold>(B)</bold> its canonical RLC equivalent circuit.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g009.tif"/>
</fig>
<sec id="s5-1">
<title>5.1 Canonical RLC equivalent circuit of DC grid with CLCB</title>
<p>Although there are several types of CLCBs, a canonical voltage source model is established to describe their external electrical characteristics, as shown in <xref ref-type="fig" rid="F6">Figure 6A</xref>. The difference lies in the fact that each type of CLCB has its own voltage source expression.</p>
<p>For example, assuming a short-circuit fault occurs at the middle position of the overhead line OL<sub>12</sub> between VSC<sub>1</sub> and VSC<sub>2</sub>, the CLCB<sub>a</sub> and CLCB<sub>b</sub> will respond, but the other CLCBs still remain silent. By employing the aforementioned rule, a canonical RLC equivalent circuit of the DC grid with CLCB can be established and illustrated in <xref ref-type="fig" rid="F9">Figure 9B</xref>. Compared with <xref ref-type="fig" rid="F1">Figures 1A,B</xref> new voltage source, <italic>u</italic>
<sub>
<italic>CLCBa</italic>
</sub>, is inserted between node n<sub>1</sub> and the fault point n<sub>0</sub>. It expresses that the CLCB<sub>a</sub> connected with VSC1 has responded to the short circuit fault. And <italic>u</italic>
<sub>
<italic>CLCBb</italic>
</sub> is done in the same way as the CLCB<sub>a</sub>.</p>
</sec>
<sec id="s5-2">
<title>5.2 Normal form of differential equation set</title>
<p>The systematic formulation approach of the differential equation set developed in <xref ref-type="sec" rid="s2">Section 2</xref> is also suitable for writing the differential equation set for the canonical RLC equivalent circuit shown in <xref ref-type="fig" rid="F9">Figure 9B</xref>. The systematic formulation of the differential equation set is a universal approach that is independent of the topology of the DC grid.</p>
<p>In the analysis of the circuit shown in <xref ref-type="fig" rid="F9">Figure 9B</xref>, if we chose the loop current like that of <xref ref-type="fig" rid="F1">Figure 1B</xref>, then the loop currents rector <bold>
<italic>i</italic>
</bold>, the capacitance voltages rector <bold>
<italic>u</italic>
</bold> and its current rector <bold>
<italic>i</italic>
</bold>
<sub>
<italic>c</italic>
</sub> are identical, also expressed by Eq. <xref ref-type="disp-formula" rid="e1">1</xref>. Therefore, the incidence matrix <bold>A</bold> and <bold>R</bold> as well as <bold>L</bold> parameter matrix have an identical form. However, compared to <xref ref-type="fig" rid="F1">Figure 1B</xref>, two new voltage sources <italic>u</italic>
<sub>CLCBa</sub> and <italic>u</italic>
<sub>CLCBb</sub> are added in <xref ref-type="fig" rid="F9">Figure 9B</xref>. As a result, voltage source rector should be modified as<disp-formula id="e26">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(26)</label>
</disp-formula>where <italic>u</italic>
<sub>CLCB1</sub> &#x3d; <italic>u</italic>
<sub>CLCBa</sub> &#x2260; 0 and <italic>u</italic>
<sub>CLCB2</sub> &#x3d; <italic>u</italic>
<sub>CLCBb</sub> &#x2260; 0, the CLCB<sub>a</sub> and CLCB<sub>b</sub> have responded since a short fault occured on the overhand line OL<sub>12</sub>. Additionally, <italic>u</italic>
<sub>CLCB3</sub> &#x3d; <italic>u</italic>
<sub>CLCB4</sub> &#x3d; <italic>u</italic>
<sub>CLCB5</sub> &#x3d; 0, implying that the other CLCBs keep silent because no short fault has occurred on the overhead line OL<sub>13</sub>, OL<sub>24</sub> and OL<sub>34</sub> illustrated in <xref ref-type="fig" rid="F9">Figure 9A</xref>.</p>
<p>Hence, the differential equation set shown Eq. <xref ref-type="disp-formula" rid="e7">7</xref> can be modified as<disp-formula id="e27">
<mml:math id="m28">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold">A</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">R</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold">L</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold">P</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(27)</label>
</disp-formula>
</p>
<p>This formula is called the normal form of the differential equation set, which will be used to calculate the short circuit fault current of DC grid with CLCB in this paper.</p>
<p>The normal form of differential equation set shares the following several merits.<list list-type="simple">
<list-item>
<p>1) The normal form is more general because it is independent of the DC grid&#x2019;s topology and the SCPE categories.</p>
</list-item>
<list-item>
<p>2) The differential equation sets, and their parameters can be easily written down by employing the systematic formulation approach of the differential equation set developed in <xref ref-type="sec" rid="s2">Section 2</xref>. This allows us to avoid the large number of matrix computations required when writing equation sets directly using KVL and KCL.</p>
</list-item>
<list-item>
<p>3) By using MATLAB, the solution can be obtained quickly because the RLC equivalent circuit is a linear-non-time variable, and the voltage source <italic>u</italic>
<sub>CLCB</sub>(t), as shown in Eqs <xref ref-type="disp-formula" rid="e18">18</xref>, <xref ref-type="disp-formula" rid="e21">21</xref>, <xref ref-type="disp-formula" rid="e25">25</xref>, is a piecewise linear expression. So, the normal form is a piecewise linear differential equation set.</p>
</list-item>
</list>
</p>
</sec>
</sec>
<sec id="s6">
<title>6 Validation and results analysis</title>
<sec id="s6-1">
<title>6.1 Validation</title>
<p>This section will validate the short circuit fault current calculation method of the DC grid with CLCB proposed in <xref ref-type="sec" rid="s5">Section 5</xref> by employing the four-terminal DC grid shown in <xref ref-type="fig" rid="F9">Figure 9A</xref>. The equivalent model of pole-to-ground short-circuit fault is illustrated in <xref ref-type="fig" rid="F9">Figure 9B</xref>. The parameters of VSC are shown in <xref ref-type="table" rid="T1">Table 1</xref>. The equivalent inductance and equivalent resistance of the overhead line are 1&#xa0;mH/km and 0.01&#xa0;&#x3a9;/km respectively. The DC reactor at both ends of the overhead line is 150&#xa0;mH. About the parameters of the three typical CLCBs: The rated voltage of a single MOA module for all three typical CLCBs is 40&#xa0;kA. The sub-module capacitor of the Capacitance-CLCB is 240&#xa0;uF. The coupling inductance pair of the Inductance-CLCB is 100&#xa0;mH, and the coupling factor is 0.9.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Parameters of MMC-VSC</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">VSC station</th>
<th align="left">1</th>
<th align="left">2</th>
<th align="left">3</th>
<th align="left">4</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Rated DC voltage U<sub>dc</sub>
</td>
<td align="left">&#xb1;500&#xa0;kV</td>
<td align="left">&#xb1;500&#xa0;kV</td>
<td align="left">&#xb1;500&#xa0;kV</td>
<td align="left">&#xb1;500&#xa0;kV</td>
</tr>
<tr>
<td rowspan="2" align="left">Control strategy</td>
<td align="left">Q &#x3d; 0MVar</td>
<td align="left">Q &#x3d; 0MVar</td>
<td align="left">Q &#x3d; 0MVar</td>
<td align="left">Q &#x3d; 0MVar</td>
</tr>
<tr>
<td align="left">
<italic>p</italic> &#x3d; 750&#xa0;MW</td>
<td align="left">U<sub>dc</sub> &#x3d; 500&#xa0;kV</td>
<td align="left">
<italic>p</italic> &#x3d; 750&#xa0;MW</td>
<td align="left">
<italic>p</italic> &#x3d; &#x2212;1000&#xa0;MW</td>
</tr>
<tr>
<td align="left">N<sub>sm</sub>
</td>
<td align="left">200</td>
<td align="left">200</td>
<td align="left">200</td>
<td align="left">200</td>
</tr>
<tr>
<td align="left">R<sub>arm</sub> (m&#x3a9;)</td>
<td align="left">1.361</td>
<td align="left">1.361</td>
<td align="left">0.908</td>
<td align="left">1.361</td>
</tr>
<tr>
<td align="left">L<sub>arm</sub> (mH)</td>
<td align="left">100</td>
<td align="left">100</td>
<td align="left">50</td>
<td align="left">50</td>
</tr>
<tr>
<td align="left">C<sub>arm</sub> (&#x3bc;F)</td>
<td align="left">15000</td>
<td align="left">10000</td>
<td align="left">15000</td>
<td align="left">15000</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In order to accurately calculate the short circuit fault current, the four-terminal DC grid with three different types of CLCB shown in <xref ref-type="fig" rid="F9">Figure 9A</xref> is established in EMT simulation tool PSCAD. The simulation results are illustrated in <xref ref-type="fig" rid="F10">Figure 10</xref> by the solid lines and its VSC in this DC grid model is based on the EMT equivalent model of half-bridge MMC proposed in (<xref ref-type="bibr" rid="B6">CIGRE WG B4.57, 2014</xref>; <xref ref-type="bibr" rid="B7">Gnanarathna et al., 2011</xref>). This simulation tool is generally considered to have the ability to accurately estimate the performance of an HVDC grid in normal or fault conditions, making the EMT simulation results a reliable reference criterion. In short circuit fault current calculation, getting the fault current curves is the most concerned. Therefore, <xref ref-type="fig" rid="F10">Figure 10</xref> shows solely the fault current <italic>i</italic>
<sub>1F</sub> and <italic>i</italic>
<sub>2F</sub> curves. However, EMT simulation would require a lot of computing resources and be expensive and time-consuming.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Fault current results of MTDC grid with different typical CLCB from EMT simulation and proposed fast-computing method; <bold>(A)</bold> MTDC grid with Inductance-CLCB <bold>(B)</bold> MTDC grid with MOA-CLCB <bold>(C)</bold> MTDC grid with Capacitance-CLCB.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g010.tif"/>
</fig>
<p>In order to improve simulating efficiency, the short circuit fault current calculation equation is derived in <xref ref-type="sec" rid="s5">Section 5</xref> as expressed in Eq. <xref ref-type="disp-formula" rid="e27">27</xref>. By applying the systematic formulation proposed in <xref ref-type="sec" rid="s2">Section 2</xref>, the parameter matrices of the equivalent circuit shown in <xref ref-type="fig" rid="F9">Figure 9B</xref> can be obtained. The incidence matrix <bold>A</bold> is also shown in <xref ref-type="fig" rid="F2">Figure 2B</xref>, and Eq. <xref ref-type="disp-formula" rid="e6">6</xref> can be used to compute the <bold>P</bold> matrix. <bold>R</bold> and <bold>L</bold> matrices are<disp-formula id="equ1">
<mml:math id="m29">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold">R</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>10</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>20</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>24</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>34</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>13</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="bold">L</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>10</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>20</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>24</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>34</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mn>13</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
</disp-formula>
</p>
<p>Using the MATLAB program to solve the short circuit fault current calculation equation set Eq. <xref ref-type="disp-formula" rid="e27">27</xref>, we can obtain the fault current <italic>i</italic>
<sub>1F</sub> and <italic>i</italic>
<sub>2F</sub> curves illustrated in <xref ref-type="fig" rid="F10">Figure 10</xref> by the dotted lines. For the inductance-CLCB, the formula of Eq. <xref ref-type="disp-formula" rid="e18">18</xref> is utilized to express <italic>u</italic>
<sub>CLCBa</sub> and <italic>u</italic>
<sub>CLCBb</sub> in the voltage source rector of Eq. <xref ref-type="disp-formula" rid="e27">27</xref>, and the simulation results are shown in <xref ref-type="fig" rid="F10">Figure 10A</xref>. The Eq. <xref ref-type="disp-formula" rid="e21">21</xref> and Eq. <xref ref-type="disp-formula" rid="e25">25</xref> are used to express the voltage source rector of MOA-CLCB and capacitance-CLCB, respectively, and the simulation results are shown in <xref ref-type="fig" rid="F10">Figures 10B,C</xref>, respectively.</p>
<p>
<xref ref-type="fig" rid="F10">Figure 10</xref> shows four curves labelled &#x2460;, &#x2461;, &#x2462; and &#x2463;. The curve &#x2460; and &#x2461; are the fault current <italic>i</italic>
<sub>1F</sub> and <italic>i</italic>
<sub>2F</sub>, respectively. It can be observed that the simulated results from MATLAB, indicated by the dotted lines, almost agreed with that of conventional EMT simulation shown in the solid lines. The curves &#x2462; and &#x2463; are error curves, illustrating that the maximum computing error is less than 5%. Therefore, these results confirm that the proposed model and formulas are accurate enough to can meet the requirement of routine engineering analysis and design.</p>
<p>In order to verify that the accuracy of the proposed DC fault current fast-computing method is not affected by changes in test conditions, such as the location of fault occurrence, we selected a four-terminal DC network with an inductive CLCB as the test object, as shown in <xref ref-type="fig" rid="F9">Figure 9</xref>. Several rounds of the general EMT simulation and the proposed DC fault current fast-computing method were performed respectively on PSCAD and MATLAB under the same test conditions. However, the fault location was selected as a variable. The compared errors between the EMT simulating results and the proposed method calculating results are shown in <xref ref-type="fig" rid="F11">Figure 11</xref>. In this case, the 2-D graph shown in <xref ref-type="fig" rid="F10">Figure 10</xref> becomes a 3-D graph in <xref ref-type="fig" rid="F11">Figure 11</xref>.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>The Impact of fault location on the compared error of the proposed DC fault current fast-computing method: <bold>(A)</bold> fault current i1F <bold>(B)</bold> fault current i2F.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g011.tif"/>
</fig>
<p>In these 3D graphs, the <italic>x-</italic>axis indicates the relative location of the short-circuit fault point, denoted as n<sub>F</sub>, on the fault line. Let <italic>x</italic> be a ratio of <italic>L</italic>
<sub>fault</sub>/<italic>L</italic>
<sub>F-line</sub>, where <italic>L</italic>
<sub>fault</sub> represents the distance from the short circuit fault point <italic>n</italic>
<sub>F</sub> to VSC1, and <italic>L</italic>
<sub>F-line</sub> is the total length of the overhead line OL<sub>12</sub>. The <italic>y</italic>-axis shows the time, and the vertical axes <italic>z</italic> shows the compared error between the proposed method and general EMT. Specifically, <xref ref-type="fig" rid="F11">Figures 11A,B</xref> represent the errors fault current <italic>i</italic>
<sub>1F</sub> and <italic>i</italic>
<sub>2F</sub>, respectively. As shown in <xref ref-type="fig" rid="F11">Figure 11</xref>, it can be observed that the maximum computing error is less than 5%, regardless of the location of the short circuit fault occurs on the line. Therefore, these results further confirm the accuracy of the proposed model and formulas with good stability and universality.</p>
</sec>
<sec id="s6-2">
<title>6.2 Computing efficiency</title>
<p>The simulation platform used in this paper was an AMD Ryzen7 4800H 2.90&#xa0;GHz CPU with 16&#xa0;GB of RAM and a 64-bit Windows 10 Operating System. The time step is 10&#xa0;&#xb5;s. For the PSCAD EMT simulation platform, it takes 0.9&#xa0;s to start up. So, for the PSCAD EMT simulation and MATLAB calculation platform, let a dc short-circuit fault occurs at t &#x3d; 1.0&#xa0;s and fault time duration equals 10&#xa0;ms. The time-consuming and efficiency improvements are listed in <xref ref-type="table" rid="T2">Table 2</xref>. It can be seen that the proposed short circuit current calculation method is much more efficient than the conventional EMT simulation and has relatively higher accuracy. The computing efficiency has been improved at least three hundred times, and the accuracy meets engineering analysis and design requirements.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Time-consuming and efficiency improving.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left"/>
<th align="left">MATLAB (ms)</th>
<th align="left">PSCAD/EMTDCT (ms)</th>
<th align="left">Efficiency</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">MOA &#x2013;CLCB</td>
<td align="left">2218</td>
<td align="left">986406</td>
<td align="left">445</td>
</tr>
<tr>
<td align="left">Capacitance-CLCB</td>
<td align="left">2375</td>
<td align="left">987484</td>
<td align="left">415</td>
</tr>
<tr>
<td align="left">Inductance-CLCB</td>
<td align="left">2781</td>
<td align="left">940625</td>
<td align="left">338</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s7">
<title>7 Performance comparison of three -typical CLCBs</title>
<p>Using the conventional EMT simulation tool would take a long time to estimate the performance of a DC grid with SCPE since it is a complex dynamic system. However, the proposed short circuit current calculation method offers greater computing efficiency and relatively higher accuracy. It is a powerful tool for one to analyses the performances of DC grid with SCPE. For DC grid projects that require the installation of SCPE equipment, this method can help design and quickly screen and evaluate the topology, configuration scheme and timing logic of the SCPEs.</p>
<p>To demonstrate this function, we analyze and compare the current limiting capabilities of three typical CLCBs under varying initial fault current values and fault locations and draw some significant conclusions. This will be done using the method described at the end of this paper.</p>
<sec id="s7-1">
<title>7.1 Comparison of current limiting capability</title>
<p>At first, since ABB-DCCB is a typical hybrid DC circuit breaker and has been widely used in HVDC grids (<xref ref-type="bibr" rid="B22">Skarby and Steiger, 2013</xref>), it is applied to the four-terminal DC grid shown in <xref ref-type="fig" rid="F9">Figure 9A</xref>. The EMT platform is utilized as a simulation tool to estimate its performance. The short circuit fault current simulation results are shown in <xref ref-type="fig" rid="F12">Figure 12</xref> by dotted lines, labelled as &#x201c;i1F-DCCB&#x201d; and &#x201c;i2F-DCCB,&#x201d; which is used as a reference criterion in this section.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Comparison of current limiting capability of three CLCB <bold>(A)</bold> Inductance-CLCB <bold>(B)</bold> MOA-CLCB <bold>(C)</bold> Capacitance-CLCB.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g012.tif"/>
</fig>
<p>In order to quickly compare the performance of the three typical CLCBs with the reference criterion, the proposed short circuit current calculation method is used to predict the fault currents. The simulation results are illustrated in <xref ref-type="fig" rid="F12">Figures 12A&#x2013;C</xref> by the solid lines, respectively.</p>
<p>According to <xref ref-type="fig" rid="F12">Figure 12</xref>, it can be seen that the three typical CLCBs have a better fault current limiting ability than the conventional ABB-DCCB, and MOA-CLCB demonstrates the best performance. However, the other CLCBs take a longer fault isolation time than the conventional ABB-DCCB.</p>
</sec>
<sec id="s7-2">
<title>7.2 Impact of the normal current and fault location on CLCB performance</title>
<p>In the four-terminal DC network shown in <xref ref-type="fig" rid="F9">Figure 9</xref>, the fault location <italic>n</italic>
<sub>F</sub> and normal current value are important parameters that will profoundly impact CLCB performance. Using the proposed short-circuit current calculation method, the performances of the three typical CLCBs can be estimated, and the simulation results are illustrated in <xref ref-type="fig" rid="F13">Figure 13</xref>, which are three-dimensional graphs (3D graphs). The short-circuit fault current peak amplitude and fault isolation time of ABB-DCCB are also taken as standard values to evaluate the performance of the three typical CLCBs, such as peak fault current reduction value, peak current limiting ratio, and fault isolating time.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>3D performance graph of three typical CLCB; <bold>(A&#x2013;C)</bold> The Impact of the normal current and fault location on the performance of MOA-CLCB. <bold>(D&#x2013;F)</bold>The Impact of the normal current and fault location on the performance of capacitance-CLCB. <bold>(G&#x2013;I)</bold> The Impact of the normal current and fault location on the performance of inductance-CLCB.</p>
</caption>
<graphic xlink:href="fenrg-12-1366283-g013.tif"/>
</fig>
<p>In these 3D graphs, the <italic>x-</italic>axis indicates the relative location of the short-circuit fault point <italic>n</italic>
<sub>F</sub> in the fault line. The <italic>y</italic>-axis shows the normal current, and the vertical axes <italic>z</italic> are the relative values of fault current peak amplitude, peak fault current limiting ratio and relative values of fault isolating time, respectively.</p>
<p>
<xref ref-type="fig" rid="F13">Figures 13A&#x2013;C</xref> are the relative value of fault current peak amplitude, peak fault current limiting ratio and fault isolating time of the DC grid with MOA-CLCB, respectively. From <xref ref-type="fig" rid="F13">Figure 13A</xref>, it can be observed that the relative value of the peak current amplitude depends on the fault point location and is almost independent of the normal value. However, the peak current limiting rate shown in <xref ref-type="fig" rid="F13">Figure 13B</xref> is correlated with the fault point location and the normal value. As shown in <xref ref-type="fig" rid="F13">Figure 13C</xref>, the isolating fault time of the MOA-CLCB is always shorter than that of ABB-DCCB. Consequently, MOA-CLCB demonstrates the best robustness performance among the three typical CLCBs.</p>
<p>
<xref ref-type="fig" rid="F13">Figures 13D&#x2013;F</xref> demonstrate the performance of the capacitance-CLCB. The performance depends on both the short-circuit fault location and the normal value, but the short-circuit fault location has a more significant impact on the performance of capacitance-CLCB. When x &#x3e; 0.8 and y &#x3c; 0.5&#xa0;kA, capacitance-CLCB will lose its current limiting ability.</p>
<p>
<xref ref-type="fig" rid="F13">Figures 13G&#x2013;I</xref> show the performance of the inductance-CLCB, which is similar to that of the capacitance-CLCB. However, its current limiting effect is better than that of the capacitance-CLCB. It can be observed that there exists an optimal current limiting area, x &#x3c; 0.3. In this area, the peak current limit rate exceeds 40%, and the relative value of peak current is greater than 5&#xa0;kA. So, its current limiting effect is far better than the other two CLCBs in the optimal current limiting area. However, it is important to note that the inductance-CLCB achieves its best current limiting effect at the expense of a longer fault isolating time.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s8">
<title>8 Conclusion</title>
<p>The proposed DC fault current fast-computing method is much more efficient than the conventional EMT simulation and has a relative higher accuracy. When compared to the conventional time-domain numerical DC fault current calculation method of MMC-HVDC grid, the proposed DC fault current fast-computing method offers three improvements.</p>
<p>Firstly, by applying modern electrical circuit theory, the proposed systematic formulation makes it easy to write the normal form of the differential equation set instead of writing the equation set using KVL, KCL, and VAR. This approach avoids complex and cumbersome manual matrix calculations, making the proposed calculation method suitable for handling large-scale MMC- HVDC grids.</p>
<p>Secondly, the proposed Y-Delta transformation in the s-domain can eliminate virtual nodes in short circuit fault RLC equivalent circuit of the complex MMC-HVDC grid. This progress makes the proposed calculation method suitable for handling complex structure MMC-HVDC grids.</p>
<p>Finally, to make SCPE become a linear circuit, a canonical voltage source model of SCPE is proposed. In this paper, three types of CLCB are taken as examples to introduce how to establish the canonical voltage source model of SCPE. This progress makes the proposed calculation method suitable for handling MMC-HVDC grids with SCPE and expands its application field.</p>
<p>To comparing with the conventional PSCAD/EMTDC, the computing efficiency is improved at least about three hundreds time and the accuracy can meet requirement of engineering analysis and design. In summary, the proposed method can be applied to large-scale MMC-HVDC grids with complex structures and provides a powerful tool for analysing fault currents and evaluating the performance of short-circuit protection devices.</p>
<p>It is first time in this paper to estimate and compare the performances of three typical CLCBs, obtaining some significant conclusions as the followings. 1). MOA-CLCB shows the best robustness performance in three typical CLCBs. 2). There is a lost current limiting area for capacitance-CLCB. 3). For the inductance-CLCB, there is an optimal current limiting area, in where its current limiting effect is far better than the other two CLCBs, but it is at expense of a longer fault isolating time.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s9">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s10">
<title>Author contributions</title>
<p>XZ: Conceptualization, Data curation, Formal Analysis, Visualization, Writing&#x2013;original draft, Writing&#x2013;review and editing. XY: Project administration, Supervision, Writing&#x2013;review and editing. CZ: Resources, Writing&#x2013;review and editing, Investigation.</p>
</sec>
<sec sec-type="funding-information" id="s11">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.</p>
</sec>
<sec sec-type="COI-statement" id="s12">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s13">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s14">
<title>Abbreviations</title>
<p>A, incidence matrix; <italic>C</italic>, capacitance; <italic>C</italic>
<sub>m</sub>, current limiting capacitance; <italic>D</italic>
<sub>m</sub>, freewheeling diode; <italic>i</italic>, current; <italic>i</italic>
<sub>
<italic>c</italic>
</sub>, capacitance current; <italic>i</italic>
<sub>&#x394;</sub>, current of Delta equivalent circuit; <italic>L</italic>, inductance; <italic>L</italic>
<sub>&#x394;</sub>, inductance of Delta equivalent circuit; MOA, metal-oxide arresters; <italic>n</italic>
<sub>0</sub>, fault point; <italic>R</italic>, resistance; <italic>R</italic>
<sub>&#x2206;</sub>, resistance of Delta equivalent circuit; <italic>T</italic>, IGBTs; <italic>T</italic>m, IGBTs of current-limiting circuit; <italic>U</italic>
<sub>r</sub>, rated voltage of MOA; UFD, ultra-fast disconnectors; <italic>U</italic>
<sub>ref</sub>, reference voltage of MOA; <italic>u</italic>, voltage; <italic>u</italic>
<sub>moa</sub>, voltage of MOA; <italic>u</italic>
<sub>
<italic>CLCB</italic>
</sub>, voltage of CLCB; <italic>u</italic>
<sub>c</sub>, capacitance voltage; VSC, voltage-source converter.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ahmad</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gong</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Nadeem</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A hybrid circuit breaker with fault current limiter circuit in a VSC-HVDC application</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>7</volume>, <fpage>43</fpage>. <pub-id pub-id-type="doi">10.1186/s41601-022-00264-9</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Alexander</surname>
<given-names>C. K.</given-names>
</name>
<name>
<surname>Sadiku</surname>
<given-names>M. N. O.</given-names>
</name>
</person-group> (<year>2013</year>). <source>Fundamentals of electric circuits</source>. <publisher-loc>New York, New York, United States</publisher-loc>: <publisher-name>McGraw-Hill</publisher-name>.</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>An</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Research and application on multi-terminal and DC grids based on VSC-HVDC technology in China</article-title>. <source>High. Volt.</source> <volume>2</volume>, <fpage>1</fpage>&#x2013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1049/hve.2017.0010</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Balabanian</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Bickart</surname>
<given-names>T. A.</given-names>
</name>
</person-group> (<year>1969</year>). <source>Electrical network theory</source>. <publisher-loc>Hoboken, New Jersey, United States</publisher-loc>: <publisher-name>John Wiley &#x26; Sons</publisher-name>.</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Belda</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Plet</surname>
<given-names>C. A.</given-names>
</name>
<name>
<surname>Smeets</surname>
<given-names>R. P. P.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Analysis of faults in multiterminal HVDC grid for definition of test requirements of HVDC circuit breakers</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>33</volume>, <fpage>403</fpage>&#x2013;<lpage>411</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2017.2716369</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="web">
<collab>CIGRE WG B4.57</collab> (<year>2014</year>). <article-title>Guide for the development of models for hvdc converters in a hvdc grid</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://e-cigre.org/publication/604-guide-for-the-development-of-models-for-hvdc-converters-in-a-hvdc-grid">https://e-cigre.org/publication/604-guide-for-the-development-of-models-for-hvdc-converters-in-a-hvdc-grid</ext-link>
</comment> (<comment>Accessed March 30, 2024</comment>).</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gnanarathna</surname>
<given-names>U. N.</given-names>
</name>
<name>
<surname>Gole</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Jayasinghe</surname>
<given-names>R. P.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Efficient modeling of modular multilevel HVDC converters (MMC) on electromagnetic transient simulation programs</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>26</volume>, <fpage>316</fpage>&#x2013;<lpage>324</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2010.2060737</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hedayati</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Jovcic</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Low voltage prototype design, fabrication and testing of ultra-fast disconnector (UFD) for hybrid DC CB</article-title>. <source>CIGRE Winn. 2017 Colloq.</source>,</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hedayati</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Jovcic</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Reducing peak current and energy dissipation in hybrid HVDC CBs using disconnector voltage control</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>33</volume>, <fpage>2030</fpage>&#x2013;<lpage>2038</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2018.2812713</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Heidary</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Bigdeli</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Rouzbehi</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Controllable reactor based hybrid HVDC breaker</article-title>. <source>High. Volt.</source> <volume>5</volume>, <fpage>543</fpage>&#x2013;<lpage>548</lpage>. <pub-id pub-id-type="doi">10.1049/hve.2019.0354</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Transient fault identification method for bipolar short-circuit fault on MMC-HVDC overhead lines based on hybrid HVDC breaker</article-title>. <source>High. Volt.</source> <volume>6</volume>, <fpage>881</fpage>&#x2013;<lpage>893</lpage>. <pub-id pub-id-type="doi">10.1049/hve2.12091</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2017b</year>). <article-title>A Pole-to-Pole short-circuit Fault Current calculation method for DC grids</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>32</volume>, <fpage>4943</fpage>&#x2013;<lpage>4953</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2017.2682110</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Holliday</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Page</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Finney</surname>
<given-names>S. J.</given-names>
</name>
<name>
<surname>Williams</surname>
<given-names>B. W.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Continuous operation of radial multiterminal HVDC systems under DC fault</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>31</volume>, <fpage>351</fpage>&#x2013;<lpage>361</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2015.2471089</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>DC fault current approximation and fault clearing methods for hybrid LCC-VSC HVDC networks</article-title>. <source>Int. J. Electr. Power &#x26; Energy Syst.</source> <volume>143</volume>, <fpage>108467</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2022.108467</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2017a</year>). <article-title>DC fault analysis for modular multilevel converter-based system</article-title>. <source>J. Mod. Power Syst. Clean Energy</source> <volume>5</volume>, <fpage>275</fpage>&#x2013;<lpage>282</lpage>. <pub-id pub-id-type="doi">10.1007/s40565-015-0174-3</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Martinez</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Durbak</surname>
<given-names>D. W.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Parameter determination for modeling systems transients&#x2014;Part V: surge arresters IEEE PES task force on data for modeling system transients of IEEE PES working group on modeling and analysis of system transients using digital simulation (general systems subcommittee)</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>20</volume>, <fpage>2073</fpage>&#x2013;<lpage>2078</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2005.848771</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mohammadi</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Rouzbehi</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Hajian</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Niayesh</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Gharehpetian</surname>
<given-names>G. B.</given-names>
</name>
<name>
<surname>Saad</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>HVDC circuit breakers: a comprehensive review</article-title>. <source>IEEE Trans. Power Electron.</source> <volume>36</volume>, <fpage>13726</fpage>&#x2013;<lpage>13739</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2021.3073895</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ning</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Xiang</surname>
<given-names>C. U. I.</given-names>
</name>
<name>
<surname>Jiangjiang</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>Q. I.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Calculation method for line-to-ground short-circuit currents of VSC-HVDC grid</article-title>. <source>Electr. POWER Constr.</source> <volume>40</volume>, <fpage>119</fpage>&#x2013;<lpage>127</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-7229.2019.04.014</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Saad</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Peralta</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Denneti&#xe8;re</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Mahseredjian</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Jatskevich</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Martinez</surname>
<given-names>J. A.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>Dynamic averaged and simplified models for MMC-based HVDC transmission systems</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>28</volume>, <fpage>1723</fpage>&#x2013;<lpage>1730</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2013.2251912</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Safaei</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zolfaghari</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gilvanejad</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gharehpetian</surname>
<given-names>G. B.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>A survey on fault current limiters: development and technical aspects</article-title>. <source>Int. J. Electr. Power &#x26; Energy Syst.</source> <volume>118</volume>, <fpage>105729</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2019.105729</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Single-ended protection method for hybrid HVDC transmission line based on transient voltage characteristic frequency band</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>8</volume>, <fpage>26</fpage>. <pub-id pub-id-type="doi">10.1186/s41601-023-00301-1</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Skarby</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Steiger</surname>
<given-names>U.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>An ultra-fast disconnecting switch for a hybrid HVDC breaker&#x2014;a technical breakthrough</article-title>. <source>Proc. CIGR&#xc9; Sess.</source>, <fpage>1</fpage>&#x2013;<lpage>9</lpage>.</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Saeedifard</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Meliopoulos</surname>
<given-names>A. P. S.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Reducing the fault-transient magnitudes in multiterminal HVdc grids by sequential tripping of hybrid circuit breaker modules</article-title>. <source>IEEE Trans. Industrial Electron.</source> <volume>66</volume>, <fpage>7290</fpage>&#x2013;<lpage>7299</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2018.2881941</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stepanov</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Mahseredjian</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Karaagac</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>Saad</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Adaptive modular multilevel converter model for electromagnetic transient simulations</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>36</volume>, <fpage>803</fpage>&#x2013;<lpage>813</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2020.2993502</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Xin</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>An inductively coupled HVDC current limiting circuit breaker</article-title>. <source>Proc. CSEE</source>, <fpage>1731</fpage>&#x2013;<lpage>1741</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.190980</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Rong</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Technical assessment on self-charging mechanical HVDC circuit breaker</article-title>. <source>IEEE Trans. Industrial Electron.</source> <volume>69</volume>, <fpage>3622</fpage>&#x2013;<lpage>3630</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2021.3070523</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Rong</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Yi</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Damping HVDC circuit breaker with current commutation and limiting integrated</article-title>. <source>IEEE Trans. Industrial Electron.</source> <volume>67</volume>, <fpage>10433</fpage>&#x2013;<lpage>10441</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2019.2962475</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>The enhanced DC Fault Current calculation method of MMC-HVDC grid with FCLs</article-title>. <source>IEEE J. Emerg. Sel. Top. Power Electron.</source> <volume>7</volume>, <fpage>1758</fpage>&#x2013;<lpage>1767</lpage>. <pub-id pub-id-type="doi">10.1109/JESTPE.2018.2888931</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Comprehensive understanding of DC pole-to-pole fault and its protection for modular multilevel converters</article-title>. <source>High. Volt.</source> <volume>3</volume>, <fpage>246</fpage>&#x2013;<lpage>254</lpage>. <pub-id pub-id-type="doi">10.1049/hve.2018.5045</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ye</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Efficient estimation and characteristic analysis of short-circuit currents for MMC-MTDC grids</article-title>. <source>IEEE Trans. Industrial Electron.</source> <volume>68</volume>, <fpage>258</fpage>&#x2013;<lpage>269</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2020.2965433</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhan</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Insulation structure design and electric field simulation of 500 kV isolation energy supply transformer for HVDC breaker</article-title>. <source>High. Volt.</source> <volume>7</volume>, <fpage>185</fpage>&#x2013;<lpage>196</lpage>. <pub-id pub-id-type="doi">10.1049/hve2.12102</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhuo</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A natural commutation current topology of hybrid HVDC circuit breaker integrated with limiting fault current</article-title>. <source>IET Generation, Transm. Distribution</source> <volume>17</volume>, <fpage>1509</fpage>&#x2013;<lpage>1524</lpage>. <pub-id pub-id-type="doi">10.1049/gtd2.12760</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>