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.

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

Correct VAV maximum air flow fraction during reheat for heating dominated use cases #10763

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

Conversation

rraustad
Copy link
Contributor

@rraustad rraustad commented Sep 24, 2024

Pull request overview

NOTE: ENHANCEMENTS MUST FOLLOW A SUBMISSION PROCESS INCLUDING A FEATURE PROPOSAL AND DESIGN DOCUMENT PRIOR TO SUBMITTING CODE

Pull Request Author

Add to this list or remove from it as applicable. This is a simple templated set of guidelines.

  • Title of PR should be user-synopsis style (clearly understandable in a standalone changelog context)
  • Label the PR with at least one of: Defect, Refactoring, NewFeature, Performance, and/or DoNoPublish
  • Pull requests that impact EnergyPlus code must also include unit tests to cover enhancement or defect repair
  • Author should provide a "walkthrough" of relevant code changes using a GitHub code review comment process
  • If any diffs are expected, author must demonstrate they are justified using plots and descriptions
  • If changes fix a defect, the fix should be demonstrated in plots and descriptions
  • If any defect files are updated to a more recent version, upload new versions here or on DevSupport
  • If IDD requires transition, transition source, rules, ExpandObjects, and IDFs must be updated, and add IDDChange label
  • If structural output changes, add to output rules file and add OutputChange label
  • If adding/removing any LaTeX docs or figures, update that document's CMakeLists file dependencies

Reviewer

This will not be exhaustively relevant to every PR.

  • Perform a Code Review on GitHub
  • If branch is behind develop, merge develop and build locally to check for side effects of the merge
  • If defect, verify by running develop branch and reproducing defect, then running PR and reproducing fix
  • If feature, test running new feature, try creative ways to break it
  • CI status: all green or justified
  • Check that performance is not impacted (CI Linux results include performance check)
  • Run Unit Test(s) locally
  • Check any new function arguments for performance impacts
  • Verify IDF naming conventions and styles, memos and notes and defaults
  • If new idf included, locally check the err file and other outputs

@rraustad rraustad added the Defect Includes code to repair a defect in EnergyPlus label Sep 24, 2024
@rraustad
Copy link
Contributor Author

The defect file results are now:

Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE1-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE2-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 0.79773
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE3-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE4-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 0.69789
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE5-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000

@rraustad rraustad changed the title Correct VAV maximum air flow fraction during reheat Correct VAV maximum air flow fraction during reheat for heating dominated use cases Sep 24, 2024
Real64 heatingMaxFlow = (this->DamperHeatingAction == Action::ReverseWithLimits)
? state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesHeatVolFlowMax
: state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesHeatVolFlow;
MaxAirVolFlowRateDes = max(state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesCoolVolFlow, heatingMaxFlow);
Copy link
Contributor Author

Choose a reason for hiding this comment

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

If ReverseWithLimits uses a different heating air flow rate from the zone design heating air flow rate then that flow rate should be used as the value for the VAV box heating air flow rate.

Copy link
Contributor Author

@rraustad rraustad Sep 24, 2024

Choose a reason for hiding this comment

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

Well, I was expecting DesHeatVolFlowMax to always be less than DesHeatVolFlow. For 5ZoneBoilerOutsideAirReset that's not the case. The diffs show that Zone 5 TU max air flow rate increased, I was only expecting decreases in max air flow rate for heating dominated zones. Zone cooling design air flow is greater than zone heating design air flow, but DesHeatVolFlowMax is larger than both of these.

image

image

// set the zone maximum heating supply air flow rate. This will be used for autosizing VAV terminal unit
// max heating flow rates
zsFinalSizing.DesHeatVolFlowMax = max(zsFinalSizing.DesHeatMaxAirFlow,
                                      zsFinalSizing.DesHeatMaxAirFlow2,
                                      max(zsFinalSizing.DesCoolVolFlow, zsFinalSizing.DesHeatVolFlow) * zsFinalSizing.DesHeatMaxAirFlowFrac);

