基于格子Boltzmann方程的大涡模拟对湍流时空关联性的研究

将格子Boltzmann方程和大涡模拟(LBE-LES)相结合,提出适应于格子Boltzmann方法(LBM)的涡黏性亚格子尺度模型,开展均匀各向同性湍流时空关联性的研究.采用D3Q19格式计算湍流的三维能谱、湍动能耗散率和其它高阶统计量,与实验和直接数值模拟结果的比较表明,该模型比传统涡黏模型有明显改进.考察了不同亚格子模型预测湍流频率波数能量谱的能力,结果表明,尺度涡产生的横扫作用是造成小尺度涡时间去关联的主要因素,不同波数的频率能量谱之间有一定的相似性,横扫速度是描述湍流频率波数能量谱的特征量.     

两种亚网格湍流模型的旋流扩散火焰大涡模拟

为研究旋流火焰结构,采用二阶矩亚网格(SOM-SGS)燃烧模型及Smagorinsky-Lilly和K方程亚网格湍流模型,对美国Sandia国家实验室测量的旋流火焰进行了大涡模拟,得到了与燃烧场速度、温度和温度脉动实测值相吻合的模拟结果。预报的瞬时温度分布云图与实际火焰的形状很相似。K方程亚网格湍流模型预报的瞬时温度分布比Smagorinsky-Lilly模型的预报结果更接近实际。燃烧火焰基本上位于回流区所在的位置,火焰在回流区被稳定。     

两种亚网格湍流模型的旋流扩散火焰大涡模拟

为研究旋流火焰结构,采用二阶矩亚网格(SOM-SGS)燃烧模型及Smagorinsky-Lilly和K方程亚网格湍流模型,对美国Sandia国家实验室测量的旋流火焰进行了大涡模拟,得到了与燃烧场速度、温度和温度脉动实测值相吻合的模拟结果。预报的瞬时温度分布云图与实际火焰的形状很相似。K方程亚网格湍流模型预报的瞬时温度分布比Smagorinsky-Lilly模型的预报结果更接近实际。燃烧火焰基本上位于回流区所在的位置,火焰在回流区被稳定。     

方柱绕流的大涡模拟

采用大涡模拟(LES)方法对雷诺数为2.2×104的方柱绕流流场进行了数值模拟.使用非交错网格的有限差分法,分别对准三维物理模型和真实三维物理模型求解不可压N-S方程,将沿流向方向方柱水平中心线上的时均速度的计算结果与实验数据进行了比较,结果表明,三维模型的模拟结果优于准三维模型的模拟结果.比较升力系数和阻力系数发现,与二维模拟(RNG)方法相比, 三维模拟的结果更加接近实验测试数据.     

制动工况下液力变矩器大涡模拟流场仿真及可视化试验验证

为了准确地揭示液力变矩器流场时空演化机理,基于计算流体动力学理论,采用5种不同的大涡模拟(LES)亚格子湍流模型(SL、WALE、WMLES、WMLES S-Omega、KET)精细仿真制动工况下液力变矩器三维流场,识别并提取涡轮内部非定常多尺度三维涡结构,籍此分析时空多尺度涡演化特征及其对流场结构演变的影响规律.基于...     

城市大气微环境大涡模拟研究

良好的空气质量是居民生活环境的重要条件之一,大气微环境研究与居民区的空气质量密切相关．大气环境由多尺度的空气运动主宰,而城市居民区空气流动属于微尺度大气流动．介绍城市居民区微尺度大气环境的现代数值方法——大涡数值模拟方法．针对居民小区复杂流动的特点,讨论如何正确应用大涡数值模拟方法,比较不同数值方法对预测城市微环境流动的影响,如：下垫面阻力元法与浸没边界法,亚格子模式的影响等．以模型小区和实际小区为例计算了风场（包括湍流统计特性）,污染物浓度等物理量分布,数值计算结果和模型小区及实际小区的实验测量数据符合良好．     

