🔬 physics

Modeling Radiative Magnetohydrodynamic Mixed Convection in Second-Grade Viscoelastic Flow over a Continuously Stretching Cylinder

This study investigates the radiative magnetohydrodynamic mixed convection of a second-grade viscoelastic fluid over a continuously stretching cylinder by addressing a unique boundary condition mismatch through a mass conservation constraint and utilizing 3D surface topologies to reveal critical non-linear coupling effects on heat and mass transfer for advanced thermal management applications.

Original authors: Abdullah Muhammad, Azhar Ali, Aftab Alam, Ibrar Hussian

Published 2026-09-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Abdullah Muhammad, Azhar Ali, Aftab Alam, Ibrar Hussian

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the world of industrial manufacturing, many materials do not behave like water or air. Instead of flowing smoothly and predictably, substances like melted plastic, thick paints, and certain biological fluids possess a unique memory. When stretched or squeezed, they resist and then slowly return to their original shape, a property known as viscoelasticity. Understanding how these complex fluids move and transfer heat is critical for industries that produce everything from synthetic fibers to medical devices. When these fluids are subjected to magnetic fields or intense heat, their behavior becomes even more intricate, creating a challenge for engineers who need to design efficient cooling systems or precise manufacturing processes. This is where the study of fluid dynamics meets the practical demands of thermal management, requiring a deep look at how heat, motion, and magnetic forces interact within these non-standard liquids.

A team of researchers has turned its attention to a specific and difficult scenario: a long, cylindrical object being continuously stretched while surrounded by a fluid that is both viscoelastic and exposed to magnetic forces and thermal radiation. While scientists have long studied how fluids flow over flat surfaces, such as a sheet of plastic being pulled, the real world often involves curved objects like wires, fibers, or pipes. The curvature of a cylinder changes how the fluid clings to the surface and how heat spreads through it, making flat-surface models inaccurate for many industrial applications. The researchers set out to model this specific situation, incorporating the effects of a magnetic field, which can slow down the flow of charged particles in the fluid, and thermal radiation, which describes how heat travels as energy waves. They also considered how the fluid generates its own heat and how its ability to conduct heat changes as it gets hotter.

To solve this complex problem, the team translated the physical laws governing the fluid's motion into a set of mathematical relationships. They faced a significant hurdle: the equations describing the fluid's movement were so complex that they had more layers of detail than the standard rules for starting and ending the calculation allowed. To fix this, the researchers introduced an extra condition based on a physically motivated asymptotic decay constraint, ensuring strict index consistency and numerical closure within the solver algorithm. They then used two different powerful computational methods to solve these equations. One method provided a step-by-step analytical solution, while the other used a numerical approach to calculate precise values. By comparing the results from both methods, they confirmed that their model was robust and accurate, with both approaches agreeing perfectly on the behavior of the fluid.

The study revealed several key insights into how this fluid behaves under different conditions. When the researchers increased the strength of the magnetic field, the fluid's speed decreased. This happens because the magnetic field exerts a force that acts like a brake on the moving fluid, creating more resistance. Conversely, when the fluid's own elastic properties were strengthened, the fluid moved faster and the layer of fluid affected by the stretching surface became thicker. The curvature of the cylinder also played a surprising role; as the cylinder became smaller and more curved, the resistance near the surface changed, altering how the fluid flowed and how heat was distributed. The researchers found that increasing the radiation parameter, which represents the intensity of heat transfer via radiation, caused the temperature of the fluid to rise significantly, as more thermal energy was absorbed. Similarly, when the fluid was allowed to generate more internal heat, the temperature profile thickened, indicating that the fluid retained more energy.

Perhaps the most significant contribution of this work is not just the specific numbers, but the way the results were visualized. Instead of looking at simple two-dimensional lines that show how one variable changes at a time, the researchers mapped the results onto three-dimensional surfaces. This approach allowed them to see how multiple factors, such as magnetic strength and radiation, interacted simultaneously. These 3D maps revealed critical thresholds where the behavior of the fluid changed dramatically, details that would have been hidden in traditional graphs. For instance, they could see exactly how the friction on the cylinder's surface and the rate of heat transfer shifted as they adjusted the magnetic field and the curvature of the cylinder together.

The findings suggest that this model is highly reliable for optimizing real-world engineering systems. The ability to predict how a viscoelastic fluid will behave under magnetic and thermal stress is valuable for designing better cooling systems for electronics, improving the efficiency of nuclear reactors, and refining the processes used in polymer extrusion and metalworking. By validating their model against established data and using dual methods to ensure accuracy, the researchers have provided a solid foundation for future work. Their work confirms that accounting for the curvature of the object and the complex interplay of magnetic and thermal forces is essential for accurate predictions in advanced thermal management systems, offering a clearer path for engineers to design more efficient and effective industrial processes.

Technical Summary: Modeling Radiative Magnetohydrodynamic Mixed Convection in Second-Grade Viscoelastic Flow over a Continuously Stretching Cylinder

Problem Statement and Motivation
This study investigates the steady, two-dimensional, laminar flow of a second-grade viscoelastic fluid over a continuously stretching cylinder. The flow is driven by mixed convection and is subjected to magnetohydrodynamic (MHD) effects, thermal radiation, and internal heat/mass generation. Unlike previous studies that predominantly focus on flat stretching sheets, this research addresses the curvature effects inherent to cylindrical geometries, which are more representative of industrial processes such as polymer extrusion, wire coating, and fiber spinning.

