Abstract The purpose of the research is to examine a tangent hyperbolic nanofluid flowing in three dimensions (3D) axisymmetrically on an unsteady rotatory stretching disk over a DarcyForchheimer porous medium. First order initial value problems (IVPs) are generated from the governing partial differential equations (PDEs) through the use of similarity transformation and linearization. The Runge-Kutta sixth order (RK6) is utilized to solve the IVP system using the shooting technique and the built-in Python software ”fsolve model10.” Articles that have already been published are used to validate the implemented approach. Graphs are used to examine how various parameters affect velocity, temperature, and concentration. Additionally, the behavior of heat, mass flux, and skin friction in response to different parameters is investigated. The study’s findings showed that as the Forchheimer number and velocity slip parameter increased, the nanofluid’s radial and tangential velocities decreased as well. As temperature and concentration slip parameters increase, correspondingly, thicker and thinner boundary layer structures are seen. Both the rates of heat and mass transfers are initiated for an increase Eckert and Prandtl numbers and demotivated for power-law index number.