The Jeffrey fluid model is capable of accurately characterizing the stress relaxation behavior of non-Newtonian fluids, which a normal viscous fluid model is unable to perform. The primary objective of this paper is to provide a comprehensive investigation into the effects of MHD and thermal radiation on the 3D Jeffery fluid flow over a permeable irregular stretching surface. The consequences of the Darcy effect, variable thickness and chemical reaction are also considered. The phenomena have been modeled as a nonlinear system of PDEs. Using similarity substitution, the modeled equations are reduced to a dimensionless system of ODEs. The parametric continuation method (PCM) is used to determine the numerical solution to the obtained sets of nonlinear differential equations. The impact of physical parameters on temperature, velocity and mass profiles are presented through Figures and Tables. It has been noticed that the energy profile magnifies with the increment of porosity term, thermal radiation and heat source term, while diminishing with the flourishing upshot of power index and Deborah number. Furthermore, the porosity term and wall thickness parameter enhance the skin friction.