Page 48 - Dynamics and Control of Nuclear Reactors
P. 48

4.4 Analytical solutions   39




                    1.35

                     1.3
                    1.25

                     1.2
                   P/P(0)
                    1.15

                     1.1

                    1.05
                      1
                       0         20        40        60        80       100       120
                                                   Time (s)
                  FIG. 4.4
                  Fractional power response to a ramp reactivity input, followed by a step change to zero
                  reactivity.





                  4.4 Analytical solutions
                  Analytical solutions for the six-delay neutron precursor group zero power model and
                  for more detailed power reactor models are too tedious for practical implementation.
                  However, an analytical solution to the point kinetics equation with one delay neutron
                  group is easily performed for a step change in reactivity. This solution is useful in
                  understanding the nature of reactor transients. Laplace transforms are a convenient
                  tool for solving linear, constant coefficient differential equations. Appendix D pro-
                  vides a brief summary of Laplace transform theory. Appendix H shows the derivation
                  of the analytical solution of the point reactor kinetics equations for a step change in
                  reactivity. The parameters used in the simulation are as follows
                     ρ ¼ +10 cents ¼ 0.1*0.0067 ¼ 0.00067
                     β ¼ 0.0067
                     Λ ¼ 0.00001 s.
                             -1
                     λ ¼ 0.08 s
                  The Laplace transform of the fractional power is given by
                                                      β
                                                 Λ  λ +  + s
                                       Ps ðÞ          Λ
                                           ¼                                     (4.9)
                                              2
                                       P 0ðÞ  Λs + β  ρ + λΛð  0  Þ s λρ 0
                  The roots of the denominator polynomial are used to solve for P(t)/P(0). The two
                  roots of this polynomial are given by.
   43   44   45   46   47   48   49   50   51   52   53