Stack effect airflow - Project itz fixes

IDEAS/Airflow/Multizone/CrackOrOperableDoor.mo

@@ -0,0 +1,176 @@
+within IDEAS.Airflow.Multizone;
+
+model CrackOrOperableDoor
+  "Door model using discretization along height coordinate"
+  extends IDEAS.Fluid.Interfaces.PartialFourPortInterface(
+    redeclare final package Medium1 = Medium,
+    redeclare final package Medium2 = Medium,
+    final allowFlowReversal1=true,
+    final allowFlowReversal2=true,
+    final m1_flow_nominal=10/3600*1.2,
+    final m2_flow_nominal=m1_flow_nominal,
+    final m1_flow_small=1E-4*abs(m1_flow_nominal),
+    final m2_flow_small=1E-4*abs(m2_flow_nominal));
+  extends IDEAS.Airflow.Multizone.BaseClasses.ErrorControl(forceErrorControlOnFlow=true); //force error control on flow rates
+
+  replaceable package Medium =
+    Modelica.Media.Interfaces.PartialMedium "Medium in the component"
+      annotation (choices(
+        choice(redeclare package Medium = IDEAS.Media.Air "Moist air")));
+
+  parameter Modelica.Units.SI.Velocity vZer=0.001
+    "Minimum velocity to prevent zero flow. Recommended: 0.001";
+  parameter Modelica.Units.SI.Length wOpe=0.9 "Width of opening"
+    annotation (Dialog(group="Geometry"));
+  parameter Modelica.Units.SI.Length hOpe=2.1 "Height of opening"
+    annotation (Dialog(group="Geometry"));
+
+  parameter BoundaryConditions.Types.InterZonalAirFlow interZonalAirFlowType
+    "Interzonal air flow type";
+  parameter Modelica.Units.SI.PressureDifference dpCloRat(
+    min=0,
+    displayUnit="Pa") = 50 "Pressure drop at rating condition of closed door"
+    annotation (Dialog(group="Rating conditions"));
+
+  parameter Modelica.Units.SI.Length h_b1 "Height at port b1";
+  parameter Modelica.Units.SI.Length h_b2 = 0 "Height at port b2";
+  parameter Modelica.Units.SI.Length h_a1 = 0 "Height at port a1";
+  parameter Modelica.Units.SI.Length h_a2  "Height at port a2";
+
+  parameter Real CDCloRat(min=0, max=1)=1
+    "Discharge coefficient at rating conditions of closed door"
+      annotation (Dialog(group="Rating conditions"));
+   parameter Modelica.Units.SI.Area A_q50 "Surface area for leakage computation (closed door)";
+   parameter Real q50(unit="m3/(h.m2)") "Surface air tightness";
+
+  parameter Modelica.Units.SI.Area LClo(min=0) = q50*A_q50*1.2/3600 *(1.2/2/dpCloRat)^mClo/CDClo
+    "Effective leakage area of closed door"
+    annotation (Dialog(group="Closed door"));
+
+  parameter Real CDOpe=0.65 "Discharge coefficient of open door"
+    annotation (Dialog(group="Open door"));
+  parameter Real CDClo=0.65 "Discharge coefficient of closed door"
+    annotation (Dialog(group="Closed door"));
+
+  parameter Real mOpe = 0.5 "Flow exponent for door of open door"
+    annotation (Dialog(group="Open door"));
+  parameter Real mClo= 0.65 "Flow exponent for crack of closed door"
+    annotation (Dialog(group="Closed door"));
+
+  parameter Integer nCom=10 "Number of compartments for the discretization";
+
+  parameter Boolean useDoor = false  "=true, to use operable door instead of a crack";
+  parameter Boolean use_y = true "=true, to use control input";
+
+  parameter Modelica.Units.SI.PressureDifference dp_turbulent(
+    min=0,
+    displayUnit="Pa") = 0.01
+    "Pressure difference where laminar and turbulent flow relation coincide. Recommended: 0.01";
+
+  Modelica.Blocks.Interfaces.RealInput y(min=0, max=1, unit="1") if useDoor and use_y
+    "Opening signal, 0=closed, 1=open"
+    annotation (Placement(transformation(extent={{-120,-10},{-100,10}}), iconTransformation(extent={{-120,-10},{-100,10}})));
+ IDEAS.Airflow.Multizone.Point_m_flow point_m_flow1(
+  redeclare package Medium = Medium, 
+  dpMea_nominal = dpCloRat, 
+  forceErrorControlOnFlow = false, 
+  mMea_flow_nominal = (if interZonalAirFlowType == IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts  then LClo/2 else LClo)*CDClo*(2*dpCloRat/1.2)^mClo,
+  m = mClo, 
+  useDefaultProperties = false) if not useDoor "Pressure drop equation" annotation(
+    Placement(visible = true, transformation(origin = {0, 60}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+ IDEAS.Airflow.Multizone.ReversibleDensityColumn col_b1(redeclare package Medium = Medium, h=h_b1) if interZonalAirFlowType == IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts "Column for port b1" annotation(
+    Placement(visible = true, transformation(origin = {0, 70}, extent = {{50, -10}, {70, 10}}, rotation = 0)));
+ IDEAS.Airflow.Multizone.ReversibleDensityColumn col_a1(redeclare package Medium = Medium, h=h_a1) if interZonalAirFlowType == IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts "Column for port a1" annotation(
+    Placement(visible = true, transformation(origin = {0, 70}, extent = {{-70, -10}, {-50, 10}}, rotation = 0)));
+ IDEAS.Airflow.Multizone.ReversibleDensityColumn col_b2(redeclare package Medium = Medium, h=h_b2) if interZonalAirFlowType == IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts "Column for port b2" annotation(
+    Placement(visible = true, transformation(origin = {0, -50}, extent = {{-70, -10}, {-50, 10}}, rotation = 0)));
+ IDEAS.Airflow.Multizone.ReversibleDensityColumn col_a2(redeclare package Medium = Medium, h=h_a2) if interZonalAirFlowType == IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts "Column for port a2" annotation(
+    Placement(visible = true, transformation(origin = {0, -50}, extent = {{50, -10}, {70, 10}}, rotation = 0)));
+ IDEAS.Airflow.Multizone.Point_m_flow point_m_flow2(
+   redeclare package Medium = Medium, 
+   dpMea_nominal = dpCloRat, 
+   forceErrorControlOnFlow = false, 
+   m = mClo,
+   mMea_flow_nominal = LClo/2*CDClo*(2*dpCloRat/1.2)^mClo,
+   useDefaultProperties = false) if not useDoor and interZonalAirFlowType == IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts "Pressure drop equation" annotation(
+    Placement(visible = true, transformation(origin = {0, -60}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+ IDEAS.Airflow.Multizone.DoorDiscretizedOperable doo(
+   redeclare package Medium = Medium,
+   final hA = 0, 
+   final hB = 0,
+   final forceErrorControlOnFlow=false,
+   dp_turbulent=dp_turbulent,
+   nCom=nCom,
+   CDOpe=CDOpe,
+   CDClo=CDClo,
+   mOpe=mOpe,
+   mClo=mClo,
+   CDCloRat=CDCloRat,
+   wOpe=wOpe,
+   hOpe=hOpe,
+   dpCloRat=dpCloRat,
+   LClo=LClo,
+   vZer=vZer) if useDoor and interZonalAirFlowType == IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts annotation(
+    Placement(visible = true, transformation(origin = {0, 0}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+ IDEAS.Fluid.Sources.Boundary_pT bou(
+   redeclare package Medium = Medium, 
+   nPorts = 2) 
+   if interZonalAirFlowType <> IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts 
+   "Sets absolute pressure when the ports are not connected externally" annotation(
+    Placement(visible = true, transformation(origin = {0, -90}, extent = {{-10, 10}, {10, -10}}, rotation = 90)));
+ Modelica.Blocks.Sources.Constant constOne(final k=1)
+   if not use_y
+   "Door constantly opened" annotation(
+    Placement(visible = true, transformation(origin = {-54, -14}, extent = {{-6, -6}, {6, 6}}, rotation = 0)));
+initial equation
+  assert( not (interZonalAirFlowType <> IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts and useDoor), 
+  "In " +getInstanceName() + ": Cannot use a controllable door unless interZonalAirFlowType == TwoPorts.");
+
+equation
+  connect(col_a1.port_a, point_m_flow1.port_a) annotation(
+    Line(points = {{-60, 80}, {-60, 84}, {-20, 84}, {-20, 60}, {-10, 60}}, color = {0, 127, 255}));
+  connect(col_b1.port_a, point_m_flow1.port_b) annotation(
+    Line(points = {{60, 80}, {60, 84}, {20, 84}, {20, 60}, {10, 60}}, color = {0, 127, 255}));
+  connect(col_b2.port_a, point_m_flow2.port_a) annotation(
+    Line(points = {{-60, -40}, {-60, -36}, {-20, -36}, {-20, -60}, {-10, -60}}, color = {0, 127, 255}));
+  connect(col_a2.port_a, point_m_flow2.port_b) annotation(
+    Line(points = {{60, -40}, {60, -36}, {20, -36}, {20, -60}, {10, -60}}, color = {0, 127, 255}));
+ connect(col_b2.port_b, port_b2) annotation(
+    Line(points = {{-60, -60}, {-100, -60}}, color = {0, 127, 255}));
+ connect(col_a2.port_b, port_a2) annotation(
+    Line(points = {{60, -60}, {100, -60}}, color = {0, 127, 255}));
+ connect(col_b1.port_b, port_b1) annotation(
+    Line(points = {{60, 60}, {100, 60}}, color = {0, 127, 255}));
+ connect(col_a1.port_b, port_a1) annotation(
+    Line(points = {{-60, 60}, {-100, 60}}, color = {0, 127, 255}));
+ connect(y, doo.y) annotation(
+    Line(points = {{-110, 0}, {-10, 0}}, color = {0, 0, 127}));
+ connect(doo.port_a1, col_a1.port_a) annotation(
+    Line(points = {{-10, 6}, {-20, 6}, {-20, 84}, {-60, 84}, {-60, 80}}, color = {0, 127, 255}));
+ connect(doo.port_b1, col_b1.port_a) annotation(
+    Line(points = {{10, 6}, {20, 6}, {20, 84}, {60, 84}, {60, 80}}, color = {0, 127, 255}));
+ connect(doo.port_b2, col_b2.port_a) annotation(
+    Line(points = {{-10, -6}, {-20, -6}, {-20, -36}, {-60, -36}, {-60, -40}}, color = {0, 127, 255}));
+ connect(doo.port_a2, col_a2.port_a) annotation(
+    Line(points = {{10, -6}, {20, -6}, {20, -36}, {60, -36}, {60, -40}}, color = {0, 127, 255}));
+ connect(bou.ports[1], port_a2) annotation(
+    Line(points = {{0, -80}, {100, -80}, {100, -60}}, color = {0, 127, 255}));
+ if  interZonalAirFlowType <> IDEAS.BoundaryConditions.Types.InterZonalAirFlow.TwoPorts then
+    connect(point_m_flow1.port_a, port_a1) annotation(
+    Line(points = {{-10, 60}, {-100, 60}}, color = {0, 127, 255}));
+    connect(point_m_flow1.port_b, port_b1) annotation(
+    Line(points = {{10, 60}, {100, 60}}, color = {0, 127, 255}));
+    connect(doo.port_b1, port_b1) annotation(
+      Line(points = {{10, 6}, {20, 6}, {20, 60}, {100, 60}}, color = {0, 127, 255}));
+    connect(doo.port_a1, port_a1) annotation(
+      Line(points = {{-10, 6}, {-20, 6}, {-20, 60}, {-100, 60}}, color = {0, 127, 255}));
+ end if;
+ connect(constOne.y, doo.y) annotation(
+    Line(points = {{-48, -14}, {-32, -14}, {-32, 0}, {-10, 0}}, color = {0, 0, 127}));
+ connect(bou.ports[2], port_b2) annotation(
+    Line(points = {{0, -80}, {-100, -80}, {-100, -60}}, color = {0, 127, 255}));
+
+annotation(
+    Diagram,
+ Icon(graphics = {Polygon(lineColor = {0, 0, 255}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, points = {{-30, -10}, {-16, -8}, {-16, -14}, {-30, -16}, {-30, -10}}), Line(points = {{-54, 48}, {-36, 48}}), Line(points = {{-54, 20}, {-36, 20}}), Rectangle(fillColor = {255, 255, 255}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-100, 100}, {100, -100}}), Line(points = {{-54, -58}, {-36, -58}}), Rectangle(lineColor = {0, 0, 255}, fillColor = {255, 128, 0}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-46, -16}, {-20, -20}}), Rectangle(lineColor = {0, 0, 255}, fillColor = {85, 75, 55}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-60, 80}, {60, -84}}), Rectangle(fillColor = {215, 215, 215}, fillPattern = FillPattern.Solid, extent = {{-54, 72}, {56, -84}}), Polygon(fillColor = {95, 95, 95}, fillPattern = FillPattern.Solid, points = {{56, 72}, {-36, 66}, {-36, -90}, {56, -84}, {56, 72}}), Rectangle(lineColor = {0, 0, 255}, fillColor = {255, 128, 0}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-100, 2}, {-46, -2}}), Polygon(visible = false, origin = {75, 50}, rotation = 360, lineColor = {0, 128, 255}, fillColor = {0, 128, 255}, fillPattern = FillPattern.Solid, points = {{-5, 10}, {25, 10}, {-5, -10}, {-5, 10}}), Text(textColor = {0, 0, 127}, extent = {{-118, 34}, {-98, 16}}, textString = "y"), Line(points = {{-54, -6}, {-36, -6}}), Line(points = {{-54, -32}, {-36, -32}}), Polygon(visible = false, origin = {-79, -50}, rotation = 360, lineColor = {0, 128, 255}, fillColor = {0, 128, 255}, fillPattern = FillPattern.Solid, points = {{10, 10}, {-20, -10}, {10, -10}, {10, 10}}), Rectangle(lineColor = {0, 0, 255}, fillColor = {255, 128, 0}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-46, 2}, {-40, -16}})}, coordinateSystem(extent = {{-100, -100}, {100, 100}})));
+end CrackOrOperableDoor;