zsFinalSizing.DesHeatMaxAirFlow2 = zsFinalSizing.DesHeatMaxAirFlowPerArea * floorArea * zoneMult;
                                 = 0.002032 * 182.49 * 1 = 0.37082

N8 , \field Maximum Flow Fraction During Reheat
     \type real
     \autosizable
     \default autosize
     \note Used only when Reheat Coil Object Type = Coil:Heating:Water and Damper Heating Action = ReverseWithLimits
     \note When autocalculating, the maximum flow fraction is set to the ratio of
     \note 0.002032 m3/s-m2 (0.4 cfm/sqft) multiplied by the zone floor area and the
     \note Maximum Air Flow Rate.
     \note This optional field limits the maximum flow allowed in reheat mode.
     \note At no time will the maximum flow rate calculated here exceed the value of
     \note Maximum Air Flow Rate.

So is DesHeatVolFlowMax just an upper limit? or is it the expected heating air flow rate?

Copy link
Contributor Author

@rraustad rraustad Sep 24, 2024

Choose a reason for hiding this comment

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

So, since DesHeatVolFlowMax looks like a limit and not a flow rate, maybe I should change this line and line 3224 to?:

Real64 heatingMaxFlow = (DesHeatVolFlowMax > DesHeatVolFlow) ? state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesHeatVolFlow :
    state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesHeatVolFlowMax);

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Based on this discussion I changed to treat DesHeatVolFlowMax as an upper limit. This should still correct the issue and result in fewer diffs.

MaxHeatAirVolFlowRateDes = state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesHeatVolFlowMax;
} else {
MaxHeatAirVolFlowRateDes = state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesHeatVolFlow;
}
Copy link
Contributor Author

Choose a reason for hiding this comment

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

Below at line 3224, ReverseWithLimits uses DesHeatVolFlowMax. That value should also be used here.

3224: MaxAirVolFlowRateDuringReheatDes = state.dataSize->TermUnitFinalZoneSizing(state.dataSize->CurTermUnitSizingNum).DesHeatVolFlowMax;

Copy link
Contributor Author

@rraustad rraustad Sep 25, 2024

Choose a reason for hiding this comment

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

As discussed above, these 2 areas of code, above (2937) and below (3231), were changed to use DesHeatVolFlowMax as a limit.

Copy link

⚠️ Regressions detected on macos-14 for commit 64c36ca

Regression Summary
  • EIO: 83
  • ESO Big Diffs: 77
  • MTR Big Diffs: 68
  • Table Big Diffs: 82
  • Table String Diffs: 76
  • ERR: 15
  • Audit: 2
  • EDD: 1

Copy link

⚠️ Regressions detected on macos-14 for commit 1979f09

Regression Summary
  • EIO: 105
  • Table Big Diffs: 105
  • ESO Small Diffs: 11
  • MTR Small Diffs: 11
  • Table String Diffs: 34
  • ESO Big Diffs: 26
  • MTR Big Diffs: 18
  • ERR: 6
  • EDD: 2

@rraustad
Copy link
Contributor Author

rraustad commented Sep 24, 2024

It turns out this does correct a reporting issue (5ZoneAutoDXVAV).

Zone Sizing Information, SPACE1-1, Cooling, 4433.93700, 4433.93700, 0.46734, 0.46734
Zone Sizing Information, SPACE1-1, Heating, 3564.64198, 3564.64198, 0.10362, 0.10362

ReverseWithLimits will operate up to the zone design heating air flow unless that flow exceeds the maximum limit. In the past the maximum limit was used to calculate the maximum reheat fraction. This change only uses the maximum reheat flow rate if it actually limits the design heating air flow rate.

In this case the maximum reheat fraction is 0.10362 m3/s / 0.46734 m3/s = 0.2217229 <- correct answer
The previous code used the maximum reheat fraction of 99.16 m2 * 0.002032 = 0.201493 m3/s / 0.46734 m3/.s = 0.431149

