The heat-conducting nature of blood is critical in the human circulatory system and features also in important thermal regulation and blood processing systems in biomedicine. Motivated by these applications, in the present investigation, the classical Graetz problem in heat transfer is extended to the case of a bio-rheological fluid model. The Quemada bio-rheological fluid model is selected since it has been shown to be accurate in mimicking physiological flows (blood) at different shear rates and hematocrits. The steady two-dimensional energy equation without viscous dissipation in stationary regime is tackled via a separation of variables approach for the isothermal wall temperature case. Following the introduction of transformation variables, the ensuing dimensionless boundary value problem is solved numerically via MATLAB based algorithm known as bvp5c (a finite difference code that implements the four-stage Lobatto IIIa collocation formula). Numerical validation is also presented against two analytical approaches namely, series solutions and Kummer function techniques. Axial conduction in terms of Péclet number is also considered. Typical values of Reynolds number and Prandtl number are used to categorize the vascular regions. The graphical representation of mean temperature, temperature gradient, and Nusselt numbers along with detail discussions are presented for the effects of Quemada non-Newtonian parameters and Péclet number. The current analysis may also have potential applications for the development of microfluidic and biofluidic devices particularly which are used in the diagnosis of diseases in addition to blood oxygenation technologies.
Keywords: Concentrated suspension; Graetz problem; Newton-Raphson method; Péclet number; Quemada fluid; S-L boundary value problem; Simpson’s rule; bvp5c.