The researchers modeled the flow of a non-Newtonian Walter’s-B nanoliquid containing nanoparticles, where Brownian motion and thermophoresis drive the transport of thermal energy and species. The Cattaneo-Christov formulation replaces the classical Fourier and Fick laws, introducing finite relaxation times that account for the memory effects in heat and mass diffusion. Gyrotactic microorganisms, which orient themselves relative to the gravitational field and swim in a directed manner, add a further layer of biological transport to the system.

Thermal radiation was included in the energy equation, contributing an additional radiative flux that augments the overall heat transfer away from the stretching boundary. The sheet was assumed to stretch linearly, and the boundary layer approximation was applied to reduce the governing partial differential equations to a set of nonlinear ordinary differential equations. The homotopy analysis method was then employed to obtain semi-analytical solutions, with the auxiliary convergence parameters tuned to ensure rapid and stable convergence of the series.

The results demonstrate that increasing the radiation parameter substantially enhances the local Nusselt number, reflecting the additional radiative flux that contributes to the thermal energy transport. Brownian motion also improves heat transfer performance by augmenting the thermal diffusion of nanoparticles. Conversely, higher values of the thermophoresis parameter reduce the Nusselt number, as the directed migration of particles away from the boundary diminishes the effective thermal gradient.

For mass transport, the Sherwood number increases gradually with the Schmidt number, indicating that the molecular diffusivity of species plays a key role in the concentration boundary layer. The Cattaneo-Christov relaxation times introduce memory effects that modify the classical concentration profiles, and the study quantifies how these finite relaxation times alter the species transport relative to the Fickian limit.

The transport of gyrotactic microorganisms was analyzed through the motile microorganism density number, which increases with both the bioconvection Lewis number and the bioconvection Péclet number. The Lewis number represents the ratio of thermal to microorganism diffusivity, and higher values augment the density gradient near the wall. The Péclet number captures the directed swimming speed of the microbes relative to their random diffusion, and increasing values enhance the accumulation of motile species near the boundary.

The study also examined the effects of the Walter’s-B fluid parameters on the skin friction coefficient, demonstrating that the non-Newtonian characteristics of the nanoliquid modify the shear stress at the stretching surface. The homotopy analysis solutions were validated against previously published results for limiting cases, with excellent agreement observed across the parameter ranges examined.

This comprehensive analysis provides a detailed framework for understanding the coupled transport of thermal energy, species, and microorganisms in non-Newtonian nanoliquids with finite relaxation times. The findings are relevant to applications in biomedical engineering, biotechnology, and materials processing, where the directed transport of biological species in complex fluids plays a critical role. The semi-analytical solutions offer a reliable foundation for future studies exploring more complex geometries and boundary conditions.

The porous medium through which the Walter’s-B nanoliquid flows is described by the Darcy-Forchheimer framework, which accounts for two distinct sources of hydraulic resistance. The classical Darcy term captures the linear viscous drag exerted by the solid matrix on the percolating fluid, while the Forchheimer term adds a quadratic contribution that becomes significant at higher pore velocities, when inertial effects within the pores can no longer be neglected. In the present formulation, increasing either the porosity parameter or the Forchheimer inertia coefficient raises the effective resistance offered by the medium, and the computed velocity profiles respond by flattening and thinning near the stretching surface. This behavior is physically intuitive: a more restrictive pore structure dissipates more of the momentum imparted by the moving boundary, confining the motion to a narrower region adjacent to the sheet.

The viscoelastic character of the base fluid acts in the opposite direction. Walter’s-B fluids possess a short memory of past deformation, and the associated viscoelastic parameter enhances the effective momentum diffusivity of the liquid. As this parameter grows, the fluid transmits the stretching action of the sheet more efficiently into the bulk, accelerating the flow and broadening the momentum boundary layer. The interplay between this elastic momentum transfer and the porous-medium drag therefore provides a tunable balance: engineers working with polymer solutions or lubricants in permeable or fibrous media can, in principle, select operating conditions that either promote or suppress boundary-layer growth depending on the desired transport outcome.

