Sublimation heat transfer occurs in a wide range of engineering processes, such as accelerated freeze drying (AFD), energy storage, and food technology. Particularly in the microwave AFD process, preservation of material with the least possible energy consumption is desirable. In connection with this, it is of interest to analyze the effect of temperature/concentration dependent heat/mass transfer properties. Given the limited literature available on sublimation, there is a general lack of physical understanding of this particular problem. The present work analyzes the nonlinear sublimation process driven by convective heat/mass transfer and evaporation of water vapor using the Legendre wavelet collocation method (LWCM). Results from the present work are shown to be in excellent agreement with the exact solution of the special case of a linear problem. Further, the present numerical technique shows good agreement with finite difference method in case of a completely nonlinear model. The model is used for a comprehensive investigation of the impact of the problem parameters, on the rate of sublimation. It is found that the sublimation rate increases with increasing values of β1 and decreasing values of β2. The impact of other dimensionless problem parameters such as Péclet numbers Pe1 and Pem, convection due to mass transfer of water vapor β, latent heat of sublimation l0 and Luikov number Lu on sublimation process is also discussed in detail. These observations offer a comprehensive theoretical and mathematical understanding of sublimation heat/mass transfer for improving the performance and efficiency of freeze-drying and related engineering processes.