TY - JOUR
T1 - Free Surface Motion of a Liquid Pool with Isothermal Sidewalls as a Benchmark for Marangoni Convection Problems
AU - Ciccotosto, Bruce E.
AU - Brooks, Caleb S.
N1 - This research was funded by the US Department of Energy, as part of a University Nuclear Leadership Program Graduate Fellowship under grant number NE0009076.
This material is based upon work supported under a Department of Energy, Office of Nuclear Energy, Integrated University Program Graduate Fellowship. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the authors and do not necessarily reflect the views of the Department of Energy Office of Nuclear Energy.
PY - 2023/10
Y1 - 2023/10
N2 - In single phase flows, benchmarks like the lid driven cavity have become recognized as fundamental tests for newly developed computational fluid dynamics, CFD, codes. For multiphase free surface flows with variable surface tension, the presently studied pool with isothermal sidewalls is suggested as it is the simplest domain where Marangoni effects can dominate. It was also chosen due to its strange sensitivity to the initial setup which is discussed at length from a chosen number of ’scenarios’. It was found that the fluid interface can reverse deformation by a change in the top boundary condition, the liquid equation of state, and the gravity level. For the top boundary condition, this reversal is due to vapor expansion within the closed volume, creating an additional convection mechanism. Not only does the interface reverse, but the peak height changes by more than an order of magnitude at the same Marangoni number. When including gravity, the peak velocity can increase significantly, but it can also cause a decrease when done in combination with a change in the top wall boundary condition. Finally, thermal expansion of the liquid phase causes the peak velocity to be reduced, with additional reductions from the gravity and top wall condition. The differences in each scenario could lead to significant errors in analyzing a practical application of Marangoni flows. Therefore, it is important to demonstrate that a new CFD code can not only resolve Marangoni convection, but also has the capability to resolve the scenario most relevant to the application at hand.
AB - In single phase flows, benchmarks like the lid driven cavity have become recognized as fundamental tests for newly developed computational fluid dynamics, CFD, codes. For multiphase free surface flows with variable surface tension, the presently studied pool with isothermal sidewalls is suggested as it is the simplest domain where Marangoni effects can dominate. It was also chosen due to its strange sensitivity to the initial setup which is discussed at length from a chosen number of ’scenarios’. It was found that the fluid interface can reverse deformation by a change in the top boundary condition, the liquid equation of state, and the gravity level. For the top boundary condition, this reversal is due to vapor expansion within the closed volume, creating an additional convection mechanism. Not only does the interface reverse, but the peak height changes by more than an order of magnitude at the same Marangoni number. When including gravity, the peak velocity can increase significantly, but it can also cause a decrease when done in combination with a change in the top wall boundary condition. Finally, thermal expansion of the liquid phase causes the peak velocity to be reduced, with additional reductions from the gravity and top wall condition. The differences in each scenario could lead to significant errors in analyzing a practical application of Marangoni flows. Therefore, it is important to demonstrate that a new CFD code can not only resolve Marangoni convection, but also has the capability to resolve the scenario most relevant to the application at hand.
KW - marangoni
KW - CFD
KW - multiphase
UR - http://www.scopus.com/inward/record.url?scp=85174062799&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85174062799&partnerID=8YFLogxK
U2 - 10.3390/en16196824
DO - 10.3390/en16196824
M3 - Article
SN - 1996-1073
VL - 16
JO - Energies
JF - Energies
IS - 19
M1 - 6824
ER -