Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Corrects crash in DOAS fan sizing #10816

Open
wants to merge 3 commits into
base: develop
Choose a base branch
from
Open

Corrects crash in DOAS fan sizing #10816

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
6bf19d9
Correct crash in DOAS fan sizing
rraustad Nov 7, 2024
0d248a0
Add unit test
rraustad Nov 7, 2024
b249e6f
Revise OA system fan sizing methodology
rraustad Nov 8, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 0 additions & 2 deletions src/EnergyPlus/AirLoopHVACDOAS.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -994,13 +994,11 @@ namespace AirLoopHVACDOAS {
this->GetDesignDayConditions(state);

if (this->m_FanIndex > 0 && this->m_FanTypeNum == SimAirServingZones::CompType::Fan_System_Object) {
state.dataFans->fans(this->m_FanIndex)->maxAirFlowRate = sizingMassFlow / state.dataEnvrn->StdRhoAir;
state.dataLoopNodes->Node(this->m_FanInletNodeNum).MassFlowRateMaxAvail = sizingMassFlow;
state.dataLoopNodes->Node(this->m_FanOutletNodeNum).MassFlowRateMaxAvail = sizingMassFlow;
state.dataLoopNodes->Node(this->m_FanOutletNodeNum).MassFlowRateMax = sizingMassFlow;
}
if (this->m_FanIndex > 0 && this->m_FanTypeNum == SimAirServingZones::CompType::Fan_ComponentModel) {
state.dataFans->fans(this->m_FanIndex)->maxAirFlowRate = sizingMassFlow / state.dataEnvrn->StdRhoAir;
state.dataFans->fans(this->m_FanIndex)->minAirFlowRate = 0.0;
state.dataFans->fans(this->m_FanIndex)->maxAirMassFlowRate = sizingMassFlow;
state.dataLoopNodes->Node(this->m_FanInletNodeNum).MassFlowRateMaxAvail = sizingMassFlow;
Expand Down
126 changes: 68 additions & 58 deletions src/EnergyPlus/Autosizing/SystemAirFlowSizing.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -777,24 +777,57 @@ Real64 SystemAirFlowSizer::size(EnergyPlusData &state, Real64 _originalValue, bo
}

} else {
auto &sysSizPeakDDNum = state.dataSize->SysSizPeakDDNum(this->curSysNum);
if (this->curOASysNum > 0) {
// should a parent object ever override DOAS or OASystem flow rate (e.g., using oaSysEqSizing().AirVolFlow)?
if (outsideAirSys(this->curOASysNum).AirLoopDOASNum > -1) {
this->autoSizedValue =
this->airloopDOAS[outsideAirSys(this->curOASysNum).AirLoopDOASNum].SizingMassFlow / state.dataEnvrn->StdRhoAir;
} else {
this->autoSizedValue = this->finalSysSizing(this->curSysNum).DesOutAirVolFlow;
}
} else {
Copy link
Contributor Author

@rraustad rraustad Nov 7, 2024 •
edited
Loading

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The code change is from here to line 780 above. The below is just git confusion.

auto &sysSizPeakDDNum = state.dataSize->SysSizPeakDDNum(this->curSysNum);

if (this->airLoopSysFlag) {
auto const &unitarysysEqSizing = this->unitarySysEqSizing(this->curSysNum);
if (unitarysysEqSizing.CoolingAirFlow && unitarysysEqSizing.HeatingAirFlow) {
this->autoSizedValue = std::max(unitarysysEqSizing.CoolingAirVolFlow, unitarysysEqSizing.HeatingAirVolFlow);
if (this->autoSizedValue == unitarysysEqSizing.CoolingAirVolFlow) {
if (sysSizPeakDDNum.CoolFlowPeakDD > 0 && sysSizPeakDDNum.CoolFlowPeakDD <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(sysSizPeakDDNum.CoolFlowPeakDD);
DDNameFanPeak = desDayInput.Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, sysSizPeakDDNum.TimeStepAtCoolFlowPk(sysSizPeakDDNum.CoolFlowPeakDD)));
}

if (this->airLoopSysFlag) {
auto const &unitarysysEqSizing = this->unitarySysEqSizing(this->curSysNum);
if (unitarysysEqSizing.CoolingAirFlow && unitarysysEqSizing.HeatingAirFlow) {
this->autoSizedValue = std::max(unitarysysEqSizing.CoolingAirVolFlow, unitarysysEqSizing.HeatingAirVolFlow);
if (this->autoSizedValue == unitarysysEqSizing.CoolingAirVolFlow) {
} else if (this->autoSizedValue == unitarysysEqSizing.HeatingAirVolFlow) {
if (this->finalSysSizing(this->curSysNum).HeatDDNum > 0 &&
this->finalSysSizing(this->curSysNum).HeatDDNum <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(this->finalSysSizing(this->curSysNum).HeatDDNum);
DDNameFanPeak = desDayInput.Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, this->finalSysSizing(this->curSysNum).SysHeatAirTimeStepPk));
}
}
} else if (unitarysysEqSizing.CoolingAirFlow) {
this->autoSizedValue = unitarysysEqSizing.CoolingAirVolFlow;
if (sysSizPeakDDNum.CoolFlowPeakDD > 0 && sysSizPeakDDNum.CoolFlowPeakDD <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(sysSizPeakDDNum.CoolFlowPeakDD);
DDNameFanPeak = desDayInput.Title;
DDNameFanPeak = state.dataWeather->DesDayInput(sysSizPeakDDNum.CoolFlowPeakDD).Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, sysSizPeakDDNum.TimeStepAtCoolFlowPk(sysSizPeakDDNum.CoolFlowPeakDD)));
}

} else if (this->autoSizedValue == unitarysysEqSizing.HeatingAirVolFlow) {
} else if (unitarysysEqSizing.HeatingAirFlow) {
this->autoSizedValue = unitarysysEqSizing.HeatingAirVolFlow;
if (this->finalSysSizing(this->curSysNum).HeatDDNum > 0 &&
this->finalSysSizing(this->curSysNum).HeatDDNum <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(this->finalSysSizing(this->curSysNum).HeatDDNum);
Expand All @@ -805,31 +838,32 @@ Real64 SystemAirFlowSizer::size(EnergyPlusData &state, Real64 _originalValue, bo
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, this->finalSysSizing(this->curSysNum).SysHeatAirTimeStepPk));
}
}
} else if (unitarysysEqSizing.CoolingAirFlow) {
this->autoSizedValue = unitarysysEqSizing.CoolingAirVolFlow;
if (sysSizPeakDDNum.CoolFlowPeakDD > 0 && sysSizPeakDDNum.CoolFlowPeakDD <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(sysSizPeakDDNum.CoolFlowPeakDD);
DDNameFanPeak = state.dataWeather->DesDayInput(sysSizPeakDDNum.CoolFlowPeakDD).Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, sysSizPeakDDNum.TimeStepAtCoolFlowPk(sysSizPeakDDNum.CoolFlowPeakDD)));
}
} else if (unitarysysEqSizing.HeatingAirFlow) {
this->autoSizedValue = unitarysysEqSizing.HeatingAirVolFlow;
if (this->finalSysSizing(this->curSysNum).HeatDDNum > 0 &&
this->finalSysSizing(this->curSysNum).HeatDDNum <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(this->finalSysSizing(this->curSysNum).HeatDDNum);
DDNameFanPeak = desDayInput.Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, this->finalSysSizing(this->curSysNum).SysHeatAirTimeStepPk));
}

} else {
this->autoSizedValue = this->finalSysSizing(this->curSysNum).DesMainVolFlow;
if (this->autoSizedValue == this->finalSysSizing(this->curSysNum).DesHeatVolFlow) {
if (this->finalSysSizing(this->curSysNum).HeatDDNum > 0 &&
this->finalSysSizing(this->curSysNum).HeatDDNum <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(this->finalSysSizing(this->curSysNum).HeatDDNum);
DDNameFanPeak = desDayInput.Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, this->finalSysSizing(this->curSysNum).SysHeatAirTimeStepPk));
}
} else if (this->autoSizedValue == this->finalSysSizing(this->curSysNum).DesCoolVolFlow) {
if (sysSizPeakDDNum.CoolFlowPeakDD > 0 && sysSizPeakDDNum.CoolFlowPeakDD <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(sysSizPeakDDNum.CoolFlowPeakDD);
DDNameFanPeak = desDayInput.Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, sysSizPeakDDNum.TimeStepAtCoolFlowPk(sysSizPeakDDNum.CoolFlowPeakDD)));
}
}
}
} else {
this->autoSizedValue = this->finalSysSizing(this->curSysNum).DesMainVolFlow;
if (this->autoSizedValue == this->finalSysSizing(this->curSysNum).DesHeatVolFlow) {
Expand All @@ -855,32 +889,8 @@ Real64 SystemAirFlowSizer::size(EnergyPlusData &state, Real64 _originalValue, bo
}
}
}
} else {
this->autoSizedValue = this->finalSysSizing(this->curSysNum).DesMainVolFlow;
if (this->autoSizedValue == this->finalSysSizing(this->curSysNum).DesHeatVolFlow) {
if (this->finalSysSizing(this->curSysNum).HeatDDNum > 0 &&
this->finalSysSizing(this->curSysNum).HeatDDNum <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(this->finalSysSizing(this->curSysNum).HeatDDNum);
DDNameFanPeak = desDayInput.Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, this->finalSysSizing(this->curSysNum).SysHeatAirTimeStepPk));
}
} else if (this->autoSizedValue == this->finalSysSizing(this->curSysNum).DesCoolVolFlow) {
if (sysSizPeakDDNum.CoolFlowPeakDD > 0 && sysSizPeakDDNum.CoolFlowPeakDD <= state.dataEnvrn->TotDesDays) {
auto &desDayInput = state.dataWeather->DesDayInput(sysSizPeakDDNum.CoolFlowPeakDD);
DDNameFanPeak = desDayInput.Title;
dateTimeFanPeak = format("{}/{} {}",
desDayInput.Month,
desDayInput.DayOfMonth,
state.dataRptCoilSelection->coilSelectionReportObj->getTimeText(
state, sysSizPeakDDNum.TimeStepAtCoolFlowPk(sysSizPeakDDNum.CoolFlowPeakDD)));
}
}
if (this->dataFractionUsedForSizing > 0.0) this->autoSizedValue = this->autoSizedValue * this->dataFractionUsedForSizing;
}
if (this->dataFractionUsedForSizing > 0.0) this->autoSizedValue = this->autoSizedValue * this->dataFractionUsedForSizing;
}
}
} else if (this->dataNonZoneNonAirloopValue > 0) {
Expand Down
19 changes: 19 additions & 0 deletions tst/EnergyPlus/unit/Autosizing/SystemAirFlowSizing.unit.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -960,6 +960,25 @@ TEST_F(AutoSizingFixture, SystemAirFlowSizingGauntlet)
eiooutput = std::string(" Component Sizing Information, Coil:Heating:Water, MyWaterCoil, Design Size Maximum Flow Rate [m3/s], 3.50000\n"
" Component Sizing Information, Coil:Heating:Water, MyWaterCoil, User-Specified Maximum Flow Rate [m3/s], 2.20000\n");
EXPECT_TRUE(compare_eio_stream(eiooutput, true));

// Test 58 - DOAS fan air flow rate autosizing
state->dataSize->HRFlowSizingFlag = false;
state->dataSize->CurOASysNum = 1;
state->dataAirLoop->OutsideAirSys(1).AirLoopDOASNum = 0;
AirLoopHVACDOAS::AirLoopDOAS thisDOAS;
state->dataAirLoopHVACDOAS->airloopDOAS.push_back(thisDOAS);
// use 0.53 m3/s multiplied by StdRhoAir = 1.2
state->dataAirLoopHVACDOAS->airloopDOAS[state->dataAirLoop->OutsideAirSys(1).AirLoopDOASNum].SizingMassFlow = 0.53 * 1.2;
inputValue = DataSizing::AutoSize;

// do sizing
sizer.wasAutoSized = false;
printFlag = true;
sizer.initializeWithinEP(*this->state, "Fan:SystemModel", "MyDOASFan", printFlag, routineName);
sizedValue = sizer.size(*this->state, inputValue, errorsFound);
EXPECT_ENUM_EQ(AutoSizingResultType::NoError, sizer.errorType); // cumulative of previous calls
EXPECT_TRUE(sizer.wasAutoSized);
EXPECT_NEAR(0.53, sizedValue, 0.01); // auto-sized value
}

} // namespace EnergyPlus
Loading