Stratification emerges as one of the most consequential features of the model. In the thermal case, the ambient temperature is assumed to vary linearly with distance from the surface rather than remaining uniform. This layered temperature environment reduces the wall-to-fluid temperature difference that drives conductive and radiative heat exchange, and the computed temperature profiles respond by compressing toward the sheet. The thermal penetration depth shrinks as the thermal stratification parameter increases, which is the mirror image of the widening produced by stronger radiation, Brownian motion, or thermophoresis. An analogous mechanism operates for the solutal field: a stratified ambient concentration diminishes the driving concentration gradient, so nanoparticle concentration distributions contract as solutal stratification intensifies. Chemical reaction effects reinforce this contraction, since reactive consumption of the diffusing species depletes the concentration field more rapidly than diffusion alone would allow.

The bioconvective component of the model introduces a third stratification mechanism acting on the density of motile microorganisms. Gyrotactic algae and similar self-propelled cells accumulate in response to gradients in their own population density, and when the ambient microbial density is itself stratified, the driving gradient for bioconvection is weakened. The study shows that increasing the microbial stratification factor lowers the local population density of motile organisms near the surface. The bioconvection Péclet number, which measures the ratio of directed swimming speed to the random diffusivity of the cells, acts in the same direction: larger Péclet numbers reduce the local density number, indicating that stronger directed swimming disperses the accumulation that would otherwise form near the boundary. These results underscore that bioconvection intensity is not fixed by biology alone but can be steered through the ambient density environment imposed on the system.

The wall transport coefficients tie the field solutions back to quantities of direct engineering interest. The skin-friction coefficient, which quantifies the shear stress exerted on the stretching surface, increases with the Darcy resistance, the Forchheimer inertia, and the viscoelastic parameter. This reflects improved interfacial momentum transfer: a more resistive medium or a more elastically responsive fluid both transmit greater traction to the wall. For heat transfer, the Nusselt number grows with the intensity of radiative effects, consistent with the additional radiative flux leaving the surface, but falls as thermophoresis strengthens, because nanoparticles drifting away from the hot wall carry less thermal energy across the near-wall gradient. The motile microorganism density number rises with the bioconvection Lewis number, confirming that the diffusivity of the microorganisms relative to other transport mechanisms is central to sustaining the bioconvection flux at the boundary.

Methodologically, the reliance on the homotopy analysis method deserves emphasis. Unlike purely numerical schemes, this semi-analytical technique expresses the solution as a convergent series whose structure is controlled by an auxiliary parameter, denoted by the convergence-control symbol h. By plotting h-curves for each governing equation, the authors identified the ranges of this parameter within which the series converges reliably, and residual error tables provided quantitative confirmation that successive approximations were approaching a stable solution. This combination of graphical and tabular convergence verification is particularly valuable for systems of strongly coupled nonlinear equations, where naive series expansions may diverge or converge to spurious branches.

The simultaneous treatment of thermal, solutal, and microbial stratification within a single framework is relatively rare in the literature on stretching-sheet flows, and the present results clarify how the three layering mechanisms interact with the distinct transport physics of a viscoelastic nanofluid. Because stratification suppresses the driving gradients for heat, species, and microorganisms alike, while radiation, Brownian diffusion, and thermophoresis amplify them, the model effectively maps out a competition between gradient-eroding and gradient-enhancing influences. This mapping offers practical guidance for systems such as solar ponds, thermal energy storage tanks, chemical reactors, and bioreactors, where layered density environments are common and where controlling the near-wall transport of heat, dissolved species, and suspended biological cells simultaneously is essential to performance.

Subject of Research: Semi-analytical transport of heat, mass, and gyrotactic microorganisms in a Walter's-B nanoliquid with Cattaneo-Christov double diffusion over a stretching sheet

