Skip to content

Fix exterior convection heat transfer coefficient for specified temperature boundary condition #443

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

Merged
merged 8 commits into from
Oct 19, 2023
Merged

Fix exterior convection heat transfer coefficient for specified temperature boundary condition #443

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
c8177c7
Include sb_hcFrc in sb_hxa calculation.
nathancoliver Oct 17, 2023
b660d9a
Add comment to consider adding exterior wind speed user input.
nathancoliver Oct 17, 2023
1e950e0
Use sb_HCAmbient() to calculate sb_hcNat.
nathancoliver Oct 17, 2023
66bc93b
Define sfExHcMult and set sfExHcModel=Input when sfExCnd=SpecifiedT.
nathancoliver Oct 17, 2023
ede178e
Update files.
nathancoliver Oct 17, 2023
8d8ab28
Explicitly set sb_hc_Frc = 0 when sfExCnd = SPECIFIEDT and sfExHcMode…
nathancoliver Oct 19, 2023
ddf468a
Revert "Update files."
nathancoliver Oct 19, 2023
b568d56
Update reference test files.
nathancoliver Oct 19, 2023
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
6 changes: 3 additions & 3 deletions src/cncult3.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2704,9 +2704,8 @@ void SBC::sb_SetCoeffs( // set convective and radiant coefficients
else if (sb_pXS->sfExCnd==C_EXCNDCH_SPECT)
{ sb_txa = sb_pXS->sfExT;
sb_txr = sb_txa;
sb_hcNat = sb_pXS->uX;
sb_hcFrc = 0.f;
sb_hxa = sb_hcMult * sb_hcNat;
sb_HCAmbient();
sb_hxa = sb_hcMult * (sb_hcNat + sb_hcFrc);
sb_hxr = 0.f;
}
else if (sb_pXS->sfExCnd == C_EXCNDCH_GROUND)
Expand Down Expand Up @@ -2939,6 +2938,7 @@ x printf( "Hit\n");

case C_CONVMODELCH_INPUT:
sb_hcNat = 1.f;
sb_hcFrc = 0.f;
break;

default:
Expand Down
4 changes: 2 additions & 2 deletions test/1ZAttic.cse
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -201,8 +201,8 @@ ZONE Z1
sfExT=50 //10 C ground temperature
sfInAbs=INABSORP
sfInHcModel = Unified //floor inside film coefficient
sfExH=100000 //floor outside film coefficient
sfExHcModel = Unified
sfExHcMult=100000 //floor outside film coefficient
sfExHcModel = INPUT

SURFACE North_Wall
sfType=Wall
Expand Down
8 changes: 4 additions & 4 deletions test/2ZAttic.cse
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -201,8 +201,8 @@ ZONE Z1
sfExT=50 //10 C ground temperature
sfInAbs=INABSORP
sfInHcModel = Unified //floor inside film coefficient
sfExH=100000 //floor outside film coefficient
sfExHcModel = Unified
sfExHcMult=100000 //floor outside film coefficient
sfExHcModel = INPUT

SURFACE North_Wall
sfType=Wall
Expand Down Expand Up @@ -390,8 +390,8 @@ ZONE Z2
sfExT=50 //10 C ground temperature
sfInAbs=INABSORP
sfInHcModel = Unified //floor inside film coefficient
sfExH=100000 //floor outside film coefficient
sfExHcModel = Unified
sfExHcMult=100000 //floor outside film coefficient
sfExHcModel = INPUT

SURFACE North_Wall
sfType=Wall
Expand Down
3 changes: 2 additions & 1 deletion test/CRTerminalAH.cse
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -191,7 +191,8 @@ ZONE Z1
sfExT=50 //10 C ground temperature
sfInAbs=INABSORP
sfInHcModel = Unified //floor inside film coefficient
sfExH=100000 //floor outside film coefficient
sfExHcMult=100000 //floor outside film coefficient
sfExHcModel = INPUT
sfModel = SFMODELX

SURFACE North_Wall
Expand Down
16 changes: 8 additions & 8 deletions test/Zonetest4W.cse
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,7 @@ MATERIAL "stucco" matCond=0.4167 matSpHt=0.2 matDens=116. // JA

//Constructions

CONSTRUCTION WCav13 // �W.13.2X4 BAW 9/2
CONSTRUCTION WCav13 // �W.13.2X4 BAW 9/2
layer lrMat="Gypbd" lrThk=0.25/12.
layer lrMat="Gypbd" lrThk=0.25/12.
layer lrMat="R11" lrThk=3.5/(12.*2)
Expand Down Expand Up @@ -153,16 +153,16 @@ Zone Single znArea=100 znCAir=18.3 infAC=0 infShld=5 infStories=1 infELA=0 znVol
znQMxH = 1000000 znQMxC = -500000 //
GAIN gnPower= select( $hour >12,0.,default 0.0) gnFrRad=0 //

SURFACE "North" sfType=Wall sfCon=WCav13 sfArea=10 sfAzm=Front sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExH=1.5 //
SURFACE "South" sfType=Wall sfCon=WCav13 sfArea=100 sfAzm=Back sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExH=1.5 //
SURFACE "North" sfType=Wall sfCon=WCav13 sfArea=10 sfAzm=Front sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExHcMult=1.5 sfExHcModel=INPUT//
SURFACE "South" sfType=Wall sfCon=WCav13 sfArea=100 sfAzm=Back sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExHcMult=1.5 sfExHcModel=INPUT//

#redefine WinU 0.4 // Window NFRC U-factor
WINDOW wnHeight=10; wnWidth=5 wnGt=GT wnSMSO=0.; wnSMSC=0.; wnU=(1./((1./WinU)-0.85)) wnInH=1.5 // BAW 10/17
WINDOW wnHeight=10; wnWidth=5 wnGt=GT wnSMSO=0.; wnSMSC=0.; wnU=(1./((1./WinU)-0.85)) wnExHcMult=1.5 wnExHcModel=INPUT // BAW 10/17

SURFACE "East" sfType=Wall sfCon=WCav13 sfArea=100 sfAzm=Back sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExH=1.5
SURFACE "West" sfType=Wall sfCon=WCav13 sfArea=100 sfAzm=right sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExH=1.5 //
SURFACE "Ceil" sfType=Ceiling sfCon=WCav13 sfArea=100 sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExH=1.5 //
SURFACE "Floor" sfType=Floor sfCon=WCav13 sfArea=100 sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExH=1.5 //
SURFACE "East" sfType=Wall sfCon=WCav13 sfArea=100 sfAzm=Back sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExHcMult=1.5 sfExHcModel=INPUT
SURFACE "West" sfType=Wall sfCon=WCav13 sfArea=100 sfAzm=right sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExHcMult=1.5 sfExHcModel=INPUT//
SURFACE "Ceil" sfType=Ceiling sfCon=WCav13 sfArea=100 sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExHcMult=1.5 sfExHcModel=INPUT//
SURFACE "Floor" sfType=Floor sfCon=WCav13 sfArea=100 sfExAbs=0. sfModel = forward_difference sfExCnd=SPECIFIEDT sfExT=70 sfExHcMult=1.5 sfExHcModel=INPUT//

// SURFACE "IExSlab" UseType "IntExpSlab" sfArea=100

Expand Down
14 changes: 7 additions & 7 deletions test/actest1bL.cse
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -99,7 +99,7 @@ MATERIAL "Crawl" matThk=1. matCond=1./6. matSpHt=.24 ma
MATERIAL "SoftWood" matCond=0.98/12 matSpHt=.39 matDens=35. // 2009 HOF 26.6 Douglas fir/larch
MATERIAL "Concrete" matCond=1. matSpHt=0.2 matDens=144.
//MATERIAL "Sheath" matThk=1./12. matCond=1./(12*Rval) matSpHt=.35 matDens=1.5 // R4 from 2009 HOF 26.6
//MATERIAL "BlownInsul" matCond=0.032051282/RCQ matSpHt=.2 matDens=0.5 // Certainteed InsulSafe�SP Fiber Glass R=2.6/in for R30-38
//MATERIAL "BlownInsul" matCond=0.032051282/RCQ matSpHt=.2 matDens=0.5 // Certainteed InsulSafe�SP Fiber Glass R=2.6/in for R30-38
MATERIAL "BlownInsul" matCond=0.032051282/RCQ matSpHt=.2 matDens=0.5 matCondCT =0.00418 // variable k


Expand Down Expand Up @@ -673,13 +673,13 @@ Zone Single znArea=1000 infAC=0 infShld=4 infStories=1 infELA=0 znVol = 9000 znE
znRSYS = MySys //BAW 120312
//znQMxHRated = 1000000 znQMxCRated = -1000000 znQMxH =1000000 znQMxC = -1000000
znModel = CZM znTH = 70 znTD = 74 znTC = 78
SURFACE Wall1 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0.
SURFACE Wall2 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0.
SURFACE Wall3 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0.
SURFACE Wall4 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0.
SURFACE Wall1 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0. sfExHcModel = INPUT sfExHcMult = 2.64
SURFACE Wall2 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0. sfExHcModel = INPUT sfExHcMult = 2.64
SURFACE Wall3 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0. sfExHcModel = INPUT sfExHcMult = 2.64
SURFACE Wall4 sfType=Wall sfArea=300 sfCon=CondSurf; sfAZM=0 sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0. sfExHcModel = INPUT sfExHcMult = 2.64

SURFACE Floor sfType=floor sfArea=900 sfCon=CondSurf; sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0.
SURFACE Ceil sfType=ceiling sfArea=900 sfCon=CondSurf; sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0.
SURFACE Floor sfType=floor sfArea=900 sfCon=CondSurf; sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0. sfExHcModel = INPUT sfExHcMult = 2.64
SURFACE Ceil sfType=ceiling sfArea=900 sfCon=CondSurf; sfExCnd=SPECIFIEDT sfModel = FD sfLThkF=LFact sfInHcModel=Input sfExAbs=0 sfExT=OutT sfexEpsLW=0. sfExHcModel = INPUT sfExHcMult = 2.64
IZXFER CZtran izNVTYPE = AirNetExt izZN1 = Single izALo = 1 izCpr=0. izExp=.65 izHD =1

Run
Loading