气固两相湍流大涡模拟中DF模型有效性研究

doi:10.11949/0438-1157.20250474

摘要/Abstract

摘要：

两相湍流大涡模拟方法中缺失的亚网格尺度湍流信息会对离散相运动特性的预测造成较大误差。迄今为止，颗粒亚格子模型均是基于单向耦合条件构建的。为了提高气固两相湍流大涡模拟的模拟精度，结合基于小波离散误差的自适应网格算法，搭建了欧拉-拉格朗日框架双向耦合条件下的气固两相湍流模拟平台。在该平台下对衰减均匀各向同性湍流进行模拟，结果表明，基于差分过滤（DF）的颗粒所见流体亚格子速度模型的加入使中小惯性范围内的颗粒对径向分布函数相对误差减小约28%，径向相对速度相对误差减小约13%，能在高载荷下以较低的计算成本补偿离散相模拟精度，为该模型在大型工业反应器及流化床等化工设备实际工程应用的数值模拟中的发展及优化奠定了理论基础。

关键词: 两相流, 湍流, 数值模拟, 颗粒亚格子模型, 双向耦合

Abstract:

In large-eddy simulations (LES) of gas-solid two phase turbulent flows, the absence of subgrid-scale (SGS) turbulence information introduces significant inaccuracies in predicting the motion characteristics of the dispersed phase. To date, all particle subgrid models are based on one-way coupling conditions. With the adaptive mesh algorithm based on wavelet discrete errors, this study addresses this limitation by developing a novel gas-solid turbulence large eddy simulation platform under two-way coupled conditions. The subgrid velocity model seen by particles based on differential filter (DF) was then applied to the simulation of decaying gas-solid homogeneous isotropic turbulence. Comparative analysis of particle-pair statistics in pre-/post-model implementation simulations demonstrates that the proposed approach effectively restores dispersed-phase accuracy while validating its efficacy under high-volume-loading two-way coupling conditions. The inclusion of the DF model achieves enhanced simulation accuracy for the dispersed phase while maintaining relatively low computational costs. It reduces the relative error of the particle pairs radial distribution function with moderate inertia by about 28%, and the relative error of the radial relative velocity by about 13%. This study establishes critical theoretical foundations for advancing engineering applications in industrial particulate systems, particularly enhancing predictive capabilities for fluidized bed dynamics and multiphase reactor optimization.

Key words: two-phase flow, turbulence flow, numerical simulation, particle sub-grid model, two-way coupled

中图分类号:

O 375.5

谭昭轶, 栗晶, 覃羿翔, 柳朝晖. 气固两相湍流大涡模拟中DF模型有效性研究[J]. 化工学报, 2025, 76(11): 5677-5686.

Zhaoyi TAN, Jing LI, Yixiang QIN, Zhaohui LIU. Study on the effectiveness of DF model in large eddy simulations of gas-solid two phase turbulence[J]. CIESC Journal, 2025, 76(11): 5677-5686.

图/表 10

表1 衰减湍流初始流场特征统计量

Table 1 Characteristic statistics of the initial decaying turbulence

特征统计量	数值
动力黏度ν	0.004
Kolmogorov尺度 $η$	0.016
耗散率ε	0.945
湍动能κ	2.389
泰勒Reynolds数Re_λ	100.3

表1 衰减湍流初始流场特征统计量

Table 1 Characteristic statistics of the initial decaying turbulence

特征统计量	数值
动力黏度ν	0.004
Kolmogorov尺度 $η$	0.016
耗散率ε	0.945
湍动能κ	2.389
泰勒Reynolds数Re_λ	100.3

表2 AMR算法参数及网格设置

Table 2 AMR algorithm parameters and mesh settings

参数	DNS	LES（0.4）	LES（0.6）
最高网格分辨率N_max	512	256	256
最低网格分辨率N_min	64	32	32
小波离散误差δ	0.1	0.4	0.6

图1 能谱图及耗散谱图

Fig.1 Energy spectrum and dissipation spectrum

图2 颗粒对径向分布函数及其相对误差

Fig.2 Particle pair radial distribution function and its relative errors

图3 颗粒对径向相对速度及其相对误差

Fig.3 Particle pair radial relative velocity and its relative errors

图4 颗粒对碰撞核函数及其相对误差

Fig.4 Particle pair collision kernel and its relative errors

图5 不同网格密度颗粒对径向分布函数随颗粒对间距比的变化曲线

Fig.5 Distribution of the particle pair g(r) as a function of r/0.25η for different mesh densities

图6 不同网格密度颗粒对径向相对速度随颗粒对间距比的变化曲线

Fig.6 Distribution of the particle pair <|wr |> as a function of r/0.25η for different mesh densities

图7 不同网格密度颗粒对碰撞核函数随颗粒对间距比的变化曲线

Fig.7 Distribution of the particle pair βcol as a function of r/0.25η for different mesh densities

图8 惯性颗粒分布及其背景流场涡量图

Fig.8 Instantaneous snapshot of inertial particle distribution and background flow field vorticity

参考文献 30

