model EpsilonNTUZ "Test model for the functions epsilon_ntuZ and ntu_epsilonZ" extends Modelica.Icons.Example; import f = Buildings.Fluid.Types.HeatExchangerFlowRegime; Real Z[6] "Ratio of capacity flow rates"; Real epsilon[6] "Heat exchanger effectiveness"; Real eps[6] "Heat exchanger effectiveness"; Real ntu[6] "Number of transfer units"; Real diff[6] "Difference in results"; equation for conf in {Integer(f.ParallelFlow), Integer(f.CounterFlow), Integer(f.CrossFlowUnmixed), Integer(f.CrossFlowCMinMixedCMaxUnmixed), Integer(f.CrossFlowCMinUnmixedCMaxMixed), Integer(f.ConstantTemperaturePhaseChange)} loop Z[conf] = abs(cos(time)); epsilon[conf] = 0.01 + 0.98*abs(sin(time)) * 1/(1+Z[conf]); ntu[conf] = ntu_epsilonZ(epsilon[conf], Z[conf], conf); eps[conf] = epsilon_ntuZ(ntu[conf], Z[conf], conf); diff[conf] = epsilon[conf] - eps[conf]; end for; end EpsilonNTUZ;

model HACoilInside "Test model for HACoilInside" extends Modelica.Icons.Example; Modelica.Blocks.Sources.Sine sine( freqHz=0.1, amplitude=0.063, offset=0.063) "Mass flow rate"; Modelica.Blocks.Sources.Sine sine1( amplitude=10, freqHz=0.1, offset=273.15 + 50) "Temperature"; Buildings.Fluid.HeatExchangers.BaseClasses.HACoilInside hASin( hA_nominal=13*(0.5+1)/0.5, m_flow_nominal=0.063) "Model for convection coefficient"; equation connect(sine1.y, hASin.T); connect(sine.y, hASin.m_flow); end HACoilInside;

model HADryCoil "Test model for dry coil convection coefficient" extends Modelica.Icons.Example; Buildings.Fluid.HeatExchangers.BaseClasses.HADryCoil hADryCoi( m_flow_nominal_w=0.1, UA_nominal=10000, m_flow_nominal_a=2.5) "Calculates an hA value for a dry coil"; Modelica.Blocks.Sources.Ramp TWat( offset=293.15, duration=20, startTime=40, height=40) "Temperature of the water"; Modelica.Blocks.Sources.Ramp TAir( height=-20, offset=293.15, duration=20, startTime=60) "Temperature of the air"; Modelica.Blocks.Sources.Ramp mWat_flow( height=0.1, duration=20, offset=0) "Water mass flow rate"; Modelica.Blocks.Sources.Ramp mAir_flow( duration=20, offset=0, height=2.5, startTime=20) "Air mass flow rate"; Buildings.Fluid.HeatExchangers.BaseClasses.HADryCoil hADryCoiAir( m_flow_nominal_w=0.1, UA_nominal=10000, m_flow_nominal_a=2.5, waterSideFlowDependent=false, waterSideTemperatureDependent=false) "Calculates an hA value for a dry coil"; Buildings.Fluid.HeatExchangers.BaseClasses.HADryCoil hADryCoiWat( m_flow_nominal_w=0.1, UA_nominal=10000, m_flow_nominal_a=2.5, airSideFlowDependent=false, airSideTemperatureDependent=false) "Calculates an hA value for a dry coil"; equation connect(TWat.y, hADryCoi.T_1); connect(hADryCoi.T_2, TAir.y); connect(mAir_flow.y, hADryCoi.m2_flow); connect(mWat_flow.y, hADryCoi.m1_flow); connect(TWat.y, hADryCoiAir.T_1); connect(hADryCoiAir.T_2, TAir.y); connect(mAir_flow.y, hADryCoiAir.m2_flow); connect(mWat_flow.y, hADryCoiAir.m1_flow); connect(TWat.y, hADryCoiWat.T_1); connect(hADryCoiWat.T_2, TAir.y); connect(mAir_flow.y, hADryCoiWat.m2_flow); connect(mWat_flow.y, hADryCoiWat.m1_flow); end HADryCoil;

