Multiple steady states are investigated for natural convection of fluids in a square enclosure with non-isothermally hot bottom wall, isothermally cold side walls, and thermally insulated top wall. A robust computation scheme involving steady-state governing equations has been developed to compute the steady states as a function of Rayleigh number () for two different Prandtl numbers (Pr = 0.026 and 0.1). Penalty Galerkin finite element method with Newton–Raphson solver is employed for the solution of the governing equations, while the solution branches are initiated by varying initial guess to the Newton–Raphson solver. In this context, a dual-perturbation scheme involving perturbations of the boundary conditions and various process parameters has been designed leading to the rich spectrum of the symmetric and asymmetric solution branches for the current symmetric problem. It is found that multiple steady states occur beyond a critical value of Ra, which depends on the magnitude of Pr. In addition to the basic solution branch (corresponding to the solutions obtained via uniform initial guesses), nineteen additional solution branches (six symmetric and thirteen asymmetric) are obtained for Pr = 0.026, while four additional solution branches (two symmetric and two asymmetric) are obtained for Pr = 0.1. The solution branches are associated with a wide spectrum of flow structures (24 distinct types for Pr = 0.026 excluding the reflection symmetric mirror images of the asymmetric solutions), which are reported for the first time. The flow structures lead to various heating scenarios within the enclosure resulting in a significant variation of heat transfer rates (more than 50%). The current results are important for the practical applications. The spectrum of the possible scenarios revealed in this work can be pivotal to design the optimal processes based on the process requirement (targeted heating or enhanced heating rates).
Skip Nav Destination
Article navigation
October 2021
Research Article|
October 05 2021
Analysis of multiple steady states for natural convection of Newtonian fluids in a square enclosure
Special Collection:
Tribute to Frank M. White on his 88th Anniversary
Madhuchhanda Bhattacharya;
Madhuchhanda Bhattacharya
a)
1
Indian Institute of Technology Madras Campus, C2-5-4C
, Delhi Avenue, Chennai 600036, India
a)Authors to whom correspondence should be addressed: madhuchhanda_bhattacharya@yahoo.co.in and tanmay@iitm.ac.in
Search for other works by this author on:
Tanmay Basak
Tanmay Basak
a)
1
Indian Institute of Technology Madras Campus, C2-5-4C
, Delhi Avenue, Chennai 600036, India
2
Department of Chemical Engineering, Indian Institute of Technology
Madras, Chennai 600036, India
a)Authors to whom correspondence should be addressed: madhuchhanda_bhattacharya@yahoo.co.in and tanmay@iitm.ac.in
Search for other works by this author on:
a)Authors to whom correspondence should be addressed: madhuchhanda_bhattacharya@yahoo.co.in and tanmay@iitm.ac.in
Note: This paper is part of the special topic, Tribute to Frank M. White on his 88th Anniversary.
Physics of Fluids 33, 103605 (2021)
Article history
Received:
July 28 2021
Accepted:
September 08 2021
Citation
Madhuchhanda Bhattacharya, Tanmay Basak; Analysis of multiple steady states for natural convection of Newtonian fluids in a square enclosure. Physics of Fluids 1 October 2021; 33 (10): 103605. https://doi.org/10.1063/5.0064899
Download citation file:
Sign in
Don't already have an account? Register
Sign In
You could not be signed in. Please check your credentials and make sure you have an active account and try again.
Sign in via your Institution
Sign in via your InstitutionPay-Per-View Access
$40.00