Line data Source code
1 : // EnergyPlus, Copyright (c) 1996-2023, The Board of Trustees of the University of Illinois,
2 : // The Regents of the University of California, through Lawrence Berkeley National Laboratory
3 : // (subject to receipt of any required approvals from the U.S. Dept. of Energy), Oak Ridge
4 : // National Laboratory, managed by UT-Battelle, Alliance for Sustainable Energy, LLC, and other
5 : // contributors. All rights reserved.
6 : //
7 : // NOTICE: This Software was developed under funding from the U.S. Department of Energy and the
8 : // U.S. Government consequently retains certain rights. As such, the U.S. Government has been
9 : // granted for itself and others acting on its behalf a paid-up, nonexclusive, irrevocable,
10 : // worldwide license in the Software to reproduce, distribute copies to the public, prepare
11 : // derivative works, and perform publicly and display publicly, and to permit others to do so.
12 : //
13 : // Redistribution and use in source and binary forms, with or without modification, are permitted
14 : // provided that the following conditions are met:
15 : //
16 : // (1) Redistributions of source code must retain the above copyright notice, this list of
17 : // conditions and the following disclaimer.
18 : //
19 : // (2) Redistributions in binary form must reproduce the above copyright notice, this list of
20 : // conditions and the following disclaimer in the documentation and/or other materials
21 : // provided with the distribution.
22 : //
23 : // (3) Neither the name of the University of California, Lawrence Berkeley National Laboratory,
24 : // the University of Illinois, U.S. Dept. of Energy nor the names of its contributors may be
25 : // used to endorse or promote products derived from this software without specific prior
26 : // written permission.
27 : //
28 : // (4) Use of EnergyPlus(TM) Name. If Licensee (i) distributes the software in stand-alone form
29 : // without changes from the version obtained under this License, or (ii) Licensee makes a
30 : // reference solely to the software portion of its product, Licensee must refer to the
31 : // software as "EnergyPlus version X" software, where "X" is the version number Licensee
32 : // obtained under this License and may not use a different name for the software. Except as
33 : // specifically required in this Section (4), Licensee shall not use in a company name, a
34 : // product name, in advertising, publicity, or other promotional activities any name, trade
35 : // name, trademark, logo, or other designation of "EnergyPlus", "E+", "e+" or confusingly
36 : // similar designation, without the U.S. Department of Energy's prior written consent.
37 : //
38 : // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
39 : // IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY
40 : // AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
41 : // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
42 : // CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
43 : // SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
44 : // THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
45 : // OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
46 : // POSSIBILITY OF SUCH DAMAGE.
47 :
48 : #include "AirflowNetwork/Properties.hpp"
49 : #include <EnergyPlus/Data/EnergyPlusData.hh>
50 : #include <EnergyPlus/General.hh>
51 :
52 : namespace EnergyPlus {
53 :
54 : namespace AirflowNetwork {
55 :
56 894 : AirState::AirState(Real64 const density)
57 894 : : temperature(20.0), humidity_ratio(0.0), density(density), sqrt_density(sqrt(density)), viscosity(AIRDYNAMICVISCOSITY_CONSTEXPR(20.0))
58 : {
59 894 : }
60 :
61 0 : AirState::AirState()
62 0 : : temperature(20.0), humidity_ratio(0.0), density(AIRDENSITY_CONSTEXPR(101325.0, 20.0, 0.0)),
63 0 : sqrt_density(std::sqrt(AIRDENSITY_CONSTEXPR(101325.0, 20.0, 0.0))), viscosity(AIRDYNAMICVISCOSITY_CONSTEXPR(20.0))
64 : {
65 0 : }
66 :
67 171458880 : Real64 AirProperties::density(Real64 P, // Barometric pressure
68 : Real64 T, // Temperature in Celsius
69 : Real64 W // Humidity ratio
70 :
71 : )
72 : {
73 171458880 : return Psychrometrics::PsyRhoAirFnPbTdbW(m_state, P, T, W);
74 : }
75 :
76 18827897 : Real64 AirProperties::thermal_conductivity(Real64 T // Temperature in Celsius
77 : )
78 : {
79 : // Dry air thermal conductivity {W/m-K}
80 : // Correlated over the range -20C to 70C
81 : // Reference Cengel & Ghajar, Heat and Mass Transfer. 5th ed.
82 :
83 18827897 : Real64 const LowerLimit = -20;
84 18827897 : Real64 const UpperLimit = 70;
85 :
86 18827897 : Real64 const a = 0.02364;
87 18827897 : Real64 const b = 0.0000754772569209165;
88 18827897 : Real64 const c = -2.40977632412045e-8;
89 :
90 18827897 : if (T < LowerLimit) {
91 333204 : if (lowerLimitErrIdx == 0) {
92 2 : ShowWarningMessage(m_state, "Air temperature below lower limit of -20C for conductivity calculation");
93 : }
94 999612 : ShowRecurringWarningErrorAtEnd(m_state,
95 666408 : format("Air temperature below lower limit of -20C for conductivity calculation. Air temperature of {:.1R} "
96 : "used for conductivity calculation.",
97 333204 : LowerLimit),
98 : lowerLimitErrIdx);
99 333204 : T = LowerLimit;
100 18494693 : } else if (T > UpperLimit) {
101 87 : if (upperLimitErrIdx == 0) {
102 1 : ShowWarningMessage(m_state, "Air temperature above upper limit of 70C for conductivity calculation");
103 : }
104 261 : ShowRecurringWarningErrorAtEnd(m_state,
105 174 : format("Air temperature above upper limit of 70C for conductivity calculation. Air temperature of {:.1R} "
106 : "used for conductivity calculation.",
107 87 : UpperLimit),
108 : upperLimitErrIdx);
109 87 : T = UpperLimit;
110 : }
111 :
112 18827897 : return a + b * T + c * pow_2(T);
113 : }
114 :
115 15465009 : Real64 AirProperties::dynamic_viscosity(Real64 T // Temperature in Celsius
116 : )
117 : {
118 15465009 : return 1.71432e-5 + 4.828e-8 * T;
119 : }
120 :
121 15464926 : Real64 AirProperties::kinematic_viscosity(Real64 P, // Barometric pressure
122 : Real64 T, // Temperature in Celsius
123 : Real64 W // Humidity ratio
124 :
125 : )
126 : {
127 : // Dry air kinematic viscosity {m2/s}
128 : // Correlated over the range -20C to 70C
129 : // Reference Cengel & Ghajar, Heat and Mass Transfer. 5th ed.
130 :
131 15464926 : Real64 const LowerLimit = -20;
132 15464926 : Real64 const UpperLimit = 70;
133 :
134 15464926 : if (T < LowerLimit) {
135 646902 : T = LowerLimit;
136 14818024 : } else if (T > UpperLimit) {
137 0 : T = UpperLimit;
138 : }
139 :
140 15464926 : return dynamic_viscosity(T) / Psychrometrics::PsyRhoAirFnPbTdbW(m_state, P, T, W);
141 : }
142 :
143 7732463 : Real64 AirProperties::thermal_diffusivity(Real64 P, // Barometric pressure
144 : Real64 T, // Temperature in Celsius
145 : Real64 W // Humidity ratio
146 : )
147 : {
148 : // Dry air thermal diffusivity {-}
149 : // Correlated over the range -20C to 70C
150 : // Reference Cengel & Ghajar, Heat and Mass Transfer. 5th ed.
151 :
152 7732463 : Real64 const LowerLimit = -20;
153 7732463 : Real64 const UpperLimit = 70;
154 :
155 7732463 : if (T < LowerLimit) {
156 0 : T = LowerLimit;
157 7732463 : } else if (T > UpperLimit) {
158 0 : T = UpperLimit;
159 : }
160 :
161 7732463 : return thermal_conductivity(T) / (AIRCP(W) * Psychrometrics::PsyRhoAirFnPbTdbW(m_state, P, T, W));
162 : }
163 :
164 7732463 : Real64 AirProperties::prandtl_number(Real64 P, // Barometric pressure
165 : Real64 T, // Temperature in Celsius
166 : Real64 W // Humidity ratio
167 : )
168 : {
169 : // Dry air Prandtl number {-}
170 : // Correlated over the range -20C to 70C
171 : // Reference Cengel & Ghajar, Heat and Mass Transfer. 5th ed.
172 :
173 7732463 : Real64 const LowerLimit = -20;
174 7732463 : Real64 const UpperLimit = 70;
175 :
176 7732463 : if (T < LowerLimit) {
177 646902 : T = LowerLimit;
178 7085561 : } else if (T > UpperLimit) {
179 0 : T = UpperLimit;
180 : }
181 :
182 7732463 : return kinematic_viscosity(P, T, W) / thermal_diffusivity(P, T, W);
183 : }
184 :
185 : } // namespace AirflowNetwork
186 :
187 2313 : } // namespace EnergyPlus
|