论文部分内容阅读
将少量柔性长链高分子聚合物或某些表面活性剂加入水或者有机溶剂中可以明显减小湍流流动中的摩擦阻力,这种现象被称为添加剂湍流减阻效应。为了使这一现象更好地在实际工业系统中得到应用,亟需对其流动特性和减阻机理进行研究。由于实验方法不能获得粘弹性流体湍流减阻流动中的分子变形场和弹性应力场信息,因而大多采用数值模拟方法。大涡数值模拟(LES)因其在相同工况下计算量小于直接数值模拟(DNS)、可获得的信息量比雷诺平均数值模拟(RANS)多而备受关注。应用LES研究特定湍流流动时,关键在于所选用的亚格子模型是否适用于该种湍流流动,几乎所有已发表的亚格子模型只针对于牛顿流体湍流流动,而粘弹性流体湍流减阻流动LES研究刚刚起步,目前只有应用时间近似解卷积模型(TADM)尝试对粘弹性流体湍流减阻流动进行数值模拟。为了更好地应用LES研究粘弹性流体湍流减阻流动,需开发更多可靠、有效的亚格子模型。本文将从这一角度出发,结合粘弹性流体湍流减阻机理,构建针对粘弹性流体湍流减阻流动的新型亚格子模型,这对粘弹性流体湍流减阻流动LES研究有重要的理论意义和学术价值。从粘弹性流体湍流减阻机理的角度出发,基于反映粘弹性效应对湍流相干结构影响的思想,创新性地建立了耦合空间过滤和时间过滤的新型亚格子模型,命名为MCT (Mixed subgrid-scale model based on Coherent structures and Temporal approximate deconvolution)。选用粘弹性流体强迫各向同性湍流和槽道湍流作为验证对象,将低雷诺数(Re)下两种湍流流动的LES结果与对应的DNS结果进行比较,同时将较高Re下的LES槽道湍流结果与实验结果进行对比,结果表明本文所提出的MCT亚格子模型能够预测无壁面效应和有壁面效应的粘弹性流体湍流减阻流动,达到了本文的预期目的。在采用MCT和TADM亚格子模型数值模拟较高Re下的粘弹性流体强迫各向同性湍流时,发现TADM亚格子模型表现出过度耗散,而MCT亚格子模型并无此现象,表明MCT亚格子模型在计算较高Re下的湍流减阻流动方面具有优越性,这为LES研究高Re下粘弹性流体湍流减阻流动奠定了基础。采用MCT亚格子模型建立了较高Re下的粘弹性流体强迫各向同性湍流的LES数据库,并分析了其流动特性。结果表明粘弹性的存在明显改变了湍流的流动特性,涡量和拟涡能明显减小,小尺度涡结构数目明显减少,随着粘弹性效应的增强减阻率增大且粘弹性流体分子拉伸更为剧烈。同时,采用一维和二维小波变换进一步探究粘弹性效应对湍流多尺度特性的影响,发现湍流流动中的间歇性随着尺度的减小而增强,在粘弹性流体中不仅整个流动区域的涡结构以及间歇性受到抑制,而且局部涡结构对间歇性的贡献减弱。采用MCT亚格子模型对较高Re下的粘弹性流体槽道湍流进行数值模拟,建立了较高Re下的LES数据库,并结合槽道湍流实验数据库分析其流动特性及减阻机理。结果表明在粘弹性流体中平均速度场、湍流脉动强度、雷诺剪切应力、低速条带结构等均发生变化,且粘弹性流体中的过渡层有向主流区扩展的趋势。根据湍流流动阻力贡献、湍动能贡献、基于小波变换的标度律以及间歇性的结果,很好地解释了粘弹性流体湍流减阻机理。从整体角度出发,粘弹性流体中湍流贡献的减小程度比引入的粘弹性贡献大,从而使湍流流动的摩擦系数减小,进而产生减阻效果;从粘弹性效应对壁面法向不同区域的影响角度出发,粘弹性效应主要抑制了过渡层中从线性底层由湍流猝发事件引起的上抛运动而来的相干结构。