modelFrancis

Model of the Francis turbine

Extends from Icons.Turbine (Turbine icon), OpenHPL.Interfaces.TwoContacts (Model of two connectors), OpenHPL.Interfaces.TurbineContacts (Model of turbine connectors).

Information

Francis Turbine Model

This is the Francis turbine model that gives possibilities for proper modelling of the Francis turbine. The mechanistic model is based on Euler equations for the Francis turbine.

Besides hydraulic input and output, there are input as the control signal for the valve opening and also output as the turbine shaft power and input as angular velocity.

Figure: Key quantities in the Francis turbine model showing inlet (1) and outlet (2) sections.

Euler Turbine Equations

The shaft power produced by the turbine is given by:

$$ \dot{W}_s = \dot{m}\omega \left(R_1\frac{\dot{V}}{A_1}\cot{\alpha_1} - R_2\left(\omega R_2 + \frac{\dot{V}}{A_2}\cot{\beta_2}\right)\right) $$

where ṁ and V̇ are mass and volumetric flow rates, ω is angular velocity, R1 and R2 are inlet and outlet radii, A1 and A2 are cross-sectional areas, α1 is inlet guide vane angle, and β2 is outlet blade angle.

Total Work and Efficiency

The total work rate is:

$$ \dot{W}_t = \dot{W}_s + \dot{W}_{ft} + \Delta p_v \dot{V} $$

where Ẇft represents various friction losses (shock, whirl, wall friction), and ΔpvV̇ accounts for guide vane pressure drop. Turbine efficiency \(\eta = \dot{W}_s / \dot{W}_t\).

Turbine Design Algorithm

There is also available the runner design algorithm that can define all geometrical parameters based on the nominal parameters (net head, flow rate, power, speed). The algorithm determines: outlet blade angle β2, runner radii R1 and R2, runner width w1, and inlet blade angle β1.

The turbine losses coefficients (k_ft1, k_ft2, k_ft3) can be also defined automatically. However, if some dynamic data from real turbine is available it is better to tune these parameters a bit more and use the defined values as a starting point.

Guide Vane Actuation

A model for servo that runs the guide vane opening is also available. A guide vane opening model relates actuator position Y to guide vane angle α1 through geometric relationships. Furthermore it is possible to automatically generate all needed parameters for the servo, or simply specify them.

Low Load Performance

This mechanistic turbine model does not work really well for low loads (<10% guide vane opening). However there is parameters that could be tuned for low load regimes. These are u_min and k_ft4.

More Information

More info about the mechanistic turbine model can be found in: [Vytvytskyi2018] and [Vytvytskyi2019]. Additional details about the servo (and turbine model) are in: Resources/Documents/Turbines_model.pdf.

Parameters

