In this article, conservation laws, temperature dependent viscosity, thermal conductivity and mass diffusion models are used for mathematical modeling of flow of Casson fluid induced by nonlinear axisymmetric stretching sheet exposed to non-uniform magnetic field. Mathematical models are solved by finite element method. Convergence is ensured and grid independent analysis is performed. Several numerical simulations are executed to investigate the behavior of velocity, temperature and concentration under the variation of parameters. The impact of viscous dissipation, Joule heating, variation of viscosity and thermal conductivity is investigated through numerical simulations. The developed computer program is tested by comparing with computed results. The intensity of applied magnetic field and Joule heating are directly proportional to each other and an increase in the intensity of magnetic field causes an increase in the temperature. A rise in viscosity due to a rise in temperature has remarkable impact on both velocity and temperature. However, an increase in thermal conductivity (due to rise in temperature) results a significant increase in the temperature as conduction of heat from hot surface into the fluid speeds up. The wall heat flux increases when intensity of the magnetic field is increased.