model HANaturalCylinder "Test model for natural convection outside of a coil" extends Modelica.Icons.Example; Buildings.Fluid.HeatExchangers.BaseClasses.HANaturalCylinder hANatCyl(redeclare package Medium = Buildings.Media.Water, ChaLen=0.1, hA_nominal=10, TFlu_nominal=293.15, TSur_nominal=313.15) "Calculates an hA value for natural convection around a cylinder"; Modelica.Blocks.Sources.Ramp TSur( duration=100, height=50, offset=293.15) "Temperature of the pipe surface"; Modelica.Blocks.Sources.Ramp TFlu( duration=100, startTime=150, height=50, offset=293.15) "Temperature of the surrounding fluid"; equation connect(TSur.y, hANatCyl.TSur); connect(TFlu.y, hANatCyl.TFlu); end HANaturalCylinder;

model HexElementLatent "Model that tests the basic element that is used to built heat exchanger models" extends Modelica.Icons.Example; package Medium_W = Buildings.Media.Water; package Medium_A = Buildings.Media.Air; Buildings.Fluid.Sources.Boundary_pT sin_2( redeclare package Medium = Medium_A, use_p_in=true, use_T_in=true, T=288.15, nPorts=1); Modelica.Blocks.Sources.Ramp PIn( height=20, duration=300, startTime=300, offset=101324.95); Buildings.Fluid.Sources.Boundary_pT sou_2( redeclare package Medium = Medium_A, use_p_in=true, use_T_in=true, T=283.15, nPorts=1); Modelica.Blocks.Sources.Ramp TWat( startTime=1, height=4, duration=300, offset=303.15) "Water temperature"; Modelica.Blocks.Sources.Constant TDb(k=278.15) "Drybulb temperature"; Modelica.Blocks.Sources.Constant POut(k=101325); Buildings.Fluid.Sources.Boundary_pT sin_1( redeclare package Medium = Medium_W, use_p_in=true, T=288.15, nPorts=1); Buildings.Fluid.Sources.Boundary_pT sou_1( redeclare package Medium = Medium_W, use_T_in=true, p=101340, T=293.15, nPorts=1); Buildings.Fluid.FixedResistances.PressureDrop res_22( m_flow_nominal=5, dp_nominal=5, redeclare package Medium = Medium_A); Buildings.Fluid.FixedResistances.PressureDrop res_12( m_flow_nominal=5, dp_nominal=5, redeclare package Medium = Medium_W); Buildings.Fluid.HeatExchangers.BaseClasses.HexElementLatent hex( m1_flow_nominal=5, m2_flow_nominal=5, UA_nominal=9999, redeclare package Medium1 = Medium_W, redeclare package Medium2 = Medium_A, dp1_nominal=10, dp2_nominal=10, energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState); Modelica.Blocks.Sources.Constant TDb1(k=303.15) "Drybulb temperature"; Modelica.Blocks.Sources.Constant hACon(k=10000) "Convective heat transfer"; Sensors.EnthalpyFlowRate senEntFlo( redeclare package Medium = Medium_W, m_flow_nominal=5, tau=0, initType=Modelica.Blocks.Types.Init.SteadyState) "Enthalpy flow rate sensor"; Sensors.EnthalpyFlowRate senEntFlo1( redeclare package Medium = Medium_W, m_flow_nominal=5, tau=0, initType=Modelica.Blocks.Types.Init.SteadyState) "Enthalpy flow rate sensor"; Sensors.EnthalpyFlowRate senEntFlo2( m_flow_nominal=5, tau=0, initType=Modelica.Blocks.Types.Init.SteadyState, redeclare package Medium = Medium_A) "Enthalpy flow rate sensor"; Sensors.EnthalpyFlowRate senEntFlo3( m_flow_nominal=5, tau=0, initType=Modelica.Blocks.Types.Init.SteadyState, redeclare package Medium = Medium_A) "Enthalpy flow rate sensor"; Modelica.Blocks.Math.MultiSum netEne(nu=4, k={1,-1,1,-1}) "Net energy balance (needs to be zero)"; equation connect(TDb.y, sou_1.T_in); connect(POut.y, sin_1.p_in); connect(PIn.y, sou_2.p_in); connect(TWat.y, sou_2.T_in); connect(POut.y, sin_2.p_in); connect(TDb1.y, sin_2.T_in); connect(hACon.y, hex.Gc_1); connect(hACon.y, hex.Gc_2); connect(sin_1.ports[1], res_12.port_b); connect(sin_2.ports[1], res_22.port_b); connect(senEntFlo.port_b, hex.port_a1); connect(senEntFlo.port_a, sou_1.ports[1]); connect(res_12.port_a, senEntFlo1.port_b); connect(hex.port_b1, senEntFlo1.port_a); connect(res_22.port_a, senEntFlo2.port_a); connect(senEntFlo2.port_b, hex.port_b2); connect(hex.port_a2, senEntFlo3.port_a); connect(senEntFlo3.port_b, sou_2.ports[1]); connect(netEne.u[1], senEntFlo.H_flow); connect(senEntFlo1.H_flow, netEne.u[2]); connect(senEntFlo2.H_flow, netEne.u[3]); connect(senEntFlo3.H_flow, netEne.u[4]); end HexElementLatent;

