4.2. Reactor Power

At any time t, the local power production at position r is assumed to be given by the space-time separated function:

(4.2-1)Q(r,t)=ψt(t)S(r)

where ψt(t) is the dimensionless, normalized power amplitude and S(r) is the steady-state reactor power in watts being produced in an axial node at location r. In terms of input quantities, S(r) is given by the product of POW and PSHAPE. Initially, the power amplitude has a value of unity and S(r) is normalized to the total steady-state reactor power. Appendix 4.1 contains a description of the internal normalization of PSHAPE performed by SAS4A/SASSYS‑1. The time-dependent power amplitude is assumed to be made up of the sum of two components:

(4.2-2)ψt(t)=ψf(t)+ψh(t)

where ψh(t) comes from the decay of fission and capture products. These two components have been separated to allow the simulation of both short- and long-term transients.

The direct fission component of the power amplitude is given by

(4.2-3)ψf(t)=ψf(0)ϕ(t)

where ϕ(t) is the dimensionless, normalized fission power amplitude given by the point reactor kinetics model:

(4.2-4)˙ϕ(t)=ϕ(t)δk(t)βΛ+iλiCi(t)

with the initial condition ϕ(0)=1.

In Eq. (4.2-4), δk(t) is the net reactivity, β is the total effective delayed-neutron fraction, Λ is the effective prompt neutron generation time, and λi is the decay constant for the delayed-neutron precursor isotope whose normalized population is Ci(t). The physical interpretation of the terms in the point reactor kinetics equation is made by Henry [4-2] and also by Bell and Glasstone [4-3].