This study introduces an efficient and accurate two-stage explicit computational scheme for solving partial differential equations (PDEs) containing first-order time derivatives. The suggested method is a modification of the classical Runge–Kutta scheme that introduces a new first-stage formulation. This minimizes numerical error with moderate step sizes while preserving the stability region of the classical method. Spatial discretization is performed using a sixth-order compact finite-difference scheme to obtain high-resolution solutions. The analysis of stability and convergence is strictly determined for both scalar and system forms of convection–diffusion-type equations. To illustrate the suitability of the method, a dimensionless mathematical model of the unsteady, incompressible, laminar flow of a Prandtl-type non-Newtonian nanofluid over a Riga plate is considered, accounting for viscous dissipation, thermophoresis, Brownian motion, and a magnetic field. Here, the Prandtl ternary nanofluid is defined as a non-Newtonian nanofluid that follows the Prandtl rheological model, and it exhibits three critical transport phenomena: heat conduction, viscous dissipation, and nanoparticle diffusion. Representative values of the Prandtl number Pr = 3 and Reynolds number Re = 5 are used to perform the simulation, and other parameters, including but not limited to the Hartmann number Ha, Williamson number We, thermophoresis Nt and Brownian motion Nb, are varied to evaluate the flow behavior. Moreover, an artificial neural network (ANN)-developed surrogate model is used to calculate the skin friction coefficient and the local Sherwood number, using five input parameters: the Reynolds number, Prandtl number, Schmidt number, Brownian motion parameter, and thermophoresis parameter. The governing partial differential equations yield high-fidelity numerical data used to train the surrogate model. The data is split into 80% for training, 10% for validation, and 10% for testing. The ANN is tested using regression analysis and error histograms, which demonstrate high accuracy and generalization capacity. Numerical simulation combined with AI-based prediction is a cost-efficient method for real-time estimation of complex non-Newtonian nanofluid systems.
Nawaz et al. (Fri,) studied this question.