IDEAS/Airflow/Multizone/Examples/TrickleVentIDEAS.mo

@@ -0,0 +1,53 @@
+within IDEAS.Airflow.Multizone.Examples;
+model TrickleVentIDEAS
+  "Model with a trickle vent modelled using the models with flow based on tabulated data"
+  extends IDEAS.Airflow.Multizone.Examples.TrickleVent(west(nPorts=2), east(nPorts=2));
+
+  IDEAS.Airflow.Multizone.TrickleVent vent(
+    redeclare package Medium = Medium,
+    dp_nominal = 10, 
+    m_flow_nominal = 0.02614, 
+    use_y = true) 
+    "Analytic trickle vent implementation"  annotation(
+    Placement(visible = true, transformation(origin = {20, -70}, extent = {{-10, 10}, {10, -10}}, rotation = 0)));
+  Modelica.Blocks.Sources.Ramp ramp(
+    duration = 1e6, 
+    height = -1, 
+    offset = 1, 
+    startTime = 1592000)  
+    "Step control signal" annotation(
+    Placement(visible = true, transformation(origin = {-90, -88}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+equation
+   connect(vent.port_a, east.ports[2]) annotation(
+    Line(points = {{10, -70}, {-25, -70}, {-25, -30}, {-30, -30}}, color = {0, 127, 255}));
+  connect(vent.port_b, west.ports[2]) annotation(
+    Line(points = {{30, -70}, {65, -70}, {65, -30}, {70, -30}}, color = {0, 127, 255}));
+  connect(ramp.y, vent.y) annotation(
+    Line(points = {{-79, -88}, {20.5, -88}, {20.5, -82}, {20, -82}}, color = {0, 0, 127}));
+  annotation (__Dymola_Commands(file="modelica://IDEAS/Resources/Scripts/Dymola/Airflow/Multizone/Examples/TrickleVentIDEAS.mos"
+        "Simulate and plot"),
+        experiment(
+      StopTime=2592000,
+      Interval=600,
+      Tolerance=1e-06),
+    Documentation(info="<html>
+<p>
+This model illustrates the use of the models
+<a href=\"modelica://IDEAS.Airflow.Multizone.TrickleVent\">
+IDEAS.Airflow.Multizone.Table_V_flow</a>, 
+which is an analytic alternative to the table implementation of
+<a href=\"modelica://IDEAS.Airflow.Multizone.Table_m_flow\">
+IDEAS.Airflow.Multizone.Table_m_flow</a> for
+
+modelling self regulating inlet vents.
+</p>
+</html>", revisions="<html>
+<ul>
+<li>
+July 9, 2023 by Filip Jorissen:<br/>
+First implementation.
+</li>
+</ul>
+</html>"),
+ Diagram);
+end TrickleVentIDEAS;

IDEAS/Airflow/Multizone/Examples/package.order

@@ -12,4 +12,5 @@ PressurizationData
 ReverseBuoyancy
 ReverseBuoyancy3Zones
 TrickleVent
+TrickleVentIDEAS
 ZonalFlow

IDEAS/Airflow/Multizone/ReversibleDensityColumn.mo

@@ -0,0 +1,42 @@
+within IDEAS.Airflow.Multizone;
+model ReversibleDensityColumn
+  "Vertical shaft with no friction and no storage of heat and mass, reversible because it can handle negative column heights"
+
+  extends IDEAS.Airflow.Multizone.MediumColumn(
+    h(min=-Modelica.Constants.inf), 
+    final densitySelection = IDEAS.Airflow.Multizone.Types.densitySelection.fromBottom);
+    // by convention, port_b must be connected to a zone instead of a flow element
+    // h is allowed to be negative to accomodate for this convention
+
+  annotation (
+Icon(graphics={
+    Line(
+      points={{0,100},{0,-100},{0,-98}}),
+    Text(origin = {-126, 2},lineColor = {0, 0, 127}, extent = {{24, -78}, {106, -100}}, textString = "Zone/Amb"),
+    Text(origin = {-130, 4}, lineColor = {0, 0, 127}, extent = {{32, 104}, {98, 70}}, textString = "FlowElem"),
+    Text(lineColor = {0, 0, 127}, extent = {{36, 26}, {88, -10}}, textString = "h=%h"),
+    Rectangle(fillColor = {255, 0, 0}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-16, 80}, {16, -80}}),
+    Rectangle(visible = false, fillColor = {85, 170, 255}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-16, 80}, {16, 0}}),
+    Rectangle(visible = false, fillColor = {85, 170, 255}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-16, 80}, {16, 54}}),
+    Rectangle(fillColor = {85, 170, 255}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-16, 0}, {16, -82}}),
+    Rectangle(visible = false, fillColor = {85, 170, 255}, pattern = LinePattern.None, fillPattern = FillPattern.Solid, extent = {{-16, -55}, {16, -80}})}),
+  defaultComponentName="col",
+  Documentation(info="<html>
+<p>
+This model describes the pressure difference of a vertical medium
+column. It can be used to model the pressure difference caused by
+stack effect.
+It is a variation on IDEAS.Airflow.Multizone.MediumColumn.
+</p>
+</html>",
+      revisions="<html>
+<ul>
+<li>
+January 19, 2022, by Klaas De Jonge:<br/>
+Adapted IDEAS.Airflow.Multizone.MediumColumn to obtain the current model where input of
+ h can be negative and cleaned out the model as the density should always be set by port_b. 
+This makes port_a not nececarilly always the top port.
+</li>
+</ul>
+</html>"));
+end ReversibleDensityColumn;

IDEAS/Airflow/Multizone/TrickleVent.mo

@@ -5,15 +5,21 @@ model TrickleVent "Self regulating trickle vent"
     final from_dp=true,
     final homotopyInitialization=true);
 
