错位六弯叶桨搅拌假塑性流体流场宏观不稳定性数值模拟
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
错位六弯叶桨搅拌假塑性流体流场宏观不稳定性数值模拟栾德玉;魏星;陈一鸣
【摘要】采用分离涡模型,对错位六弯叶(6PBT)搅拌槽内流场结构和宏观不稳定性(MI)进行了数值模拟.工作介质分别选用去离子水和不同质量分数的黄原胶水溶液,并将水的速度场分布与PIV实验结果进行了比较.通过采集监测点的速度时间序列,结合MATLAB软件编程,计算得到流场宏观不稳定频率的变化.结果表明:流场结构和速度矢量的计算值与PIV实验数据吻合较好,分离涡模型的计算结果可靠;提高转速,槽内假塑性流体MI频率峰值增大,脉动强度提高,当转速达到225 r·min-1时,MI频率特征消失,频谱图呈现谱带现象,意味着流场进入混沌;流体的流变性对MI没有影响,MI现象是流体流动的共有特征.
【期刊名称】《化工学报》
【年(卷),期】2018(069)005
【总页数】8页(P1999-2006)
【关键词】错位六弯叶搅拌器;假塑性流体;宏观不稳定性;分离涡模型
【作者】栾德玉;魏星;陈一鸣
【作者单位】青岛科技大学机电工程学院,山东青岛266061;青岛科技大学机电工程学院,山东青岛266061;青岛科技大学机电工程学院,山东青岛266061
【正文语种】中文
【中图分类】TQ027.2
引言
搅拌槽内流体的流动过程存在一个大尺度、低频率的流型变化,称之为宏观不稳定现象(macro instability, MI),在近桨区的速度频谱图中表现为一个低频显著峰值。
宏观不稳定是槽内流体运动能量耗散与流场大涡运移的综合表现[1],它的存在有
助于增强流体的湍动程度,提高混合效率,对传质、传热性能和局部气含量也有影响[2-3]。
因此,研究这一现象对于深入了解流动形态和流体混合机理具有重要意义。
近年来,国内外学者采用不同的实验方法,相继开展了对搅拌槽宏观不稳定现象的研究[4-8],取得的成果表明:宏观不稳定频率与搅拌槽的尺寸、搅拌桨型式
及流体物性参数等因素无关,在不同Reynolds数范围内,宏观不稳定频率有一定的差异。
假塑性流体在工业生产中发挥着十分重要的作用[9-10],具有剪切变稀化的流变特性,会在搅拌槽内形成不同于牛顿流体的独特流场结构,从而改变混合效率[11-12],因此假塑性流体的混合一直是搅拌技术的难点之一[13-15]。
在搅拌槽内提高流体流动的混沌程度,是有效提高假塑性流体混合效率的根本途径[16-19],错位
六弯叶搅拌器(6PBT)可产生不对称流场结构,达到提高流场混沌程度的目的[20-22]。
为此本文采用分离涡数值模拟(DES)方法,结合MATLAB软件编程,研究
6PBT桨流场结构和宏观不稳定性特征,探索假塑性流体搅拌流场诱发混沌混合的有效方法。
1 几何模型与网格划分
平底圆柱形搅拌槽内径 T = 210 mm,内壁均布四块挡板,挡板宽度 w = T/10,液位高度 H = T;搅拌轴直径 d = 16 mm,搅拌桨直径 D = T/2,桨叶宽度 b = T/10,厚度 2 mm,后弯角为30°;桨叶中心离底高度 C=T/3。
坐标原点位于桨
叶中心。
P1点是速度采样点,距轴心和槽底距离分别为 57 mm和 70 mm,如图
1所示。
图1 搅拌槽及搅拌器结构示意图Fig.1 Schematic of mixing tank and agitator 利用 Gambit ( Fluent Inc. ) 软件,对整个计算域采用四面体网格(tetrahedral)离散,在桨叶和挡板区域进行局部网格加密,以桨叶附近的速度大小及搅拌器功率准数的改变量均不超过 3%为条件,来进行网格无关性验证,最终网格单元数确定为1309900个左右,其中最大和最小网格尺寸分别为3 mm和0.5 mm。
2 计算方法
2.1 数值模型
采用质量分数分别为0.5%、0.75%、1.0% 的黄原胶水溶液,这是典型的假塑性
流体,使用 DV3T流变仪,在恒定温度30℃下测得流变学参数,如表1所示。
表1 不同浓度黄原胶水溶液的流变学参数Table 1 Rheological parameters of different concentration of xanthan gum solutionsConcentration wt/%(mass) n K/(Pa· sn) ρ/(kg·m−3)0.5 0.63 0.32 9980.75 0.41 1.9 9951.0 0.3 5.1 992
黄原胶水溶液处在湍流状态时,其黏度模型符合power law模型,表观黏度η 的表达式为
根据 Metzner-Otto 关系式[23],平均剪切速率与搅拌转速的关系为
剪切速率常数的表达式为[24]
则表观黏度η 可表示为
湍流场模拟采用 DES 模型,该模型能够很好地捕捉流场的湍流特征和宏观不稳定
性[25-26]。
速度场的计算采用动网格法(dynamic mesh),将整个搅拌流体域作
为一体,整个区域的网格运动,该方法非常适合处理非稳态的流动问题,其网格的迭代速率大约是滑移网格法迭代速率的三倍[27]。
2.2 模拟方法
动量方程按二阶上迎风格式离散,压力速度耦合项采用SIMPLE算法,为了保证
模拟精度,将各变量的收敛残差设为10−5。
时间步长的最佳大小为搅拌器旋转周期的 1/120[28],为此模拟计算的时间步长设定为0.001 s。
黄原胶水溶液的黏度
模型为幂律(power law)模型,近壁区域的流动计算采用标准壁面函数法,槽内壁面定义为无滑移边界条件(no slip),液面处采用对称边界条件(symmetric)。
在
建立的黏度模型中,n设置为1.0,K设置为0.001 Pa·sn,即可通过假塑性流体
模型来计算牛顿流体(水)的流动[23]。
3 计算模型验证
3.1 PIV 实验
本实验所用的PIV测试系统主要包括激光发射器、跨帧CCD相机、同步器、计算机以及图像数据处理系统,激光发射器与同步器均为 Dyntec Dynamic A/S公司
生产,CCD相机为尼康公司生产,像素为2048×2048。
搅拌装置主要由搅拌桨、搅拌槽、方形槽、驱动电机,扭矩传感器等部件组成。
如图2所示。
图2 PIV实验设备布置Fig.2 PIV experimental equipment arrangement
PIV 实验选取示踪粒子密度为ρP = 998 kg⋅ m−3的聚苯乙烯微球(polystyrene),粒子直径 dP =1~5 μm,具有良好的成像可见性和很好的跟随性。
片光厚度调整至 1 mm 左右,激光脉冲间隔为220 μs,所有工况的速度采样频率设为 10 Hz。
利用粒子图像分析系统软件MicroVec V2.3,对采集的速度矢量图像进行分析,
用于验证模拟结果的可靠性。
实验所用搅拌槽和错位六弯叶桨的结构及尺寸大小与数值模型完全一致,介质选用去离子水,在转速N =60 r⋅min−1条件下进行,此
时 Re = 11005,流动处于湍流状态。
3.2 宏观流场结构
从图3可以看出,二者流场结构基本吻合,说明所建立的分离涡模型能够较好地反映搅拌流场特征。
桨叶的上下端都存在涡心,但涡形明显不一致,表明错位桨可产生不对称流场结构。
图3 纵截面宏观流场结构的速度矢量图Fig.3 Distribution of velocity field in PIV experiment
3.3 时均速度分布
选取搅拌槽内不同高度截面的速度矢量分布,将 PIV 实验与分离涡模拟结果进行对比分析。
图4~图6分别显示了纵向截面内速度大小、轴向和径向速度的分布。
可以看出,模拟计算值与 PIV 实验数据吻合较好,曲线的变化趋势是一致的,说明利用分离涡模型建立的模拟算法计算流场的速度分布是可行的,该方法能够捕捉到流场的流动特征。
图4 不同高度速度大小沿径向的变化Fig.4 Velocity magnitude vary from radial direction at different high level
图5 不同高度轴向速度沿径向的变化Fig.5 Axial velocity vary from radial direction at different high level
图6 不同高度径向速度沿径向的变化Fig.6 Radial velocity vary from radial direction at different high level
4 结果与讨论
4.1 流场宏观结构
以质量分数wt = 0.5% 的黄原胶水溶液为例,在不同转速下的宏观流场结构分别如图7所示。
可以看出,不同转速下流场的涡心A与C,B与D都不对称,左右涡形有很大差异,这说明流场的混合隔离区不复存在,特别是转速为225
r·min−1时,流场出现了从主涡分离出的两个小涡,这将更有效地将能量耗散出去,使得混合更好。
从流场分布来看,表明基于分离涡模型的模拟方法可以精确地捕捉搅拌槽内非稳态的流动特征。
4.2 速度时间序列
对质量分数为 0.5% 的黄原胶水溶液径向速度取值,采样时间为200 s,采样频率为10 Hz,速度时间序列如图8所示。
不同转速下,槽内流体的径向速度总是围绕一个值在上下波动,表明槽内流体的流动处于非稳定状态。
随着转速的增加,流体的速度波动频率不断加大,表现出湍流脉动程度增强的特征。
4.3 宏观不稳定性分析
图7 不同转速下宏观流场结构Fig.7 Macroscopic flow field at different speeds
4.3.1 转速对宏观不稳定性的影响不同转速下搅拌槽内MI频率如图9所示。
可以看出,在较低转速N = 60 r·min−1下,脉动强度较低,MI只有单一频率峰值f = 0.198 Hz,对应于无量纲频率f* =f/N =0.198。
当转速为125 r·min−1时,MI频率在整个频段内的幅值都随之增大,此时槽内流体的脉动程度明显增加,出现两个MI频率特征峰值,分别为f1=0.291 Hz和f2= 0.436 Hz,对应于无量纲频率分别为f*(1)=0.140和f*(2)=0.209,这表明槽内流体处于较为明显的不稳定状态,这种不稳定状态是由于流场中大涡的拟序结构运移造成的,通过控制流场的拟序结构可以改变MI频率的分布。
当转速提高到180 r·min−1时,MI频率特征峰值为 f = 0.686 Hz,对应于无量纲频率为f* = 0.229,流场的宏观不稳定特征得到进一步增强。
在转速增加到225 r·min−1时,流场的MI频率特征消失,呈现出多尺度结构特征,表明流场的拟序结构发生了明显变化,这必然导致流场进入混沌状态。
比较不同工况下流场的MI频率变化可知,提高搅拌转速,整个频段内的频率幅值
随着增加,脉动强度都相应提高,峰值也随之增大,这说明槽内流体的能量不断增加,湍动程度增强。
转速继续增加,MI频率特征消失,呈现多尺度谱带现象,此
时槽内流体脉动随机性加强,意味着流场进入混沌,这与樊建华等[29]和 Galletti 等[30]研究涡轮搅拌器流场得到的结论是一致的。
4.3.2 流变性对宏观不稳定性的影响图10和图11为转速N =125 r·min−1时,不同流变指数流体的速度时间序列和MI频谱图,可以看出,径向速度的波动情况大致相同,不同流变指数下的MI频率峰值基本相同,分布特征也很相近,流变性对流场的MI基本没有影响,这表明MI是流体流动所具有的共性。
5 结论
本文采用分离涡模型,对错位六弯叶搅拌槽内黄原胶水溶液的流场特性进行研究,分析了搅拌转速和流体流变性对MI的影响,得到如下结论。
图8 不同转速下径向速度时间序列Fig.8 Radial velocity time series at different speeds
(1)随转速的增加,假塑性流体的流场中会同时呈现多种涡形,这种大涡运移和小涡扩散的共同作用,将加速能量耗散,同时表明分离涡模型可以精确捕捉搅拌槽内假塑性流体非稳态的流动特征。
图9 不同转速下的MI频谱图Fig.9 Macroscopic instability spectrum at different speeds
(2)提高搅拌转速,可使槽内假塑性流体MI频率峰值增大,脉动强度相应提高。
当转速达到225 r·min−1时,流体的MI频率特征消失,出现了谱带现象,表明
流场拟序结构发生了明显变化,将呈现多尺度结构特征,这必然导致流场进入混沌。
图10 不同流变性的速度时间序列Fig.10 Velocity time series of fluid at different rheology
图11 不同流变性的MI频谱图Fig.11 Macroscopic instability spectrum of
fluid at different rheology
(3)不同流变性流体的MI频率特征峰值基本相同,流变性对流场的MI没有影响,表明MI现象是流体流动的共有特征。
符号说明
B——形状系数
b——桨叶宽度,mm
C——桨叶中心离底高度,mm
D——搅拌桨直径,mm
d——搅拌轴直径,mm
dP——粒子直径,μm
f——MI单一频率峰值,Hz
f*——无量纲频率
f1,f2——分别为两个MI频率特征峰值,Hz
H——液体高度,mm
K——稠度系数,Pa·sn
KS——Metzner-Otto 常数,对于径向流搅拌器,KS=11.5
N ——搅拌桨转速,r·min−1
n——流变指数
Re——Reynolds数
T——搅拌槽内径,mm
t——时间,s
W——挡板宽度,mm
wt——质量浓度,%
γ——平均剪切速率,s−1
η——表观黏度,Pa·s
ρ——密度,kg·m−3
ρP——示踪粒子的密度,kg·m−3
References
【相关文献】
[1] HARTMANN H, DERKSEN J, VAN DEN AKKER H E A, et al.Macroinstability uncovered in
a Rushton turbine stirred tank by means of LES[J]. American Institute of Chem. Eng., 2004, 50(10): 2383-2393.
[2] HASAL P, MONTES J L, BOISSON H C, et al. Macro-instabilities of velocity field in stirred vessel: detection and analysis[J]. Chem. Eng.Sci., 2000, 55(2): 391-401.
[3] DUCCI A, DOULGERAKIS Z, YIANNESKIS M. Decomposition of flow structures in stirred reactors and implications for mixing enhancement [J]. Ind. Eng. Chem. Res., 2008, 47(10): 3664-3676.
[4] PAGLIANTI A, LIU Z, MONTANTE G, et al. Effect of macroinstabilities in single- and multiple-impeller stirred tanks[J]. Ind.Eng .Chem. Res., 2008, 47(14): 4944-4952.
[5] FAN J, RAO Q, WANG Y, et al. Spatio-temporal stirred vessel via digital particle image velocimetry macro-instability analysis of DPIV[J]. Chem. Eng .Sci., 2004, 59(6): 1863-1873.
[6] 樊建华, 饶麒, 王运东. 搅拌槽内流场脉动的频谱分析研究[J]. 高校化学工程学报,2004, 18(3): 287-292.FAN J H, RAO Q, WANG Y D, et al. Spectral analysis of the velocity fluctuations in
a mechanically stirred tank[J]. J. Chem.Eng. Chin. Univ.,2004, 18(3): 287-292.
[7] GALLETTI C, PAGLIANTI A, YIANNESKIS M, et al. Observations on the significance of instabilities turbulence and intermittent motions on fluid mixing processes in stirred reactors[J]. Chem. Eng .Sci., 2005,60(4): 2317-2331.
[8] GALLETTI C, LEE K C, PAGLIANTI A, et al. Reynolds number and impeller diameter effects on instabilities in stirred vessels[J]. Chem.Eng. Sci., 2004, 50(3): 2050-2063.
[9] EIN-MOZAFFARI F, UPRETI S R. Using ultrasonic Doppler velocimetry and CFD modeling to investigate the mixing of non-Newtonian fluids possessing yield stress[J]. Chem. Eng. Res. Des.,2009, 87(4): 515-523.
[10] ADAMS L W, BARIGOU M. CFD analysis of caverns and pseudocaverns developed during mixing of non-Newtonian fluids[J]. Chem.Eng. Res. Des., 2007, 85(5): 598-604.
[11] 栾德玉, 周慎杰, 陈诵英. 错位六弯叶桨搅拌假塑性流体的洞穴变化[J].机械工程学报, 2012,
48(16): 152-157.LUAN D Y, ZHOU S J, CHEN S Y, et al. Cavern development of pseudoplastic fluids stirred by impeller of perturbed six-bent-bladed turbine[J]. Journal of Mechanical Engineering, 2012, 48(16): 152-157.
[12] XIAO Q, YANG N, ZHU J, et al. Modeling of cavern formation in yield stress fluids in stirred tanks[J]. AIChE J., 2014, 60(8): 3057-3070.
[13] PAKZAD L, EIN-MOZAFFARI F, CHAN P. Using electrical resistance tomography and computational fluid dynamics modeling to study the formation of cavern in the mixing of pseudoplastic fluids possessing yield stress[J]. Chem. Eng. Sci., 2008, 63(9): 2508-2522. [14] GALINDOT E, NIENOW A W. Mixing of highly viscous simulated xanthan fermentation broths with the lightnin A-315 impeller[J].Biotechnology Progress, 2010, 8(3): 233-239. [15] SOSSA-ECHEVERRIA J, TAGHIPOUR F. Computational simulation of mixing flow of shear thinning non-Newtonian fluids with various impellers in a stirred tank[J]. Chem. Eng. Process., 2015, 93: 66-78.
[16] XU B, GILCHRIST J F. Shear migration and chaotic mixing of particle suspensions in a time-periodic lid-driven cavity[J]. Physics of Fluids,2010, 22(5): 1-7.
[17] XU B, GILCHRIST J F. Shear migration and chaotic mixing of particle suspensions in a time-periodic lid-driven cavity[J]. Physics of Fluids,2010, 22(5): 1-7.
[18] CUDAK M, KARCZ J. Local momentum transfer process in a wall region of an agitated vessel equipped with an eccentric impeller[J]. Ind.Eng. Chem. Res., 2011, 50(7): 4140-4149.
[19] 栾德玉, 周慎杰, 陈诵英. 错位六弯叶搅拌槽内假塑性流体的混合特性[J]. 高校化学工程学报, 2012, 26(5): 787-792.LUAN D Y, ZHOU S J, CHEN S Y, et a1. Mixing characteristics of pseudoplastic fluid in a stirred tank with the stirrer composed of perturbed six-bent-bladed turbine[J]. Journal of Chemical Engineering of Chinese Universities,2012, 26(5): 787-792.
[20] 栾德玉, 周慎杰, 陈诵英, 等. 错位六弯叶在假塑性流体中的混沌搅拌特性[J]. 化学工程, 2011, 39(9): 41-46.LUAN D Y, ZHOU S J, CHEN S Y, et al. Chaotic agitation characteristics of shifted 6-bent-blade impeller in pseudoplastic fluid[J]. Chemical Engineering(China),2011, 39(9): 41-46.
[21] 栾德玉, 张盛峰, 郑深晓, 等. 基于流固耦合的错位桨搅拌假塑性流体动力学特性[J]. 化工学报, 2017, 68(6): 2328-2335.LUAN D Y, ZHANG S F, ZHENG S X, et al. Dynamic characteristics of impeller of perturbed six-bent-bladed turbine in pseudoplastic fluid based on fluid-structure interaction[J]. CIESC Journal,2017, 68(6):2328-2335.
[22] LUAN D Y, CHEN Q, ZHOU S J. Numerical simulation and analysis of power consumption and Metzner-Otto constant for impeller of 6PBT[J]. Chin. J. Mech. Eng., 2014, 27(3): 635-640.
[23] METZNER A B, OTTO R E. Agitation of non-Newtonian fluids[J].AIChE J., 1957, 3(1): 3-
11.
[24] HOUARI A. Effect of some parameters on the performance of anchor impellers for stirring shear-thinning fluids in a cylindrical vessel[J].Journal of Hydrodynamics, 2016,
28(4): 669-675.
[25] JAYARAJU, BROUNS S T. Large eddy and detached eddy simulations of fluid flow and particle deposition in a human mouth-throat[J].Journal of Aerosol Science, 2008, 39(10): 862-875.
[26] 杨锋苓. 偏心搅拌槽内宏观不稳定性的分离涡模拟及实验研究[D].济南: 山东大学, 2010.YANG
F L. Detached eddy simulation and PIV study of the macroinstability in eccentrically stirred tanks[D]. Jinan: Shandong University,2010.
[27] 黄思, 杨富翔, 郭京, 等. 运用三维动网格技术模拟计算离心泵非定常流动[J]. 科技导报, 2013, 31(24): 33-36.HUANG S, YANG F X, GUO J, et al. Numerical simulation of 3D unsteady flow in centrifugal pump by dynamic mesh technique[J].Science&Technology Review, 2013,
31(24): 33-36.
[28] LI M Z, GRAEME W, DEREK W, et a1. Scale up study of retreat curve impeller stirred tanks using LDA measurements and CFD simulation[J].Chemical Engineering Journa1, 2005, 108(1): 81-90.
[29] 樊建华, 饶麒, 王运东. 搅拌槽内流场脉动的频谱分析研究[J]. 高校化学工程学报, 2004, 18(3): 287-292.FAN J H, RAO Q, WANG Y D. Spectral analysis of the velocity fluctuations in a mechanically stirred tank[J]. Journal of Chemical Engineering of Chinese Universities, 2004, 18(3): 287-292.
[30] GALLETTI C, PAGLIANTI A, LEE K C, et al. Reynolds number and impeller diameter effects on instabilities in stirred vessels[J]. AIChE J.,2004, 50(9): 2050-2063.。