Line data Source code
1 : // EnergyPlus, Copyright (c) 1996-2024, 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 : // C++ Headers
49 : #include <cmath>
50 :
51 : // EnergyPlus Headers
52 : #include <EnergyPlus/Autosizing/HeatingWaterDesCoilLoadUsedForUASizing.hh>
53 : #include <EnergyPlus/Data/EnergyPlusData.hh>
54 : #include <EnergyPlus/DataEnvironment.hh>
55 : #include <EnergyPlus/FluidProperties.hh>
56 : #include <EnergyPlus/Psychrometrics.hh>
57 :
58 : namespace EnergyPlus {
59 :
60 2388 : Real64 HeatingWaterDesCoilLoadUsedForUASizer::size(EnergyPlusData &state, Real64 _originalValue, bool &errorsFound)
61 : {
62 2388 : if (!this->checkInitialized(state, errorsFound)) {
63 0 : return 0.0;
64 : }
65 2388 : this->preSize(state, _originalValue);
66 :
67 2388 : if (this->curZoneEqNum > 0) {
68 1996 : if (!this->wasAutoSized && !this->sizingDesRunThisZone) {
69 0 : this->autoSizedValue = _originalValue;
70 : } else {
71 1996 : if (this->termUnitSingDuct && (this->curTermUnitSizingNum > 0)) {
72 3664 : Real64 const Cp = FluidProperties::GetSpecificHeatGlycol(state,
73 1832 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidName,
74 : Constant::HWInitConvTemp,
75 1832 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidIndex,
76 : this->callingRoutine);
77 3664 : Real64 const rho = FluidProperties::GetDensityGlycol(state,
78 1832 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidName,
79 : Constant::HWInitConvTemp,
80 1832 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidIndex,
81 : this->callingRoutine);
82 1832 : this->autoSizedValue = this->dataWaterFlowUsedForSizing * this->dataWaterCoilSizHeatDeltaT * Cp * rho;
83 1832 : state.dataRptCoilSelection->coilSelectionReportObj->setCoilReheatMultiplier(state, this->compName, this->compType, 1.0);
84 1996 : } else if ((this->termUnitPIU || this->termUnitIU) && (this->curTermUnitSizingNum > 0)) {
85 120 : Real64 const Cp = FluidProperties::GetSpecificHeatGlycol(state,
86 60 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidName,
87 : Constant::HWInitConvTemp,
88 60 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidIndex,
89 : this->callingRoutine);
90 120 : Real64 const rho = FluidProperties::GetDensityGlycol(state,
91 60 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidName,
92 : Constant::HWInitConvTemp,
93 60 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidIndex,
94 : this->callingRoutine);
95 120 : this->autoSizedValue = this->dataWaterFlowUsedForSizing * this->dataWaterCoilSizHeatDeltaT * Cp * rho *
96 60 : this->termUnitSizing(this->curTermUnitSizingNum).ReheatLoadMult;
97 164 : } else if (this->zoneEqFanCoil || this->zoneEqUnitHeater) {
98 166 : Real64 const Cp = FluidProperties::GetSpecificHeatGlycol(state,
99 83 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidName,
100 : Constant::HWInitConvTemp,
101 83 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidIndex,
102 : this->callingRoutine);
103 166 : Real64 const rho = FluidProperties::GetDensityGlycol(state,
104 83 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidName,
105 : Constant::HWInitConvTemp,
106 83 : state.dataPlnt->PlantLoop(this->dataWaterLoopNum).FluidIndex,
107 : this->callingRoutine);
108 83 : this->autoSizedValue = this->dataWaterFlowUsedForSizing * this->dataWaterCoilSizHeatDeltaT * Cp * rho;
109 83 : state.dataRptCoilSelection->coilSelectionReportObj->setCoilReheatMultiplier(state, this->compName, this->compType, 1.0);
110 83 : } else {
111 21 : Real64 DesMassFlow = 0.0;
112 21 : if (this->zoneEqSizing(this->curZoneEqNum).SystemAirFlow) {
113 0 : DesMassFlow = this->zoneEqSizing(this->curZoneEqNum).AirVolFlow * state.dataEnvrn->StdRhoAir;
114 21 : } else if (this->zoneEqSizing(this->curZoneEqNum).HeatingAirFlow) {
115 12 : DesMassFlow = this->zoneEqSizing(this->curZoneEqNum).HeatingAirVolFlow * state.dataEnvrn->StdRhoAir;
116 : } else {
117 9 : DesMassFlow = this->finalZoneSizing(this->curZoneEqNum).DesHeatMassFlow;
118 : }
119 : Real64 const CoilInTemp =
120 21 : this->setHeatCoilInletTempForZoneEqSizing(setOAFracForZoneEqSizing(state, DesMassFlow, this->zoneEqSizing(this->curZoneEqNum)),
121 21 : this->zoneEqSizing(this->curZoneEqNum),
122 21 : this->finalZoneSizing(this->curZoneEqNum));
123 : // Real64 CoilInHumRat =
124 : // this->setHeatCoilInletHumRatForZoneEqSizing(setOAFracForZoneEqSizing(DesMassFlow, this->zoneEqSizing(this->curZoneEqNum)),
125 : // this->zoneEqSizing(this->curZoneEqNum),
126 : // this->finalZoneSizing(this->curZoneEqNum));
127 21 : Real64 const CoilOutTemp = this->finalZoneSizing(this->curZoneEqNum).HeatDesTemp;
128 21 : Real64 const CoilOutHumRat = this->finalZoneSizing(this->curZoneEqNum).HeatDesHumRat;
129 21 : this->autoSizedValue = Psychrometrics::PsyCpAirFnW(CoilOutHumRat) * DesMassFlow * (CoilOutTemp - CoilInTemp);
130 : }
131 : }
132 392 : } else if (this->curSysNum > 0) {
133 392 : if (!this->wasAutoSized && !this->sizingDesRunThisAirSys) {
134 0 : this->autoSizedValue = _originalValue;
135 : } else {
136 392 : Real64 OutAirFrac = 1.0;
137 392 : if (this->curOASysNum > 0) {
138 31 : OutAirFrac = 1.0;
139 361 : } else if (this->finalSysSizing(this->curSysNum).HeatOAOption == this->minOA) {
140 345 : if (this->dataAirFlowUsedForSizing > 0.0) {
141 345 : OutAirFrac = this->finalSysSizing(this->curSysNum).DesOutAirVolFlow / this->dataAirFlowUsedForSizing;
142 345 : OutAirFrac = min(1.0, max(0.0, OutAirFrac));
143 : } else {
144 0 : OutAirFrac = 1.0;
145 : }
146 : } else {
147 16 : OutAirFrac = 1.0;
148 : }
149 392 : Real64 CoilInTemp = 0.0;
150 392 : if (this->curOASysNum == 0 && this->primaryAirSystem(this->curSysNum).NumOAHeatCoils > 0) {
151 28 : CoilInTemp = OutAirFrac * this->finalSysSizing(this->curSysNum).PreheatTemp +
152 28 : (1.0 - OutAirFrac) * this->finalSysSizing(this->curSysNum).HeatRetTemp;
153 364 : } else if (this->curOASysNum > 0 && this->outsideAirSys(this->curOASysNum).AirLoopDOASNum > -1) {
154 1 : CoilInTemp = this->airloopDOAS[this->outsideAirSys(this->curOASysNum).AirLoopDOASNum].HeatOutTemp;
155 : } else {
156 363 : CoilInTemp = OutAirFrac * this->finalSysSizing(this->curSysNum).HeatOutTemp +
157 363 : (1.0 - OutAirFrac) * this->finalSysSizing(this->curSysNum).HeatRetTemp;
158 : }
159 : // coil load
160 392 : Real64 const CpAirStd = Psychrometrics::PsyCpAirFnW(0.0);
161 392 : if (this->curOASysNum > 0) {
162 31 : if (this->dataDesicRegCoil) {
163 0 : this->autoSizedValue = CpAirStd * state.dataEnvrn->StdRhoAir * this->dataAirFlowUsedForSizing *
164 0 : (this->dataDesOutletAirTemp - this->dataDesInletAirTemp);
165 31 : } else if (this->outsideAirSys(this->curOASysNum).AirLoopDOASNum > -1) {
166 1 : this->autoSizedValue = CpAirStd * state.dataEnvrn->StdRhoAir * this->dataAirFlowUsedForSizing *
167 1 : (this->airloopDOAS[this->outsideAirSys(this->curOASysNum).AirLoopDOASNum].PreheatTemp - CoilInTemp);
168 : } else {
169 30 : this->autoSizedValue = CpAirStd * state.dataEnvrn->StdRhoAir * this->dataAirFlowUsedForSizing *
170 30 : (this->finalSysSizing(this->curSysNum).PreheatTemp - CoilInTemp);
171 : }
172 : } else {
173 361 : if (this->finalSysSizing(this->curSysNum).HeatingCapMethod == DataSizing::FractionOfAutosizedHeatingCapacity) {
174 1 : this->dataFracOfAutosizedHeatingCapacity = this->finalSysSizing(this->curSysNum).FractionOfAutosizedHeatingCapacity;
175 : }
176 361 : if (this->dataDesicRegCoil) {
177 0 : this->autoSizedValue = CpAirStd * state.dataEnvrn->StdRhoAir * this->dataAirFlowUsedForSizing *
178 0 : (this->dataDesOutletAirTemp - this->dataDesInletAirTemp);
179 : } else {
180 361 : this->autoSizedValue = CpAirStd * state.dataEnvrn->StdRhoAir * this->dataAirFlowUsedForSizing *
181 361 : (this->finalSysSizing(this->curSysNum).HeatSupTemp - CoilInTemp);
182 : }
183 : }
184 : }
185 : }
186 : // heating coil can't have negative capacity
187 2388 : this->autoSizedValue = std::max(0.0, this->autoSizedValue) * this->dataHeatSizeRatio * this->dataFracOfAutosizedHeatingCapacity;
188 2388 : if (this->overrideSizeString) {
189 2388 : if (this->isEpJSON) {
190 0 : this->sizingString = "water_heating_design_coil_load_for_ua_sizing";
191 : }
192 : }
193 2388 : this->selectSizerOutput(state, errorsFound);
194 2388 : if (this->isCoilReportObject && this->curSysNum <= this->numPrimaryAirSys) {
195 2387 : Real64 constexpr FanCoolLoad = 0.0;
196 2387 : Real64 constexpr TotCapTempModFac = 1.0;
197 2387 : Real64 constexpr DXFlowPerCapMinRatio = 1.0;
198 2387 : Real64 constexpr DXFlowPerCapMaxRatio = 1.0;
199 2387 : state.dataRptCoilSelection->coilSelectionReportObj->setCoilHeatingCapacity(state,
200 2387 : this->compName,
201 2387 : this->compType,
202 : this->autoSizedValue,
203 2387 : this->wasAutoSized,
204 : this->curSysNum,
205 : this->curZoneEqNum,
206 : this->curOASysNum,
207 : FanCoolLoad,
208 : TotCapTempModFac,
209 : DXFlowPerCapMinRatio,
210 : DXFlowPerCapMaxRatio);
211 : }
212 2388 : return this->autoSizedValue;
213 : }
214 :
215 : } // namespace EnergyPlus
|
