The aim of this paper is to study discretizations of convection-diffusion-reaction equations using splitting methods. Estimates for the physical splitting errors and the numerical splitting errors are established. These estimates lead to the selection of optimal sequences and coupling of physical phenomena and adequate use of implicitness and explicitness. Numerical simulations of two chemical industry problems are included.