Scribd
Upload
93 views16 pages

NHT Assignment

This program calculates pressure drop along bare nuclear fuel bundles and bundles with spacer grids using correlations from various sources. It takes user input for the geometry, material properties, and flow conditions of the bundle. For bare bundles, it calculates pressure drop based on correlations for square or hexagonal geometry. For bundles with spacer grids, it adds calculations of grid parameters and uses different correlations to account for the grids.

Uploaded by

Engr Saif Ullah
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
93 views16 pages

NHT Assignment

This program calculates pressure drop along bare nuclear fuel bundles and bundles with spacer grids using correlations from various sources. It takes user input for the geometry, material properties, and flow conditions of the bundle. For bare bundles, it calculates pressure drop based on correlations for square or hexagonal geometry. For bundles with spacer grids, it adds calculations of grid parameters and uses different correlations to account for the grids.

Uploaded by

Engr Saif Ullah
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

ASSIGNMENT 01

TO DEVELOP FORTRAN PROGRAM FOR CALCULATING PRESSURE DROP

SUBJECT: NUCLEAR HEAT TRANSPORT

SUBMITTED BY: SAIF ULLAH

SUBMITTED TO: DR WASEEM SIDDIQUE

CLASS # 18

DEPARTMENT OF NUCLEAR ENGINEERING


