Research Article | | Peer-Reviewed

Homotopy Perturbation Method for Analyzing the Effect of Viscous Dissipation on Steady Natural Convection Couette Flow with Convective Boundary Conditions

Received: 7 November 2024     Accepted: 28 November 2024     Published: 23 December 2024
Views:       Downloads:
Abstract

This research article presents an analytical study of convective flow in a vertical channel with convective boundary conditions. Because of the nonlinear nature of the governing energy and momentum equations, the homotopy perturbation method was employed. The effects of various physical parameters on temperature and velocity profiles are illustrated in Figures 2 to 9, and a comparison table is provided to validate the results. Notably, both temperature and velocity distributions increased with higher viscous dissipation. Furthermore, the velocity profile decreased with an increase in the Biot number, while the temperature profile adjacent to the plate increased as the Biot number grew. Shear stress also exhibited an upward trend with rising viscous dissipation. Finally, an increase in the Grashof number and Biot number is found to elevate the skin friction on both plates. The mean temperature is higher when air is used as the working fluid compared to mercury. To validate this study, the temperature and velocity results were compared with previously published work, showing excellent agreement. This confirms the efficiency of the Homotopy Perturbation Method in solving coupled and nonlinear system of differential equations. Additionally, it was observed that both temperature and velocity increase with a rise in the Prandtl number, attributed to the dominance of momentum diffusivity over thermal diffusivity.

Published in International Journal of Fluid Mechanics & Thermal Sciences (Volume 10, Issue 3)
DOI 10.11648/j.ijfmts.20241003.11
Page(s) 45-56
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2024. Published by Science Publishing Group

Keywords

Convective Boundary Condition, Couette Flow, Natural Convection, Viscous Dissipation, Homotopy Perturbation Method

