Mixed convection in four-sided lid-driven sinusoidally heated porous cavity using stream function-vorticity formulation

This study presents the mixed convection inside a four-sided lid-driven square porous cavity whose right wall is maintained at a sinusoidal temperature condition, the left wall of the cavity is maintained at a cold temperature, while the top and the bottom walls are adiabatic. We have discussed two different cases depending upon the direction of the moving walls. Brinkmann-extended Darcy model is represented in terms of ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi $$\end{document} and ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document} using the stream function-vorticity formulation to simulate the momentum transfer in the porous medium. This formulation is used to solve the governing equations as a coupled system of equations which consists of the field variables, vorticity (ξ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\xi )$$\end{document}, stream function (ψ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\psi )$$\end{document}, and temperature (T). The velocity components (u, v) are derived from the stream function (ψ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\psi )$$\end{document} whereas the average Nusselt number is derived from temperature. The stability and consistency of the applied numerical scheme to the considered problem has been proven by matrix method. The numerical results are investigated by ranging the various dimensionless numbers such as Grashof number (103≤Gr≤105)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(10^3 \le { {Gr}} \le 10^5)$$\end{document}, Darcy number (10-1≤Da≤10-5)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(10^{-1} \le { {Da}} \le 10^{-5})$$\end{document}, Reynolds number (10≤Re≤1000)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(10 \le { {Re}} \le 1000)$$\end{document} and keeping the Prandtl number (Pr=0.7)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({ {Pr}}=0.7)$$\end{document} fixed.