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