+  parameter Real l(min=1e-10, max=1) = 0.0001
+    "Valve leakage, l=Kv(y=0)/Kv(y=1)";
   parameter Real deltaM(min=1E-6) = 0.3
     "Fraction of nominal mass flow rate where transition to turbulent occurs"
        annotation(Evaluate=true,
                   Dialog(group = "Transition to laminar",
                          enable = not linearized));
-
+  parameter Boolean use_y = false 
+    "=true, to enable control input"
+    annotation(Evaluate=true, Dialog(group="Control"));
   final parameter Real k = if computeFlowResistance then
         m_flow_nominal_pos / sqrt(dp_nominal_pos) else 0
     "Flow coefficient, k=m_flow/sqrt(dp), with unit=(kg.m)^(1/2)";
+  Modelica.Blocks.Interfaces.RealInput y(min=0,max=1) if use_y "Control input for trickle vent Kv value" annotation(
+    Placement(visible = true, transformation(origin = {0, 120}, extent = {{-20, -20}, {20, 20}}, rotation = -90), iconTransformation(origin = {0, 120}, extent = {{-20, -20}, {20, 20}}, rotation = -90)));
 protected
   final parameter Boolean computeFlowResistance=(dp_nominal_pos > Modelica.Constants.eps)
     "Flag to enable/disable computation of flow resistance"
