<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Archiving and Interchange DTD v2.3 20070202//EN" "archivearticle.dtd">
<article article-type="methods-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">901708</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2022.901708</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Methods</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Black Start and Voltage Establishment Strategy for PMSG-Based Wind Turbine</article-title>
<alt-title alt-title-type="left-running-head">Ji et al.</alt-title>
<alt-title alt-title-type="right-running-head">Black Start Using Fixed Load</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Ji</surname>
<given-names>Xiaotong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Dan</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Hu</surname>
<given-names>Pan</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1747899/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Cao</surname>
<given-names>Kan</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jiang</surname>
<given-names>Kezheng</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Huang</surname>
<given-names>Guohang</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1690717/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Huang</surname>
<given-names>Sheng</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1374826/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>State Grid Hubei Electric Power CO., LTD.</institution>, <addr-line>Wuhan</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Hubei Electric Power Research Institute</institution>, <addr-line>Wuhan</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>College of Electrical and Information Engineering</institution>, <institution>Hunan University</institution>, <addr-line>Changsha</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1259467/overview">Liansong Xiong</ext-link>, Nanjing Institute of Technology (NJIT), China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1737753/overview">Pengda Wang</ext-link>, Technical University of Denmark, Denmark</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1708972/overview">Weiyu Bao</ext-link>, Shandong University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Guohang Huang, <email>mrhung33@outlook.com</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Process and Energy Systems Engineering, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>11</day>
<month>05</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>901708</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>03</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>04</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Ji, Liu, Hu, Cao, Jiang, Huang and Huang.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Ji, Liu, Hu, Cao, Jiang, Huang and Huang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>During the black start, backup ac power sources such as diesel generators can offer line-side voltage reference for wind turbine and keep line-side converter of wind turbine work properly so that the dc-link capacitor voltage within converter can be established without overcharge. This study proposes a black start control strategy and line-side voltage establishment method for PMSG-based wind turbines with no ac power source. Unlike the traditional control strategy of full power grid-connected converters, the dc-link voltage within back-to-back full power converters of wind turbines can be controlled by generator-side converters, and the line-side voltage can be established by line-side converters with the help of fixed loads. The mechanical power can be balanced by pitch angle control, and the power unbalance between mechanical power and electrical power will be reflected in the rotor speed of PMSG. By this method, a single wind turbine can establish the line-side voltage with no extra backup ac power source, offering voltage reference for the other wind turbines during the black start of a wind power plant.</p>
</abstract>
<kwd-group>
<kwd>black start</kwd>
<kwd>PMSG based</kwd>
<kwd>wind turbine</kwd>
<kwd>voltage establishment</kwd>
<kwd>power balance</kwd>
<kwd>voltage reference</kwd>
</kwd-group>
<contract-num rid="cn001">4000-202122070A-0-0-00</contract-num>
<contract-sponsor id="cn001">Science and Technology Project of State Grid<named-content content-type="fundref-id">10.13039/501100013096</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Black start is an inevitable starting procedure during the progress of grid restoration after a blackout happens (<xref ref-type="bibr" rid="B4">Chou et al., 2013</xref>; <xref ref-type="bibr" rid="B6">Erdiwansyah et al., 2021</xref>). Power sources with high reliability, such as thermal generators and hydro generators (<xref ref-type="bibr" rid="B3">Benato et al., 2019</xref>; <xref ref-type="bibr" rid="B9">Lindstrom, 1990</xref>), are usually used as backup power sources during black start. With the huge progress in virtual inertial control of wind power systems, wind turbines now can work as black start power sources (<xref ref-type="bibr" rid="B11">Pape and Kazerani, 2020</xref>). However, a complete black start scheme for wind turbines still requires reliable power supplies and stable voltage references, which are normally offered by backup diesel generators (<xref ref-type="bibr" rid="B17">Tang et al., 2017</xref>).</p>
<p>Generally, the first step of the black start procedure in a wind turbine is the startup of electrical subsystems. In order to boot up the electrical subsystems such as secondary circuits, pitch and yaw actuators, and cooling systems, the deployment of batteries or energy storage systems (<xref ref-type="bibr" rid="B10">Liu et al., 2021</xref>) is one of the best and most necessary solutions (<xref ref-type="bibr" rid="B15">Satpathy et al., 2014</xref>). After the startup of subsystems, the dc-link voltage must be established to make sure the power can be transferred from the generator to the line-side. The establishment of dc-link voltage is usually accomplished by a line-side converter or dc-link energy storage system. In the former case, the line-side voltage must be built by a backup motor before the grid line-converter can operate as expected, and a DFIG grid-forming technique using torque synchronization and voltage droop in the generator-side converter to form transient ac voltage has been proposed (<xref ref-type="bibr" rid="B12">Rodriguez-Amenedo et al., 2021</xref>), in which the dc-link voltage is maintained by the line-side converter; in the latter case, the capacity of energy storage system needs to be much larger since it will not only provide power to the subsystem but also help build the dc-link voltage (<xref ref-type="bibr" rid="B18">Xu et al., 2012</xref>; <xref ref-type="bibr" rid="B5">Deng et al., 2017</xref>; <xref ref-type="bibr" rid="B16">Sun et al., 2018</xref>).</p>
<p>When it comes to the circumstances with no extra back generators, a unique black start strategy for wind power plants using a centralized high-voltage dc-link converter has been proposed (<xref ref-type="bibr" rid="B14">Sakamuri et al., 2019</xref>; <xref ref-type="bibr" rid="B13">Sabor&#xed;o-Romano et al., 2019</xref>). By transmitting power from wind turbines to centralized high-voltage dc-link, the high centralized dc-link voltage can be established by the onshore high-voltage dc converter. However, the dc-link voltage of each single wind turbine which is supposed to be controlled by the generator-side converter in this study is considered constant. An active power-sharing method of wind turbine using the generator-side converter to control dc-link voltage and the grid-side converter to control active power is proposed (<xref ref-type="bibr" rid="B7">Fathabadi, 2017</xref>). The power balance is achieved between two wind turbines with the power-sharing strategy and one common local load, but the power flow between the grid and generator of one single turbine needs to be balanced during the effect of the power-sharing strategy with two turbines. A novel centralized control strategy based on the look-up table to ensure optimal power sharing is proposed (<xref ref-type="bibr" rid="B1">Alavi and Ghazi, 2022</xref>), using the optimal results from the centralized controller and data in the look-up table to achieve the optimal active power-sharing between multiple turbines. Though the optimal control command is obtained from the available wind power, it can still be a good reference when it comes to power balance with loads. In order to deal with the challenge caused by irregular, nonlinear, and nonstationary characteristics of wind power, an uncertainty modeling method based on the prediction of wind power is proposed (<xref ref-type="bibr" rid="B19">Yan et al., 2021</xref>), which, on the other hand, provides the possibility of grid-forming by wind power under variable wind speed and multiple loads.</p>
<p>From the recent research, two main challenges of black start control of wind turbines can be summarized: one is the establishment of dc and ac link voltage, and the other is the balance of power between generator and load. When a blackout happens with no backup generator, there is no grid line voltage, and dc-link voltage cannot be built by the line-side converter at the beginning. So, using the generator-side converter to build dc-link voltage is a better option. After the dc-link voltage rises and stays at a certain value, power begins to flow into the line-side converter, and the line-side voltage starts to rise on local load, which leads to the problem of grid voltage control. Since the volt-ampere characteristic on fixed load under certain power input is also fixed, the line-side voltage forming problem can be turned into a power control problem for wind turbines, which is the main idea of this study.</p>
<p>The main contribution of this study is as follows: instead of backup diesel generators, fixed three-phase local loads are used to establish line-side voltage, which is much better for costs; the proposed power balance method is proposed to maintain the power balance between mechanical power input and electrical power consumption; the first wind turbine using the proposed black start control strategy in this study can be used as line-side voltage reference during the black start of a wind power plant, and the black start of the whole wind power plant can be completed with no extra backup generator. However, the power supply of subsystems in the first turbine still needs to be satisfied by power storage or battery.</p>
<p>In the proposed black start method for wind turbines, the dc-link voltage is established by the generator-side converter, and the line-side voltage is built by the line-side converter with power balancing control of wind turbine and fixed load in ac line. The crowbar circuit is used to maintain the stability of rotor speed and dc-link voltage. By this method, the black start of a single wind turbine can be accomplished, and the line-side voltage can be established, offering voltage reference for the other wind turbines during the black start period of a wind power plant. The main content of this study is organized as follows. In <xref ref-type="sec" rid="s2">Section 2</xref>, the main idea of this study and the topology of a wind turbine system during black start are briefly introduced. In <xref ref-type="sec" rid="s3">Section 3</xref>, the black start control strategies of power converters and the power balance method in a wind turbine are detailed. In <xref ref-type="sec" rid="s4">Section 4</xref>, 3 sets of simulations are carried out to testify to the validity of the proposed control strategy. In <xref ref-type="sec" rid="s5">Section 5</xref>, the conclusion of the proposed method in this study is given and discussed.</p>
</sec>
<sec id="s2">
<title>2 System Description</title>
<p>When a blackout happens, the line-side voltage falls to zero. The wind turbine needs to restart with no grid voltage reference, and the power from the generator needs to be consumed with the local load. The main topology of the wind turbine based on PMSG during the black start procedure is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>, the line-side converter is connected to the fixed load, and there is also an extra crowbar circuit within the dc-link.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Topology of the wind turbine during a black start.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g001.tif"/>
</fig>
<p>Since there is no grid voltage during the whole black start procedure, the traditional MPPT control strategy of the full power converter is no longer suitable, the dc-link voltage can be controlled by the generator-side converter, and the line-side converter can be used to control the frequency, phase, and amplitude of the line-side voltage.</p>
<p>The local three-phase fixed load is designed for both power consumption and voltage establishment. With the exact information about the impedance characteristic of load and target amplitude and frequency of line-side voltage, the power demand of line-side voltage establishment can be calculated, which can be used to control the mechanical power capture of the wind turbine. During the black start, the wind turbine cannot operate at full power because the extra power cannot be absorbed by the grid. The power balance between the generator and load must be achieved by pitch angle control, as the pitch angle is the only controllable variable left to adjust mechanical power during black start. When the power balance is achieved steadily, the voltage of the line-side will be established and maintained at the designed value.</p>
<p>However, considering the power loss within the system, the power command calculated based on load and target line-side voltage cannot match the actual power need, which will cause the fluctuation in the rotor speed. In this study, the power command is designed higher than the calculated power, and the extra power is absorbed using the designed crowbar circuit to maintain the stability of rotor speed in a certain range.</p>
</sec>
<sec id="s3">
<title>3 Modeling and Control</title>
<p>In this section, the detailed model and control strategy of both the generator-side converter, the line-side converter, the wind turbine power capture method, and the crowbar circuit for a black start is presented.</p>
<sec id="s3-1">
<title>3.1 System Modeling</title>
<p>As in <xref ref-type="fig" rid="F1">Figure 1</xref>, the whole system contains a PMSG, a back-to-back full power converter, a crowbar circuit, an LCL filter, and a fixed load. The mathematical model of the system can be described based on <xref ref-type="fig" rid="F1">Figure 1</xref>. The mechanical model of a PMSG-based wind turbine can be described as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi mathvariant="italic">&#x3c1;A</mml:mi>
<mml:msubsup>
<mml:mi mathvariant="italic">V</mml:mi>
<mml:mi mathvariant="italic">w</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mi mathvariant="italic">p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="italic">&#x3bb;,&#x3b2;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mi mathvariant="italic">p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="italic">&#x3bb;</mml:mi>
<mml:mi mathvariant="normal">,</mml:mi>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mi mathvariant="italic">&#x3bb;</mml:mi>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi mathvariant="italic">e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3bb;</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.0068</mml:mn>
<mml:mi mathvariant="italic">&#x3bb;,</mml:mi>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3bb;</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mi mathvariant="italic">&#x3bb;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.008</mml:mn>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>0.035</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3bb;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">V</mml:mi>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e1">Eqs 1</xref>&#x2013;<xref ref-type="disp-formula" rid="e4">4</xref> (<xref ref-type="bibr" rid="B8">Li et al., 2012</xref>; <xref ref-type="bibr" rid="B7">Fathabadi, 2017</xref>), the variable <inline-formula id="inf1">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the mechanical power input, <inline-formula id="inf2">
<mml:math id="m6">
<mml:mi>&#x3c1;</mml:mi>
</mml:math>
</inline-formula> is the air density, <inline-formula id="inf3">
<mml:math id="m7">
<mml:mi>A</mml:mi>
</mml:math>
</inline-formula>is the rotation area of a wind turbine, <inline-formula id="inf4">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the wind speed, <inline-formula id="inf5">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the wind power coefficient, <inline-formula id="inf6">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the rotating speed of wind turbine and also can represent the rotor speed of PMSG, and <inline-formula id="inf7">
<mml:math id="m11">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula> is the tip speed ratio and <inline-formula id="inf8">
<mml:math id="m12">
<mml:mi>&#x3b2;</mml:mi>
</mml:math>
</inline-formula> is the pitch angle. The constant values from <inline-formula id="inf9">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> to <inline-formula id="inf10">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are 0.5173, 116, 0.4, 5, and 21.</p>
<p>As described in <xref ref-type="disp-formula" rid="e1">equations (1)</xref>-<xref ref-type="disp-formula" rid="e4">(4)</xref>, it is clear that wind power capture is related to the pitch angle and rotor speed.</p>
<p>The power flow within the wind power system can be written as<disp-formula id="e5">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="italic">J</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">dc</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">dc</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">loss</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e5">Equation 5</xref>, <inline-formula id="inf11">
<mml:math id="m18">
<mml:mi>J</mml:mi>
</mml:math>
</inline-formula> is the inertial coefficient of PMSG.</p>
<p>During the black start, when dc-link voltage is stable, the electrical power <inline-formula id="inf12">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> consists of load power <inline-formula id="inf13">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and power loss <inline-formula id="inf14">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. The power unbalance between <inline-formula id="inf15">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf16">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be reflected on <inline-formula id="inf17">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
<sec id="s3-2">
<title>3.2 Control of Back-To-Back Converter</title>
<p>The stability of dc-link voltage is quite important during the operation of the wind turbine. As is mentioned above, the line-side converter is not capable of maintaining the dc-link voltage during the initial transient non-voltage state. Instead, the generator-side converter can be used to control the dc-link voltage, the current equation of PMSG is represented as follows (<xref ref-type="bibr" rid="B8">Li et al., 2012</xref>):<disp-formula id="e8">
<mml:math id="m25">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">e</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mi mathvariant="italic">q</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">q</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">u</mml:mi>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mi mathvariant="italic">q</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">q</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">q</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">e</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">e</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c6;</mml:mi>
<mml:mi mathvariant="italic">f</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">u</mml:mi>
<mml:mi mathvariant="italic">q</mml:mi>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e8">Equation 8</xref>, <inline-formula id="inf18">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf19">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the d-q axis inductance of PMSG, <inline-formula id="inf20">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf21">
<mml:math id="m29">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the d-q axis current component of PMSG, <inline-formula id="inf22">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the electrical angle frequency, <inline-formula id="inf23">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf24">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the d-q axis terminal voltage component , and <inline-formula id="inf25">
<mml:math id="m33">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c6;</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the permanent magnet flux linkage of PMSG. Based on <xref ref-type="disp-formula" rid="e8">Equation 8</xref>, the control scheme of the generator-side converter is given in <xref ref-type="fig" rid="F2">Figure 2</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Control strategy of the generator-side converter.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g002.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F2">Figure 2</xref>, the generator-side converter uses a dc-link voltage outer loop and current inner loop to achieve the stability of dc-link voltage, the d-axis current component is controlled to 0 to achieve the maximum torque. The transfer diagram of the generator-side converter is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. In <xref ref-type="fig" rid="F3">Figure 3</xref>, <inline-formula id="inf26">
<mml:math id="m34">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf27">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represent the gain of the PI controller, <inline-formula id="inf28">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the gain of the generator-side converter which is normally regarded as an inertial link. <inline-formula id="inf29">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the stator winding, <inline-formula id="inf30">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the effect of permanent magnet flux linkage, and <inline-formula id="inf31">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>6</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the inertia coefficient of PMSG and wind turbine.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Transfer diagram of the generator-side converter.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g003.tif"/>
</fig>
<p>The dc-link voltage in <xref ref-type="fig" rid="F1">Figure 1</xref> can be described as<disp-formula id="e9">
<mml:math id="m40">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">dc</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">c</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">dc</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">cb</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>From <xref ref-type="disp-formula" rid="e9">Equation 9</xref> it is clear that the changes in dc-link voltage are related to the control of line-side and crowbar control. When the dc-link voltage is relatively stable, the power generated can flow into the line-side, and the power flow can be described using <xref ref-type="disp-formula" rid="e5">Eqs (5)</xref> and <xref ref-type="disp-formula" rid="e6">(6)</xref>.</p>
<p>The line-side converter is used to control the frequency and amplitude of output voltage, offering voltage reference for the other turbines. The current equation is given as follows (<xref ref-type="bibr" rid="B2">Ashourianjozdani et al., 2018</xref>):<disp-formula id="e10">
<mml:math id="m41">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gd</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gd</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gd</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gq</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gq</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">u</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gd</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">e</mml:mi>
<mml:mi mathvariant="italic">d</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gq</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gq</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gq</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gd</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gd</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">u</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">gq</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">e</mml:mi>
<mml:mi mathvariant="italic">q</mml:mi>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e10">Equation 10</xref>, <inline-formula id="inf32">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf33">
<mml:math id="m43">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the d-q axis line-side inductance, <inline-formula id="inf34">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf35">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the d-q axis component of line-side output current, <inline-formula id="inf36">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the line-side frequency which is controlled to 50Hz, <inline-formula id="inf37">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf38">
<mml:math id="m48">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the d-q axis component of the line-side terminal voltage, and the line-side voltage source <inline-formula id="inf39">
<mml:math id="m49">
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf40">
<mml:math id="m50">
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are set to 0 because there is no grid voltage during the whole black start procedure.</p>
<p>Based on <xref ref-type="disp-formula" rid="e10">Equation 10</xref>, the control diagram can be obtained in <xref ref-type="fig" rid="F4">Figure 4</xref>. As in <xref ref-type="fig" rid="F4">Figure 4</xref>, the line-side converter uses phase-to-ground voltage amplitude outer loop and current inner loop, the q-axis component is set to zero to improve the power factor. The transfer diagram of the line-side converter is shown in <xref ref-type="fig" rid="F5">Figure 5</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Control strategy of the line-side converter.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Transfer diagram of the line-side converter.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g005.tif"/>
</fig>
<p>In <xref ref-type="fig" rid="F5">Figure 5</xref>, <inline-formula id="inf41">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf42">
<mml:math id="m52">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>8</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represent the gain of PI controller, <inline-formula id="inf43">
<mml:math id="m53">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>9</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the transfer function of the grid-side converter which is also treated as an inertia link, <inline-formula id="inf44">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the power loss in line before the capacitor branch of LCL filter, and <inline-formula id="inf45">
<mml:math id="m55">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mn>11</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the fixed load designed for voltage establishment.</p>
<p>Ignoring the voltage of the filter capacitor, the transfer function of <xref ref-type="fig" rid="F5">Figure 5</xref> can be obtained as follows:<disp-formula id="e11">
<mml:math id="m56">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">amref</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">am</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mn>8</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mn>9</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mrow>
<mml:mn>11</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mn>8</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mn>9</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mn>8</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mn>9</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">G</mml:mi>
<mml:mrow>
<mml:mn>11</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>
<xref ref-type="disp-formula" rid="e11">Eq. 11</xref> can be regarded as a phase shift, and the amplitude of load voltage will be controlled as the command value. Because there is no line voltage reference, the PLL is no longer needed, and the phase signal <inline-formula id="inf46">
<mml:math id="m57">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> which is used to convert the coordinate is calculated based on the target frequency.</p>
</sec>
<sec id="s3-3">
<title>3.3 Voltage Establishment and Power Balance Control</title>
<p>With the control structure built by the back-to-back converter, the amplitude and frequency of voltage on load can be controlled to a certain designed value, which can offer voltage reference during the black start to multiple wind turbines. In order to maintain the voltage level on load, the mechanical power input must be high enough to make sure the power flow within the converter can satisfy the power demand on load. When the amplitude and frequency of load voltage are determined, the mathematical relationship between voltage, load, and three-phase power on load can be represented as<disp-formula id="e12">
<mml:math id="m58">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">S</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:msubsup>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">amp</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">Z</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e12">Equation 12</xref>, <inline-formula id="inf47">
<mml:math id="m59">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the apparent power on load, <inline-formula id="inf48">
<mml:math id="m60">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the phase-to-ground voltage amplitude on load, and <inline-formula id="inf49">
<mml:math id="m61">
<mml:mrow>
<mml:msub>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the impedance of the load. With a fixed resistive load, the power under a certain voltage level can be calculated as<disp-formula id="e13">
<mml:math id="m62">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:msubsup>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">amp</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>Also, the space state equation of the line-side in <xref ref-type="fig" rid="F4">Figure 4</xref> can be obtained as<disp-formula id="e14">
<mml:math id="m63">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:mi mathvariant="italic">L</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">FC</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">FC</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">FC</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mi mathvariant="italic">F</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">FC</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e14">Equation 14</xref>, <inline-formula id="inf50">
<mml:math id="m64">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf51">
<mml:math id="m65">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represent the <inline-formula id="inf52">
<mml:math id="m66">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf53">
<mml:math id="m67">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> in <xref ref-type="fig" rid="F4">Figure 4</xref>, the power loss of the line-side can be obtained from <xref ref-type="disp-formula" rid="e14">Equation 14</xref> as<disp-formula id="e15">
<mml:math id="m68">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">lossL</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mi mathvariant="italic">g</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi mathvariant="italic">R</mml:mi>
<mml:mi mathvariant="italic">L</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">L</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">FC</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mi mathvariant="italic">F</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">FC</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">dt</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
</p>
<p>Unlike power loss of line-side, power loss that exists in the winding of PMSG is hard to quantify. In order to simplify the analysis in this study, the power loss in PMSG is symbolized by <inline-formula id="inf54">
<mml:math id="m69">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, therefore the power loss and electrical power can be described using <xref ref-type="disp-formula" rid="e16">the</xref> following equation:<disp-formula id="e16">
<mml:math id="m70">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">loss</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">lossL</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">lossG</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">loss</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">load</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
</p>
<p>The constant voltage on load means stable power consumption on the line-side. The load voltage can maintain stability as long as <inline-formula id="inf55">
<mml:math id="m71">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be supplied by the generator. The line-side voltage can be supported when <inline-formula id="inf56">
<mml:math id="m72">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> meets the need of <inline-formula id="inf57">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> in <xref ref-type="disp-formula" rid="e16">Equation 16</xref>; otherwise, the rotor speed of PMSG will change because of the power unbalance until the system loses its control stability.</p>
<p>In order to make sure that <inline-formula id="inf58">
<mml:math id="m74">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> meets <inline-formula id="inf59">
<mml:math id="m75">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the power capture of a wind turbine is the only controllable link left if there is no power storage. As in <xref ref-type="disp-formula" rid="e1">Eqs 1</xref>, <xref ref-type="disp-formula" rid="e2">2</xref>, <inline-formula id="inf60">
<mml:math id="m76">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be controlled by pitch angle control and rotor speed control. If the power loss in PMSG can be ignored, the ideal power balance can be described with the following:<disp-formula id="e17">
<mml:math id="m77">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">mi</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>A</mml:mi>
<mml:msubsup>
<mml:mi mathvariant="italic">V</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">wi</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi mathvariant="italic">C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">pi</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3bb;</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:mi mathvariant="normal">,</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">ei</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">lossi</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">loadi</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">mi</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">ei</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e17">Equation 17</xref>, the letter subscript <inline-formula id="inf61">
<mml:math id="m78">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> means the ideal situation. Under the ideal situation, the balance between <inline-formula id="inf62">
<mml:math id="m79">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf63">
<mml:math id="m80">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is achieved by pitch control of the wind turbine, the wind power coefficient <inline-formula id="inf64">
<mml:math id="m81">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is used to describe <inline-formula id="inf65">
<mml:math id="m82">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> based on the wind speed <inline-formula id="inf66">
<mml:math id="m83">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, initial rotor speed <inline-formula id="inf67">
<mml:math id="m84">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and the data in the look-up table in the theoretical analysis. The rotor speed <inline-formula id="inf68">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> will maintain steady because there is no delay in pitch control under an ideal situation, and the power balance can be achieved as long as the system operates.</p>
<p>However, considering the power loss in PMSG and the delay in pitch control, it is difficult to maintain the power balance as the ideal situation, as a result <inline-formula id="inf69">
<mml:math id="m86">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> will fluctuate with power unbalance, damaging the stability of the system. In case of that, a power balance control strategy using pitch control and crowbar circuit to maintain the relative stability of rotor speed is proposed.</p>
<p>When wind speed reaches the cut-in speed, the brake of the wind turbine should be released to make sure there is enough rotational kinetic energy saved in the rotor, otherwise, the mechanical power will be smaller than electrical power, and the capacitor charging and load voltage establishment will not be completed. After the initial rotor kinetic energy satisfies the electrical power need, pitch angle control starts to react to the power command to make sure wind power captured in <xref ref-type="disp-formula" rid="e1">Equation 1</xref> is higher than <inline-formula id="inf70">
<mml:math id="m87">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> in <xref ref-type="disp-formula" rid="e17">Equation 17</xref>. The updated <inline-formula id="inf71">
<mml:math id="m88">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> value is obtained based on the transient real-time rotor speed, tip speed ratio, and look-up table.</p>
<p>Based on <xref ref-type="disp-formula" rid="e5">Equation 5</xref>, the rotor speed will increase because of the power unbalance. A designed crowbar circuit is used to consume the extra power, the power consumed in the crowbar circuit can be described using <xref ref-type="disp-formula" rid="e18">the</xref> following equation:<disp-formula id="e18">
<mml:math id="m89">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">cb</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="italic">U</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">dc</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">Z</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">cb</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">cb</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>o</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>the power of the crowbar circuit must satisfy the following constraint:<disp-formula id="e19">
<mml:math id="m90">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">lossi</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">loadi</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">cb</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">lossi</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">loadi</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>
</p>
<p>When the power flow satisfies <xref ref-type="disp-formula" rid="e19">Equation 19</xref>, the rotor speed can be stabled within a certain range. The control of the crowbar circuit based on the detected rotor speed is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. In <xref ref-type="fig" rid="F6">Figure 6</xref>, <inline-formula id="inf72">
<mml:math id="m91">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf73">
<mml:math id="m92">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the threshold values.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Control strategy of the crowbar circuit.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g006.tif"/>
</fig>
<p>When there is a sudden disturbance in the system, the pitch angle control reacts to the power disturbance at first. When the mechanical power input is adjusted by pitch control to balance with the electrical power, the crowbar circuit is then activated to adjust the rotor speed after the power balance is relatively accomplished. The change in rotor speed caused by power imbalance during power regulation by pitch control can be compensated by a crowbar. In other words, the crowbar circuit works as a secondary regulator. However, the crowbar circuit in this work is necessary because the change in rotor speed during the power regulation by pitch control must be regulated so that the control stability issues can be avoided.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Results</title>
<sec id="s4-1">
<title>4.1 Case 1: Wind Speed Fixed Situation</title>
<p>In this section, the simulation results in the wind turbine system shown in <xref ref-type="fig" rid="F1">Figure 1</xref> are given. The designed phase-to-phase line-side voltage RMS is 690&#xa0;V, and the load value is fixed so that the power on the load changes with the designed load voltage. The results are obtained using the MATLAB/Simulink.</p>
<p>First, the dc-link voltage is set to 1100&#xa0;V at t &#x3d; 0s, 1200&#xa0;V at t &#x3d; 3&#xa0;s, then back to 1100&#xa0;V at t &#x3d; 6&#xa0;s to verify the control of generator-side. The result is shown in <xref ref-type="fig" rid="F7">Figure 7</xref>.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Change of the dc-link voltage.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g007.tif"/>
</fig>
<p>As in <xref ref-type="fig" rid="F7">Figure 7</xref>, the dc-link voltage can follow the control reference as expected. The wind speed is set to 10&#xa0;m/s, the initial rotor speed of the wind turbine is set to 1.5&#xa0;rad/s because in reality the brake of the wind turbine should be released and the initial rotor kinetic energy should be high enough to make sure the PMSG will not stall against the electrical power and lose the operation stability during black start. The dc-link voltage, line-side voltage, line-side frequency, and power are shown in <xref ref-type="fig" rid="F8">Figure 8</xref>.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>The dc-link voltage, line-side voltage, and frequency in case 1. <bold>(A)</bold> The dc-link voltage <bold>(B)</bold> Detailed three-phase line-side voltage <bold>(C)</bold> Line-side voltage amplitude <bold>(D)</bold> line-side frequency.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g008.tif"/>
</fig>
<p>As in <xref ref-type="fig" rid="F8">Figure 8A</xref>, the dc-link voltage is controlled to 1100&#xa0;V, and the specifics of the line-side voltage are presented in <xref ref-type="fig" rid="F8">Figure 8B&#x2013;D</xref>. As in <xref ref-type="fig" rid="F8">Figure 8C</xref>, the phase-to-ground voltage amplitude of load voltage is controlled to 563&#xa0;V, and the frequency in <xref ref-type="fig" rid="F8">Figure 8D</xref> is maintained as 50&#xa0;Hz. When the crowbar circuit is activated, the power flow satisfies <xref ref-type="disp-formula" rid="e11">Equation 11</xref> which causes the periodic fluctuation in electrical power.</p>
<p>The changes in the rotor speed and power are given in <xref ref-type="fig" rid="F9">Figures 9A,B</xref>. During t &#x3d; 0&#x2013;0.3&#xa0;s, <inline-formula id="inf74">
<mml:math id="m93">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> rises to a transient peak level and returns to stability very quickly because the line-side voltage and dc-link voltage are under control in a short time. The rotor speed increases at first and then decreases because the initial rotor speed of 1.5&#xa0;rad/s is at the right side of the maximum power point so that <inline-formula id="inf75">
<mml:math id="m94">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> decreases when rotor speed increases. When <inline-formula id="inf76">
<mml:math id="m95">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> falls lower than<inline-formula id="inf77">
<mml:math id="m96">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the rotor speed <inline-formula id="inf78">
<mml:math id="m97">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> will decrease until it is lower than the threshold values<inline-formula id="inf79">
<mml:math id="m98">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. Then the crowbar circuit will start to work periodically. The status of the crowbar circuit is presented in <xref ref-type="fig" rid="F9">Figure 9C</xref>. Because the rotor speed is higher than the trigger threshold value<inline-formula id="inf80">
<mml:math id="m99">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the crowbar IGBT is kept open at the beginning till the rotor speed decreases to<inline-formula id="inf81">
<mml:math id="m100">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, which explains the reduction in <inline-formula id="inf82">
<mml:math id="m101">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> at about 1.9&#xa0;s.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>The rotor speed, electrical power, and crowbar status in case 1. <bold>(A)</bold> The rotor speed <bold>(B)</bold> Electrical power Pe <bold>(C)</bold> Crowbar status.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g009.tif"/>
</fig>
</sec>
<sec id="s4-2">
<title>4.2 Case 2: Variable Line-Side Voltage</title>
<p>In this case, the reference value of phase-to-ground voltage amplitude of load is set to 563&#xa0;V at t &#x3d; 0&#xa0;s, 600V at t &#x3d; 5&#xa0;s, and 500&#xa0;V at t &#x3d; 10&#xa0;s. However, the mechanical power command is not adjusted with the change of line-side voltage so that the system dynamics can be shown more independently. The status of load phase-to-ground voltage is shown in <xref ref-type="fig" rid="F10">Figure 10A</xref>, the changes in rotor speed and power are shown in <xref ref-type="fig" rid="F10">Figures 10B,C</xref>, and the status of the crowbar circuit is given in <xref ref-type="fig" rid="F10">Figure 10D</xref>. The crowbar circuit is not activated until about t &#x3d; 5&#xa0;s when the rotor speed is higher than the threshold. Then, <inline-formula id="inf83">
<mml:math id="m102">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is higher than <inline-formula id="inf84">
<mml:math id="m103">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> after t &#x3d; 5&#xa0;s because the line-side voltage rises, which makes the rotor speed continue to drop and the crowbar circuit, therefore, is deactivated. At t &#x3d; 10&#xa0;s, the line-side voltage drops and <inline-formula id="inf85">
<mml:math id="m104">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> decreases, which makes the rotor speed rises till the crowbar circuit is activated periodically.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>The line-side voltage amplitude, rotor speed, electrical power, and crowbar status in case 2. <bold>(A)</bold> Line-side voltage amplitude <bold>(B)</bold> Rotor speed <bold>(C)</bold> Electrical power Pe <bold>(D)</bold> Crowbar status.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g010.tif"/>
</fig>
</sec>
<sec id="s4-3">
<title>4.3 Case 3: Change of Wind Speed</title>
<p>In this case, the wind speed changes from 10&#xa0;m/s to 7&#xa0;m/s at t &#x3d; 3&#xa0;s, and then changes back to 10&#xa0;m/s at t &#x3d; 6&#xa0;s, the change of rotor speed is shown in <xref ref-type="fig" rid="F11">Figure 11</xref>. As in <xref ref-type="fig" rid="F11">Figure 11</xref>, when wind speed changes at t &#x3d; 3&#xa0;s and t &#x3d; 6&#xa0;s, the delay of pitch action causes a transient power gap between <inline-formula id="inf86">
<mml:math id="m105">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and<inline-formula id="inf87">
<mml:math id="m106">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, thus rotor speed changes with power. After wind speed returns to 10&#xa0;m/s, the crowbar circuit can still work as a rotor speed stabilizer.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Change of the rotor speed in case 3.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g011.tif"/>
</fig>
<p>The dc-link voltage, line-side voltage, and power are presented in <xref ref-type="fig" rid="F12">Figure 12</xref>. In <xref ref-type="fig" rid="F12">Figure 12</xref>, it is clear that instead of changes in wind speed, the changes in <inline-formula id="inf88">
<mml:math id="m107">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are directly related to the changes in <inline-formula id="inf89">
<mml:math id="m108">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf90">
<mml:math id="m109">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, which improves that as long as the electrical power can be offered by the wind power input and rotor kinetic energy storage, the voltage of dc-link and line-side can be maintained by converter without being affected by change of rotor speed. Instead, the rotor speed can be treated as a variable of the state which reflects the power unbalance between <inline-formula id="inf91">
<mml:math id="m110">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf92">
<mml:math id="m111">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>The DC-link voltage, line-side voltage, power, and crowbar status in case 3. <bold>(A)</bold> The dc-link voltage <bold>(B)</bold> Detailed three-phase line-side voltage <bold>(C)</bold> Electrical power Pe <bold>(D)</bold> Crowbar status.</p>
</caption>
<graphic xlink:href="fenrg-10-901708-g012.tif"/>
</fig>
</sec>
</sec>
<sec id="s5">
<title>5 Conclusion</title>
<p>In this study, a control strategy of the black start and line-side voltage establishment for PMSG-based wind power system is proposed. The dc-link voltage is stabled using the generator-side converter, and the line-side voltage is controlled by the line-side converter. The pitch angle is used to balance mechanical power and electrical power, and an extra crowbar circuit is applied to maintain the rotor speed. The proposed black start control strategy and the line-side voltage establishment method should be carried out in the following order:</p>
<p>
<statement content-type="step" id="Step_1">
<label>Step 1</label>
<p>Calculating the load power based on the desired line-side voltage level and local load configuration.</p>
</statement>
</p>
<p>
<statement content-type="step" id="list0">
<label>Step 2</label>
<p>Releasing the brake device of the wind turbine when wind speed reaches the cut-in speed, initiating the controller of the generator-side converter, the line-side converter, and the crowbar circuit after the rotor speed reaches a point under which the initial rotor kinetic energy can match the load power need in step 1.</p>
</statement>
</p>
<p>
<statement content-type="step" id="list1">
<label>Step 3</label>
<p>: Adjusting the pitch angle using a power closed loop to make sure the mechanical power input and the electrical power demand match the description in <xref ref-type="disp-formula" rid="e19">Eq. 19</xref>. At this point, the rotor speed will rise and then be controlled periodically by the crowbar circuit.</p>
</statement>
</p>
<p>
<statement content-type="step" id="list2">
<label>Step 4</label>
<p>When the line-side voltage reaches the designed value, the other wind turbines can follow up their procedures using the line-side voltage of the first wind turbine as a reference, which will be testified in our next work.</p>
<p>By the proposed control strategy, the problems about the stability of the dc-link voltage and the establishment of line-side voltage can be solved. The simulation results show the validity of the proposed method, and rotor speed can be used to check the unbalance between the electrical power and mechanical power. Once the line-side voltage is established, more wind turbines can start up using the voltage reference built by the one using the proposed method in this study.</p>
</statement>
</p>
</sec>
</body>
<back>
<sec id="s6">
<title>Data Availability Statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>XJ: designing of the main control strategy. DL: modifying of the topology. PH: construction of simulation. KC: designing of the simulation control group. KJ: writing of the original manuscript. GH: revision of the manuscript. SH: optimization of the control strategy.</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of Interest</title>
<p>Author XJ is employed by State Grid Hubei Electric Power CO., LTD.; Authors DL, PH, KC, and KJ are employed by Electric Power Research Institute, Hubei Electric Power Company, State Grid; authors GH and SH are employed by the College of Electrical and Information Engineering, Hunan University. This study received funding from the science and technology project of the State Grid Corporation of China, project number 4000-202122070A-0-0-00. The funder had involvement in the study design, collection, and data analysis. All authors declare no other competing interests.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alavi</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Ghazi</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A Novel Control Strategy Based on a Look-Up Table for Optimal Operation of MTDC Systems in Post-contingency Conditions</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>7</volume>, <fpage>4</fpage>. <pub-id pub-id-type="doi">10.1186/s41601-022-00224-3</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ashourianjozdani</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lopes</surname>
<given-names>L. A. C.</given-names>
</name>
<name>
<surname>Pillay</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Power Electronic Converter Based PMSG Emulator: A Testbed for Renewable Energy Experiments</article-title>. <source>IEEE Trans. Ind. Appl.</source> <volume>54</volume>, <fpage>3626</fpage>&#x2013;<lpage>3636</lpage>. <pub-id pub-id-type="doi">10.1109/TIA.2018.2819618</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benato</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Bruno</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Sessa</surname>
<given-names>S. D.</given-names>
</name>
<name>
<surname>Giannuzzi</surname>
<given-names>G. M.</given-names>
</name>
<name>
<surname>Ortolano</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Pedrazzoli</surname>
<given-names>G.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>A Novel Modeling for Assessing Frequency Behavior during a Hydro-To-Thermal Plant Black Start Restoration Test</article-title>. <source>IEEE Access</source> <volume>7</volume>, <fpage>47317</fpage>&#x2013;<lpage>47328</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2019.2909321</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chou</surname>
<given-names>Y.-T.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.-W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.-J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>C.-C.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>C.-C.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Development of a Black Start Decision Supporting System for Isolated Power Systems</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>28</volume>, <fpage>2202</fpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2013.2237792</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deng</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Control Strategy of Wind Turbine Based on Permanent Magnet Synchronous Generator and Energy Storage for Stand-Alone Systems</article-title>. <source>Chin. J. Electr. Eng.</source> <volume>3</volume>, <fpage>51</fpage>&#x2013;<lpage>62</lpage>. <pub-id pub-id-type="doi">10.23919/CJEE.2017.7961322</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Erdiwansyah</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mahidin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Husin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Nasaruddin</surname>
<given-names>fnm.</given-names>
</name>
<name>
<surname>Zaki</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Muhibbuddin</surname>
<given-names>fnm.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A Critical Review of the Integration of Renewable Energy Sources with Various Technologies</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>6</volume>, <fpage>3</fpage>. <pub-id pub-id-type="doi">10.1186/s41601-021-00181-3</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fathabadi</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Novel Maximum Electrical and Mechanical Power Tracking Controllers for Wind Energy Conversion Systems</article-title>. <source>IEEE J. Emerg. Sel. Top. Power Electron.</source> <volume>5</volume>, <fpage>1739</fpage>&#x2013;<lpage>1745</lpage>. <pub-id pub-id-type="doi">10.1109/JESTPE.2017.2727978</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Haskew</surname>
<given-names>T. A.</given-names>
</name>
<name>
<surname>Swatloski</surname>
<given-names>R. P.</given-names>
</name>
<name>
<surname>Gathings</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Optimal and Direct-Current Vector Control of Direct-Driven PMSG Wind Turbines</article-title>. <source>IEEE Trans. Power Electron.</source> <volume>27</volume>, <fpage>2325</fpage>&#x2013;<lpage>2337</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2011.2174254</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lindstrom</surname>
<given-names>R. R.</given-names>
</name>
</person-group> (<year>1990</year>). <article-title>Simulation and Field Tests of the Black Start of A Large Coal-Fired Generating Station Utilizing Small Remote Hydro Generation</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>5</volume>, <fpage>162</fpage>&#x2013;<lpage>168</lpage>. <pub-id pub-id-type="doi">10.1109/59.49101</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>T. E.</given-names>
</name>
<name>
<surname>Schomer</surname>
<given-names>I.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Operational Optimization of a Building-Level Integrated Energy System Considering Additional Potential Benefits of Energy Storage</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>6</volume>, <fpage>4</fpage>. <pub-id pub-id-type="doi">10.1186/s41601-021-00184-0</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pape</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kazerani</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Turbine Startup and Shutdown in Wind Farms Featuring Partial Power Processing Converters</article-title>. <source>IEEE Open J. Power Energy</source> <volume>7</volume>, <fpage>254</fpage>&#x2013;<lpage>264</lpage>. <pub-id pub-id-type="doi">10.1109/OAJPE.2020.3006352</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodriguez-Amenedo</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>Gomez</surname>
<given-names>S. A.</given-names>
</name>
<name>
<surname>Martinez</surname>
<given-names>J. C.</given-names>
</name>
<name>
<surname>Alonso-Martinez</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Black-Start Capability of DFIG Wind Turbines through a Grid-Forming Control Based on the Rotor Flux Orientation</article-title>. <source>IEEE Access</source> <volume>9</volume>, <fpage>142910</fpage>&#x2013;<lpage>142924</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2021.3120478</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Saborio-Romano</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Bidadfar</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sakamuri</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Goksu</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Cutululis</surname>
<given-names>N. A.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>Primary Frequency Response from Offshore Wind Farms Connected to HVdc via Diode Rectifiers</article-title>,&#x201d; in <conf-name>Proceeding of the 2019 IEEE Milan PowerTech</conf-name>, <conf-loc>Milan, Italy</conf-loc>, <conf-date>June 2019</conf-date> (<publisher-name>IEEE</publisher-name>). <pub-id pub-id-type="doi">10.1109/PTC.2019.8810907</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Sakamuri</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Bidadfar</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Saborio-Romano</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Jain</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Cutululis</surname>
<given-names>N. A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Black Start by HVdc-Connected Offshore Wind Power Plants</article-title>. <conf-name>Proceeding of the 45th Annual Conference of the IEEE Industrial Electronics Society</conf-name>, <conf-loc>Lisbon, Portugal</conf-loc>, <conf-date>Oct. 2019</conf-date> (<publisher-name>IEEE</publisher-name>). <pub-id pub-id-type="doi">10.1109/IECON.2019.8927615</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Satpathy</surname>
<given-names>A. S.</given-names>
</name>
<name>
<surname>Kishore</surname>
<given-names>N. K.</given-names>
</name>
<name>
<surname>Kastha</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Sahoo</surname>
<given-names>N. C.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Control Scheme for a Stand-Alone Wind Energy Conversion System</article-title>. <source>IEEE Trans. Energy Convers.</source> <volume>29</volume>, <fpage>418</fpage>&#x2013;<lpage>425</lpage>. <pub-id pub-id-type="doi">10.1109/TEC.2014.2303203</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Application of Type 3 Wind Turbines for System Restoration</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>33</volume>, <fpage>3040</fpage>&#x2013;<lpage>3051</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2017.2762009</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Frequency Control Strategy for Black Starts via PMSG-Based Wind Power Generation</article-title>. <source>Energies</source> <volume>10</volume>, <fpage>358</fpage>. <pub-id pub-id-type="doi">10.3390/en10030358</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Morrow</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Coordinated DC Voltage Control of Wind Turbine with Embedded Energy Storage System</article-title>. <source>IEEE Trans. Energy Convers.</source> <volume>27</volume>, <fpage>1036</fpage>&#x2013;<lpage>1045</lpage>. <pub-id pub-id-type="doi">10.1109/TEC.2012.2220361</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yan</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Uncertainty Modeling of Wind Power Frequency Regulation Potential Considering Distributed Characteristics of Forecast Errors</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>6</volume>, <fpage>22</fpage>. <pub-id pub-id-type="doi">10.1186/s41601-021-00200-3</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>