mcnp从入门到精通讲座(7)

message: datapath=D:\\transport\\data\\mcnpxs

title -- Problem 1: cylinder along z axis using surfaces --- c

c ************************************************************ c

c !!! REMEMBER TO CHANGE THE MESSAGE BLOCK FOR YOUR SET-UP !!! c

c ************************************************************ c

c *** cell descriptions *** c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c density in g/cm^3 (-ve) c

20 0 +1:+5:-10 imp:n=0 $ outside world c

1 cz 5 $ infinite cyl along z-axis rad 5 cm c

5 pz +5 $ top of cylinder at 5cm

10 pz -5 $ bottom of cylinder at -5 cm c

c *** Data Cards after blank line***

c *** Source Card*** c

sdef $ default source: 14 MeV point neutron at (0,0,0) c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 using weight fractions (-ve) c

c ** Print first 50 histories print table 110 *** print 110 $ add table 110 to minimum print-out c

c *** time and nps cut cards ***

nps 60 $ nps 60 - stop simulation after 60 source particles ctme 1 $ ctme 1 - stop simulation after 1 min (60 sec) c

c *** Blank line terminator ***

End of File

例二：

message: datapath=D:\\transport\\data\\mcnpxs

title -- Ex 1b cylinder along z axis, in water using surfaces --- c c

c ************************************************************ c

c !!! REMEMBER TO CHANGE THE MESSAGE BLOCK FOR YOUR SET-UP !!! c

c ************************************************************ c

c *** cell descriptions *** c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c

c Use of the not operator - use with a cell ie #10 means not cell 10 c 15 2 -1 -15 +20 -25 +30 -35 +40 #10 imp:n=1 $ box but not cylinder c

c Expanded definition use of : operator, use of parenthesis and continuation

15 2 -1 (+1 -15 20 -25 30 -35 40):(-15 20 -25 30 -35 5) :(-15 20 -25 30 40 -10)

imp:n=1 $box but not cylinder

20 0 +15:-20:+25:-30:+35:-40 imp:n=0 $ outside world c

1 cz 5 $ infinite cyl along z-axis rad 5 cm c

5 pz +5 $ top of cylinder at 5 cm

10 pz -5 $ bottom of cylinder at -5 cm c

15 px +10 $ +ve x extent of warer 20 px -10 $ -ve x extent of water 25 py +10 $ +ve y extent of water 30 py -10 $ -ve y extent of water 35 pz +10 $ +ve z extent of water 40 pz -10 $ -ve z extent of water c

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238

m2 1001 +1 8016 +2 $ water (using atome fractions (+ve)) c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File 例三

message: datapath=D:\\transport\\data\\mcnpxs

title --Ex 1c cylinder in z-y plane rot. about x in water using surfaces -- c

c ************************************************************ c

c !!! REMEMBER TO CHANGE THE MESSAGE BLOCK FOR YOUR SET-UP !!! c

c ************************************************************ c

c *** cell descriptions *** c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c

c Use of the not operator - use with a cell ie #10 means not cell 10 c 15 2 -1 -15 +20 -25 +30 -35 +40 #10 imp:n=1 $ box but not cylinder c

c Expanded definition use of : operator, use of parenthesis and continuation

15 2 -1 (+1 -15 20 -25 30 -35 40):(-15 20 -25 30 -35 5) :(-15 20 -25 30 40 -10)

imp:n=1 $box but not cylinder

20 0 +15:-20:+25:-30:+35:-40 imp:n=0 $ outside world c

1 1 cz 5 $ infinite cyl along z-axis rad 5 cm

c

5 1 pz +5 $ top of cylinder at 5 cm

10 1 pz -5 $ bottom of cylinder at -5 cm c

15 px +10 $ +ve x extent of water 20 px -10 $ -ve x extent of water 25 py +10 $ +ve y extent of water 30 py -10 $ -ve y extent of water 35 pz +10 $ +ve z extent of water 40 pz -10 $ -ve z extent of water c

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 m2 1001 +1 8016 +2 $ water c

c *** Rotation matrix for cylinder *** c

c o1 o2 o3 b1 b2 b3 b4 b5 b6 b7 b8 b9 M c x y z xx' yx' zx' xy' yy' zy' xz' yz' zz' +1|-1

c Tr1 0 0 0 1 0 0 0 .7071 -.7071 0 .7071 .7071 1 $ cosines of angles *tr1 0 0 0 0 90 90 90 315 135 90 45 45 1 $ angles in degrees c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File

例四：

message: datapath=D:\\transport\\data\\mcnpxs

title -- cylinder along z axis, in box using surfaces --- c c

c *** cell descriptions ***

c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c

c Use of the not operator - use with a cell ie #10 means not cell 10 c 15 2 -1 -15 +20 -25 +30 -35 +40 #10 imp:n=1 $ H2O but not cylinder c

c Expanded deffinition use of : opperator, use of paranthesis and continuation

15 2 -1 (+1 -15 20 -25 30 -35 40):(-15 20 -25 30 -35 5) :(-15 20 -25 30 40 -10)

imp:n=1 $ H2O but not cylinder c

20 3 -2.7 -45 +50 +60 -55 -65 +35 imp:n=1 $ Al top 25 3 -2.7 -45 +50 +60 -55 -40 +70 imp:n=1 $ Al bottom 30 3 -2.7 -15 +20 +60 -30 -35 +40 imp:n=1 $ -y side 35 3 -2.7 -15 +20 +25 -55 -35 +40 imp:n=1 $ +y side 40 3 -2.7 -45 +15 +60 -55 -35 +40 imp:n=1 $ +x side 45 3 -2.7 -20 +50 +60 -55 -35 +40 imp:n=1 $ -x side 50 0 +45:-50:+55:-60:+65:-70 imp:n=0 $ outside world c

1 cz 5 $ infinite cyl along z-axis rad 5 cm c

5 pz +5 $ top of cylinder at 5 cm

10 pz -5 $ bottom of cylinder at -5 cm c

15 px +10 $ +ve x extent of warer 20 px -10 $ -ve x extent of water 25 py +10 $ +ve y extent of water 30 py -10 $ -ve y extent of water 35 pz +10 $ +ve z extent of water 40 pz -10 $ -ve z extent of water 45 px +12 $ +ve x extent of Al 50 px -12 $ -ve x extent of Al 55 py +12 $ +ve y extent of Al 60 py -12 $ -ve y extent of Al 65 pz +12 $ +ve z extent of Al 70 pz -12 $ -ve z extent of Al c

c *** Source Card*** c

sdef $ default source c

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库mcnp从入门到精通讲座(7)在线全文阅读。