Line data Source code
1 : // EnergyPlus, Copyright (c) 1996-2025, 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 315 : AirState::AirState(Real64 const density)
57 315 : : temperature(20.0), humidity_ratio(0.0), density(density), sqrt_density(sqrt(density)), viscosity(AIRDYNAMICVISCOSITY_CONSTEXPR(20.0))
58 : {
59 315 : }
60 :
61 8 : AirState::AirState()
62 8 : : temperature(20.0), humidity_ratio(0.0), density(AIRDENSITY_CONSTEXPR(101325.0, 20.0, 0.0)),
63 8 : sqrt_density(std::sqrt(AIRDENSITY_CONSTEXPR(101325.0, 20.0, 0.0))), viscosity(AIRDYNAMICVISCOSITY_CONSTEXPR(20.0))
64 : {
65 8 : }
66 :
67 5299181 : Real64 AirProperties::density(Real64 P, // Barometric pressure
68 : Real64 T, // Temperature in Celsius
69 : Real64 W // Humidity ratio
70 :
71 : )
72 : {
73 5299181 : return Psychrometrics::PsyRhoAirFnPbTdbW(m_state, P, T, W);
74 : }
75 :
76 798713 : 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 798713 : Real64 const LowerLimit = -20;
84 798713 : Real64 const UpperLimit = 70;
85 :
86 798713 : Real64 const a = 0.02364;
87 798713 : Real64 const b = 0.0000754772569209165;
88 798713 : Real64 const c = -2.40977632412045e-8;
89 :
90 798713 : if (T < LowerLimit) {
91 1 : if (lowerLimitErrIdx == 0) {
92 3 : ShowWarningMessage(m_state, "Air temperature below lower limit of -20C for conductivity calculation");
93 : }
94 7 : ShowRecurringWarningErrorAtEnd(m_state,
95 2 : format("Air temperature below lower limit of -20C for conductivity calculation. Air temperature of {:.1R} "
96 : "used for conductivity calculation.",
97 : LowerLimit),
98 1 : lowerLimitErrIdx);
99 1 : T = LowerLimit;
100 798712 : } else if (T > UpperLimit) {
101 1 : if (upperLimitErrIdx == 0) {
102 3 : ShowWarningMessage(m_state, "Air temperature above upper limit of 70C for conductivity calculation");
103 : }
104 7 : ShowRecurringWarningErrorAtEnd(m_state,
105 2 : format("Air temperature above upper limit of 70C for conductivity calculation. Air temperature of {:.1R} "
106 : "used for conductivity calculation.",
107 : UpperLimit),
108 1 : upperLimitErrIdx);
109 1 : T = UpperLimit;
110 : }
111 :
112 1597426 : return a + b * T + c * pow_2(T);
113 : }
114 :
115 798748 : Real64 AirProperties::dynamic_viscosity(Real64 T // Temperature in Celsius
116 : )
117 : {
118 798748 : return 1.71432e-5 + 4.828e-8 * T;
119 : }
120 :
121 798702 : 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 798702 : Real64 const LowerLimit = -20;
132 798702 : Real64 const UpperLimit = 70;
133 :
134 798702 : if (T < LowerLimit) {
135 1 : T = LowerLimit;
136 798701 : } else if (T > UpperLimit) {
137 1 : T = UpperLimit;
138 : }
139 :
140 798702 : return dynamic_viscosity(T) / Psychrometrics::PsyRhoAirFnPbTdbW(m_state, P, T, W);
141 : }
142 :
143 399359 : 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 399359 : Real64 const LowerLimit = -20;
153 399359 : Real64 const UpperLimit = 70;
154 :
155 399359 : if (T < LowerLimit) {
156 1 : T = LowerLimit;
157 399358 : } else if (T > UpperLimit) {
158 1 : T = UpperLimit;
159 : }
160 :
161 399359 : return thermal_conductivity(T) / (AIRCP(W) * Psychrometrics::PsyRhoAirFnPbTdbW(m_state, P, T, W));
162 : }
163 :
164 399351 : 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 399351 : Real64 const LowerLimit = -20;
174 399351 : Real64 const UpperLimit = 70;
175 :
176 399351 : if (T < LowerLimit) {
177 1 : T = LowerLimit;
178 399350 : } else if (T > UpperLimit) {
179 1 : T = UpperLimit;
180 : }
181 :
182 399351 : return kinematic_viscosity(P, T, W) / thermal_diffusivity(P, T, W);
183 : }
184 :
185 : } // namespace AirflowNetwork
186 :
187 : } // namespace EnergyPlus
|
