model FlowControlled_dp "Fan or pump with ideally controlled head dp as input signal" extends Buildings.Fluid.Movers.BaseClasses.PartialFlowMachine( final preVar=Buildings.Fluid.Movers.BaseClasses.Types.PrescribedVariable.PressureDifference, final computePowerUsingSimilarityLaws=per.havePressureCurve, preSou(dp_start=dp_start, control_dp= not prescribeSystemPressure), final stageInputs(each final unit="Pa") = heads, final constInput(final unit="Pa") = constantHead, final _m_flow_nominal = m_flow_nominal, filter( final y_start=dp_start, u(final unit="Pa"), y(final unit="Pa"), x(each nominal=dp_nominal), u_nominal=dp_nominal), eff( per( final pressure= if per.havePressureCurve then per.pressure else Buildings.Fluid.Movers.BaseClasses.Characteristics.flowParameters( V_flow = {i/(nOri-1)*2.0*m_flow_nominal/rho_default for i in 0:(nOri-1)}, dp = {i/(nOri-1)*2.0*dp_nominal for i in (nOri-1):-1:0}), final etaHydMet= if (per.etaHydMet == Buildings.Fluid.Movers.BaseClasses.Types.HydraulicEfficiencyMethod.Power_VolumeFlowRate or per.etaHydMet == Buildings.Fluid.Movers.BaseClasses.Types.HydraulicEfficiencyMethod.EulerNumber) and not per.havePressureCurve then Buildings.Fluid.Movers.BaseClasses.Types.HydraulicEfficiencyMethod.NotProvided else per.etaHydMet, final etaMotMet= if (per.etaMotMet == Buildings.Fluid.Movers.BaseClasses.Types.MotorEfficiencyMethod.Efficiency_MotorPartLoadRatio or per.etaMotMet == Buildings.Fluid.Movers.BaseClasses.Types.MotorEfficiencyMethod.GenericCurve) and (not per.haveWMot_nominal and not per.havePressureCurve) then Buildings.Fluid.Movers.BaseClasses.Types.MotorEfficiencyMethod.NotProvided else per.etaMotMet), r_N(start=if abs(dp_nominal) > 1E-8 then dp_start/dp_nominal else 0))); parameter Modelica.Units.SI.PressureDifference dp_start( min=0, displayUnit="Pa") = 0 "Initial value of pressure raise"; parameter Modelica.Units.SI.MassFlowRate m_flow_nominal( final min=Modelica.Constants.small) "Nominal mass flow rate"; // For air, we set dp_nominal = 600 as default, for water we set 10000 parameter Modelica.Units.SI.PressureDifference dp_nominal( final min=Modelica.Constants.small, displayUnit="Pa") = if rho_default < 500 then 500 else 10000 "Nominal pressure raise, used to normalized the filter if use_inputFilter=true, to set default values of constantHead and heads, and and for default pressure curve if not specified in record per"; parameter Modelica.Units.SI.PressureDifference constantHead( min=0, displayUnit="Pa") = dp_nominal "Constant pump head, used when inputType=Constant"; // By default, set heads proportional to sqrt(speed/speed_nominal) parameter Modelica.Units.SI.PressureDifference[:] heads( each min=0, each displayUnit="Pa") = dp_nominal*{(per.speeds[i]/per.speeds[end])^2 for i in 1:size(per.speeds, 1)} "Vector of head set points, used when inputType=Stages"; parameter Boolean prescribeSystemPressure = false "=true, to control mover such that pressure difference is obtained across two remote points in system"; Modelica.Blocks.Interfaces.RealInput dpMea( final quantity="PressureDifference", final displayUnit="Pa", final unit="Pa")=gain.u if prescribeSystemPressure "Measurement of pressure difference between two points where the set point should be obtained"; Modelica.Blocks.Interfaces.RealInput dp_in(final unit="Pa") if inputType == Buildings.Fluid.Types.InputType.Continuous "Prescribed pressure rise"; Modelica.Blocks.Interfaces.RealOutput dp_actual(final unit="Pa") "Pressure difference between the mover inlet and outlet"; protected Modelica.Blocks.Math.Gain gain(final k=-1); equation assert(inputSwitch.u >= -1E-3, "Pressure set point for mover cannot be negative. Obtained dp = " + String(inputSwitch.u)); if use_inputFilter then connect(filter.y, gain.u); else connect(inputSwitch.y, gain.u); end if; connect(inputSwitch.u, dp_in); connect(preSou.dp_in, gain.y); connect(senRelPre.p_rel, dp_actual); end FlowControlled_dp;