[1]	Xiong Y, Cai R W, Li J, et al. Structural subgrid scale model based on wavelet filter for large eddy simulation of particle-laden turbulence[J]. International Journal of Multiphase Flow, 2021, 143: 103730.
[2]	Zhou Z D, Wang S Z, Yang X L, et al. A structural subgrid-scale model for the collision-related statistics of inertial particles in large-eddy simulations of isotropic turbulent flows[J]. Physics of Fluids, 2020, 32(9): 095103.
[3]	Smagorinsky J. General circulation experiments with the primitive equations[J]. Monthly Weather Review, 1963, 91(3): 99-164.
[4]	Wang Q Z, Squires K D. Large eddy simulation of particle-laden turbulent channel flow[J]. Physics of Fluids, 1996, 8(5): 1207-1223.
[5]	Fukagata K, Zahrai S, Bark F H. Dynamics of Brownian particles in a turbulent channel flow[J]. Heat and Mass Transfer, 2004, 40(9): 715-726.
[6]	Amiri A E, Hannani S K, Mashayek F. Large-eddy simulation of heavy-particle transport in turbulent channel flow[J]. Numerical Heat Transfer, Part B: Fundamentals, 2006, 50(4): 285-313.
[7]	Simonin O, Zaichik L I, Alipchenkov V M, et al. Connection between two statistical approaches for the modelling of particle velocity and concentration distributions in turbulent flow: the mesoscopic Eulerian formalism and the two-point probability density function method[J]. Physics of Fluids, 2006, 18(12): 125107-125109.
[8]	Fede P, Simonin O, Villedieu P, et al. Stochastic modeling of the turbulent subgrid fluid velocity along inertial particle trajectories[C]//Proceedings of the Summer Program. Center for Turbulence Research, 2006: 247-258.
[9]	Vinkovic I, Aguirre C, Ayrault M, et al. Large-eddy simulation of the dispersion of solid particles in a turbulent boundary layer[J]. Boundary-Layer Meteorology, 2006, 121(2): 283-311.
[10]	Pozorski J, Apte S V. Filtered particle tracking in isotropic turbulence and stochastic modeling of subgrid-scale dispersion[J]. International Journal of Multiphase Flow, 2009, 35(2): 118-128.
[11]	Jin G D, He G W, Wang L P, et al. Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence[J]. International Journal of Multiphase Flow, 2010, 36(5): 432-437.
[12]	Jin G D, He G W. A nonlinear model for the subgrid timescale experienced by heavy particles in large eddy simulation of isotropic turbulence with a stochastic differential equation[J]. New Journal of Physics, 2013, 15(3): 035011.
[13]	Shotorban B, Mashayek F. A stochastic model for particle motion in large-eddy simulation[J]. Journal of Turbulence, 2006, 7: N18.
[14]	Pesmazoglou I, Kempf A M, Navarro-Martinez S. A dynamic model for the Lagrangian stochastic dispersion coefficient[J]. Physics of Fluids, 2013, 25(12): 125108.
[15]	Innocenti A, Marchioli C, Chibbaro S. Lagrangian filtered density function for LES-based stochastic modelling of turbulent particle-laden flows[J]. Physics of Fluids, 2016, 28(11): 115106.
[16]	Fung J C H, Hunt J C R, Malik N A, et al. Kinematic simulation of homogeneous turbulence by unsteady random Fourier modes[J]. Journal of Fluid Mechanics, 1992, 236: 281-318.
[17]	Ray B, Collins L R. A subgrid model for clustering of high-inertia particles in large-eddy simulations of turbulence[J]. Journal of Turbulence, 2014, 15(6): 366-385.
[18]	Kuerten J G M. Subgrid modeling in particle-laden channel flow[J]. Physics of Fluids, 2006, 18(2): 025108.
[19]	Shotorban B B. Modeling of Subgrid-Scale Effects on Particles in Large-Eddy Simulation of Turbulent Two-Phase Flows[M]. Chicago: University of Illinois at Chicago, 2005.
[20]	Shotorban B, Zhang K K Q, Mashayek F. Improvement of particle concentration prediction in large-eddy simulation by defiltering[J]. International Journal of Heat and Mass Transfer, 2007, 50(19/20): 3728-3739.
[21]	Park G I, Bassenne M, Urzay J, et al. A simple dynamic subgrid-scale model for LES of particle-laden turbulence[J]. Physical Review Fluids, 2017, 2(4): 044301.
[22]	Bassenne M, Esmaily M, Livescu D, et al. A dynamic spectrally enriched subgrid-scale model for preferential concentration in particle-laden turbulence[J]. International Journal of Multiphase Flow, 2019, 116: 270-280.
[23]	Vreman A W. An eddy-viscosity subgrid-scale model for turbulent shear flow: algebraic theory and applications[J]. Physics of Fluids, 2004, 16(10): 3670-3681.
[24]	Maxey M R, Riley J J. Equation of motion for a small rigid sphere in a nonuniform flow[J]. Physics of Fluids, 1983, 26(4): 883-889.
[25]	Auton T R, Hunt J C R, Prud'Homme M. The force exerted on a body in inviscid unsteady non-uniform rotational flow[J]. Journal of Fluid Mechanics, 1988, 197: 241-257.
[26]	Boivin M, Simonin O, Squires K D. Direct numerical simulation of turbulence modulation by particles in isotropic turbulence[J]. Journal of Fluid Mechanics, 1998, 375(1): 235-263.
[27]	罗坤. 气固两相自由剪切流动的直接数值模拟和实验研究[D]. 杭州: 浙江大学, 2005.
	Luo K. Direct numerical simulation and experimental study on gas-solid two-phase free shear flow[D]. Hangzhou: Zhejiang University, 2005.
[28]	Sweldens W. The lifting scheme: a construction of second generation wavelets[J]. SIAM Journal on Mathematical Analysis, 1998, 29(2): 511-546.
[29]	van Hooft J A, Popinet S, van Heerwaarden C C, et al. Towards adaptive grids for atmospheric boundary-layer simulations[J]. Boundary-Layer Meteorology, 2018, 167(3): 421-443.
[30]	Pope S B. Turbulent Flows[M]. Cambridge: Cambridge University Press, 2000.