Article Title: Thermal and mass management of Walter-B nanoliquid with heat radiation and the Cattaneo–Christov double diffusion model incorporating microorganisms: convergence solution for stratified effect

Article References: Ali, F., Hussain, I., Zafar, S. S., & Ahmed, M. F. (2026). Thermal and mass management of Walter-B nanoliquid with heat radiation and the Cattaneo–Christov double diffusion model incorporating microorganisms: convergence solution for stratified effect. Frontiers in Thermal Engineering, 6, Article 1864538. https://doi.org/10.3389/fther.2026.1864538

Image Credits: AI Generated

DOI: 10.3389/fther.2026.1864538

Keywords: Walter's-B nanoliquid, Cattaneo-Christov double diffusion, gyrotactic microorganisms, stretching sheet, homotopy analysis, thermal radiation, Brownian motion, thermophoresis, bioconvection, Nusselt number, Thermal, mass

The team modeled the working fluid as a Casson fluid, a mathematical description famously suited to blood. Unlike water or air, which begin to flow the instant any force is applied, a Casson fluid behaves like a solid until the applied stress crosses a critical threshold known as the yield stress. Only after that barrier is breached does it begin to move, and even then the relationship between stress and deformation is nonlinear. This yield-stress behavior captures a real feature of blood, which resists deformation in small vessels in ways that ordinary Newtonian models cannot reproduce. The Casson parameter in the model expresses the ratio of yield stress to plastic viscosity, and it turns out to be one of the most powerful levers controlling how fast fluid moves near the vessel wall.

The geometry at the heart of the study is deliberately idealized. Rather than simulating a full network of branching arterioles, venules, and capillaries, the researchers represented a segment of micro-vessel as a curved sheet that stretches exponentially along its length. This choice isolates two effects the authors consider critical for near-wall transport: the curvature of the cross-section, which alters pressure and shear stress distributions, and the distending motion of the wall, which mimics the pulsing and deformation that real vessels experience with every heartbeat. The authors are careful to stress that the model is not a physiologically resolved simulation of the microvascular system. It is a boundary-layer formulation designed to reveal qualitative trends, not to compute precise drug-delivery efficiencies.

Into this curved, stretching environment the researchers packed an unusually complete set of physics. The fluid is electrically conducting and subjected to a transverse magnetic field, which generates a Lorentz force that opposes the motion. The channel walls bound a porous medium described by both Darcy resistance, which dominates at low speeds, and Forchheimer inertia, which adds extra momentum loss when the fluid pushes faster through the porous structure. Heat transport includes thermal radiation treated through the Rosseland approximation, viscous dissipation captured by the Brinkman number, and the cross-diffusion effects of Brownian motion and thermophoresis that characterize nanoscale particle transport. The concentration field incorporates Arrhenius activation energy, the concept proposed by Svante Arrhenius in 1889, which sets the minimum energy barrier molecules must overcome before a chemical reaction can proceed. Finally, the model tracks motile microorganisms — living, swimming carriers — whose collective drifting, known as bioconvection, is included as a stand-in for biological drug-delivery vehicles navigating toward hypoxic tumor regions.

Solving the resulting equations is no small feat. The governing partial differential equations for momentum, energy, concentration, and microorganism density are coupled and nonlinear, and the curvature terms make them more complex than their flat-plate counterparts. The authors applied similarity transformations to collapse the partial differential equations into ordinary differential equations, then handed the system to MATLAB’s BVP4C boundary-value solver. The higher-order equations were rewritten as a system of ten first-order equations, and the solver iterated until the boundary conditions were satisfied to a precision of one part in a million. As a check on accuracy, the team compared their computed skin friction coefficients against previously published results for a special case with no magnetic field and effectively infinite Casson parameter, finding agreement to within roughly one part in ten thousand across a range of curvature values.

