Dynamic simulation of a solar-powered solid absorption ice maker(II). The finite difference method

In order to verify the correctness of the model of reaction fronts presented in part I of this paper, the finite difference method was used to solve the unsteady heat conduction equations inside the reactor. After discreting the unsteady heat conduction equations in time coordinate, a system of ordinary differential equations was obtained, then it was solved using a general software package DALI which can simulate the system described by one order ordinary differential equations and algebraic equations. Due to the axisymmetry of reactor area and boundary conditions, only 1/4 of reactor area needs to be considered. Computation results show a good agreement between these two methods. The model of reaction fronts saves a lot of computation time, so it can be used to guide the design of the solar-powered solid absorption ice maker.