<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "https://jats.nlm.nih.gov/nlm-dtd/publishing/3.0/journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" xml:lang="en" dtd-version="3.0" article-type="research-article">
  <front>
    <journal-meta><journal-id journal-id-type="publisher">WES</journal-id><journal-title-group>
    <journal-title>Wind Energy Science</journal-title>
    <abbrev-journal-title abbrev-type="publisher">WES</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Wind Energ. Sci.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">2366-7451</issn><publisher>
    <publisher-name>Copernicus Publications</publisher-name>
    <publisher-loc>Göttingen, Germany</publisher-loc>
  </publisher></journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.5194/wes-11-2037-2026</article-id><title-group><article-title>Experimental validation and extension of a blade element momentum model for counter-rotation dual-rotor wind turbine with double rotational armature design</article-title><alt-title>Experimental validation and extension of a BEM model</alt-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Adema</surname><given-names>Niels Cornelis</given-names></name>
          <email>ni.c.adema@pl.hanze.nl</email>
        <ext-link>https://orcid.org/0000-0002-6592-8606</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1 aff2">
          <name><surname>Gil-Vernet Pagonabarraga</surname><given-names>Josep</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>Swart Ranshuysen</surname><given-names>Wouter</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>de Ruijter</surname><given-names>Arjen</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1 aff3">
          <name><surname>Schepers</surname><given-names>Gerard</given-names></name>
          
        </contrib>
        <aff id="aff1"><label>1</label><institution>Entrance – Centre of Expertise Energy part of Hanze UAS, Groningen, 9747AA, the Netherlands</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>Institute of Engineering, section Mechanical Engineering, Hanze UAS, Groningen, 9747AS, the Netherlands</institution>
        </aff>
        <aff id="aff3"><label>3</label><institution>TNO Energy Transition and Materials, 2628 CK, Delft, the Netherlands</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Niels Cornelis Adema (ni.c.adema@pl.hanze.nl)</corresp></author-notes><pub-date><day>10</day><month>June</month><year>2026</year></pub-date>
      
      <volume>11</volume>
      <issue>6</issue>
      <fpage>2037</fpage><lpage>2051</lpage>
      <history>
        <date date-type="received"><day>9</day><month>November</month><year>2025</year></date>
           <date date-type="rev-request"><day>19</day><month>November</month><year>2025</year></date>
           <date date-type="rev-recd"><day>18</day><month>April</month><year>2026</year></date>
           <date date-type="accepted"><day>25</day><month>May</month><year>2026</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2026 Niels Cornelis Adema et al.</copyright-statement>
        <copyright-year>2026</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026.html">This article is available from https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026.html</self-uri><self-uri xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026.pdf">The full text article is available as a PDF file from https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026.pdf</self-uri>
      <abstract><title>Abstract</title>

      <p id="d2e130">Small wind turbines face significant challenges in achieving commercial viability due to lower efficiency and higher energy costs compared to utility-scale systems and competing renewable technologies. Counter-rotating dual-rotor wind turbines (CR-DRWTs) with dual-rotational armature configurations offer a potential pathway for efficiency improvements through doubled direct-drive power and minimal mechanical complexity, suitable for urban applications. This study presents a detailed experimental investigation of a 1.6 m diameter CR-DRWT through wind tunnel testing at the Centre Scientifique et Technique du Bâtiment (CSTB) in Nantes, France, conducted at wind speeds between 4 and 15 m s<sup>−1</sup>. An improved turbine design is tested with enhanced instrumentation, including independent measurements of rotor rotational speed and blade pitch angle, enabling a detailed characterization of rotor interaction, operating behaviour, and power performance. The turbine achieved a maximum electrical power output of 1014 W and a peak measured power coefficient of 0.33 while demonstrating reliable self-starting at wind speeds as low as 3.5 m s<sup>−1</sup>. To interpret and generalize the experimental findings, an existing blade element momentum (BEM) model for dual-rotor systems is extended to explicitly represent torque balance and power production in a dual-rotational armature configuration. The extended model shows good agreement with experimental trends and is further applied in a sweep optimization. The optimized configuration predicts a theoretical maximum power coefficient of 0.56, highlighting substantial remaining performance potential. By combining wind tunnel measurements with a validated BEM model extension, this work provides unique reference data and supports the development of mechanically simple CR-DRWTs for future small-scale wind energy systems.</p>
  </abstract>
    
<funding-group>
<award-group id="gs1">
<funding-source>Nationaal Regieorgaan Praktijkgericht Onderzoek SIA</funding-source>
<award-id>RAAK.MKB17.007</award-id>
</award-group>
</funding-group>
</article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d2e166">While modern wind turbines have become the largest rotating machines on earth with further upscaling planned, renewed interest in small wind turbines (SWTs) is fostered through local energy transition and smart grid development. SWTs have traditionally lacked the aerodynamic refinement of larger turbines, resulting in lower efficiency, lower capacity factors, and higher energy costs  (Bianchini et al., 2022). Still, SWTs serve diverse applications worldwide, including power generation for households, industrial centres, farms, and isolated communities; hybrid energy systems for remote monitoring and telecommunications; and direct energy services like water pumping, desalination, and purification (Chagas et al., 2020). While SWTs can outperform photovoltaic (PV) systems in annual power generation at specific locations, investments require site-specific wind resource assessments, and support schemes must avoid subsidizing low-potential areas (Jurasz et al., 2025). Studies have also shown potential for urban wind applications (Bereziartua-Gonzalez et al., 2025; Calautit and Johnstone, 2023). Despite these opportunities, SWTs face significant economic challenges. Their energy costs typically exceed both residential electricity prices and those of competing technologies like rooftop solar and utility-scale wind farms. These elevated costs stem from limited development compared to large-scale systems, as well as disproportionately high expenses for electrical connection, resource assessment, and installation (Simic et al., 2013). To achieve commercial viability, SWTs must either achieve substantial cost reductions, for example, through targeted policy incentives (Jurasz et al., 2025), or must significantly improve their energy capture capability (Bianchini et al., 2022).</p>
      <p id="d2e169">Dual-rotor wind turbine (DRWT) configurations represent one potential pathway for efficiency gains, though this approach introduces trade-offs such as higher manufacturing costs, increased structural weight, and greater mechanical complexity. DRWTs can be categorized into co-rotating systems (CO-DRWTs), where both rotors rotate in the same clockwise direction, or counter-rotating systems (CR-DRWTs) where one rotor rotates clockwise and the other counter-clockwise. Several mechanical configurations exist for dual-rotor wind turbines (DRWTs). Multiple experiments are performed on individual wind turbines placed close together or a set-up with two rotors, each with a separate generator. To design a complete single DRWT system, one approach employs a bevel gear system that connects both rotors to a single generator shaft (Jung et al., 2005; Schepers et al., 2024). Alternatively, a double rotational armature generator features one rotor connected to the generator stator and the other to the generator rotor (Adema et al., 2025; Booker et al., 2010; Mitulet et al., 2015). This latter configuration offers notable advantages, including doubled direct-drive power, minimal starting torque, high electrical and mechanical efficiency, compact design, and reduced noise and vibration – characteristics that make it particularly suitable for small-scale urban wind applications.</p>
      <p id="d2e172">Newman (1986) extended the classical 1D Betz theory for single-actuator discs to configurations with multiple discs, demonstrating that a dual-rotor system can theoretically achieve a power coefficient (CP) of <inline-formula><mml:math id="M3" display="inline"><mml:mrow><mml:mn mathvariant="normal">8</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">75</mml:mn></mml:mrow></mml:math></inline-formula> <inline-formula><mml:math id="M4" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 0.64, with an induction factor of <inline-formula><mml:math id="M5" display="inline"><mml:mrow><mml:mn mathvariant="normal">3</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">5</mml:mn></mml:mrow></mml:math></inline-formula> at the second rotor. These values exceed those of a single rotor operating at the Betz limit (CP <inline-formula><mml:math id="M6" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 0.59, induction factor <inline-formula><mml:math id="M7" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> <inline-formula><mml:math id="M8" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>/</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:math></inline-formula>). Through smoke visualization experiments in a wind tunnel using porous discs at varying separation distances, Newman recommended a minimum spacing of 0.5 <inline-formula><mml:math id="M9" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> (where <inline-formula><mml:math id="M10" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> is the rotor diameter) between discs to minimize flow curvature and non-uniformity, thereby maintaining the validity of the 1D assumptions. More recently, Sundararaju et al. (2017) investigated CR-DRWTs with equal-diameter rotors. They found that the maximum achievable CP reaches 0.814 when the rotors are separated by an axial distance of 2.8 <inline-formula><mml:math id="M11" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula>. Their work also revealed that reducing the axial spacing toward zero causes the power coefficient to decrease progressively toward the single-rotor Betz limit.</p>
      <p id="d2e254">A limited amount of field- and lab-scale testing on DRWTs has been performed in the last few years. A field test by Jung et al. (2005) showed a CP of approximately 0.5 for a 30 kW CR-DRWT with asymmetric rotors (5.5 m front, 11 m rear) connected via bevel gears to a planetary gearbox, with simulations indicating optimal performance at 0.5 <inline-formula><mml:math id="M12" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> spacing. Small-scale wind tunnel experiments on DRWTs are more common and have produced varying results.  Habash et al. (2011) and Mitulet et al. (2015) both reported approximately 60 % increases in energy production compared to single-rotor configurations, testing turbines with rotor diameters of 23 cm and 2.5 m, respectively, at spacings ranging from 0.3 to 2.3  and 0.4 <inline-formula><mml:math id="M13" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula>. The first used an individual generator per rotor, and the latter used a double-rotational armature configuration.  Zhao et al. (2020) found relative improvements in the CP of 5.3 %–28.9 % with a 0.55 m radius turbine at optimal 0.3 <inline-formula><mml:math id="M14" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> spacing and individual generators per rotor. The CP values found were 0.34–0.41 versus 0.31–0.35 for single-rotor operation at 8–14 m s<sup>−1</sup> wind speed. A wind tunnel test including particle image velocimetry (PIV) measurements comparing co- and counter-rotating DRWT systems conducted by Ozbay et al. (2014) found a 60 % power increase for a counter-rotating DRWT and a 48 % increase for a co-rotating DRWT compared to a single-rotor configuration. Individual turbine models with 1.27 m radius at different spacings (0.7–6.5 <inline-formula><mml:math id="M16" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula>) were tested by Yuan et al. (2014). It was found that counter-rotating produced 20 % more power compared to co-rotating in the near wake, which decreased to only 4 % at 5.0 <inline-formula><mml:math id="M17" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula>. Wang et al. (2018) tested asymmetric rotor designs (0.28 m front, 0.15 m rear, with 0.25 <inline-formula><mml:math id="M18" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> spacing) and observed smaller gains of 7.2 % for counter-rotating and 1.8 % for co-rotating configurations. At 3.5 <inline-formula><mml:math id="M19" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula>, Mühle et al. (2017) found a significant improvement of 2 % in power production by operating the upstream rotor counter rotating.</p>
      <p id="d2e320">On the contrary, no power increase was found for a bevel gear CR-DRWT system with 1.6 m rotors at 0.64 <inline-formula><mml:math id="M20" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> spacing by Schepers et al. (2024). However, more recently, Adema et al. (2025) tested the same turbine with a double-rotational armature design instead of bevel gears and observed a 10 % CP increase, highlighting how mechanical and electrical design may influence the observed performance of DRWT systems.</p>
      <p id="d2e330">Modelling of DRWT systems can be categorized in two primary approaches. Blade element momentum (BEM) and computational fluid dynamics (CFD) offer complementary insights into DRWT performance with the BEM method commonly used as the preferred tool (Peng et al., 2025). Several BEM models have been developed for DRWT systems. Lee et al. (2012) proposed a modified BEM model for CR-DRWTs in which the rear rotor operates in the fully developed wake of the front rotor, prescribing the axial velocity deficit and neglecting tangential effects and downstream-to-upstream interactions. Yin et al. (2022) improved physical fidelity by coupling a BE/M–T model with free-wake lifting-line data, capturing both axial and tangential inter-rotor velocity interferences, albeit with increased modelling complexity. Amoretti et al. (2023) instead developed a configurable BEM-based approach that retains the dominant upstream wake effects while deliberately neglecting rear-to-front feedback. This offers a robust and computationally efficient framework that is particularly well suited for CR-DRWT performance prediction and parametric design studies and that will be used in this study. Experimental validation is crucial to increase the accuracy of these methods. Both Amoretti et al. (2023) and  Yin et al. (2022) developed BEM models dedicated to DRWT systems, with the first finding an increase in CP of 10.6 % for a CR-DRWT with respect to a single rotor and with the latter finding a 5 % CP increase. Additionally, several studies modelled DRWT systems using CFD simulations. Wang et al. (2022) modelled a diverse range of DRWT systems using CFD-RANS (co- and counter-rotating rotors of varied sizes, including configurations with equally but also non-equally sized rotors). Compared to single-rotor systems, an increase in performance is found for all configurations. The benefit in terms of performance in a counter-rotating system compared to in a co-rotating system is less conclusive. A net benefit of 7 % CP is found by Rosenberg et al. (2014) using RANS by modelling the NREL 5 MW rotor in combination with a 25 % size secondary rotor at 0.2 <inline-formula><mml:math id="M21" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> and a 4.6 % benefit using large-eddy simulations. Another CFD simulation on the NREL 5 MW rotor with an additional rotor of 180 kW at 0.1 <inline-formula><mml:math id="M22" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> downstream of the main rotor resulted in additional performance of 1.74 % despite both rotors having slight efficiency decreases (Peng et al., 2025). For equal-sized rotors at 0.5 <inline-formula><mml:math id="M23" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> spacing, a peak CP of 0.53 is found by Koehuan et al. (2017) using CFD simulations.</p>
      <p id="d2e354">The modelling of DRWT systems has furthermore received growing attention as a means to interpret experiments and to explore operational parameters. Building on the work of Amoretti et al. (2023), the present study adopts and extends their model to represent a counter-rotating DRWT equipped with a double-rotational armature generator, using experimental data as the primary reference for validation and model refinement. The core contribution of this study lies in the presentation and analysis of experimental wind tunnel measurements obtained from a counter-rotating dual-rotor wind turbine tested at the Centre Scientifique et Technique du Bâtiment (CSTB) in Nantes, France, within the framework of the 2025 International Small Wind Turbine Competition. The tested turbine represents an improved design iteration of the system presented by Adema et al. (2025) and incorporates additional instrumentation, including independent rotational speed and pitch measurements for both rotors. These measurements enable a more detailed characterization of rotor interactions, operating conditions, and power performance that has, so far, been scarcely documented for CR-DRWTs with double-rotational armature design. In addition to reporting experimental results, the study employs an extended BEM model to interpret the measurements and to perform an operation optimization of the pitch angle settings and the design tip speed ratios. By grounding the modelling and optimization directly in experimental evidence, this work provides new reference data and valuable insights that support the future development and optimization of compact counter-rotating wind turbines for small-scale applications.</p>
      <p id="d2e357">The remainder of the paper is structured as follows: Sect. 2.1 presents the design of the CR-DRWT, after which the wind tunnel set-up and procedure are explained in Sect. 2.2. The BEM model, as well as adaptations for a double-rotational armature design, is shown in Sect. 2.3. An optimization algorithm for operational parameters is laid out in Sect. 2.4. The results of the wind tunnel test, BEM model, and optimization are shown in Sect. 3, after which the results are discussed in Sect. 4. Conclusions are presented in Sect. 5.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Methods and materials</title>
      <p id="d2e368">This chapter describes the method and materials used in the study. The mechanical, electrical, and aerodynamic designs of the CR-DRWT are presented, after which the test procedure in the wind tunnel is explained. The BEM model used is explained, including proposed extensions to model a dual-rotational armature design. Finally, an optimization strategy is presented to determine optimal operating conditions for future designs.</p>
