Add moisture source and transport mechanisms to HAMT

idd/Energy+.idd.in

@@ -8903,6 +8903,64 @@ ConstructionProperty:InternalHeatSource,
        \note 0.0 means in line with the tubing, 1.0 means at the midpoint between two adjacent pipes
        \note this field is ignored for 1-D calculations
 
+ConstructionProperty:InternalMoistureSource,
+      \memo HeatBalanceAlgorithm = CombinedHeatAndMoistureFiniteElement solution algorithm only.
+      \memo Adds moisture source to selected layer inside the component.
+      \memo Has no effect with other HeatBalanceAlgorithm solution algorithms
+  A1, \field Name
+      \required-field
+      \type alpha
+      \reference InternalMoistureSourceNames
+  A2, \field Construction Name
+      \type object-list
+      \object-list ConstructionNames
+  N1, \field Moisture Source Present In Layer Number
+      \required-field
+      \type integer
+      \minimum 1
+      \maximum 10
+      \note Layer number to which the moisture source should be applied (from outside to inside)
+      \note refers to the list of layers specified for the construction applied to the above surface
+      \note If a source occurs on a layer surface, the layer should be split into two layers (thin where source occurs) and the source added to the thin layer
+  A3, \field Source Type
+      \required-field
+      \type choice
+      \key UserDefined
+      \key StackAndOverPressure
+      \key AirflowNetwork
+      \note Type of Moisture Source Calculation
+      \note UserDefined: User defined input of air flow through component
+      \note StackAndOverPressure: Calculation of dynamic air flow through component due to stack effect and ventilation overpressure (according to Kuenzel, Zirkelbach and Schafazcek 2012)
+      \note AirflowNetwork: Calculation of dynamic air flow through component with AirflowNetwork model
+  N2, \field Air Flow Rate
+      \units m/s
+      \type real
+      \minimum -0.001
+      \maximum 0.001
+      \note Air flux density in m3 per m2 and second.
+      \note Only required for Source Type UserDefined.
+  N3, \field Stack Height
+      \units m
+      \type real
+      \minimum 0
+      \maximum 100
+      \note Height of the connected airspace in the building envelope element
+      \note Only required for Source Type StackAndOverPressure.
+  N4, \field Component Air Permeance
+      \units m/s
+      \type real
+      \minimum 0
+      \maximum 1
+      \note Moisture specific air permeance of the component in m3 per m2, hour and Pascal.
+      \note Only required for Source Type StackAndOverPressure.
+  N5; \field Mechanical Ventilation Overpressure
+      \units Pa
+      \type real
+      \minimum 0
+      \maximum 100
+      \note Constant mechanical ventilation overpressure in Pascal.
+      \note Only required for Source Type StackAndOverPressure.
+
 Construction:AirBoundary,
        \memo Indicates an open boundary between two zones. It may be used for base surfaces and fenestration surfaces.
        \memo The two adjacent zones are grouped together for solar, daylighting and radiant exchange.

src/EnergyPlus/HeatBalanceHAMTManager.hh

