adding faulty mode to CCG specific components (#279)

helenev-met · web-flow · commit 986a9c969e47 · 2022-09-13T17:00:12.000+02:00
* adding faulty mode to CCG specific components

* no more max value for filter_Kfr fouling

* again : changing Kfr ;p
diff --git a/MetroscopeModelingLibrary/MultiFluid/HeatExchangers/Evaporator.mo b/MetroscopeModelingLibrary/MultiFluid/HeatExchangers/Evaporator.mo
@@ -31,6 +31,10 @@ model Evaporator
     Units.Temperature T_cold_out;
     Units.Temperature T_hot_out;
 
+    // Failure modes
+    parameter Boolean faulty = false;
+    Units.Percentage fouling(min = 0, max=100); // Fouling percentage
+
     // Initialization parameters
     parameter Units.MassFlowRate Q_cold_0 = 500;
     parameter Units.MassFlowRate Q_hot_0 = 50;
@@ -77,6 +81,11 @@ model Evaporator
   WaterSteam.Connectors.Outlet C_cold_out annotation (Placement(transformation(extent={{-40,60},{-20,80}}), iconTransformation(extent={{-40,60},{-20,80}})));
 
 equation
+  // Failure modes
+  if not faulty then
+    fouling = 0;
+  end if;
+
   // Definitions
   Q_cold = cold_side_heating.Q;
   Q_hot = hot_side_vaporising.Q;
@@ -115,7 +124,7 @@ equation
   cold_side_vaporising.h_out = x_steam_out * h_vap_sat + (1-x_steam_out)*h_liq_sat;
 
   HX_vaporising.W = W_vaporising;
-  HX_vaporising.Kth = Kth;
+  HX_vaporising.Kth = Kth*(1-fouling/100);
   HX_vaporising.S = S_vaporising;
   HX_vaporising.Q_cold = Q_cold;
   HX_vaporising.Q_hot = Q_hot;
diff --git a/MetroscopeModelingLibrary/MultiFluid/HeatExchangers/FuelHeater.mo b/MetroscopeModelingLibrary/MultiFluid/HeatExchangers/FuelHeater.mo
@@ -25,6 +25,10 @@ model FuelHeater
   Units.Temperature T_hot_in;
   Units.Temperature T_hot_out;
 
+  // Failure modes
+  parameter Boolean faulty = false;
+  Units.Percentage fouling(min = 0, max=100); // Fouling percentage
+
   // Initialization parameters
   parameter Units.MassFlowRate Q_cold_0 = 500;
   parameter Units.MassFlowRate Q_hot_0 = 50;
@@ -57,6 +61,12 @@ model FuelHeater
   Fuel.Pipes.Pipe cold_side_pipe annotation (Placement(transformation(extent={{-52,-10},{-32,10}})));
   Fuel.BaseClasses.IsoPFlowModel cold_side annotation (Placement(transformation(extent={{0,-10},{20,10}})));
 equation
+
+  // Failure modes
+  if not faulty then
+    fouling = 0;
+  end if;
+
     // Definitions
   Q_cold = cold_side.Q;
   Q_hot = hot_side.Q;
@@ -78,7 +88,7 @@ equation
   // Power Exchange
   HX.W = W;
   HX.S = S;
-  HX.Kth = Kth;
+  HX.Kth = Kth*(1-fouling/100);
   HX.Q_cold = Q_cold;
   HX.Q_hot = Q_hot;
   HX.T_cold_in = T_cold_in;
diff --git a/MetroscopeModelingLibrary/Partial/HeatExchangers/hrsg_monophasic_HX.mo b/MetroscopeModelingLibrary/Partial/HeatExchangers/hrsg_monophasic_HX.mo
@@ -23,6 +23,10 @@ partial model hrsg_monophasic_HX
   Units.Temperature T_cold_out;
   Units.Temperature T_hot_out;
 
+    // Failure modes
+  parameter Boolean faulty = false;
+  Units.Percentage fouling(min = 0, max=100); // Fouling percentage
+
     // Initialization parameters
   parameter Units.MassFlowRate Q_cold_0 = 500;
   parameter Units.MassFlowRate Q_hot_0 = 50;
@@ -69,6 +73,11 @@ protected
   MetroscopeModelingLibrary.Media.FlueGasesMedium.ThermodynamicState state_hot_out; // estimation of the flue gases outlet thermodynamic state
 
 equation
+  // Failure modes
+  if not faulty then
+    fouling = 0;
+  end if;
+
   // Definitions
   Q_cold =cold_side.Q;
   Q_hot =hot_side.Q;
@@ -89,7 +98,7 @@ equation
 
   // Power Exchange
   HX.W = W;
-  HX.Kth =  Kth;
+  HX.Kth =  Kth*(1-fouling/100);
   HX.S = S;
   HX.Q_cold = Q_cold;
   HX.Q_hot = Q_hot;
diff --git a/MetroscopeModelingLibrary/Partial/Pipes/Pipe.mo b/MetroscopeModelingLibrary/Partial/Pipes/Pipe.mo
@@ -9,8 +9,19 @@ partial model Pipe
   Inputs.InputDifferentialHeight delta_z(nominal=5) "Height difference between outlet and inlet";
   Units.DifferentialPressure DP_f "Singular pressure loss";
   Units.DifferentialPressure DP_z "Singular pressure loss";
+
+  // Failure modes
+  parameter Boolean faulty = false;
+  Real fouling; // Fouling coefficient
+
 equation
-  DP_f = -Kfr * Q^2 / rho;
+
+  // Failure modes
+  if not faulty then
+    fouling = 0;
+  end if;
+
+  DP_f = - (1+ fouling)  * Kfr * Q^2 / rho;
   DP_z = -rho * Constants.g * delta_z;
 
   DP = DP_f + DP_z;
diff --git a/MetroscopeModelingLibrary/Tests/FlueGases/Pipes/Filter_faulty.mo b/MetroscopeModelingLibrary/Tests/FlueGases/Pipes/Filter_faulty.mo
@@ -0,0 +1,15 @@
+within MetroscopeModelingLibrary.Tests.FlueGases.Pipes;
+model Filter_faulty
+  extends Filter_direct(filter(faulty=true));
+
+  Real Failure_fouling(start=0);
+
+equation
+
+  // Failure input
+  Failure_fouling = 0 + 1*time;
+
+  // Failure definition
+  filter.fouling = Failure_fouling;
+  annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram(coordinateSystem(preserveAspectRatio=false)));
+end Filter_faulty;
diff --git a/MetroscopeModelingLibrary/Tests/FlueGases/Pipes/package.order b/MetroscopeModelingLibrary/Tests/FlueGases/Pipes/package.order
@@ -2,3 +2,4 @@ Pipe_direct
 Pipe_reverse
 Filter_direct
 Filter_reverse