TypeNameDefaultDescription
Given data
BooleanGivenDatatrueIf checked the user specifies whole set of parameters. Otherwise, the design algorithm will be used for missed parameters
Guide vane
Booleandp_v_conditionfalseIf checked then it includes the pressure drop through the guide vane (leave it unchecked - this doesn't work well)
SI.RadiusR_v_2.89/2Radius of the guide vane suspension circle
SI.Lengthw_v_w_1_Width of the guide vane suspension circle
Realk_fv0Friction loss coefficient from turbine entrance and across guide vanes
RealReduction0.2Reduction the given formula for the guide vane pressure drop
Servo
BooleanGivenServoDatatrueIf checked the user specifies parameters for the servo. Otherwise, the design algorithm will be used for missed parameters
SI.Lengthr_v_1.1Radius of servo circle
SI.Lengthr_Y_1.2Radius to servo connection
SI.LengthR_Y_3Radius to servo
Realu_start_2.28Servo position with zero flow
Realu_end_2.4Servo position with full flow
Nominal parameters
SI.HeightH_n460Nominal head
SI.VolumeFlowRateVdot_n24.3Nominal flow
SI.PowerP_n103e6Nominal power
Modelica.Units.NonSI.AngularVelocity_rpmn_n500Nominal turbine speed
Runner
SI.RadiusR_1_2.63/2Radius of the turbine blade inlet
SI.RadiusR_2_1.55/2Radius of the turbine blade outlet
SI.Lengthw_1_0.2Width of the turbine/blades inlet
Modelica.Units.NonSI.Angle_degbeta1_110Turbine inlet blade angle
Modelica.Units.NonSI.Angle_degbeta2_162.5Turbine outlet blade angle
Guide vane › Scroll case
SI.DiameterD_i1.632Diameter of the inlet pipe
Losses in runner
BooleanGiven_lossestrueFriction shock loss coefficient due to shock
Realk_ft1_7e5Friction shock loss coefficient due to shock
Realk_ft2_0Hydraulic friction loss coefficient due to effluent whirl
Realk_ft3_1.63e4Friction loss coefficient due to wall friction
Parameters for low load
Realk_ft4k_ft3_*100Friction loss coefficient that is used for low load (under u_min)
Realu_min0.03Control signal value under which the moodel used k_f4 friction term to balance the model
Condition
BooleanWaterCompressfalseIf checked the water is compressible in the penstock
I/O › Outputs
Booleanenable_P_out (from TurbineContacts)falseIf checked, get a connector for the output power

Connectors

TypeNameDefaultDescription
Contact_ii (from TwoContacts)Inlet contact (positive design flow direction is from i to o)
Contact_oo (from TwoContacts)Outlet contact (positive design flow direction is from i to o)
Modelica.Blocks.Interfaces.RealInputu_t (from TurbineContacts)[Guide vane|nozzle] opening of the turbine(=1: completely open, =0: completely closed)
Modelica.Blocks.Interfaces.RealOutputP_out (from TurbineContacts)Mechanical Output power
Modelica.Blocks.Interfaces.RealInputw_inInput angular velocity from the generator

Components

TypeNameDefaultDescription
DatadataUsing standard class with constants
SI.Pressurep_r1Runner inlet pressure
SI.Pressuredp_trTurbine pressure drop
SI.Pressuredp_rRunner pressure drop
SI.Pressurep_tr2Turbine outlet pressure
SI.Pressuredp_vGuide vane pressure drop
SI.AreaA_1Runner inlet cross section
SI.AreaA_0Turbine inlet cross section
SI.AreaA_vGuide vane cross section
SI.AreaA_2Runner outlet cross section
SI.EnergyFlowRateWdot_sShaft power
SI.EnergyFlowRateWdot_ftTotal runner losses
SI.EnergyFlowRateW_t1Euler first term
SI.EnergyFlowRateW_t2Euler second term
SI.EnergyFlowRateWdot_ft_sShock losses
SI.EnergyFlowRateWdot_ft_wWhirl losses
SI.EnergyFlowRateWdot_ft_lFriction losses
SI.EnergyFlowRateWdot_tTotal power
SI.VolumeFlowRateVdotFlow rate
SI.MassFlowRatemdotMass flow rate
SI.AngularVelocityww_inAngular velocity
SI.Velocityu_2Outlet reference velocity
SI.Velocityc_m2Outlet meridional velocity
SI.Velocityc_m1Inlet meridional velocity
SI.Velocityu_1Inlet reference velocity
SI.Velocityc_u1Inlet tangential velocity
Modelica.Units.NonSI.Angle_degbeta1Inlet blade angle
Modelica.Units.NonSI.Angle_degbeta2Outlet blade angle
Modelica.Units.NonSI.Angle_deg_beta1
SI.Anglealpha1Inlet guide vane angle
SI.AnglephiOne of servo angles
SI.AnglepsiOne of servo angles
SI.AnglethetaOne of servo angles
SI.AngledthetaOne of servo angles
Realk_ft1
Realk_ft2
Realk_ft3
Realcot_a1
Realcot_a2
Realcot_b1
Realcot_b2
Realcot_g1
Realsin_a1
Realcoef
SI.Lengthl1.8*sqrt((R_v - r_v)^2/2)Servo term
SI.LengthdServo term
SI.LengthR_1Inlet runner radius
SI.LengthR_2Outlet runner radius
SI.LengthR_vGuide vane radius
SI.Lengthw_1Inlet runner width/geight
SI.Lengthw_vGuide vane width/geight
SI.Lengthr_YServo term
SI.LengthR_YServo term
SI.Lengthr_vServo term
RealYServo position
RealY0sqrt(R_Y^2 - r_Y^2)Initial servo position
Realtheta0Modelica.Math.acos(r_Y/R_Y)Initial servo angle
Reald0_2l*(r_v^2 - R_v^2)/(l - R_v)Initial servo term d
Realtheta_0theta0 - Modelica.Math.acos((r_v^2 + R_v^2 - d0_2)/(2*r_v*R_v))Servo angle for fully close guide vane
Realu_endServo position for fully open guide vane
Realu_startServo position for fully close guide vane
RealW_t2_nEuler second term, nominal
RealW_t1_nEuler first term, nominal
RealWdot_t_nTotal power, nominal
Realcot_a1_nCotant nominal alpha
RealVdot_n_Vdot_n/0.99Flow rate for fully open guide vane
Reald_nNominal servo term
Realtheta_nServo angle for fully open guide vane
SI.Anglealpha1_nNominal inlet guide vane angle