12.07.2015 Views

Monte Carlo Particle Transport Methods: Neutron and Photon - gnssn

Monte Carlo Particle Transport Methods: Neutron and Photon - gnssn

Monte Carlo Particle Transport Methods: Neutron and Photon - gnssn

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

176 <strong>Monte</strong> <strong>Carlo</strong> <strong>Particle</strong> <strong>Transport</strong> <strong>Methods</strong>: <strong>Neutron</strong> <strong>and</strong> <strong>Photon</strong> CalculationsComparison of Equations (5.79) <strong>and</strong> (5.86) shows that the two equations are identical if thefollowing two conditions hold.1. The statistical weight of a particle emerging from a collision is independent of thenumber of secondaries in the collision <strong>and</strong> is the same in both the physical <strong>and</strong>hypothetical games. Denoting this weight by W", the condition is formulated asw: = w; = w (5.87)2. The contribution functions in the two games are related asJdP"C(P',P")W"f,(P',P") = c a(P')W a f a(P') + c s(P')|dP"C s(P',P'')W''f s(P',P'')+ 6 f(P') 2 nq„(P') fdP"C„(P',P")WT n(P',P") (5.88)H= I •>Condition 2 requires that the expected score contributed by a particle entering a collision atP be equal in the two games. We have thus proven the following.Theorem 5.3 — Under conditions 1 <strong>and</strong> 2 the hypothetical nonmultiplying game resultsin the same expected score as the (possibly multiplying) physical game.•For the sake of illustration let us consider a couple of cases when conditions 1 <strong>and</strong> 2are satisfied. Obviously, if the games are analog (with respect to the physical <strong>and</strong> hypotheticaltransport processes, respectively), then all the statistical weights are unity <strong>and</strong> Equation(5.87) is automatically met. The rules of generating the statistical weights in an unbiasednonanalog game will be derived in Section 5.V.B; nevertheless, it is heuristically obviouswithout any derivation that if in the nonanalog game only the transition kernel t is differentfrom the analog kernel, but the nonanalog <strong>and</strong> analog collision kernels are identical, thenthe statistical weights of the particles are not changed in a collision [except for the multiplicationby c(P') in the hypothetical game] <strong>and</strong> Equation (5.87) is again satisfied. Finally,if the nonanalog <strong>and</strong> analog kernels differ only in the mean number of secondaries percollision, i.e., if Equation (5.83) can be rewritten asC(P',F') = c(P')C,(P',P") =c(P')C(P',P")/c(P)then the postcollision weight is not affected by the type of the collision <strong>and</strong> condition 1holds. Note that in the three examples above, Equation (5.87) is satisfied in such a way thatthe pre- <strong>and</strong> postcollision weights are related asW" = W" = a(P')W' (5,89)with some function a(P') independent of the postcollision coordinates P". This relation holdsin most practical cases <strong>and</strong> if so, the postcollision coordinates P" are selected from the analogdensity, which also means that there is no need to change the statistical weight of the particleswhen they emerge from a collision. The rigorous proof of the reasoning above follows fromTheorem 5.8 in Section 5.V.B.The following examples demonstrate that condition 2 is also met in most practical cases.Let us first consider the case when the contribution functions f rdo not depend on the

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!