基于脉动压力时均值的除尘器气液耦合研究
基于湍流脉动压力的波数—频率谱预报流噪声
基于湍流脉动压力的波数—频率谱预报流噪声伍宏亮;周其斗;吕晓军;孟庆昌【摘要】[目的]根据Lighthill声类比方程及其发展理论,可以将壁面湍流脉动压力的波数—频率谱作为声源项来预报流噪声,且分析湍流脉动压力的波数—频率谱有助于了解湍流结构的时空关联特性.[方法]以NACA 0012翼型为例,采用大涡模拟(LES)方法进行流场仿真计算,然后通过Fourier变换得到壁面湍流脉动压力波数—频率谱的数值解,并与Corcos的平板湍流边界层脉动压力波数—频率谱模型进行比较;在此基础上,将该波数—频率谱作为声源输入,代入Goldstein版本的声类比方程中预报辐射噪声,并与软件计算的流噪声结果以及Brooks试验拟合结果进行比较.[结果]结果发现:小曲率变化的NACA 0012翼型表面的波数—频率谱具有与平板表面相似的一般特性;在中、低频段采用该方法预报的流噪声结果与Brooks试验结果拟合更好.[结论]所得结果表明开展波数—频率谱研究是有必要的,将其作为主要声源项来预报亚声速下产生的流噪声是合理的.%[Objectives]According to the Lighthill acoustic analogy equation and its development theory, it is feasible to analyze the wavenumber-frequency spectrum of turbulent wall pressure fluctuations,then make it an acoustic source in order to predict flow noise. Moreover, the study of the wavenumber-frequency spectrum is useful for understanding the temporal and spatial characteristics of turbulent structures.[Methods]Taking the NACA 0012 airfoil,which was studied by Brooks,as an example,we employ the Large Eddy Simulation (LES)method to calculate the flow field and obtain a numerical solution of the wavenumber-frequency spectrum via the Fourier transform. On this basis,we take the wavenumber-frequency spectrum as an input conditionfor predicting the radiated noise using the acoustic analogue equation of the Goldstein version. At the same time,acoustic software is used to calculate the flow noise. Comparing these two sets of results with Brooks' empirical formula,the sound pressure level is found to be within the same order of magnitude.[Results]The results show that the spectrum on an airfoil surface with a small curvature change is comparable with the Corcos spectrum model on a flat plate,and their general characteristics are similar. Finally,we conclude that the forecast results of the method in this paper accord better with Brooks' experimental results at low and medium frequencies. [Conclusions]This shows that it is necessary to carry out the study of wavenumber-frequency spectra,and it is reasonable to make it the main sound source in order to predict flow noise produced at subsonic speed.【期刊名称】《中国舰船研究》【年(卷),期】2017(012)006【总页数】7页(P36-42)【关键词】波数—频率谱;Fourier变换;流噪声;声类比方程【作者】伍宏亮;周其斗;吕晓军;孟庆昌【作者单位】海军工程大学舰船工程系,湖北武汉 430033;海军工程大学舰船工程系,湖北武汉 430033;海军工程大学舰船工程系,湖北武汉 430033;海军工程大学理学院,湖北武汉 430033【正文语种】中文【中图分类】U661.44湍流脉动压力是湍流非定常特性的重要表征,也是流体诱发结构振动、产生噪声的重要来源。
3天了解CFD——湍流、多相流的方程推导和数值解法(上)
如果为大量固定质量的流体微元他就是有限控制体,少量流体微团(但足够表现流体宏观量)的就是流动微元(fluid
element),但是体积和面积会变化。之所以守恒和不守恒是因为空间位置固定的微元或者控制体可以应用质量守恒
方程来推导,空间位置随流因为质量固定不能应用与质量守恒因此成为非守恒形式。
1.2 动量方程
+
�������������
+
������������������������ ������������
�������������������������
������������������������������������������������
−
������������
������������������������������������������������������������������������
�������������������������
������������������������������������������������
−
������������
������������������������������������������������������������������������
+
������������
∙
∇�
+
�������������������������������������
+
∇
∙
(������������������������)�
=
������������
������������������������ ������������
管道流固耦合振动及声传播的研究现状及展望
第5卷第2期船舶力学V ol.5N o.2 2001年4月Journal of Ship M echanics A pr.2001文章编号:1007-7294(2001)02-0082-09管道流固耦合振动及声传播的研究现状及展望刘忠族,孙玉东,吴有生(中国船舶科学研究中心,无锡214082)摘要:本文对管道振动及声传播领域的研究现状和水平分类作了比较全面的综述,并对今后开展潜艇管道振动、声传播及噪声控制的研究进行了展望。
关键词:管道系统;流固耦合作用;管道振动;声传播中图分类号:O353.1,O422文献标识码:A1引言管道振动、噪声及其控制技术的研究有着广泛的应用背景。
对于工程上的管道系统,其动力学及声学分析非常复杂,不同的流体与结构物理模型的组合可派生出不同的动力学问题。
管道系统涉及的问题很多,如水锤,流固耦合振动及稳定性,管道声传播等。
水锤直接影响到管道的安全性,同时也会激发管道系统振动。
当管内流体速度增加到某一临界值时,系统将发生失稳,这种失稳是有害的,它可以使管道系统的薄弱环节产生破坏;另外,颤振失稳也会辐射噪声。
管道振动及噪声研究领域有许多理论问题没有解决,例如:复杂空间管系振动传递特性分析,管道系统对激励的响应及衰减特性,液压系统的动态特性分析以及管道声传播及辐射特性等。
管道振动噪声研究具有重要的学术价值和工程意义。
2管道振动及声传播的研究现状2.1管道振动及噪声源的描述管道噪声的含义是:管壁结构的振动和管内流体中声波相互作用,并沿管壁和管内流体传播与辐射的过程。
作用于管道的激励有两种:力激励与声激励。
力激励主要是指泵作用于管壁与流体的动态力,声激励是指泵的水动力噪声与阀门等处的湍流噪声。
湍流噪声是由控制阀和流量调节器产生的,在控制阀和流量调节器处,通流面积突然缩小或扩大,由此产生湍流,在下游伴随着压力脉动。
这种脉动压力在流体中向下游传递,从而形成平面波而成为湍流噪声[1]。
对湍流噪声的机理已有许多研究。
基于流体-结构耦合振动的液压脉动滤波器试验研究
‘
=l
(8)
式 中 ,ao为柱 塞泵 平均 流量 ,m3/s。 由式 (8)可 知 ,瞬 时流 量可 以分 成 常量 和变 量 两部 分 。试验 中采用 7柱 塞泵 ,根 据 柱塞 泵 的流量
特性 ,取基 波和 一 次谐波 来逼 近 原 曲线 ,其瞬 时流 量 的傅立 叶级 数 展开表 达 式为
能器 、滤 波器 、可 调 节流 阀和各 种传 感 器等 。柱 塞 泵 的压 力脉 动是 由流量 脉动 的基 频及 其 谐波 组成 的 周 期信 号 ,它 的瞬 时流 量 Q可 以表 达成 】
sin( ̄b+Zon)sin( + )
Z0
O=职 ta】1 ∑sin[ ̄b+(i—Oe]=SR ̄tany
变 成 结构 振动 而被 消耗 ,流 体 压力 脉动 幅值 明显 减 小 ,压力 输 出 曲线更 加平 缓 。
衰
戋
变频调速 电机
、囊
l 压力表 蓄
专矩转速 I 7、\ ] 能
传感器 l
器
压1 力 传 感
压 l力 【传 I感
器
传 递 矩
H 为 厂
=
成 结 构 (结 构振 动体 )的振 动 ,消 除液 压 系统 中流 量和 压 力 的周 期性 波动 ,吸收 压力 冲击 ,降低 因液
压 振动 造 成 的机 械 振动 和 噪声 。
2 试验 系统及滤波器工作原理
2.1 试验 系 统 液 压 系统压 力 脉动 测试 系统 原理 如 图 2所 示 ,包 括变 频调 速器 电机 、柱 塞 式变量 泵 、溢流 阀、蓄
关 键 词 :泵源回路;压力脉动衰减;传递矩阵;耦合振动;试验研究
中图分 类号 :U664.84;TB53
基于CFD
第8卷㊀第6期2023年11月气体物理PHYSICSOFGASESVol.8㊀No.6Nov.2023㊀㊀DOI:10.19527/j.cnki.2096 ̄1642.1088基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析李泳德ꎬ㊀郭㊀力ꎬ㊀季㊀辰(中国航天空气动力技术研究院ꎬ北京100074)CharacterizationofTransonicAerodynamicDampingofRocketsBasedonCFD/CSDCouplingLIYong ̄deꎬ㊀GUOLiꎬ㊀JIChen(ChinaAcademyofAerospaceAerodynamicsꎬBeijing100074ꎬChina)摘㊀要:随着新型大推力火箭的发展ꎬ弯曲模态频率的不断降低ꎬ以及流动分离和跨声速飞行时产生的激波震荡等因素ꎬ其在跨声速飞行过程中更容易出现非定常振动发散ꎮ文章以某带助推的运载火箭模型为研究对象ꎬ通过数值计算获取火箭强迫振动时的气动阻尼ꎬ并对影响火箭气动阻尼的因素进行了分析ꎮ包括结构节点位置㊁振动振幅大小㊁脉动压力等ꎮ研究表明:助推主要起到增大气动阻尼的作用ꎻ前节点主要影响收缩段的气动阻尼ꎻ振动振幅大小和脉动压力对气动阻尼的影响可忽略不计ꎮ关键词:气动阻尼ꎻ数值计算ꎻ跨声速ꎻ气动弹性ꎻ运载火箭㊀㊀㊀收稿日期:2023 ̄09 ̄25ꎻ修回日期:2023 ̄10 ̄23第一作者简介:李泳德(1995 ̄)㊀男ꎬ工学硕士ꎬ助理工程师ꎬ主要研究方向为气动弹性分析ꎮE ̄mail:562064169@qq.com通信作者简介:季辰(1982 ̄)㊀男ꎬ工学博士ꎬ研究员ꎬ主要研究方向为气动弹性力学ꎮE ̄mail:jichen167@hotmail.com中图分类号:V475.1㊀㊀文献标志码:AAbstract:Withthedevelopmentofnewhigh ̄thrustrocketsꎬthedecreasingfrequencyofthebendingmodesoftherocketꎬaswellasthefactorssuchasflowseparationandshockoscillationsgeneratedduringtransonicflightmakeitmorepronetonon ̄constantvibration.Inthispaperꎬalaunchvehiclemodelwithboostwastakenastheresearchobjectꎬandtheaerody ̄namicdampingoftherocketduringforcedvibrationwasobtainedthroughnumericalcalculation.Thefactorsaffectingtheaerodynamicdampingoftherocketwereanalyzedꎬincludingthepositionofstructuralnodesꎬthemagnitudeofvibrationam ̄plitudeꎬpulsatingpressureandsoon.Thestudyshowsthattheboostmainlyplaystheroleofincreasingaerodynamicdamp ̄ingandthefrontnodemainlyaffectstheaerodynamicdampingofthecontractionsection.Thevibrationamplitudesizeandthepulsatingpressurehaveanegligibleeffectontheaerodynamicdamping.Keywords:aerodynamicdampingꎻnumericalcalculationꎻtransonicꎻaeroelasticityꎻlaunchvehicle引㊀言通常情况下人们认为气动力对火箭的振动起到阻尼作用ꎬ即气动阻尼为正值ꎮ然而随着大推力火箭发展ꎬ火箭的长细比逐渐加大ꎬ导致弯曲刚度越来越小ꎬ同时为了满足有效载荷的外形要求ꎬ火箭头部整流罩尺寸不断加大ꎬ后续箱体的直径却保持不变ꎬ形成了典型的锤头体外形ꎮ国内外大量的火箭研制经验表明[1 ̄9]ꎬ对于此类锤头体外形火箭的气动设计ꎬ必须要进行动态气动载荷与动态气弹稳定性分析ꎬ否则设计的疏忽可能会导致火箭结构出现毁灭性的破坏进而导致发射失败ꎮ目前常用的衡量气弹稳定性的方法是通过风洞试验来获取气动阻尼系数ꎮ早在1963年ꎬ美国国家航空航天局Ames研究中心(NASAAmesRe ̄searchCenter)采用半刚性模型开展试验研究[10]ꎬ获取火箭头部的气动阻尼来评估其稳定性ꎬ但这只能用来模拟火箭弯曲振型前节点之前部分的结构动力学特性ꎮ直到兰利研究中心(NASALangleyResearchCenter)开发了全弹性模型气动阻尼试验气体物理2023年㊀第8卷技术ꎬ其可以模拟整体的结构动力学特性以及气动外形ꎬ并应用于多款运载火箭研制[11 ̄15]ꎮ国内ꎬ中国航天空气动力技术研究院对气动阻尼问题开展过较多的研究[16 ̄20]ꎬ从模型设计方法㊁模型制作工艺㊁试验机构设计和数据处理等诸多方面ꎬ逐步改进实现了从半刚性模型到全弹性模型的过渡ꎬ并在多个型号上得到验证ꎮ然而通过风洞试验研究气动弹性问题ꎬ技术难度大ꎬ试验成本高ꎬ同时几乎不可能开展全尺寸试验ꎮ因此通过数值计算的方法开展相关研究是另一种重要的手段ꎮ刘子强等[21]实现了通过数值计算确定气动阻尼系数的技术和方法ꎬ并与试验结果进行对比ꎬ证实了该方法的可靠性ꎮ冉景洪等[22]通过模态数据结合准定常理论的方法分析了减阻杆加后体这一弹性结构的气动阻尼ꎬ结果表明减阻杆造成的分离流会对后体的气动阻尼系数产生影响ꎮ朱剑等[23]针对新一代捆绑式运载火箭发展了非结构网格下的气动阻尼计算方法ꎬ并分析了攻角㊁Mach数等参数对气动阻尼的影响ꎮ本文在之前的计算方法[23]的基础上采用IDDES模型ꎬ考虑脉动压力的影响ꎬ通过强迫振动的方式ꎬ针对捆绑式运载火箭的某一特定模态进行数值计算仿真ꎬ研究前节点位置ꎬ振动振幅ꎬ脉动压力等参数对气动阻尼的影响规律ꎮ1㊀计算方法图1为本文所用的捆绑式运载火箭的计算模型ꎬ是典型的锤头体结构ꎮ在跨声速阶段ꎬ其头部会产生激波造成激波边界层干扰ꎬ而在锤头体外形的过渡段会出现气流分离ꎮ为探究各部分气动阻尼的变化ꎬ将整个箭体分为头部㊁过渡段㊁弹身3个部分ꎮ图1㊀表面网格及区域划分Fig.1㊀Surfacegridandregiondivision1.1㊀流场仿真模型本文分别用Reynolds平均法(Reynolds ̄averagedNavier ̄StokesꎬRANS)和改进的延迟分离涡模拟(improveddelayeddetached ̄eddysimulationꎬID ̄DES)[24 ̄25]进行计算ꎬ在RANS方程中ꎬ将变量分为平均值和波动值两部分ꎬ对于速度分量有ui=ui+uᶄi其中ꎬi=1ꎬ2ꎬ3ꎬui和uᶄi分别代表平均量和波动量ꎬ对于压强和其他标量也采用类似的形式ꎬ将这种形式代入连续性方程和动量方程中ꎬ并写成张量形式∂ρ∂t+∂∂xi(ρui)=0(1)∂∂t(ρui)+∂∂xj(ρuiuj)=∂p∂xi+∂∂xjμ∂ui∂xj+∂uj∂xi-23δij∂uk∂xkæèçöø÷éëêêùûúú+∂∂xj(-ρuᶄiuᶄj)(2)其中ꎬiꎬjꎬk可分别取1ꎬ2ꎬ3ꎻρ是密度ꎻt是时间ꎻ当i=j时δij取0ꎬ否则取1ꎮ式(1)㊁(2)是RANS方程ꎬ由方程可知RANS方法将湍流脉动对平均流动的作用模化为Reynolds应力项即-ρuᶄiuᶄjꎬ之后采用湍流模型进行封闭ꎬ本文采用的湍流模型为SSTk ̄ω模型ꎬ其输运方程为∂∂t(ρk)+∂∂xi(ρkui)=∂∂xjΓk∂k∂xjæèçöø÷+Gk-Yk∂∂t(ρω)+∂∂xi(ρωui)=∂∂xjΓω∂ω∂xjæèçöø÷+Gω-Yω其中ꎬk和ω分别代表湍流动能和湍流耗散率ꎬΓk和Γω分别代表k和ω的有效扩散系数ꎬGk和Gω分别代表k和ω的生成率ꎬYk和Yω分别代表k和ω的耗散率ꎮ因此RANS方法只能计算大尺度的平均流动ꎬ本文采用IDDES方法计算脉动压力对气动阻尼的影响ꎮIDDES方法是由分离涡模拟(detached ̄eddysimulationꎬDES)方法改进而来ꎬ其本质思想与DES方法相同ꎬ是想以网格尺度和模型中的特征尺度隐式划分RANS和大涡模拟(large ̄eddysimulationꎬLES)区域ꎬ使其既能处理RANS方法无法得到的脉动场ꎬ也能降低LES方法在模拟高Reynolds数流动时所需的计算资源ꎮ区别在于当边界层较厚或者分离区域较窄时ꎬDES方法会出现如模型应力损耗(modeledstressdepletionꎬMSD)ꎬ网格诱导分离(grid ̄inducedseparationꎬGIS)以及对数层不匹配(logarithmic ̄layermismatchꎬLLM)问题[24]ꎬ而IDDES模型通过改良计算区域划分ꎬ结合延迟分离涡模拟(delayeddetached ̄eddysimulationꎬDDES)和03第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析壁面模型大涡模拟(wall ̄modeledlarge ̄eddysimula ̄tionꎬWMLES)ꎬ定义新的长度尺度解决了这些问题ꎬ具体公式详见文献[25]ꎮ流场网格如图2㊁图3所示ꎬ边界层采用棱柱层结构ꎬ并调整第1层网格高度使得y+小于1ꎬ远场部分采用六面体结构网格ꎬ与边界层的过渡层采用非结构网格ꎮ整体网格单元数量为4.2ˑ106ꎮ图2㊀y方向截面网格示意图Fig.2㊀Schematicdiagramofcross ̄sectionalgridinthey ̄direction图3㊀x方向截面网格示意图Fig.3㊀Schematicdiagramofcross ̄sectionalgridinthex ̄direction物面边界条件为无滑移壁面条件ꎬ远场采用压力远场边界条件ꎬ湍流模型采用SSTk ̄ω模型ꎬ采用密度基求解ꎬ气体黏性采用Sutherland定律ꎬ空间离散采用2阶迎风格式ꎬ对流通量采用Roe格式ꎮ1.2 结构分析模型结构与流场耦合分析过程中ꎬ结构部分可以采用模态方法描述ꎮ结构模态可以通过有限元方法与结构模态试验方法获得ꎮ本文采用有限元分析结果获得的模态ꎬ图4所示为结构的前3阶模态ꎬ本文只分析计算结果中气动阻尼最小的第2阶模态ꎮ(a)f=1.200Hz(b)f=2.460Hz(c)f=2.957Hz图4㊀结构的前3阶模态Fig.4㊀Firstthreemodesofthestructure由于火箭结构外形简单ꎬ一般不考虑其扭转影响ꎬ因此可以将其简化为简单的梁模型ꎬ这样就可以给出其模态振动方程q㊆i+2biωiq˙i+ω2iqi=fi(3)式中ꎬqi为第i阶模态的广义位移ꎬbi为第i阶模态的结构阻尼系数ꎬωi为第i阶模态的固有频率ꎬ13气体物理2023年㊀第8卷fi为第i阶模态下质量归一化的广义气动力ꎮ若将fi按照Taylor展开并略去高阶项ꎬ可以将其转化为气动阻尼项与气动刚度项的形式ꎬ则式(3)可写为q㊆i+2(bi+Bi)ωiq˙i+(Ki+1)ω2iqi=0(4)式中ꎬBi为气动阻尼系数ꎬKi为气动刚度系数ꎬ研究表明[26]ꎬ气动刚度相对于结构刚度为小量可以忽略不计ꎬ而在计算中结构阻尼往往设置为0ꎬ因此气动阻尼可以直接反映其气弹稳定性ꎮ1.3㊀气动阻尼分析原理气动阻尼的分析可以采用强迫振动或者自由振动的方式进行ꎬ这两种方法获得的时域数据不同ꎬ提取气动阻尼的方式也不同ꎮ强迫振动方法初始演化过程较短ꎬ因此计算量较小ꎬ同时能够分析某一种振动形式的气动阻尼ꎬ明确该振动形式是收敛还是发散ꎮ分析过程中能够获得不同部位与部件的气动阻尼ꎮ但是对于多模态相互作用引起的发散(例如颤振)较难预测ꎮ自由振动方法需要一定的自由演化时间才能够对时域数据进行分析ꎬ不过自由振动方法能够获得最能够吸收能量的模态及其振动频率ꎮ对于本研究所关注的问题ꎬ气动载荷对结构振动的过程中气动阻尼的影响较大ꎬ而对气动刚度与气动质量影响较小ꎬ即结构的固有振动频率受到来流的影响较小ꎬ其稳定性问题主要由气动阻尼的正㊁负引起ꎬ所以采用强迫振动方法分析ꎮ强迫振动下结构做简谐模态振动qi(t)=Asin(ωit)式中ꎬA表示振动的振幅ꎬ将其代入计算气动力的公式中[21]并做正交积分可得Bi=ʏl0Bx(x)dx=-1MiAω2iTʏl0ʏt0+Tt0G(xꎬt)cos(ωit)dtdx(5)式中ꎬMi为第i阶模态的模态质量ꎬT为整数倍周期ꎬG为广义气动力ꎮ根据式(5)便可以得到局部或分区域的气动阻尼ꎮ1.4㊀耦合计算流程首先进行模态分析ꎬ以确定结构的模态频率与振型ꎬ用以设计强迫振动的频率和振幅ꎮ非定常流场计算前先进行定常流场计算ꎬ来加快非定常计算的演化速度并增强收敛性ꎬ结构节点位移通过径向基函数(RBF)插值方法[27]映射到气动网格节点上ꎬ来进行网格的变形ꎬ这里径向基函数选用WendlandC2ꎬ如下所示φ(x)=(1-x)4(4x+1)最后将计算出来的广义力提取出来ꎬ截取演化完毕的整数倍周期ꎬ进行气动阻尼计算ꎮ耦合计算流程图如图5所示ꎮ图5㊀耦合计算流程图Fig.5㊀Flowchartofcoupledcalculation2㊀结果分析与讨论2.1㊀流场分析结果计算的来流Mach数范围为0.7~1.2ꎮ其中中截面的压力分布如图6所示ꎮ可以看出在头部出现了膨胀波以及跨声速激波ꎬ在过渡段存在流动分离ꎬ随着Mach数的增大ꎬ头部低压区域逐渐扩张ꎬ并且能明显看到ꎬ在流动再附的位置产生了再附激波ꎮ(a)Ma=0.7023第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析(b)Ma=0.75(c)Ma=0.80(d)Ma=0.85(e)Ma=0.88(f)Ma=0.90(g)Ma=0.92(h)Ma=0.96(i)Ma=0.9833气体物理2023年㊀第8卷(j)Ma=1.00(k)Ma=1.05(l)Ma=1.10图6㊀不同Mach数下的中截面压力分布Fig.6㊀PressuredistributioninthemiddlesectionatdifferentMachnumbers2.2 气动阻尼分布通过上述流场分析ꎬ可以看出火箭不同部位流动结构并不相同ꎬ在头部与箭身上ꎬ流动主要为附着流动ꎬ而在过渡段会出现较为复杂的波系结构以及流动分离ꎮ针对不同的流动结构随流向站位x的变化ꎬ设该位置上广义力与广义位移的相位差为φ(x)ꎬ并且简谐振动没有引入其他模态的广义力ꎬ则广义力的表达式为G(xꎬt)=Fgen sin[ωt+φ(x)]+F0(6)其中ꎬFgen为广义力的振动幅度ꎬF0为广义力的常数偏移量ꎮ将式(6)代入到式(5)中得到B(x)=-FgenMAω2Tʏt0+Tt0sin[ωt+φ(x)]cos(ωt)dt其中ꎬ广义力的常数偏移量F0的积分为0ꎬ因此省略ꎮ通过将等式中的正弦函数部分进行和差化积得到B(x)=-FgenMAω2Tʏt0+Tt0sin(ωt)cos[φ(x)]cos(ωt)dt+[ʏt0+Tt0sin[φ(x)]cos(ωt)cos(ωt)dt](7)式(7)中第1部分在整个周期中的积分为0ꎬ只有第2部分保留ꎬ因此得到B(x)=-Fgensin[φ(x)]MAω2Tʏt0+Tt0cos2(ωt)dt(8)式(8)中积分部分恒为正值ꎬ决定整个气动阻尼的部分只有相位角φ(x)的正弦值sin[φ(x)]ꎬ为了能够更加直观地获得相位角与气动阻尼B之间的关系ꎬ须将符号转化为对应的正弦函数转角ꎬ根据正弦关系ꎬ此转角为πꎬ因此得到B(x)=-Fgen(x)sin[φ(x)+π]MAω2Tʏt0+Tt0cos2(ωt)dt(9)图7为气动阻尼变化曲线ꎬ可以看出随着Mach数的增大ꎬ整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ图7㊀有助推时气动阻尼变化曲线Fig.7㊀Aerodynamicdampingchangecurvewithboost根据式(9)ꎬ得到相位角与气动阻尼B之间的关系为:当φ(x)ɪ(-πꎬ0)时ꎬ相位角滞后ꎬ气动阻尼B为负值ꎻ当φ(x)ɪ(0ꎬπ)ꎬ相位角提前ꎬ43第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析气动阻尼B为正值ꎻ为当φ(x)=0时ꎬ无相位角差别ꎬ气动阻尼B为0ꎮ在过渡段上ꎬ复杂的波系结构以及流动分离ꎬ使得气动力与结构位移之间会出现较为明显的迟滞现象ꎬ从而导致相位角φ(x)ɪ(-πꎬ0)ꎬ由此在过渡段上产生了负的气动阻尼ꎮ计算过程中的广义力与广义位移随时间变化曲线如图8所示ꎬ可以看出所有工况计算结果都表现良好ꎬ需要注意的是在非定常计算初期ꎬ演化的不完全导致广义力存在一些突变异常的结果ꎬ计算气动阻尼时须剔除ꎬ选择后面演化完全的周期ꎮ本文计算了9个周期ꎬ剔除了第1个周期出现的错误结果ꎬ采用后8个周期进行气动阻尼分析ꎮ强迫运动振幅为芯级直径的0.5%ꎮ(a)Ma=0.70㊀㊀㊀(b)Ma=0.75(c)Ma=0.80㊀㊀㊀(d)Ma=0.85(e)Ma=0.88㊀㊀㊀(f)Ma=0.9053气体物理2023年㊀第8卷(g)Ma=0.92㊀㊀㊀(h)Ma=0.96(i)Ma=0.98㊀㊀㊀(j)Ma=1.00(k)Ma=1.05㊀㊀㊀(l)Ma=1.10图8㊀不同工况下的广义力与广义位移随时间变化曲线Fig.8㊀Timedependentcurvesofgeneralizedforceandgeneralizeddisplacementunderdifferentoperatingconditions2.3㊀气动阻尼影响因素2.3.1㊀有无助推对气动阻尼的影响捆绑式运载火箭相比于传统的运载火箭ꎬ最大的区别就是在尾部四周捆绑了助推器ꎬ使得其流场特性变得复杂ꎬ因此须分析其对气动阻尼的影响ꎮ图7㊁图9分别为有无助推时气动阻尼变化曲线ꎬ可以看出随着Mach数的增大整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过63第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ对比两个图可知ꎬ助推主要起增大气动阻尼的作用ꎮ还可以看出有无助推情况下头部的气动阻尼变化很小ꎬ意味着在箭体尾部施加控制很难影响到头部的气动阻尼ꎬ特别是在超声速流场中ꎮ图9㊀无助推时气动阻尼变化曲线Fig.9㊀Aerodynamicdampingchangecurvewithoutboost2.3.2㊀前节点位置影响为了考察前节点位置变化对气动阻尼的影响ꎬ在保持振动频率不变㊁头部最大振型位置与振幅不变的条件下移动前节点ꎬ变化后的振型如图10所示ꎮ(a)Frontnodeafterthetransitionregion(b)Frontnodeinthetransitionregion(c)Frontnodebeforethetransitionregion图10㊀前节点变化后的振型Fig.10㊀Vibrationmodeafterthechangeofformernode根据对计算结果的分析分别获得了不同前节点位置的整体气动阻尼对比与过渡段气动阻尼对比ꎬ如图11㊁图12所示ꎬ可以看出前节点位置的改变并没有影响整体气动阻尼随Mach数增大而增大的趋势ꎬ且前节点在过渡段上与过渡段前的整体气动阻尼相差不大ꎬ而前节点在过渡段后的整体气动阻尼要高于另两种情况ꎬ因此过渡段与头部放在同一侧有助于提高气动阻尼ꎮ过渡段的气动阻尼会随着前节点的变化发生剧烈改变ꎬ前节点在过渡段前后随Mach数增大的变化规律相反ꎬ节点前后的振动相位变化导致不同节点位置过渡段的振动相位不同ꎬ进而导致气动阻尼发生变化ꎮ图11㊀不同节点位置的整体气动阻尼Fig.11㊀Overallaerodynamicdampingatdifferentnodepositions图12㊀不同节点位置的过渡段气动阻尼Fig.12㊀Aerodynamicdampingofthetransitionregionatdifferentnodepositions2.3.3㊀强迫振动振幅大小对气动阻尼的影响为了考察强迫振动振幅大小对气动阻尼的影响ꎬ在保证流场结构不发生改变的前提下ꎬ振动振幅分别为原来的一半和两倍ꎬ根据工程经验ꎬ如果振幅超过芯级直径的5%ꎬ则须考虑流场结构改变所造成的影响ꎮ图13㊁图14分别为不同振幅下的整体与头部气动阻尼ꎮ73气体物理2023年㊀第8卷图13㊀不同振幅下整体气动阻尼Fig.13㊀Overallaerodynamicdampingatdifferentamplitudes图14㊀不同振幅下头部气动阻尼Fig.14㊀Aerodynamicdampingoftheheadregionatdifferentamplitudes可以发现改变振幅无论是对整体气动阻尼还是头部气动阻尼来说变化都很小ꎬ这意味着气动阻尼的大小主要取决于气动力与结构振动的相位差ꎬ不依赖于振动幅度的大小ꎮ2.3.4㊀脉动压力对气动阻尼的影响为了模拟出脉动压力的影响ꎬ采用IDDES方法对火箭气动阻尼进行计算ꎬ计算来流Mach数为0.92ꎬ计算过程中的广义力与广义位移如图15所示ꎬ相较于图8可以看出广义力随时间变化曲线并不光滑ꎬ脉动压力的存在导致广义力由多个频率叠加而成ꎮ由于第2阶模态的频率为2.46Hzꎬ而由分离流㊁激波振荡等引起的脉动压力频率往往远大于此频率ꎬ因此这里选择3.5Hz为分界ꎬ将高于3.5Hz的部分视为由抖振脉动压力引起的广义力ꎬ低于3.5Hz的部分视为强迫振动引起的广义力ꎬ通过低通滤波把高于3.5Hz的广义力滤掉ꎬ可以获得由强迫振动引起的广义力与广义位移变化曲线ꎬ如图16所示ꎬ通过此广义力计算的气动阻尼为2.08ɢꎮ同样地ꎬ进行高通滤波将低于3.5Hz的广义力滤掉ꎬ可以获得由抖振脉动压力引起的气动阻尼为(2.94ˑ10-3)ɢꎬ由此得到脉动压力引起的气动阻尼变化为0.14%ꎬ可以忽略不计ꎮ同时使用RANS方法计算的气动阻尼为2.07ɢꎬ与IDDES的计算结果相比误差约为(2.94ˑ10-3+2.08-2.07)/2.07ʈ0.48%ꎬ这说明针对气动阻尼的模拟ꎬ抖振引起的脉动压力对气动阻尼的计算结果影响很小ꎬ起主要作用的还是广义力的变化ꎬ该变化由强迫振动引起的结构边界变化所导致ꎮ图15㊀基于IDDES的广义力与广义位移变化曲线Fig.15㊀VariationcuresofgeneralizedforceandgeneralizeddisplacementbasedonIDDES图16㊀滤波后的广义力与广义位移变化曲线Fig.16㊀Variationcuresofgeneralizedforceandgeneralizeddisplacementvariationcurveafterfiltering3㊀结论本文通过数值计算方法研究了火箭的气动阻尼特性ꎮ根据流动特征分析与理论推导ꎬ发现火箭过渡段几何外形的收缩导致该区域出现复杂的分离与激波结构ꎬ从而造成了气动力相对于结构振动83第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析相位的滞后ꎬ导致了该区域为气动负阻尼ꎬ即气动不稳定性的主要来源ꎮ在此机理的基础上ꎬ分析了前节点位置㊁振动振幅㊁脉动压力等因素对气动阻尼的影响规律ꎮ可以得出以下结论:1)助推增加了正阻尼区域的面积ꎬ从而相对于没有助推的构型起到了增加气动阻尼的作用ꎮ2)前节点位置的改变对过渡段气动阻尼影响很大ꎬ节点前后的振动方向相反ꎬ导致节点在过渡段前后的气动阻尼变化规律也截然相反ꎬ将过渡段与头部区域放在节点的同一侧有助于增加气动阻尼ꎮ3)在不改变流场结构的前提下ꎬ改变振动的振幅ꎬ气动力也会产生相应幅度的变化ꎬ因此结构振幅对气动阻尼的影响可忽略不计ꎮ4)高频部分的广义力对气动阻尼的贡献很小ꎬ即结构振动引起的广义力变化对气动阻尼起主要作用ꎬ而脉动压力对计算气动阻尼影响不大ꎬ可忽略不计ꎮ参考文献(References)[1]㊀CoeCF.Steadyandfluctuatingpressuresattransonicspeedsontwospace ̄vehiclepayloadshape[R].NASATMX ̄503ꎬ1961.[2]ColeSRJrꎬHenningTLꎬRaineyAG.NASAspaceve ̄hicledesigncriteria[R].NASASP ̄8001(REV)ꎬ1964. [3]EricssonLEꎬRedingJP.Analysisofflowseparationeffectsonthedynamicsofalargespacebooster[J].Jour ̄nalofSpacecraftandRocketsꎬ1965ꎬ2(4):481 ̄490. [4]RedingJPꎬEricssonLE.Effectofaeroelasticconsidera ̄tionsonseasat ̄Apayloadshrouddesign[J].JournalofSpacecraftandRocketsꎬ1981ꎬ18(3):241 ̄247. [5]程镇煌.宇航飞行器跨音速气动弹性问题探讨[J].上海航天ꎬ1997(6):16 ̄21.ChengZH.Discussionontransonicpneumaticelasticityofspacecraft[J].AerospaceShanghaiꎬ1997(6):16 ̄21(inChinese).[6]倪嘉敏.我国运载火箭气动设计回顾[C].近代空气动力学研讨会论文集ꎬ2005.NiJM.ReviewofaerodynamicdesignofChinaᶄslaunchvehicle[C].Proceedingsofthemodernaerodynamicssymposiumꎬ2005(inChinese).[7]臧涛成ꎬ胡焕性.大长细比弹箭弹性效应研究综述[J].弹道学报ꎬ1999ꎬ11(3):89 ̄93ꎬ96.ZangTCꎬHuHX.Areviewofgreatslendernessratioprojectileelasticeffectresearch[J].JournalofBallisticsꎬ1999ꎬ11(3):89 ̄93ꎬ96(inChinese). [8]吴志刚ꎬ杨超.细长体弹箭的气动弹性问题与研究方法[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.WuZGꎬYangC.Aeroelasticityproblemsandresearchmethodsofslenderbodymissiles[C].9thNationalSym ̄posiumonAeroelasticityꎬ2005(inChinese). [9]张贺ꎬ黄晓鹏.弹性细长旋转弹箭运动稳定性问题的研究进展[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.ZhangHꎬHuangXP.Researchprogressonthestabilityofelasticslenderrotatingprojectiles[C].9thNationalSymposiumonAeroelasticityꎬ2005(inChinese). [10]ColeHAJr.Dynamicresponseofhammerheadlaunchvehiclestotransonicbuffeting[R].NASATND ̄1982ꎬ1963.[11]BartelsREꎬWiesemanCDꎬMineckRE.ComputationalaeroelasticanalysisoftheAreslaunchvehicleduringas ̄cent[R].AIAA2010 ̄4374ꎬ2010.[12]AzevedoJLF.Aeroelasticanalysisoflaunchvehiclesintransonicflight[J].JournalofSpacecraftandRocketsꎬ1989ꎬ26(1):14 ̄23.[13]SinclairAꎬFlowersG.Low ̄orderaeroelasticmodeloflaunch ̄vehicledynamics[R].AIAA2010 ̄7725ꎬ2010. [14]DotsonKW.Transientcouplingoflaunchvehiclebendingresponseswithaerodynamicflowstatevariations[J].JournalofSpacecraftandRocketsꎬ2001ꎬ38(1):97 ̄104.[15]ColeSRꎬHenningTL.Buffetresponseofahammerheadlaunchvehiclewind ̄tunnelmodel[J].JournalofSpacecraftandRocketsꎬ1992ꎬ29(3):379 ̄385.[16]崔尔杰.流固耦合力学研究与应用进展[C].钱学森科学贡献暨学术思想研讨会论文集ꎬ2001.CuiEJ.Researchandapplicationprogressoffluid ̄struc ̄tureinteractionmechanics[C].ProceedingsofSeminarofQianXuesenScientificContributionsandAcademicThoughtsꎬ2001(inChinese).[17]冯明溪ꎬ王志安.火箭跨音速动导数和抖振实验[J].宇航学报ꎬ1987(1):55 ̄63.FengMXꎬWangZA.Experimentsoftransonicderiva ̄tivesandbuffetingofrocket[J].JournalofAstronauticsꎬ1987(1):55 ̄63(inChinese).[18]白葵ꎬ冯明溪.弹性模型实验技术[J].流体力学实验与测量ꎬ1999ꎬ13(1):38 ̄42.BaiKꎬFengMX.Aeroelasticmodelandthebuffetex ̄perimentaltechnique[J].ExperimentsandMeasurementsinFluidMechanicsꎬ1999ꎬ13(1):38 ̄42(inChinese). [19]JiCꎬRanJHꎬLiFꎬetal.Theaerodynamicdamping93气体物理2023年㊀第8卷testofelasticlaunchvehiclemodelintransonicflow[C].Proceedingsofthe64thInternationalAstronauticalCon ̄gressꎬ2013.[20]季辰ꎬ吴彦森ꎬ何岗ꎬ等.运载火箭气动阻尼风洞试验研究[C].第十二届全国空气弹性学术交流会论文集ꎬ2011.JiCꎬWuYSꎬHeGꎬetal.Experimentalstudyonaero ̄dynamicdampingwindtunneloflaunchvehicle[C].Pro ̄ceedingsofthe12thNationalAeroelasticityConferenceꎬ2011(inChinese).[21]刘子强ꎬ白葵ꎬ毛国良ꎬ等.锤头体弹性振动跨音速气动阻尼系数的确定[J].宇航学报ꎬ2002ꎬ23(6):1 ̄7.LiuZQꎬBaiKꎬMaoGLꎬetal.Thedeterminationofaerodynamicdampingonhammerheadlaunchvehiclesattransonicspeeds[J].JournalofAstronauticsꎬ2002ꎬ23(6):1 ̄7(inChinese).[22]冉景洪ꎬ刘子强ꎬ胡静ꎬ等.减阻杆气动阻尼研究[J].力学学报ꎬ2014ꎬ46(4):636 ̄641.RanJHꎬLiuZQꎬHuJꎬetal.Researchofaero ̄dampingforbluntwithspike[J].ChineseJournalofThe ̄oreticalandAppliedMechanicsꎬ2014ꎬ46(4):636 ̄641(inChinese).[23]朱剑ꎬ冉景洪ꎬ吴彦森ꎬ等.捆绑式运载火箭的气动阻尼数值计算方法[C].第十三届全国空气弹性学术交流会论文集.哈尔滨:中国力学学会ꎬ中国空气动力学会ꎬ2013.ZhuJꎬRanJHꎬWuYSꎬetal.Numericalcalculationmethodforaerodynamicdampingofbundlelaunchvehicles[C].Proceedingsofthe13thNationalAeroelasticityCon ̄ference.Harbin:ChineseSocietyofTheoreticalandAp ̄pliedMechanicsꎬChinaAeromechanicsSocietyꎬ2013(inChinese).[24]SpalartPR.Detached ̄eddysimulation[J].AnnualReviewofFluidMechanicsꎬ2009ꎬ41:181 ̄202. [25]GritskevichMSꎬGarbarukAVꎬSchützeJꎬetal.Devel ̄opmentofDDESandIDDESformulationsforthek ̄ωshearstresstransportmodel[J].FlowꎬTurbulenceandCombustionꎬ2012ꎬ88(3):431 ̄449.[26]季辰ꎬ吴彦森ꎬ侯英昱ꎬ等.捆绑式运载火箭跨声速气动阻尼特性试验研究[J].实验流体力学ꎬ2020ꎬ34(6):24 ̄31.JiCꎬWuYSꎬHouYYꎬetal.Experimentalstudyofaerodynamicdampingcharacteristicsofalaunchvehiclewithboostersintransonicflow[J].JournalofExperimentsinFluidMechanicsꎬ2020ꎬ34(6):24 ̄31(inChinese). [27]AllenCꎬRendallTCS.UnifiedapproachtoCFD ̄CSDinterpolationandmeshmotionusingradialbasisfunctions[R].AIAA2007 ̄3804ꎬ2007.04。
土木建筑工程:结构力学考试答案(题库版)
土木建筑工程:结构力学考试答案(题库版)1、问答题潜体或浮体在重力G和浮力P的作用,会出现哪三种情况?正确答案:1)重力大于浮力,则物体下沉至底;2)重力等于浮力,则物体可在任一水深处维持平衡;3)重力小于浮力,则(江南博哥)物体浮出液体表面,直至液体下部分所排开的液体重量等于物体重量为止。
2、名词解释牛顿流体正确答案:在一维剪切流动情况下,当有剪切应力于定温下施加到两个相距dr 的流体平行层面并以相对速度dv运动,剪切应力与剪切速率成线性关系的流体称为牛顿流体.3、名词解释拉格朗日法正确答案:通过描述每一质点的运动达到了解流体运动的方法称为拉格朗日法。
4、问答题为什么流体静压强的方向必垂直作用面的内法线?正确答案:由于流体在静止时,不能承受拉力和切力,所以流体静压强的方向必然是沿着作用面的内法线方向。
5、填空题流体紊流运动的特征是(),处理方法是()。
正确答案:脉动现行;时均法6、单选在一个储水箱的侧面上、下安装有两只水银U形管测压计(如图二),当箱顶部压强p0=1个大气压时,两测压计水银柱高之差△h=h1-h2=760mm(Hg),如果顶部再压入一部分空气,使p0=2个大气压时。
则△h应为()A.△h=-760mm(Hg)B.△h=0mm(Hg)C.△h=760mm(Hg)D.△h=1520mm(Hg)正确答案:C7、填空题作用于曲面上的水静压力的铅直分力等于()。
正确答案:其压力体内的水重8、判断题粘滞性是引起流体运动能量损失的根本原因。
正确答案:对9、判断题当AB杆件刚度系数时,杆件的B端为定向支座。
正确答案:错10、问答题什么是水力半径?什么是当量直径?正确答案:水力半径R的定义为过流断面面积和湿周之比,即;当量直径。
11、名词解释水力光滑管与水力粗糙管.正确答案:流体在管内作紊流流动时,用符号△表示管壁绝对粗糙度,δ0表示粘性底层的厚度,则当δ0>△时,叫此时的管路为水力光滑管。
基于PumpLinx的变量燃料泵流动特性仿真研究
基于PumpLinx的变量燃料泵流动特性仿真研究王凯;孙涛;史小锋;伊寅;李永东;黄艳芬【摘要】To reveal the flow characteristics of the variable fuel pump for a torpedo, this paper proposes a hydraulic simulation method of the static variable cylinder valve axial piston pump by using the rotary cylinder plunger pump template in the professional pump computational fluid dynamics(CFD) simulation software PumpLinx, and further sim-ulates and analyzes the flow, pressure and cavitation of the pump. Simulation results show that this pump does not pro-duce cavitation, but there exists obvious pulse offlux with unsmooth transition. This research may benefit the improve-ment and design of the fuel pump with variable valve angle.%为了研究鱼雷变量燃料泵的流动特性, 文中提出了一种使用标准模版实现静缸式柱塞泵的研究方法, 通过3D建模软件UG与计算流体力学(CFD)仿真软件PumpLinx 完成对柱塞泵的仿真.在此基础上, 进一步对泵的流量、压力和空化进行了仿真和分析.仿真结果表明, 变量泵不会产生空化现象, 但流量脉动较大且高低压过渡不平稳, 这些结论对变量泵的改进设计有一定的参考价值.【期刊名称】《鱼雷技术》【年(卷),期】2015(023)006【总页数】5页(P444-448)【关键词】鱼雷;变量燃料泵;静缸式【作者】王凯;孙涛;史小锋;伊寅;李永东;黄艳芬【作者单位】中国船舶重工集团公司第705研究所,陕西西安, 710075;水下信息与控制重点实验室,陕西西安, 710075;中国船舶重工集团公司第705研究所,陕西西安, 710075;中国船舶重工集团公司第705研究所,陕西西安, 710075;中国船舶重工集团公司第705研究所,陕西西安, 710075;中国船舶重工集团公司第705研究所,陕西西安, 710075;北京海基嘉盛科技有限公司,北京, 610041【正文语种】中文【中图分类】TJ630.32动力系统是鱼雷的重要组成部分,它对鱼雷的航速、航深、可靠性等都有着重要影响[1]。
高效节能低污染脉动燃烧技术介绍及其应用
特 征 , 要关 注 热 声耦 合 起 振机 理 及相 应 的 主 声场 和 燃 烧 特 征 ; 迫 脉 动 燃烧 则侧 重 描 强
述 燃烧 对 声场 的响 应 行 为 。 一方 面 强迫 脉 动
它 的特 点 在 于 不 仅 可 以燃 用 气 、 燃 料 , 液 而 且还可 以方 便地使 用 固体燃 料 j 。 () ro n t : 4Mac n e 型 这种 类 型 的燃 烧器 最
排 放量 相 比, 脉动 燃烧 器 的排 放量 可 减小 到
分热 能 自发 地 转 化 成 声 能 , 生 自激 励 持 产
续 稳 定 的脉 动 , 由于 没 有 机 械 运 动 部 件 , 组
件 使 用 寿 命 相 对 较 长 , 点 是 噪 声 较 大 【。 缺 4 】 程 与脉动 压力 的准 则关 系 , 数学 表达式 为 :
根据 脉 动 燃烧 器 内声 场 发 生 方 式 分 为 自激 励 式脉 动燃 烧 脉
1 自激 励 脉 动 式 : 用 热 声转 换 将 部 、 利
流 的脉 动 削弱 了管 壁边 界 层 的 阻力 作用 , 使 热 量传递 加快 【。 3 】 () 3 降低污 染物 的排放 量 ( 别是 NO ) 特 x 充 分 的 燃 烧 使 得 排 气 中 的 C、H、C O 及 HC的含 量 都 降 低 , 由于脉 动 燃烧 的温 且 度 一般 较 低 , 得 热 力 型 的 N0x的 排 放 更 使 低, 常规锅 炉 5 与 8~ 1 8 3 mL / L的 N Ox
率 , 快 了化 学 反应 速 度 , 加 燃烧 充分 , 使得 脉
3脉 动燃 烧器 的分 类
动燃烧 器具有 很 高的燃烧 强度 【。 2 】
() 2 具有极 高 的传热 效率 脉动 燃 烧 室 中 的 声共 振 引起 了较 大 的 流 体振 动 , 附加 在 这 些振 动 上 的 是湍 流速 度 振 动 , 流 的作 用 增 强 了 传 热 。 湍 同时 由于 气
流-固耦合共振式液压滤波器性能研究
1
2
因而就不 可避免 地要产 生流量 脉动 与压力 脉动沿 管 路 传播 。介质脉 动所产 生 的交 变力使 管路及 它所 连 接 的 附件产 生振 动 , 强烈 的 振动 会使 管 路 附件 以及 它们 的连接 部分 松动 或破 裂 , 引起 系统 结 构 的振 并 动 。为 了对 液压 系统 中的振 动和 噪声 进 行抑 制 , 人 们研制了各种管路流体滤波器 , 常用的滤波器有共
中 图分 类 号 : H 3 T 17 文献 标 识 码 : A 文章 编 号 -17 — 9 4( 00)5 0 1—0 6 2 80 2 1 0 — 0 0 0 3
引 言
在 各 种流 体 管 路 系统 中 , 如液 压 系统 、 油 系 燃 统 、 却水 系统 等 , 冷 由于要周 期 性地 泵送 工 作介 质 ,
第 5期 ( 第 4 总 2期 ) 21 0 0年 9月
流 体钴幼与 控副
Fl i P w r Trns iso a d u d o e a m s i n n Co to nr】
N . S r l o 2 o5(e a . ) i N4
S p 2 1 e.00
流 一固耦 合 共振 式液 压滤 波器 性 能研 究
振型 、 容腔型 、 阻性 型[ 1 ] 。 共振 型液 压滤波器 是利用 共振来 吸收 能量 的原 理 而达 到脉动 消减 的 目的 。 这种 滤波器 直流压 损小 , 结 构简 单 , 在其 共振 频率 范 围 内衰减 效 果好 。容腔 型液 压滤波 器主要 减小有 效容积模 量或 增加有 效容
插入 损失 表示 为
见 :2 l 1 丝 ) ot 。 ) ( g
p2 1
管道 作腼ຫໍສະໝຸດ 积 为 () 8 (
离心泵蜗壳内压力脉动特性数值分析
离心泵蜗壳内压力脉动特性数值分析刘厚林;杜辉;董亮;吴贤芳;刘东喜【摘要】为揭示离心泵蜗壳流道内的压力脉动变化规律,采用雷诺时均方法(RANS),对3种工况下的离心泵内部三维非定常湍流流场进行数值计算,分析同一蜗壳断面不同位置以及沿蜗壳周向不同点的压力脉动特性.结果表明:蜗壳流道内具有非常明显的压力脉动,在各种工况下压力脉动的主频均是叶片通过频率;同一蜗壳断面上的压力脉动从蜗壳底部到蜗壳背面先减小后增大,蜗壳底部监测点的高频脉动成分较多;沿蜗壳周向,随着圆周角的增大,压力脉动减弱,隔舌附近压力脉动幅度最大,且高频脉动成分明显增加.%In order to reveal the variation laws of the pressure fluctuation in the volute of a centrifugal pump, the three-dimensional unsteady flow fields in the centrifugal pump under three different conditions are numerically solved by means of the Reynolds-averaged Navier-Stokes method. The pressure fluctuations at different points of the same section of the volute and different points along the circumferential direction of the volute are analyzed. The results show that the pressure fluctuations in the volute are obvious. The dominant frequencies under three conditions are found to be the ones that the blades pass. The pressure fluctuations at the same section decrease first and then increase from the bottom to the rear of the volute, and high-frequency fluctuations at the monitoring points on the bottom of the volute are evident. Along the circumferential direction, with the increase of circular angles, the pressure fluctuations decrease. The amplitude of pressure fluctuations isthe largest near the tongue where high-frequency fluctuations dramatically increase.【期刊名称】《水利水电科技进展》【年(卷),期】2013(033)001【总页数】5页(P18-21,32)【关键词】离心泵;蜗壳隔舌;压力脉动;雷诺时均法;数值分析【作者】刘厚林;杜辉;董亮;吴贤芳;刘东喜【作者单位】江苏大学流体机械工程技术研究中心,江苏镇江212013;江苏大学流体机械工程技术研究中心,江苏镇江212013;江苏大学流体机械工程技术研究中心,江苏镇江212013;江苏大学流体机械工程技术研究中心,江苏镇江212013;江苏大学流体机械工程技术研究中心,江苏镇江212013【正文语种】中文【中图分类】TV136+.2离心泵的空间非对称结构使其内部流动呈现出复杂的非定常特性,这种特性使泵在产生静态压力分量的同时还会产生动态压力分量,也就是压力脉动[1]。
脉动热管可视化实验研究进展
CHEMICAL INDUSTRY AND ENGINEERING PROGRESS 2018年第37卷第8期·2880·化 工 进展脉动热管可视化实验研究进展孙潇1,韩东阳1,焦波2,甘智华1,3,4(1浙江大学制冷与低温研究所,浙江 杭州 310027;2哈尔滨理工大学荣成校区机械工程系,山东 荣成 264300;3浙江省制冷与低温技术重点实验室,浙江 杭州 310027;4国家冷冻冷藏设备质量检验中心(河南),河南 商丘476800)摘要:脉动热管是一种结构简单、可用于微小空间、高热流密度条件下的传热元件。
其优良的传热特性源于内部工质的蒸发、冷凝相变以及两相流动等热-水动力学特性。
可视化实验是研究工质流动特性的常用方法,对于理解脉动热管的传热机理非常重要。
本文介绍了5种可视化技术手段,简述了近二十年来脉动热管可视化实验观察到的气液相变和多种流型,分析了流型的影响因素,并总结了流动特性和传热性能之间的对应关系。
指出目前可视化实验观察到的结果可以定性地解释一些传热现象,却难以定量地揭示脉动热管的传热规律。
未来可以将可视化实验中测得的气液塞运动速度、液膜厚度和相分布等信息利用起来,考虑管壁材质的影响,建立更准确的气液两相流模型,用来设计、评估和预测脉动热管的性能。
关键词:脉动热管;可视化;气液两相流;相变;传热中图分类号:TK172.4 文献标志码:A 文章编号:1000–6613(2018)08–2880–12 DOI :10.16085/j.issn.1000-6613.2017-2168Research progress on visualization of pulsating heat pipesSUN Xiao 1, HAN Dongyang 1, JIAO Bo 2, GAN Zhihua 1, 3, 4(1Institute of Refrigeration and Cryogenics, Zhejiang University, Hangzhou 310027, Zhejiang, China ;2Department of Mechanical Engineering, Rongcheng Campus, Harbin University of Science and Technology, Rongcheng 264300, Shandong, China ;3Key Laboratory of Refrigeration and Cryogenic Technology of Zhejiang Province, Hangzhou 310027,Zhejiang, China ;4National Quality Inspection Center of Refrigeration Equipment (Henan), Shangqiu 476800,Henan, China )Abstract :Pulsating heat pipe (PHP )is a heat transfer device with simple structure capable of transferring heat effectively in comparatively small space. Visualization experiment is a common method to study flow characteristics of working fluid, which is very important to understand thermo-hydrodynamics of PHP. Previous studies have shown that bubble flow, slug flow, annular flow or semi-annular flow are common flow patterns. Phenomena such as nucleation boiling, coalescence of bubbles, formation of vapor plug and liquid slug were observed in the closed loop PHP. In this paper, five visualization techniques were introduced. Gas-liquid phase transition and flow patterns observed in the visualization experiments of PHPs in recent two decades were briefly introduced. The influencing factors of flow pattern were analyzed, and the correspondence between flow characteristics and heat transfer performance were summarized. At last, future developments were suggested.Key words : pulsating heat pipe; visualization; gas-liquid flow; phase change; heat transfer管。
布袋除尘器扩容流场均匀性研究与应用
科技与创新┃Science and Technology&Innovation ·162·2023年第24期文章编号:2095-6835(2023)24-0162-02布袋除尘器扩容流场均匀性研究与应用吴连营(江苏徐矿综合利用发电有限公司,江苏徐州221000)摘要:某300MW循环流化床锅炉因燃料的改变,机组烟气中含尘激增,造成布袋除尘器在实际运行中系统阻力增加,在提出通过除尘器一仓室扩容改造方案后,为了验证方案的可行性,了解扩容后除尘器内部流场分布以及除尘器风门风速,利用CFD(Computational Fluid Dynamics,计算流体动力学)研究方法和FLUENT软件设计,对除尘器进行建模分析,提出增加灰斗进风导流板方法,消除流场不均、风速过高滤袋磨损现象,验证改造项目的可靠性,为技术改造提供理论依据。
关键词:300MW流化床锅炉;布袋除尘器;扩容改造;导流板中图分类号:TK223.1文献标志码:A DOI:10.15913/ki.kjycx.2023.24.048某锅炉是由东方锅炉(集团)有限公司制造的型号为DG1025/17.5-Ⅱ2的300MW亚临界循环流化床锅炉。
锅炉设置2台布袋除尘器,每台布袋除尘器滤袋4032只,每台锅炉共计8064只,采用压缩空气离线低压脉冲清灰,设计运行阻力1300Pa。
该公司近几年加大煤泥和污泥掺烧比例,由于燃料的改变,机组烟气中含尘激增,造成布袋除尘器在实际运行中系统阻力增加,经常在1600~2500Pa。
另外,滤袋底部磨损严重,造成除尘器出口烟尘排放质量浓度升高,使布置在除尘器后的低温省煤器磨损、积灰加剧,给机组的安全运行带来隐患。
为解决滤袋磨损问题,电厂利用CFD研究方法,对除尘器流场分布情况和烟气流速进行分析研究,提出增加灰斗进风导流板方法,消除流场不均、风速过高滤袋磨损现象,验证改造项目的可靠性,为技术改造提供理论依据。
脉动热管技术的研究进展_朱晓琼_201508
抽真空管 冷 凝 段
注液管 液塞
绝 热 段
蒸 发 段
单向阀
气泡
闭式回路型
带单向阀的 闭式回路型
弯头
图 1 脉动热管示意图
开路型
2 脉动热管的实验研究现状
2.1 脉动热管的可视化研究 通过可视化实验观察脉动热管启动和稳定运行时的 工质流态是一种常用的研究方法。Tong 等人[5-7]发现在 热管蒸发段可以观察到核态沸腾现象及气泡的膨胀、集 合等行为,且在加热量达到一定程度后工质在热管内会 形成定向循环流动。Yang 等[8]发现在充液率为 40%~ 70%时,低加热功率时只有液弹的振荡运动,高加热功率 时逐渐变为循环流动。Xu 等人[7]发现工质的种类会影响 气泡在脉动热管内的速度及位置,屈健等[9]则发现管径、 倾角及弯头数都会影响脉动热管内工质的流动状态,从 而影响换热。 曲伟等 在 [10] 半可视化脉动热管实验台发现,管内工 质的流型与加热功率有关,功率较小时为间歇振动,功率 较大时为单向脉动流动。同时发现热阻与加热器功率、 热管倾角及不凝性气体含量有关。杨蔚原等人[6]通过可 视化实验发现充液率、温差、对流换热条件、加热段与冷 却段的长度比例、管内表面状况都对脉动热管的稳定性 有影响。曹小林等人[11]在可视化实验中发现脉动热管存 在传热极限,且在最佳充液体率为 50%和最佳倾角为 50° 时的传热极限最高,当热流密度较大时,通道形状对传热 性能的影响较小。 一些学者也采用了特殊的探测技术对脉动热管进行 了可视化研究。Borgmeyer 等 采 [12] 用中子成像技术对铜 管式脉动热管进行了可视化研究,发现加热功率会影响 热管内的工质脉动,随着加热功率的增大,工质脉动会逐 渐趋于稳定。李志宏等人 通过 [13-14] 电容成像方法,测量 了脉动热管内两相流的液膜厚度,发现脉动热管流型呈 现出弹状流、混合流和环状流等。 总体而言,脉动热管内的工质流动及流型比较复杂, 在不同工况下会出现多种流型,伴随着气泡的生成、膨 胀、聚合、凝结与消失,流型与传热性能密切相关,随着加
基于CFD的加工中心吸尘系统管路仿真与优化
计算与设计流量 之差 /(m3/h)
计算与设 计流量误
差 /%
排气口 1 12029.1 12000.0
29.1
0.24
排气口 2 6025.5
6000.0
25.5
0.42
吸气 1
5278.8
5000.0
278.8
5.58
吸气 2
5568.8
3600.0
1968.8
54.69
吸气 3 吸气 4 吸气 5 吸气 6 吸气 7
吸尘装置一般由吸尘罩、吸尘管路、吸尘器和风机等 部分构成。许多科技工作者对吸尘装置开展了一系列数值 模拟和优化研究,如林钰珍等人通过数值模拟对铣床除尘 系统吸风罩进行了优化设计,发现当吸风罩张角为 60°时, 吸尘效果最佳;王沣浩等人对影响吸尘器效果的除尘速度、 进风口尺寸以及风口到滤网的距离等参数进行了优化分析; 张亚蕊人等对常规滤筒内部增加锥体结构进行了数值模拟 仿真,发现改进后的滤筒清灰效果明显提高,并且延长了 滤筒寿命 [3-5]。目前,技术人员对除尘装置的研究主要集中 在除尘器与吸尘罩方面,而对吸尘管路的研究较少,所以 本文对吸尘管路管径进行优化具有一定的应用价值。
吸尘系统会引起空气流动,当吸尘速度稳定时,流动 状态与时间无关,可以当作稳态来计算。经计算得,马赫 数为 0.051,远小于 0.3,可认为空气是不可压缩流体。然 后对于湍流数值模拟,采用标准 k-ε 两方程模型;近壁面 采用壁面函数法,估算距壁面第一层网格高度值为 1.2mm, 且边界层网格数设为 5 层。
对此加工中心的吸尘系统的几何模型进行简化,重点 细化吸气口所包围区域的几何轮廓,尽可能与实际几何模 型保持一致,减小模型简化对气流运动所造成的影响。分 析模型的网格共有 770 万个混合网格单元,网格模型如
输气工高级考试(试卷编号151)
输气工高级考试(试卷编号151)1.[单选题]工业设备上的铭牌所标的额定电压、电流均指( )。
A)平均值B)有效值C)最大值D)最小值答案:B解析:无2.[单选题]当压力变送器发生无输出故障时,下面解决办法中哪项是错误的( )。
A)把电源极性接正确B)必须保证供给变送器的电源电压小于12VC)更换表头D)把电源线接在电源接线端子上答案:B解析:无3.[单选题]测量脉动压力时,最大工作压力不超过量程的1/2;测量高压力时,最大工作力不应超过量程的A)2/3B)3/4C)3/5D)4/5答案:C解析:无4.[单选题]自动控制系统安全措施主要是防止雷击等自然灾害损坏系统和人的()造成的损坏。
A)误操作B)不安全行为C)故意破坏D)技术低下答案:B解析:无5.[单选题]输气管道中气流温度的变化,取决于( )和与周围的热交换。
A)起点气流温度B)环境温度C)运动的物理条件解析:无6.[单选题]输气管道的压降比是( )。
A)进口压力/出口压力B)出口压力/进口压力C)实际压力降/理论压力降D)理论压力降/实际压力降答案:C解析:无7.[单选题]()具有承载能力大、传动平稳、使用寿命长等特点。
A)直齿圆柱齿轮B)斜齿圆柱齿轮C)圆锥齿轮D)伞齿轮答案:B解析:8.[单选题]材料在受力后是否发生破坏,取决于构件的应力是否超过( )。
A)抗拉强度B)极限抗力C)剪切应力D)屈服极限答案:B解析:无9.[单选题]甲醇是( )的物质。
A)易燃但无毒B)对视神经无伤害,喝人甲醇后都不会失明C)易燃且有毒D)易燃,高度挥发的液体但不能与水相混溶答案:C解析:无10.[单选题]油浸变压器在正常情况下为使绝缘油不致过速氧化,上层油温不不宜超过()。
A)75℃B)85℃C)95℃D)105℃答案:B11.[单选题]设备信息管理模块的主要功能是利用计算机硬件、软件、网络设备、通信设备以及其他自动化设备,进行设备信息的收集、传输、加工、储存、更新和()。
搅拌釜生物反应器中污泥絮凝曝气的CFD多相种群平衡耦合模型的联合仿真
搅拌釜生物反应器中污泥絮凝曝气的CFD多相种群平衡耦合模型的联合仿真摘要:曝气搅拌槽生物反应器通常用于通过微生物环境转化有机物和去除营养物质。
基于气(气)固(泥)群体平衡模型(pbm)的多相流计算流体动力学CFD(Computational Fluid Dynamic CFD)在不同时间尺度上模拟复杂的多物理现象。
联合模拟耦合CFD-PBM 动力学模型,捕捉生物质生长动力学和反应器水力动力学对整个过程的影响。
絮体粒径分布是至关重要的,因为它决定了最终絮体结构的不同等级,分别用作商业好氧废水处理。
耦合群体平衡方程采用连续偏微分方程的非均相离散方法求解污泥絮体,气相采用标准矩法求解,二者同时求解同时求解,影响生物质动力学的增长。
此外,在全瞬态CFD分析中,还研究了搅拌器转速和空气流速对生物量动态增长的影响。
通过用户定义函数(UDF),模型中使用了生物量增长率的动力学,作为氧饱和度和空气/流量平衡的函数。
多相系统方法学中的多物理综合模拟(co-simulation)是研究生物反应器中生物量生长各个参数影响的一个有价值的工具。
引言:群体平衡模型(PBM)在生物制药、污水处理厂、食品加工、生物反应器、细菌生长培养等领域有着广泛的应用(Zhang,2009;ANSYS Theory guide,2015;Qian Li,2017)。
在生物化学工业中,种群平衡模型对颗粒或气泡尺寸分布(PSD)的确定起着重要作用。
(Ding和Biggs,2006)使用离散PBM模型检验了活性污泥絮凝,确定了剪切和颗粒大小对碰撞效率的破碎率系数影响。
平均流速梯度增大,絮体粒径减小。
过程是周期性的,絮体的大小随着时间的推移而变化,聚集和分解是唯一的过程机制。
在不假设悬浮液的流变依赖性与碎裂、异凝、异聚、颗粒吸附或核化现象之间的额外相互作用的情况下,简化了文献中的种群平衡模型的大小(Heath等人,2003;Oshinowo等人,2016;Chen等人,2004)。
平板壁面湍流脉动压力及其波数-频率谱的大涡模拟计算分析研究
平板壁面湍流脉动压力及其波数-频率谱的大涡模拟计算分析研究张晓龙;张楠;吴宝山【摘要】Turbulent wall pressure fluctuations beneath turbulent boundary layers are important sources of flow noise. The computation of wall pressure fluctuation and its wavenumber-frequency spectrum is a hot topic in the field of flow-acoustic coupling. It is necessary to carry out corresponding research. In this pa-per, wall pressure fluctuation and its wavenumber-frequency spectrum are computed using large eddy sim-ulation (LES) with approriate sub-grid scale model, grid number and discretization methods. The results are compared with the experiment of Abraham and discussed in detail. Firstly, some fundamentals of the nu-merical simulation are presented, including the philosophy of LES, formulations of sub-grid scale models, discretization methods and boundary conditions, etc. Secondly, the rectangular test section of Abraham’s experiment and its computational domain are depicted. Thirdly, the scaling of turbulent wall pressure fluc-tuations and its variations due to different free stream velocities are discussed. The wavenumber-frequency spectra of turbulent wall pressure fluctuations are computed by taking FFT of the correlation function of the simulation data, the computed spectra of the wall pressure fluctuations are compared with those of Abra- ham’s experiment and analyzed qualitatively and quantitatively. Finally, comparison of typical theoretical models ofwavenumber-frequency spectra is made based on computational and Abraham’s experimental results. Groundwork i s made for further research in turbulent wall pressure fluctuation and its wavenumberfrequency spectrum.%壁面湍流脉动压力是重要的流噪声声源,对壁面湍流脉动压力及其波数-频率谱进行数值计算是流声耦合领域的重要课题,开展相应的研究十分必要。
中低温脉动热管传热性能的模拟研究
中低温脉动热管传热性能的模拟研究马文统;陈曦;曹广亮;武飞【摘要】Pulsating heat pipe(PHP)is a high-efficiency heat transfer device with simple structure and excellent heat-transfer performance. In order to study the pulsating heat pipe which can be used in cryogenic refrigerator at the liquid R508B zone,a three-dimensional numerical model was established and numerically analyzed based on the VOF method in the paper. The filling ratio(30%,50%,and 70%)and input power(20 W,40 W,60 W,80 W,100 W and 120 W) were applied in present numerical investigation. The distributions of vapor-liquid in PHP after the initial of working fluid and the temperature distribution were simulated. The influence of filling ratio and input power effect on the performance of heat transfer was discussed. The simulation results show that the characteristics of the flow in the R508B zone are simi-lar with the room temperature PHP,when the heat pipe working at the stabilization stage,the fluctuation of temperature is periodic.The heat transfer performance of pulsating heat pipe is changed with the increase of heating power,when the fill-ing ratio is low,the temperature difference and the equivalent thermal resistance of cold and hot part are relatively small in the low heating power,the equivalent thermal resistance will increase in the high heating power. When the filling ratio is high,the temperature difference and the equivalent thermal resistance of cold and hot part are relatively big in the low heating power,and the equivalent thermal resistance will decrease with theincreasing of heating power.%脉动热管是一种新型传热元件,具有结构简单,传热性能突出的优点。
基于流固耦合的风力机气动噪声影响研究
基于流固耦合的风力机气动噪声影响研究郭茂丰;张立茹;王占洋;李得银【摘要】考虑风力机叶片与空气的流固耦合作用,基于ANSYS workbench工作平台,采用双向流固耦合的方法,模拟预测风力机的气动噪声,并与额定工况下的实验数据对比.结果发现:耦合作用下风力机气动噪声增大,且耦合模拟得到的气动噪声声压级与实验值更为接近,证明计算模型的准确性;风轮后的辐射声最大声压级在叶片径向0.57R~0.71R位置,风力机叶片与空气的流固耦合作用,增大了辐射声的声压级,而对于辐射声的传播规律影响很小;耦合作用下随尖速比的增加,相同位置气动噪声的声压级呈现缓慢增大的趋势,不同尖速比下气动噪声的声压级随轴向距离的增加变化规律大致相同,均呈不断减小的趋势.【期刊名称】《可再生能源》【年(卷),期】2018(036)008【总页数】6页(P1238-1243)【关键词】水平轴风力机;流固耦合;气动噪声;声压级【作者】郭茂丰;张立茹;王占洋;李得银【作者单位】北方联合电力有限责任公司呼和浩特热电厂, 内蒙古呼和浩特010030;内蒙古工业大学能源与动力工程学院, 内蒙古呼和浩特 010051;内蒙古工业大学能源与动力工程学院, 内蒙古呼和浩特 010051;内蒙古工业大学能源与动力工程学院, 内蒙古呼和浩特 010051【正文语种】中文【中图分类】TK830 引言由于风力发电装机容量的大幅增加及风电场规模的不断扩大,风力发电机组日益临近居民区,由此引发的问题也逐渐突显,风力机的气动噪声成为风电场设计及运行人员研究的重点[1]。
风力机在复杂多变的风工况下运行,使得风轮承受气动力、离心力等交变载荷的作用,这些载荷使得风力机叶片发生变形,同时这些变形反过来又会改变风的运动,从而造成风力机流场的变化,因流场与声场相互关联,进而引起气动噪声的变化。
因此,考虑风力机流场与结构场的相互耦合,分析在流固耦合作用下的气动噪声显得尤为重要。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第36卷第8期煤 炭 学 报Vol.36 No.8 2011年8月JOURNAL OF CHINA COAL SOCIETYAug. 2011 文章编号:0253-9993(2011)08-1412-05基于脉动压力时均值的除尘器气液耦合研究李小川1,2,胡亚非1,蒋炎炎1,熊建军1(1.中国矿业大学化工学院,江苏徐州 221008;2.煤炭加工与高效洁净利用教育部重点实验室,江苏徐州 221008)摘 要:以节流型自激式水幕除尘器为研究对象,通过紊流状况下除尘器脉动压力时均值的测量,对除尘器内部气-液耦合状况进行了研究,得出了接触腔的高度直接影响除尘器气-水两相耦合程度进而影响除尘器除尘效果,接触腔合适高度为300mm ;接触腔压力沿流线方向逐渐降低,随着节流强度α值增大接触腔中压力整体下降,漩涡区的涡强也随之增大;随α值变化气-水耦合形式主要表现为水流增加节流阻力㊁激发除尘水滴㊁改变气流流向以及3种形式共同作用,对于400mm×300mm (宽×高)的除尘器,当α值在250/300~235/300时为3种形式共同作用,具有较好的除尘效果㊂关键词:气液耦合;节流型自激式水幕除尘;时均压力;节流强度中图分类号:TD714.45 文献标志码:A收稿日期:2010-11-06 责任编辑:许书阁 基金项目:国家高技术研究发展计划(863)资助项目(2007AA05Z339) 作者简介:李小川(1984 ),男,四川大英人,博士研究生㊂E-mail:xiaochuanli2008@Study of gas⁃liquid coupling in dust separator based ontime⁃averaged value of pulsatile pressureLI Xiao⁃chuan 1,2,HU Ya⁃fei 1,JIANG Yan⁃yan 1,XIONG Jian⁃jun 1(1.School of Chemical Engineering and Technology ,China University of Mining and Technology ,Xuzhou 221008,China ;2.Key Laboratory of Coal Process⁃ing and Efficient Utilization ,Ministry of Education ,Xuzhou 221008,China )Abstract :Study object is the throttle⁃type self⁃excitation dust collection.Through the measurement of time⁃averaged value of pulsatile pressure on turbulence flow,gas⁃liquid coupling status in dust separator was studied.The primary conclusions are made as follows:the dedusting effect is influenced by the height of the catch⁃cavity because the gas⁃liquid coupling status is changed with the height of the catch⁃cavity.And the suitable height of catch⁃cavity is 300mm.Pressure of catch⁃cavity reduce gradually along with the direction of flown line.As the throttle strength value αaggrandizing,the whole pressure value of catch⁃cavity drops,and the vortex intensity produced in whirlpool area is also reduced.As the changing of αvalue,the water injected in the catch⁃cavity will influence the gas⁃liquid coupling of catch⁃cavity.There are several influence formats which are increasing throttle resistance,exciting dust removal drip,transferring the airflow direction and the interaction of all.If it has the size of 400mm×300mm,the dust collection has the third influence format and with good dedusting effect when the αvalue between 250/300~235/300.Key words :gas⁃liquid coupling;throttle⁃type self⁃excitation dust collecting;time⁃averaged pressure value;throttlestrength 在煤炭开采㊁洗选㊁利用等环节产生大量粉尘,并有工序多㊁尘源分散㊁粉尘分散度高等特点,严重威胁到安全生产和操作人员的身体健康[1-5]㊂一些早期科研工作者以选煤系统和煤矿掘进工作面为背景,对湿式除尘技术的除尘机理㊁除尘效率等进行了相关研究和应用探讨[6-8]㊂近年来,科研工作者在气-液两相流的理论与工程实践方面取得了不少研究成果,包括气-液两相流气含率㊁流型判定㊁流量测定㊁绕流压第8期李小川等:基于脉动压力时均值的除尘器气液耦合研究降特性㊁气-液两相流的影响因素㊁数值模拟等[9-15],然而除尘器内部气-液耦合状况等并未得到深入研究㊂笔者以节流型自激式水幕除尘器为对象,通过测量紊流状况下除尘器脉动压力时均值,对其内部气-液耦合状况进行研究,以期更深入地了解该类除尘器的特性㊂1 除尘器原理及压力脉动特性节流型自激式水幕除尘器结构原理如图1所示,含尘气流在通风机吸力作用下经进气口进入除尘器,气流经节流引导装置加速后在节流口对倾斜流道中自由流动的除尘水流形成冲击,产生高速运动的除尘水滴,大量水滴以一定角度向接触腔运动形成除尘水幕,捕捉固相粉尘,最后洁净气体经排气口排出[16]㊂图1 节流型自激式水幕除尘器结构Fig.1 Structure chart of the dust separator1.1 压力脉动特性除尘器内部气流处于紊流状态,在理想状况下,高速紊流气流与水流连续均匀的对流,随着补给水的流入,除尘水位置势能具有阶跃特性,即当气流对水流的冲击力正好等于水的位能时最高处的水流微粒处于平衡状态,只有当水位能增加到一定程度后这种平衡才被打破,水微粒克服气流冲击,迅速溃流到下游,此时流失的水多于处于平衡状态后增加的水量,气流的冲力大于水流压力,气流开始蓄势,水重新增加,此时除尘器内部压力可简化为周期性脉动压力,即p 1=2A (ft -[ft ])-A(1)式中,A 为振幅;方括号为高斯函数;f 为频率,Hz;t 为时间,s㊂另外,由于电力波动㊁风机运行不稳定和外流场干扰等带入压力随机波动,记为p 2㊂因此,除尘器的脉动压力可表示为周期性脉动压力和随机波动之和,即f (p )=p 1+p 2(2) 在较短时间内,周期性脉动是对除尘器内部压力的主要影响因数,在充分避免干扰的情况下,电力波动等随机波动在较长时间内变化趋势才较为明显,由此,通过对d t 时间内连续测量的压力值进行积分并取其时间平均,即可得到表征此段时间的时均压力值,即p -=[∫tf (p )d t ]/t =p -1+[∫tp 2d t ]/t(3)1.2 节流强度节流强度α是节流口加速含尘气流能力的度量,用节流引导装置始端水平高度B 与末端水平高度b 表示,即α=(B -b )/B(4)2 实验方案及数据实验1:两台横截面积相同,截面尺寸分别为(宽×高)400mm×300mm(定义为I 号)和200mm×600mm(定义为II 号)的除尘器,设定适当的除尘风机运行频率,实时测量不同α值时除尘器流量㊁压力损耗,并观察水幕效果,结果见表1㊂表1 不同节流强度条件下除尘器水幕效果㊁风机流量和压力损耗Table 1 Water screen impression ,flow rate and press expense of dust separator on different throttle strength condition横截面节流强度/(mm㊃mm -1)入口风速/(m㊃s -1)流量/(m 3㊃h -1)压阻/Pa 水幕效果270/3004.311861.92851水滴细小且较少,幕状流动限于接触腔中下部I 号250/3004.932129.76734水幕较好的充满接触腔,气㊁水接触充分220/3005.492371.68680水滴减少,流线清晰,幕状流动明显减弱200/3006.082626.56632水滴产生较少570/6002.531092.961190水滴细小且较少,幕状流动限于接触腔下部II 号550/6003.941702.08843幕状流动限于接触腔中下部,气㊁水分层520/6004.762056.32721水滴减少,流线清晰,幕状流动限于较小区域500/6005.132216.16703水滴进一步减少,幕状流动限于较小区域 实验2:在空气介质和气-水介质条件下,对I 号除尘器接触腔内部压力分布进行测量,接触腔压力测3141煤 炭 学 报2011年第36卷孔布置如图2所示,试验结果如图3㊁4所示㊂图2 除尘器压力测孔布置位置Fig.2 Installing place of pressure sensor3 结果分析与讨论3.1 接触腔高度对气-水耦合程度的影响由表1可得,随着α值减小,除尘器流量增加,阻力降低,气-水耦合程度先增加后降低(接触腔中水幕到达的高度先升高后降低),当α减小到一定值后,气流对来流基本不产生冲击作用,气-水为单相耦合㊂不同除尘器气-水耦合最佳程度的α值不同,当I 号除尘器α=250/300时水幕较好的充满接触腔顶部,气㊁水充分耦合,当II 号除尘器α=550/600时水幕达到最佳状态,但幕状流动仅限于接触腔中下部,气-水单相耦合㊂相比II 号,I 号除尘器在同一板高条件下流量较大㊁阻力较小㊂由此可知,接触腔高度直接影响除尘器气-水两相耦合程度进而影响除尘器除尘效果,接触腔高度很低时,除尘器处理量受到限制;接触腔高度很高时,水图3 不同α值和不同介质的压力对比Fig.3 Pressure comparison on different αvalues and different medium conditions幕不能到达接触腔上部,大量含尘气流绕开水幕排出,部分粉尘无法捕捉;只有接触腔高度合适(约300mm)才能使水幕到达接触腔上部,气-水充分耦合,且达到较大处理量和较小阻力㊂3.2 除尘器内部气-水耦合规律及影响因素除尘器内部介质性质㊁节流强度等都不同程度地影响着除尘器内部气-水耦合规律,对这些规律的研究有利于更深入地了解除尘器的特性㊂3.2.1 除尘腔内部压力分布规律实验2得到的不同α值下接触腔各点压力分布如图3所示,不同介质下接触腔压力随α值变化如图4所示㊂对同一α值接触腔在空气介质和气-水介质的压力分布趋势相似,压力分别在点2㊁5㊁8有极小值,在3㊁7点有极大值㊂接触腔中压力具有递变趋势,5㊁6㊁7点压力有逐渐升高的趋势,1㊁2㊁3㊁4点压力有逐渐降低的趋势,但在2点压力明显偏低于趋势线㊂部分点在两种不同介质下有较大偏离,对空气介质,随着α值增大各点压力都有所升高但变化较为平缓;气-水介质各点压力变化不同,当α较小时各点压力变化趋势与空气介质时接近,随着α值增大,1㊁2㊁9点压力升高,5点降低,4㊁6㊁7㊁8点先降低后有微弱上升趋势,3点先升高后降低,且随着α值增大,5㊁6及3㊁4㊁7㊁8㊁9点压力分别先后降到空气介质时压力以下,5点随着α值增大压力变化幅度最大㊂3.2.2 节流强度对接触腔中压力分布的影响由空气介质条件下接触腔中压力随α变化情况可以看出,节流强度直接影响接触腔中压力分布㊂由图4(a)可以看出,随着α值减小,接触腔中各点的压力都升高,主要原因是节流引导装置的节流阻力降4141第8期李小川等:基于脉动压力时均值的除尘器气液耦合研究图4 各点压力随α值变化情况Fig.4 Pressure change αvalue低,除尘风机负载阻力下降㊂除5㊁9点外,各点压力p (Pa)与α近似呈线性关系,可表示为p =-kB (1-α)(5)式中,k 为比例系数㊂同一流线上的各点压力有升高趋势的原因可由伯努利方程解释,即p a ρg +v 2a2g =p b ρg +v 2b 2g(6)其中,p a ㊁p b 和v a ㊁v b 分别为流线上a ㊁b 两点的压力(Pa)和流速(m /s)㊂当a 点取节流口时,压力最小,速度最大,且随着α的增大,压力增大,速度降低,沿流线前进,b 点位于容积较大的接触腔,由此流体质点速度降低,压力升高㊂对于2㊁5㊁9点压力明显高于附近其它点,可从除尘腔结构方面考虑为这3点分别由于不连续过渡和结构死区造成的气流漩涡,涡强为Г=v 2/2r(7)其中,v 为气体质点流速,m /s;r 为质点离涡流核心区的距离,m㊂将式(7)代入式(6),可得,离涡流核心区越远,压力越大,反之越远㊂因此当α越小,相同r 值处v 和Г值较大,该处压力相对较低,反之较高㊂3.2.3 介质变化对接触腔中压力分布的影响接触腔的压力分布在气-水介质条件下较空气介质下有明显的变化㊂当α=270/300时,各点负压都有较大升高,但压力整体趋势较空气介质时变化不大,主要原因是α值较小时节流引导装置本身的节流阻力较大,当在接触腔中加入水时,沿倾斜流道下流的水流对气流的阻碍作用非常明显,极大的增加节流阻力,造成接触腔中负压整体较高;同时气流流量降低,造成对水的冲击力较小,激起水滴也较少,对接触腔中的流动状态影响较小㊂随α值不断增大,接触腔内压力受水介质影响先减弱后增强㊂当α=250/300和235/300时,各点负压向空气介质时靠近,并有个别点压力低于空气介质时的压力,主要原因是随着α值增大节流阻力变小,相比节流阻力降低水流对气流的阻碍作用对接触腔压力的影响减弱;同时气流流量升高,对水流的冲击作用增强,激起水滴增多,接触腔中有环量Г增大,局部压力(点4㊁5㊁6)降低幅度增强㊂当α=220/300时,2㊁3㊁4㊁8㊁9点压力进一步接近空气介质时,由此说明节流阻力减小对接触腔压力变化的影响削弱,而5㊁6㊁7点压力进一步降低,说明水流对接触腔压力影响增强㊂另外由于接触腔中激起水滴较少,说明气流对倾斜板上流动的水流冲击时动能传递较少,而更多的是水流对气流的变向作用,从而使5点气流速度较大,Г较大㊂以上现象说明气-水耦合形式表现为:当α较大时,水主要增加节流阻力来影响接触腔中的压力场;当α较小时,水流只起改变气流流向作用来影响接触腔压力场;当α在以上两者之间时,水流主要以改变气流流向和激起水滴来影响流场㊂因此只有当α值取合适时,才能保证气-水耦合良好,达到较好的除尘效果,即产生相当数量的水滴又一定程度上改变气流流向㊂4 结 论(1)得出了接触腔的高度直接影响除尘器气-水两相耦合程度进而影响除尘器除尘效果,高度过高水滴较多分布在接触腔下部,气-水耦合程度低,除尘效率不高,接触腔最佳高度为300mm㊂(2)随着节流强度α值增大接触腔中压力整体下降,接触腔压力沿流线方向逐渐降低,压力与节流强度近似呈线性关系,接触腔结构的不连续过渡和结构死区将产生漩涡,漩涡强度随α值增大而增大㊂(3)随α值变化气-水耦合表现为4种形式:增加节流阻力㊁激发除尘水滴㊁改变气流流向与3种形式共同作用,对于400mm×300mm 的除尘器,当α值在250/300~235/300时为3种形式共同作用,即产生相当数量高速运动的水滴,又一定程度上改变气流流向,达到较好的除尘效果㊂参考文献:[1] 刘 建,姚海飞,魏传光,等.掘进工作面湿式离心除尘器的结5141煤 炭 学 报2011年第36卷构优化及数值模拟[J].煤炭学报,2010,35(3):424-428.Liu Jian,Yao Haifei,Wei Chuanguang,et al.The structural optimiza⁃tion and numerical simulation of wet centrifugal precipitator at work⁃ing face of diving airway[J].Journal of China Coal Society,2010,35(3):424-428.[2] 姚国芳,李 勇,王广胜.煤矿巷道掘进的综合降尘与防尘[J].煤炭技术,2009,28(12):95-97.Yao Guofang,Li Yong,Wang Guangsheng.Integrated dust fall and dust proof of roadway excavation in coal mine[J].Coal Technology, 2009,28(12):95-97.[3] Tasneem Abbasi,Abbasi S A.Dust explosions⁃cases,causes,conse⁃quences,and control[J].Journal of Hazardous Materials,2007,140: 7-44.[4] Petavratzi E,Kingman S,Lowndes I.Particulates from mining opera⁃tions:a review of sources,effects and regulations[J].Minerals Engi⁃neering,2005,18:1183-1199.[5] 甘正旺,许振良.洁净煤技术及其发展前景[J].辽宁工程技术大学学报,2005,24(4):253-255.Gan Zhengwang,Xu Zhenliang.Clean coal technology and develop⁃ment prospect[J].Journal of Liaoning Technical University,2005, 24(4):253-255.[6] 陈维民.掘进工作面自激式除尘器除尘机理的研究[J].中国矿业大学学报,1993,22(3):67-73.Chen Weimin.The mechanism study of dust capture of self⁃impinge⁃ment scrubber in heading face[J].Journal of China University of Mining&Technology,1993,22(3):67-73.[7] 蒋仲安,金龙哲,陈立武,等.掘进巷道粉尘控制技术的研究[J].中国安全科学学报,1999,9(1):11-15.Jiang Zhong’an,Jin Longzhe,Chen Liwu,et al.Study of dust control technology on tunnelling[J].China Safety Science Journal,1999,9(1):11-15.[8] Yu Yong,Zhou Lixing,Wang Baoguo,et al.A USM⁃Θtwo⁃phase tur⁃bulence model for simulating dense gas⁃particle flow[J].Acta Me⁃chanica Sinica,2005,21:228-234.[9] 陈 斌,傅宇晨,郭烈锦,等.水平管束间气液两相流局部含气率分布的实验研究[J].化工学报,2003,54(3):316-320.Chen Bin,Fu Yuchen,Guo Liejin,et al.Experimential investigationof distribution of void fraction between horizontal tube bundle[J].Journal of Chemical Industry and Engineering(China),2003,54(3):316-320.[10] 张金红,阎昌琪,孙中宁.摇摆状态下水平管内气液两相流流型转换研究[J].哈尔滨工程大学学报,2008,29(10):1050-1053.Zhang Jinhong,Yan Changqi,Sun Zhongning.Investigation of flowpattern transition for gas⁃liquid two⁃phase flow in horizontal rollingpipes[J].Journal of Harbin Engineering University,2008,29(10):1050-1053.[11] 孙自强,张宏建.气液两相流水平绕流柱体的动态压降特性[J].高校化学工程学报,2008,22(3):378-383.Sun Ziqiang,Zhang Hongjian.Dynamic pressure drop property ofgas⁃liquid two⁃phase flow around a horizontal cylinder[J].Journalof Chemical Engineering of Chinese Universities,2008,22(3):378-383.[12] 阎昌琪,于凯秋,栾 锋,等.摇摆对气-液两相流流型及空泡份额的影响[J].核动力工程,2008,29(4):35-38,49.Yan Changqi,Yu Kaiqiu,Luan Feng,et al.Rolling effects on two⁃phase flow pattern and void fraction[J].Nuclear Power Engineer⁃ing,2008,29(4):35-38,49.[13] Karl Bühler.Pattern formation in rotating fluids[J].Journal ofThermal Science,2009(2):109-118.[14] 梁法春,王 栋,林宗虎.基于管壁取样的气液两相流量测量[J].西安交通大学学报,2008,42(1):52-55.Liang Fachun,Wang Dong,Lin Zonghu.Gas⁃liquid two⁃phase flowmetering with pipe wall sampler[J].Journal of Xi’an Jiaotong Uni⁃versity,2008,42(1):52-55.[15] You Changfu,Qi Haiying,Xu Xuchang.Lift force on rotating sphereat low reynolds numbers and high rotational speeds[J].Acta Me⁃chanica Sinica,2003,19(4):300-307.[16] 李小川,胡亚非,吴晓冰.石墨机加工车间的粉尘处理[J].中国科技论文在线,2010(12):952-956.Li Xiaochuan,Hu Yafei,Wu Xiaobing.Dust removing for lathingprocess of graphite⁃produce[J].Sciencepaper Online,2010(12):952-956.6141。