4. METHOD OF SOLUTION
The finite difference method is employed, with up to 100 mesh points in the active part of the channel and 10 more points in the riser. The program is intended for forced circulation and highly pressurized systems, for which the pressure drop in the channel may be considered negligible compared to the general pressure level, so that the coolant fluid properties may be assumed to be independent of space and time. The driving pressure may be taken as a quadratic function of the mass flow, to simulate pump characteristics, or given as a time table. The power distribution in the heating element is given as a fixed arbitrary shape, while the power level is any tabulated function of time.