This study introduces a modified computational scheme for handling linear and nonlinear fractal time-dependent partial differential equations. The method is constructed using three different stages that provide third-order accuracy in the fractal time variable. The stability of the approach is examined using scalar fractal models and Fourier analysis, while convergence is established for coupled convection–diffusion systems. The numerical algorithm is applied to analyze the mixed convective flow of a Carreau–Yasuda non-Newtonian fluid over stationary and oscillating plates under the influence of viscous dissipation and magnetic field effects. For spatial discretization, the incompressible continuity equation is handled by a first-order difference scheme, whereas higher-order compact schemes are implemented for the momentum, thermal, and concentration equations. The numerical findings show that increasing the Weissenberg number and magnetic field inclination reduces the velocity distribution. An accuracy assessment against existing numerical techniques demonstrates that the present method yields smaller computational errors, particularly when central difference schemes are used in space. In addition, a surrogate machine learning model is developed to predict the skin friction coefficient and local Nusselt number using Reynolds, Weissenberg, Prandtl, and Eckert numbers as input features. The predictive capability of the model is validated through Parity plots, bar charts for sensitivity analysis, scatter visualization, and Taylor Diagrams, confirming strong agreement with the numerical results.
Nawaz et al. (Thu,) studied this question.