180°方形弯管中湍流流动的大涡模拟

把大涡模拟应用于复杂几何形状的流域是一个重要研究课题。对１８０°方形弯管中的湍流流动作大涡模拟,以网格节点为中心的交错网格作为离散网格,用有限体积方法和中心差分进行空间离散,用Ａｄａｍｓ－Ｂａｓｈｆｏｒｔｈ格式进行时间离散,用简单的Ｓｍａｇｏｒｉｎｓｋｙ模型来作为次网格湍流模型。用多重网格方法求解压力方程以提高计算效率。数值模拟结果与实验结果取得了很好的一致。     

大涡模拟在水利中的研究进展

大涡模拟作为一种研究湍流流动的有效手段，在国际上已经得到了广泛的应用．本文在介绍了大涡模拟的发展历程及其理论方法的基础上，就国内外大涡模拟在水利中的最新研究成果作了较为详细的阐述和系统的分析，指出了其中应该注意的问题和今后的发展方向。     

基于改进大涡模拟模型的风荷载验证

为了获得与近地大气边界层实际脉动风场特征相一致的三维脉动风速,文章首先应用DSRFG(Discretizing and Synthesizing of Random Flow Generation)方法,模拟生成计算流体动力学中大涡模拟算法所需的模拟入口边界脉动风速时程,然后采用改进的动态一方程模型作为大涡模拟的亚格子模型,计算长宽高为1∶1∶3的高层单体矩形建筑在0度风向角下的风荷载分布及特性,并与日本TPU风洞试验进行了对比.结果表明:1当高层建筑长宽比(B/D)=1时,扭矩的能量主要来自涡旋脱落而不是入口湍流的影响;2DSRFG结合改进的动态一方程模型能够较为准确地模拟得到高层建筑表面的平均及脉动风荷载特征.此方法有望应用到实际高层建筑风荷载评估中.     

二维大涡数值模拟的算例及其分析

根据文[1]中给出的二维大涡模拟方程的离散形式及数值解法,编制了一套数值计算程序,并通过计算二维后阶湍流流动和双流道式污水泵叶轮内部流场的实例计算,验证了文[1]中推导出的方程的正确性、可解性,以及提供的数值解法的可行性和可靠性.经文中的实例计算,并把计算结果与实验实测结果相比较,令人十分满意.从而表明:二维大涡数值模拟既保持了原有三维形式的优点,又大大简化了三维求解对计算条件的苛刻要求,有望使大涡数值实现广泛的工程应用.     

离心泵叶轮内部清水湍流的动态大涡模拟

提出旋转坐标下的单相液体湍流的二阶双系数动态亚格子应力大涡模拟模型并加以验证。考虑亚格子应力的对称性和量纲一致性 ,在基于应变率张量和旋转率张量的不可约量及 Smagorinsky模型和亥姆霍兹定理的基础上 ,提出亚格子应力应表示为可解的应变率张量和旋转率张量的函数。在贴体坐标系下 ,利用交错网格系统和有限体积法离散湍流控制方程 ,采用 SIMPL EC方法求解方程。水泵叶轮中湍流的流速分布规律和压力分布等计算结果都与试验结果基本吻合 ,从而证明了此模型的可行性 ,以及计算方法和程序的可靠性     

二维大涡数值模拟方程及其解法

采用高斯型滤波函数对Navier Stokes方程进行滤波处理 ,并引入亚格子雷诺应力模型推导出了在水力机械内部流场计算中实用的二维大涡模拟方程及其在物理平面和计算平面上的离散形式 ,并介绍了方程的数值解法     

矩形射流流动的大涡模拟