model HexElementLatentLoop "Model that tests the basic element that is used to built heat exchanger models" extends Modelica.Icons.Example; package Medium_W = Buildings.Media.Water; package Medium_A = Buildings.Media.Air; parameter Modelica.SIunits.MassFlowRate mW_flow_nominal = 0.1 "Water mass flow rate"; parameter Modelica.SIunits.MassFlowRate mA_flow_nominal = 0.14 "Air mass flow rate"; Buildings.Fluid.Sources.Boundary_pT sin_2( redeclare package Medium = Medium_A, use_p_in=false, use_T_in=false, nPorts=1) "Pressure boundary condition"; Buildings.Fluid.Sources.Boundary_pT sin_1( redeclare package Medium = Medium_W, nPorts=1) "Pressure boundary condition"; Sources.MassFlowSource_T sou( redeclare package Medium = Medium_W, m_flow=5, nPorts=1, use_T_in=true) "Mass flow source"; Buildings.Fluid.HeatExchangers.BaseClasses.HexElementLatent hex( redeclare package Medium1 = Medium_W, redeclare package Medium2 = Medium_A, dp2_nominal=5, UA_nominal=1E3, m1_flow_nominal=mW_flow_nominal, m2_flow_nominal=mA_flow_nominal, dp1_nominal=0, show_T=true, energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState) "Hex element"; Modelica.Blocks.Sources.Constant hACon(k=70) "Convective heat transfer"; Movers.FlowControlled_m_flow fan( redeclare package Medium = Medium_A, m_flow_nominal=mA_flow_nominal, inputType=Buildings.Fluid.Types.InputType.Constant, addPowerToMedium=false, nominalValuesDefineDefaultPressureCurve=true, dp_nominal=5, use_inputFilter=false, energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial) "Fan"; Buildings.Utilities.Psychrometrics.pW_TDewPoi TDewPoi "Model to compute the water vapor pressure at the dew point"; public Buildings.Utilities.Psychrometrics.X_pW humRatPre(use_p_in=false) "Model to convert water vapor pressure into humidity ratio"; Modelica.Blocks.Sources.Constant TWat(k=273.15 + 5) "Water inlet temperature"; Buildings.Fluid.HeatExchangers.PrescribedOutlet hea( redeclare package Medium = Medium_A, m_flow_nominal=mA_flow_nominal, dp_nominal=0, energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState, use_X_wSet=false) "Heater"; Modelica.Blocks.Sources.Constant TAir(k=273.15 + 20) "Air inlet temperature"; Modelica.Blocks.Sources.RealExpression TSur(y=hex.masExc.TSur) "Surface temperature of heat exchanger"; equation connect(hACon.y, hex.Gc_1); connect(hACon.y, hex.Gc_2); connect(sou.ports[1], hex.port_a1); connect(sin_1.ports[1], hex.port_b1); connect(fan.port_a, hex.port_b2); connect(sin_2.ports[1], fan.port_a); connect(TDewPoi.p_w, humRatPre.p_w); connect(TWat.y, sou.T_in); connect(fan.port_b, hea.port_a); connect(hea.port_b, hex.port_a2); connect(TAir.y, hea.TSet); connect(TSur.y, TDewPoi.T); end HexElementLatentLoop;

