Finite Difference Modeling of Time Fractal Impact on Unsteady Magneto-hydrodynamic Darcy–Forchheimer Flow in Non-Newtonian Nanofluids with the q-Derivative

Abstract

This contribution addresses a fractal numerical scheme that can be employed for handling fractal time-dependent parabolic equations. The numerical scheme presented in this contribution can be used to discretize integer order and fractal derivatives in a given differential equation. Therefore, the scheme and results can be used for both cases. The proposed finite difference scheme is based on two stages. Fractal time derivatives are discretized by employing the proposed approach. For the scalar convection–diffusion equation, we derive the stability condition of the proposed fractal scheme. Using a nonlinear chemical reaction, the approach is also used to solve the Quantum Calculus model of a Williamson nanofluid’s unsteady Darcy–Forchheimer flow over flat and oscillatory sheets. The findings indicate a negative correlation between the velocity profile and the porosity parameter and inertia coefficient, with an increase in these factors resulting in a drop in the velocity profile. Additionally, the fractal scheme under consideration is being compared to the fractal Crank–Nicolson method, revealing that the proposed scheme exhibits a superior convergence speed compared to the fractal Crank–Nicolson method. Several problems involving the motion of non-Newtonian nanofluids through magnetic fields and porous media can be investigated with the help of the proposed numerical scheme. This research has implications for developing more efficient heat transfer and energy conversion devices based on nanofluids.