- A_refReference area [m^2]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Reference area [m^2]
- D_wheelAny physical dimension of the turbine as a measure of the turbine's size, usually the rotor diameter [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Any physical dimension of the turbine as a measure of the turbine's size, usually the rotor diameter [m]
- head_coefficientHead coefficient vs flow coefficient function [-]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Head coefficient vs flow coefficient function [-]
- inertia_coeffTurbine inertia coefficients [kg-m^2]
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Turbine inertia coefficients [kg-m^2]
- inertia_constTurbine inertia constant [kg-m^2]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Turbine inertia constant [kg-m^2]
- inletTurbine inlet
C++ Type:BoundaryName
Unit:(no unit assumed)
Controllable:No
Description:Turbine inlet
- omega_ratedRated turbine speed [rad/s]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Rated turbine speed [rad/s]
- outletTurbine outlet
C++ Type:BoundaryName
Unit:(no unit assumed)
Controllable:No
Description:Turbine outlet
- positionSpatial position of the center of the junction [m]
C++ Type:libMesh::Point
Unit:(no unit assumed)
Controllable:No
Description:Spatial position of the center of the junction [m]
- power_coefficientPower coefficient vs flow coefficient function [-]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Power coefficient vs flow coefficient function [-]
- speed_cr_ITurbine speed threshold for inertia [-]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Turbine speed threshold for inertia [-]
- speed_cr_frTurbine speed threshold for friction [-]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Turbine speed threshold for friction [-]
- tau_fr_coeffTurbine friction coefficients [N-m]
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Turbine friction coefficients [N-m]
- tau_fr_constTurbine friction constant [N-m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Turbine friction constant [N-m]
- volumeVolume of the junction [m^3]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Volume of the junction [m^3]
ShaftConnectedTurbine1Phase
The 1-phase shaft-connected turbine component is modeled as a volume junction with source terms added to the momentum and energy equations due to the turbine's torque and pressure drop.
The turbine model relies on input parameters, upstream flow conditions, and rotational speed to calculate the turbine pressure drop, driving torque, friction torque, power, and moment of inertia. Turbine pressure drop and net torque (sum of driving and friction torques) are used in source terms to affect flow conditions coming out of the turbine. Net torque and moment of inertia are sent to the connected Shaft
component and affect the rotational speed.
Usage
A ShaftConnectedTurbine1Phase
component must be connected to a Shaft component, which controls the turbine rotational speed. The user must specify which flow channel boundaries are connected to the turbine with the "inlet" and "outlet" parameters.
"A_ref" should generally be assigned a value in between the "A" values of the connected inlet
or outlet
. A_ref
can act as a means to scale the turbine pressure drop. A_ref
can also be used in combination with "K" to apply an additional pressure drop due to form loss. The turbine volume, "volume", should be a reasonable estimate of the inlet to outlet fluid volume of the modeled turbine.
The user assigns values for the rated turbine speed, "omega_rated", and the turbine wheel diameter, "D_wheel". D_wheel
can be any physical dimension of the turbine as a measure of the turbine's size, usually the rotor diameter. These rated values, along with the user-supplied functions "head_coefficient" and "power_coefficient", define the turbine performance model. Further discussion on these parameters, their defining equations, and the input data used in them is found in Turbine performance data.
The user inputs values for "tau_fr_coeff", "tau_fr_const", and "speed_cr_fr" that are used to compute friction torque as defined in Friction torque.
The turbine moment of inertia affects how the turbine accelerates. The user inputs values for "inertia_coeff", "inertia_const", and "speed_cr_I" that are used to compute the turbine moment of inertia as defined in Moment of inertia.
Input Parameters
- K0Form loss factor [-]
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:Yes
Description:Form loss factor [-]
- apply_velocity_scalingFalseSet to true to apply the scaling to the normal velocity. See documentation for more information.
Default:False
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Set to true to apply the scaling to the normal velocity. See documentation for more information.
- initial_TInitial temperature [K]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Initial temperature [K]
- initial_pInitial pressure [Pa]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Initial pressure [Pa]
- initial_vel_xInitial velocity in x-direction [m/s]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Initial velocity in x-direction [m/s]
- initial_vel_yInitial velocity in y-direction [m/s]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Initial velocity in y-direction [m/s]
- initial_vel_zInitial velocity in z-direction [m/s]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Initial velocity in z-direction [m/s]
- scaling_factor_rhoEV1Scaling factor for rho*E*V [-]
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for rho*E*V [-]
- scaling_factor_rhoV1Scaling factor for rho*V [-]
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for rho*V [-]
- scaling_factor_rhouV1Scaling factor for rho*u*V [-]
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for rho*u*V [-]
- scaling_factor_rhovV1Scaling factor for rho*v*V [-]
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for rho*v*V [-]
- scaling_factor_rhowV1Scaling factor for rho*w*V [-]
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for rho*w*V [-]
- use_scalar_variablesTrueTrue if the junction variables are scalar variables
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:True if the junction variables are scalar variables
Optional Parameters
- control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
- enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Unit:(no unit assumed)
Controllable:No
Description:Set the enabled status of the MooseObject.
Advanced Parameters
Formulation
The conservation of mass, momentum, and energy equations for the turbine volume are similar to the equations used by VolumeJunction1Phase but add the turbine momentum and energy source terms,
(1)and
(2)where
is the pressure drop from turbine inlet to outlet,
is the cross-sectional area of the turbine,
is the orientation of the flow channel connected to the turbine outlet,
is the driving torque,
is the friction torque, and
is the shaft speed.
Turbine performance data
The turbine model is based on dimensional analysis using Buckingham's Pi Theorem Anon (2021). "The Affinity Laws were derived using the dimensional analysis methods developed by Buckingham to reduce the complex performance map into distinct dimensionless curves which are more useful output parameters. For a single-phase gas turbine, the turbine's performance parameters are non-dimensionalized to evaluate the effects of flow rate (represented by the flow coefficient, ), fluid density (), rotational speed () and turbine wheel diameter () on the power input represented by the head coefficient () and on useful turbine work represented by the power output coefficient () Patil et al. (2019)."
Values are obtained by means of an empirical turbine performance map. The independent variables are flow coefficient () and shaft speed (). The dependent variables are the head coefficient () and power coefficient ().
The flow coefficient is computed as
(3)where is volumetric flow rate at the turbine inlet.
The head coefficient is computed as
(4)The power coefficient is computed as
(5)where is the fluid density within the turbine volume.
Pressure drop () and driving torque () are calculated using the input turbine performance data, dimensional coefficients defined above, and the two equations below.
(6)(7)Friction torque
Friction torque always resists the direction of motion, so the sign of the turbine friction torque depends on the sign of the connected shaft speed. Positive shaft speed results in negative friction torque while negative shaft speed results in positive friction torque.
The friction torque magnitude, , is a function of shaft speed, , and four input parameters. If the ratio of shaft speed to "omega_rated", , is less than "speed_cr_fr", friction torque magnitude equals "tau_fr_const",
otherwise, is a function of shaft speed and friction torque coefficients, "tau_fr_coeff",
Moment of inertia
The turbine moment of inertia, , is a function of shaft speed, , and four input parameters. If the ratio of shaft speed to "omega_rated", , is less than "speed_cr_I", turbine inertia equals "inertia_const",
otherwise, is a function of shaft speed and inertia coefficients, "inertia_coeff",
References
- Anon.
Basic principles of turbomachines: lecture 3.
2021.
\url https://nptel.ac.in/courses/112/104/112104117/.[BibTeX]
- A. Patil, S. Sundar, Y. Wang, M. Solom, K.V. Kirkland, and G. Morrison.
Characterization of steam impulse turbine for two-phase flow.
Int. J. Heat Fluid Flow, pages, 2019.
doi:10.1016/j.ijheatfluidflow.2019.108439.[BibTeX]