model HexElementSensible "Model that tests the basic element that is used to built heat exchanger models" extends Modelica.Icons.Example; package Medium_W = Buildings.Media.Water; package Medium_A = Buildings.Media.Air; Buildings.Fluid.Sources.Boundary_pT sin_2( redeclare package Medium = Medium_A, use_p_in=true, use_T_in=true, T=288.15, nPorts=1); Modelica.Blocks.Sources.Ramp PIn( height=20, duration=300, startTime=300, offset=101325); Buildings.Fluid.Sources.Boundary_pT sou_2( redeclare package Medium = Medium_A, use_p_in=true, use_T_in=true, T=283.15, nPorts=1); Modelica.Blocks.Sources.Ramp TWat( startTime=1, height=4, duration=300, offset=303.15) "Water temperature"; Modelica.Blocks.Sources.Constant TDb(k=278.15) "Drybulb temperature"; Modelica.Blocks.Sources.Constant POut(k=101325); Buildings.Fluid.Sources.Boundary_pT sin_1( redeclare package Medium = Medium_W, use_p_in=true, T=288.15, nPorts=1); Buildings.Fluid.Sources.Boundary_pT sou_1( redeclare package Medium = Medium_W, use_T_in=true, nPorts=1, p=101340, T=293.15); Buildings.Fluid.FixedResistances.PressureDrop res_22( m_flow_nominal=5, dp_nominal=5, redeclare package Medium = Medium_A); Buildings.Fluid.FixedResistances.PressureDrop res_12( m_flow_nominal=5, dp_nominal=5, redeclare package Medium = Medium_W); Buildings.Fluid.HeatExchangers.BaseClasses.HexElementSensible hex( m1_flow_nominal=5, m2_flow_nominal=5, UA_nominal=9999, redeclare package Medium1 = Medium_W, redeclare package Medium2 = Medium_A, dp1_nominal=5, dp2_nominal=5, energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial); Modelica.Blocks.Sources.Constant TDb1(k=303.15) "Drybulb temperature"; Buildings.Fluid.FixedResistances.PressureDrop res_11( m_flow_nominal=5, dp_nominal=5, redeclare package Medium = Medium_W); Buildings.Fluid.FixedResistances.PressureDrop res_21( m_flow_nominal=5, dp_nominal=5, redeclare package Medium = Medium_A); Modelica.Blocks.Sources.Constant hACon(k=10000) "Convective heat transfer"; equation connect(TDb.y, sou_1.T_in); connect(POut.y, sin_1.p_in); connect(PIn.y, sou_2.p_in); connect(TWat.y, sou_2.T_in); connect(POut.y, sin_2.p_in); connect(hex.port_b1, res_12.port_a); connect(res_22.port_a, hex.port_b2); connect(TDb1.y, sin_2.T_in); connect(res_11.port_b, hex.port_a1); connect(hex.port_a2, res_21.port_b); connect(hACon.y, hex.Gc_1); connect(hACon.y, hex.Gc_2); connect(sou_2.ports[1], res_21.port_a); connect(sin_1.ports[1], res_12.port_b); connect(sou_1.ports[1], res_11.port_a); connect(sin_2.ports[1], res_22.port_b); end HexElementSensible;

model MassExchange "Test model for latent heat exchange" extends Modelica.Icons.Example; package Medium = Buildings.Media.Air; Buildings.Fluid.HeatExchangers.BaseClasses.MassExchange masExc(redeclare package Medium = Medium) "Model for mass exchange"; Modelica.Blocks.Sources.Ramp TSur( duration=1, height=20, offset=273.15 + 5) "Surface temperature"; Modelica.Blocks.Sources.Constant X_w(k=0.01) "Humidity mass fraction in medium"; Modelica.Blocks.Sources.Constant Gc(k=1) "Sensible convective thermal conductance"; equation connect(TSur.y, masExc.TSur); connect(Gc.y, masExc.Gc); connect(X_w.y, masExc.XInf); end MassExchange;

model RayleighNumber "Test model for RayleighNumber" extends Modelica.Icons.Example; Buildings.Fluid.HeatExchangers.BaseClasses.RayleighNumber rayleighNumber(ChaLen= 0.1, redeclare package Medium = Buildings.Media.Water); Modelica.Blocks.Sources.Sine TSur( amplitude=50, freqHz=1/50, offset=273.15 + 50) "Temperature of the hot surface"; Modelica.Blocks.Sources.Sine TFlu( amplitude=50, freqHz=1/25, offset=273.15 + 50) "Temperature of the surrounding fluid"; equation connect(TSur.y, rayleighNumber.TSur); connect(TFlu.y, rayleighNumber.TFlu); end RayleighNumber;