The velocity results carry the study’s most direct message for drug delivery. Increasing the magnetic parameter strengthens the Lorentz force, which decelerates the fluid and keeps drug particles lingering close to the vessel wall rather than sweeping away downstream. Raising the Darcy permeability parameter — which in this formulation represents greater resistance from the porous medium — likewise diminishes the velocity profile and thins the hydrodynamic boundary layer. The Forchheimer parameter, capturing inertial drag at higher flow speeds, produces the same retarding trend. Even the Casson parameter itself reduces velocity as it increases, meaning that changes in the yield-stress character of blood can alter how long particles remain near a capillary wall. In the drug-delivery framing, all of this resistance is potentially useful: a slower, wall-hugging flow gives pharmaceutical particles more opportunity to adhere to or penetrate the intended target tissue, particularly in tumor regions where distribution occurs through permeable biological material.

The thermal findings are equally consequential. The Brinkman number, which measures heat generated by viscous friction relative to heat conducted away, raises the temperature field as it increases — more fluid friction means more built-in warmth. The thermal Biot number, the ratio of surface heat-transfer resistance to internal fluid resistance, also elevates temperature when it grows, indicating stronger thermal exchange at the wall and, by extension, enhanced interaction between drug particles and the vessel surface. Thermal radiation amplifies the effect further, intensifying the temperature throughout the boundary layer, with the strongest influence closest to the curved stretching surface. Since drug absorption and penetration at a target site depend on local thermal conditions, these parameters describe how the flow itself could shape the environment a drug particle encounters.

On the mass-transport side, the activation energy parameter widens the concentration profile, meaning that a higher energy barrier for chemical reactions leaves more of the dissolved drug intact in the fluid rather than consumed by reaction. Brownian motion and thermophoresis work in the opposite direction, narrowing the concentration profile while boosting the Sherwood number, the dimensionless measure of mass transfer from the surface into the fluid. Thermophoresis, in particular, drives particles from warmer toward cooler regions, a mechanism the authors note could actively steer medication through temperature gradients. Meanwhile, the bioconvective parameters govern the microorganisms: a higher bio-Schmidt number suppresses the spreading of the swimming carriers, confining them near the wall where they are easier to localize, while the bio-convective Peclet number balances convective sweeping against random diffusive dispersal, with higher values reducing microorganism density in the outer field.

The authors are candid about the limits of what they have built. Real microvascular networks branch and taper asymmetrically, generating secondary flows that a two-dimensional curved sheet cannot capture. Physiological blood flow is pulsatile and unsteady, driven by the cardiac cycle, whereas the model assumes steady conditions. Blood itself is treated as a single-phase homogeneous Casson fluid, ignoring the discrete nature of red blood cells and the cell-free layer that forms near vessel walls. The team proposes in vitro validation using biomimetic curved microfluidic channels and blood-analogue fluids as the necessary next step, along with extending the framework to fully three-dimensional, pulsatile formulations. Until then, the model should be read as a mechanistic map rather than a clinical prescription.

Even with those caveats, the study offers a genuinely useful synthesis. By combining, for the first time in a single boundary-layer model, an exponentially curved stretching geometry with Darcy–Forchheimer porous resistance, Arrhenius activation energy, magnetohydrodynamics, and oxytactic bioconvection, the work identifies permeability and the Casson parameter as the most influential levers on near-wall particle residence. That insight could eventually inform the design of stents, catheters, and other curved microvascular devices engineered to influence how long therapeutic particles linger at a treatment site. The message is elegant in its simplicity: in the narrow, curving world of the microvasculature, the physics of resistance — magnetic, porous, and yield-stress alike — may be the very tool that keeps a drug where it belongs.

Subject of Research: Numerical boundary-layer modeling of magnetohydrodynamic Casson fluid flow with Arrhenius activation energy and bioconvection over an exponentially curved stretching sheet for targeted drug delivery

Article Title: Numerical analysis of Casson fluid flow past an exponentially curved stretching sheet with arrhenius activation energy: implications for targeted drug delivery

