<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">EE</journal-id>
<journal-id journal-id-type="nlm-ta">EE</journal-id>
<journal-id journal-id-type="publisher-id">EE</journal-id>
<journal-title-group>
<journal-title>Energy Engineering</journal-title>
</journal-title-group>
<issn pub-type="epub">1546-0118</issn>
<issn pub-type="ppub">0199-8595</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">28965</article-id>
<article-id pub-id-type="doi">10.32604/ee.2023.028965</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Model Predictive Control Strategy of Multi-Port Interline DC Power Flow Controller</article-title>
<alt-title alt-title-type="left-running-head">Model Predictive Control Strategy of Multi-Port Interline DC Power Flow Controller</alt-title>
<alt-title alt-title-type="right-running-head">Model Predictive Control Strategy of Multi-Port Interline DC Power Flow Controller</alt-title>
</title-group>
<contrib-group>
<contrib id="author-1" contrib-type="author">
<name name-style="western"><surname>Wang</surname><given-names>He</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western"><surname>Xu</surname><given-names>Xiangsheng</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-3" contrib-type="author">
<name name-style="western"><surname>Shen</surname><given-names>Guanye</given-names></name><xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-4" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Jing</surname><given-names>Bian</given-names></name><xref ref-type="aff" rid="aff-1">1</xref><email>bj_jjj@163.com</email></contrib>
<aff id="aff-1"><label>1</label><institution>Key Laboratory of Modern Power System Simulation and Control &#x0026; Renewable Energy Technology, Ministry of Education (Northeast Electric Power University)</institution>, <addr-line>Jilin, 132012</addr-line>, <country>China</country></aff>
<aff id="aff-2"><label>2</label><institution>Marketing Department, State Grid Jilin Electric Power Co., Ltd.</institution>, <addr-line>Changchun, 130061</addr-line>, <country>China</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Bian Jing. Email: <email>bj_jjj@163.com</email></corresp>
</author-notes>
<pub-date date-type="collection" publication-format="electronic">
<year>2023</year></pub-date>
<pub-date date-type="pub" publication-format="electronic"><day>28</day>
<month>9</month>
<year>2023</year></pub-date>
<volume>120</volume>
<issue>10</issue>
<fpage>2251</fpage>
<lpage>2272</lpage>
<history>
<date date-type="received">
<day>19</day>
<month>1</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>05</day>
<month>5</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2023 Wang et al.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Wang et al.</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_EE_28965.pdf"></self-uri>
<abstract>
<p>There are issues with flexible DC transmission system such as a lack of control freedom over power flow. In order to tackle these issues, a DC power flow controller (DCPFC) is incorporated into a multi-terminal, flexible DC power grid. In recent years, a multi-port DC power flow controller based on a modular multi-level converter has become a focal point of research due to its simple structure and robust scalability. This work proposes a model predictive control (MPC) strategy for multi-port interline DC power flow controllers in order to improve their steady-state dynamic performance. Initially, the mathematical model of a multi-terminal DC power grid with a multi-port interline DC power flow controller is developed, and the relationship between each regulated variable and control variable is determined; The power flow controller is then discretized, and the cost function and weight factor are built with numerous control objectives. Sub module sorting method and nearest level approximation modulation regulate the power flow controller; Lastly, the MATLAB/Simulink simulation platform is used to verify the correctness of the established mathematical model and the control performance of the suggested MPC strategy. Finally, it is demonstrated that the control strategy possesses the benefits of robust dynamic performance, multi-objective control, and a simple structure.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>DC power flow controller</kwd>
<kwd>model predictive control</kwd>
<kwd>modular multi-level converter</kwd>
<kwd>control strategy</kwd>
<kwd>dynamic performance</kwd>
</kwd-group>
<funding-group>
<award-group id="awg1">
<funding-source>National Natural Science Foundation of China</funding-source>
<award-id>52177074</award-id>
</award-group>
</funding-group>
</article-meta>
</front>
<body>
<sec id="s1">
<label>1</label>
<title>Introduction</title>
<p>Flexible DC transmission technology with the advantages of independent control of active and reactive power and no phase change failure has been recognized as an effective solution to the problem of renewable energy grid connection and electricity consumption [<xref ref-type="bibr" rid="ref-1">1</xref>,<xref ref-type="bibr" rid="ref-2">2</xref>] due to the rapid development of new energy generation. Flexible DC grids have increased grid redundancy and enhanced transmission reliability than a standard DC grid. However, it can pose challenges for the power system, such as challenging current regulation and sophisticated fault protection. If there are more than N-1 lines in a DC grid with N converter stations, the converter stations cannot manage the currents on all lines individually. In order to improve the flexibility of DC grid power flow control, a further type of power electronics-based device, the DCPFC [<xref ref-type="bibr" rid="ref-3">3</xref>], is introduced into the DC grid.</p>
<p>Among the many DCPFC topologies, interline DCPFC (interline direct current power flow controller, IDCPFC) does not require additional power supply and has a greater advantage in terms of engineering feasibility and technical economy. The first IDCPFC topology with dual H-bridges back-to-back was proposed in reference [<xref ref-type="bibr" rid="ref-4">4</xref>], and the loss of the dual H-bridge structure is about 2.4% of that of the varactor topology at 1 kA DC current. Compared with the dual H-bridge interline DC current controller, the single capacitor interline DC current controller proposed in reference [<xref ref-type="bibr" rid="ref-5">5</xref>] was able to save 33% of power electronics in the unidirectional current condition, but still does not solve the current ripple problem caused by individual capacitor charging and discharging. Reference [<xref ref-type="bibr" rid="ref-6">6</xref>] also proposes a power electronics saving DCPFC topology and designed an hysteresis band current control (HCC) based control system for its proposed DCPFC. Reference [<xref ref-type="bibr" rid="ref-7">7</xref>] extended the dual H-bridge interline DC current controller to multiple transmission lines with fewer power electronics used, but the interline DC current controller with a single capacitor has very limited expansion capability. Reference [<xref ref-type="bibr" rid="ref-8">8</xref>] compared different power flow control methods and uses a dual H-bridge back-to-back IDCPFC topology in reference [<xref ref-type="bibr" rid="ref-4">4</xref>] to study the role of DC power flow controller, where power flow current control is achieved by connecting a single capacitor to 2 lines at high frequency.</p>
<p>In recent years, DCPFCs based on modular multi-level converters (MMC) have attracted a great deal of interest. MMCs offer numerous advantages, including adaptability, controllability, superior power quality, and high power transmission capacity. Numerous writers [<xref ref-type="bibr" rid="ref-9">9</xref>&#x2013;<xref ref-type="bibr" rid="ref-12">12</xref>] have proposed their very own MMC-based DCPFCs by combining different MMC characteristics. By adjusting the number and direction of sub-module inputs in a single bridge arm, which has the benefits of a simple structure and easy extension into multi-ports, these DCPFCs make use of the operating characteristics of MMC sub-modules to meet the requirements of power quality and control effects.</p>
<p>Concerning the DCPFC control strategy, the bulk of previously published articles [<xref ref-type="bibr" rid="ref-13">13</xref>&#x2013;<xref ref-type="bibr" rid="ref-15">15</xref>] have established control strategies for their proposed DCPFC topologies; however, they are all standard linear control strategies based on the proportional-integral (PI) control technique. Using the PI control methodology, reference [<xref ref-type="bibr" rid="ref-16">16</xref>] developed and analyzes the mathematical model of a three-terminal DC grid with a DCPFC placed between two lines. Reference [<xref ref-type="bibr" rid="ref-17">17</xref>] proposed an expandable N-interline DCPFC and a control method based on proportional-integral-derivative (PID) control. Other publications have created control strategies for DCPFCs operating under specific conditions. Specifically, reference [<xref ref-type="bibr" rid="ref-18">18</xref>] examined a successful control strategy and the damping control mechanism of the DCPFC under sub-synchronous oscillations. The problem of a shift in the control degree of freedom caused by a non-transient disconnection defect in the DC ring grid is addressed in the study [<xref ref-type="bibr" rid="ref-19">19</xref>] by the creation of a control degree switching method for DCPFC. Reference [<xref ref-type="bibr" rid="ref-20">20</xref>] provided a control strategy for coordinated management of DCPFC and DC circuit breakers in order to reduce fault currents during DC faults.</p>
<p>Numerous investigations on DCPFC topologies and control strategies have been conducted in recent years. MPC is a nonlinear predictive control system that predicts controlled quantities and completes control using a predictive model of future states and online cyclic optimization calculations. In the area of process control, it has been successfully used in [<xref ref-type="bibr" rid="ref-21">21</xref>,<xref ref-type="bibr" rid="ref-22">22</xref>]. Power electronics and transmission control have seen a lot of study on finite control set model predictive control (FCS-MPC) over the past ten years. For the first time, FCS-MPC is utilized in [<xref ref-type="bibr" rid="ref-23">23</xref>] to control the MMC-HVDC system. MPC is also used in [<xref ref-type="bibr" rid="ref-24">24</xref>] to regulate the MMC&#x2019;s AC output current, and loop control is added to complete multi-objective control. In order to lower the computational cost of the MPC and increase the computational speed, a voltage sequencing algorithm is added to the control strategy [<xref ref-type="bibr" rid="ref-25">25</xref>]. To further reduce the computational cost when utilizing the voltage sequencing technique, reference [<xref ref-type="bibr" rid="ref-26">26</xref>] substituted a finite control set of output level numbers for the conventional effective switching state finite control set. In reference [<xref ref-type="bibr" rid="ref-27">27</xref>], to guarantee the output voltage&#x2019;s consistency, the number of bridge arm input modules is only allowed to fluctuate within a small range throughout the subsequent sampling time. This prevents the output voltage from changing too quickly.</p>
<p>In conclusion, DCPFC&#x2019;s steady-state control strategy is mostly built on the typical linear control strategy based on PI control, while MPC is not yet implemented. In order to improve the dynamic performance of the MMC-based DCPFC, a MPC strategy is suggested for the novel multi-port interline power flow controller (MI-PFC), established in reference [<xref ref-type="bibr" rid="ref-28">28</xref>], in this paper. On the basis of the original topology and power flow control technique, a mathematical model of the multi-terminal DC grid comprising the MI-PFC is developed, and a MPC strategy for the MI-PFC with optimum steady-state operating performance is subsequently devised. In order to test the performance of the control strategy, the MI-PFC electromagnetic transient model is constructed using the MPC strategy and the MATLAB/Simulink simulation tool.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>MI-PFC Topology</title>
<p><xref ref-type="fig" rid="fig-1">Fig. 1</xref> depicts the topology of the MI-PFC employed in the article. Full-bridge submodule forms the foundation for the MMC bridge arm. The primary components of the DC current controller are the equivalent resistor R<sub>1</sub>, capacitor C<sub>1</sub>, bridge arm reactor L<sub>1</sub>, and DC reactor L<sub>dc</sub>. The AC and DC routes of the MI-PFC, which consist of capacitor C<sub>1</sub> and DC reactor L<sub>dc</sub>, keep the AC and DC currents separate. The MMC bridge arm can be utilized as a modular device and is therefore easily expandable, allowing the number of sub-modules and bridge arms to be altered as needed to fulfill power quality specifications.</p>
<fig id="fig-1">
<label>Figure 1</label>
<caption>
<title>Topology of MI-PFC</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-1.tif"/>
</fig>
</sec>
<sec id="s3">
<label>3</label>
<title>MI-PFC Working Principle and Mathematical Model</title>
<p>MPC requires a discrete mathematical model to foresee future states. Consequently, it must be designed with MI-PFC-specific architecture and actual operating conditions in mind. After introducing the topology of MI-PFC in the preceding section, this section establishes the dynamic mathematical model of MI-PFC by analyzing its operating principle, then provides the discrete mathematical model required in MPC by discretizing it, and finally derives the necessary equation for prediction.</p>
<p>In a four-terminal DC grid, the three-port MI-PFC is installed at the site of converter station 1 as illustrated in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>, where <italic>I</italic><sub><italic>ij</italic></sub> represents the line current, <italic>R</italic><sub><italic>ij</italic></sub> represents the DC line resistance, <italic>U</italic><sub><italic>i</italic></sub> represents the voltage at the converter station, and <italic>I</italic><sub><italic>i</italic></sub> represents the outlet current.</p>
<fig id="fig-2">
<label>Figure 2</label>
<caption>
<title>Equivalent circuit diagram of power flow control</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-2.tif"/>
</fig>
<p>The interline DC current controller must have both the basic function of controlling current and the ability to balance its own power. Consequently, the MI-PFC bridge arm voltage is given by <xref ref-type="disp-formula" rid="eqn-1">Eq. (1)</xref>.</p>
<p><disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn><mml:mrow><mml:mtext>DC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn><mml:mrow><mml:mtext>AC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mtext>DC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mtext>AC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn><mml:mrow><mml:mtext>DC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn><mml:mrow><mml:mtext>AC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B8;</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>U</italic><sub>M<italic>i</italic>DC</sub> is the amplitude of the DC component of the voltage at both ends of the MMC and <italic>U</italic><sub>M<italic>i</italic>AC</sub> is the amplitude of the AC component of the voltage at both ends of the MMC; the former is responsible for controlling the grid current and the latter for balancing its own power; <italic>&#x03B8;</italic><sub><italic>i</italic></sub> is the initial phase angle of the AC component (<italic>i</italic> &#x003D; 1&#x007E;3).</p>
<p>A few kilovolts of voltage drop in a DC grid line is negligible compared to typical system-level values of several hundred kilovolts. Consequently, a minimal number of <italic>U</italic><sub>M1DC</sub>&#x2013;<italic>U</italic><sub>M3DC</sub> submodules of the MMC full-bridge variety can perform the role of controlling DC grid currents. The DC grid depicted in <xref ref-type="fig" rid="fig-2">Fig. 2</xref> is controlled by master-slave control, with converter station 1 adopting a constant voltage control strategy and the other three converter stations using a constant power control strategy.</p>
<p>Therefore, the currents <italic>I</italic><sub>2</sub>&#x007E;<italic>I</italic><sub>4</sub> at the converter station can be considered known quantities. <xref ref-type="fig" rid="fig-2">Fig. 2</xref> illustrates the relationship between DC voltage and line current according to Kirchhoff&#x2019;s voltage law.</p>
<p><disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn><mml:mrow><mml:mtext>DC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mtext>DC</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>34</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mn>34</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mtext>DC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn><mml:mrow><mml:mtext>DC</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>34</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>34</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>I</italic><sub>12</sub>, <italic>I</italic><sub>13</sub>, <italic>I</italic><sub>14</sub>, <italic>I</italic><sub>23</sub>, and <italic>I</italic><sub>34</sub> are the line currents as shown in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>; <italic>I</italic><sub>M1</sub>, <italic>I</italic><sub>M2</sub>, and <italic>I</italic><sub>M3</sub> are the bridge arm currents as shown in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>; <italic>I</italic><sub>2</sub>, <italic>I</italic><sub>3</sub>, and <italic>I</italic><sub>4</sub> are the converter station currents as shown in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>; <italic>R</italic><sub>12</sub>, <italic>R</italic><sub>13</sub>, <italic>R</italic><sub>14</sub>, <italic>R</italic><sub>23</sub>, and <italic>R</italic><sub>34</sub> are the line resistances in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>.</p>
<p>In addition, the entire power of the MMC bridge arm of the MI-PFC must be balanced, and the AC power of the MI-PFC can only be exchanged inside it so that it can achieve its intended function of power balancing. To generate matched bridge arm currents, the MMC bridge arm voltage consists of both DC and AC components, as stated in the following <xref ref-type="disp-formula" rid="eqn-3">Eq. (3)</xref>. In the MI-PFC topology, C<sub>1</sub> is responsible for &#x201C;blocking DC current and passing AC current,&#x201D; whereas L<sub>dc</sub> is responsible for &#x201C;blocking AC current and passing DC current.&#x201D; The existence of the DC reactor L<sub>dc</sub> prevents the AC current of the MI-PFC from migrating to the DC grid line, causing it to flow only within the reactor.</p>
<p><disp-formula id="eqn-3"><label>(3)</label><mml:math id="mml-eqn-3" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn><mml:mrow><mml:mtext>AC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03C6;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mtext>AC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03C6;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn><mml:mrow><mml:mtext>AC</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03C6;</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>I</italic><sub>M<italic>i</italic>AC</sub> is the AC current fundamental amplitude of the bridge arm current, and <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:msub><mml:mi>&#x03C6;</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:math></inline-formula> is the initial phase angle of the AC current (<italic>i</italic> &#x003D; 1&#x007E;3).</p>
<p>In each bridge arm of the MI-PFC, there are <italic>N</italic> full-bridge sub-modules, and it is assumed that each sub-module receives a voltage of <italic>U</italic><sub>M<italic>i</italic>C</sub>. The following will serve to establish the relationship between the number of sub-modules, the voltage of each individual sub-module, and the voltage of each bridge arm.</p>
<p><disp-formula id="eqn-4"><label>(4)</label><mml:math id="mml-eqn-4" display="block"><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>N</mml:mi><mml:mo>&#x2217;</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:math></disp-formula>where <italic>U</italic><sub>M<italic>i</italic></sub> is the bridge arm voltage, <italic>N</italic> is the number of submodules, and <italic>U</italic><sub>M<italic>i</italic>C</sub> is the submodule capacitance voltage (<italic>i</italic> &#x003D; 1&#x007E;3).</p>
<p>The bridge arm containing <italic>N</italic> sub-modules is analogous to a whole capacitor, and the voltage on each bridge arm is determined by</p>
<p><disp-formula id="eqn-5"><label>(5)</label><mml:math id="mml-eqn-5" display="block"><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mrow><mml:mtext>eq</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula>where <italic>U</italic><sub>M<italic>i</italic></sub> is the bridge arm voltage, <italic>I</italic><sub>M<italic>i</italic></sub> is the bridge arm current and C<sub>eq</sub> is the equivalent capacitor, and its value can be expressed as</p>
<p><disp-formula id="eqn-6"><label>(6)</label><mml:math id="mml-eqn-6" display="block"><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mrow><mml:mtext>eq</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mi>C</mml:mi><mml:mi>N</mml:mi></mml:mfrac></mml:math></disp-formula>where <italic>C</italic> is the capacitor value of each submodule.</p>
<p>The KCL equation for the capacitor between the two bridge arms is given as</p>
<p><disp-formula id="eqn-7"><label>(7)</label><mml:math id="mml-eqn-7" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>U</italic><sub>C1</sub> and <italic>U</italic><sub>C2</sub> are the voltages of the capacitors <italic>C</italic><sub>1</sub> and <italic>C</italic><sub>2</sub> between the bridge arms.</p>
<p>The mathematical model presented in the paper depicts station 1&#x2019;s constant voltage converter as being equivalent to the DC voltage source <italic>U</italic><sub>dc</sub>. The current source shunt capacitor is equated to the remaining three converter stations. The equivalent KCL equations for these three converter stations are listed below:</p>
<p><disp-formula id="eqn-8"><label>(8)</label><mml:math id="mml-eqn-8" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>34</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>34</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>U</italic><sub>2</sub>, <italic>U</italic><sub>3</sub>, and <italic>U</italic><sub>4</sub> and refer to the voltage of converter station 2, 3, and 4, respectively. <italic>C</italic><sub>S2</sub>, <italic>C</italic><sub>S3</sub>, and <italic>C</italic><sub>S4</sub> are the equivalent capacitors of the converter stations.</p>
<p>Following is how the final dynamic mathematical model of MI-PFC is obtained:</p>
<p><disp-formula id="eqn-9"><label>(9)</label><mml:math id="mml-eqn-9" display="block"><mml:mi mathvariant="bold-italic">X</mml:mi><mml:mo mathvariant="bold">=</mml:mo><mml:mrow><mml:mi mathvariant="bold-italic">A</mml:mi><mml:mi mathvariant="bold-italic">X</mml:mi></mml:mrow><mml:mo mathvariant="bold">+</mml:mo><mml:mrow><mml:mi mathvariant="bold-italic">B</mml:mi><mml:mi mathvariant="bold-italic">U</mml:mi></mml:mrow></mml:math></disp-formula></p>
<p><disp-formula id="eqn-10"><label>(10)</label><mml:math id="mml-eqn-10" display="block"><mml:mrow><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>X</mml:mi></mml:mstyle><mml:mi>T</mml:mi></mml:msup><mml:mo>=</mml:mo><mml:mtable><mml:mtr><mml:mtd><mml:mrow><mml:mo stretchy='false'>[</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mtext>M</mml:mtext><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mtext>M</mml:mtext><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>12</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>13</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>&#x22EF;</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>23</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>34</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mtext>C</mml:mtext><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mtext>C</mml:mtext><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo>&#x22EF;</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mtext>M</mml:mtext><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mtext>M</mml:mtext><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mtext>M</mml:mtext><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mo stretchy='false'>]</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:math></disp-formula></p>
<p><disp-formula id="eqn-11"><label>(11)</label><mml:math id="mml-eqn-11" display="block"><mml:msup><mml:mrow><mml:mi mathvariant="bold">U</mml:mi></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mtable columnalign="left left left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>dc</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mn>4</mml:mn></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula>where A and B are matrices of coefficients.</p>
<p>MPC requires a discrete MI-PFC dynamic mathematical model; hence, the aforementioned dynamic mathematical model is discretized using the forward Euler method, as indicated in <xref ref-type="disp-formula" rid="eqn-12">Eq. (12)</xref>.</p>
<p><disp-formula id="eqn-12"><label>(12)</label><mml:math id="mml-eqn-12" display="block"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>&#x2248;</mml:mo><mml:mfrac><mml:mrow><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>x</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula>where x(<italic>k</italic> &#x002B; 1) is the value of the variable at the next sampling time, x(<italic>k</italic>) is the value of the variable at the current sampling time, and <italic>T</italic><sub>s</sub> is the sampling cycle time.</p>
</sec>
<sec id="s4">
<label>4</label>
<title>MPC Strategy for MI-PFC</title>
<sec id="s4_1">
<label>4.1</label>
<title>The Overall Control Structure</title>
<p>The completion of MI-PFC&#x2019;s primary core function is dependent on the MMC full-bridge sub-module; thus, the initial control strategy also employs the direct current control method based on the PI control strategy, which is a standard linear control strategy. The direct current control method is separated into an outer-loop controller and an inner-loop current controller, requiring the setting of proportional and integral coefficients for the two control loops, resulting in a more complicated control structure. Moreover, when the control target is quickly altered, the linear control features of PI control would result in poor dynamic performance. As a nonlinear optimal control strategy, MPC has many benefits, including simple control structure, multi-objective control, and good dynamic performance. Consequently, with the rapid development of digital processors, FCS-MPC has been widely studied in the fields of power electronics and power drives. The principle of MPC is to obtain the predicted value of the controlled quantity at the next sampling time using a model predictive method under the action of a limited combination of switch states, thereby obtaining the optimal switch group state combination.</p>
<p>Therefore, the paper chooses model predictive control instead of conventional PI-based control for current and voltage control in MI-PFC. The overall structure of the proposed model predictive control strategy is described in detail below.</p>
<p><xref ref-type="fig" rid="fig-3">Fig. 3</xref> illustrates that the MI-PFC MPC strategy encompasses five essential components, comprising a discrete prediction model, reference value computation, model prediction, rolling optimization, sub-module voltage ordering, and nearest level approximation modulation. Included in the calculation of the reference value are the bridge arm voltage command and line current command reference values. The discrete MI-PFC mathematical model computes the bridge arm voltage and line current prediction values at the next sampling instant, inputs them into the defined cost function, and then performs rolling optimization to determine the ideal number of input submodules. In order to limit the rate of voltage change so that the bridge arm voltage does not fluctuate drastically and abruptly, the MPC strategy proposed in the paper limits the number of submodule inputs to the range &#x2212;1, 0, and 1, i.e., only up to one submodule can be increased, decreased, or held constant.</p>
<fig id="fig-3">
<label>Figure 3</label>
<caption>
<title>Block diagram of MPC strategy for MI-PFC</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-3.tif"/>
</fig>
<p>After establishing the ideal number of input submodules, <italic>N</italic>, submodule voltage sorting is used to balance the submodule capacitor voltage. Due to the low voltage level of MI-PFC, a small number of full-bridge submodules are utilized. In the sorting method, traditional bubble sorting is utilized, which compares the voltage of each submodule to the others and selects the submodule with the lowest voltage while charging and the submodule with the greatest voltage when discharging. Thus, the required sub-module can be chosen so that the voltage in the full-bridge type sub-module within the bridge arm can be kept as consistent as feasible, and the computing work of the MPC strategy can be minimized to some degree. The mathematical model is given as</p>
<p><disp-formula id="eqn-13"><label>(13)</label><mml:math id="mml-eqn-13" display="block"><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mo movablelimits="true" form="prefix">max</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mtable columnalign="left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>4</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:mtd></mml:mtr></mml:mtable><mml:mo>}</mml:mo></mml:mrow><mml:mrow><mml:mtext>&#xA0;discharging</mml:mtext></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mo movablelimits="true" form="prefix">min</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>1</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>2</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mtable columnalign="left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>3</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mn>4</mml:mn><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:mtd></mml:mtr></mml:mtable><mml:mo>}</mml:mo></mml:mrow><mml:mrow><mml:mtext>&#xA0;charging</mml:mtext></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>U</italic><sub>M<italic>i</italic>C</sub> is the sub-module capacitor voltage.</p>
<p>Using nearest level modulation (NLM), the final modulation signal is created for input into the IGBT of the submodule. <italic>U</italic><sub>s</sub>(<italic>t</italic>) is supposed to represent the instantaneous value of the modulating wave, and <italic>U</italic><sub>C</sub> represents the output voltage level span, where <italic>U</italic><sub>C</sub> is equal to the average value of the sub-module capacitor voltage for a well-balanced sub-module voltage. After each sampling interval, the reference wave <italic>U</italic><sub>s</sub>(<italic>t</italic>) is sampled, and <italic>N</italic> sub-modules are inserted into the bridge arm between each sampling period. The expression of the relationship is as follows:</p>
<p><disp-formula id="eqn-14"><label>(14)</label><mml:math id="mml-eqn-14" display="block"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mtext>round</mml:mtext></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mfrac><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mi>S</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mfrac><mml:mo>]</mml:mo></mml:mrow></mml:math></disp-formula>where the round (<italic>x</italic>) function returns the integer that is closest to <italic>x</italic>.</p>
<p>Since the MI-PFC used in the paper relies on the forward and reverse inputs of full-bridge submodules to complete the bidirectional control of the DC grid power flow, the calculated number of submodule inputs ranges from <italic>N</italic>, where a submodule input number greater than zero indicates the forward input of submodules and a submodule input number less than zero indicates the reverse input of submodules.</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Reference Value Computation</title>
<p>MI-PFC generates amplitude <italic>U</italic><sub>M<italic>i</italic>DC</sub> and <italic>U</italic><sub>M<italic>i</italic>AC</sub> DC voltage components and AC voltage components, the former of which is responsible for achieving the objective of managing grid currents and the latter of which is responsible for achieving the goal of balancing its own power.</p>
<p>In addition, <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref> depicts the DC line currents in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>.</p>
<p><disp-formula id="eqn-15"><label>(15)</label><mml:math id="mml-eqn-15" display="block"><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula>where <italic>P</italic><sub><italic>ij</italic></sub> is the line power, <italic>U</italic><sub><italic>i</italic></sub> is the converter station voltage, and <italic>I</italic><sub><italic>ij</italic></sub> is the line current (<italic>i</italic> &#x003D; 1&#x007E;3, <italic>j</italic> &#x003D; 1&#x007E;4).</p>
<p>It can be inferred that the value of <italic>U<sub>i</sub></italic> in <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref> remains consistent, enabling the tidal current controller to implement a mechanism that governs the line current with the intention of regulating the tidal current. Consequently, the line current and the bridge arm voltage are the two control objectives of the MI-PFC, and their respective reference values must be calculated. In <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, <italic>I</italic><sub><italic>ij</italic>ref</sub> is the DC line current reference value that is determined by the power flow control objective. Using the known line resistance and the constant power converter station current, i.e., <italic>U</italic><sub>M<italic>i</italic>ref</sub> in <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, it is possible to determine the bridge arm voltage reference value while adjusting the power flow current. Its specific computation can be taken from the literature on steady-state control strategies [<xref ref-type="bibr" rid="ref-28">28</xref>].</p>
<fig id="fig-4">
<label>Figure 4</label>
<caption>
<title>Sub-module input quantity calculation method</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-4.tif"/>
</fig>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Constructing the Cost Function</title>
<p>The two control objectives of the MI-PFC MPC strategy are bridge arm voltage and line current, and these two control quantities are included in the design of the cost function. Incorporating a weighting component into the multi-objective control problem further integrates it. Cost functions are created in accordance with the aforementioned control objectives.</p>
<p>(1) Bridge arm voltage</p>
<p>According to the discrete mathematical model of MI-PFC presented in the preceding section, the bridge arm voltage prediction equation can be derived as</p>
<p><disp-formula id="eqn-16"><label>(16)</label><mml:math id="mml-eqn-16" display="block"><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mrow><mml:mtext>eq</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2217;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:math></disp-formula>where <italic>U</italic><sub>M<italic>i</italic></sub>(<italic>k</italic> &#x002B; 1) is the predicted value of the bridge arm voltage at the next sampling time, <italic>U</italic><sub>M<italic>i</italic></sub>(<italic>k</italic>) is the actual value of the bridge arm voltage at this time, and <italic>T</italic><sub>s</sub> is the sampling time.</p>
<p>Therefore, the cost function for the control objective of the bridge arm voltage can be designed as follows:</p>
<p><disp-formula id="eqn-17"><label>(17)</label><mml:math id="mml-eqn-17" display="block"><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi><mml:mrow><mml:mtext>ref</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>U</italic><sub>M<italic>i</italic>ref</sub> is the bridge arm voltage reference value.</p>
<p>(2) Line current</p>
<p>The power distribution of a DC grid is determined by the direction and magnitude of the current on each line. Consequently, line current is the ultimate aim of DC current control. In the initial DC current control strategy, the bridge arm voltage command reference value is derived from the line current command reference value. The MPC strategy for MI-PFC directly uses the line current as a control target, which will manage the DC grid power flow more directly and more quickly, and is vital for enhancing the MI-PFC strategy&#x2019;s dynamic performance.</p>
<p>According to the discrete mathematical model of MI-PFC, the following line current prediction equation can be derived:</p>
<p><disp-formula id="eqn-18"><label>(18)</label><mml:math id="mml-eqn-18" display="block"><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac><mml:mo>&#x2217;</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:math></disp-formula>where <italic>I</italic><sub><italic>ij</italic></sub>(<italic>k</italic> &#x002B; 1) represents the anticipated value of the line current at the next sampling time, <italic>I</italic><sub><italic>ij</italic></sub>(<italic>k</italic>) represents the actual value of the line current at this time, and <italic>T</italic><sub>s</sub> represents the sampling time.</p>
<p>Consequently, the cost function for the line current control goal can be constructed as follows:</p>
<p><disp-formula id="eqn-19"><label>(19)</label><mml:math id="mml-eqn-19" display="block"><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext mathvariant="italic">ijref</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>I</italic><sub><italic>ij</italic>ref</sub> is the line current reference value.</p>
<p>Combining the aforementioned two control objectives, the overall cost function is created according to <xref ref-type="disp-formula" rid="eqn-20">Eq. (20)</xref>, where the weight factor must be found through trial generation and <italic>&#x03BB;</italic><sub>1</sub> &#x003D; 0.01, <italic>&#x03BB;</italic><sub>2</sub> &#x003D; 0.2 are assumed in the study.</p>
<p><disp-formula id="eqn-20"><label>(20)</label><mml:math id="mml-eqn-20" display="block"><mml:mi>J</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>M</mml:mtext></mml:mrow><mml:mi>i</mml:mi><mml:mrow><mml:mtext>ref</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>i</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi><mml:mrow><mml:mtext>&#xA0;ref&#xA0;</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>&#x03BB;</italic><sub>1</sub>, <italic>&#x03BB;</italic><sub>2</sub> are the weight factors.</p>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Simulation Analysis</title>
<sec id="s5_1">
<label>5.1</label>
<title>Model Validation</title>
<p>The dynamic mathematical model of the four-terminal DC grid comprising MI-PFC established in the previous publication is simulated in the form of a MATLAB program, and the electromagnetic transient model of MI-PFC employing a PI-based technique is constructed in MATLAB/Simulink. The specific parameters are shown in <xref ref-type="table" rid="table-1">Tables 1</xref>&#x2013;<xref ref-type="table" rid="table-4">4</xref>. The following simulations are done for the MI-PFC during normal operation. The length of lines L<sub>12</sub>/L<sub>14</sub> is 2 km, and the length of lines L<sub>13</sub>/L<sub>23</sub>/L<sub>34</sub> is 1 km.</p>
<table-wrap id="table-1">
<label>Table 1</label>
<caption>
<title>DC grid component parameters</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Items</th>
<th>Values</th>
</tr>
</thead>
<tbody>
<tr>
<td>DC line resistance</td>
<td>0.01 &#x03A9;/km</td>
</tr>
<tr>
<td>Bridge arm inductors L<sub>0</sub></td>
<td>53 mH</td>
</tr>
<tr>
<td>DC reactor L<sub>dc</sub></td>
<td>0.02 H</td>
</tr>
<tr>
<td>Sub-module capacitor C<sub>0</sub></td>
<td>7 mF</td>
</tr>
<tr>
<td>Bridge arm equivalent resistance R<sub>0</sub></td>
<td>1 &#x03A9;</td>
</tr>
<tr>
<td>Number of sub-modules <italic>N</italic><sub>MMC</sub></td>
<td>10</td>
</tr>
</tbody>
</table>
</table-wrap><table-wrap id="table-2">
<label>Table 2</label>
<caption>
<title>Converter station parameters</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Items</th>
<th>Values</th>
</tr>
</thead>
<tbody>
<tr>
<td>Converter station 1 equivalent capacitor C<sub>S1</sub></td>
<td>142.9 &#x03BC;F</td>
</tr>
<tr>
<td>Converter station 2 equivalent capacitor C<sub>S2</sub></td>
<td>142.9 &#x03BC;F</td>
</tr>
<tr>
<td>Converter station 3 equivalent capacitor C<sub>S3</sub></td>
<td>142.9 &#x03BC;F</td>
</tr>
</tbody>
</table>
</table-wrap><table-wrap id="table-3">
<label>Table 3</label>
<caption>
<title>Converter station control parameters</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Items</th>
<th>Values</th>
</tr>
</thead>
<tbody>
<tr>
<td>Voltage of converter station 1</td>
<td>200 kV</td>
</tr>
<tr>
<td>Capacity of converter station 2</td>
<td>100 MW</td>
</tr>
<tr>
<td>Capacity of converter station 3</td>
<td>400 MW</td>
</tr>
<tr>
<td>Capacity of converter station 4</td>
<td>200 MW</td>
</tr>
</tbody>
</table>
</table-wrap><table-wrap id="table-4">
<label>Table 4</label>
<caption>
<title>MI-PFC parameters</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Items</th>
<th>Values</th>
</tr>
</thead>
<tbody>
<tr>
<td>Bridge arm equivalent resistance R<sub>1</sub></td>
<td>0.01 &#x03A9;</td>
</tr>
<tr>
<td>Bridge arm inductors L<sub>1</sub></td>
<td>0.5 mH</td>
</tr>
<tr>
<td>Sub-module capacitor C<sub>sm</sub></td>
<td>1 mF</td>
</tr>
<tr>
<td>Number of sub-modules <italic>N</italic><sub>MMC</sub></td>
<td>4</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The model validation was performed by comparing the dynamic mathematical model of MI-PFC and the electromagnetic transient model under operating condition 1 as shown in <xref ref-type="table" rid="table-5">Table 5</xref> to verify the correctness and validity of the dynamic mathematical model used in the paper.</p>
<table-wrap id="table-5">
<label>Table 5</label>
<caption>
<title>Control parameters of condition 1</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th></th>
<th><italic>U</italic><sub>M<italic>i</italic>DC</sub>/kV</th>
<th><italic>U</italic><sub>M<italic>i</italic>AC</sub>/kV</th>
<th><italic>&#x03B8;</italic><sub><italic>i</italic></sub>/degree</th>
</tr>
</thead>
<tbody>
<tr>
<td><italic>i</italic> &#x003D; 1</td>
<td>&#x2212;1.29</td>
<td>4.129</td>
<td>0</td>
</tr>
<tr>
<td><italic>i</italic> &#x003D; 2</td>
<td>&#x2212;0.79</td>
<td>4.629</td>
<td>19</td>
</tr>
<tr>
<td><italic>i</italic> &#x003D; 3</td>
<td>1.21</td>
<td>&#x2212;4.209</td>
<td>71</td>
</tr>
</tbody>
</table>
</table-wrap>
<p><xref ref-type="fig" rid="fig-5">Figs. 5a</xref> and <xref ref-type="fig" rid="fig-6">6b</xref> depict the currents of each line of MI-PFC for the two models of four-terminal DC grids simulated by the MATLAB software, whereas <xref ref-type="fig" rid="fig-5">Figs. 5c</xref> and <xref ref-type="fig" rid="fig-5">5d</xref> depict the voltages of the converter station simulated by the two models. <italic>U</italic><sub>M1C</sub>, <italic>U</italic><sub>M2C</sub>, and <italic>U</italic><sub>M3C</sub> are the dynamic mathematical models simulated by the submodule capacitor voltages simulated by the dynamic mathematical model in <xref ref-type="fig" rid="fig-5">Figs. 5e</xref> and <xref ref-type="fig" rid="fig-5">5f</xref>. The comparison in <xref ref-type="fig" rid="fig-5">Fig. 5</xref> demonstrates that the DC current, converter station voltage, and sub-module capacitor voltage simulated by the dynamic mathematical model are highly consistent with the MI-PFC electromagnetic transient model employing the original PI control strategy, and the error is due to the fact that the mathematical model disregards the high harmonic components of the internal AC voltage of the MI-PFC. The above analysis demonstrates the correctness and precision of the article&#x2019;s dynamic mathematical model of MI-PFC.</p>
<fig id="fig-5">
<label>Figure 5</label>
<caption>
<title>Model verification simulation diagram</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-5.tif"/>
</fig><fig id="fig-6">
<label>Figure 6</label>
<caption>
<title>Simulation comparison of line current under condition 1</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-6.tif"/>
</fig>
</sec>
<sec id="s5_2">
<label>5.2</label>
<title>Verification of the Performance of the MI-PFC Model Control Strategy</title>
<p>To further demonstrate the benefits of the MPC strategy for MI-PFC, an electromagnetic transient model of the MI-PFC is constructed and simulated in MATLAB/Simulink. The parameters utilized in the preceding mathematical model are maintained.</p>
<p>(1) Condition 1: The current on control line L<sub>13</sub> falls as the current on control line L<sub>14</sub> rises.</p>
<p>Before MI-PFC is activated, the DC grid is stabilized with the following line currents: <italic>I</italic><sub>12</sub> &#x003D; 0.767 kA, <italic>I</italic><sub>13</sub> &#x003D; 1.8 kA, <italic>I</italic><sub>14</sub> &#x003D; 0.933 kA, <italic>I</italic><sub>34</sub> &#x003D; &#x2212;0.067 kA, and <italic>I</italic><sub>23</sub> &#x003D; &#x2212;0.267 kA. To maintain the consistency of the MPC strategy and the sampling frequency of nearest level approximation modulation, the sampling period for the main circuit parameters and the operation parameters of the Dinghai converter station of the Zhoushan multi-terminal flexible DC transmission project is set to 100 &#x03BC;s. The key control parameters in the simulation system are shown in <xref ref-type="table" rid="table-6">Table 6</xref>. <italic>I</italic><sub>13</sub> &#x003D; 1.5 kA and <italic>I</italic><sub>14</sub> &#x003D; 1.5 kA are the power flow current control goals for condition 1 at <italic>t</italic> &#x003D; 1 s. According to the operating principle described in the preceding study, the bridge arm voltage reference value is as shown in <xref ref-type="table" rid="table-5">Table 5</xref>.</p>
<table-wrap id="table-6">
<label>Table 6</label>
<caption>
<title>Key control parameters</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Parameters</th>
<th>Value</th>
</tr>
</thead>
<tbody>
<tr>
<td>Weighting factor &#x03BB;<sub>1</sub></td>
<td>0.01</td>
</tr>
<tr>
<td>Weighting factor &#x03BB;<sub>2</sub></td>
<td>0.2</td>
</tr>
<tr>
<td>Sampling time <italic>T</italic><sub>S</sub></td>
<td>100 &#x03BC;s</td>
</tr>
</tbody>
</table>
</table-wrap>
<p><xref ref-type="fig" rid="fig-6">Fig. 6a</xref> depicts the distribution of DC grid current. During this procedure, the power flow controller adjusts <italic>I</italic><sub>13</sub> and <italic>I</italic><sub>14</sub> to the desired values in approximately 0.3 s. <xref ref-type="fig" rid="fig-6">Fig. 6b</xref> depicts the line currents regulated by the MPC strategy technique described in the article, in which this power flow controller brings <italic>I</italic><sub>13</sub> and <italic>I</italic><sub>14</sub> to the desired control value in around 0.2 s. From <xref ref-type="fig" rid="fig-6">Fig. 6b</xref>, it can be observed that the DC current with the MPC strategy follows the reference value more quickly and steadily, and its fluctuations are lower than those with the original control strategy. Therefore, when the DC current controller is activated, the effect on the current of each line of the DC grid can be enhanced.</p>
<p>Compared to the bridge arm voltages <italic>U</italic><sub>M1</sub>, <italic>U</italic><sub>M2</sub>, and <italic>U</italic><sub>M3</sub> in <xref ref-type="fig" rid="fig-7">Fig. 7a</xref> with PI control, the bridge arm voltages <italic>U</italic><sub>M1</sub>, <italic>U</italic><sub>M2</sub>, and <italic>U</italic><sub>M3</sub> in <xref ref-type="fig" rid="fig-7">Fig. 7b</xref> with MPC strategy demonstrate that the MPC strategy designed in this paper provides superior control dynamic performance for the bridge arm voltages. As illustrated in <xref ref-type="fig" rid="fig-7">Fig. 7a</xref>, the bridge arm voltage under PI controlled condition 1 begins to vary significantly at 1.01 s, whereas in <xref ref-type="fig" rid="fig-7">Fig. 7b</xref>, the bridge arm voltage under MPC control strategy of condition 1 begins to change significantly only at 1.002 s.</p>
<fig id="fig-7">
<label>Figure 7</label>
<caption>
<title>Simulation comparison of bridge arm voltage under condition 1</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-7.tif"/>
</fig>
<p>In accordance with the paper&#x2019;s MPC strategy, the submodule capacitor voltage can follow the reference value more quickly based on the simulation comparison in <xref ref-type="fig" rid="fig-8">Fig. 8</xref>. In conjunction with the sub-module balancing algorithm, it can maintain the balance and stability of the sub-module capacitor voltage in the MI-PFC bridge arm. Similarly, the peak value of <italic>U</italic><sub>M3C</sub> in <xref ref-type="fig" rid="fig-8">Fig. 8b</xref> is lower than that in <xref ref-type="fig" rid="fig-8">Fig. 8a</xref> due to the more accurate following ability of submodule voltage brought about by the MPC strategy, indicating that the overcharging problem at the initial charging of submodule capacitor under PI control can be mitigated to some degree.</p>
<fig id="fig-8">
<label>Figure 8</label>
<caption>
<title>Simulation comparison of sub-module capacitor under condition 1</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-8.tif"/>
</fig>
<p>(2) Condition 2: The current on line L<sub>13</sub> decreases and the current on line L<sub>14</sub> increases.</p>
<p>After the DC grid is stabilized, the power flow control target for MI-PFC at <italic>t</italic> &#x003D; 1 s in condition 2 is set at <italic>I</italic><sub>13</sub> &#x003D; 1 kA and <italic>I</italic><sub>14</sub> &#x003D; 0.5 kA, and the sampling time is also 100 &#x03BC;s. The reference values of the bridge arm voltage in Condition 2 are shown in <xref ref-type="table" rid="table-7">Table 7</xref>.</p>
<table-wrap id="table-7">
<label>Table 7</label>
<caption>
<title>Control parameters of condition 2</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th></th>
<th><italic>U</italic><sub>M<italic>i</italic>DC</sub>/kV</th>
<th><italic>U</italic><sub>M<italic>i</italic>AC</sub>/kV</th>
<th><italic>&#x03B8;</italic><sub><italic>i</italic></sub>/degree</th>
</tr>
</thead>
<tbody>
<tr>
<td><italic>i</italic> &#x003D; 1</td>
<td>2</td>
<td>&#x2212;6.614</td>
<td>0</td>
</tr>
<tr>
<td><italic>i</italic> &#x003D; 2</td>
<td>&#x2212;2.5</td>
<td>&#x2212;6.114</td>
<td>124</td>
</tr>
<tr>
<td><italic>i</italic> &#x003D; 3</td>
<td>&#x2212;3</td>
<td>&#x2212;5.614</td>
<td>103</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In <xref ref-type="fig" rid="fig-9">Fig. 9b</xref>, <italic>I</italic><sub>13</sub> and <italic>I</italic><sub>14</sub> are the DC line currents of condition 2 with the MPC strategy, and it can be seen that <italic>I</italic><sub>13</sub> and <italic>I</italic><sub>14</sub> can be reduced to the target value in <xref ref-type="fig" rid="fig-9">Fig. 9b</xref>. In terms of the dynamic performance of the line current control in condition 2, the original PI control can be completed in about 0.3 s to follow the line current reference value, while the MPC strategy can be completed in about 0.2 s to basically follow the line current reference value. At the same time, the fluctuation of the current magnitude will be reduced.</p>
<fig id="fig-9">
<label>Figure 9</label>
<caption>
<title>Simulation comparison of line current under condition 2</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-9.tif"/>
</fig>
<p><xref ref-type="fig" rid="fig-10">Fig. 10b</xref> depicts the bridge arm voltages in the DC grid under the MPC strategy, while <xref ref-type="fig" rid="fig-10">Fig. 10a</xref> depicts the bridge arm voltages with the original PI control strategy. Comparing the two figures reveals, similar to condition 1, that the bridge arm voltage under the MPC strategy responds approximately 0.098 s faster than the bridge arm voltage under PI control.</p>
<fig id="fig-10">
<label>Figure 10</label>
<caption>
<title>Simulation comparison of bridge arm voltage under condition 2</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-10.tif"/>
</fig>
<p>In <xref ref-type="fig" rid="fig-11">Fig. 11b</xref>, the submodule capacitor voltages <italic>U</italic><sub>M1C</sub>, <italic>U</italic><sub>M2C</sub>, and <italic>U</italic><sub>M3C</sub> under the MPC strategy have a faster control target tracking speed than the submodule capacitor voltages <italic>U</italic><sub>M1C</sub>, <italic>U</italic><sub>M2C,</sub> and <italic>U</italic><sub>M3C</sub> with the original PI control strategy. Comparing the peak value of <italic>U</italic><sub>M1C</sub> in <xref ref-type="fig" rid="fig-11">Figs. 11b</xref> and <xref ref-type="fig" rid="fig-11">11a</xref>, there is a reduction utilizing the MPC strategy designed in the study, demonstrating that the control strategy can mitigate the problem of submodule capacitor overcharging to some extent.</p>
<fig id="fig-11">
<label>Figure 11</label>
<caption>
<title>Simulation comparison of sub-module capacitor under condition 2</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-11.tif"/>
</fig>
<p>The simulation analysis of the aforementioned two operating conditions demonstrates that the MPC strategy designed in this paper can control line current, bridge arm voltage, and sub-module voltage effectively under different operating conditions and has a faster dynamic response and more accurate control effect than the original control strategy for the designed control target.</p>
</sec>
<sec id="s5_3">
<label>5.3</label>
<title>Precharge Method for MI-PFC</title>
<p>Incorporating an initial energy for the storage element in the MI-PFC can effectively minimize the intensity of the transient process, which subsequently helps reduce the voltage shock during the start-up phase [<xref ref-type="bibr" rid="ref-29">29</xref>]. Therefore, this research endeavors to include a precharge method in the start-up process by following the subsequent steps:</p>
<p>(1) Firstly, the sub-module will be blocked so that the DC line current charges the sub-module of the MMC bridge arm. The sub-module capacitors will have an uneven charging rate due to the varying DC line current. Therefore, to avoid over-voltage, the sub-module will be unblocked once the capacitor of the largest sub-module in all bridge arms reaches 85% of the rated value. The sub-module capacitor will then store energy.</p>
<p>(2) Secondly, the sub-module will be removed to enable the DC line current to magnetize the inductor L<sub>1</sub>. During the initial stage, the charging current will decrease as the capacitor voltage rises, leading to an inability to magnetize the bridge arm inductor. However, after removing the sub-module, the DC line current starts to magnetize the inductor, with <italic>I</italic><sub>M1</sub> &#x2212; <italic>I</italic><sub>M<italic>i</italic></sub> indicating the magnetization is complete, and the MI-PFC is ready for power flow control.</p>
<p>As shown in <xref ref-type="fig" rid="fig-12">Figs. 12</xref> and <xref ref-type="fig" rid="fig-13">13</xref>, simulation experiments were conducted to demonstrate the effect of the precharge method on the control of sub-module capacitor voltage fluctuations. The simulation results show that the voltage fluctuation of small module capacitors can be reduced by adding the precharge method, and the parameters of the capacitor and IGBT can be more easily selected.</p>
<fig id="fig-12">
<label>Figure 12</label>
<caption>
<title>Precharge sub-module capacitor voltage of condition 1</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-12.tif"/>
</fig><fig id="fig-13">
<label>Figure 13</label>
<caption>
<title>Precharge sub-module capacitor voltage of condition 2</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_28965-fig-13.tif"/>
</fig>
</sec>
</sec>
<sec id="s6">
<label>6</label>
<title>Conclusion</title>
<p>The article proposes a MPC strategy for the MI-PFC that has the benefits of good dynamic performance, a simple control structure, and the ability to implement multi-objective control. The subsequent conclusions are reached:</p>
<p>(1) The MPC strategy presented in this study can replace the original PI-based voltage control and current control structures. It has two control objectives: DC line current and bridge arm voltage.</p>
<p>(2) The dynamic mathematical model of a multi-terminal DC grid with a MI-PFC can accurately represent the relationship between each controlled quantity and the control quantity, thereby achieving the prediction output of voltage and current commands and completing the model prediction link of the paper&#x2019;s control strategy.</p>
<p>(3) The MPC strategy created in this study has a certain degree of dynamic performance improvement and optimization compared to the original control strategy; it also has greater control accuracy and simplifies the construction of the control strategy for the MI-PFC. The model predictive control strategy designed in the article can improve the DC line current control speed by 30%, and the control time to complete the line current in the four-terminal DC network used in the MPC strategy in the article is 0.2 s, which saves 0.1 s compared to PI control. Under the PI control strategy, the bridge arm voltage starts to change significantly at t &#x003D; 1.01 s, while with the MPC control strategy, the bridge arm voltage starts to change significantly at t &#x003D; 1.002 s the bridge arm voltage strategy starts to change significantly, which is about a 20% improvement in dynamic response speed.</p>
<p>(4) The paper utilizes a model predictive control approach to enhance the control performance of MI-PFC. A cost function is employed to normalize the two control objectives; however, certain limitations still exist. The cost function used in the study relies on the traditional trial generation method, which may be challenging and time-consuming in determining the process. Additionally, the capacitor voltage sorting technique utilized in the research is based on the traditional bubble sorting algorithm, which requires a substantial amount of computation. Voltage sorting algorithm used is the traditional bubble sorting, which requires a large amount of calculation.</p>
<p>(5) Moreover, while MI-PFC is a unique inter-line DCPFC topology, the model predictive control strategy developed may not be suitable for other inter-line DC current controllers. Therefore, proposed extensions to more interline DC current controllers utilizing general models will be explored in future studies to enhance the scope of this research.</p>
</sec>
</body>
<back>
<ack><p>None.</p></ack>
<sec><title>Funding Statement</title>
<p>The research is funded by National Natural Science Foundation of China (52177074).</p>
</sec>
<sec><title>Author Contributions</title>
<p>The authors confirm contribution to the paper as follows: study conception and design: He Wang, Jing Bian; data collection: Guanye Shen; analysis and interpretation of results: He Wang, Xiangsheng Xu, Jing Bian; draft manuscript preparation: Xiangsheng Xu, Jing Bian. All authors reviewed the results and approved the final version of the manuscript.</p>
</sec>
<sec sec-type="data-availability"><title>Availability of Data and Materials</title>
<p>The authors confirm that the data supporting the findings of this study are available within the article. And the additional data that support the findings of this study are available on request from the corresponding author, upon reasonable request.</p>
</sec>
<sec sec-type="COI-statement"><title>Conflicts of Interest</title>
<p>The authors declare that they have no conflicts of interest to report regarding the present study.</p>
</sec>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1"><label>1.</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><surname>Tang</surname>, <given-names>G. F.</given-names></string-name></person-group> (<year>2010</year>). <source>High voltage DC transmission technology based on voltage source converter</source>. <publisher-loc>China</publisher-loc>: <publisher-name>China Electric Power Press</publisher-name>.</mixed-citation></ref>
<ref id="ref-2"><label>2.</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><surname>Zhao</surname>, <given-names>C. Y.</given-names></string-name></person-group> (<year>2014</year>). <source>Modeling and simulation techniques for flexible DC power transmission</source>. <publisher-loc>China</publisher-loc>: <publisher-name>China Electric Power Press</publisher-name>.</mixed-citation></ref>
<ref id="ref-3"><label>3.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Gomis-Bellmunt</surname>, <given-names>O.</given-names></string-name>, <string-name><surname>Sau-Bassols</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Prieto-Araujo</surname>, <given-names>E.</given-names></string-name>, <string-name><surname>Cheah-Mane</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Flexible converters for meshed HVDC grids: From flexible AC transmission systems (FACTS) to flexible DC grids</article-title>. <source>IEEE Transactions on Power Delivery</source><italic>,</italic> <volume>35</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>2</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2016.11.233</pub-id></mixed-citation></ref>
<ref id="ref-4"><label>4.</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><surname>Barker</surname>, <given-names>C. D.</given-names></string-name>, <string-name><surname>Whitehouse</surname>, <given-names>R. S.</given-names></string-name></person-group> (<year>2012</year>). <article-title>A current flow controller for use in HVDC grids</article-title>. <conf-name>10th IET International Conference on AC and DC Power Transmission (ACDC 2012)</conf-name>, pp. <fpage>1</fpage>&#x2013;<lpage>5</lpage>. <publisher-loc>Birmingham, UK</publisher-loc>, <publisher-name>IET</publisher-name>. <pub-id pub-id-type="doi">10.1049/cp.2012.1973</pub-id></mixed-citation></ref>
<ref id="ref-5"><label>5.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Sau-Bassols</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Prieto-Araujo</surname>, <given-names>E.</given-names></string-name>, <string-name><surname>Gomis-Bellmunt</surname>, <given-names>O.</given-names></string-name></person-group> (<year>2015</year>). <article-title>Modelling and control of an interline current flow controller for meshed HVDC grids</article-title>. <source>IEEE Transactions on Power Delivery</source><italic>,</italic> <volume>32</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>11</fpage>&#x2013;<lpage>22</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2015.2513160</pub-id></mixed-citation></ref>
<ref id="ref-6"><label>6.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Diab</surname>, <given-names>H. Y.</given-names></string-name>, <string-name><surname>Marei</surname>, <given-names>M. I.</given-names></string-name>, <string-name><surname>Tennakoon</surname>, <given-names>S. B.</given-names></string-name></person-group> (<year>2016</year>). <article-title>Reduced switch count topology of current flow control apparatus for MTDC grids</article-title>. <source>Journal of Power Electronics</source><italic>,</italic> <volume>16</volume><italic>(</italic><issue>5</issue><italic>),</italic> <fpage>1743</fpage>&#x2013;<lpage>1751</lpage>. <pub-id pub-id-type="doi">10.6113/JPE.2016.16.5.1743</pub-id></mixed-citation></ref>
<ref id="ref-7"><label>7.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Diab</surname>, <given-names>H. Y.</given-names></string-name>, <string-name><surname>Marei</surname>, <given-names>M. I.</given-names></string-name>, <string-name><surname>Tennakoon</surname>, <given-names>S. B.</given-names></string-name></person-group> (<year>2016</year>). <article-title>Operation and control of an insulated gate bipolar transistor-based current controlling device for power flow applications in multi-terminal high-voltage direct current grids</article-title>. <source>IET Power Electronics</source><italic>,</italic> <volume>10</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>305</fpage>&#x2013;<lpage>315</lpage>. <pub-id pub-id-type="doi">10.1049/iet-pel.2015.0525</pub-id></mixed-citation></ref>
<ref id="ref-8"><label>8.</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><surname>Diab</surname>, <given-names>H.</given-names></string-name>, <string-name><surname>Tennakoon</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Gould</surname>, <given-names>C.</given-names></string-name>, <string-name><surname>Marei</surname>, <given-names>M. I.</given-names></string-name></person-group> (<year>2015</year>). <article-title>An investigation of power flow control methods in multi terminal high voltage DC grids</article-title>. <conf-name>2015 50th International Universities Power Engineering Conference (UPEC)</conf-name>, pp. <fpage>1</fpage>&#x2013;<lpage>5</lpage>. <publisher-loc>Stoke on Trent, UK</publisher-loc>, <publisher-name>IEEE</publisher-name>. <pub-id pub-id-type="doi">10.1109/UPEC.2015.7339839</pub-id></mixed-citation></ref>
<ref id="ref-9"><label>9.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname>, <given-names>H.</given-names></string-name>, <string-name><surname>Bian</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>G. Q.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>Z. H.</given-names></string-name>, <string-name><surname>Yang</surname>, <given-names>Y.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2017</year>). <article-title>Multi-port DC power flow controller for flexible DC power grid</article-title>. <source>Automation of Electric Power Systems</source><italic>,</italic> <volume>41</volume><italic>(</italic><issue>22</issue><italic>),</italic> <fpage>7</fpage>. <pub-id pub-id-type="doi">10.7500/AEPS20170105004</pub-id></mixed-citation></ref>
<ref id="ref-10"><label>10.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname>, <given-names>G. Q.</given-names></string-name>, <string-name><surname>Bian</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>H.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>Z. H.</given-names></string-name></person-group> (<year>2020</year>). <article-title>A circulating current interline DC power flow controller for DC grid</article-title>. <source>Transactions of China Electrotechnical Society</source><italic>,</italic> <volume>35</volume><italic>(</italic><issue>5</issue><italic>),</italic> <fpage>10</fpage>. <pub-id pub-id-type="doi">10.19595/j.cnki.1000-6753.tces.190299</pub-id></mixed-citation></ref>
<ref id="ref-11"><label>11.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Fang</surname>, <given-names>W. X.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>J. W.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>J. C.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>J. Q.</given-names></string-name>, <string-name><surname>Shi</surname>, <given-names>G.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2021</year>). <article-title>A multiport embedded DC power flow controller for modular multilevel converter</article-title>. <source>Proceedings of the CSEE</source><italic>,</italic> <volume>41</volume><italic>(</italic><issue>7</issue><italic>),</italic> <fpage>2484</fpage>&#x2013;<lpage>2495</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2021.3126983</pub-id></mixed-citation></ref>
<ref id="ref-12"><label>12.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liu</surname>, <given-names>Z. C.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>J. W.</given-names></string-name>, <string-name><surname>Shi</surname>, <given-names>G.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>J. Q.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>Y. X.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2022</year>). <article-title>An embedded DC power flow controller and its control strategy for flexible high voltage DC power system</article-title>. <source>Proceedings of the CSEE</source><italic>,</italic> <volume>42</volume><italic>(</italic><issue>19</issue><italic>),</italic> <fpage>7170</fpage>&#x2013;<lpage>7182</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.211458</pub-id></mixed-citation></ref>
<ref id="ref-13"><label>13.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wu</surname>, <given-names>W.</given-names></string-name>, <string-name><surname>Wu</surname>, <given-names>X. Z.</given-names></string-name>, <string-name><surname>Jing</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Wei</surname>, <given-names>M. H.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>Y.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Studies of a novel DC power flow controller and its control strategy</article-title>. <source>Proceedings of the CSEE</source><italic>,</italic> <volume>39</volume><italic>(</italic><issue>13</issue><italic>),</italic> <fpage>3744</fpage>&#x2013;<lpage>3757</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.181365</pub-id></mixed-citation></ref>
<ref id="ref-14"><label>14.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhong</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Zhu</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Chi</surname>, <given-names>Y. Liu</given-names></string-name>, <string-name><surname>S.</surname> <given-names>Q.</given-names></string-name>, <string-name><surname>Cai</surname>, <given-names>X.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Composite DC power flow controller</article-title>. <source>IEEE Transactions on Power Electronics</source><italic>,</italic> <volume>35</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>3530</fpage>&#x2013;<lpage>3542</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2019.2936773</pub-id></mixed-citation></ref>
<ref id="ref-15"><label>15.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Xu</surname>, <given-names>Q.</given-names></string-name>, <string-name><surname>Huang</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Chu</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Shuai</surname>, <given-names>Z.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2021</year>). <article-title>Analysis and control of modular multi-terminal DC power flow controller with fault current limiting function</article-title>. <source>Journal of Modern Power System and Clean Energy</source><italic>,</italic> <volume>9</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>1375</fpage>&#x2013;<lpage>1385</lpage>. <pub-id pub-id-type="doi">10.35833/MPCE.2020.000876</pub-id></mixed-citation></ref>
<ref id="ref-16"><label>16.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname>, <given-names>P.</given-names></string-name>, <string-name><surname>Feng</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Liu</surname>, <given-names>P.</given-names></string-name>, <string-name><surname>Jiang</surname>, <given-names>N.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>X. P.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Nyquist stability analysis and capacitance selection for DC current flow controllers in meshed multi-terminal HVDC grids</article-title>. <source>CSEE Journal of Power and Energy Systems</source><italic>,</italic> <volume>7</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>114</fpage>&#x2013;<lpage>127</lpage>. <pub-id pub-id-type="doi">10.17775/CSEEJPES.2020.00320</pub-id></mixed-citation></ref>
<ref id="ref-17"><label>17.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhong</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Zhu</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>H.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2020</year>). <article-title>Modular interline DC power flow controller</article-title>. <source>IEEE Transactions on Power Electronics</source><italic>,</italic> <volume>35</volume><italic>(</italic><issue>11</issue><italic>),</italic> <fpage>11707</fpage>&#x2013;<lpage>11719</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2020.2989197</pub-id></mixed-citation></ref>
<ref id="ref-18"><label>18.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wu</surname>, <given-names>W.</given-names></string-name>, <string-name><surname>Wu</surname>, <given-names>X. Z.</given-names></string-name>, <string-name><surname>Zhao</surname>, <given-names>Y. M.</given-names></string-name>, <string-name><surname>Jing</surname>, <given-names>L.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Active damping control of multi-port DC power follow controller for suppressing sub-synchronous oscillation</article-title>. <source>Power System Technology</source><italic>,</italic> <volume>45</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>1388</fpage>&#x2013;<lpage>1400</lpage>. <pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2020.0328</pub-id></mixed-citation></ref>
<ref id="ref-19"><label>19.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname>, <given-names>P. Y.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>J. Y.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>X. P.</given-names></string-name>, <string-name><surname>Guo</surname>, <given-names>C. Y.</given-names></string-name></person-group> (<year>2022</year>). <article-title>Characteristics and control &#x0026; protection strategy of a multi-line CFC integrated M2TDC grid under non-instantaneous open-circuit faults</article-title>. <source>Power System Technology</source><italic>,</italic> <volume>46</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>11</fpage>. <pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2021.0154</pub-id></mixed-citation></ref>
<ref id="ref-20"><label>20.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liao</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>N.</given-names></string-name>, <string-name><surname>Qin</surname>, <given-names>Z.</given-names></string-name>, <string-name><surname>Purgat</surname>, <given-names>P.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>Q.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2021</year>). <article-title>Coordination control of power flow controller and hybrid DC circuit breaker in MVDC distribution networks</article-title>. <source>Journal of Modern Power System and Clean Energy</source><italic>,</italic> <volume>9</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>1257</fpage>&#x2013;<lpage>1268</lpage>. <pub-id pub-id-type="doi">10.35833/MPCE.2021.000299</pub-id></mixed-citation></ref>
<ref id="ref-21"><label>21.</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><surname>Camacho</surname>, <given-names>E. F.</given-names></string-name>, <string-name><surname>Bordons</surname>, <given-names>C.</given-names></string-name></person-group> (<year>2004</year>). <source>Model predictive control</source>. <publisher-loc>London, UK</publisher-loc>: <publisher-name>Springer</publisher-name>.</mixed-citation></ref>
<ref id="ref-22"><label>22.</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><surname>Xi</surname>, <given-names>Y. G.</given-names></string-name></person-group> (<year>2013</year>). <source>Predictive control</source>. <publisher-loc>China</publisher-loc>: <publisher-name>Nation Defense Industry Press</publisher-name>.</mixed-citation></ref>
<ref id="ref-23"><label>23.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Qin</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Saeedifard</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2012</year>). <article-title>Predictive control of a modular multilevel converter for a back-to-back HVDC system</article-title>. <source>IEEE Transactions on Power Delivery</source><italic>,</italic> <volume>27</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>1538</fpage>&#x2013;<lpage>1547</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRD.2012.2191577</pub-id></mixed-citation></ref>
<ref id="ref-24"><label>24.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Bocker</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Freudenberg</surname>, <given-names>B.</given-names></string-name>, <string-name><surname>The</surname>, <given-names>A.</given-names></string-name>, <string-name><surname>Dieckerhoff</surname>, <given-names>S.</given-names></string-name></person-group> (<year>2014</year>). <article-title>Experimental comparison of model predictive control and cascaded control of the modular multilevel converter</article-title>. <source>IEEE Transactions on Power Electronics</source><italic>,</italic> <volume>30</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>422</fpage>&#x2013;<lpage>430</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2014.2309438</pub-id></mixed-citation></ref>
<ref id="ref-25"><label>25.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zheng</surname>, <given-names>G.</given-names></string-name>, <string-name><surname>Peng</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Yuan</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Wu</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Guo</surname>, <given-names>G.</given-names></string-name></person-group> (<year>2016</year>). <article-title>Design and experimental evaluation of fast model predictive control for modular multilevel converters</article-title>. <source>IEEE Transactions on Industrial Electronics</source><italic>,</italic> <volume>63</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>3845</fpage>&#x2013;<lpage>3856</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2015.2497254</pub-id></mixed-citation></ref>
<ref id="ref-26"><label>26.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liang</surname>, <given-names>Y. Y.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>T.</given-names></string-name>, <string-name><surname>Liu</surname>, <given-names>J. Z.</given-names></string-name>, <string-name><surname>Yang</surname>, <given-names>Q. X.</given-names></string-name>, <string-name><surname>Mei</surname>, <given-names>H. M.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2016</year>). <article-title>The applications of the model predictive control for MMC-HVDC</article-title>. <source>Transactions of China Electrotechnical Society</source><italic>,</italic> <volume>31</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>128</fpage>&#x2013;<lpage>138</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-6753.2016.01.016</pub-id></mixed-citation></ref>
<ref id="ref-27"><label>27.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhu</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Fu</surname>, <given-names>X. W.</given-names></string-name>, <string-name><surname>Hu</surname>, <given-names>X. B.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>X. R.</given-names></string-name></person-group> (<year>2014</year>). <article-title>Model predictive control of modular multilevel converter for HVDC system</article-title>. <source>Power System Protection and Control</source><italic>,</italic> <volume>42</volume><italic>(</italic><issue>16</issue><italic>),</italic> <fpage>1</fpage>&#x2013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.7667/j.issn.1674-3415.2014.16.001</pub-id></mixed-citation></ref>
<ref id="ref-28"><label>28.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Bian</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>G. Q.</given-names></string-name>, <string-name><surname>Jin</surname>, <given-names>R. K.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>H.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Multiport interline DC power flow controller and control strategy optimization</article-title>. <source>Proceedings of the CSEE</source><italic>,</italic> <volume>40</volume><italic>(</italic><issue>12</issue><italic>),</italic> <fpage>3980</fpage>&#x2013;<lpage>3990</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.190564</pub-id></mixed-citation></ref>
<ref id="ref-29"><label>29.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tian</surname>, <given-names>K.</given-names></string-name>, <string-name><surname>Wu</surname>, <given-names>B.</given-names></string-name>, <string-name><surname>Du</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Xu</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Cheng</surname>, <given-names>Z.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2015</year>). <article-title>A simple and cost-effective precharge method for modular multilevel converters by using a low-voltage DC source</article-title>. <source>IEEE Transactions on Power Electronics</source><italic>,</italic> <volume>31</volume><italic>(</italic><issue>7</issue><italic>),</italic> <fpage>5321</fpage>&#x2013;<lpage>5329</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2015.2484222</pub-id></mixed-citation></ref>
</ref-list>
</back></article>