@@ -23,13 +29,18 @@ protected
     then if from_dp then k^2/m_flow_nominal_pos else m_flow_nominal_pos/k^2
     else 0
     "Precomputed coefficient to avoid division by parameter";
+    Modelica.Blocks.Interfaces.RealInput y_internal "Internal variable for conditional variables";
 initial equation
  if computeFlowResistance then
    assert(m_flow_turbulent > 0, "m_flow_turbulent must be bigger than zero.");
  end if;
 
  assert(m_flow_nominal_pos > 0, "m_flow_nominal_pos must be non-zero. Check parameters.");
 equation
+  connect(y, y_internal);
+  if not use_y then
+    y_internal = 1;
+  end if;
   // Pressure drop calculation
   if computeFlowResistance then
     if linearized then
@@ -40,7 +51,7 @@ equation
             m_flow_nominal,
             IDEAS.Fluid.BaseClasses.FlowModels.basicFlowFunction_dp(
               dp=dp,
-              k=k,
+              k=k*(l + (1-l)*y_internal),
               m_flow_turbulent=m_flow_turbulent),
             m_flow_nominal/10);
     end if; // linearized
@@ -56,12 +67,20 @@ The positive mass flow rate is limited to <code>m_flow_nominal</code>
 at a pressure difference of <code>dp_nominal</code>.
 For negative pressure differences the mass flow rate is not limited.
 </p>
+<p>
+An optional control input can be enabled, which reduces the flow rate for the same
+pressure difference, but which does not affect the maximum flow rate.
+</p>
 </html>", revisions="<html>
 <ul>
 <li>
+July 9, 2023, by Filip Jorissen:<br/>
+Add control input.
+</li>
+<li>
 September 21, 2021, by Filip Jorissen:<br/>
 First implementation.
 </li>
 </ul>
 </html>"));
-end TrickleVent;
+end TrickleVent;

IDEAS/Airflow/Multizone/package.order

@@ -11,6 +11,7 @@ MediumColumnDynamic
 Orifice
 Point_m_flow
 Points_m_flow
+ReversibleDensityColumn
 Table_V_flow
 Table_m_flow
 ZonalFlow_ACS
@@ -20,3 +21,4 @@ Examples
 Validation
 BaseClasses
 TrickleVent
+CrackOrOperableDoor