Article References: Pisho, M. J., Shankar, G., Siva, E., & Loganathan, K. (2026). Numerical analysis of Casson fluid flow past an exponentially curved stretching sheet with arrhenius activation energy: implications for targeted drug delivery. Results in Engineering, 32, Article 113049. https://doi.org/10.1016/j.rineng.2026.113049

Image Credits: AI Generated

DOI: 10.1016/j.rineng.2026.113049

Keywords: Casson fluid, targeted drug delivery, magnetohydrodynamics, Arrhenius activation energy, bioconvection, porous medium, Darcy-Forchheimer, stretching sheet, boundary layer, thermal radiation, microvasculature, BVP4C

The study’s central innovation lies in its sheer scope. Rather than examining one or two phenomena in isolation, the researchers constructed a single mathematical framework that couples magnetohydrodynamics, porous media, thermal radiation, viscous dissipation, Soret and Dufour cross-diffusion, the Buongiorno nanofluid model, multiple velocity and thermal slip conditions, and bioconvection driven by gyrotactic microorganisms. The working fluid itself is a dual material: a Casson fluid, which behaves like a plastic solid until a critical yield stress is exceeded, blended with micropolar characteristics that allow suspended particles to spin independently of the bulk motion. This micro-rotation matters enormously in real materials, because the internal structure of blood, polymer melts, and many slurries directly influences how efficiently they carry heat.

Mathematically, the team began with the full Casson rheological equation, in which the Cauchy stress tensor depends on whether the product of the deformation rate components exceeds a critical threshold. They then wrote the conservation equations for momentum, angular momentum, energy, nanoparticle concentration, and microorganism density for flow over a vertical linearly stretching sheet, where the wall velocity grows linearly with distance from the origin. Buoyancy enters the momentum balance through three separate Grashof numbers, one each for thermal, solutal, and microorganism density differences, meaning the fluid can be driven simultaneously by heat, by chemical concentration gradients, and by the swimming of bacteria. A transverse magnetic field of strength B0 imposes a Lorentz retarding force, while the porous matrix adds both Darcy resistance and Forchheimer inertial drag.

To render this formidable system tractable, the researchers applied classical similarity transformations that collapse the governing partial differential equations into a set of coupled nonlinear ordinary differential equations for five dimensionless profiles: the axial velocity, the micro-rotation or angular velocity, the temperature, the nanoparticle concentration, and the density of motile microorganisms. Each equation bristles with dimensionless groups, from the Casson parameter and magnetic parameter to the Eckert number measuring viscous heating, the Brownian motion and thermophoresis parameters from Buongiorno’s model, the Lewis and Schmidt numbers governing diffusion, and the Peclet number describing how strongly microorganisms swim relative to how fast they diffuse. Five distinct slip parameters modify the boundary conditions at the wall, allowing velocity, temperature, concentration, and microorganism density to jump discontinuously at the fluid-solid interface, a realistic feature at microscale that almost no previous study has treated all at once.

The numerical solution was carried out with MATLAB’s bvp5c solver, which implements a three-stage Lobatto IIIA collocation formula specially designed for boundary value problems. The higher-order equations were first reduced to a system of first-order equations, and a shooting method was used to satisfy the conditions at both the wall and the far field. To verify the code, the team compared their computed skin friction values against previously published results for the micropolar Casson problem, finding agreement to six or more significant figures across magnetic parameter values ranging from 1 to 1000. This benchmark, they argue, establishes the reliability of the framework before it was extended to the full multi-physics case.

The machine learning component is where the work takes a distinctly modern turn. The researchers generated training data from the bvp5c solutions and built a multilayer perceptron artificial neural network trained with the Levenberg-Marquardt backpropagation algorithm, a method prized for its rapid and stable convergence on nonlinear least-squares problems. Eight input physical parameters fed the network, which predicted all five flow profiles simultaneously. The data were split into training, validation, and testing subsets, and the team tracked mean squared error across epochs, observing steady decline until optimal validation performance was reached after as few as nine epochs in some scenarios. Error histograms showed errors tightly clustered around zero, and regression plots placed every data point from every subset on the equality line, with coefficients of determination equal to one. Absolute errors between the network predictions and the numerical reference solutions fell between ten to the minus four and ten to the minus five, a level of fidelity the authors say confirms the network as a trustworthy surrogate for the full solver.

