A model is developed to simulate the solidification of the steel shell in the mold region of the continuous casting process. Conduction-dominated temperature fields in the mold, mold flux, steel shell, and molten steel regions are determined through the development of an evolution equation for the solidifying front. This equation is derived in the limit of small aspect ratio, mold width to height, using asymptotic methods. These results are coupled with a lubrication-theory model for the mold flux region. This model assumes a temperature-dependent viscosity for the mold flux and allows for solidification of the flux at temperatures below a critical value. System response to changing casting speeds, superheat, mold wall temperatures, and mold flux properties is investigated.