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 : #ifndef HeatBalanceIntRadExchange_hh_INCLUDED
49 : #define HeatBalanceIntRadExchange_hh_INCLUDED
50 :
51 : // ObjexxFCL Headers
52 : #include <ObjexxFCL/Array1S.hh>
53 : #include <ObjexxFCL/Array2A.hh>
54 : #include <ObjexxFCL/Array2S.hh>
55 : #include <ObjexxFCL/Optional.hh>
56 :
57 : // EnergyPlus Headers
58 : #include <EnergyPlus/Data/BaseData.hh>
59 : #include <EnergyPlus/EnergyPlus.hh>
60 :
61 : namespace EnergyPlus {
62 :
63 : // Forward declarations
64 : struct EnergyPlusData;
65 :
66 : #define EP_HBIRE_SEQ
67 :
68 : namespace HeatBalanceIntRadExchange {
69 :
70 : void CalcInteriorRadExchange(EnergyPlusData &state,
71 : Array1S<Real64> const SurfaceTemp, // Current surface temperatures
72 : int const SurfIterations, // Number of iterations in calling subroutine
73 : Array1D<Real64> &NetLWRadToSurf, // Net long wavelength radiant exchange from other surfaces
74 : ObjexxFCL::Optional_int_const ZoneToResimulate = _, // if passed in, then only calculate for this zone
75 : std::string_view CalledFrom = "");
76 :
77 : void UpdateMovableInsulationFlag(EnergyPlusData &state,
78 : bool &MovableInsulationChange, // set to true if there is a change in the movable insulation state
79 : int const SurfNum // surface number of surface being investigated
80 : );
81 :
82 : void InitInteriorRadExchange(EnergyPlusData &state);
83 :
84 : void InitSolarViewFactors(EnergyPlusData &state);
85 :
86 : void AlignInputViewFactors(EnergyPlusData &state,
87 : std::string const &cCurrentModuleObject, // Object type
88 : bool &ErrorsFound // True when errors are found
89 : );
90 :
91 : void GetInputViewFactors(EnergyPlusData &state,
92 : std::string const &EnclosureName, // Needed to check for user input view factors.
93 : int const N, // NUMBER OF SURFACES
94 : Array2A<Real64> F, // USER INPUT DIRECT VIEW FACTOR MATRIX (N X N)
95 : const Array1D_int &SPtr, // pointer to actual surface number
96 : bool &NoUserInputF, // Flag signifying no input F's for this
97 : bool &ErrorsFound // True when errors are found in number of fields vs max args
98 : );
99 :
100 : void GetInputViewFactorsbyName(EnergyPlusData &state,
101 : std::string const &ZoneName, // Needed to check for user input view factors.
102 : int const N, // NUMBER OF SURFACES
103 : Array2A<Real64> F, // USER INPUT DIRECT VIEW FACTOR MATRIX (N X N)
104 : const Array1D_int &SPtr, // pointer to actual surface number
105 : bool &NoUserInputF, // Flag signifying no input F's for this
106 : bool &ErrorsFound // True when errors are found in number of fields vs max args
107 : );
108 :
109 : void CalcApproximateViewFactors(EnergyPlusData &state,
110 : int const N, // NUMBER OF SURFACES
111 : const Array1D<Real64> &A, // AREA VECTOR- ASSUMED,BE N ELEMENTS LONG
112 : const Array1D<Real64> &Azimuth, // Facing angle of the surface (in degrees)
113 : const Array1D<Real64> &Tilt, // Tilt angle of the surface (in degrees)
114 : Array2A<Real64> F, // APPROXIMATE DIRECT VIEW FACTOR MATRIX (N X N)
115 : const Array1D_int &SPtr // pointer to REAL(r64) surface number (for error message)
116 : );
117 :
118 : void FixViewFactors(EnergyPlusData &state,
119 : int const N, // NUMBER OF SURFACES
120 : const Array1D<Real64> &A, // AREA VECTOR- ASSUMED,BE N ELEMENTS LONG
121 : Array2A<Real64> F, // APPROXIMATE DIRECT VIEW FACTOR MATRIX (N X N)
122 : std::string &enclName, // Name of Enclosure being fixed
123 : std::vector<int> const &zoneNums, // Zones which are part of this enclosure
124 : Real64 &OriginalCheckValue, // check of SUM(F) - N
125 : Real64 &FixedCheckValue, // check after fixed of SUM(F) - N
126 : Real64 &FinalCheckValue, // the one to go with
127 : int &NumIterations, // number of iterations to fixed
128 : Real64 &RowSum, // RowSum of Fixed
129 : bool const anyIntMassInZone // are there any surfaces in the zone that are thermal mass
130 : );
131 :
132 : bool DoesZoneHaveInternalMass(EnergyPlusData &state,
133 : int const numZoneSurfaces, // Number of surfaces in the zone
134 : const Array1D_int &surfPointer // Pointers to the surfaces in the zone
135 : );
136 :
137 : void CalcScriptF(EnergyPlusData &state,
138 : int const N, // Number of surfaces
139 : Array1D<Real64> const &A, // AREA VECTOR- ASSUMED,BE N ELEMENTS LONG
140 : Array2<Real64> const &F, // DIRECT VIEW FACTOR MATRIX (N X N)
141 : Array1D<Real64> &EMISS, // VECTOR OF SURFACE EMISSIVITIES
142 : Array2<Real64> &ScriptF // MATRIX OF SCRIPT F FACTORS (N X N) //Tuned Transposed
143 : );
144 :
145 : void CalcFMRT(EnergyPlusData &state,
146 : int const N, // Number of surfaces
147 : Array1D<Real64> const &A, // AREA VECTOR- ASSUMED,BE N ELEMENTS LONG
148 : Array1D<Real64> &FMRT // VECTOR OF MEAN RADIANT TEMPERATURE "VIEW FACTORS"
149 : );
150 :
151 : void CalcFp(int const N, // Number of surfaces
152 : Array1D<Real64> const &EMISS, // VECTOR OF SURFACE EMISSIVITIES
153 : Array1D<Real64> const &FMRT, // VECTOR OF MEAN RADIANT TEMPERATURE "VIEW FACTORS"
154 : Array1D<Real64> &Fp // VECTOR OF OPPENHEIM RESISTNACE VALUES
155 : );
156 :
157 : void CalcMatrixInverse(Array2<Real64> &A, // Matrix: Gets reduced to L\U form
158 : Array2<Real64> &I // Returned as inverse matrix
159 : );
160 :
161 : int GetRadiantSystemSurface(EnergyPlusData &state,
162 : std::string const &cCurrentModuleObject, // Calling Object type
163 : std::string const &RadSysName, // Calling Object name
164 : int const RadSysZoneNum, // Radiant system zone number
165 : std::string const &SurfaceName, // Referenced surface name
166 : bool &ErrorsFound // True when errors are found
167 : );
168 :
169 : } // namespace HeatBalanceIntRadExchange
170 :
171 : struct HeatBalanceIntRadExchgData : BaseGlobalStruct
172 : {
173 :
174 : int MaxNumOfRadEnclosureSurfs = 0; // Max saved to get large enough space for SurfaceTempInKto4th
175 : bool CarrollMethod = false; // Use Carroll MRT method
176 : bool CalcInteriorRadExchangefirstTime = true; // Logical flag for one-time initializations
177 :
178 : // variables added as part of strategy to reduce calculation time - Glazer 2011-04-22
179 : Array1D<Real64> SurfaceTempRad;
180 : Array1D<Real64> SurfaceTempInKto4th;
181 : Array1D<Real64> SurfaceEmiss;
182 : bool ViewFactorReport = false; // Flag to output view factor report in eio file
183 : int LargestSurf = 0;
184 :
185 2126 : void init_constant_state([[maybe_unused]] EnergyPlusData &state) override
186 : {
187 2126 : }
188 :
189 1152 : void init_state([[maybe_unused]] EnergyPlusData &state) override
190 : {
191 1152 : }
192 :
193 2100 : void clear_state() override
194 : {
195 2100 : this->MaxNumOfRadEnclosureSurfs = 0;
196 2100 : this->CarrollMethod = false;
197 2100 : this->CalcInteriorRadExchangefirstTime = true;
198 2100 : this->SurfaceTempRad.deallocate();
199 2100 : this->SurfaceTempInKto4th.deallocate();
200 2100 : this->SurfaceEmiss.deallocate();
201 2100 : this->ViewFactorReport = false;
202 2100 : this->LargestSurf = 0;
203 2100 : }
204 : };
205 :
206 : } // namespace EnergyPlus
207 :
208 : #endif
|
