Abstract: Based on the Oceanic Regional Circulation and Tidal Model (ORCTM) and its ice-sea coupling module, this paper establishes a numerical model of sea ice in the Bohai Sea by comprehensively considering the dynamic and thermodynamic processes of sea ice generation and elimination. Firstly, the simulation results of barotropic tide and sea surface temperature are evaluated to verify the reliability of the model. Secondly, sensitivity experiments are carried out with three key sea ice thermodynamic parameters, revealing the effect of critical ice thickness on the density of nascent ice. The sea ice area increases as critical ice thickness decreases, and the optimal value of critical ice thickness for the Bohai Sea is 0.1 m. Significant effects of sensible and latent heat flux coefficients on the heat fluxes at the sea-air and air-ice interfaces are found. The relationship between salinity and freezing point in the Bohai Sea is modeled in conjunction with previous studies. Finally, a typical sea ice process in the Bohai Sea in the winter of 2009/2010 is simulated using adjusted model parameter configuration and the Japan Meteorological Agency JRA55-TL319 reanalysis data as meteorological forcing. The results show that the simulated sea ice coverage is consistent with satellite remote sensing data in the same period, indicating that the established sea ice numerical model can well reproduce the evolution of the Bohai sea ice generation and elimination, and can be used in the medium- and long-term numerical simulation of the Bohai Sea ice.