@@ -81,7 +81,10 @@ namespace HeatBalanceHAMTManager {
         // Members
         int matid;        // Material Id Number
         int sid;          // Surface Id Number
+        int layer_id;     // Layer Id Number
+        int source_id;    // Internal moisture source ID Number
         Real64 Qadds;     // Additional sources of heat
+        Real64 Wadds;     // Additional sources of water
         Real64 density;   // Density
         Real64 wthermalc; // Moisture Dependant Thermal Conductivity
         Real64 spech;     // Specific Heat capacity
@@ -112,14 +115,65 @@ namespace HeatBalanceHAMTManager {
 
         // Default Constructor
         subcell()
-            : matid(-1), sid(-1), Qadds(0.0), density(-1.0), wthermalc(0.0), spech(0.0), htc(-1.0), vtc(-1.0), mu(-1.0), volume(0.0), temp(0.0),
-              tempp1(0.0), tempp2(0.0), wreport(0.0), water(0.0), vp(0.0), vpp1(0.0), vpsat(0.0), rh(0.1), rhp1(0.1), rhp2(0.1), rhp(10.0),
-              dwdphi(-1.0), dw(-1.0), origin(3, 0.0), length(3, 0.0), overlap(6, 0.0), dist(6, 0.0), adjs(6, 0), adjsl(6, 0)
+            : matid(-1), sid(-1), layer_id(-1), source_id(-1), Qadds(0.0), Wadds(0.0), density(-1.0), wthermalc(0.0), spech(0.0), htc(-1.0),
+              vtc(-1.0), mu(-1.0), volume(0.0), temp(0.0), tempp1(0.0), tempp2(0.0), wreport(0.0), water(0.0), vp(0.0), vpp1(0.0), vpsat(0.0),
+              rh(0.1), rhp1(0.1), rhp2(0.1), rhp(10.0), dwdphi(-1.0), dw(-1.0), origin(3, 0.0), length(3, 0.0), overlap(6, 0.0), dist(6, 0.0),
+              adjs(6, 0), adjsl(6, 0)
         {
         }
     };
 
-    void ManageHeatBalHAMT(EnergyPlusData &state, int const SurfNum, Real64 &SurfTempInTmp, Real64 &TempSurfOutTmp);
+    // structure internal moisture source
+    struct InternalMoistureSource
+    {
+        int id;         // Internal moisture source id
+        int surface_id; // Surface id
+        int layer_id;   // Layer id
+
+        // int SourceType; // (1 = air flow input; 2 = infiltration model (Kuenzel); 3 = multi zone air flow model)
+        enum class Type
+        {
+            UserDefined = 1,
+            StackAndOverPressure,
+            AirflowNetwork
+        };
+        Type type;
+
+        // SourceType == 1: Air Flow input
+        Real64 moist_airflow_input; // moist air flow input for source type 1 in m3 / m2 s
+
+        // Source Type == 2: Air infiltration model
+        Real64 stack_height;                        // Stack height in m to create pressure difference for source type 2
+        Real64 component_air_permeance;             // Component air permeance in
+        Real64 mechanical_ventilation_overpressure; // Additional mechanical ventilation overpressure in Pa
+        Real64 delta_pressure;                      // Pressure difference over component in Pa
+
+        // Source Type == 3: Air flow network
+        int afn_id; // Id of the component element in the airflow network
+
+        Real64 moist_airflow;
+
+        InternalMoistureSource(int id = -1,
+                               int surface_id = -1,
+                               int layer_id = -1,
+                               Type type = Type::UserDefined,
+                               Real64 moist_airflow_input = 0.0,
+                               Real64 stack_height = 0.0,
+                               Real64 component_air_permeance = 0.0,
+                               Real64 mechanical_ventilation_overpressure = 0.0,
+                               Real64 delta_pressure = 0.0,
+                               int afn_id = 0,
+                               Real64 moist_airflow = 0.0)
+            : id(id), surface_id(surface_id), layer_id(layer_id), type(type), moist_airflow_input(moist_airflow_input), stack_height(stack_height),
+              component_air_permeance(component_air_permeance), mechanical_ventilation_overpressure(mechanical_ventilation_overpressure),
+              delta_pressure(delta_pressure), afn_id(afn_id), moist_airflow(moist_airflow)
+        {
+        }
+    };
+
+    // Functions
+
+    void ManageHeatBalHAMT(EnergyPlusData &state, int const SurfNum, Real64 &TempSurfInTmp, Real64 &TempSurfOutTmp);
 
     void GetHeatBalHAMTInput(EnergyPlusData &state);
 
@@ -136,6 +190,8 @@ namespace HeatBalanceHAMTManager {
 
     Real64 WVDC(Real64 const Temperature, Real64 const ambp);
 
+    Real64 SatAbsHum(EnergyPlusData &state, Real64 const Temperature);
+
     //                                 COPYRIGHT NOTICE
 
     //     Portions Copyright (c) University College London 2007.  All rights
@@ -186,6 +242,7 @@ struct HeatBalHAMTMgrData : BaseGlobalStruct
     bool latswitch = false;   // latent heat switch,
     bool rainswitch = false;  // rain switch,
     Array1D<HeatBalanceHAMTManager::subcell> cells;
+    Array1D<HeatBalanceHAMTManager::InternalMoistureSource> sources;
     bool OneTimeFlag = true;
     int qvpErrCount = 0;
     int qvpErrReport = 0;

src/EnergyPlus/RoomAirModelAirflowNetwork.cc

@@ -52,6 +52,7 @@
 
 // EnergyPlus Headers
 #include <AirflowNetwork/Elements.hpp>
+#include <EnergyPlus/AirflowNetworkBalanceManager.hh>
 #include <EnergyPlus/BaseboardElectric.hh>
 #include <EnergyPlus/BaseboardRadiator.hh>
 #include <EnergyPlus/Data/EnergyPlusData.hh>