基于DNS和POD方法的粘弹性湍流减阻流动特征研究

来源 :中国力学大会2011暨钱学森诞辰100周年纪念大会 | 被引量 : 0次 | 上传用户:tnzx911
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  本文将直接数值模拟(DNS)方法和本征正交分解(POD)技术结合起来研究不同弹性减阻流体流动的平均特性和脉动特征结构,提取了速度和附加变形(与减阻剂带来的附加应力有关)的脉动特征结构,对含不同能量的典型湍流结构和附加变形结构进行了分析。结果显示减阻剂主要对大尺度含能结构产生作用,使得大尺度结构的能量比重增大,小尺度结构的能量比重减小,使减阻流动比牛顿流体流动表现出更强的各向异性。随着流体弹性的增强,减阻流体的主流脉动大尺度结构在流向伸长,在壁面方向梯度减小,交替出现了一系列顺时针和反时针流向大涡,流动的有序性得以增强;脉动流线由无规则缠绕逐渐变为沿流向和垂向的大尺度旋转和沿展向的对称分布;脉动涡结构增大,在流向被拉长,在垂向移向槽道内部,在展向被拉宽;脉动变形结构由壁面向流动区域内部不断迁移整合。这些说明减阻流体弹性的增强和减阻率的提高与大尺度螺旋运动和内部附加变形网络结构的增大有关,此时脉动结构分布更加有序,能量缓冲作用更大。
其他文献
Bypass转捩是工程中常见的现象。目前对其过程的描述是,自由流中的低频扰动进入边界层并激发条带结构,在条带结构下出现不稳定的高频扰动,这些扰动增长并导致湍斑出现,最后湍斑合并并达到湍流状态。但是这对breakdown过程的描述还不够详细。由于条带结构下的平均流剖面存在拐点,其不稳定性远大于Blasius剖面,因而导致Bypass转捩breakdown过程的关键因素可能和平均流剖面的稳定性有关。本
基于有限体积法采用雷诺应力湍流模型并结合Cross本构方程数值模拟高分子对Reyleigh-Bénard对流在软(Ra=5×105~4×107)、硬(Ra>4×107)湍流区热输运的影响。模拟结果显示:添加高分子可加强湍流Reyleigh-Bénard对流热输运,相对于软湍流区在硬湍流区的有效作用范围较大;相应地大尺度环流也随之充满整个腔体,封闭腔体内对角线上两个对称小涡也随着高分子浓度和Ra数的
线性渐近稳定流场是指一类在受到小扰动后经过足够长的时间能够消除干扰影响恢复原稳定状态的流动现象。最近,越来越多的研究人员开始对线性渐近稳定流场受到干扰后短时间内的流场演化过程展开研究。本文采用直接数值模拟方法计算低雷诺数(Re=100)下二维周期性旋转圆柱绕流流场,并将该非定常结果作为基本流场进行稳定性分析,发现在给定的转速范围内,旋转圆柱扰流的尾迹演化是一个线性渐近稳定或弱不稳定的过程。研究结果
本文采用Mortar谱元法和多处理器并行计算技术模拟了Kelvin-Helmhtz界面不稳定性湍流的混合发展过程,通过对混合层动量厚度、能谱和总动能的计算,评估了Kelvin-Helmhtz混合层的演化机理。计算结果表明:三维Mortar谱元法可以有效混合层流动的湍流混合和演化,能够捕捉到涡的合并现象和大涡到小涡的级联过程。
本文采用数值模拟的方法研究了展向周期变化槽道流的全局稳定性。首先用线性稳定性求得的扰动特征函数作为初始值,计算扰动随时间的演化。经过充分的发展,得到了全局稳定性的特征解。给出了不同b下的计算结果以及全局稳定性为中性时,Re与展向不均匀程度的关系,及其对应的特征函数、x-z平面上扰动的等值线图以及全局稳定性下的中性曲线。
Burnett方程是Boltzmann方程的二阶近似,能够描述轻微偏离热力学平衡时的稀薄气体流动。但是Burnett方程在小扰动下不稳定,这是限制Burnett方程广泛应用的一个重要原因。本文通过线性小扰动理论,首次得到了二维常规Burnett方程的稳定性特征方程,并把稳定性方程的解表示在复平面上,得到了二维稳定性特征曲线。通过扰动增长系数和扰动波数的关系,得到二维常规Burnett方程的临界Kn
结合理论分析和数值模拟对流体-微粒两相湍流的相互作用和减阻机制进行研究,采用点-力耦合模型的负载颗粒槽道湍流的直接数值模拟,结合拉格朗日轨道模型,通过模拟两相间的双向耦合刻画其相互作用的动力学过程,获得各阶次的统计数据,分析了微米级颗粒在壁湍流中的运动以及对湍流所产生的调制作用。同时,从含颗粒作用外力项的Navier-Stokes方程出发,推导颗粒加载湍流的阻力系数公式,研究颗粒的不同物理参数如直
应用热线测速技术,对沟槽壁面平板湍流边界层的减阻机理进行了实验测量研究。用热线风速仪和X形双丝热线探针测量了风洞中沟槽壁面及光滑壁面平板湍流边界层在不同法向位置的瞬时流向u、展向速度分量w的时间序列信号。用流向速度分量时间序列信号的多尺度子波系数辨识壁湍流多尺度相干结构,用条件采样和相位平均技术提取了壁湍流多尺度相干结构喷射和扫掠时流向速度分量u、展向速度分量w和瞬时雷诺应力分量uw的条件相位平均
采用充分发展的槽道湍流直接数值模拟数据库,系统地研究了壁面可测信号(流向脉动摩擦应力τwx、展向脉动摩擦应力τwz和脉动压力pw)和近壁区流向涡的内在关系。发现τwx与下游流向涡密切相关,而τwz则反映了正上方流向涡对壁面的影响。pw和近壁区流向涡的关系较为复杂:在探测点上游,高压区(pw>0)对应着流体的下扫,低压区(pw<0)则对应着流体的上抛;在探测点,近壁区流向涡和正下方低压区存在强相关性
在25℃,1~300s-1剪切速率范围内,通过对70、100、200、300和500 ppm CTAC/NaSal表面活性剂溶液在瞬态模式下的流变测试,研究了表面活性剂稀溶液的流变性质。研究表明,表面活性剂稀溶液在剪切作用下会形成剪切诱导结构(SIS),导致了溶液表观粘度的增大,且SIS开始形成和完全形成所需时间均与剪切速率成负相关,符合幂指数为-1.0的幂函数关系。表面活性剂稀溶液中,各向异性转