An assessment on the performance of sub-grid scale models of large eddy simulation in modeling bubbly flows

• Published in: Powder technology Vol. 374; pp. 470 - 481
• Main Authors: Liu, Zhongqiu, Wu, Yingdong, Li, Baokuan
• Format: Journal Article
• Language: English
• Published: Lausanne Elsevier B.V 01.09.2020; Elsevier BV
• Subjects: Adaptiveness; Bubbly flow; Computational fluid dynamics; Dynamic models; Energy transfer; High frequencies; Kinetic energy; Large eddy simulation; Reynolds stress; Scale models; SGS models; SGS turbulent dissipation; SGS turbulent kinetic energy; Turbulence; Turbulent flow; Void fraction; Vortices; SGS turbulent dissipation; Large eddy simulation; Bubbly flow; SGS turbulent kinetic energy; SGS models
• ISSN: 0032-5910; 1873-328X
• DOI: 10.1016/j.powtec.2020.07.055

An Euler-Euler large eddy simulation (LES) of transient gas-liquid turbulence flow in a bubble column was conducted to investigate the adaptiveness of various sub-grid scale (SGS) models. The performance of three SGS models, respectively the standard Smagorinsky model, the wall-adapting local eddy (WALE) model and the Dynamic model, on the predictions of bubbly flow was studied. Comparison between SGS models for time-averaged and transient velocities, gas void fraction, power spectrum density (PSD), SGS turbulent dissipation, SGS turbulent kinetic energy, and Reynolds stress is conducted. The transient multi-scale irregular vortex structures, fluctuation characteristics of liquid motion and spiral floating motion of bubbles can be well captured by LES. The Dynamic model and WALE model with CW = 0.325 give reasonably better agreement with the time-averaged experimental data. The transient fluctuation behavior of gas-liquid turbulent flow can be correctly captured by the current Euler-Euler LES model, whatever the SGS model. The classical −5/3 law of PSD of vertical liquid velocity is hold properly in the low frequency region. The slope of PSD is around −25/3 at the high frequency region. The Smagorinsky model and WALE model are sensitive to the model constant, at least for the energy transfer of high frequency. The molecular dissipation is negligible compared to the SGS turbulent dissipation. The Dynamic model gives higher predictions for fluctuation amplitude and time-averaged values of both SGS turbulent dissipation rate profile and SGS turbulent kinetic energy profile than other SGS models. A certain calculation time is needed to reach the dynamic stable state in LES, and the specific time is related to the scale of the calculation model. [Display omitted] •Adaptiveness of various SGS models in LES for two-phase turbulent flow is studied.•First to study SGS turbulent dissipation, SGS turbulent kinetic energy, and Reynolds stress tensor for bubbly flow.•A certain calculation time is needed to reach dynamic stability process in LES.