为探讨矩形射流中的拟序运动特征,用大涡模拟方法(LES)对三维矩形射流流动过程进行了数值模拟。统计平均的流向速度横向分布曲线在发展区某些位置上呈现"马鞍型"分布,速度峰值偏离中心轴线,与圆射流显著不同,这与已有的实验和直接数值模拟(DNS)结果吻合良好。计算结果表明,在发展区湍流脉动最强,且应力的横向分布不均匀。随着流向距离的增大,速度功率谱(PSD)曲线呈现峰值向低频移动的现象。速度的PSD分布表明:矩形射流中拟序运动主要位于射流的核心区末端和发展区,射流发展区中的拟序运动具有明显的时间准周期性。     

基于大涡模拟的液滴蒸发液相传热过程

建立了燃油液滴蒸发瞬态气液两相流大涡模拟流体体积(VOF)模型,综合考虑了气液两相流流场,气液相界面的温度梯度和浓度梯度等因素,研究了气液相温度梯度和雷诺数对液滴液相传热的影响.研究结果表明:随着温度梯度的增大,液滴温度显著提高,温度梯度越大,液滴温度增加幅度越大;雷诺数不仅能够影响液滴周围流场的分布,而且还引起液滴内部温度场分布的变化;随着雷诺数的增大,液滴内部加热区域逐渐均匀,雷诺数越大,液滴的温度越高.     

液力缓速器瞬态两相流动大涡模拟及性能预测

为了全面掌握液力缓速器各个流动单元内速度场和压力场的分布特性,提取了全流道几何模型作为计算区域,利用CFD平台的大涡模拟法和多流动区域耦合计算的滑动网格法对液力缓速器内部气液两相流动进行了三维瞬态数值模拟,将混合模型与欧拉模型交替运用在其多相流模型中,获得了不同充液率下液力缓速器内流场结构的变化及两相体积分数分布情况,分析了流场内二次流、脱流及涡旋的产生机理,并计算了缓速器的外特性.结果表明:数值计算结果与试验结果吻合很好,误差在8％以内;流场计算十分准确,运用的大涡模拟方法可以有效地模拟液力缓速器流场内的真实液流结构,其结果可以用来指导液力缓速器的设计及其结构优化.     

基于动力拟序涡模型的水波三维大涡模拟

利用大涡模拟理论数值模拟了规则波和不规则波的传播过程,采用动力拟序涡黏性亚格子应力模型进行波高及速度场的模拟.计算结果表明,模型对水波的波高及由波产生的流体速度分布模拟结果与解析值符合良好,且精度较高.对比分析了人工海绵层及Sommerfeld辐射型耦合边界条件与仅采用Sommerfeld辐射型边界条件对水波数值模拟精度的影响.结果表明动力拟序涡粘性模型可以准确地模拟水波传播特性.     

液力缓速器部分充液流场大涡模拟及特性预测

利用CFD平台的大涡模拟法和滑动网格法对液力缓速器内部气液两相流动进行三维瞬态数值模拟,将混合模型与欧拉模型交替运用在其多相流模型中,得到速度场和压力场的分布特性,总结不同充液率下液力缓速器内流场结构的变化及两相体积分布情况.基于流场数值解计算了液力缓速器不同工况下的制动转矩,进而预测其性能,性能预测结果与试验结果吻合较好,误差在8％以内,表明流场计算是较准确的,其结果可用于液力缓速器的设计及其结构优化.     

平面自由湍射流拟序结构的大涡模拟研究

对空间发展的平面不可压缩湍射流拟序结构的非定常演化过程进行了大涡模拟。采用标志物浓度等值线分布示踪平面射流气相流动 ,Reynolds数为 1130 0 ,模拟结果再现了平面射流中 Kelvin- Helm holtz不稳定性的发生以及展向大尺度涡的卷起、合并、破碎过程。捕获到了射流拟序结构剧烈相互作用的“偶极子”和“三极子”现象。在拟序结构的配对过程中 ,大尺度旋涡的尾迹混合形成了流向“发卡”型拟序涡结构。平面射流拟序结构的时空演化过程所呈现出的卡门涡街特征主要集中于初始段的后部和过渡段的前端。数值流场显示结果和染线法的实验结果吻合得很好