A multi-dimensional freeze drying analysis program, which took simultaneous heat and mass transfer, sublimation of ice and motion of sublimation front into full account, was developed using finite volume method with fixed grid. The effect due to the motion of sublimation interface was modeled by an enthalpy formulation streamlined for the freeze drying problems. The efficiency and accuracy of the program was validated by solving one- and two-dimensional freeze drying problems frequently encountered in industrial processes. Finally, the freeze drying processes of cylinder and slab objects with permeable side surfaces were simulated, which demonstrated the capability of the present analysis program in solving multi-dimensional freeze drying problems with complex sublimation interface configurations.