model FlowControlled_m_flow "Fan or pump with ideally controlled mass flow rate as input signal" extends Buildings.Fluid.Movers.BaseClasses.PartialFlowMachine( final preVar=Buildings.Fluid.Movers.BaseClasses.Types.PrescribedVariable.FlowRate, final computePowerUsingSimilarityLaws=per.havePressureCurve, final stageInputs(each final unit="kg/s")=massFlowRates, final constInput(final unit="kg/s")=constantMassFlowRate, final _m_flow_nominal = m_flow_nominal, filter( final y_start=m_flow_start, u(final unit="kg/s"), y(final unit="kg/s"), x(each nominal=m_flow_nominal), u_nominal=m_flow_nominal), eff( per( final pressure = if per.havePressureCurve then per.pressure else Buildings.Fluid.Movers.BaseClasses.Characteristics.flowParameters( V_flow = {i/(nOri-1)*2.0*m_flow_nominal/rho_default for i in 0:(nOri-1)}, dp = {i/(nOri-1)*2.0*dp_nominal for i in (nOri-1):-1:0}), final etaHydMet= if (per.etaHydMet == Buildings.Fluid.Movers.BaseClasses.Types.HydraulicEfficiencyMethod.Power_VolumeFlowRate or per.etaHydMet == Buildings.Fluid.Movers.BaseClasses.Types.HydraulicEfficiencyMethod.EulerNumber) and not per.havePressureCurve then Buildings.Fluid.Movers.BaseClasses.Types.HydraulicEfficiencyMethod.NotProvided else per.etaHydMet, final etaMotMet= if (per.etaMotMet == Buildings.Fluid.Movers.BaseClasses.Types.MotorEfficiencyMethod.Efficiency_MotorPartLoadRatio or per.etaMotMet == Buildings.Fluid.Movers.BaseClasses.Types.MotorEfficiencyMethod.GenericCurve) and (not per.haveWMot_nominal and not per.havePressureCurve) then Buildings.Fluid.Movers.BaseClasses.Types.MotorEfficiencyMethod.NotProvided else per.etaMotMet), r_N(start=if abs(m_flow_nominal) > 1E-8 then m_flow_start/m_flow_nominal else 0)), preSou(m_flow_start=m_flow_start)); parameter Modelica.Units.SI.MassFlowRate m_flow_nominal( final min=Modelica.Constants.small) "Nominal mass flow rate"; // For air, we set dp_nominal = 600 as default, for water we set 10000 parameter Modelica.Units.SI.PressureDifference dp_nominal( final min=Modelica.Constants.small, displayUnit="Pa") = if rho_default < 500 then 500 else 10000 "Nominal pressure raise, used for default pressure curve if not specified in record per"; parameter Modelica.Units.SI.MassFlowRate m_flow_start(min=0) = 0 "Initial value of mass flow rate"; parameter Modelica.Units.SI.MassFlowRate constantMassFlowRate=m_flow_nominal "Constant pump mass flow rate, used when inputType=Constant"; // By default, set massFlowRates proportional to (speed/speed_nominal) parameter Modelica.Units.SI.MassFlowRate[:] massFlowRates=m_flow_nominal*{per.speeds[ i]/per.speeds[end] for i in 1:size(per.speeds, 1)} "Vector of mass flow rate set points, used when inputType=Stage"; parameter Modelica.Units.SI.Pressure dpMax( min=0, displayUnit="Pa") = 2*max(eff.per.pressure.dp) "Maximum pressure allowed to operate the model, if exceeded, the simulation stops with an error"; Modelica.Blocks.Interfaces.RealInput m_flow_in( final unit="kg/s", nominal=m_flow_nominal) if inputType == Buildings.Fluid.Types.InputType.Continuous "Prescribed mass flow rate"; Modelica.Blocks.Interfaces.RealOutput m_flow_actual( final unit="kg/s", nominal=m_flow_nominal) "Actual mass flow rate"; equation assert(-dp <= dpMax, "In " + getInstanceName() + ": Model operates with head -dp = " + String(-dp) + " Pa, exceeding the pressure allowed by the parameter " + getInstanceName() + ".dpMax. This can happen if the model forces a high mass flow rate through a closed actuator, or if the performance record is unreasonable. Please verify your model, and consider using one of the other pump or fan models."); if use_inputFilter then connect(filter.y, m_flow_actual); connect(filter.y, preSou.m_flow_in); else connect(inputSwitch.y, m_flow_actual); connect(inputSwitch.y, preSou.m_flow_in); end if; connect(inputSwitch.u, m_flow_in); end FlowControlled_m_flow;

model SpeedControlled_y "Fan or pump with ideally controlled normalized speed y as input signal" extends Buildings.Fluid.Movers.BaseClasses.PartialFlowMachine( final preVar=Buildings.Fluid.Movers.BaseClasses.Types.PrescribedVariable.Speed, final nominalValuesDefineDefaultPressureCurve=false, final computePowerUsingSimilarityLaws=true, final stageInputs(each final unit="1") = per.speeds, final constInput(final unit="1") = per.constantSpeed, filter( final y_start=y_start, u(final unit="1"), y(final unit="1")), eff( per(final pressure = per.pressure, final etaHydMet = per.etaHydMet, final etaMotMet = per.etaMotMet), r_N(start=y_start))); parameter Real y_start(min=0, max=1, unit="1")=0 "Initial value of speed"; Modelica.Blocks.Interfaces.RealInput y( unit="1") if inputType == Buildings.Fluid.Types.InputType.Continuous "Constant normalized rotational speed"; protected Modelica.Blocks.Math.Gain gain(final k=-1) "Pressure gain"; initial equation assert(per.havePressureCurve, "SpeedControlled_y requires to set the pressure vs. flow rate curve in record 'per'."); equation connect(inputSwitch.u, y); connect(eff.dp, gain.u); connect(gain.y, preSou.dp_in); if use_inputFilter then connect(filter.y, eff.y_in); else connect(inputSwitch.y, eff.y_in); end if; end SpeedControlled_y;