The physics that emerges from the simulations is rich and often counterintuitive. Increasing the magnetic parameter slows the fluid dramatically, because the Lorentz force opposes motion and thins the momentum boundary layer, a result with direct implications for controlling cooling rates in metallurgical and manufacturing processes. A stronger porous matrix resistance similarly drains momentum from the flow. By contrast, all three buoyancy mechanisms accelerate the fluid: raising the thermal, solutal, or microorganism Grashof number curves the velocity profile upward as buoyant plumes near the heated, concentrated, or bioactively populated wall drive stronger convection. The spin gradient parameter sharpens angular velocity near the wall, while the microinertia parameter suppresses it, reflecting the resistance of suspended elements to changing their rotational state.

On the thermal side, the Eckert number raises fluid temperature by converting kinetic energy into heat through viscous dissipation, thickening the thermal boundary layer. Both Brownian motion, which randomizes nanoparticle trajectories and intensifies fluid-particle collisions, and thermophoresis, which drives nanoparticles from hot regions toward cold ones, elevate the temperature profile. Yet the two mechanisms split on concentration: Brownian motion flattens the nanoparticle concentration profile by spreading particles more uniformly, while thermophoresis piles particles up near the cold far field and steepens the concentration boundary layer. For the microorganisms, a higher Peclet number means swimming transport overwhelms diffusion, pushing cells away from the wall and lowering their density in the boundary layer. Slip effects act as universal suppressors: velocity, thermal, solutal, and microorganism slip each weaken the transfer of momentum, heat, mass, or cells from the stretching surface into the fluid, thinning the corresponding boundary layers.

The engineering quantities tell a similarly coherent story. The Forchheimer number, which quantifies nonlinear inertial drag in the porous medium, increases skin friction by steepening the velocity gradient at the wall. A heat source thickens the thermal boundary layer and reduces the Nusselt number, while a heat sink does the opposite, extracting heat and boosting heat transfer rates. Higher Schmidt numbers thin the concentration boundary layer and raise the Sherwood number, and higher Lewis numbers reduce microorganism density at the surface. Together, these trends offer designers of thermal systems a tunable menu: magnetic fields and porous matrices to tame runaway flows, heat sinks to maximize cooling, and slip engineering to modulate transport at microstructured walls.

The authors are candid about the limits of the work. The model is laminar, steady, and incompressible, and no experimental validation or real-time physical implementation was attempted. Future extensions, they suggest, could incorporate more complex geometries, variable material properties, and more advanced machine learning architectures to push predictive accuracy higher while cutting computational cost. Even so, the study stands as a striking demonstration that neural networks can now serve as near-perfect surrogates for some of the most forbidding equation systems in fluid mechanics, potentially allowing engineers to explore design spaces that brute-force numerical simulation would take orders of magnitude longer to map. For industries that live and die by heat transfer, from electronics cooling to biomedical devices, that capability may prove as consequential as the fluid physics itself.

Subject of Research: Numerical and machine learning analysis of MHD Casson-micropolar nanofluid flow with bioconvection over a stretching sheet

Article Title: Numerical and machine learning investigation of dual-material fluid flow behavior under thermal and solutal rheological mechanisms

Article References: Muhammad, S., Khan, M., Muhammad, Rehman, G., Atif, H. M., & Hussain, S. M. (2026). Numerical and machine learning investigation of dual-material fluid flow behavior under thermal and solutal rheological mechanisms. Results in Engineering, 32, Article 113022. https://doi.org/10.1016/j.rineng.2026.113022

Image Credits: AI Generated

DOI: 10.1016/j.rineng.2026.113022

