introduction to helios t. k. kim argonne national laboratory
TRANSCRIPT
Introduction to HELIOSIntroduction to HELIOS
T. K. KimT. K. Kim
Argonne National LaboratoryArgonne National Laboratory
2Argonne National LaboratoryArgonne National Laboratory
Overview of HELIOSOverview of HELIOS
• Expert inputExpert input– Described by parametersDescribed by parameters– Not sensitive to geometry or composition dataNot sensitive to geometry or composition data
• Short inputShort input– Assign value to parametersAssign value to parameters– Sensitive to geometry and composition dataSensitive to geometry and composition data
HELIOSHELIOSGeneral 2D geometryGeneral 2D geometry
subgroup resonance methodsubgroup resonance methodcurrent coupling collision probability methodcurrent coupling collision probability method
HERMES file (data base file)
AURORAAURORAInput processorInput processor
ZENTIHZENTIHOutput processorOutput processor
3Argonne National LaboratoryArgonne National Laboratory
Flow Diagram of AURORA-HELIOS-Flow Diagram of AURORA-HELIOS-ZENITHZENITH
AURORAAURORAExpert inputExpert input
AURORAAURORA
AURORA SETsAURORA SETs(HERMES file)(HERMES file)
AURORAAURORAShort inputShort input
AURORAAURORA
AURORA SETsAURORA SETs(HERMES file)(HERMES file)
HELIOSHELIOS
ZENITHZENITHExpert inputExpert input
AURORAAURORA
ZENITH SETsZENITH SETs(HERMES file)(HERMES file)
ZENTIHZENTIH
ORIONORION
OutputOutput(ASCII)(ASCII)
PicturePicture(postscript)(postscript)
createcreate
createcreate
updateupdate
updateupdate
createcreate
createcreate
4Argonne National LaboratoryArgonne National Laboratory
MethodologiesMethodologies
• CPM in structure CPM in structure
– STR consists of region and circular geometry STR consists of region and circular geometry (CCS)(CCS)
– Region is defined by nodesRegion is defined by nodes
• CCPM between structures (CNX)CCPM between structures (CNX)
– Assign current coupling orderAssign current coupling order
• Boundary condition (BDRY)Boundary condition (BDRY)
– Reflective or Albedo conditions Reflective or Albedo conditions
5Argonne National LaboratoryArgonne National Laboratory
BWR assemblyBWR assembly
6Argonne National LaboratoryArgonne National Laboratory
AURORA Expert input - 1 AURORA Expert input - 1
• Create geometry data set Create geometry data set &ADD = SET('SBWR200.set'/AURORA;SBWR200; &Geometry)&ADD = SET('SBWR200.set'/AURORA;SBWR200; &Geometry)
$p = PAR("$FuelPitch")$p = PAR("$FuelPitch")
$rf = PAR("$FuelRadius")$rf = PAR("$FuelRadius")
$gw = PAR("$InnerGapThickness")$gw = PAR("$InnerGapThickness")
$DREF = PAR("$DENH-$AverageVoidFraction*($DENH-$DSSH)") $DREF = PAR("$DENH-$AverageVoidFraction*($DENH-$DSSH)") ! reference density!! reference density!
$DV000 = PAR("$DENH-0.00*($DENH-$DSSH)") $DV000 = PAR("$DENH-0.00*($DENH-$DSSH)") ! 0% void Density [g/cc]!! 0% void Density [g/cc]!
'COO-REF' = MAT(NB / 'COO-REF' = MAT(NB / ! Reference void + reference ppm ! ! Reference void + reference ppm !
$DREF/1001, 11.19 ; 8016, 88.81 ; 5000, $shrf) $DREF/1001, 11.19 ; 8016, 88.81 ; 5000, $shrf)
Fpin = CCS("0.90*$rf", $rf, $rc // Fuel1,Fuel2, Clad )Fpin = CCS("0.90*$rf", $rf, $rc // Fuel1,Fuel2, Clad )
$Fcell1 = PAR(("-$p2","-$p2") ("-$p2",$p2) ($p2,$p2) ($p2,"-$p2") $Fcell1 = PAR(("-$p2","-$p2") ("-$p2",$p2) ($p2,$p2) ($p2,"-$p2") ! 1- 4 !! 1- 4 !
("-$p2",0) (0,$p2) ($p2,0) (0,"-$p2") ("-$p2",0) (0,$p2) ($p2,0) (0,"-$p2") ! 5- 8 !! 5- 8 !
("-$rc",0) (0,$rc) ($rc,0) (0,"-$rc") ("-$rc",0) (0,$rc) ($rc,0) (0,"-$rc") ! 9-12 !! 9-12 !
/ 4,Cool /Fpin(0,0)/ 1,5,9,12,8,Cool; 5,2,6,10,9,Cool; 6,3,7,11,10,Cool)/ 4,Cool /Fpin(0,0)/ 1,5,9,12,8,Cool; 5,2,6,10,9,Cool; 6,3,7,11,10,Cool)
FUELPART = CNX($PinLoading/FUELPART = CNX($PinLoading/
( 1,3,4)$k( 2,2,1) / ( 2,3,4)$k( 3,2,1) / ( 1,3,4)$k( 2,2,1) / ( 2,3,4)$k( 3,2,1) / ! 1st row !! 1st row !
( 3,3,4)$k( 4,2,1) / ( 4,3,4)$k( 5,2,1) / ( 3,3,4)$k( 4,2,1) / ( 4,3,4)$k( 5,2,1) /
( 5,3,4)$k( 6,2,1) / ( 6,3,4)$k( 7,2,1) /( 5,3,4)$k( 6,2,1) / ( 6,3,4)$k( 7,2,1) /
( 7,3,4)$k( 8,2,1) /( 7,3,4)$k( 8,2,1) /
…………… ……………..
System = CNX(CROD,FUELPART,WALL/ System = CNX(CROD,FUELPART,WALL/
( 1-2,1)( 1-2,4)$k( 2-1,2)( 2-1,3) /( 1-2,1)( 1-2,4)$k( 2-1,2)( 2-1,3) /
( 1-7,3)( 1-7,4)$k( 3-1,2)( 3-1,1) )( 1-7,3)( 1-7,4)$k( 3-1,2)( 3-1,1) )
System = BDRY((1-1,2,2)1(0)) ! Specular boundary condition !System = BDRY((1-1,2,2)1(0)) ! Specular boundary condition !
7Argonne National LaboratoryArgonne National Laboratory
AURORA Expert input - 2 AURORA Expert input - 2
• Create overlay data set Create overlay data set &ADD = SET('SBWR200.set'/ AURORA;SBWR200; &Overlays) &ADD = SET('SBWR200.set'/ AURORA;SBWR200; &Overlays)
ovldNormal ovldNormal = OVLD(1/*-*-** ) = OVLD(1/*-*-** )
ovltHot ovltHot = OVLT( = OVLT( $TMOH /*-*-** / $TMOH /*-*-** /
$TCLH /*-*-*-Clad /$TCLH /*-*-*-Clad /
$TFUH /*-*-*-(Fuel1,Fuel2))$TFUH /*-*-*-(Fuel1,Fuel2))
ovlmFuel ovlmFuel = OVLM($FuelComp) ! Fuel material overlay != OVLM($FuelComp) ! Fuel material overlay !
ovsdBase ovsdBase = OVSD(ovldNormal)= OVSD(ovldNormal)
ovstBase ovstBase = OVST(ovltHot) = OVST(ovltHot)
ovsmBase ovsmBase = OVSM(ovlmNonFuel, ovlmFuel ) = OVSM(ovlmNonFuel, ovlmFuel )
statBase statBase = STAT(ovsmBase , ovsdBase, ovstBase , $HFPpd) = STAT(ovsmBase , ovsdBase, ovstBase , $HFPpd)
pathBase pathBase = PATH(/CG,(statBase),-50,-150,-500,-1000,-11000/10,-12500,= PATH(/CG,(statBase),-50,-150,-500,-1000,-11000/10,-12500,
-15000,-60000/18 )-15000,-60000/18 )
8Argonne National LaboratoryArgonne National Laboratory
AURORA Expert input - 3 AURORA Expert input - 3
• Create output data set Create output data set &ADD = SET('SBWR200.set'/ AURORA;SBWR200; &Output) &ADD = SET('SBWR200.set'/ AURORA;SBWR200; &Output)
ng1 = GROUP(N/0 ) ! one group !ng1 = GROUP(N/0 ) ! one group !
AllArea = AREA( <*-*-**> )AllArea = AREA( <*-*-**> )
AllFuel = AREA( <*-*-*-(Fuel1,Fuel2)>)AllFuel = AREA( <*-*-*-(Fuel1,Fuel2)>)
XSFuelMap XSFuelMap = MACRO(ng1, FuelMap / bu,kf,ab,fi)= MACRO(ng1, FuelMap / bu,kf,ab,fi)
XSMicAll XSMicAll = MICRO(ngs, AllArea / 1,53635,54635,62649/ab)= MICRO(ngs, AllArea / 1,53635,54635,62649/ab)
ESide ESide = FACE((3- 3,3,4)(3- 4,3,4)(3- 5,3,4)(3- 6,3,4)(3- 7,3,4)= FACE((3- 3,3,4)(3- 4,3,4)(3- 5,3,4)(3- 6,3,4)(3- 7,3,4)
(3- 8,3,4)(3- 9,3,4)(3-11,3,4)(3-13,3,4)(3-23,3,4))(3- 8,3,4)(3- 9,3,4)(3-11,3,4)(3-13,3,4)(3-23,3,4))
NECorner NECorner = FACE((3- 3,2,4))= FACE((3- 3,2,4))
SideCurrents = CUR(ngs,ESide ,NSide ,WSide ,SSide )SideCurrents = CUR(ngs,ESide ,NSide ,WSide ,SSide )
CornerCurrents = CUR(ngs,NECorner,NWCorner,SECorner,SWCorner)CornerCurrents = CUR(ngs,NECorner,NWCorner,SECorner,SWCorner)
9Argonne National LaboratoryArgonne National Laboratory
Aurora Short InputAurora Short Input
• Call library file and create HERME fileCall library file and create HERME file'SBWR200' = CASE('C:\TKKIM\HELIOS\Library\xslib035-1.5' /'SBWR200.hrf'/'SBWR 200 Short 'SBWR200' = CASE('C:\TKKIM\HELIOS\Library\xslib035-1.5' /'SBWR200.hrf'/'SBWR 200 Short
Input')Input')
• Assign values to parametersAssign values to parameters&Geometry = SET('SBWR200.set'/AURORA;SBWR200)&Geometry = SET('SBWR200.set'/AURORA;SBWR200)
&Overlays = SET('SBWR200.set'/AURORA;SBWR200)&Overlays = SET('SBWR200.set'/AURORA;SBWR200)
&Output = SET('SBWR200.set'/AURORA;SBWR200)&Output = SET('SBWR200.set'/AURORA;SBWR200)
$k = PAR(4) $k = PAR(4) ! Current coupling !! Current coupling !
$HFPpd = PAR("21.00") $HFPpd = PAR("21.00") ! HFP Power density [W/gU] ! ! HFP Power density [W/gU] !
$FuelPitch = PAR("1.6200") $FuelPitch = PAR("1.6200") ! Fuel pitch [cm] !! Fuel pitch [cm] !
$FuelRadius = PAR("0.5220") $FuelRadius = PAR("0.5220") ! Fuel rdaius [cm] !! Fuel rdaius [cm] !
'SBWR200' = RUN()'SBWR200' = RUN()
10Argonne National LaboratoryArgonne National Laboratory
Zenith Expert InputZenith Expert Input
• Call library file and create HERME fileCall library file and create HERME file
&ADD = SET('SBWR200.set'/ ZENITH; &StateData) &ADD = SET('SBWR200.set'/ ZENITH; &StateData)
KINF KINF = SEL(kinf /STATUS/C1;/$calp)= SEL(kinf /STATUS/C1;/$calp)
KINFB KINFB = SEL(kinfb/STATUS/C1;/$calp)= SEL(kinfb/STATUS/C1;/$calp)
TotalVolume TotalVolume = SEL(vo /MICRO/C1;XSMicAll )= SEL(vo /MICRO/C1;XSMicAll )
dr dr = SEL(dr /MACRO/C1;XSMacAll /$calp )= SEL(dr /MACRO/C1;XSMacAll /$calp )
tr tr = SEL(tr /MACRO/C1;XSMacAll /$calp )= SEL(tr /MACRO/C1;XSMacAll /$calp )
jps jps = SEL(jp /CUR /C1;SideCurrents /$calp ) = SEL(jp /CUR /C1;SideCurrents /$calp )
jms jms = SEL(jm /CUR /C1;SideCurrents /$calp ) = SEL(jm /CUR /C1;SideCurrents /$calp )
sideADF sideADF = FOR(2*(jps+jms)/fx)= FOR(2*(jps+jms)/fx)
DiffusionCoefficient DiffusionCoefficient = FOR(dr/(3*tr))= FOR(dr/(3*tr))
TransXs TransXs = FOR(1/(3*DiffusionCoefficient))= FOR(1/(3*DiffusionCoefficient))
11Argonne National LaboratoryArgonne National Laboratory
Zenith Short InputZenith Short Input
• Call library file and create HERME fileCall library file and create HERME file
&ADD = SET('SBWR200.set'/ ZENITH; &StateData) &ADD = SET('SBWR200.set'/ ZENITH; &StateData)
KINF KINF = SEL(kinf /STATUS/C1;/$calp)= SEL(kinf /STATUS/C1;/$calp)
KINFB KINFB = SEL(kinfb/STATUS/C1;/$calp)= SEL(kinfb/STATUS/C1;/$calp)
TotalVolume TotalVolume = SEL(vo /MICRO/C1;XSMicAll )= SEL(vo /MICRO/C1;XSMicAll )
dr dr = SEL(dr /MACRO/C1;XSMacAll /$calp )= SEL(dr /MACRO/C1;XSMacAll /$calp )
tr tr = SEL(tr /MACRO/C1;XSMacAll /$calp )= SEL(tr /MACRO/C1;XSMacAll /$calp )
jps jps = SEL(jp /CUR /C1;SideCurrents /$calp ) = SEL(jp /CUR /C1;SideCurrents /$calp )
jms jms = SEL(jm /CUR /C1;SideCurrents /$calp ) = SEL(jm /CUR /C1;SideCurrents /$calp )
sideADF sideADF = FOR(2*(jps+jms)/fx)= FOR(2*(jps+jms)/fx)
DiffusionCoefficient DiffusionCoefficient = FOR(dr/(3*tr))= FOR(dr/(3*tr))
TransXs TransXs = FOR(1/(3*DiffusionCoefficient))= FOR(1/(3*DiffusionCoefficient))
12Argonne National LaboratoryArgonne National Laboratory
Zenith Expert InputZenith Expert Input
BEGIN ('SBWR200 Fuel Calculation by T.K.Kim'/0.5;0.5/LR;BT)BEGIN ('SBWR200 Fuel Calculation by T.K.Kim'/0.5;0.5/LR;BT)
%OPTIONS(LINPAG=500, MXAUX=130000, %OPTIONS(LINPAG=500, MXAUX=130000, MXLNAM=10000,MXCON=30000,MXSET=1500)MXLNAM=10000,MXCON=30000,MXSET=1500)
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
! <<< Base Depletion Case >>>! <<< Base Depletion Case >>>
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
%'Case_Switch' = 0 %'Case_Switch' = 0 ! Base depletion edit flag.! Base depletion edit flag.
%'Form_Function' = 1 %'Form_Function' = 1 ! Form Function edit flag, 0/1=no/yes ! Form Function edit flag, 0/1=no/yes
%'Power_Map' = 1 %'Power_Map' = 1 ! Power map edit flag, 0/1=no/yes ! Power map edit flag, 0/1=no/yes
%'Fluxr_Map' = 0 %'Fluxr_Map' = 0 ! Flux map edit flag, 0/1=no/yes ! Flux map edit flag, 0/1=no/yes
%'Thorium_Fuel' = 0 %'Thorium_Fuel' = 0 ! 0/1 =non Thorium fuel / Thorium fuel ! 0/1 =non Thorium fuel / Thorium fuel
%'HM_number_density' = 1 %'HM_number_density' = 1 ! Heavy meatl number ddensity falg 0/1=no/yes! Heavy meatl number ddensity falg 0/1=no/yes
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
! -------------------------------------------------------------- Run Zenith --! -------------------------------------------------------------- Run Zenith --
! ---------------------------------------------------------------------------- ! ----------------------------------------------------------------------------
$Hfile = PAR('SBWR200.hrf') ! Hermes format file name $Hfile = PAR('SBWR200.hrf') ! Hermes format file name
$Hcase = PAR('SBWR200') ! Case name $Hcase = PAR('SBWR200') ! Case name
C1 = IMP(HELIOS;$Hcase/$Hfile) ! HELIOS Import pointC1 = IMP(HELIOS;$Hcase/$Hfile) ! HELIOS Import point
&RunZenith = SET('SBWR200.set'/ ZENITH) &RunZenith = SET('SBWR200.set'/ ZENITH)
END()END()