1. Introduction
The study of natural convection in Couette flow has gained significant attention due to its applications in various fields, including engineering, geothermal energy, the petroleum industry, physics, nuclear power plants, and lubrication technology. Momohjimoh et al. observed a reduction in velocity boundary layer thickness with increasing Biot and Prandtl numbers. Ajibade et al. noted that greater internal conductive resistance diminishes buoyancy forces, resulting in reduced fluid thickness, temperature, and velocity. Agunbiade and Oyekunle reported that higher viscous dissipation enhances fluid temperature and velocity. Similarly, Tafida and Ajibade highlighted the positive influence of viscous dissipation on these parameters. Ajibade and Tafida found that velocity boundary layer thickness enhances with larger Biot numbers.
Astanina et al. demonstrated that variable viscosity in the working fluid effectively controls heat transfer and fluid flow. Hussain et al. confirmed that the rate of heat transfer improves with higher Prandtl numbers. Ajibade et al. revealed that hydrodynamic and thermodynamic distributions improve with increased viscous dissipation, while Ajibade and Tafida observed that viscous dissipation raises fluid temperature but tends to reduce velocity. Finally, Ajibade and Princely concluded that greater viscous dissipation leads to higher fluid temperature and velocity within the channel. Mehdi and Kourosh concluded that growing the generalized Prandtl number increases heat transfer while reducing skin friction, and also decreases the thickness of the temperature boundary layer. Muawiya and Tafida observed that the mass flux increases with higher viscous dissipation (Ec). Isa et al. reported that skin friction rises with an increase in the Grashof number, while the Nusselt number decreases. Hamza noted that the heat source parameter increases with a rise in the Biot number. Masthanaiah et al. indicated that higher values of viscous dissipation result in greater temperature generation. Tafida et al. concluded that velocity and temperature distribution develop as viscous dissipation increases.
Ajibade et al. observed that the velocity profile decreases adjacent to the cold plate as viscous dissipation increases. Liu et al. highlighted the superiority of their method over conventional approaches, emphasizing its flexibility in selecting the initial approximation. Abker successfully provided solutions for both linear and nonlinear equations by the method of homotopy perturbation (HPM). Abu-Zeid found that temperature rises with an increase in the Brinkman number. Adamu concluded that HPM is a robust and effective tool for solving coupled differential equations. Ajibade and Tafida demonstrated the effectiveness of their HPM-based solution, achieving excellent results. Nasution reported that HPM is suitable for solving nonlinear differential equation systems, as evidenced by simulations where HPM results closely matched numerical simulations using the 4th-order Runge-Kutta method. The main objective of this research is to apply the homotopy perturbation method (HPM) to examine the impact of viscous dissipation on steady natural convection Couette flow under convective boundary conditions. The governing equations for this flow are coupled and nonlinear, making them suitable for numerical schemes or approximate solution methods. Among these, the perturbation method is effective; however, its solutions are limited to small perturbation parameters. To address this limitation, the homotopy perturbation method (HPM) is introduced as an alternative. HPM offers rapid convergence, often requiring only a limited terms in the series solution to achieve high correctness. Given the nonlinearity and coupling of the governing equations, HPM is employed in this study to solve the energy and momentum equations.
2. Methodology
We examine a steady flow of an incompressible fluid with viscous dissipation occurring between two vertical parallel plates positioned at specific locations. and , where the hot and cold walls are maintained at uniform temperatures and , respectively. It is assumed that the flow occur in the direction, vertically upward along the plates, while the axis is perpendicular to the plates, as depicted in Figure 1. Given the infinite length of the plates, both the velocity and temperature fields are assumed to depend solely on the spatial coordinate. .
Figure 1. Schematic diagram of the problem.
The flow behavior and heat transfer in a vertical channel based on the Boussinesq’s approximation, are described by the following governing equations:
(1)
(2)
The corresponding boundary conditions for these equations are:
at
at (3)
As the quantities involved have different dimensions, we introduce appropriate dimensionless variables to transform the governing equations and boundary conditions into a dimensionless form. The dimensionless quantities used in equations (1) – (2) and boundary condition (3) are:
(4)
Using these definitions, equations (1) to (3) can be rewritten in the following non-dimensional forms:
(5)
(6)
The non-dimensional boundary conditions are:
at
at(7)
Method of Solution
To address the problem with the Homotopy Perturbation Method, we formulate a convex homotopy for the energy and momentum equations. As a result, the equations can be expressed in the following form without requiring an initial approximation:
(8)
such that
(9)
By putting equation (9) into equation (8), we have
(10)
By relating the coefficient of , we obtain
,(11)
(12)
(13)
(14)
The boundary conditions are transformed to:
(15)
Similarly, for equation (6), there is no initial approximation; thus, it can be transformed into:
(16)
such that
(17)
By substituting equation (17) into equation (16), we obtain:
(18)
By relating the coefficient of , we obtain
,(19)
(20)
(21)
(22)
The boundary conditions for these equations are transformed to:
(23)
Thus, by solving equations (11) and (19) and applying the boundary conditions and and we obtain the solutions
(24)
(25)
Thus, by solving equations (12) and (20) and applying the boundary conditions and and we obtain the solutions
(26)
(27)
Thus, by solving equations (13) and (21) and applying the boundary conditions and and we obtain the solutions
(28)
(29)
Thus, by solving equations (14) and (22) and applying the boundary conditions and and we obtain the solutions
(30)
(31)
Where,
The rate of heat transfer and skin friction and are expressed as follows:
To find the mass flux, , we have
To obtain the mean temperature,
,
(32)
3. Results and Discussion
This section presents a discussion of the results obtained. The analytical solutions for temperature and velocity are displayed in graphs and analyzed for various values of the monitoring parameters. Moreover, the calculated numerical values for the rate of heat transfer, skin friction, mass flux, and mean temperature are provided in tables for further evaluation.
The temperature and velocity profiles for different values of Eckert number (Ec) is displayed in Figures 2 and 3. Both distributions increase with higher Eckert numbers, as the fluid's kinetic energy surpasses the boundary layer enthalpy difference. This intensifies convection currents, reducing fluid density and resulting in higher temperature and velocity profiles within the channel. It is also noted that the effect of the Eckert number on the temperature distribution is more pronounced near the cold plate.
The effect of the thermal buoyancy ratio (Gr) is shown in Figures 4 and 5. It is observed that the temperature distribution increases near the heated plate as the Grashof number increases, while a decrease occurs near the cold plate with rising Grashof number. Additionally, the Grashof number has no noticeable effect on the temperature distribution at the center of the channel. In Figure 5, the velocity distribution increases across the channel with higher Gr values, which is attributed to the buoyancy forces overpowering the viscous forces. This causes fluid particles in more energetic regions to move toward less energetic areas.
Figures 6 and 7 illustrate the effects of the Prandtl number (Pr) on the fluid's temperature and velocity distributions. It is observed that both temperature and velocity increase as the Prandtl number rises. This is due to the dominance of momentum diffusivity over thermal diffusivity. Physically, as the momentum diffusivity rate surpasses the thermal diffusivity, the fluid's boundary layer thickness increases, which enhances the convection currents within the channel, leading to higher temperature and velocity distributions.
Figures 8 and 9 illustrate the effects of the Biot number (Bi) on the velocity and temperature distributions, respectively. Figure 8 shows that temperature distribution increases nearby the heated plate as the Biot number rises, while it decreases near the cold plate with an increasing Biot number. In Figure 9, the velocity profile is observed to decrease as the Biot number increases. This behavior is due to the dominance of internal conduction resistance over external convection resistance, which affects heat transfer to the fluid, consequently reducing its temperature and weakening fluid motion.
Table 1 presents the computed numerical values for the rate of heat transfer at both plates. It is evident from the table that, the rate of heat transfer decreases at both plates as viscous dissipation increases. Furthermore, an increase in both the Prandtl number and viscous dissipation results in a higher rate of heat transfer.
Table 2 shows the calculated numerical values of skin friction at both plates. It is observed that an increase in the Grashof number and Biot number inclines to raise the skin friction at both plates.
Table 3 presents the mass flux values. Mass flux increases as the Biot number and Grashof number rise, while it decreases with an increase in the Prandtl number.
Table 4 shows the mean temperature values. The mean temperature is higher when the working fluid is air compared to mercury. Additionally, the mean temperature increases with rising Grashof and Biot numbers.
Validation of results
To validate this study, we compared our results for temperature and velocity with those of Muawiya and Tafida , finding excellent agreement, as presented in Table 5. This confirms that the Homotopy Perturbation Method is an effective approach for solving coupled and nonlinear differential equations. Additionally, to further substantiate our findings, we extended the work of Muawiya and Tafida by introducing a different boundary condition, specifically the convective boundary condition. The comparison is outlined in Table 5.
4. Conclusion
This article examines the effect of viscous dissipation on steady natural convection Couette flow with convective boundary conditions, employing the homotopy perturbation method. The energy and momentum equations were solved, and the results were presented graphically and analyzed for various values of the governing parameters. The study concludes that:
Both temperature and velocity profiles rise with increasing viscous dissipation,
The velocity profile decreases as the Biot number increases,
A rise in the Prandtl number and viscous dissipation leads to an enhanced rate of heat transfer,
The mean temperature is greater when air is used as the working fluid compared to mercury.
Figure 2. Temperature profile for various values of .
Figure 3. Velocity profile for various values of .
Figure 4. Temperature profile for various values of .
Figure 5. Velocity profile for various values of .
Figure 6. Velocity profile for various values of .
Figure 7. Temperature profile for various values of .
Figure 8. Velocity profile for various values of .
Figure 9. Temperature profile for various values of .
Table 1. Calculated numerical values for the rate of heat transfer at both plates.