<sec id="Ch1.S2.SS1">
  <label>2.1</label><title>Design of the CR-DRWT</title>
      <p id="d2e378">The mechanical, electrical, and aerodynamic designs of the CR-DRWT used in this study are presented in earlier work for a different wind tunnel (Adema et al., 2025). The mechanical and electrical designs of the current turbine are identical to those in the earlier work. The generator is still an inline Windstream 1.4 kW permanent magnet direct current (PMDC) with the following operating conditions: 0–120 V, maximum current of 20 A for 30 min, 10 A continuous operation, and startup torque of 0.153 Nm. The overall efficiency of the generator under the contest conditions has been determined in previous work to be 70 %–90 %  (Adema et al., 2025; Jansma, 2024). An overview of the mechanical design, including the blade pitching mechanism and blade mounts, is shown in Fig. 1. Some significant changes have been made to the aerodynamic design, as well as to the sensors on the turbine.</p>

      <fig id="F1" specific-use="star"><label>Figure 1</label><caption><p id="d2e383">Overview of the mechanical lay-out of the CR-DRWT. <bold>(a)</bold> Nacelle overview, <bold>(b)</bold> pitch mechanism (in which A is the linear actuators, B is the generator, C is the generator cap, D is the disc brakes, E is the slip ring, F is the plate with external bearing, G shows the plates with internal bearing, H is the push–pull plate, I is the triangular place, and J shows the blade holders.). Taken from Adema et al. (2025).</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f01.png"/>

        </fig>

      <p id="d2e398">Numerous studies emphasize the importance of measuring blade pitch angles and rotor RPM during wind tunnel measurements of wind turbines (Adema et al., 2025; Amoretti et al., 2023; Bai et al., 2023; Bontempo and Manna, 2025; Erturk et al., 2018; Zhao et al., 2020). Thereto, additional sensors measuring RPM and pitch angles are mounted on the CR-DRWT. Rotor RPM is measured with light pulses passing through slotted discs connected to the rotor and stator and measured by two calibrated light gate sensors. During testing, only the upstream rotor RPM measurement was functional. Pitch angles are set with four Actuonix Motion P16-R linear actuators controlled by an Arduino Mega board with two L298B motor drivers. An Arduino-based control system is used to operate pitch and brake actuators, while the generator output is monitored in real time through a Python-based data acquisition interface.</p>
      <p id="d2e402">The initial pitch angle for the mounting system for the blades is changed slightly, and, therefore, the twist distribution at zero pitch changes. The corresponding chord and twist distribution are shown in Fig. 2 alongside a zoom-in of the root section of the blade and connection to the blade holders. The blades consist of SG6043 airfoils and are fabricated with a 3D-printed core laminated with two carbon fibre layers. For simplicity, the downstream rotor is mirrored with respect to the upstream rotor. The distance between rotors is fixed at 0.64 m (0.39 <inline-formula><mml:math id="M24" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula>) for this CR-DRWT. The rotor diameter for both rotors is 1.6 m, leading to a 2 m<sup>2</sup> swept frontal area.</p>

      <fig id="F2" specific-use="star"><label>Figure 2</label><caption><p id="d2e423">Aerodynamic design of both rotors of the DRWT. Left: the chord distribution. Middle: the twist distribution. Right: zoom-in on the root section of the blades.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f02.png"/>

        </fig>

</sec>
<sec id="Ch1.S2.SS2">
  <label>2.2</label><title>Wind tunnel set-up and test procedure</title>
      <p id="d2e440">The CR-DRWT is tested in the aerodynamic test section at CSTB. The section is 12 m long and 6 m wide and has a height of 5 m. The airflow is controlled to keep a uniform and constant velocity, taking into consideration air density variations. The maximum free-stream wind speed is 70 m s<sup>−1</sup>. Turbulence intensity in the empty test section is less than 1.5 % (Braud et al., 2024). The blockage ratio of the 2 m<sup>2</sup> rotor is 6.7 %, which is below the 10 % mentioned in the literature, above which corrections need to be made (Al-Obaidi and Madivaanan, 2022; Chen and Liou, 2011; Jeong et al., 2018).</p>
      <p id="d2e464">The CR-DRWT is placed on a steel frame connected to the floor of the wind tunnel; see Fig. 3. For mounting, testing, and dismantling of the turbines, a time slot of 2 h was available. Wind speeds were increased in 1 m s<sup>−1</sup> increments from the moment the turbine self-started, starting at 4 m s<sup>−1</sup> up to a maximum of 15 m s<sup>−1</sup>. The DC generator is connected to a programmable Chroma DC load. The resistance of the programmable load is increased at each wind speed until the maximum power point is found. The resistance, voltage, current, and power production are measured for 10 s with 0.1 s intervals. From the 100 data points per wind speed, average values and standard deviations are determined. The measured average power is then combined with the tunnel speed, air density, and rotor surface area to determine the power coefficient (CP) of the wind turbine.</p>

      <fig id="F3" specific-use="star"><label>Figure 3</label><caption><p id="d2e505">The CR-DRWT in the wind tunnel at CSTB.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f03.jpg"/>

        </fig>

</sec>
<sec id="Ch1.S2.SS3">
  <label>2.3</label><title>BEM Model for dual-rotational armature</title>
      <p id="d2e522">The BEM model in this study is based on the work of Amoretti et al. (2023), who developed a configurable BEM model for dual-rotor systems with two distinct rotors with individual generators. The present work extends this model: as in a dual-rotational armature configuration, both rotors are connected through a single generator. The model has the following assumptions: the flow is considered to be inviscid and incompressible, the system is in a stationary state, the retroaction from the second rotor to the first one is neglected, there is no radial speed taken into account, and the expansion of the flow at the rear of the first rotor is not considered. The last assumption is made because the expansion of two adjacent concentric annular sections would result in the intersection of the two respective flows and give complex flow velocity calculations (Amoretti et al., 2023). The model and the proposed extensions to it are presented in the following sections.</p>
