The objective of the work presented is to formulate a mathematical description of frazil ice dynamics. The formulation is to be in balance with the current knowledge of the physical processes, for example secondary nucleation. As the knowledge of some of these processes is fragmentary, this means that a conceptually simple formulation is sought. A number of processes are known to influence the supercooling rate and the frazil ice production. The present formulation accounts for the following processes: initial seeding, secondary nucleation, gravitational removal, growth due to cooling of water volume and flocculation/break up. Equations are formulated for these present considering a resolution in time and radius of particles but not in space (well-mixed jar). The equations are solved using a simple explicit numerical scheme. Preliminary results indicate that the model can be calibrated to describe the experimental results reported in the literature. It is mainly the supercooling curves that are used for comparison but some information about the crystal size distribution is also considered. It is to be noted that the model is calibrated to fit the experiments, due to the lack of detailed mathematical description of some of the physical processes. Sensitivity analysis is also used in order to establish that the model behaves according to experimental findings and expectations. The main conclusion of the study is that a fairly simple mathematical model can be formulated and calibrated, which fits the experimental data reported in the literature hitherto. It is further concluded that a resolution in radial space gives additional insight into the dynamics of the process. The evolution of the size distribution and its sensitivity to seeding and dissipation rate has been predicted with results that look physically plausible.