Table 2. Calculated numerical values for the skin friction at both plates.

Table 3. Calculated numerical values of mass flux .

Table 4. Calculated numerical values of mean temperature.

Table 5. Comparison of calculated numerical values between the present problem and that of Muawiya and Tafida .

Muawiya and Tafida

Present problem

Ec

Temperature

Velocity

Temperature

Velocity

Abbreviations

g

Acceleration to Due Gravity [ms-2]

β

Coefficient of Thermal Expansion [K-1}

h

Width of the Channel [m]

μ

Coefficient of Viscosity [Kgm-1s-1]

T*

Dimensional Fluid Temperature [K]

ν

Kinematic Viscosity [m2s-1]

Tw

Channel Wall Temperature [K]

Pr

Prandtl Number

T0

Temperature of the Ambience [K]

Gr

Grashof Number

T

Dimensionless Fluid Temperature

Ec

Eckert Number

u*

Dimensional Velocity [ms-1]

Bi

Biot Number

u

Dimensionless Velocity

cp

Specific Heat at Constant Pressure [m2s-2K-1]

U

Dimensional Velocity of the Moving Plate [ms-1]

ρ

Density of the Fluid [Kgm-3]

y*

Coordinate Perpendicular to the Plate [m]

α

Thermal Diffusivity of the Fluid [Kgm-3]