<sec id="Ch1.S2.SS3.SSS1">
  <label>2.3.1</label><title>Single-rotor BEM</title>
      <p id="d2e533">For this study, the definition of rotor, blade, and blade elements is in accordance with what is shown in Fig. 4. A rotor with radius <inline-formula><mml:math id="M31" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi mathvariant="normal">rotor</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and hub radius <inline-formula><mml:math id="M32" display="inline"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi mathvariant="normal">hub</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> rotates at angular velocity <inline-formula><mml:math id="M33" display="inline"><mml:mi mathvariant="normal">Ω</mml:mi></mml:math></inline-formula> about the axial direction <inline-formula><mml:math id="M34" display="inline"><mml:mover accent="true"><mml:mrow><mml:msub><mml:mi>e</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow><mml:mo mathvariant="normal">→</mml:mo></mml:mover></mml:math></inline-formula>. Each blade is discretized into <inline-formula><mml:math id="M35" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula> elements of length <inline-formula><mml:math id="M36" display="inline"><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:mi>r</mml:mi></mml:mrow></mml:math></inline-formula> at radial position <inline-formula><mml:math id="M37" display="inline"><mml:mi>r</mml:mi></mml:math></inline-formula>, characterized by chord length <inline-formula><mml:math id="M38" display="inline"><mml:mi>c</mml:mi></mml:math></inline-formula>, twist angle <inline-formula><mml:math id="M39" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula>, and the aerodynamic profile.</p>

      <fig id="F4" specific-use="star"><label>Figure 4</label><caption><p id="d2e620">BEM parameter definition. Left: global wind turbine; middle: blade element definition; right: flow velocities and forces seen by a blade element. Adapted from Figs. 2 to 4 from Amoretti et al. (2023).</p></caption>
            <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f04.png"/>

          </fig>

      <p id="d2e629">The axial and tangential forces on a blade element are then as in Eqs. (1) and  (2):

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M40" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E1"><mml:mtd><mml:mtext>1</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi>F</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo>⋅</mml:mo><mml:mi>c</mml:mi><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:mi>r</mml:mi><mml:mo>⋅</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mi mathvariant="normal">br</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>⋅</mml:mo><mml:mo mathsize="1.5em">[</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">L</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>cos⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">D</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mo mathsize="1.5em">]</mml:mo><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E2"><mml:mtd><mml:mtext>2</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi>F</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo>⋅</mml:mo><mml:mi>c</mml:mi><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:mi>r</mml:mi><mml:mo>⋅</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mi mathvariant="normal">br</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>⋅</mml:mo><mml:mo mathsize="1.5em">[</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">L</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">D</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>cos⁡</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mo mathsize="1.5em">]</mml:mo><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            where <inline-formula><mml:math id="M41" display="inline"><mml:mi mathvariant="italic">ρ</mml:mi></mml:math></inline-formula> is the air density; <inline-formula><mml:math id="M42" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">br</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the relative velocity seen by the blade; <inline-formula><mml:math id="M43" display="inline"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">L</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M44" display="inline"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">D</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are the airfoil lift and drag coefficients; and <inline-formula><mml:math id="M45" display="inline"><mml:mi mathvariant="normal">Φ</mml:mi></mml:math></inline-formula> is the sum of the angle of attack <inline-formula><mml:math id="M46" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula>, pitch angle, and twist angle <inline-formula><mml:math id="M47" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula>. The torque produced by a blade element is described as in Eq. (3), leading to a total rotor power as in Eq. (4), in which <inline-formula><mml:math id="M48" display="inline"><mml:mi>B</mml:mi></mml:math></inline-formula> is the number of blades of the turbine. The power coefficient (<inline-formula><mml:math id="M49" display="inline"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">p</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) can be determined as in Eq. (5), where <inline-formula><mml:math id="M50" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi mathvariant="normal">ax</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the upstream wind speed. 

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M51" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E3"><mml:mtd><mml:mtext>3</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mi mathvariant="normal">d</mml:mi><mml:mi>M</mml:mi><mml:mo>=</mml:mo><mml:mi>r</mml:mi><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:msub><mml:mi>F</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E4"><mml:mtd><mml:mtext>4</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">turbine</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>B</mml:mi><mml:mo>⋅</mml:mo><mml:msubsup><mml:mo>∫</mml:mo><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi mathvariant="normal">rotor</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi mathvariant="normal">hub</mml:mi></mml:msub></mml:mrow></mml:msubsup><mml:mi mathvariant="normal">Ω</mml:mi><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">d</mml:mi><mml:mi>M</mml:mi></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E5"><mml:mtd><mml:mtext>5</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">p</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">turbine</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">wind</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">turbine</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mfenced close="]" open="["><mml:mrow><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo>⋅</mml:mo><mml:mi mathvariant="italic">π</mml:mi><mml:mo>⋅</mml:mo><mml:msup><mml:msub><mml:mi>R</mml:mi><mml:mi mathvariant="normal">rotor</mml:mi></mml:msub><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>⋅</mml:mo><mml:msup><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi mathvariant="normal">ax</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:mrow></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            The axial and tangential velocities seen by the blade are modified by the induction coefficients <inline-formula><mml:math id="M52" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:msup><mml:mi>a</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>, leading to the relative velocity following Eq. (8).

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M54" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E6"><mml:mtd><mml:mtext>6</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">bax</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mo>-</mml:mo><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mi>a</mml:mi></mml:mrow></mml:mfenced><mml:mo>⋅</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi mathvariant="normal">ax</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E7"><mml:mtd><mml:mtext>7</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">bt</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mo>+</mml:mo><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msup><mml:mi>a</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:mfenced><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">Ω</mml:mi><mml:mo>⋅</mml:mo><mml:mi>r</mml:mi></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E8"><mml:mtd><mml:mtext>8</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">br</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:msubsup><mml:mi>V</mml:mi><mml:mi mathvariant="normal">bax</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mi mathvariant="normal">bt</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msubsup></mml:mrow></mml:msqrt></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            The induction coefficients are computed iteratively via a fixed-point algorithm presented in Amoretti et al. (2023), and, for each iteration, the coefficients are calculated using Eqs. (9) and (10). The iterative convergence criterion is set to <inline-formula><mml:math id="M55" display="inline"><mml:mrow><mml:mi mathvariant="italic">ε</mml:mi><mml:mo>=</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>.

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M56" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E9"><mml:mtd><mml:mtext>9</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mn mathvariant="normal">4</mml:mn><mml:mo>⋅</mml:mo><mml:mi>F</mml:mi><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:msup><mml:mfenced close=")" open="("><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mi>s</mml:mi><mml:mo>⋅</mml:mo><mml:mfenced close="]" open="["><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">L</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>cos⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mo>+</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">D</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E10"><mml:mtd><mml:mtext>10</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi>a</mml:mi><mml:msub><mml:msup><mml:mi/><mml:mo>′</mml:mo></mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:mstyle displaystyle="false"><mml:mfrac style="text"><mml:mrow><mml:mn mathvariant="normal">4</mml:mn><mml:mo>⋅</mml:mo><mml:mi>F</mml:mi><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mo>⋅</mml:mo><mml:mi>cos⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced></mml:mrow><mml:mrow><mml:mi>s</mml:mi><mml:mo>⋅</mml:mo><mml:mfenced open="[" close="]"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">L</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mo>-</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">D</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>cos⁡</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            In the above, <inline-formula><mml:math id="M57" display="inline"><mml:mi>s</mml:mi></mml:math></inline-formula> is the local solidity shown in Eq. (11), and <inline-formula><mml:math id="M58" display="inline"><mml:mi>F</mml:mi></mml:math></inline-formula> is Prandtl's tip loss factor (Hansen, 2015), as in Eqs. (12) and (13):

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M59" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E11"><mml:mtd><mml:mtext>11</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi>s</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mi>c</mml:mi><mml:mo>⋅</mml:mo><mml:mi>B</mml:mi></mml:mrow><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mo>⋅</mml:mo><mml:mi mathvariant="italic">π</mml:mi><mml:mo>⋅</mml:mo><mml:mi>r</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E12"><mml:mtd><mml:mtext>12</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mi mathvariant="normal">F</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">2</mml:mn><mml:mi mathvariant="italic">π</mml:mi></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">arccos</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mi>f</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E13"><mml:mtd><mml:mtext>13</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi>f</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>B</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mi mathvariant="normal">rotor</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mi>r</mml:mi><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            For high induction factors (<inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:mi>a</mml:mi><mml:mo>&gt;</mml:mo><mml:mn mathvariant="normal">0.2</mml:mn></mml:mrow></mml:math></inline-formula>), Spera's correction is applied (Hansen, 2015), and Eq. (9) is replaced by the following:

              <disp-formula id="Ch1.E14" content-type="numbered"><label>14</label><mml:math id="M61" display="block"><mml:mtable rowspacing="0.2ex" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>)</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mn mathvariant="normal">2</mml:mn></mml:mfrac></mml:mstyle><mml:mo>⋅</mml:mo><mml:mo mathsize="2.0em">[</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mo>+</mml:mo><mml:mi>K</mml:mi><mml:mo>⋅</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mi>a</mml:mi></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>-</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:mi>K</mml:mi><mml:mfenced close=")" open="("><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mi>a</mml:mi></mml:mrow></mml:mfenced><mml:mo>+</mml:mo><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:mn mathvariant="normal">4</mml:mn><mml:mfenced open="(" close=")"><mml:mrow><mml:mi>K</mml:mi><mml:mo>⋅</mml:mo><mml:msup><mml:mi>a</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:mfenced></mml:mrow></mml:msqrt><mml:mo mathsize="2.0em">]</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>

            with

              <disp-formula id="Ch1.E15" content-type="numbered"><label>15</label><mml:math id="M62" display="block"><mml:mrow><mml:mi>K</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mn mathvariant="normal">4</mml:mn><mml:mo>⋅</mml:mo><mml:mi>F</mml:mi><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:msup><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:mi>s</mml:mi><mml:mo>⋅</mml:mo><mml:mfenced open="[" close="]"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">L</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>cos⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced><mml:mo>+</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">D</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:mi>sin⁡</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="normal">Φ</mml:mi></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
</sec>
<sec id="Ch1.S2.SS3.SSS2">
  <label>2.3.2</label><title>Wake velocity evolution</title>
      <p id="d2e1711">For the second rotor, the axial and tangential velocities downstream of the upstream rotor evolve with the distance after the rotor. The velocity evolution is defined with a distance coefficient <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">distance</mml:mi></mml:msub><mml:mo>(</mml:mo><mml:mi>x</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> based on propeller stages  (Gur, 2019), in which  <inline-formula><mml:math id="M64" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> is the distance from the first rotor.

              <disp-formula id="Ch1.E16" content-type="numbered"><label>16</label><mml:math id="M65" display="block"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">distance</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mi>x</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>x</mml:mi><mml:msqrt><mml:mrow><mml:msup><mml:mi>x</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:msub><mml:mi>R</mml:mi><mml:mi mathvariant="normal">rotor</mml:mi></mml:msub><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:msqrt></mml:mfrac></mml:mstyle></mml:mrow></mml:math></disp-formula>

            The wake velocities at distance  <inline-formula><mml:math id="M66" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> from the rotor are then calculated as follows:

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M67" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E17"><mml:mtd><mml:mtext>17</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">ax</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mi>x</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi mathvariant="normal">ax</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:mrow></mml:msub><mml:mo>⋅</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">distance</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mi>x</mml:mi></mml:mfenced><mml:mo>⋅</mml:mo><mml:mi>a</mml:mi><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E18"><mml:mtd><mml:mtext>18</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">t</mml:mi></mml:msub><mml:mfenced open="(" close=")"><mml:mi>x</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:msup><mml:mi>a</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>⋅</mml:mo><mml:mi mathvariant="normal">Ω</mml:mi><mml:mo>⋅</mml:mo><mml:mi>r</mml:mi><mml:mo>⋅</mml:mo><mml:msub><mml:mi>C</mml:mi><mml:mi mathvariant="normal">distance</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mi>x</mml:mi></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula></p>
</sec>
<sec id="Ch1.S2.SS3.SSS3">
  <label>2.3.3</label><title>Dual-rotor configuration</title>
      <p id="d2e1895">For the second rotor at distance <inline-formula><mml:math id="M68" display="inline"><mml:mrow><mml:mi>x</mml:mi><mml:mo>=</mml:mo><mml:mi>d</mml:mi></mml:mrow></mml:math></inline-formula> downstream from the upstream rotor, the upstream wind velocities are now <inline-formula><mml:math id="M69" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">ax</mml:mi></mml:msub><mml:mfenced open="(" close=")"><mml:mi>d</mml:mi></mml:mfenced></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">t</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mi>d</mml:mi></mml:mfenced></mml:mrow></mml:math></inline-formula>. The axial velocity seen by the second rotor blade is calculated using Eq. (19). The tangential velocity depends on rotation direction – in this case, the counter-rotating configuration, thus following Eq. (20).

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M71" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E19"><mml:mtd><mml:mtext>19</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi mathvariant="normal">bax</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:mfenced><mml:mo>⋅</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">ax</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mi>d</mml:mi></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E20"><mml:mtd><mml:mtext>20</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi mathvariant="normal">bt</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>+</mml:mo><mml:mi>a</mml:mi><mml:msub><mml:msup><mml:mi/><mml:mo>′</mml:mo></mml:msup><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:mfenced><mml:mo>⋅</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub><mml:mo>⋅</mml:mo><mml:mi>r</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mi mathvariant="normal">t</mml:mi></mml:msub><mml:mfenced open="(" close=")"><mml:mi>d</mml:mi></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            The BEM algorithm is applied to both rotors sequentially, with the downwind rotor calculation incorporating modified inflow conditions. All other BEM equations (forces, induction factors, corrections) remain structurally identical for the second rotor. The model accepts geometrical parameters (blade geometry, rotor radii, number of blades), aerodynamic data (lift–drag polars), operational parameters (rotational speeds), environmental conditions (wind speed, air density), and configuration settings (rotor spacing, rotation direction).</p>
</sec>
<sec id="Ch1.S2.SS3.SSS4">
  <label>2.3.4</label><title>Extension to the model for dual-rotational armature design</title>
      <p id="d2e2045">In order for a CR-DRWT with a dual-rotational armature design to be in balance, both rotor torques need to be equal (Kutt et al., 2020; Li et al., 2021). For a range of prescribed tip speed ratios of the upstream rotor, valid operating tip speed ratios of the downstream rotor are determined by enforcing a torque balance condition under the wake induced by the upstream rotor. The procedure used to determine these torque-balanced operating points is described in Table 1.</p>

<table-wrap id="T1" specific-use="star"><label>Table 1</label><caption><p id="d2e2051">Algorithm to determine balance torque of CR-DRWT.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="1">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:thead>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">Algorithm torque balance CR-DRWT:</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1"><italic>for each</italic> wind speed:</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><italic>for each</italic> prescribed front rotor tip speed ratio:</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">  Compute the aerodynamic torque and rotational speed of the upstream rotor</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">  Compute the wake conditions at the downstream rotor plane induced by the upstream rotor</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">  Fix the wake conditions and evaluate the downstream rotor torque over a predefined range of rear rotor tip speed ratios</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">  Identify all intervals where the downstream rotor torque crosses the upstream rotor torque</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><italic>for each</italic> identified sign change interval:</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">     Apply bisection to find the downstream rotor tip speed ratio at which both torques are balanced</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">  Compute the downstream rotor rotational speed and power output</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">  Compute the total system power</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">  Store the torque balanced operating point</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><italic>end for</italic></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1"><italic>end for</italic></oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1"><italic>end for</italic></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Output: torque-balanced operating points and power values for all wind speeds</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

      <p id="d2e2186">The output is a list of operating points where the system is in balance according to Eq. (22). Throughout the algorithm shown in Table 1, a tolerance for torque balance of 0.05 Nm is given in order to find possible solutions. For each point of balance, the net generator RPM is the sum of the RPM of both rotors as the system is counter-rotating, as shown in Eq. (21). The power then follows Eq. (23), and the power coefficient of the CR-DRWT is according to Eq. (24), where <inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">wind</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is calculated as in Eq. (5).

                  <disp-formula specific-use="gather" content-type="numbered"><mml:math id="M73" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E21"><mml:mtd><mml:mtext>21</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mi mathvariant="normal">generator</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mtext>upstream rotor</mml:mtext></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mtext>downstream rotor</mml:mtext></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E22"><mml:mtd><mml:mtext>22</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>T</mml:mi><mml:mtext>CR-DRWT</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mtext>upstream rotor</mml:mtext></mml:msub><mml:mspace width="0.125em" linebreak="nobreak"/><mml:mo>=</mml:mo><mml:mspace linebreak="nobreak" width="0.125em"/><mml:mspace linebreak="nobreak" width="0.125em"/><mml:msub><mml:mi>T</mml:mi><mml:mtext>downstream rotor</mml:mtext></mml:msub></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E23"><mml:mtd><mml:mtext>23</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mi>P</mml:mi><mml:mtext>CR-DRWT</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mi mathvariant="normal">generator</mml:mi></mml:msub><mml:mo>⋅</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mtext>CR-DRWT</mml:mtext></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E24"><mml:mtd><mml:mtext>24</mml:mtext></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi mathvariant="normal">p</mml:mi><mml:mi mathvariant="normal">_</mml:mi><mml:mi mathvariant="normal">total</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mtext>CR-DRWT</mml:mtext></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi mathvariant="normal">wind</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula></p>