PAKISTAN INSTITUTE OF ENGINEERING AND APPLIED SCIENCES
Dated 13 Nov, 2020
(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

C THIS PROGRAM HAS BEEN DEVELOPED TO FIND THE PRESSURE DROP


C ALONG BARE BUNDLES AND WITH SPACER GRIDS IN BOTH SQUARE AND
C HEXAGONAL GEOMETRY

PROGRAM PRESSURE_DROP_LN_ET_AL_and_STORDEUR_AND REHME_CORRELATIONS

C DECLARATION OF VARIABLES

real GST,GSH,Dia,PITCH,rho,MEU,Deq,DeqG,Vv,Vg,mdot,gg,z,g,k_rehme
real area_grid_frontal,A_subch,A_sub_grid,Re_Sub,Re_grid,p,n

real PD_Ratio,PD_1,PD_2,PD_3,ai,b1i,b2i,ae,b1e,b2e,ac,b1c,b2c,FcL
real fel,cefil,fil,cifil,W,BFA,Flowarea_C,Flowarea_E,fnew

real c_form_grid,cfg,c_fric_rod,Area_grid_wet,cfgl,cfgt,k_,y
real form_a,form_b,form_c,area_rod_Wetted_grid,eps,Re_grid_stord

real Wetted_sum,Wetted_i,Wetted_E,Wetted_C,DEB,WD_Ratio,L
real Dei,Dee,Dec,sumc,sume,sumi,cbl,f,pdrop,Flowarea_I,CcFIL

character*10 bundle, bare,square

Print*,"To Find pressure drop along bare bundle(type bare)"


Print*,"To Find pressure drop in bundle with spacer grid(type spacer)"
read*, bare
if ([Link].'bare') then ! Condition for Bare or spacer grid bundle

c input data for square and hexagonal geometry bare rod bundles"
Print*, "Enter the Diameter of rod(m)"
read*,Dia
Print*,"Enter the Pitch(m)"
read*,Pitch
Print*,"Enter the Length(m)"
read*,l
Print*,"Enter the width(rod+clearance)(m)"
read*,w
Print*,"Enter the rho(SI units) based on average temp"
read*,rho
Print*,"Enter the dynamic viscosity(SI units) based on average temp"
read*,meu
Print*,"Enter the Mass flow rate(Kg/sec)"
read*,mdot
Print*,"Enter the Number of assemblies"
read*,p
Print*,"Enter the number of pins in one array"
read*,n
c Calculations for square geometry bare rod bundles
A_subch= Pitch**2-0.25*3.14*Dia**2
Deq= 4*A_subch/(3.14*Dia)
Q=((MDOT/(p*N**2)))/rho
print*, "Q",q
VV=Q/A_subch

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

print*, "vv", vv
Re_Sub=(rho*VV*Deq)/MEU
Print*, "Reynold Number", Re_Sub
PD_Ratio=Pitch/Dia
WD_Ratio=W/Dia
Print*,"PD_Ratio",PD_Ratio
if (PD_Ratio.lt.1.0) then
Print*,"Correlation Not Applicable"
go to 7
else
PD_1=1.0
pD_2=1.1
PD_3=1.5
Print*," Select one option"
Print*, "Enter sq for square or Hex for hexagonal geometry"
read*,square
if ([Link]."sq") then ! condition for square geometry
if(Re_Sub.LE.2300) then ! condition for laminar flow
Print*,"LAMINAR FLOW"
if(PD_Ratio.LE.PD_2.AND.PD_Ratio.GE.PD_1) then
Print*,"Pressure drop for square Geometry Calculations"
ai=26.37
b1i=374.2
b2i=-493.9
ae=26.18
b1e=554.5
b2e=-1480
ac=28.62
b1c=715.9
b2c=-2807
c For interior ===============================================
CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2
FiL=CiFIL/(Re_Sub)**0.12
c For edge ===============================================
CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2
FeL=CeFIL/(Re_Sub)**0.12
c For corner ===============================================
CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2
FcL=CcFIL/(Re_Sub)**0.12

Flowarea_I=((N-1)**2)*((Pitch**2)-((3.14/4)*Dia**2))
Flowarea_E=4*(N-1)*((W-(Dia/2))*Pitch-((3.14/8)*Dia**2))
Flowarea_C=4*((W-(Dia/2))**2-((3.14/16)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=((N-1)**2)*3.14*Dia
Wetted_E=4*((N-1)*(3.14/2)*Dia+Pitch)
Wetted_C=4*((3.14/4)*Dia+(2*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

SUMI=((Flowarea_I)/BFA)*((DEI/DEB))*(CiFIL/DEI)**(-1)
SUME=((Flowarea_E)/BFA)*((DEE/DEB))*(CeFIL/DEE)**(-1)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB))*(CcFIL/DEC)**(-1)

CBL=DEB*(SUMC+SUME+SUMI)**(-1.0)

F=CBL/Re_Sub

PDROP=(F*(L/Deq)*(0.5/9.81)*(VV**2))/100
Print*,"PDROP ACROSS 1 ASSEMBLY=",PDROP*n**2
Print*,"PDROP ACROSS CORE =",PDROP*p*n**2

ELSE IF(PD_Ratio.GT.PD_2.AND.PD_Ratio.LE.PD_3) THEN


ai=35.55
b1i=263.7
b2i=-190.2
ae=44.40
b1e=256.7
b2e=-267.6
ac=58.83
b1c=160.7
b2c=-203.5

CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2
CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2
CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2
Flowarea_I=((N-1)**2)*((P**2)-((3.14/4)*Dia**2))
Flowarea_E=4*(N-1)*((W-(Dia/2))*P-(3.14/8)*Dia**2)
Flowarea_C=4*((W-(Dia/2))**2-((3.14/16)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=((N-1)**2)*3.14*Dia
Wetted_E=4*(N-1)*((3.14/2)*Dia+Pitch)
Wetted_C=4*((3.14/4)*Dia+(2*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

SUMI=((Flowarea_I)/BFA)*((DEI/DEB))*(CiFIL/DEI)**(-1)
SUME=((Flowarea_E)/BFA)*((DEE/DEB))*(CeFIL/DEE)**(-1)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB))*(CcFIL/DEC)**(-1)

CBL=DEB*(SUMC+SUME+SUMI)**(-1.0)
F=CBL/(Re_Sub)**0.12
print*, "Friction Factor",f
PDROP=(F*(L/Deq)*(0.5/9.81)*(VV**2))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2


(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

END IF
ELSE
Print*,"turbulent case for square"
IF(PD_Ratio.LE.PD_2.AND.PD_Ratio.GE.PD_1) THEN
Print*,"square case"
ai=0.09423
b1i=0.5806
b2i=-1.239
ae=0.09377
b1e=0.8732
b2e=-3.341
ac=0.09755
b1c=1.127
b2c=-6.304
! center one========
CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2
CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2
CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2
Flowarea_I=((N-1)**2.0)*((Pitch**2)-((3.14/4)*Dia**2))
Flowarea_E=4*(N-1)*(((W-(Dia/2))*Pitch)-(3.14/8)*Dia**2)
Flowarea_C=4*((W-(Dia/2.0))**2-((3.14/16.0)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=((N-1)**2)*3.14*Dia
Wetted_E=4*(N-1)*((3.14/2)*Dia+Pitch)
Wetted_C=4*((3.14/4)*Dia+(2*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

SUMI=((Flowarea_I)/BFA)*((DEI/DEB)**0.11)*(CiFIL/DEI)**(-0.55)
SUME=((Flowarea_E)/BFA)*((DEE/DEB)**0.11)*(CeFIL/DEE)**(-0.55)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB)**0.11)*(CcFIL/DEC)**(-0.55)

CBL=DEB*((SUMC+SUME+SUMI)**(-1.8))
F=CBL/(Re_Sub)**0.25
print*,"f",f
PDROP=(F*(L/DEq)*(0.5/9.81)*(VV**2))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2

ELSE IF(PD_Ratio.GT.PD_2.AND.PD_Ratio.LE.PD_3) THEN


ai=0.1339
b1i=0.09059
b2i=-0.09926
ae=0.1430
b1e=0.04199
b2e=-0.04428
ac=0.1452
b1c=0.02681
b2c=-0.03411

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

Print*,"======================================"
C CENTRAL SUBCHANNEL======================================================
CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2
C EDGE SUBCHANNEL ========================================================
CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2
C CORNER SUBCHANNEL ========================================================
CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2
C EQUIVALENT DIA CALCULATIONS ==============================================
Flowarea_I=((N-1)**2)*((P**2)-((3.14/4)*Dia**2))
Flowarea_E=4*(N-1)*((W-(Dia/2))*P-(3.14/8)*Dia**2)
Flowarea_C=4*((W-(Dia/2))**2-((3.14/16)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=((N-1)**2)*3.14*Dia
Wetted_E=4*(N-1)*((3.14/2)*Dia+Pitch)
Wetted_C=4*((3.14/4)*Dia+(2*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

SUMI=((Flowarea_I)/BFA)*((DEI/DEB)**0.11)*(CiFIL/DEI)**(-0.55)
SUME=((Flowarea_E)/BFA)*((DEE/DEB)**0.11)*(CeFIL/DEE)**(-0.55)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB)**0.11)*(CcFIL/DEC)**(-0.55)

CBL=DEB*((SUMC+SUME+SUMI)**-1.8)

F=CBL/(Re_Sub)**0.25
print*,"f",f

PDROP=(F*(L/DEq)*(0.5/9.81)*(VV**2.0))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2
else
Print*,"Using MaAdams Correlation if P/D is greater than 1.5"
F=0.184*Re_Sub**-0.2
PDROP=(F*(L/DEq)*(0.5)*(9.81)*(VV**2))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2

END IF
END IF

c Hexagonal geometry Calculations started

else
Print*,"TRIANGULAR ARRAY with bare rod bundles"
IF(Re_Sub.LE.2300) THEN
Print*,"LAMINAR FLOW"
IF(PD_Ratio.LE.PD_2.AND.PD_Ratio.GE.PD_1) THEN
ai=26.00

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

b1i=888.2
b2i=-3334.9
ae=26.18
b1e=554.5
b2e=-1480
ac=26.98
b1c=1636
b2c=-10050
C CENTRAL SUBCHANNEL =======================================
CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2

C EDGE SUBCHANNEL ==========================================


CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2

C CORNER SUBCHANNEL =======================================


CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2

C EQUIVALENT DIA CALCULATIONS =============================


Flowarea_I=(6*(N)**2)*(((SQRT(3.0)*Pitch**2)/4.0)-((3.14/8)*Dia**2))
Flowarea_E=6*(N)*(((W-(Dia/2))*Pitch)-((3.14/8)*Dia**2))
Flowarea_C=6*(((1/SQRT(3.0))*((W-(Dia/2))**2))-((3.14/24)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=(6*(N)**2)*3.14*(Dia/2.0)
Wetted_E=6*(N)*(((3.14/2)*Dia)+Pitch)
Wetted_C=6*(((3.14/6)*Dia)+((2/SQRT(3.0))*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

SUMI=((Flowarea_I)/BFA)*((DEI/DEB))*(CiFIL/DEI)**(-1)
SUME=((Flowarea_E)/BFA)*((DEE/DEB))*(CeFIL/DEE)**(-1)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB))*(CcFIL/DEC)**(-1)

CBL=DEB*(SUMC+SUME+SUMI)**(-1.0)

F=CBL/(Re_Sub)**0.12

PDROP=(F*(L/Deq)*(0.5/9.81)*(VV**2))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2

ELSE IF(PD_Ratio.GT.PD_2.AND.PD_Ratio.LE.PD_3) THEN


ai=62.97
b1i=216.9
b2i=-190.2
ae=44.40
b1e=256.7
b2e=-267.6
ac=87.26

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

b1c=38.59
b2c=-55.12
C CENTRAL SUBCHANNEL ===========================================
CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2

C EDGE SUBCHANNEL ============================================


CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2

C CORNER SUBCHANNEL ===========================================


CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2

C EQUIVALENT DIA CALCULATIONS ==================================


Flowarea_I=(6*(N)**2)*(((SQRT(3.0)*Pitch**2)/4.0)-((3.14/8)*Dia**2))
Flowarea_E=6*(N)*(((W-(Dia/2))*Pitch)-((3.14/8)*Dia**2))
Flowarea_C=6*(((1/SQRT(3.0))*((W-(Dia/2))**2))-((3.14/24)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=(6*(N)**2)*3.14*(Dia/2.0)
Wetted_E=6*(N)*(((3.14/2)*Dia)+Pitch)
Wetted_C=6*(((3.14/6)*Dia)+((2/SQRT(3.0))*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

SUMI=((Flowarea_I)/BFA)*((DEI/DEB))*(CiFIL/DEI)**(-1)
SUME=((Flowarea_E)/BFA)*((DEE/DEB))*(CeFIL/DEE)**(-1)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB))*(CcFIL/DEC)**(-1)

CBL=DEB*(SUMC+SUME+SUMI)**(-1.0)

F=CBL/(Re_Sub)**0.12

PDROP=(F*(L/Deq)*(0.5/9.81)*(VV**2))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2

END IF
ELSE
Print*,"turbulent case for HEXAGON"
IF(PD_Ratio.LE.PD_2.AND.PD_Ratio.GE.PD_1) THEN
Print*,"HEXAGONE case"
ai=0.09378
b1i=1.398
b2i=-8.664
ae=0.09377
b1e=0.8732
b2e=-3.341
ac=0.1004
b1c=1.625
b2c=-11.85

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2
CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2
CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2
Flowarea_I=(6*(N)**2)*(((SQRT(3.0)*Pitch**2)/4.0)-((3.14/8)*Dia**2))
Flowarea_E=6*(N)*(((W-(Dia/2))*Pitch)-((3.14/8)*Dia**2))
Flowarea_C=6*(((1/SQRT(3.0))*((W-(Dia/2))**2))-((3.14/24)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=(6*(N)**2)*3.14*(Dia/2.0)
Wetted_E=6*(N)*(((3.14/2)*Dia)+Pitch)
Wetted_C=6*(((3.14/6)*Dia)+((2/SQRT(3.0))*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

SUMI=((Flowarea_I)/BFA)*((DEI/DEB))*(CiFIL/DEI)**(-1)
SUME=((Flowarea_E)/BFA)*((DEE/DEB))*(CeFIL/DEE)**(-1)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB))*(CcFIL/DEC)**(-1)

CBL=DEB*(SUMC+SUME+SUMI)**(-1.0)

F=CBL/(Re_Sub)**0.25

PDROP=(F*(L/Deq)*(0.5/9.81)*(VV**2))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2

ELSE IF(PD_Ratio.GT.PD_2.AND.PD_Ratio.LE.PD_3) THEN


ai=0.1458
b1i=0.03632
b2i=-0.03333
ae=0.1430
b1e=0.04199
b2e=-0.04428
ac=0.1499
b1c=0.006706
b2c=-0.009567
Print*,"======================================"
CiFIL=ai+b1i*(PD_Ratio-1)+b2i*(PD_Ratio-1)**2
CeFIL=ae+b1e*(WD_Ratio-1)+b2e*(WD_Ratio-1)**2
CcFIL=ac+b1c*(WD_Ratio-1)+b2c*(WD_Ratio-1)**2
Flowarea_I=(6*(N)**2)*(((SQRT(3.0)*Pitch**2)/4.0)-((3.14/8)*Dia**2))
Flowarea_E=6*(N)*(((W-(Dia/2))*Pitch)-((3.14/8)*Dia**2))
Flowarea_C=6*(((1/SQRT(3.0))*((W-(Dia/2))**2))-((3.14/24)*Dia**2))
BFA=Flowarea_C+Flowarea_E+Flowarea_I
Wetted_I=(6*(N)**2)*3.14*(Dia/2.0)
Wetted_E=6*(N)*(((3.14/2)*Dia)+Pitch)
Wetted_C=6*(((3.14/6)*Dia)+((2/SQRT(3.0))*(W-(Dia/2))))
Wetted_SUM=Wetted_I+Wetted_E+Wetted_C


(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

DEB=(4*BFA)/Wetted_SUM

DEI=(4*Flowarea_I/Wetted_I)
DEE=(4*Flowarea_E/Wetted_E)
DEC=(4*Flowarea_C/Wetted_C)

SUMI=((Flowarea_I)/BFA)*((DEI/DEB))*(CiFIL/DEI)**(-1)
SUME=((Flowarea_E)/BFA)*((DEE/DEB))*(CeFIL/DEE)**(-1)
SUMC=((Flowarea_C)/BFA)*((DEC/DEB))*(CcFIL/DEC)**(-1)

CBL=DEB*(SUMC+SUME+SUMI)**(-1.0)

F=CBL/(Re_Sub)**0.25
PDROP=(F*(L/Deq)*(0.5/9.81)*(VV**2))/100
Print*,"Pressure drop across one assembly in KPa=",PDROP*n**2
Print*,"Pressure drop across core in KPa =",PDROP*p*n**2
END IF
END IF
end if
end if
else
Print*, " Type Sq for Square Geometry or Hex for Hexagonal Geometry"
read*, square
if ( [Link]."sq") then

Print*, "Enter the Grid strap thickness(m)"


read*,gst
Print*,"Enter the Grid Strap height(m)"
read*,gsh
Print*, "Enter the Diameter of rod(m)"
read*,Dia
Print*,"Enter the Pitch(m)"
read*,Pitch
Print*,"Enter the rho(SI units)"
read*,rho
Print*,"Enter the dynamic viscosity(SI units)"
read*,meu
Print*,"Enter the Mass flow rate(Kg/sec)"
read*,mdot
Print*,"Enter the Number of assemblies"
read*,p
Print*,"Enter the number of pins in one array"
read*,n
Print*,"Enter the number of Grid spacer"
read*,y
area_grid_frontal= pitch**2-(pitch-gst)**2 ! area grid frontal calculated
Print*,"area_grid_frontal=",area_grid_frontal
A_subch= pitch**2-0.25*3.14*Dia**2 ! area of subchannel away Grid
Print*, "areA_subchannel=",A_subch
A_sub_grid=(pitch-gst)**2-0.25*3.14*Dia**2 ! area of subchannel at Grid
Print*, "area of sub channel at grid =", A_sub_grid
Deq= 4*A_subch/(3.14*Dia) ! Eq Diameter subchannel away Grid
Print*, "Equivalent Dia at sub channel=",Deq
DeqG=4*A_sub_grid/(3.14*Dia+6*((pitch/2.0)-gst/sqrt(3.0)))!Equivalent Dia at Grid

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

Print*, "Equivalent Diameter at Grid=",DeqG


q=(mdot/(p*n**2))/rho
vv=q/A_subch
vg=q/A_sub_grid
Re_Sub=(rho*vv*Deq)/meu
Re_grid=(rho*vg*DeqG)/meu
Print*, "Reynold number away from grid=",Re_Sub
Print*, "Reynold number at grid=",Re_grid
G=(mdot/(p*n**2))/A_subch
Gg=(mdot/(p*n**2))/A_sub_grid
c_form_grid=2.75-0.27*log10(Re_Sub)
Print*, "c_form_grid=",c_form_grid
z=((3e4)*meu)/gg
cfgl= 1.328*((Gg*(gsh-((3e4)*meu/gg))/meu)**-0.5)
cfgt= 0.523*((alog((0.06*Gg*(gsh-((3e4)*meu/gg)))/meu))**-2.0)

if([Link]) then
cfg=(cfgl*(3e4)*(meu))/(gg*gsh)+(cfgt)*((gsh-((3e4*meu)/gg))/gsh)
else
cfg=(cfgl*(3e4)*(meu))/(gg*gsh)
end if
c_fric_rod=0.184*((Re_grid)**-0.2)
Area_grid_wet=4*(pitch-gst)*gsh
area_rod_Wetted_grid=3.14*Dia*gsh
eps=(area_grid_frontal/A_subch)

form_a=( c_form_grid*(eps/(1-eps)**2))

form_b=((cfg* Area_grid_wet)/(A_subch*((1-eps)**2)))

form_c=(c_fric_rod*(area_rod_Wetted_grid))/(A_subch*((1-eps)**2))
Print*,"======================================"
Print*,"Ln [Link] Correlation for pressure drop"

Print*,"K=", form_a+form_b+form_c
Print*,"Pressure drop", (y*(form_a+form_b+form_c)*g**2)/(2*rho)
Print*,"======================================"
Print*, "DE STORDEUR'S MODEL"
Re_grid_stord= (gst*gg/2*rho)
k_= 1.7*(area_grid_frontal*A_subch**2)/(A_subch*A_sub_grid**2)
Print*,k_
Print*, "Pressure drop =",(y*k_*g**2)/(2*rho)
Print*,"======================================"
Print*, "Rehme Correlation"
k_rehme= 6.5*(area_grid_frontal/A_subch)**2
Print*, "Pressure drop=", y*k_rehme*g**2/(2*rho)

else
Print*, "Enter the Grid strap thickness(m)"
read*,gst
Print*,"Enter the Grid Strap height(m)"
read*,gsh
Print*, "Enter the Diameter of rod(m)"
read*,Dia

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

Print*,"Enter the Pitch(m)"


read*,Pitch
Print*,"Enter the Density(SI units)"
read*,rho
Print*,"Enter the dynamic viscosity(SI units)"
read*,meu
Print*,"Enter the Mass flow rate(Kg/sec)"
read*,mdot
Print*,"Enter the Number of assemblies"
read*,p
Print*,"Enter the number of pins in one assembly"
read*,n
Print*,"Enter the number of Grid spacers"
read*,y
area_grid_frontal= (3.0/2.0)*(pitch*gst-(gst**2/sqrt(3.0)))
A_subch= (((Sqrt(3.0))/4.0)*pitch**2)-((3.14*dia**2)/8.0)

A_sub_grid= A_subch-area_grid_frontal

Deq= 4*A_subch/(3.14*Dia)
DeqG=4*A_sub_grid/(3.14*Dia+6*((pitch/2.0)-gst/sqrt(3.0)))
q=(mdot/(p*n))/rho
vv=q/A_subch
vg=q/A_sub_grid
Re_Sub=(rho*vv*Deq)/meu
Re_grid=(rho*vg*DeqG)/meu
G=(mdot/(p*n))/A_subch
Gg=(mdot/(p*n))/A_sub_grid
c_form_grid=2.75-0.27*log10(Re_Sub)
z=((3e4)*meu)/gg
cfgl= 1.328*((Gg*(gsh-((3e4)*meu/gg))/meu)**-0.5)
cfgt= 0.523*((alog((0.06*Gg*(gsh-((3e4)*meu/gg)))/meu))**-2.0)

if([Link]) then
cfg=(cfgl*(3e4)*(meu))/(gg*gsh)+(cfgt)*((gsh-((3e4*meu)/gg))/gsh)
else
cfg=(cfgl*(3e4)*(meu))/(gg*gsh)
end if
c_fric_rod=0.184*((Re_grid)**-0.2)
Area_grid_wet=6*((pitch/2)-(gst/(sqrt(3.0))))*gsh
area_rod_Wetted_grid=3.14*Dia*gsh
eps=(area_grid_frontal/A_subch)

form_a=( c_form_grid*(eps/(1-eps)**2))

form_b=((cfg* Area_grid_wet)/(A_subch*((1-eps)**2)))

form_c=(c_fric_rod*(area_rod_Wetted_grid))/(A_subch*((1-eps)**2))
Print*,"======================================"
Print*,"Ln [Link] Correlation for pressure drop"

Print*,"K=", form_a+form_b+form_c
Print*,"Pressure drop in KPa ", (y*(form_a+form_b+form_c)*g**2)/(2*rho)
Print*,"======================================"

(?18&/($5(1*,1((5,1*?=(526(0(67(5?&)?1HZIROGHU?6RIWZDUH?6RIWZDUHV?J?ELQ?QKWWIRU 7KXUVGD\1RYHPEHU30

Print*, "DE STORDEUR'S MODEL"


Re_grid_stord= (gst*gg/2*rho)
k_= 1.7*(area_grid_frontal*A_subch**2)/(A_subch*A_sub_grid**2)
Print*, "Pressure drop in KPa =",(y*k_*g**2)/(2*rho)
Print*,"======================================"
Print*, "Rehme Correlation"
k_rehme= 6.5*(area_grid_frontal/A_subch)**2
Print*, "Pressure drop in KPa=", y*k_rehme*g**2/(2*rho)
end if
end if
7 continue
end


Example 9.9 Run by Program. The results have been shown in Fig 1.

 
          Figure 1 Example 9.9 square geometry with spacer grids 

 
Problem 2: Pressure Drop in a test facility with Given Conditions for Hexagonal geometry with 
spacer Grids. 

Figure 2 Test Facility Pressure drop calculations by Using program 

 
Bare Rod bundle pressure drop having square geometry

 
Figure 3 Bare rod bundle with square geometry 

You might also like