3. Mathematical Modelling of Volumetric Coupling and Delivered Power
3.1. Modelling the Spatial Extension of Volumetric Coupling Propagation
Once the physical origin of the force acting on the rotor has been established, the next logical step is to determine the spatial extension of the region perturbed by the influence of an ideal turbine. This characterisation is essential for quantifying the mobilised mass and, consequently, for estimating the intensity of the resulting force.
For the turbine rotor to effectively interact with the flow, the rotor assembly must be connected to an electric generator, which imposes a resistive load on the rotation axis. In the absence of this resistance, the rotor spins freely in synchrony with the flow, allowing particles to pass nearly rectilinearly between the blades with minimal deflection. However, when subjected to the generator's resistive load, the rotor ceases to rotate in synchrony with the stream and begins to oppose the flow. This decoupling causes a deviation from the natural particle trajectories, leading to interactions with the blade profiles and triggering the vectorial deflection of momentum. This process initiates the dynamic coupling between the flow and the rotor.
To describe the propagation of this hydrodynamic perturbation according to the principles of incompressible flow physics, it is considered that its propagation speed corresponds to the speed of pressure waves in the medium, that is, the speed of sound in water, approximately 1500 m/s [
10]. This propagation defines the volumetric coupling length
L, which conceptually represents the spatial extension of upstream mass mobilisation. It physically substantiates the proposed volumetric model as follows:
where Δ
t represents the time interval during which a portion of the flow passes through the rotor under resistive load. This same interval defines the temporal extension of the hydrodynamic disturbance generated, which propagates upstream at speed
c, thereby defining the coupling length.
Although a precise quantitative determination of Δt is not developed at this stage, its conceptual role is essential to justify the coupling between the rotor and the surrounding fluid volume.
The implications of this coupling will be addressed in the following sections. The mathematical development will demonstrate that the interaction time Δt cancels out in the final expressions for both force and power. This cancellation stems from a fundamental property of the continuous regime: although the momentum deflection occurs along a curved trajectory inside the rotor, each point along this path is at all times occupied by a fraction of the coupled volume, whose mass differentially carries the momentum being redirected.
Therefore, the total force on the rotor results from the continuous summation of the forces associated with the vectorial deflection of momentum elements, spatially distributed along the active region. This spatial integral yields the same momentum variation that would be obtained by tracking a single differential of momentum over its entire trajectory. Hence, the individual interaction time becomes irrelevant to the final result. What matters, in this regime, is the continuous rate of momentum reorientation within the coupled volume, sustained by hydrodynamic equilibrium.
3.1.1. Volumetric Attenuation and Proposal of a Realistic Efficiency Factor
As previously discussed, the hydrodynamic reaction generated by the presence of the rotor is not uniformly distributed throughout the mobilised volume. The influence exerted on the pressure field decreases with distance from the perturbation point, resulting in a distributed and attenuated response along the coupling length .
This behaviour is widely recognised in the hydrodynamic literature, which demonstrates that, in incompressible flows, only a fraction of the mobilised mass behaves as if fully coupled to the motion of the solid body. Both classical and contemporary studies, including the theory of added mass and the concept of hydrodynamic inertia, indicate that the fluid response tends to dissipate with distance, so that only part of the volume actually exerts a significant influence on the rotor [
11,
12,
13,
14].
These effects are initially represented by potential flow models, which assume instantaneous and idealised propagation [
15,
16], but are later refined by experimental evidence and numerical simulations that confirm that real coupling is limited and not absolute [
12,
13]. This behaviour is particularly well documented in naval and offshore applications, such as in the study of inertial forces acting on floating and submerged structures
11,
13].
To incorporate this physical refinement into the model, a volumetric efficiency factor, denoted , is introduced. It is defined as the fraction of the mobilised mass that effectively transfers its momentum to the turbine rotor. This factor acknowledges that the idealised model overestimates volumetric coupling and accepts that only part of the mobilised mass actually contributes to the generation of torque and power.
Based on comparative analyses with studies on added mass and on conservative estimates in the literature, a reference value is proposed as:
This value reflects the empirical ratio between the theoretically mobilised mass and the fraction that exerts significant mechanical influence on the solid body. By incorporating this factor into the subsequent equations, the model adjusts its predictions to a more realistic level, consistent with the hydrodynamic behaviour observed in real-world conditions, without compromising the integrity of the theoretical framework and in fact reinforcing its scientific soundness.
It is important to emphasise that, in order to preserve conceptual clarity and maintain a clear distinction from classical theory, the factor will be applied only at the final stage of power estimation, and not in the intermediate expressions for mass, force, or torque. This choice allows the mathematical development to be presented in its idealised form, highlighting the conceptual proposal of volumetric interaction. The attenuation factor will then be incorporated into the final power expression, at which point the cumulative effect of volumetric limitation becomes physically manifest, refining the model’s predictions to reflect a more realistic and technically defensible performance.
3.2. Definition of the Flow Velocity in the Rotor’s Working Annular Region- v″
After characterising the spatial extent of the interaction between the flow and the turbine, the next step is to determine the local flow velocity across the active region of the rotor blades, based on the free-stream velocity . This velocity is essential for estimating the tangential force and, consequently, the torque generated.
In the proposed model, the turbine is equipped with a duct that not only houses the rotor but also channels the incoming current, increasing the effective interaction area. For illustrative and modelling purposes, the duct inlet radius is considered to be 20% greater than the rotor radius, forming a hydrodynamic funnel that progressively narrows the flow cross-section. The rotor’s effective working area is treated as annular, defined by the difference between the area of a circle with radius (the rotor) and the area of the central hub, assumed to have radius .
These simplified proportions were adopted to facilitate analytical development, providing a coherent basis for estimating the acceleration induced by geometric constraints. Although illustrative, they do not represent fixed design parameters and may be adjusted based on specific engineering requirements. Ultimately, both will be incorporated into a global coefficient that encapsulates the ideal geometric and operational contributions of the turbine configuration within the theoretical formulation.
This parameterisation is now applied to ensure mathematical continuity and to derive a preliminary estimation of the flow acceleration within the rotor’s active region.
3.2.1. Acceleration by the Duct (Venturi Effect)
According to the principle of volumetric flow rate conservation, the relationship between the flow at the duct inlet (with velocity
) and the flow in the rotor region (with velocity
already accelerated due to the duct contraction) is given by:
where the inlet area of the duct is:
and the rotor's frontal area is:
Substituting these areas into Equation (5):
This relationship shows that the velocity in the rotor region increases proportionally to the reduction in cross-sectional area imposed by the duct geometry. This acceleration mechanism is consistent with the Venturi effect, in which the contraction of the conduit’s cross-sectional area promotes a higher flow velocity while maintaining a constant volumetric flow rate.
The assumed ratio of the duct inlet radius to the rotor radius, which leads to the (36/25) factor, was adopted for analytical clarity and does not reflect a definitive design constraint. It will be generalised in the formulation through a global geometric coefficient.
3.2.2. Acceleration Due to Restriction to the Effective Annular Area
In the proposed model, after the redirection of the flow towards the rotor region, it is considered that the effective interaction with the blades occurs only within the annular area between the central hub and the external edge of the rotor. This effective flow area is defined by the difference between the area of the circle with radius
R (the rotor radius) and the area of the circle corresponding to the hub with radius
The effective annular area is given by:
The continuity of volumetric flow rate between the total rotor area and the effective annular area yields the following relationship:
Rearranging to isolate
:
As in the previous section, the use of as the hub radius reflects a modelling simplification adopted to facilitate analytical derivation. This proportion will also be incorporated into the global geometric coefficient in subsequent generalisations.
3.2.3. Additional Hypothesis: Relative Hydrodynamic Confinement Due to Hydrostatic Pressure at Depth
Although the channelling effect promoted by the turbine duct geometry is conceptually based on the classical Venturi principle, it must be acknowledged that this effect is typically observed under conditions of rigid physical confinement, such as in pipes or closed conduits. In open environments, including the seabed or natural channels, there are no material barriers to restrict the lateral flow of water. This raises a legitimate question regarding the actual effectiveness of flow acceleration induced solely by the converging geometry of the duct under such conditions.
However, in submerged conditions and at greater depths, the surrounding fluid exists under high hydrostatic pressure and fully occupies the available space. Any attempt to deviate the flow laterally requires that this surrounding water mass be displaced outward, which demands additional hydrodynamic effort. This inherent resistance to radial deviation may favour the axial conduction of the flow, even in the absence of rigid walls. It is therefore proposed that the fluid environment itself may act as a relative confinement mechanism, potentially contributing to an increase in flow velocity in the rotor region.
Although this hypothesis is consistent with hydrodynamic principles, its actual contribution to turbine performance has yet to be experimentally verified. Nonetheless, it is formally incorporated into the performance estimates developed in this article, with appropriate conceptual caution. Its influence may range from negligible to substantial, as suggested by the geometric relationships previously established.
The uncertainty associated with the effectiveness of hydrodynamic confinement will be addressed at a later stage in the article, when global performance coefficients are introduced. At that point, a plausible interval of turbine performance will be proposed, based on the physical and mathematical boundaries derived from the assumption of flow with or without local acceleration.
3.3. Calculation of the Resultant Force from the Change in Momentum within the Extended Volume
Based on the previously established physical model, the total force generated by the turbine results from the change in linear momentum of the fluid mass mobilised by the upstream pressure field.
This mass corresponds to a volume of water with density
, associated with a cylindrical volume defined by the base area
(the cross-sectional area of the mobilised flow) and the length
, which represents the upstream propagation reach, as defined in Equation (3). Thus, the mobilised mass can be expressed as:
Although the flow velocity remains constant in magnitude as it passes around the turbine blades, the axial component of linear momentum is progressively redirected into the tangential direction due to the interaction with the blades. From the turbine’s perspective, this corresponds to a gradual reduction of the fluid’s axial linear momentum, from to zero, over the interaction time. This process characterises the complete conversion of the axial component into a tangential force applied to the blades, opposite to the direction of the discharged flow.
Accordingly, the resultant force can be derived by applying Newton’s second law, as defined in Equation (2), to the mass of fluid mobilised upstream, as expressed in Equation (12). This yields:
The mathematical cancellation of Δt is justified physically under the steady-state regime assumed in the model. Although the transfer of momentum associated with an individual fluid element occurs over a finite time interval, the resulting pressure field organises itself continuously along the profile of the blades. Thus, the integral of differential forces applied along the curved trajectory of the deflected flow is equivalent to the spatial integral of the contributions from all fluid elements simultaneously occupying different positions along this path at any given instant. This equivalence between temporal and spatial integration is a characteristic of continuous systems in steady state. Therefore, the resultant force is correctly described as a continuous function of fluid density, axial flow velocity, and the coupling length determined by the propagation of the perturbation.
The permanence of the stationary regime is ensured by the energy restoration mechanism described in the ERGF proposal, which maintains the continuity of the axial flow even under the resistance imposed by the rotor, without compromising the hydrodynamic balance of the system.
3.4. Calculation of Torque by Integration of the Distributed Force over the Annular Region
The total force
, as expressed in Equation (13), is distributed over the annular area defined by Equation (9), resulting in a surface force density
, that is, the tangential force per unit area acting on the rotor. Considering the cylindrical symmetry of the system, the direction of the force is tangential to the rotor and can be represented by the unit vector
defined as the angular direction orthogonal to the radius in polar coordinates. Thus, the vectorial surface force density is expressed as:
Although this density has the same dimensionality as pressure (N/m²), it should not be confused with the scalar pressure of the fluid. In the context of this model, represents the distributed tangential force with a defined direction, exerted by the blades over the effective annular area of the rotor, resulting from the hydrodynamic interaction. It is, therefore, a distributed vectorial force responsible for generating torque, rather than an isotropic pressure such as the static pressure of the fluid.
The differential area element of an annular region with radius
r and thickness
dr is given by:
The corresponding differential force acting on this ring is:
In this expression, denotes the scalar magnitude of the total tangential force distributed across the annular area, while indicates its direction in the rotational plane. The differential force , therefore, is tangential to the annular region and proportional to the local radius .
Since
and
are orthogonal in the rotor plane (respectively radial and tangential), the resulting differential torque points in the direction normal to the plane, represented by
:
The total torque is obtained by integrating this differential element over the annular region, from
to
:
Substituting the expression for the area
, as defined in Equation (6), into the force expression given in Equation (13), and then inserting this result into the torque formulation presented in Equation (18), yields the complete expression for the total torque generated by the turbine:
To define the maximum angular velocity
of the turbine under free rotational condition, it is assumed that the blade tips move at the local flow velocity
v″, as expressed in Equation (11). This assumption leads to the following relationship between the linear tip speed and the angular velocity:
This expression will be used in the next section for the calculation of the theoretical power converted by the turbine.
3.5. Estimation of Optimal Power Based on the Characteristic Operating Curve
The instantaneous mechanical power available in the turbine is given by the dot product between the torque vector
and the angular velocity vector
, resulting in a scalar quantity of power
P:
It is important to note that the torque derived in Equation (19) corresponds to the blocked-rotor condition, where the shaft is immobilised and the angular velocity is zero. This condition yields the maximum possible torque but results in no mechanical power, as ω = 0. Conversely, Equation (20) represents the free-spinning condition, in which the angular velocity reaches its maximum and the resistive torque vanishes, also leading to zero power output.
Maximum mechanical power is achieved at an intermediate operating point, where torque and angular velocity are both approximately half of their respective maximum values. This balance reflects the classical behaviour of rotating machines under optimal load conditions. By combining Equations (19) and (20), the ideal power output can be estimated as the product of the average torque and the average angular velocity.
This leads to Equation (22), in which the denominator 4 results from multiplying half of the maximum torque by half of the free rotational speed:
In the transition from Equation (19) to the power expression, the unit vector is omitted as the dot product with the angular velocity vector yields a scalar quantity, reflecting the alignment of the vectors along the axis of rotation.
Simplifying the numerical factors involved in the mathematical development yields the aggregate coefficient:
This coefficient represents the combined contribution of the turbine geometry, the force distribution over the effective working area, and the optimal operating condition. It synthesises the effect of all ideal factors involved in the development of the model.
However, it is important to highlight that the coefficient in its full form considers the total acceleration of the flow up to the effective region of the rotor, estimated as . Since power depends on the square of the local velocity over the turbine blades, this acceleration introduces a multiplicative factor of in the term .
If one adopts the conservative hypothesis of no local acceleration, that is, assuming
, the value of
K must be corrected by the same quadratic factor, resulting in a reduced coefficient:
To incorporate the physical limitations associated with the volumetric attenuation previously discussed, the volumetric efficiency factor
is introduced. The equation for the ideal power then becomes:
This expression can be rewritten in a form analogous to the classical formulation for ideal power according to Betz, by employing the conventional notation
for the cross-sectional area of the flow:
where:
is the volumetric efficiency factor;
is the geometric and dynamic coupling coefficient;
is the fluid density;
c is the propagation speed of the perturbation;
A is the cross-sectional area of the flow;
v is the flow velocity.
Although this equation exhibits quadratic dependence on the free-stream velocity v, it is important to emphasise that the inclusion of c (a velocity term) ensures dimensional consistency of power as a quantity proportional to the cube of velocity. However, unlike the classical model, in which power depends exclusively on v3, the proposed model introduces a distinct physical interpretation, in which power scales with the interaction between the flow and the volumetric propagation effects, represented by cv2.
To illustrate the implications of this formulation and highlight the discrepancy between the classical and volumetric models, the following section presents numerical examples comparing the theoretical predictions of both models for typical turbine parameters.
3.8. Comparative Example between the Classical and Volumetric Models
Consider a turbine with rotor radius , operating in a current of in seawater with density , adopting the parameters , , . The theoretical power estimates provided by the two models are:
This substantial difference should not be interpreted solely as evidence of the superiority of the volumetric formulation. It also highlights the limitations of the classical Betz model when applied to liquids. Although originally developed for idealised media such as incompressible air at low velocities, the classical model does not account for the constraints imposed by physically incompressible media, such as water.
Hence, the discrepancy observed between the two models suggests a potential limitation of classical assumptions, particularly in dense and incompressible media like water. While the Betz model is well-established in aerodynamics, its application in aquatic environments demands further investigation, especially when considering the broader hydrodynamic effects captured by the volumetric model.
This comparison remains theoretical and is based on estimated parameters, especially the volumetric efficiency coefficient which, although conservatively selected, still requires experimental confirmation to validate the predictive relevance of the model.