Unsteady, two-dimensional magnetohydrodynamic (MHD) analysis of Casson fluid flow in a porous cavity with heated cylindrical obstacles

Bai Mbye Cham,Shams-ul-Islam,M. Saleem,Shaiza Talib,Shafee Ahmad
DOI: https://doi.org/10.1063/5.0178827
IF: 1.697
2024-04-01
AIP Advances
Abstract:This research presents an analysis of entropy generation during natural convection in a porous medium using triangular heated cylindrical obstacles with equal spacing. The study consists of three cylindrical obstacles arranged in a triangular pattern. Each cylinder is uniformly spaced from its neighboring cylinders, creating equilateral triangles throughout the arrangement. All of these cylindrical obstacles are heated. The triangular arrangement guarantees an even distribution of obstacles across the experimental space. The governing equations, with entropy, are numerically solved using the finite element method. The study aims to investigate the interactions between several key elements in fluid dynamics: Casson fluid, magnetohydrodynamics, the Darcy–Forchheimer model, entropy, and natural convection. The goal is to gain insights into the individual behaviors of these elements and their interactions in combined systems. The results indicate that the Casson fluid parameter has an impact on the flow and heat transfer characteristics, while the Hartmann and Nusselt numbers exhibit control mechanisms for the intensity of natural convection and affect the patterns of isotherms, streamlines, and entropy.
materials science, multidisciplinary,nanoscience & nanotechnology,physics, applied
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to analyze the magnetohydrodynamics (MHD) behavior in the Casson fluid flow with a heated cylindrical obstacle in a porous cavity and its influence on entropy generation. Specifically, the research focuses on the following aspects: 1. **Unsteady two - dimensional magnetohydrodynamics (MHD) analysis**: Study the unsteady two - dimensional flow of Casson fluid in porous media, taking into account the magnetohydrodynamics effect. 2. **Casson fluid characteristics**: The Casson fluid is a non - Newtonian fluid, and its viscosity depends not only on the shear rate but also on the yield stress. The research explores the influence of Casson fluid parameters on fluid flow and heat transfer characteristics. 3. **Porous media model**: Use the Darcy - Forchheimer model to describe the fluid flow in porous media, taking into account the influence of linear and nonlinear resistances. 4. **Natural convection**: Study the natural convection phenomenon caused by temperature differences in the absence of external forced convection. 5. **Entropy generation analysis**: Evaluate the irreversibility and energy dissipation in the system by analyzing entropy generation, thereby optimizing the efficiency of the system. ### Specific problem statement - **Research background**: In the porous cavity, there are three heated cylindrical obstacles arranged at equal intervals, and these obstacles are evenly distributed throughout the experimental space. - **Research objective**: Solve the governing equations by numerical methods (finite element method) to study the flow and heat transfer characteristics of Casson fluid under the combined action of magnetohydrodynamics, Darcy - Forchheimer model, natural convection and entropy generation. - **Key parameters**: The key parameters considered in the research include Reynolds number (Ra), radiation parameter (Rd), Hartmann number (Ha), Eckert number (Ec), Casson fluid parameter (β), Forchheimer number (Γ), inverse Darcy number (γ) and Prandtl number (Pr). ### Mathematical model - **Continuity equation**: \[ \frac{\partial u}{\partial x}+\frac{\partial v}{\partial y} = 0 \] - **Momentum equation (x - direction)**: \[ \frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=-\frac{1}{\rho}\frac{\partial P}{\partial x}+\nu\left(1 + \frac{1}{\beta}\right)\nabla^{2}u-\nu\left(1 + \frac{1}{\beta}\right)\epsilon+\frac{Ku - F\epsilon+2\sqrt{K}|V|u-\sigma\beta_{0}^{2}\epsilon\rho u}{} \] - **Momentum equation (y - direction)**: \[ \frac{\partial v}{\partial t}+u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}=-\frac{1}{\rho}\frac{\partial P}{\partial y}+\nu\left(1 + \frac{1}{\beta}\right)\nabla^{2}v-\nu\left(1 + \frac{1}{\beta}\right)\epsilon+\frac{Kv - F\epsilon+2\sqrt{K}|V|v+g\beta(T - T_{c})}{} \] - **Energy equation**: \[ \frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y} \]