and 0.10362 m3/s / 99.16 m2 = 0.001049778 m3/s/m2 <- correct answer

image

Copy link

⚠️ Regressions detected on macos-14 for commit b275d42

Regression Summary
  • EIO: 105
  • Table Big Diffs: 105
  • ESO Small Diffs: 11
  • MTR Small Diffs: 11
  • Table String Diffs: 34
  • ESO Big Diffs: 26
  • MTR Big Diffs: 18
  • ERR: 6
  • EDD: 2

@rraustad
Copy link
Contributor Author

rraustad commented Sep 24, 2024

Defect file results after latest change to use DesHeatVolFlowMax as a limit, not a flow rate:

Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE1-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE2-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 0.60486
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE3-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE4-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 0.52735
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE5-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000

expectedMaxAirVolFractionDuringReheat = state->dataSize->TermUnitFinalZoneSizing(1).DesHeatVolFlowMax / thisSys.MaxAirVolFlowRate;
EXPECT_EQ(expectedZoneMinAirFracDes, thisSys.ZoneMinAirFracDes);
EXPECT_EQ(1.0, thisSys.MaxAirVolFractionDuringReheat);
EXPECT_EQ(expectedMaxAirVolFractionDuringReheat, thisSys.MaxAirVolFractionDuringReheat);
Copy link
Contributor Author

Choose a reason for hiding this comment

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

Added 3 new tests for maximum air flow fraction during reheat to this unit test.

@rraustad rraustad added this to the EnergyPlus 25.1 milestone Sep 25, 2024
@rraustad
Copy link
Contributor Author

I changed the defect file to use Damper Heating Action = Reverse and get expected results regarding the Maximum Flow Fraction during Reheat and the air flow rate used to size the heating coil.

Zone Sizing Information, SPACE1-1, Cooling, 1991.75827, 1991.75827, 0.16753, 0.16753
Zone Sizing Information, SPACE1-1, Heating, 3970.36898, 3970.36898, 0.26058, 0.26058
Zone Sizing Information, SPACE2-1, Cooling, 2110.05356, 2110.05356, 0.17747, 0.17747
Zone Sizing Information, SPACE2-1, Heating, 1635.59910, 1635.59910, 0.10735, 0.10735
Zone Sizing Information, SPACE3-1, Cooling, 1822.43896, 1822.43896, 0.15328, 0.15328
Zone Sizing Information, SPACE3-1, Heating, 3862.85865, 3862.85865, 0.25352, 0.25352
Zone Sizing Information, SPACE4-1, Cooling, 2411.91565, 2411.91565, 0.20286, 0.20286
Zone Sizing Information, SPACE4-1, Heating, 1630.01391, 1630.01391, 0.10698, 0.10698
Zone Sizing Information, SPACE5-1, Cooling, 1144.55309, 1144.55309, 9.62850E-002, 0.13906
Zone Sizing Information, SPACE5-1, Heating, 3259.03012, 3259.03012, 0.21389, 0.21389

Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE1-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE1-1 VAV REHEAT, Design Size Reheat Coil Sizing Air Volume Flow Rate [m3/s], 0.26058
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE2-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE2-1 VAV REHEAT, Design Size Reheat Coil Sizing Air Volume Flow Rate [m3/s], 0.10735
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE3-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE3-1 VAV REHEAT, Design Size Reheat Coil Sizing Air Volume Flow Rate [m3/s], 0.25352
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE4-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE4-1 VAV REHEAT, Design Size Reheat Coil Sizing Air Volume Flow Rate [m3/s], 0.10698
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE5-1 VAV REHEAT, Design Size Maximum Flow Fraction during Reheat [], 1.00000
Component Sizing Information, AirTerminal:SingleDuct:VAV:Reheat, SPACE5-1 VAV REHEAT, Design Size Reheat Coil Sizing Air Volume Flow Rate [m3/s], 0.21389

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Defect Includes code to repair a defect in EnergyPlus
Projects
None yet
5 participants