The governing physics involve a non-Newtonian fluid where the stress tensor includes material constants for viscoelasticity. The model incorporates the Boussinesq approximation for buoyancy, variable thermal conductivity dependent on temperature, and a radiative heat flux term represented by $16\sigma T^3_\infty / (3(\rho C_p)k) \partial^2 T / \partial r^2$ in the energy equation. A specific mathematical challenge identified in the formulation is that the resulting fourth-order nonlinear momentum ordinary differential equation (ODE) is subject to only three boundary conditions, creating a uniqueness issue.

Methodology
The authors employ a dual-approach strategy to solve the coupled system of nonlinear ODEs derived from the continuity, momentum, energy, and concentration equations via similarity transformations:

Homotopy Analysis Method (HAM): A semi-analytical technique used to obtain series solutions. The authors utilize specific initial guesses and linear operators to ensure convergence. Convergence analysis is performed using $h$ -curves to determine the allowable range of auxiliary parameters ( $h_f, h_\theta, h_\phi$ ) for the series solutions.
Numerical Solution (bvp4c): The system is solved numerically using MATLAB's bvp4c solver. To address the boundary condition mismatch (fourth-order ODE with three conditions), the authors introduce an additional asymptotic decay constraint, $f''(\infty) = 0$ , representing reasonable decay at infinity, following the approach of Sadeghy and Sharifi. This ensures the system order matches the number of boundary conditions.

The results from both methods are cross-validated to ensure accuracy. The study also benchmarks its findings against previously published classical results to verify the model's reliability.

Key Contributions and Novelty

Geometric Realism: The study shifts the focus from flat sheets to a stretching cylinder, explicitly modeling curvature effects which alter the boundary layer dynamics.
Mathematical Resolution: The paper addresses the uniqueness challenge of the fourth-order momentum ODE by implementing an asymptotic decay constraint ( $f''(\infty) = 0$ ) to close the system for numerical solution, following the approach of Sadeghy and Sharifi.
Multi-Dimensional Visualization: Moving beyond traditional 2D parametric plots, the authors map the simultaneous impacts of governing parameters using 3D surface topologies. This approach captures non-linear coupling effects on the skin friction coefficient, Nusselt number, and Sherwood number that are often obscured in one-dimensional analysis.
Comprehensive Parameter Study: The model integrates MHD, radiation, mixed convection, and variable thermal conductivity for a second-grade fluid, filling a gap in literature regarding the combined influence of these factors on cylindrical geometries.

Results and Discussion
The study analyzes the effects of various dimensionless parameters on velocity ( $f'$ ), temperature ( $\theta$ ), and concentration ( $\phi$ ) profiles:

Velocity Profile:
- Permeability ( $\phi_p$ ): Increasing the impermeability parameter $\phi_p$ decreases the velocity profile due to increased flow resistance.
- Magnetic Parameter ( $M_p$ ): Increasing $M_p$ suppresses the velocity profile due to the Lorentz force acting as a drag.
- Curvature Parameter ( $\gamma$ ): Near the cylinder surface, increasing $\gamma$ (decreasing cylinder radius) reduces the boundary layer thickness and velocity due to higher viscous resistance. Farther from the surface, the trend reverses.
- Mixed Convection ( $\lambda$ ): Increasing $\lambda$ reduces the velocity and momentum boundary layer thickness.
- Viscoelastic Parameter ( $\beta$ ): Increasing $\beta$ (associated with reduced viscosity) increases both the velocity and the momentum boundary layer thickness.
Temperature Profile:
- Curvature ( $\gamma$ ): Near the surface, higher $\gamma$ reduces temperature distribution.
- Variable Thermal Conductivity ( $\epsilon$ ): Temperature distribution decreases as $\epsilon$ decreases.
- Radiation ( $R_p$ ) and Heat Generation ( $Q_p$ ): Both parameters lead to an increase in the temperature profile and thermal boundary layer thickness.
- Prandtl Number ( $Pr$ ): Increasing $Pr$ reduces thermal diffusivity, leading to a thinner thermal boundary layer and lower temperature profiles.
Concentration Profile:
- Schmidt Number ( $Sc$ ): Increasing $Sc$ (indicating lower mass diffusivity) reduces the concentration profile and boundary layer thickness.
- Curvature ( $\gamma$ ): Similar to temperature, the concentration behavior near the surface is inversely related to the curvature parameter.
Engineering Parameters:
- Skin Friction ( $C_f$ ): The 3D analysis shows that skin friction is negatively affected by the mixed convection parameter, indicating enhanced drag.
- Nusselt Number ( $Nu$ ): Heat transfer rates increase with higher radiation parameters.
- Sherwood Number ( $Sh$ ): Mass transfer rates decrease as the Schmidt number increases.

Significance and Claims
The paper claims that the proposed model offers a robust and physically reliable framework for analyzing complex viscoelastic flows in cylindrical geometries. The primary significance lies in:

Validation: The excellent agreement between the HAM analytical solutions and the bvp4c numerical results, along with consistency with classical literature, confirms the mathematical and computational accuracy of the model.
Practical Application: The findings are presented as having potential significance for advanced thermal management systems. The authors specifically cite applications in MHD generators, nuclear and solar thermal energy systems, metallurgical processing, polymer extrusion, electronic cooling modules, and high-performance heat exchangers.
Methodological Rigor: The successful implementation of the asymptotic decay constraint to resolve the boundary condition mismatch provides a methodological reference for similar high-order ODE problems in fluid dynamics.

The study concludes that the model serves as a rigorous benchmark for future research on complex viscoelastic flows and offers predictive capabilities for optimizing engineering designs involving coupled transport phenomena.

Technical Summary: Modeling Radiative Magnetohydrodynamic Mixed Convection in Second-Grade Viscoelastic Flow over a Continuously Stretching Cylinder

More like this