Keywords: Casson fluid, micropolar fluid, magnetohydrodynamics, nanofluid, bioconvection, artificial neural network, Buongiorno model, thermal radiation, stretching sheet, boundary layer, Levenberg-Marquardt, porous medium

Casson fluids occupy a special place in rheology. Unlike simple Newtonian liquids such as water, they possess a yield stress: below a critical applied stress they behave like solids, and only above that threshold do they flow. Blood is the canonical example, but polymer melts, printing inks and industrial suspensions also fall into this class. The Casson parameter in the new model quantifies how strongly the yield-stress behavior manifests, and the researchers found that raising it lets the fluid flow more freely along the stretching sheet, accelerating both the primary streamwise velocity and, to a much lesser degree, the secondary cross-flow induced by Hall currents. The secondary component barely responds to yield-stress changes, a subtle but telling sign that Casson rheology dominates the direction aligned with the stretching motion while the magnetically driven cross-flow marches to a different physical drumbeat.

The magnetic side of the problem introduces its own drama. A transverse magnetic field applied perpendicular to the sheet generates a Lorentz force that resists fluid motion, and the team captured this through a Hall-modified magnetic parameter that accounts for the fact that, under strong fields, the classical Ohm’s law breaks down. As this parameter grows, both velocity components decline because the magnetic damping force actively opposes forward motion. Yet the same magnetic field heats the fluid through Joule heating, the electrical resistance of the conducting nanofluid converting electromagnetic work into thermal energy. The result is a striking asymmetry: magnetic fields slow the flow while warming it, a combination that designers of electromagnetic pumping systems and metallurgical casting operations must weigh carefully.

The paper’s most counterintuitive finding concerns thermal radiation. Modeled through the Rosseland diffusion approximation for optically thick media, radiation increases the effective thermal diffusivity of the fluid, thickening the thermal boundary layer and raising temperatures throughout. But the tabulated Nusselt numbers tell a paradoxical story: as the radiation parameter climbs from 0.1 to 2, the wall heat transfer rate falls from 0.42871 to 0.1945. Radiation spreads thermal energy through the bulk of the fluid, flattening the temperature gradient right at the wall, and it is that near-wall gradient that determines how much heat actually escapes the surface. For engineers designing combustion chambers, solar receivers or glass manufacturing lines, where radiative effects dominate, this means that simply increasing radiative transport can quietly undermine heat removal, and mitigation strategies such as low-emissivity coatings or enhanced forced convection become essential.

Viscous dissipation compounds the thermal challenge. Quantified by the Eckert number, it represents the irreversible conversion of mechanical energy into heat through internal friction, and it grows severe in high-speed flows, polymer extrusion and microfluidic channels. In the new calculations, raising the Eckert number from zero to 0.5 drives the Nusselt number down from 0.35639 to 0.19763, because the internally generated heat erodes the temperature difference that drives wall heat transfer. In extreme cases, dissipation can push fluid temperatures above the wall temperature itself, creating thermal overshoots that threaten material integrity. The study’s quantitative tables give designers concrete numbers for anticipating when these overshoots become dangerous.

Cross-diffusion effects, the reciprocal coupling between temperature and concentration fields, add another layer of complexity. The Soret effect describes how temperature gradients drive mass flux, while the Dufour effect describes how concentration gradients drive heat flux. In the energy equation, the Dufour term acts as an additional heat sink in the boundary layer, and the results show that increasing the Dufour number from zero to 0.5 boosts the Nusselt number from 0.29273 to 0.4523, a genuinely beneficial coupling that enhances heat removal beyond what conduction and convection alone can achieve. The Soret effect, by contrast, slightly impairs mass transfer, nudging the Sherwood number downward as it rises, while modestly elevating concentration profiles through thermodiffusion. The authors note that these effects become order-unity corrections in high-gradient environments such as chemical vapor deposition, isotope separation and geothermal energy extraction, and cannot be neglected there even though they are secondary in milder settings.

