基于双向流固耦合的错位六弯叶桨搅拌特性
发布时间:2020-07-11 02:12
【摘要】:基于ANSYS Workbench分析平台,采用双向流固耦合计算方法,对六弯叶桨(6BT)和错位六弯叶桨(6PBT)搅拌性能进行了对比分析,根据桨叶与流体之间相互耦合运动特性,考察了桨叶的等效应力和变形量的分布,探讨了流场结构对混合效率的影响。结果表明:速度矢量的计算值与PIV实验值吻合较好,基于双向流固耦合计算模型的数据结果可靠;同6BT桨相比,6PBT桨叶端部应力提高了51%,根部应力降低了22.2%,应力分布更趋均匀化,可有效提高桨叶强度,增强能量传递,而总变形量与6BT桨基本相同;6PBT桨能够产生不对称流场结构,混合时间明显降低,混合效率显著提高,其中6PBT桨的单位体积混合能只有6BT桨48%,体现出错位桨的优越性,这可为搅拌器的优化设计提供理论参考。
【图文】:
的叶端速度Utip=0.7m/s,雷诺数Re=23339,此时流动为湍流状态。搅拌轴与叶轮均采用结构钢,密度为7.85×103kg/m3,泊松比为0.3,弹性模量为2.0×1011Pa,屈服强度为2.5×108Pa。使用Workbench中ICEM对计算区域进行网格划分,采用非结构化四面体单元进行离散,对于速度梯度较大的转子区域,将网格进行加密处理,如图2所示。为验证网格无关性,以转子区的速度大小及搅拌器功率准数的改变量均不超过3%为依据,本文最终采用1106500个左右的网格单元数对流体计算域进行离散。图2搅拌槽计算网格示意2计算策略双向流固耦合计算是一个复杂的非定常过程,在流体区域,选用workbench中的FluidFlow(Fluent)模块对流体运动进行计算,采用多重参考系法(MRF)和κ-ε双方程湍流模型,叶片与流体的接触面设置为system-coupling,槽内壁面设置为deforming,所有变量的收敛残差均小于1×10-4;在固体区域,通过workbench中的Transientstructure模块对搅拌器进行瞬态动力学分析,搅拌轴在径向和轴向上固定约束,叶片面设置为fluidsolidinterface,两个区域设置相同的时间步长0.001s。最后将两个区域的模型一并导入workbench中的system-coupling模块中,将流体区域中的system-coupling面与固体区域中的fluidsolidi
的。根据已有的研究结果【14】,将监测点选在靠近液面处比较合适。为此,本文选择监测点位置坐标为(95,0,115),在计算浓度场分布时,组分收敛残差设定为10-7,时间步长设置为1×10-3s。图6六弯叶桨纵截面速度示意图7错位六弯叶桨纵截面速度示意搅拌器的混合速率时常用混合时间数Tm来表征,其表达式为:Tm=N·t95(1)式中:Tm———常用混合时间数;N———转速,r/s;t95———混合时间,s。Tm数值越小,表明搅拌器的混合速率越高,当搅拌器转速相同时,可用混合时间表征。图8为搅拌槽内示踪剂的浓度对时间的响应曲线。由图8可知,对于六弯叶桨,t95=2.7s,而错位桨的混合时间t95=1.4s,明显小于六弯叶桨。可见错位桨的混合速率更高。图8示踪剂浓度变化曲线利用单位体积混合能WV可表征混合效率,WV越小,混合效率越高。它的表达式为:WV=PV·t95(2)PV=2πNM/V(3)·48·石油化工设备技术2017年
本文编号:2749811
【图文】:
的叶端速度Utip=0.7m/s,雷诺数Re=23339,此时流动为湍流状态。搅拌轴与叶轮均采用结构钢,密度为7.85×103kg/m3,泊松比为0.3,弹性模量为2.0×1011Pa,屈服强度为2.5×108Pa。使用Workbench中ICEM对计算区域进行网格划分,采用非结构化四面体单元进行离散,对于速度梯度较大的转子区域,将网格进行加密处理,如图2所示。为验证网格无关性,以转子区的速度大小及搅拌器功率准数的改变量均不超过3%为依据,本文最终采用1106500个左右的网格单元数对流体计算域进行离散。图2搅拌槽计算网格示意2计算策略双向流固耦合计算是一个复杂的非定常过程,在流体区域,选用workbench中的FluidFlow(Fluent)模块对流体运动进行计算,采用多重参考系法(MRF)和κ-ε双方程湍流模型,叶片与流体的接触面设置为system-coupling,槽内壁面设置为deforming,所有变量的收敛残差均小于1×10-4;在固体区域,通过workbench中的Transientstructure模块对搅拌器进行瞬态动力学分析,搅拌轴在径向和轴向上固定约束,叶片面设置为fluidsolidinterface,两个区域设置相同的时间步长0.001s。最后将两个区域的模型一并导入workbench中的system-coupling模块中,将流体区域中的system-coupling面与固体区域中的fluidsolidi
的。根据已有的研究结果【14】,将监测点选在靠近液面处比较合适。为此,本文选择监测点位置坐标为(95,0,115),在计算浓度场分布时,组分收敛残差设定为10-7,时间步长设置为1×10-3s。图6六弯叶桨纵截面速度示意图7错位六弯叶桨纵截面速度示意搅拌器的混合速率时常用混合时间数Tm来表征,其表达式为:Tm=N·t95(1)式中:Tm———常用混合时间数;N———转速,r/s;t95———混合时间,s。Tm数值越小,表明搅拌器的混合速率越高,当搅拌器转速相同时,可用混合时间表征。图8为搅拌槽内示踪剂的浓度对时间的响应曲线。由图8可知,对于六弯叶桨,t95=2.7s,而错位桨的混合时间t95=1.4s,明显小于六弯叶桨。可见错位桨的混合速率更高。图8示踪剂浓度变化曲线利用单位体积混合能WV可表征混合效率,WV越小,混合效率越高。它的表达式为:WV=PV·t95(2)PV=2πNM/V(3)·48·石油化工设备技术2017年
【相似文献】
相关期刊论文 前7条
1 胡跃华;蒋诚航;闫怀磊;万先平;金志江;;流固耦合作用下固支输液管道有限元分析[J];化工机械;2012年02期
2 张杰;赵军;冯武文;;某大型塔器流固耦合的数值仿真[J];科技信息;2009年14期
3 潘海丽;张亚新;;管道流体双向流固耦合的动力学模拟分析[J];中国石油和化工标准与质量;2013年06期
4 周利剑;;液体与储液结构互相作用流固耦合的实现[J];油气田地面工程;2009年11期
5 王艳林;吴兰鹰;王自东;;新型输液90°弯头结构流固耦合特性[J];北京科技大学学报;2010年07期
6 刘延安;;流固耦合管道固有频率及参数灵敏度分析[J];机械科学与技术;2012年03期
7 ;[J];;年期
相关博士学位论文 前1条
1 石代嗯;搅拌槽内流固耦合的实验研究与数值模拟[D];北京化工大学;2014年
相关硕士学位论文 前1条
1 王恒宇;高塔流固耦合振动特性分析及结构改进研究[D];湘潭大学;2015年
本文编号:2749811
本文链接:https://www.wllwen.com/kejilunwen/lxlw/2749811.html