y

Dimensionless Co-ordinate Perpendicular to the Plate

Conflicts of Interest
The authors declare no conflicts of interest.
References
[1] Momohjimoh, S. A., Oyem, O. A., Momoh, O. S., Onojovwo, T,. Felix. Convective heat and casson nanofluid flow over a vertical plate with heat source. Fudma Journal of Sciences. (2023), 7(2); 9 – 18.
[2] Ajibade, A. O., Jeremiah, J. G., and Basant, K. J. Effect of darcy and viscous dissipation on natural convection flow in a vertical tube partially filled with porous material under convective boundary condition. International Journal of Applied Computational Mathematics. (2024), 10 – 84.
[3] Agunbiade, S. A., and Oyekunle, S. A. Effects of magnetic field on natural convective steady flow past an irregular vertical channel in the presence of viscous dissipation. FUW Trends in Science and Technology Journal. (2019), 4(3); 666 – 670.
[4] Tafida, M. K. and Ajibade, A. O. Variable viscosity and viscous dissipation effects on free convective Couette flow in a vertical channel: the homotopy perturbation method approach. (2020). The Journal of the Mathematical Association of Nigerian. 47(1), 161 – 175.
[5] Ajibade, A. O., Tafida, M. K. Viscous dissipation effect on steady natural convection Couette flow with convective boundary condition. International Journal of Non-linear Sciences and Numerical Simulation. (2022),
[6] Astanina, M. S., Mikhail, A. S. numerical study on natural convection of fluid with Temperature dependent viscosity. International Journal of Thermofluids. 2023, 100266.
[7] Hussain, A., Isah, B. Y., and Zayyanu, S. Y. Impression of some fluid flow possession on unsteady conductive free convective fluid on account to dissipative porous walls. Himalayan Journal of Applied Medical Sciences and Research. 4(7), (2022), 31-52.
[8] Ajibade, A. O., Ayuba, M. U., and Tafida, M. K. An analytical study on effects of viscous dissipation and suction/injection on a steady MHD natural convection Couette flow of heat generating/absorbing fluid. Advances in Mechanical Engineering. 13(5), (2021), 1-12.
[9] Ajibade, A. O. and Tafida, M. K. Viscous dissipation effect of on steady natural Convection Couette flow with convective boundary condition. International Journal of Non-Linear Sciences and Numerical Simulation (IJNSNS). (2022), 7: 105 – 116.
[10] Ajibade, A. O. and Princely, O. O. Effect of viscous dissipation on steady natural convection, heat and mass transfer in a vertical channel with variable viscosity and thermal conductivity. Journal of Heat and Mass Transfer Research. (2020). 7: 105 116.
[11] Mehdi, M. and Kourosh, J. Viscous dissipation effect in the free convection of Non- Newtonian fluid with heat generation or absorption effect on the vertical wavy surface. Hindawi Journal of Applied Mathematics ID: 7567981, 14 pages. (2021):
[12] Muawiya, H. U. and Tafida, M. K. Viscous dissipation effects on steady natural convection Couette flow in a vertical channel. Wudil Journal of Pure and Applied Sciences. (2019), 1(2), 169 – 179.
[13] Isa, B. U., Yale, I. D., Sarki, M. N., and Hamza, M. M. Effects of viscous dissipative fluid on Couette flow in a vertical channel due to Newtonian heating. International Journal of Science for Global Sustainability. (2024),
[14] Hamza, M. M. Navier slip and convective boundary condition effects on transient combined convection of dissipative fluid. The Transactions of the Nigerian Association Mathematical Physics. (2023), vol 19 (Jan – Dec issues).
[15] Masthanaiah, Y., Nainaru, T., Ijaz, K. M., Rushikesava, S. M., Bandar, M. F., Sherzod, S., Abdullaer, S., and Sayed, M. E. Impact of viscous dissipation and entropy generation on cold liquid via channel with porous medium by analytical analysis. Case Studies in Thermal Engineering. (2023), 47; 103059.
[16] Tafida, M. K., Ayuba, M. U., and Lawal, U. Magnetohydrodynamics free convection Couette flow and heat transfer through a vertical porous plate with heat generation and absorption effects. Dutse Journal of Pure and Applied Sciences. (2023), 9(2a); 133 - 145.
[17] Ajibade, A. O., Ayuba, M. U. and Tafida, M. K. An Analytical study on effects of viscous dissipation and suction/injection on a steady MHD convection Couette flow of heat generating/absorbing fluid. Advances in Mechanical Engineering. (2021). 13(5), 1 – 12.
[18] Liu, Z. J., Adamu, M. Y., Suleiman, E., and He, J. H. Hybridization of Homotopy Perturbation method and Laplace transformation for the partial differential equations. Thermal Science. (2017), 21; 1843 – 1846.
[19] Abker, A. H. Solution of the linear and nonlinear Schrodinger equations using Homotopy Perturbation method and Iterative method. American Journal of Engineering Research. (2017), 6(3); 107 – 114.
[20] Abou-Zeid, M. Homotopy Perturbation Method for Magnetohydrodynamics Non-Newtonian Nanofluid flow through a porous medium in eccentric annuli with peristalsis. Thermal Science. (2017), 21(5); 2069 – 2080.
[21] Adamu, A. Y. Parameterized Homotopy Perturbation method. Nonlinear Science Letter A. (2017), 8(2); 240 – 243.
[22] Ajibade, A. O. and Tafida, M. K. Effects of variability in viscosity and thermal conductivity Couette flow and heat transfer. Transactions of NAMP vol. 11, (January – June, 2020 Issue).
[23] Nasution, H. Exploring of homotopy perturbation method (HPM) for Solving Spread of Covid-19. Jambura Journal of Biomathematics. (2023), 4(2); 138 – 154.
Cite This Article
  • APA Style

    Kabir, T. M., Abdulganiyu, T. (2024). Homotopy Perturbation Method for Analyzing the Effect of Viscous Dissipation on Steady Natural Convection Couette Flow with Convective Boundary Conditions. International Journal of Fluid Mechanics & Thermal Sciences, 10(3), 45-56. https://doi.org/10.11648/j.ijfmts.20241003.11

    Copy | Download

    ACS Style

    Kabir, T. M.; Abdulganiyu, T. Homotopy Perturbation Method for Analyzing the Effect of Viscous Dissipation on Steady Natural Convection Couette Flow with Convective Boundary Conditions. Int. J. Fluid Mech. Therm. Sci. 2024, 10(3), 45-56. doi: 10.11648/j.ijfmts.20241003.11

    Copy | Download

    AMA Style

    Kabir TM, Abdulganiyu T. Homotopy Perturbation Method for Analyzing the Effect of Viscous Dissipation on Steady Natural Convection Couette Flow with Convective Boundary Conditions. Int J Fluid Mech Therm Sci. 2024;10(3):45-56. doi: 10.11648/j.ijfmts.20241003.11

    Copy | Download

  • @article{10.11648/j.ijfmts.20241003.11,
      author = {Tafida Mohammed Kabir and Tajuddeen Abdulganiyu},
      title = {Homotopy Perturbation Method for Analyzing the Effect of Viscous Dissipation on Steady Natural Convection Couette Flow with Convective Boundary Conditions
    },
      journal = {International Journal of Fluid Mechanics & Thermal Sciences},
      volume = {10},
      number = {3},
      pages = {45-56},
      doi = {10.11648/j.ijfmts.20241003.11},
      url = {https://doi.org/10.11648/j.ijfmts.20241003.11},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.ijfmts.20241003.11},
      abstract = {This research article presents an analytical study of convective flow in a vertical channel with convective boundary conditions. Because of the nonlinear nature of the governing energy and momentum equations, the homotopy perturbation method was employed. The effects of various physical parameters on temperature and velocity profiles are illustrated in Figures 2 to 9, and a comparison table is provided to validate the results. Notably, both temperature and velocity distributions increased with higher viscous dissipation. Furthermore, the velocity profile decreased with an increase in the Biot number, while the temperature profile adjacent to the plate increased as the Biot number grew. Shear stress also exhibited an upward trend with rising viscous dissipation. Finally, an increase in the Grashof number and Biot number is found to elevate the skin friction on both plates. The mean temperature is higher when air is used as the working fluid compared to mercury. To validate this study, the temperature and velocity results were compared with previously published work, showing excellent agreement. This confirms the efficiency of the Homotopy Perturbation Method in solving coupled and nonlinear system of differential equations. Additionally, it was observed that both temperature and velocity increase with a rise in the Prandtl number, attributed to the dominance of momentum diffusivity over thermal diffusivity.
    },
     year = {2024}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Homotopy Perturbation Method for Analyzing the Effect of Viscous Dissipation on Steady Natural Convection Couette Flow with Convective Boundary Conditions
    
    AU  - Tafida Mohammed Kabir
    AU  - Tajuddeen Abdulganiyu
    Y1  - 2024/12/23
    PY  - 2024
    N1  - https://doi.org/10.11648/j.ijfmts.20241003.11
    DO  - 10.11648/j.ijfmts.20241003.11
    T2  - International Journal of Fluid Mechanics & Thermal Sciences
    JF  - International Journal of Fluid Mechanics & Thermal Sciences
    JO  - International Journal of Fluid Mechanics & Thermal Sciences
    SP  - 45
    EP  - 56
    PB  - Science Publishing Group
    SN  - 2469-8113
    UR  - https://doi.org/10.11648/j.ijfmts.20241003.11
    AB  - This research article presents an analytical study of convective flow in a vertical channel with convective boundary conditions. Because of the nonlinear nature of the governing energy and momentum equations, the homotopy perturbation method was employed. The effects of various physical parameters on temperature and velocity profiles are illustrated in Figures 2 to 9, and a comparison table is provided to validate the results. Notably, both temperature and velocity distributions increased with higher viscous dissipation. Furthermore, the velocity profile decreased with an increase in the Biot number, while the temperature profile adjacent to the plate increased as the Biot number grew. Shear stress also exhibited an upward trend with rising viscous dissipation. Finally, an increase in the Grashof number and Biot number is found to elevate the skin friction on both plates. The mean temperature is higher when air is used as the working fluid compared to mercury. To validate this study, the temperature and velocity results were compared with previously published work, showing excellent agreement. This confirms the efficiency of the Homotopy Perturbation Method in solving coupled and nonlinear system of differential equations. Additionally, it was observed that both temperature and velocity increase with a rise in the Prandtl number, attributed to the dominance of momentum diffusivity over thermal diffusivity.
    
    VL  - 10
    IS  - 3
    ER  - 

    Copy | Download

Author Information