Chemical reactions with Arrhenius activation energy intensify the thermal-concentration coupling because reaction rates depend exponentially on temperature, so small thermal fluctuations can produce order-of-magnitude swings in reaction kinetics. This nonlinearity matters enormously in catalytic reactors, combustion chambers and reactive separation processes. Meanwhile, the bioconvection component of the model captures the behavior of gyrotactic microorganisms, which swim in directed ways in response to gravitational and viscous torques and migrate toward nutrient-rich regions. Their collective motion generates density gradients that induce convective flow, and the study shows that the bioconvection Lewis number and Peclet number exert dramatic control over microorganism distributions. As the Lewis number climbs from 0.1 to 2, the motile microorganism density number rises from 0.30252 to 1.0326, reflecting how low-motility organisms become sharply confined near surfaces, potentially triggering bioconvective instabilities marked by collective upwelling.

The numerical machinery behind these findings is rigorous. The team transformed the full set of coupled nonlinear partial differential equations into ordinary differential equations via similarity transformations, then solved them with a fourth-order Runge–Kutta-Fehlberg method paired with a shooting technique and Newton–Raphson iteration, converging to a tolerance of ten to the minus six. Crucially, they validated the code against an exact analytical solution available in a simplified limiting case, achieving agreement to better than one part in ten thousand across a range of Casson parameters, and confirmed that their computational domain and far-field residuals were fully converged. This kind of benchmark-anchored validation lends real credibility to the parametric tables that form the study’s practical core.

Those tables reveal a clear hierarchy of influence that could streamline engineering design. Momentum transport is governed chiefly by the Casson parameter, magnetic damping and porous-medium permeability, with the permeability parameter showing nonmonotonic behavior in skin friction as Darcy resistance first concentrates momentum near the wall and then slows the flow entirely. Thermal transport is dominated by the Prandtl number, which at a fixed height cuts dimensionless temperature by 82 percent as it rises from 0.71 to 2.0, followed by radiation and dissipation. Mass transfer hinges on the Schmidt number, which lifts the Sherwood number from 0.1968 to 0.59234 across its tested range, while bioconvection is steered by the Lewis and Peclet numbers. Slip conditions and cross-diffusion act as secondary but non-negligible corrections, important chiefly in microscale and rarefied applications where the no-slip assumption fails.

The practical reach of the work spans heat exchangers, electronic cooling, nuclear reactor thermal management, biomedical therapies, bioreactors, wastewater treatment and drug delivery systems, wherever simultaneous control of momentum, heat, mass and living microorganisms is required. In algae photobioreactors, for instance, maintaining moderate Lewis and Peclet numbers prevents excessive wall accumulation that would shade cells and form biofilms, while microbial fuel cells benefit from the opposite regime, maximizing electrode colonization. The authors caution that their model assumes constant thermophysical properties, steady forcing and classical Fourier heat conduction, and they flag unsteady effects, nanoparticle aggregation and experimental validation as the natural next steps. Even so, by quantifying how radiation, dissipation, cross-diffusion, magnetism and bioconvection compete within a single validated framework, the study offers engineers a rational map for navigating some of the most tangled multiphysics in modern thermal design.

Subject of Research: Thermal radiation and coupled transport mechanisms in bioconvective Casson nanofluid flow over a stretching sheet

Article Title: Impact of thermal radiation in bioconvection Casson nanofluid flow with cross-diffusion

Article References: Madhusudanan, V., Guerrini, L., Iranian, D., R., R., Khan, I., & Ozsahin, I. (2026). Impact of thermal radiation in bioconvection Casson nanofluid flow with cross-diffusion. Results in Chemistry, 31, Article 103888. https://doi.org/10.1016/j.rechem.2026.103888

Image Credits: AI Generated

DOI: 10.1016/j.rechem.2026.103888

Keywords: Casson nanofluid, thermal radiation, bioconvection, Soret-Dufour, magnetohydrodynamics, Hall current, viscous dissipation, stretching sheet, boundary layer, Nusselt number, gyrotactic microorganisms, heat transfer