+Filter_faulty
diff --git a/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Economiser_direct.mo b/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Economiser_direct.mo
@@ -53,7 +53,7 @@ equation
 
   connect(economiser.C_cold_in, cold_source.C_out) annotation (Line(points={{1,7},{0,7},{0,30},{10,30},{10,39}}, color={28,108,200}));
   connect(economiser.C_hot_out, hot_sink.C_in) annotation (Line(points={{5,0},{63,0}}, color={95,95,95}));
-  connect(economiser.C_cold_out, cold_sink.C_in) annotation (Line(points={{-5,-7},{-5,-41},{-12,-41}}, color={28,108,200}));
+  connect(economiser.C_cold_out, cold_sink.C_in) annotation (Line(points={{-5,7},{-5,-41},{-12,-41}},  color={28,108,200}));
   connect(economiser.C_hot_in, hot_source.C_out) annotation (Line(points={{-9,0},{-49,0}}, color={95,95,95}));
   annotation (Icon(coordinateSystem(preserveAspectRatio=false), graphics={
         Ellipse(lineColor = {75,138,73},
diff --git a/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Economiser_faulty.mo b/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Economiser_faulty.mo
@@ -0,0 +1,16 @@
+within MetroscopeModelingLibrary.Tests.Multifluid.HeatExchangers;
+model Economiser_faulty
+  extends Economiser_direct(economiser(faulty=true));
+
+  Real Failure_fouling(start=0);
+
+equation
+
+  // Failure input
+  Failure_fouling = 0 + 10*time;
+
+  // Failure definition
+  economiser.fouling = Failure_fouling;
+
+  annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram(coordinateSystem(preserveAspectRatio=false)));
+end Economiser_faulty;
diff --git a/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Evaporator_faulty.mo b/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Evaporator_faulty.mo
@@ -0,0 +1,16 @@
+within MetroscopeModelingLibrary.Tests.Multifluid.HeatExchangers;
+model Evaporator_faulty
+   extends Evaporator_direct(evaporator(faulty=true));
+
+  Real Failure_fouling(start=0);
+
+equation
+
+  // Failure input
+  Failure_fouling = 0 + 10*time;
+
+  // Failure definition
+  evaporator.fouling = Failure_fouling;
+
+  annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram(coordinateSystem(preserveAspectRatio=false)));
+end Evaporator_faulty;
diff --git a/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/FuelHeater_direct.mo b/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/FuelHeater_direct.mo
@@ -48,7 +48,8 @@ equation
 
   connect(fuelHeater.C_cold_in, cold_source.C_out) annotation (Line(points={{-35,0},{-63,0}}, color={213,213,0}));
   connect(fuelHeater.C_cold_out, cold_sink.C_in) annotation (Line(points={{35,0},{65,0}}, color={213,213,0}));
-  connect(fuelHeater.C_hot_out, hot_sink.C_in) annotation (Line(points={{-15,33.6},{-14,33.6},{-14,57}}, color={28,108,200}));
-  connect(fuelHeater.C_hot_in, hot_source.C_out) annotation (Line(points={{15,33.6},{16,33.6},{16,57}}, color={28,108,200}));
+  connect(fuelHeater.C_hot_out, hot_sink.C_in) annotation (Line(points={{-20,-33.6},{-14,-33.6},{-14,57}},
+                                                                                                         color={28,108,200}));
+  connect(fuelHeater.C_hot_in, hot_source.C_out) annotation (Line(points={{20,33.6},{16,33.6},{16,57}}, color={28,108,200}));
   annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram(coordinateSystem(preserveAspectRatio=false)));
 end FuelHeater_direct;
diff --git a/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/FuelHeater_faulty.mo b/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/FuelHeater_faulty.mo
@@ -0,0 +1,15 @@
+within MetroscopeModelingLibrary.Tests.Multifluid.HeatExchangers;
+model FuelHeater_faulty
+  extends FuelHeater_direct(fuelHeater(faulty=true));
+
+  Real Failure_fouling(start=0);
+
+equation
+
+  // Failure input
+  Failure_fouling = 0 + 10*time;
+
+  // Failure definition
+  fuelHeater.fouling = Failure_fouling;
+  annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram(coordinateSystem(preserveAspectRatio=false)));
+end FuelHeater_faulty;
diff --git a/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Superheater_faulty.mo b/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/Superheater_faulty.mo
@@ -0,0 +1,15 @@
+within MetroscopeModelingLibrary.Tests.Multifluid.HeatExchangers;
+model Superheater_faulty
+    extends Superheater_direct(superheater(faulty=true));
+
+  Real Failure_fouling(start=0);
+
+equation
+
+  // Failure input
+  Failure_fouling = 0 + 10*time;
+
+  // Failure definition
+  superheater.fouling = Failure_fouling;
+  annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram(coordinateSystem(preserveAspectRatio=false)));
+end Superheater_faulty;
diff --git a/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/package.order b/MetroscopeModelingLibrary/Tests/Multifluid/HeatExchangers/package.order
@@ -1,9 +1,13 @@
 Economiser_direct
 Economiser_reverse
