组合隔板抑制液舱晃荡机理的数值方法研究.docx

该【组合隔板抑制液舱晃荡机理的数值方法研究 】是由【科技星球】上传分享，文档一共【18】页，该文档可以免费在线阅读，需要了解更多关于【组合隔板抑制液舱晃荡机理的数值方法研究 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。组合隔板抑制液舱晃荡机理的数值方法研究*??罗鑫,董胜(中国海洋大学工程学院,山东青岛266100)液体在晃荡过程中会周期性的冲击舱壁,特别是当外部激励频率接近液舱一阶自然频率时,会产生极端冲击压力,威胁液舱结构的安全。有效的抑制手段成为液舱晃荡研究领域的热点。由于液舱通常是随运输工具的生产而制造安装的,考虑在舱内加入构件而抑制晃荡的方法最为实用。安装内嵌隔板是目前公认的最为有效的方法,一方面,隔板的加入可以在一定程度上改变液舱固有频率,使之远离外部激励频率;另一方面,隔板的存在可以产生阻尼,将随载体运动而固定输入的能量耗散掉,达到稳定液面、减小舱壁冲击压力的目的。但是,隔板的加入不仅会增加舱体自重、缩小载液容积,还存在受流体冲击而损坏的风险,并且这个风险会随隔板的体积增大而升高。因此,将内嵌隔板的总体积控制在一定范围内,并使其发挥最佳抑制晃荡效果至关重要。内嵌隔板的有效设计不仅依赖于对流体动力特性的详细了解,还需要理解隔板抑制晃荡的机理,并有效运用。隔板的布置形式需要考虑很多因素,如安装位置、形状、尺寸,数量、开孔处理以及组合方式,同样,外部激励属性如振幅、频率以及运动方式和液舱固有属性如载液率、几何形状、尺寸以及液体特性等也应被考虑在设计工作之内[1]。目前,内嵌式隔板通常为水平隔板、垂直隔板或二者的有机结合,国内外学者通过理论解析法、模型实验法以及数值分析法对不同隔板的抑制晃荡效果进行了大量的研究[2]。解析法可以从理论高度为模型试验和数值模拟提供指导。Abramson[3]最先开始研究液舱内嵌隔板对舱壁冲击载荷的消减问题;Evans和Mclver[4]用解析法研究了矩形液舱中垂直隔板对其共振频率的影响,发现靠近自由液面处的隔板可显著改变共振频率,而底部安装的隔板作用不明显;Maleki和Ziyaeifar[5]通过求解拉普拉斯微分方程给出了液舱中水平隔板和垂直隔板抑制晃荡的理论阻尼模型并指出环形隔板抑制晃荡的效果更好;Faltinsen和Timokha[6]提出了一个精确的理论近似解用于模拟带有垂直隔板的矩形液舱,发现舱内的垂直隔板会显著的改变液舱共振频率。在模型试验方面,Akyildiz和Unal[7]进行了矩形液舱内晃荡引起的压力分布的模型试验,并设置了不同载液率和隔板布置形式的工况,发现在舱底中心处安装垂直隔板可有效削弱舱壁所受压力,其原因在于流体经过垂直隔板时可在附近产生剪切层并耗散能量;Panigrahy等[8]同样考虑了安装隔板和无隔板情况下对舱壁压力和液面爬高的影响并进行了一系列试验,并得出环形隔板比普通水平隔板更有利于减小作用于舱壁上的压力;Sauret等[9]通过试验研究了晃荡过程中气泡的耗能作用并得出仅舱壁附近的气泡可有效耗能的结论;Nayak和Biswal[10]进行了多种构件抑制晃荡的效果对比试验,试验结果显示自由液面处悬挂式隔板能最大程度地减小晃荡幅度。数值模拟可以解决理论上无法考虑的复杂情况,还可以低成本的处理模型试验中难以解决的问题。Cho和Lee[11]使用数值方法研究了二维液舱水平受迫运动下隔板对液面爬高的抑制作用,并指出布置内嵌隔板对于抑制晃荡的重要性;基于线性化势流理论,Cho等[12]通过添加自由液面下工阻尼,给出了横向受迫运动下二维带隔板液舱内共振特性的数值分析,对比了不同载液深度、隔板数量及安装位置对抑制晃荡效果的影响;Delorme等[13]使用SPH法和模型试验研究了浅水晃荡的冲击压力载荷,并发现极端冲击不是发生于外部激励频率等于一阶自然频率处,而是稍大值处;Liu和Lin[14]通过求解Navier-Stolces(N-S)方程研究了内嵌隔板对于液舱晃荡的影响,并得出垂直隔板抑制晃荡的效果优于水平隔板的结论;Biswal和Bhattacharyya则考量了流体与舱内弹性隔板之间的耦合关系,发现流体和结构响应可以通过内嵌隔板有效控制[15];Sygulski[16]从数值角度分析了三维带有任意隔板的液舱的自然频率和模态;Akyildiz[17]用数值模拟了矩形液舱中垂直隔板高度对液舱晃荡的影响,研究发现垂直隔板抑制晃荡的机理主要源于其水动力效应而与其高度关系不大;Wu等[18];Hasheminejad和Mohammadi[19]基于线性势流理论数值模拟了圆形液舱中隔板抑制晃荡的效果,结果显示隔板安装在不同位置会改变液舱的晃荡频率及动压力场;Ebrahimian[20]发现隔板的安装位置比隔板尺寸更能影响晃荡的自然频率。上述研究表明,单一隔板抑制晃荡的效果有限,多隔板组合的方式可能会产生更好的抑制晃荡作用,而布置方式是组合隔板抑制晃荡效果的重要影响因素。为此,必须深入了解垂直隔板与水平隔板的抑制晃荡机理,并分析其组合效应。鉴于此,本文采用具有高阶精度的ConstrainedInterpolationPorfile法(简称CIP法)建立数值模型,并以基于VOF法优化格式捕捉自由液面,通过数值计算,分析了单一隔板的抑制晃荡机理,并且对比了不同隔板组合方式抑制晃荡的效果。、动量守恒和能量守恒定律,本研究以气体与水体为介质的二维、两相不可压缩的黏性流动问题,忽略温度变化的影响,控制方程如下:(1)(2)忽略结构的弹性,引入运动坐标系以简化壁面边界条件,所有网格点均在随体运动坐标系下进行,将计算区域内的固体视为密度无限大的流体,转化为多相流问题。固体边界条件定义为界面无滑移。为捕捉锐利的自由液面,在每个计算网格上引入体积函数φm(m=1,2,3),以VOF法的思想将计算区域内的不同介质区分开,体积函数满足如下方程:(3)式中:m=1表示网格为液体,m=2表示网格为气液混合,m=3表示网格为气体;Ωm表示被水体或空气占据的计算区域;φm表示水或气在控制单元中所占比值。式(3)的最左端为全微分,其物理意义为随着流体质点运动,该质点的φm不随时间变化。该方法将气体和水体作为一种流体模拟,自动满足界面处切应力动力条件,不会使水气交界面发生突变。自由面可作为一种内交界面通过运动界面追踪方法计算得到。,最早由Takewaki等[21]为求解双曲型偏微分方程而提出。其特征是使用三次内插值进行N-S方程对流项的离散求解,最初被称为CubicInterpolatedPseudo-Particle法。其原理是用已知时间层的函数值及其导数,通过三次内插值近似,进而时间推进式的求解下一时间层网格单元内部变量的信息,简洁地构造出三阶精度差分。CIP法用以离散求解N-S方程的对流项,故下面以一维对流方程为例,简要说明CIP法的原理。(4)式中:f为函数曲线;传播速度u为常数;x为矩形波传播方向。通常所采用的一阶迎风差分格式由两个相邻的网格点构造而成,线性插值的使用会导致网格内的剖面信息丢失,引起较大的数值耗散;而为了避免此问题而使用的高阶迎风差分格式需要使用多个网格点,但此时不仅会增加运算量,更会造成峰突。而CIP法采用一种独特的三阶内插值方式,利用变量及其空间导数,只使用两个网格点构造出高阶差分格式,更为准确的描述网格内信息。对式(4)求空间导数可得(5)式中:g=?f/?x。为了简化起见,将对流速度u设为常数,式右边项为0。这样,两个方程就有了相同的形式。不妨假设u>0,在迎风向网格单元[xi-1,xi]内n时刻剖面函数可以近似为(6)式中:ai、bi、ci、di分别为待定系数。在n+1时刻的单元剖面函数fn+1可以通过将n时刻的剖面函数fn平移-uΔt得到,函数f和g的时间演变可以通过下面的拉格朗日变化得到(7)(8)(9)CIP法基于两个节点信息使用三次插值函数得到未知网格点信息,格式紧凑,与其他传统的高阶迎风格式相比,使用较少节点从而提高解算效率,同时具备高阶精度。求解N-S方程采用分步算法对动量方程进行时间积分,其数值计算分为两个过程:对流项求解过程和非对流项求解过程,后者又分为两步:非对流步(I)处理方程的耗散及其他,非对流步(II)处理方程的压力耦合,得到最终流场信息。在液舱晃荡问题的研究中,如何捕捉自由液面是一个难点。本文引入体积函数,使用THINC格式捕捉自由液面。该格式是由Hirt[22]和Xiao[23]等基于VOF法的原理提出的,可精准捕捉到锐利的自由液面,可清晰地观察到液舱晃荡过程中冲击舱壁、液滴飞溅等现象。,这种高精度、低耗散的有限差分法十分适合用于液舱晃荡现象的模拟。Hu等[24]使用CIP法对无隔板液舱晃荡中产生的强烈液面变形进行了模拟;Zhu等[25]验证了CIP法在模拟流固耦合中处理大变形运动界面的优秀能力;Kishev[26]使用CIP法模拟了剧烈的液舱晃现象;Chen[27]将惯性坐标系与CIP法结合模拟了液舱晃荡现象,并以THINC格式捕捉到更为锐利的界面,通过与Xue[28]和Kishev的实验及数值结果对比,证明了CIP法能精确的模拟液舱晃荡现象,并进行了网格依赖性验证。本文将使用CIP法模拟带有不同内嵌隔板的液舱在晃荡过程中的水动力特征,通过对舱壁处的压力监测,分析不同组合形式的隔板抑制晃荡的效果。,其中,x表示容器长度方向尺寸,y表示液面至容器底的距离。液舱长L=,高H=,载液深度固定为h=。[29]研究表明,载液率h/~,舱内液体的晃荡最为剧烈,本文所选工况在此区间内,适合用于隔板抑制晃荡效果的研究。液舱横向受迫,外部激励可表示为η=A·sin(ωt)。(10)式中:η表示自由液面位置;A为振幅,;ω为外部激励频率。对于给定的液舱,其自然频率可按下式计算:(11)式中:n为阶数。当外激频率接近一阶自然频率ω1=,舱内液体发生共振,此时舱壁处发生最大冲击载荷,故以下工况均采取此外激频率。在自由液面与两侧壁交汇处设置测点,监测压力及液面高度变化。计算工况的振荡周期约为1s,计算域均布90000个网格点,空间步长dx=dy=2×10-3,时间步长dt=1×10-4。运动开始后,液面自静止状态起动先向左侧壁冲击,液体沿着舱壁不断爬升,在约1/4周期时达到最高点并开始回落,1/2周期时回复至相对平静状态,后1/2周期反向冲击,形态及过程与前者大致相同,在一个周期内对左右舱壁分别产生一个局部冲击压力的峰值,如此往复,会在某个周期中产生极大值。图2截取了上述过程中冲击两壁的瞬间,可清晰的捕捉到水舌、液滴飞溅等。。可以看出,在外激频率接近液舱一阶自然频率时,壁面出现较大动压载荷,冲击压力周期性的呈现双峰形态,首峰由高速流体冲击壁面产生,次峰由流体冲击后回落产生,呈现“教堂穹顶”形态。,根据现实条件从安装成本、性价比以及风险防控角度出发,内嵌构件的总体积是有限的,所以应使总体积固定的隔板通过组合的方式发挥最佳抑制效果。为达到此目的,先对单一隔板抑制晃荡的机理及特性进行研究。所有涉及到的隔板均采取统一尺寸,,。先考虑底部安装的垂直隔板对液舱晃荡的影响,将该工况编号为case1,垂直隔板安装于舱底中线处,相应的,将无隔板情况编号为case0。图4为case1隔板布置示意图,晃荡运动开始后,液面的周期性运动规律与case0大致相同,但由于隔板效应,流体流经时会产生涡旋耗散能量,达到抑制晃荡的目的。。可以看出液面开始冲击左侧壁的过程中,在垂直隔板左侧逐步产生了一个逆时针涡旋,在液面爬高达到峰值并开始回落时有收缩及向右传播的趋势,随后左侧的逆时针涡旋逐渐消融,隔板右侧开始凝聚一个顺时针涡旋,在经历了冲击右侧壁并回落后又向隔板左侧传播,如此往复,伴随着涡旋的不断产生,系统总能量大量耗散。对比case0与case1的两侧壁压力时历曲线。,压力峰值失去双峰特性,而双峰形态系动压载荷猛烈冲击舱壁所致。另外,晃荡过程中周期性呈现的压力值所对应频率也稍有改变,可能与内嵌隔板的调频作用有关,在今后的研究中将进一步找出其相关性。,对液舱晃荡过程中周期性的流体运动有了基本了解,观察到垂直隔板附近涡旋的产生与变化并对内嵌隔板的工作机理有了基本