The non-orthogonal boundary-fitted coordinate transformation method is applied to the solution of steady three-dimensional conservation equations of mass, momentum, energy and speciescontinuity to obtain the laminar velocity, temperature and concentration fields for simulation of polymerization of styrene in arbitrary cross-sectional duct reactors. Variable physical properties (except for speci...