+Economiser_faulty
 Superheater_direct
 Superheater_reverse
+Superheater_faulty
 Evaporator_direct
 Evaporator_reverse
+Evaporator_faulty
 FuelHeater_direct
 FuelHeater_reverse
+FuelHeater_faulty
 HXmoistAirWater_reverse
diff --git a/MetroscopeModelingLibrary/Tests/WaterSteam/Machines/Pump_reverse.mo b/MetroscopeModelingLibrary/Tests/WaterSteam/Machines/Pump_reverse.mo
@@ -3,9 +3,9 @@ model Pump_reverse
   extends MetroscopeModelingLibrary.Icons.Tests.WaterSteamTestIcon;
 
   // Boundary conditions
-  input Units.Pressure source_P(start=2e5);
-  input Units.Temperature source_T(start=20 + 273.15);
-  input Units.NegativeMassFlowRate source_Q(start=-100);
+  input Units.Pressure source_P(start=2.3e5);
+  input Units.Temperature source_T(start=7.8 + 273.15);
+  input Units.NegativeMassFlowRate source_Q(start=-30);
 
   // Component parameters
   parameter Real pump_VRot = 1400;
@@ -22,8 +22,8 @@ model Pump_reverse
   output Real pump_b3;
 
   // Calibration inputs
-  input Units.Pressure pump_P_out(start=6e5);
-  input Units.Temperature pump_T_out(start=20 + 273.15);
+  input Units.Pressure pump_P_out(start=2.31e5);
+  input Units.Temperature pump_T_out(start=7.801+ 273.15);
 
   .MetroscopeModelingLibrary.WaterSteam.Machines.Pump pump annotation (Placement(transformation(extent={{-10,-10},{10,10}}, origin={-30,0})));
   .MetroscopeModelingLibrary.WaterSteam.BoundaryConditions.Source source annotation (Placement(transformation(extent={{-10,-10},{10,10}}, origin={-70,0})));
