Dynamic simulation of non-Newtonian boundary layer flow: An enhanced exponential time integrator approach with spatially and temporally variable heat sources

Scientific inquiry into effective numerical methods for modelling complex physical processes has led to the investigation of fluid dynamics, mainly when non-Newtonian properties and complex heat sources are involved. This paper presents an enhanced exponential time integrator approach to dynamically simulate non-Newtonian boundary layer flow with spatially and temporally varying heat sources. We propose an explicit scheme with second-order accuracy in time, demonstrated to be stable through Fourier series analysis, for solving time-dependent partial differential equations (PDEs). Utilizing this scheme, we construct and solve dimensionless PDEs representing the flow of Williamson fluid under the influence of space- and temperature-dependent heat sources. The scheme discretizes the continuity equation of incompressible fluid and Navier–Stokes, energy, and concentration equations using the central difference in space. Our analysis illuminates how factors affect velocity, temperature, and concentration profiles. Specifically, we observe a rise in temperature profile with enhanced coefficients of space and temperature terms in the heat source. Non-Newtonian behaviours and geographical/temporal variations in heat sources are critical factors influencing overall dynamics. The novelty of our work lies in developing an explicit exponential integrator approach, offering stability and second-order accuracy, for solving time-dependent PDEs in non-Newtonian boundary layer flow with variable heat sources. Our results provide valuable quantitative insights for understanding and controlling complex fluid dynamics phenomena. By addressing these challenges, our study advances numerical techniques for modelling real-world systems with implications for various engineering and scientific applications.