</sec>
</sec>
<sec id="Ch1.S2.SS4">
  <label>2.4</label><title>Optimization of extended CR-DRWT BEM model</title>
      <p id="d2e2331">The ISWTC competition of 2025 included a Weibull distribution (<inline-formula><mml:math id="M74" display="inline"><mml:mrow><mml:mi>A</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">7.1</mml:mn></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M75" display="inline"><mml:mrow><mml:mi>k</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">2.4</mml:mn></mml:mrow></mml:math></inline-formula>) to calculate the annual energy production (AEP). Thereto, the turbine is optimized for performance under these conditions. A sweep optimization strategy to identify the pitch angle combination that maximizes the annual energy production (AEP) of the dual-rotor wind turbine is employed. The upstream and downstream rotor pitch angles are each varied independently across a grid from 0 to 10° in steps of 1°, yielding 121 unique pitch combinations to be evaluated in full. For each combination, the aerodynamic performance is computed across a range of wind speeds from 4 to 15 m s<sup>−1</sup>, with a 1 m s<sup>−1</sup> interval, using the extended BEM model. For each wind speed, only the operating point yielding the highest overall power output is selected. The AEP for each pitch combination is then calculated using the Weibull parameters. The pitch combination producing the highest total AEP across all wind speeds is then identified as the global optimum. Finally, the BEM will then be simulated with the theoretically optimal pitch and tip speed ratio (TSR) settings to assess the maximum performance of the CR-DRWT. This will then be compared to the wind tunnel measurements as well. </p>
</sec>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Results</title>
      <p id="d2e2392">This chapter present the results of the wind tunnel test and BEM model, including optimization. First, the results of the wind tunnel are presented together with the BEM model validation.</p>
<sec id="Ch1.S3.SS1">
  <label>3.1</label><title>Airfoil data</title>
      <p id="d2e2402">The aerodynamic data to describe the airfoil performances were generated with the 2D airfoil solver XFoil in QBlade  (QBlade, 2026; Drela, 1989). Lift and drag coefficients for angles of attack from <inline-formula><mml:math id="M78" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>10  to 25° are calculated with a transition amplification ratio of <inline-formula><mml:math id="M79" display="inline"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="normal">crit</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">9</mml:mn></mml:mrow></mml:math></inline-formula>, corresponding to clean wind tunnel conditions. Figure 5 displays the aerodynamic data. The curve is then extrapolated to 360° (Montgomerie, 2004). For the BEM code, the aerodynamic data are interpolated according to the individual Reynolds number for each blade element.</p>

      <fig id="F5"><label>Figure 5</label><caption><p id="d2e2429">Profile aerodynamics of SG6043 airfoil at <inline-formula><mml:math id="M80" display="inline"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi mathvariant="normal">crit</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">9</mml:mn></mml:mrow></mml:math></inline-formula>.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f05.png"/>

        </fig>

</sec>
<sec id="Ch1.S3.SS2">
  <label>3.2</label><title>Wind tunnel test</title>
      <p id="d2e2461">Air pressure in the wind tunnel at the time of the experiment was 1018 hPa, average humidity was 51.2 %, and average temperature was 25.5 °C. The density is then calculated to be 1.18 kg m<sup>−3</sup>. The wind tunnel speed is determined through the measurement of the dynamic pressure and air density. Rotor RPM was measured in 10 RPM increments. The data were acquired through a laptop interface connected with the turbine, and measurements were sampled at discrete wind speed intervals rather than being continuously logged due to the streaming nature of the sensor output. The generator RPM is determined through the voltage and current measurements, as well as the generator characteristics. At 6 m s<sup>−1</sup>, optimal operational pitch settings were determined by manually changing the pitch angles and resistance values of the DC load to reach maximum power production. The settings for the collective pitch angles were 8.6° for the upstream rotor and 1.7° for the downstream rotor throughout the experiment. The results for the wind tunnel test are presented in Table 2. The turbine self-starts around 3.5 m s<sup>−1</sup> and starts producing power at 4 m s<sup>−1</sup>. A startup procedure is developed where the downstream rotor is allowed to operate first by breaking the upstream rotor and consecutively allowing the upstream rotor to start only once the downwind rotor reached an operating point. Power, voltage, and current output are measured for 10 s with a 0.1 s sampling rate. Mean values are calculated, as well as the standard deviation. Maximum power production was 1013.79 <inline-formula><mml:math id="M85" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 8.58 W at 15.00 m s<sup>−1</sup> wind tunnel speed. The maximum RPM of the generator was 2150 RPM. The rear rotor was generally rotating slower than the front rotor at lower wind speeds, while, at higher wind speeds, the RPM became more similar for both rotors, and the rear rotor seemed to stabilize. Net generator RPM did keep increasing with wind speed.</p>

<table-wrap id="T2" specific-use="star"><label>Table 2</label><caption><p id="d2e2535">Measurement results from CSTB wind tunnel tests.</p></caption><oasis:table frame="topbot"><oasis:tgroup cols="8">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="right"/>
     <oasis:colspec colnum="3" colname="col3" align="right"/>
     <oasis:colspec colnum="4" colname="col4" align="right"/>
     <oasis:colspec colnum="5" colname="col5" align="right"/>
     <oasis:colspec colnum="6" colname="col6" align="right"/>
     <oasis:colspec colnum="7" colname="col7" align="right"/>
     <oasis:colspec colnum="8" colname="col8" align="right"/>
     <oasis:thead>
       <oasis:row>
         <oasis:entry colname="col1">Wind tunnel</oasis:entry>
         <oasis:entry colname="col2">Voltage</oasis:entry>
         <oasis:entry colname="col3">Current</oasis:entry>
         <oasis:entry colname="col4">Power</oasis:entry>
         <oasis:entry colname="col5">CP</oasis:entry>
         <oasis:entry colname="col6">Downstream</oasis:entry>
         <oasis:entry colname="col7">Upstream</oasis:entry>
         <oasis:entry colname="col8">Generator</oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">speed [m s<sup>−1</sup>]</oasis:entry>
         <oasis:entry colname="col2">[V]</oasis:entry>
         <oasis:entry colname="col3">[A]</oasis:entry>
         <oasis:entry colname="col4">[W]</oasis:entry>
         <oasis:entry colname="col5">[–]</oasis:entry>
         <oasis:entry colname="col6">rotor [RPM]</oasis:entry>
         <oasis:entry colname="col7">rotor [RPM]</oasis:entry>
         <oasis:entry colname="col8">[RPM]</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">4.14</oasis:entry>
         <oasis:entry colname="col2">22.71 <inline-formula><mml:math id="M88" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.54</oasis:entry>
         <oasis:entry colname="col3">0.42 <inline-formula><mml:math id="M89" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.01</oasis:entry>
         <oasis:entry colname="col4">9.49 <inline-formula><mml:math id="M90" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 1.63</oasis:entry>
         <oasis:entry colname="col5">0.11 <inline-formula><mml:math id="M91" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.02</oasis:entry>
         <oasis:entry colname="col6">60</oasis:entry>
         <oasis:entry colname="col7">340</oasis:entry>
         <oasis:entry colname="col8">400</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">5.27</oasis:entry>
         <oasis:entry colname="col2">31.57 <inline-formula><mml:math id="M92" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.33</oasis:entry>
         <oasis:entry colname="col3">1.05 <inline-formula><mml:math id="M93" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.01</oasis:entry>
         <oasis:entry colname="col4">33.31 <inline-formula><mml:math id="M94" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 4.00</oasis:entry>
         <oasis:entry colname="col5">0.19 <inline-formula><mml:math id="M95" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.02</oasis:entry>
         <oasis:entry colname="col6">260</oasis:entry>
         <oasis:entry colname="col7">390</oasis:entry>
         <oasis:entry colname="col8">650</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">5.90</oasis:entry>
         <oasis:entry colname="col2">39.38 <inline-formula><mml:math id="M96" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.29</oasis:entry>
         <oasis:entry colname="col3">1.43 <inline-formula><mml:math id="M97" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.03</oasis:entry>
         <oasis:entry colname="col4">56.33 <inline-formula><mml:math id="M98" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 1.61</oasis:entry>
         <oasis:entry colname="col5">0.23 <inline-formula><mml:math id="M99" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.01</oasis:entry>
         <oasis:entry colname="col6">250</oasis:entry>
         <oasis:entry colname="col7">540</oasis:entry>
         <oasis:entry colname="col8">790</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">7.10</oasis:entry>
         <oasis:entry colname="col2">38.38 <inline-formula><mml:math id="M100" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.05</oasis:entry>
         <oasis:entry colname="col3">3.48 <inline-formula><mml:math id="M101" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.01</oasis:entry>
         <oasis:entry colname="col4">133.58 <inline-formula><mml:math id="M102" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 1.27</oasis:entry>
         <oasis:entry colname="col5">0.32 <inline-formula><mml:math id="M103" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">210</oasis:entry>
         <oasis:entry colname="col7">630</oasis:entry>
         <oasis:entry colname="col8">840</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">8.16</oasis:entry>
         <oasis:entry colname="col2">41.75 <inline-formula><mml:math id="M104" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.22</oasis:entry>
         <oasis:entry colname="col3">5.01 <inline-formula><mml:math id="M105" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.02</oasis:entry>
         <oasis:entry colname="col4">209.25 <inline-formula><mml:math id="M106" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 1.99</oasis:entry>
         <oasis:entry colname="col5">0.33 <inline-formula><mml:math id="M107" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">245</oasis:entry>
         <oasis:entry colname="col7">705</oasis:entry>
         <oasis:entry colname="col8">950</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">9.14</oasis:entry>
         <oasis:entry colname="col2">49.86 <inline-formula><mml:math id="M108" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.22</oasis:entry>
         <oasis:entry colname="col3">6.04 <inline-formula><mml:math id="M109" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.01</oasis:entry>
         <oasis:entry colname="col4">301.05 <inline-formula><mml:math id="M110" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 2.20</oasis:entry>
         <oasis:entry colname="col5">0.33 <inline-formula><mml:math id="M111" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">360</oasis:entry>
         <oasis:entry colname="col7">830</oasis:entry>
         <oasis:entry colname="col8">1190</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">10.05</oasis:entry>
         <oasis:entry colname="col2">57.54 <inline-formula><mml:math id="M112" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.19</oasis:entry>
         <oasis:entry colname="col3">6.97 <inline-formula><mml:math id="M113" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.01</oasis:entry>
         <oasis:entry colname="col4">401.31 <inline-formula><mml:math id="M114" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 1.50</oasis:entry>
         <oasis:entry colname="col5">0.33 <inline-formula><mml:math id="M115" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">400</oasis:entry>
         <oasis:entry colname="col7">940</oasis:entry>
         <oasis:entry colname="col8">1340</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">11.15</oasis:entry>
         <oasis:entry colname="col2">66.45 <inline-formula><mml:math id="M116" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.15</oasis:entry>
         <oasis:entry colname="col3">8.07 <inline-formula><mml:math id="M117" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.05</oasis:entry>
         <oasis:entry colname="col4">535.91 <inline-formula><mml:math id="M118" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 3.40</oasis:entry>
         <oasis:entry colname="col5">0.33 <inline-formula><mml:math id="M119" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">520</oasis:entry>
         <oasis:entry colname="col7">1000</oasis:entry>
         <oasis:entry colname="col8">1520</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">12.18</oasis:entry>
         <oasis:entry colname="col2">73.51 <inline-formula><mml:math id="M120" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.23</oasis:entry>
         <oasis:entry colname="col3">8.91 <inline-formula><mml:math id="M121" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.07</oasis:entry>
         <oasis:entry colname="col4">654.75 <inline-formula><mml:math id="M122" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 5.59</oasis:entry>
         <oasis:entry colname="col5">0.31 <inline-formula><mml:math id="M123" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">620</oasis:entry>
         <oasis:entry colname="col7">1050</oasis:entry>
         <oasis:entry colname="col8">1670</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">13.04</oasis:entry>
         <oasis:entry colname="col2">80.12 <inline-formula><mml:math id="M124" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.22</oasis:entry>
         <oasis:entry colname="col3">9.74 <inline-formula><mml:math id="M125" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.07</oasis:entry>
         <oasis:entry colname="col4">780.29 <inline-formula><mml:math id="M126" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 5.54</oasis:entry>
         <oasis:entry colname="col5">0.30 <inline-formula><mml:math id="M127" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">800</oasis:entry>
         <oasis:entry colname="col7">1050</oasis:entry>
         <oasis:entry colname="col8">1850</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">13.93</oasis:entry>
         <oasis:entry colname="col2">85.91 <inline-formula><mml:math id="M128" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.32</oasis:entry>
         <oasis:entry colname="col3">10.37 <inline-formula><mml:math id="M129" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.08</oasis:entry>
         <oasis:entry colname="col4">890.98 <inline-formula><mml:math id="M130" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 5.54</oasis:entry>
         <oasis:entry colname="col5">0.28 <inline-formula><mml:math id="M131" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">960</oasis:entry>
         <oasis:entry colname="col7">1050</oasis:entry>
         <oasis:entry colname="col8">2010</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">15.00</oasis:entry>
         <oasis:entry colname="col2">91.72 <inline-formula><mml:math id="M132" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.31</oasis:entry>
         <oasis:entry colname="col3">11.05 <inline-formula><mml:math id="M133" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.07</oasis:entry>
         <oasis:entry colname="col4">1013.79 <inline-formula><mml:math id="M134" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 8.58</oasis:entry>
         <oasis:entry colname="col5">0.25 <inline-formula><mml:math id="M135" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula> 0.00</oasis:entry>
         <oasis:entry colname="col6">1100</oasis:entry>
         <oasis:entry colname="col7">1050</oasis:entry>
         <oasis:entry colname="col8">2150</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup></oasis:table></table-wrap>

      <p id="d2e3308">A simulation using the BEM code from Sect. 2.3, with identical inputs for pitch during the experiment, has been performed. The corresponding pair of rotational speeds where torque balance is present and which are closest in power output to the measurement is selected. The results for the power production, the CP, and RPMs are shown in Figs. 6,  7, and  8, together with the CSTB wind tunnel results. The BEM code imposes a 70 %–90 % electrical efficiency of the generator on the aerodynamic power. This efficiency of the generator has been determined for similar operational conditions in previous studies  (Adema et al., 2025; Jansma, 2024). The CP values for the downstream rotor are calculated with the reduced wind speed after the upstream rotor following Eq. (17).</p>

      <fig id="F6" specific-use="star"><label>Figure 6</label><caption><p id="d2e3314">Results of power production during the wind tunnel test compared to BEM results. The figure shows both rotors individually, the CR-DRWT, and the wind tunnel test data.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f06.png"/>

        </fig>

      <fig id="F7" specific-use="star"><label>Figure 7</label><caption><p id="d2e3325">Results of the coefficient of performance during the wind tunnel test compared to BEM results. The figure shows both rotors individually, the CR-DRWT, and the wind tunnel test data.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f07.png"/>

        </fig>

      <fig id="F8" specific-use="star"><label>Figure 8</label><caption><p id="d2e3336">Results of rotational speeds during the wind tunnel test compared to BEM results. The figure shows both rotors individually, the CR-DRWT, and the wind tunnel test data.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f08.png"/>

        </fig>

      <p id="d2e3345">Overall, the BEM model follows the outcome of the wind tunnel tests well with respect to power output and CP values. Both rotors operate at similar CP values. However, the modelled RPM values start to deviate from around 11 m s<sup>−1</sup> onwards. The measurements show a stable upstream rotor RPM of 1050, but the model requires a faster rotating upstream rotor to enforce torque balance at that power point. The model, however, does follow the operational RPM of the generator, indicating that the inferred RPM values from the generator characteristic were correct. The discussion in the next section further elaborates on the observed deviation in RPM measurements.</p>
</sec>
<sec id="Ch1.S3.SS3">
  <label>3.3</label><title>Optimization</title>
      <p id="d2e3368">The results of the optimization sweep are presented in Fig. 9. The plot shows the cumulative results of 121 pairs of pitch combinations. The optimization favours when the upstream rotor operates at a high pitch angle and when the downstream rotor operates at a low pitch angle. The best configuration is an upstream rotor at pitch 9° and a downstream rotor at 0°, yielding a possible AEP of 2150 kWh. The results of the extended BEM model with the optimized pitch settings are shown in Figs. 10 and 11. The maximum aerodynamic power output of the CR-DRWT at 15 m s<sup>−1</sup> is 2251 W, with a maximum CP of 0.56 at 14 and 15 m s<sup>−1</sup>. This maximum aerodynamic power is above the rated capacity of 1400 W of the current generator. The corresponding TSR values are calculated using the free-stream wind speed for the upstream rotor and the reduced axial velocity for the downstream rotor. The mean TSR values across the wind speed range for the upstream and downstream rotor are 3.2 and 6.6, respectively, leading to a maximum RPM of the generator at 15 m s<sup>−1</sup> of 1603 RPM, which is significantly below the measured 2150 RPM in the wind tunnel test at the same wind speed. Compared to the wind tunnel tests at CSTB, the optimized CR-DRWT could outperform the current turbine. The current plots show only aerodynamic performance. However, even including 70 %–90 % electrical losses, the optimized turbine outperforms the current measurements. Also, Fig. 11 shows that the optimization maximizes the performance of the downstream rotor, with a much higher CP for the downstream rotor for all wind speeds compared to the CP for the upstream rotor.</p>

      <fig id="F9" specific-use="star"><label>Figure 9</label><caption><p id="d2e3409">Results of optimization.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f09.png"/>

        </fig>

      <fig id="F10" specific-use="star"><label>Figure 10</label><caption><p id="d2e3420">Power production of optimized CR-DRWT.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f10.png"/>

        </fig>

      <fig id="F11" specific-use="star"><label>Figure 11</label><caption><p id="d2e3432">Power coefficient of optimized CR-DRWT.</p></caption>
          <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f11.png"/>

        </fig>

</sec>
</sec>
<sec id="Ch1.S4">
  <label>4</label><title>Discussion</title>
      <p id="d2e3450">The current iteration of the CR-DRWT shows improvement with respect to earlier work (Adema et al., 2025), especially as precise RPM and pitch measurements are performed. Absolute CP values of around 0.5 are found in both field studies and modelling of CR-DRWT systems (Amoretti et al., 2023; Jung et al., 2005; Koehuan et al., 2017). The current experimental results reach a maximum CP of only 0.33. The findings of this study are more in line with earlier wind tunnel experiments on a CR-DRWT of comparable size (Mitulet et al., 2015; Zhao et al., 2020). There is therefore still potential to increase the efficiency of turbines with limited rotor size. The optimization does show theoretical values reaching a CP of 0.56, confirming the possible improvements.</p>
      <p id="d2e3453">While the extended BEM model follows the experiment measurements well, the deviation between the modelled and measured upstream rotor RPM values above 11 m s<sup>−1</sup> warrants further consideration. The experimental measurements at the upstream rotor plateau at approximately 1050 RPM, and the model predicts a continuously increasing RPM to keep torque balance. A possible explanation lies in the measurement system employed during the wind tunnel campaign. The upstream rotor RPM was measured using a light gate sensor reading optical pulses through a slotted disc. It is suspected that a limitation in the pulse-counting implementation code effectively imposed an artificial ceiling on the recorded RPM values. This would cause the sensor output to remain at a fixed reading even as the true rotational speed continued to increase. The generator RPM, derived from the generator characteristic curve, is therefore considered to be the more reliable indicator of actual generator RPM. This interpretation is further supported by the fact that the total generator RPM coming from the BEM model matches well with the derived generator RPM.</p>
      <p id="d2e3468">The theoretical optimization shows the downstream rotor operating with high CP values and limiting the upstream rotor efficiency. This raises questions about if a single rotor could achieve better AEP performance than the optimized CR-DRWT. Thereto, the BEM model is run with an isolated upstream rotor, effectively modelling single-rotor operation, and the optimal operation point is found by sweeping across a range of pitch angles and TSR values and maximizing AEP. The optimal aerodynamic performance of a single rotor is found at 0° pitch and a TSR of 6.5. The results are shown in Fig. 10. A single rotor can reach equally good power and efficiency. The CR-DRWT has slightly higher AEP with 2322 kWh compared to 2093 kWh for the single rotor. However, the operational point of both systems is different, with the CR-DRWT reaching maximum power at 1603 RPM at the generator and 1164 RPM for the single rotor, highlighting that, for both configurations, a different generator is likely needed to reach optimal operation. What is worth mentioning is that Adema et al. (2025) concluded that the current direct-drive configuration is not optimized for single-rotor operation. The same turbine with 2.5° pitch for the upstream rotor and a stationary downstream rotor only reached a maximum CP in single-rotor operation of 0.23. For the current small wind turbine generator, the additional costs for a second rotor may be lower than those for the addition of a gearbox, highlighting the potential of DRWT systems in compact simplified designs (Booker et al., 2010).</p>

      <fig id="F12" specific-use="star"><label>Figure 12</label><caption><p id="d2e3474">Comparison between theoretical optimal power production and power coefficient of CR-DRWT and single rotor.</p></caption>
        <graphic xlink:href="https://wes.copernicus.org/articles/11/2037/2026/wes-11-2037-2026-f12.png"/>

      </fig>

      <p id="d2e3483">In most existing wind tunnel experiments, Reynolds numbers are much lower than those encountered in actual situations, limiting their ability to replicate real-world conditions. This limits the ability to model and predict the performance of DRWT's under real-world conditions (Hollands et al., 2020). The same holds for the current experiment with Reynolds numbers between 25 000 and 200 000. Attention needs to be paid to possible inaccuracies in predicting airfoil characteristics at such low numbers.</p>
      <p id="d2e3486">An attempt to use a BEM model incorporating both axial and tangential influences on the inflow at the downstream rotor (Amoretti et al., 2023), as well as an extension for a dual-rotational armature design, has shown good agreement with the current findings, though some deviations remain. An additional CFD analysis of the current turbine configuration presented in this study will be a valuable addition in understanding the flow field around both rotors. Knowing the detailed flow field around the downstream rotor, a detailed aerodynamic design can be performed to further optimize the CR-DRWT. Also, such an analysis may reveal unknown (three-dimensional) aerodynamic effects not captured in current BEM model.</p>
      <p id="d2e3489">Finally, calculation methods have been proposed for optimizing the downstream rotor geometry for DRWT designs (Bontempo and Manna, 2025; Wang et al., 2018). Currently, for simplicity, the downstream rotor is mirrored with respect to the upstream rotor. Using the BEM model presented in this work, a detailed blade design using these design methods can be performed to further increase the performance of dual-rotor wind turbines.</p>
</sec>
<sec id="Ch1.S5" sec-type="conclusions">
  <label>5</label><title>Conclusions</title>
      <p id="d2e3501">This study investigates the performance of a CR-DRWT with a dual-rotational armature configuration through wind tunnel testing at CSTB in Nantes, France. The research extends previous work by incorporating enhanced measurement capabilities including RPM and pitch angle sensors, the development of an extended BEM model for dual-rotational armature systems, and implementation of a sweep optimization strategy to identify optimal operating parameters. Wind tunnel tests were conducted at wind speeds ranging from 4 to 15 m s<sup>−1</sup>, with a 1.6 m diameter rotor and two rotors with 0.39 <inline-formula><mml:math id="M142" display="inline"><mml:mi>D</mml:mi></mml:math></inline-formula> spacing. The CR-DRWT achieved a maximum power output of 1014 W. The extended BEM model is validated against experimental data. The main conclusions are outlined below.</p>
      <p id="d2e3523">The extended BEM mode for dual-rotational armature design showed good overall agreement with wind tunnel measurements for power production, efficiency, and generator RPM, and an optimization strategy revealed optimum operational parameters. The upstream rotor RPM measured by the light gate sensor was limited at approximately 1050 RPM above 11 m s<sup>−1</sup>, likely due to a pulse-counting limitation in the data acquisition code, as the derived generator RPM continued to rise, in line with model predictions across the full wind speed range.</p>
      <p id="d2e3538">Experimental wind tunnel testing achieved a maximum CP of 0.33, which is higher than previous iterations of the same turbine but below the theoretically optimized value of 0.56, indicating substantial room for improvement through optimal pitch angle, tip speed ratios for both rotors, and generator choice or redesign.</p>
      <p id="d2e3541">The sweep optimization identified an upstream rotor pitch of 9° and downstream rotor pitch of 0° as the optimal configuration for this CR-DRWT, with a theoretical AEP of 2150 kWh. The optimization consistently favoured a high-efficiency downstream rotor operation at high CP values, while the upstream rotor operates with reduced efficiency, indicating that shared power production across both rotors is essential for maximizing overall system performance.</p>
      <p id="d2e3545">The CR-DRWT with dual-rotational armature design demonstrated reliable operation with self-starting capability at 3.5 m s<sup>−1</sup>. The dual-rotational armature configuration eliminates the need for gearboxes, potentially offering a more cost-effective and mechanically simpler solution for compact small wind turbine designs, particularly suited for urban applications requiring low noise and vibration.</p>
      <p id="d2e3560">To advance the understanding and performance of CR-DRWT systems, detailed computational fluid dynamics simulations should be performed to capture 3D flow effects not modelled by the current BEM approach. Experiments at higher Reynolds numbers closer to real-world operational conditions are needed to improve performance prediction accuracy and to validate airfoil characteristics. The development of blade geometries specifically optimized for the downstream rotor, rather than using mirrored configurations, should be pursued using the validated BEM model and advanced design methods from the recent literature. Field testing under real atmospheric conditions is essential to assess the practical performance, durability, and economic viability of the optimized CR-DRWT configuration.</p>
</sec>

      
      </body>
    <back><notes notes-type="codedataavailability"><title>Code and data availability</title>

      <p id="d2e3567">The models and data used in this paper can be requested from the corresponding author upon request.</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d2e3573">NA: conceptualization, formal analysis, writing (original draft preparation). JGV: investigation, resources, writing (original draft preparation). WS: investigation, resources, writing (review and editing). AdR: funding acquisition, supervision, writing (review and editing). GS: writing (review and editing).</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d2e3579">The contact author has declared that none of the authors has any competing interests.</p>
  </notes><notes notes-type="disclaimer"><title>Disclaimer</title>

      <p id="d2e3585">Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.</p>
  </notes><ack><title>Acknowledgements</title><p id="d2e3591">The authors of this paper would like to express their sincere gratitude to Callum Morgan, Sjoerd Kluvers, Wessel Smeenge, Bram Roorda, and Marijn Verbeek, who helped design and build the current iteration of the turbine used in this work; their efforts cannot go unmentioned. Second, the staff at CSTB is thanked for their assistance and guidance during the wind tunnel measurements. Finally, the authors would like to thank Caroline Braud for her effort in facilitating the hosting of the ISWTC at CSTB in 2025.</p></ack><notes notes-type="financialsupport"><title>Financial support</title>

      <p id="d2e3596">This research has been supported by the Nationaal Regieorgaan Praktijkgericht Onderzoek SIA (grant no. RAAK.MKB17.007).</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d2e3602">This paper was edited by Jonathan Whale and reviewed by two anonymous referees.</p>
  </notes><ref-list>
    <title>References</title>

      <ref id="bib1.bib1"><label>1</label><mixed-citation>Adema, N., Swart-Ranshuysen, W., De Ruijter, A., and Schepers, G.: Wind Tunnel Test of Counter-Rotating Dual Rotor Wind Turbine With Double Rotational Armature Design, Wind Energy, 28, <ext-link xlink:href="https://doi.org/10.1002/we.70039" ext-link-type="DOI">10.1002/we.70039</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bib2"><label>2</label><mixed-citation>Al-Obaidi, A. S. M. and Madivaanan, G.: Investigation of the Blockage Correction to Improve the Accuracy of Taylor's Low-Speed Wind Tunnel, J. Phys. Conf. Ser., 2222, 012008, <ext-link xlink:href="https://doi.org/10.1088/1742-6596/2222/1/012008" ext-link-type="DOI">10.1088/1742-6596/2222/1/012008</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib3"><label>3</label><mixed-citation>Amoretti, T., Huet, F., Garambois, P., and Roucoules, L.: Configurable dual rotor wind turbine model based on BEM method: Co-rotating and counter-rotating comparison, Energ. Convers. Manage., 293, 117461, <ext-link xlink:href="https://doi.org/10.1016/j.enconman.2023.117461" ext-link-type="DOI">10.1016/j.enconman.2023.117461</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib4"><label>4</label><mixed-citation>Bai, H., Wang, N., and Wan, D.: Numerical study of aerodynamic performance of horizontal axis dual-rotor wind turbine under atmospheric boundary layers, Ocean Eng., 280, 114944, <ext-link xlink:href="https://doi.org/10.1016/j.oceaneng.2023.114944" ext-link-type="DOI">10.1016/j.oceaneng.2023.114944</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib5"><label>5</label><mixed-citation>Bereziartua-Gonzalez, L., Retegi, A., and Ukar, O.: Human-centered integration of small wind turbines in urban environments: a semi-systematic review from an industrial design perspective, Front. Sustain. Cities, 7, <ext-link xlink:href="https://doi.org/10.3389/frsc.2025.1561894" ext-link-type="DOI">10.3389/frsc.2025.1561894</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bib6"><label>6</label><mixed-citation>Bianchini, A., Bangga, G., Baring-Gould, I., Croce, A., Cruz, J. I., Damiani, R., Erfort, G., Simao Ferreira, C., Infield, D., Nayeri, C. N., Pechlivanoglou, G., Runacres, M., Schepers, G., Summerville, B., Wood, D., and Orrell, A.: Current status and grand challenges for small wind turbine technology, Wind Energ. Sci., 7, 2003–2037, <ext-link xlink:href="https://doi.org/10.5194/wes-7-2003-2022" ext-link-type="DOI">10.5194/wes-7-2003-2022</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib7"><label>7</label><mixed-citation>Bontempo, R. and Manna, M.: Optimum design of contra-rotating wind turbines with adjacent rotors, Energ. Conv. Manage., 324, 119267, <ext-link xlink:href="https://doi.org/10.1016/j.enconman.2024.119267" ext-link-type="DOI">10.1016/j.enconman.2024.119267</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bib8"><label>8</label><mixed-citation>Booker, J. D., Mellor, P. H., Wrobel, R., and Drury, D.: A compact, high efficiency contra-rotating generator suitable for wind turbines in the urban environment, Renew. Energ., 35, 2027–2033, <ext-link xlink:href="https://doi.org/10.1016/j.renene.2010.02.003" ext-link-type="DOI">10.1016/j.renene.2010.02.003</ext-link>, 2010.</mixed-citation></ref>
      <ref id="bib1.bib9"><label>9</label><mixed-citation>Braud, C., Podvin, B., and Deparday, J.: Study of the wall pressure variations on the stall inception of a thick cambered profile at high Reynolds number, Phys. Rev. Fluids, 9, 014605, <ext-link xlink:href="https://doi.org/10.1103/PhysRevFluids.9.014605" ext-link-type="DOI">10.1103/PhysRevFluids.9.014605</ext-link>, 2024.</mixed-citation></ref>
      <ref id="bib1.bib10"><label>10</label><mixed-citation>Calautit, K. and Johnstone, C.: State-of-the-art review of micro to small-scale wind energy harvesting technologies for building integration, Energ. Conv. Manage.,  20, 100457, <ext-link xlink:href="https://doi.org/10.1016/j.ecmx.2023.100457" ext-link-type="DOI">10.1016/j.ecmx.2023.100457</ext-link>, 2023.</mixed-citation></ref>
      <ref id="bib1.bib11"><label>11</label><mixed-citation>Chagas, C. C. M., Pereira, M. G., Rosa, L. P., da Silva, N. F., Freitas, M. A. V., and Hunt, J. D.: From Megawatts to Kilowatts: A Review of Small Wind Turbine Applications, Lessons From The US to Brazil, Sustainability, 12, 2760, <ext-link xlink:href="https://doi.org/10.3390/su12072760" ext-link-type="DOI">10.3390/su12072760</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib12"><label>12</label><mixed-citation>Chen, T. Y. and Liou, L. R.: Blockage corrections in wind tunnel tests of small horizontal-axis wind turbines, Exp. Therm. Fluid Sci., 35, 565–569, <ext-link xlink:href="https://doi.org/10.1016/j.expthermflusci.2010.12.005" ext-link-type="DOI">10.1016/j.expthermflusci.2010.12.005</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib13"><label>13</label><mixed-citation>Drela, M.: XFOIL: An Analysis and Design System for Low Reynolds Number Airfoils, in: Low Reynolds Number Aerodynamics, 1–12, <ext-link xlink:href="https://doi.org/10.1007/978-3-642-84010-4_1" ext-link-type="DOI">10.1007/978-3-642-84010-4_1</ext-link>, 1989.</mixed-citation></ref>
      <ref id="bib1.bib14"><label>14</label><mixed-citation>Erturk, E., Sivrioglu, S., and Bolat, F. C.: Analysis Model of a Small Scale Counter-Rotating Dual Rotor Wind Turbine with Double Rotational Generator Armature, International Journal of Renewable Energy Research (IJRER), 8, 1849–1858, 2018.</mixed-citation></ref>
      <ref id="bib1.bib15"><label>15</label><mixed-citation>Gur, O.: Extending Blade-Element Model to Contra-Rotating Configuration, IOP Conf. Ser., Mater. Sci. Eng., 638, 012001, <ext-link xlink:href="https://doi.org/10.1088/1757-899X/638/1/012001" ext-link-type="DOI">10.1088/1757-899X/638/1/012001</ext-link>, 2019.</mixed-citation></ref>
      <ref id="bib1.bib16"><label>16</label><mixed-citation>Habash, R. W. Y., Groza, V., Yang, Y., Blouin, C., and Guillemette, P.: Performance of a Contrarotating Small Wind Energy Converter, International Scholarly Research Notices, 2011, 828739, <ext-link xlink:href="https://doi.org/10.5402/2011/828739" ext-link-type="DOI">10.5402/2011/828739</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bib17"><label>17</label><mixed-citation>Hansen, M.: Aerodynamics of Wind Turbines, 3rd Edn., Routledge, London, 188 pp., <ext-link xlink:href="https://doi.org/10.4324/9781315769981" ext-link-type="DOI">10.4324/9781315769981</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib18"><label>18</label><mixed-citation>Hollands, E. O., He, C., and Gan, L.: A particle image velocimetry study of dual-rotor counter-rotating wind turbine near wake, J. Vis., 23, 425–435, <ext-link xlink:href="https://doi.org/10.1007/s12650-020-00643-0" ext-link-type="DOI">10.1007/s12650-020-00643-0</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib19"><label>19</label><mixed-citation>Jansma, H.: Development and Optimization of a Generator Test Rig for Small Wind Turbines: Improving Performance and Efficiency, MSc Thesis, Hanze UAS, Groningen, 74 pp., 2024.</mixed-citation></ref>
      <ref id="bib1.bib20"><label>20</label><mixed-citation>Jeong, H., Lee, S., and Kwon, S.-D.: Blockage corrections for wind tunnel tests conducted on a Darrieus wind turbine, J. Wind Eng. Ind. Aerod., 179, 229–239, <ext-link xlink:href="https://doi.org/10.1016/j.jweia.2018.06.002" ext-link-type="DOI">10.1016/j.jweia.2018.06.002</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib21"><label>21</label><mixed-citation>Jung, S. N., No, T.-S., and Ryu, K.-W.: Aerodynamic performance predictionof a 30 kW counter-rotating wind turbine system, Renew. Energ., 30, 631–644, <ext-link xlink:href="https://doi.org/10.1016/j.renene.2004.07.005" ext-link-type="DOI">10.1016/j.renene.2004.07.005</ext-link>, 2005.</mixed-citation></ref>
      <ref id="bib1.bib22"><label>22</label><mixed-citation>Jurasz, J., Bochenek, B., Wieczorek, J., Jaczewski, A., Kies, A., and Figurski, M.: Energy potential and economic viability of small-scale wind turbines, Energy, 322, 135608, <ext-link xlink:href="https://doi.org/10.1016/j.energy.2025.135608" ext-link-type="DOI">10.1016/j.energy.2025.135608</ext-link>, 2025.</mixed-citation></ref>
      <ref id="bib1.bib23"><label>23</label><mixed-citation>Koehuan, V. A., Sugiyono, and Kamal, S.: Investigation of Counter-Rotating Wind Turbine Performance using Computational Fluid Dynamics Simulation, IOP Conf. Ser. Mater. Sci. Eng., 267, 012034, <ext-link xlink:href="https://doi.org/10.1088/1757-899X/267/1/012034" ext-link-type="DOI">10.1088/1757-899X/267/1/012034</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib24"><label>24</label><mixed-citation>Kutt, F., Blecharz, K., and Karkosiński, D.: Axial-Flux Permanent-Magnet Dual-Rotor Generator for a Counter-Rotating Wind Turbine, Energies, 13, 2833, <ext-link xlink:href="https://doi.org/10.3390/en13112833" ext-link-type="DOI">10.3390/en13112833</ext-link>, 2020.</mixed-citation></ref>
      <ref id="bib1.bib25"><label>25</label><mixed-citation>Lee, S., Kim, H., Son, E., and Lee, S.: Effects of design parameters on aerodynamic performance of a counter-rotating wind turbine, Renew. Energ., 42, 140–144, <ext-link xlink:href="https://doi.org/10.1016/j.renene.2011.08.046" ext-link-type="DOI">10.1016/j.renene.2011.08.046</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bib26"><label>26</label><mixed-citation>Li, Y., Zhang, J., Li, Z., Yang, P., and Wang, H.: Design and verification of a novel double rotor without stator wind turbine generation system, Energy Reports, 7, 161–172, <ext-link xlink:href="https://doi.org/10.1016/j.egyr.2021.10.041" ext-link-type="DOI">10.1016/j.egyr.2021.10.041</ext-link>, 2021.</mixed-citation></ref>
      <ref id="bib1.bib27"><label>27</label><mixed-citation>Mitulet, L.-A., Oprina, G., Chihaia, R.-A., Nicolaie, S., Nedelcu, A., and Popescu, M.: Wind Tunnel Testing for a New Experimental Model of Counter-rotating Wind Turbine, Procedia Engineer., 100, 1141–1149, <ext-link xlink:href="https://doi.org/10.1016/j.proeng.2015.01.477" ext-link-type="DOI">10.1016/j.proeng.2015.01.477</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib28"><label>28</label><mixed-citation>Montgomerie, B.: Methods for Root Effects, Tip Effects and Extending the Angle of Attack Range to <inline-formula><mml:math id="M145" display="inline"><mml:mo>±</mml:mo></mml:math></inline-formula>180<inline-formula><mml:math id="M146" display="inline"><mml:mi mathvariant="italic">°</mml:mi></mml:math></inline-formula>, with Application to Aerodynamics for Blades on Wind Turbines and Propellers, Swedish Defence Research Agency, <uri>https://www.foi.se/rest-api/report/FOI-R--1305--SE</uri> (last access: 3 June 2026), 2004.</mixed-citation></ref>
      <ref id="bib1.bib29"><label>29</label><mixed-citation>Mühle, F., Adaramola, M. S., and Sætran, L.: The effect of rotational direction on the wake of a wind turbine rotor – a comparison study of aligned co- and counter rotating turbine arrays, Energy Proced., 137, 238–245, <ext-link xlink:href="https://doi.org/10.1016/j.egypro.2017.10.346" ext-link-type="DOI">10.1016/j.egypro.2017.10.346</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib30"><label>30</label><mixed-citation>Newman, B. G.: Multiple actuator-disc theory for wind turbines,J. Wind Eng. Ind. Aerod., 24, 215–225, <ext-link xlink:href="https://doi.org/10.1016/0167-6105(86)90023-1" ext-link-type="DOI">10.1016/0167-6105(86)90023-1</ext-link>, 1986.</mixed-citation></ref>
      <ref id="bib1.bib31"><label>31</label><mixed-citation>Ozbay, A., Tian, W., and Hu, H.: An Experimental Investigation on the Aeromechanics and Near Wake Characteristics of Dual-Rotor Wind Turbines (DRWTs), in: 32nd ASME Wind Energy Symposium, American Institute of Aeronautics and Astronautics, <ext-link xlink:href="https://doi.org/10.2514/6.2014-1085" ext-link-type="DOI">10.2514/6.2014-1085</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bib32"><label>32</label><mixed-citation>Peng, X., Duan, L., Li, G., Jin, Y., and Han, Z.: Interference between main and auxiliary rotors in floating dual-rotor wind turbines under stationary and surge conditions, Ocean Eng., 322, 120462, <ext-link xlink:href="https://doi.org/10.1016/j.oceaneng.2025.120462" ext-link-type="DOI">10.1016/j.oceaneng.2025.120462</ext-link>, 2025. </mixed-citation></ref>
      <ref id="bib1.bib33"><label>33</label><mixed-citation>QBlade: QBlade – Next Generation Wind Turbine Simulation, <uri>https://qblade.org/</uri>, last access: 13 April 2026.</mixed-citation></ref>
      <ref id="bib1.bib34"><label>34</label><mixed-citation>Rosenberg, A., Selvaraj, S., and Sharma, A.: A Novel Dual-Rotor Turbine for Increased Wind Energy Capture, J. Phys. Conf. Ser., 524, 012078, <ext-link xlink:href="https://doi.org/10.1088/1742-6596/524/1/012078" ext-link-type="DOI">10.1088/1742-6596/524/1/012078</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bib35"><label>35</label><mixed-citation>Schepers, J. G., Adema, N. C., Lipian, M., Kulak, M., Shahid, A., Best, A., Bendre, T., Gallicchio, I., Elsabbagh, A., Mostafa, A., Kim, T., Mikkelsen, R., Gaunaa, M., Teuwen, J. J. E., Rudolf, R. T., Wood, D., and Holierhoek, J. G.: Lessons learned from 10 years of wind tunnel tests on small wind turbines designed by students, J. Phys. Conf. Ser., 2767, 072009, <ext-link xlink:href="https://doi.org/10.1088/1742-6596/2767/7/072009" ext-link-type="DOI">10.1088/1742-6596/2767/7/072009</ext-link>, 2024.</mixed-citation></ref>
      <ref id="bib1.bib36"><label>36</label><mixed-citation>Simic, Z., Havelka, J. G., and Bozicevic Vrhovcak, M.: Small wind turbines – A unique segment of the wind power market, Renew. Energ., 50, 1027–1036, 2013.</mixed-citation></ref>
      <ref id="bib1.bib37"><label>37</label><mixed-citation>Sundararaju, H., Lo, K. H., Metcalfe, R., and Wang, S. S.: Aerodynamics and CFD analysis of equal size dual-rotor wind turbine, J. Renew. Sustain. Ener., 9, 043305, <ext-link xlink:href="https://doi.org/10.1063/1.4999500" ext-link-type="DOI">10.1063/1.4999500</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib38"><label>38</label><mixed-citation>Wang, K., Liu, T., Wan, Y., Ong, M. C., and Wu, T.: Numerical Investigation on Aerodynamic Characteristics of Dual-Rotor Wind Turbines, J. Mar. Sci. Eng., 10, 1887, <ext-link xlink:href="https://doi.org/10.3390/jmse10121887" ext-link-type="DOI">10.3390/jmse10121887</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib39"><label>39</label><mixed-citation>Wang, Z., Ozbay, A., Tian, W., and Hu, H.: An experimental study on the aerodynamic performances and wake characteristics of an innovative dual-rotor wind turbine, Energy, 147, 94–109, <ext-link xlink:href="https://doi.org/10.1016/j.energy.2018.01.020" ext-link-type="DOI">10.1016/j.energy.2018.01.020</ext-link>, 2018.</mixed-citation></ref>
      <ref id="bib1.bib40"><label>40</label><mixed-citation>Yin, F. F., Chen, J. J., Li, X. K., Ye, Z. L., Tang, W., Shen, X., and Guo, X. J.: A blade element momentum model for dual-rotor wind turbines considering inter-rotor velocity interferences, J. Phys. Conf. Ser., 2265, 042058, <ext-link xlink:href="https://doi.org/10.1088/1742-6596/2265/4/042058" ext-link-type="DOI">10.1088/1742-6596/2265/4/042058</ext-link>, 2022.</mixed-citation></ref>
      <ref id="bib1.bib41"><label>41</label><mixed-citation>Yuan, W., Tian, W., Ozbay, A., and Hu, H.: An experimental study on the effects of relative rotation direction on the wake interferences among tandem wind turbines, Sci. China Phys. Mech. Astron., 57, 935–949, <ext-link xlink:href="https://doi.org/10.1007/s11433-014-5429-x" ext-link-type="DOI">10.1007/s11433-014-5429-x</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bib42"><label>42</label><mixed-citation>Zhao, X., Zhou, P., Liang, X., and Gao, S.: The aerodynamic coupling design and wind tunnel test of contra-rotating wind turbines, Renew. Energ., 146, 1–8, <ext-link xlink:href="https://doi.org/10.1016/j.renene.2019.06.118" ext-link-type="DOI">10.1016/j.renene.2019.06.118</ext-link>, 2020.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>Experimental validation and extension of a blade element momentum model for counter-rotation dual-rotor wind turbine with double rotational armature design</article-title-html>
<abstract-html/>
<ref-html id="bib1.bib1"><label>1</label><mixed-citation>
      Adema, N., Swart-Ranshuysen, W., De Ruijter, A., and Schepers, G.: Wind
Tunnel Test of Counter-Rotating Dual Rotor Wind Turbine With Double
Rotational Armature Design, Wind Energy, 28,
<a href="https://doi.org/10.1002/we.70039" target="_blank">https://doi.org/10.1002/we.70039</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>2</label><mixed-citation>
      Al-Obaidi, A. S. M. and Madivaanan, G.: Investigation of the Blockage
Correction to Improve the Accuracy of Taylor's Low-Speed Wind Tunnel, J.
Phys. Conf. Ser., 2222, 012008,
<a href="https://doi.org/10.1088/1742-6596/2222/1/012008" target="_blank">https://doi.org/10.1088/1742-6596/2222/1/012008</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>3</label><mixed-citation>
      Amoretti, T., Huet, F., Garambois, P., and Roucoules, L.: Configurable dual
rotor wind turbine model based on BEM method: Co-rotating and
counter-rotating comparison, Energ. Convers. Manage., 293, 117461,
<a href="https://doi.org/10.1016/j.enconman.2023.117461" target="_blank">https://doi.org/10.1016/j.enconman.2023.117461</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>4</label><mixed-citation>
      Bai, H., Wang, N., and Wan, D.: Numerical study of aerodynamic performance
of horizontal axis dual-rotor wind turbine under atmospheric boundary
layers, Ocean Eng., 280, 114944,
<a href="https://doi.org/10.1016/j.oceaneng.2023.114944" target="_blank">https://doi.org/10.1016/j.oceaneng.2023.114944</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>5</label><mixed-citation>
      Bereziartua-Gonzalez, L., Retegi, A., and Ukar, O.: Human-centered
integration of small wind turbines in urban environments: a semi-systematic
review from an industrial design perspective, Front. Sustain. Cities, 7,
<a href="https://doi.org/10.3389/frsc.2025.1561894" target="_blank">https://doi.org/10.3389/frsc.2025.1561894</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>6</label><mixed-citation>
      Bianchini, A., Bangga, G., Baring-Gould, I., Croce, A., Cruz, J. I., Damiani, R., Erfort, G., Simao Ferreira, C., Infield, D., Nayeri, C. N., Pechlivanoglou, G., Runacres, M., Schepers, G., Summerville, B., Wood, D., and Orrell, A.: Current status and grand challenges for small wind turbine technology, Wind Energ. Sci., 7, 2003–2037, <a href="https://doi.org/10.5194/wes-7-2003-2022" target="_blank">https://doi.org/10.5194/wes-7-2003-2022</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>7</label><mixed-citation>
      Bontempo, R. and Manna, M.: Optimum design of contra-rotating wind turbines
with adjacent rotors, Energ. Conv. Manage., 324, 119267,
<a href="https://doi.org/10.1016/j.enconman.2024.119267" target="_blank">https://doi.org/10.1016/j.enconman.2024.119267</a>, 2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>8</label><mixed-citation>
      Booker, J. D., Mellor, P. H., Wrobel, R., and Drury, D.: A compact, high
efficiency contra-rotating generator suitable for wind turbines in the urban
environment, Renew. Energ., 35, 2027–2033,
<a href="https://doi.org/10.1016/j.renene.2010.02.003" target="_blank">https://doi.org/10.1016/j.renene.2010.02.003</a>, 2010.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>9</label><mixed-citation>
      Braud, C., Podvin, B., and Deparday, J.: Study of the wall pressure
variations on the stall inception of a thick cambered profile at high
Reynolds number, Phys. Rev. Fluids, 9, 014605,
<a href="https://doi.org/10.1103/PhysRevFluids.9.014605" target="_blank">https://doi.org/10.1103/PhysRevFluids.9.014605</a>, 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>10</label><mixed-citation>
      Calautit, K. and Johnstone, C.: State-of-the-art review of micro to
small-scale wind energy harvesting technologies for building integration,
Energ. Conv. Manage.,  20, 100457,
<a href="https://doi.org/10.1016/j.ecmx.2023.100457" target="_blank">https://doi.org/10.1016/j.ecmx.2023.100457</a>, 2023.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>11</label><mixed-citation>
      Chagas, C. C. M., Pereira, M. G., Rosa, L. P., da Silva, N. F., Freitas, M.
A. V., and Hunt, J. D.: From Megawatts to Kilowatts: A Review of Small Wind
Turbine Applications, Lessons From The US to Brazil, Sustainability, 12,
2760, <a href="https://doi.org/10.3390/su12072760" target="_blank">https://doi.org/10.3390/su12072760</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>12</label><mixed-citation>
      Chen, T. Y. and Liou, L. R.: Blockage corrections in wind tunnel tests of
small horizontal-axis wind turbines, Exp. Therm. Fluid Sci.,
35, 565–569, <a href="https://doi.org/10.1016/j.expthermflusci.2010.12.005" target="_blank">https://doi.org/10.1016/j.expthermflusci.2010.12.005</a>, 2011.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>13</label><mixed-citation>
      Drela, M.: XFOIL: An Analysis and Design System for Low Reynolds Number
Airfoils, in: Low Reynolds Number Aerodynamics, 1–12,
<a href="https://doi.org/10.1007/978-3-642-84010-4_1" target="_blank">https://doi.org/10.1007/978-3-642-84010-4_1</a>, 1989.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>14</label><mixed-citation>
      Erturk, E., Sivrioglu, S., and Bolat, F. C.: Analysis Model of a Small Scale
Counter-Rotating Dual Rotor Wind Turbine with Double Rotational Generator
Armature, International Journal of Renewable Energy Research (IJRER), 8,
1849–1858, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>15</label><mixed-citation>
      Gur, O.: Extending Blade-Element Model to Contra-Rotating Configuration, IOP Conf. Ser., Mater. Sci. Eng., 638, 012001,
<a href="https://doi.org/10.1088/1757-899X/638/1/012001" target="_blank">https://doi.org/10.1088/1757-899X/638/1/012001</a>, 2019.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>16</label><mixed-citation>
      Habash, R. W. Y., Groza, V., Yang, Y., Blouin, C., and Guillemette, P.:
Performance of a Contrarotating Small Wind Energy Converter, International
Scholarly Research Notices, 2011, 828739,
<a href="https://doi.org/10.5402/2011/828739" target="_blank">https://doi.org/10.5402/2011/828739</a>, 2011.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>17</label><mixed-citation>
      Hansen, M.: Aerodynamics of Wind Turbines, 3rd Edn., Routledge, London, 188 pp., <a href="https://doi.org/10.4324/9781315769981" target="_blank">https://doi.org/10.4324/9781315769981</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>18</label><mixed-citation>
      Hollands, E. O., He, C., and Gan, L.: A particle image velocimetry study of
dual-rotor counter-rotating wind turbine near wake, J. Vis., 23, 425–435,
<a href="https://doi.org/10.1007/s12650-020-00643-0" target="_blank">https://doi.org/10.1007/s12650-020-00643-0</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>19</label><mixed-citation>
      Jansma, H.: Development and Optimization of a Generator Test Rig for Small Wind Turbines: Improving Performance and Efficiency, MSc Thesis, Hanze UAS, Groningen,
74 pp., 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>20</label><mixed-citation>
      Jeong, H., Lee, S., and Kwon, S.-D.: Blockage corrections for wind tunnel
tests conducted on a Darrieus wind turbine, J. Wind Eng. Ind. Aerod., 179, 229–239,
<a href="https://doi.org/10.1016/j.jweia.2018.06.002" target="_blank">https://doi.org/10.1016/j.jweia.2018.06.002</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>21</label><mixed-citation>
      Jung, S. N., No, T.-S., and Ryu, K.-W.: Aerodynamic performance predictionof a 30 kW counter-rotating wind turbine system, Renew. Energ., 30,
631–644, <a href="https://doi.org/10.1016/j.renene.2004.07.005" target="_blank">https://doi.org/10.1016/j.renene.2004.07.005</a>, 2005.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>22</label><mixed-citation>
      Jurasz, J., Bochenek, B., Wieczorek, J., Jaczewski, A., Kies, A., and
Figurski, M.: Energy potential and economic viability of small-scale wind
turbines, Energy, 322, 135608, <a href="https://doi.org/10.1016/j.energy.2025.135608" target="_blank">https://doi.org/10.1016/j.energy.2025.135608</a>,
2025.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>23</label><mixed-citation>
      Koehuan, V. A., Sugiyono, and Kamal, S.: Investigation of Counter-Rotating
Wind Turbine Performance using Computational Fluid Dynamics Simulation, IOP
Conf. Ser. Mater. Sci. Eng., 267, 012034,
<a href="https://doi.org/10.1088/1757-899X/267/1/012034" target="_blank">https://doi.org/10.1088/1757-899X/267/1/012034</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>24</label><mixed-citation>
      Kutt, F., Blecharz, K., and Karkosiński, D.: Axial-Flux Permanent-Magnet
Dual-Rotor Generator for a Counter-Rotating Wind Turbine, Energies, 13,
2833, <a href="https://doi.org/10.3390/en13112833" target="_blank">https://doi.org/10.3390/en13112833</a>, 2020.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>25</label><mixed-citation>
      Lee, S., Kim, H., Son, E., and Lee, S.: Effects of design parameters on
aerodynamic performance of a counter-rotating wind turbine, Renew.
Energ., 42, 140–144, <a href="https://doi.org/10.1016/j.renene.2011.08.046" target="_blank">https://doi.org/10.1016/j.renene.2011.08.046</a>, 2012.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>26</label><mixed-citation>
      Li, Y., Zhang, J., Li, Z., Yang, P., and Wang, H.: Design and verification
of a novel double rotor without stator wind turbine generation system,
Energy Reports, 7, 161–172, <a href="https://doi.org/10.1016/j.egyr.2021.10.041" target="_blank">https://doi.org/10.1016/j.egyr.2021.10.041</a>,
2021.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>27</label><mixed-citation>
      Mitulet, L.-A., Oprina, G., Chihaia, R.-A., Nicolaie, S., Nedelcu, A., and
Popescu, M.: Wind Tunnel Testing for a New Experimental Model of
Counter-rotating Wind Turbine, Procedia Engineer., 100, 1141–1149,
<a href="https://doi.org/10.1016/j.proeng.2015.01.477" target="_blank">https://doi.org/10.1016/j.proeng.2015.01.477</a>, 2015.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>28</label><mixed-citation>
      Montgomerie, B.: Methods for Root Effects, Tip Effects and Extending the
Angle of Attack Range to ±180<i>°</i>, with Application to
Aerodynamics for Blades on Wind Turbines and Propellers, Swedish Defence
Research Agency, <a href="https://www.foi.se/rest-api/report/FOI-R-1305-SE" target="_blank"/> (last access: 3 June 2026), 2004.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>29</label><mixed-citation>
      Mühle, F., Adaramola, M. S., and Sætran, L.: The effect of
rotational direction on the wake of a wind turbine rotor – a comparison
study of aligned co- and counter rotating turbine arrays, Energy Proced.,
137, 238–245, <a href="https://doi.org/10.1016/j.egypro.2017.10.346" target="_blank">https://doi.org/10.1016/j.egypro.2017.10.346</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>30</label><mixed-citation>
      Newman, B. G.: Multiple actuator-disc theory for wind turbines,J. Wind Eng. Ind. Aerod., 24, 215–225,
<a href="https://doi.org/10.1016/0167-6105(86)90023-1" target="_blank">https://doi.org/10.1016/0167-6105(86)90023-1</a>, 1986.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib31"><label>31</label><mixed-citation>
      Ozbay, A., Tian, W., and Hu, H.: An Experimental Investigation on the
Aeromechanics and Near Wake Characteristics of Dual-Rotor Wind Turbines
(DRWTs), in: 32nd ASME Wind Energy Symposium, American Institute of
Aeronautics and Astronautics, <a href="https://doi.org/10.2514/6.2014-1085" target="_blank">https://doi.org/10.2514/6.2014-1085</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib32"><label>32</label><mixed-citation>
      Peng, X., Duan, L., Li, G., Jin, Y., and Han, Z.: Interference between main
and auxiliary rotors in floating dual-rotor wind turbines under stationary
and surge conditions, Ocean Eng., 322, 120462,
<a href="https://doi.org/10.1016/j.oceaneng.2025.120462" target="_blank">https://doi.org/10.1016/j.oceaneng.2025.120462</a>, 2025.


    </mixed-citation></ref-html>
<ref-html id="bib1.bib33"><label>33</label><mixed-citation>
      QBlade: QBlade – Next Generation Wind Turbine Simulation, <a href="https://qblade.org/" target="_blank"/>, last access: 13 April 2026.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib34"><label>34</label><mixed-citation>
      Rosenberg, A., Selvaraj, S., and Sharma, A.: A Novel Dual-Rotor Turbine for
Increased Wind Energy Capture, J. Phys. Conf. Ser., 524, 012078,
<a href="https://doi.org/10.1088/1742-6596/524/1/012078" target="_blank">https://doi.org/10.1088/1742-6596/524/1/012078</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib35"><label>35</label><mixed-citation>
      Schepers, J. G., Adema, N. C., Lipian, M., Kulak, M., Shahid, A., Best, A.,
Bendre, T., Gallicchio, I., Elsabbagh, A., Mostafa, A., Kim, T., Mikkelsen,
R., Gaunaa, M., Teuwen, J. J. E., Rudolf, R. T., Wood, D., and Holierhoek,
J. G.: Lessons learned from 10 years of wind tunnel tests on small wind
turbines designed by students, J. Phys. Conf. Ser., 2767, 072009,
<a href="https://doi.org/10.1088/1742-6596/2767/7/072009" target="_blank">https://doi.org/10.1088/1742-6596/2767/7/072009</a>, 2024.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib36"><label>36</label><mixed-citation>
      Simic, Z., Havelka, J. G., and Bozicevic Vrhovcak, M.: Small wind turbines
– A unique segment of the wind power market, Renew. Energ., 50,
1027–1036, 2013.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib37"><label>37</label><mixed-citation>
      Sundararaju, H., Lo, K. H., Metcalfe, R., and Wang, S. S.: Aerodynamics and
CFD analysis of equal size dual-rotor wind turbine, J. Renew.
Sustain. Ener., 9, 043305, <a href="https://doi.org/10.1063/1.4999500" target="_blank">https://doi.org/10.1063/1.4999500</a>, 2017.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib38"><label>38</label><mixed-citation>
      Wang, K., Liu, T., Wan, Y., Ong, M. C., and Wu, T.: Numerical Investigation
on Aerodynamic Characteristics of Dual-Rotor Wind Turbines, J.
Mar. Sci. Eng., 10, 1887,
<a href="https://doi.org/10.3390/jmse10121887" target="_blank">https://doi.org/10.3390/jmse10121887</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib39"><label>39</label><mixed-citation>
      Wang, Z., Ozbay, A., Tian, W., and Hu, H.: An experimental study on the
aerodynamic performances and wake characteristics of an innovative
dual-rotor wind turbine, Energy, 147, 94–109,
<a href="https://doi.org/10.1016/j.energy.2018.01.020" target="_blank">https://doi.org/10.1016/j.energy.2018.01.020</a>, 2018.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib40"><label>40</label><mixed-citation>
      Yin, F. F., Chen, J. J., Li, X. K., Ye, Z. L., Tang, W., Shen, X., and Guo,
X. J.: A blade element momentum model for dual-rotor wind turbines
considering inter-rotor velocity interferences, J. Phys. Conf. Ser., 2265,
042058, <a href="https://doi.org/10.1088/1742-6596/2265/4/042058" target="_blank">https://doi.org/10.1088/1742-6596/2265/4/042058</a>, 2022.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib41"><label>41</label><mixed-citation>
      Yuan, W., Tian, W., Ozbay, A., and Hu, H.: An experimental study on the
effects of relative rotation direction on the wake interferences among
tandem wind turbines, Sci. China Phys. Mech. Astron., 57, 935–949,
<a href="https://doi.org/10.1007/s11433-014-5429-x" target="_blank">https://doi.org/10.1007/s11433-014-5429-x</a>, 2014.

    </mixed-citation></ref-html>
<ref-html id="bib1.bib42"><label>42</label><mixed-citation>
      Zhao, X., Zhou, P., Liang, X., and Gao, S.: The aerodynamic coupling design and wind tunnel test of contra-rotating wind turbines, Renew. Energ.,
146, 1–8, <a href="https://doi.org/10.1016/j.renene.2019.06.118" target="_blank">https://doi.org/10.1016/j.renene.2019.06.118</a>, 2020.

    </mixed-citation></ref-html>--></article>
