TY - JOUR
T1 - Side wall boundary effect on the Rayleigh–Taylor instability
AU - Yang, Junxiang
AU - Lee, Hyun Geun
AU - Kim, Junseok
N1 - Publisher Copyright:
© 2020 Elsevier Masson SAS
PY - 2021/1/1
Y1 - 2021/1/1
N2 - Rayleigh–Taylor instability (RTI) is a common phenomenon in daily life and in industrial application. Despite the long history of investigating the RTI with periodic boundary condition, there are only few studies for the RTI with no-slip side wall boundary condition. In fact, the RTI with the side wall boundary condition is more realistic physics. In this work, we investigate the side wall boundary effect on the RTI using a phase-field model. No-slip conditions are employed on all boundaries. The governing equations are the coupled Navier–Stokes (NS) and Cahn–Hilliard (CH) system. We temporally solve the NS equation using Chorin's projection method and the CH equation using Eyre's nonlinear splitting scheme. A multigrid method is employed to solve the discrete system of CH equation. By various numerical experiments, we find that the side wall will delays the evolution of bubble and spike, causes a longer interface length and complex interface shape. Besides that, some other parameters, such as Reynolds number and aspect ratio of domain etc., can also cause significantly different phase evolutions.
AB - Rayleigh–Taylor instability (RTI) is a common phenomenon in daily life and in industrial application. Despite the long history of investigating the RTI with periodic boundary condition, there are only few studies for the RTI with no-slip side wall boundary condition. In fact, the RTI with the side wall boundary condition is more realistic physics. In this work, we investigate the side wall boundary effect on the RTI using a phase-field model. No-slip conditions are employed on all boundaries. The governing equations are the coupled Navier–Stokes (NS) and Cahn–Hilliard (CH) system. We temporally solve the NS equation using Chorin's projection method and the CH equation using Eyre's nonlinear splitting scheme. A multigrid method is employed to solve the discrete system of CH equation. By various numerical experiments, we find that the side wall will delays the evolution of bubble and spike, causes a longer interface length and complex interface shape. Besides that, some other parameters, such as Reynolds number and aspect ratio of domain etc., can also cause significantly different phase evolutions.
KW - Phase-field model
KW - Projection method
KW - Rayleigh–Taylor instability
KW - Side wall boundary
UR - http://www.scopus.com/inward/record.url?scp=85094318049&partnerID=8YFLogxK
U2 - 10.1016/j.euromechflu.2020.10.001
DO - 10.1016/j.euromechflu.2020.10.001
M3 - Article
AN - SCOPUS:85094318049
SN - 0997-7546
VL - 85
SP - 361
EP - 374
JO - European Journal of Mechanics, B/Fluids
JF - European Journal of Mechanics, B/Fluids
ER -