大涡模拟二维平面尾迹湍流拟序结构

首先采用大涡模拟方法对二维平面尾迹流进行了数值模拟,分别分析了整个流场的流动特点和近壁处展向大涡结构的形成、发展、脱落和破碎过程.然后定量计算了流场下游不同断面上流向时均速度分布、流向时均湍流强度分布以及流向相均速度的分布情况,并与Yang等最近所做的实验进行了对比分析.数值结果表明大涡模拟方法可以较好地模拟二维平面尾迹流.     

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

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

平面射流拟序结构的离散涡模拟

采用离散涡方法对平面射流流场进行了数值模拟。显示了平面射流中的大尺度涡结构的卷起、旋转和合并过程,模拟结果揭示了平面射流中的大尺度拟序结构的时空演化特征和发展变化规律。得到了定点的瞬时速度和流场的时均速度分布,时均速度场的计算结果和实验值吻合,瞬时速度分布揭示了平面射流拟序结构的非定常、不稳定的非线性特征。     

湍流横向射流大涡模拟及大尺度结构的涡环模型

采用大涡模拟方法数值模拟了流向椭圆喷嘴的湍流横向射流,重点研究其旋涡结构的产生、发展等动力学演化过程．结果表明：文献中所报道的横向射流基本涡结构,如反向旋转涡对、前缘涡、后缘涡、悬涡、肾涡、反肾涡等并非独立的物理实体,而是分别对应于新发现的横向射流基本涡结构——起始于喷嘴的三维拉伸涡环的局部结构．     

湍流横向射流的大涡模拟及其涡结构特性

采用大涡模拟方法数值模拟了流向和展向椭圆喷嘴的湍流横向射流, 重点研究了其中旋涡结构的产生、发展等动力学演化过程. 结果表明文献中所报道的横向射流基本涡结构, 如反向旋转涡对、前缘涡、后缘涡、悬涡、肾涡、反肾涡等等是分别对应于新发现的横向射流中的基本涡结构——起始于喷嘴的三维拉伸涡环的局部结构, 因此, 在湍流横向射流中真正占主导作用的是拉伸、扭曲、沿展向摆动和沿流向扭动的三维涡环. 研究还发现: 涡环的脱落频率比流场信号分析得到的脉动频率小得多.     

矩形射流流动的大涡模拟

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

基于尺度相似假设的大涡模拟动力方法

当前的大涡模拟动力方法对G erm ano等式中的两个尺度的亚格子应力均使用相同的亚格子基模型。这样为求得基模型系数,必须假定该系数在测试滤波操作下保持不变。但实际上该系数存在剧烈的脉动,与假定并不一致。该文提出的方法对以上两个应力分别应用亚格子基模型和尺度相似假设。这样就无需对系数分布做出假定并保证了方法的内在一致性。用该方法对R e为104的方腔流动进行了大涡模拟。结果表明该方法对流动统计量的预测有所改进,并可有效消除亚格子模型系数的奇异性。     

基于大涡模拟的螺旋桨梢涡数值分析

为研究螺旋桨的梢涡脱落规律,基于STAR-CCM+软件平台,采用重叠网格技术和大涡模拟方法,通过求解纳维-斯托克斯方程,对螺旋桨的速度场和涡量场进行了详细的研究,并且比较了不同进速系数下螺旋桨梢涡形态的差别.数值研究结果表明:基于正交网格的大涡模拟能够较好地模拟螺旋桨的梢涡形态,低进速系数下螺旋桨的梢涡强度增大并且显得比较紊乱;高进速系数下螺旋桨梢涡则呈现出比较规则的螺旋涡管形状.     

近床面风沙流的颗粒拟流体大涡模拟分析

通过气相大涡模拟与基于颗粒动理学理论的颗粒相湍流拟流体模型进行耦合,对近床面风沙气固两相三维运动进行了数值计算。相同来流风速下,对跃移层内沙相速度、沙相浓度及沙粒脉动均方根速度沿近床面高度变化进行了统计分析,对比了模拟参数结果与相关实验值之间的差异。结果表明:所建计算模型能较好地揭示近床面风沙流特性,低沙相浓度运动模拟中2种曵力系数模型计算得到的各个宏观参数统计量趋于一致;跃移层内沙相平均流向速度随近床面高度的增加呈递增趋势,沙相浓度随近床面高度的增加呈递减趋势,沙相流向脉动均方根速度在跃移底层呈现单一峰值状。模拟结果获得了与相关实验结论相类似的变化规律,得出准确选取沙粒碰撞弹性恢复系数是准确揭示沙相速度脉动强弱的关键。     

可压缩混合层流动的三维大涡模拟

为了加深对复杂流动混合现象的认识,为化工设备设计提供理论基础,用大涡模拟方法对可压缩气体混合层进行了三维数值模拟.计算结果显示了混合层流动的三维瞬态发展过程以及流动过程中涡的产生、生长、破裂以及涡列的形成与各涡的合并过程,揭示了可压混合层三维演变过程中同样具有自相似性,计算结果与实验结果相符.通过对混合层初始状态条件下不同扰动和不同初始速度比的计算表明:随机扰动量对混合层的失稳演变过程影响最大,而大混合层速度比同样可使混合层失稳加剧并提前.     

方柱绕流的大涡模拟

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

大涡模拟预报可靠性分析

针对明渠道湍流流动,研究了计算网格、计算域大小、亚格子（SGS）模型和雷诺数等因素对大涡模拟（LES）预报精度的影响。研究表明,即使简单外形的湍流流动,网格尺寸大小及匹配仍较大地影响大涡模拟计算结果,因此具体使用中应予以充分重视。     

获得空间发展槽道流动大涡模拟入口条件的最不稳定波方法

在空间发展流动的大涡模拟中,构造有效的入口扰动是关键。该文提出了利用Fourier分析获得流动最不稳定波的方法来确定入口扰动。结果表明,与同均方差的G auss分布扰动相比,最不稳定波扰动能够显著缩短流动的起涡距离。同时,要达到同样的起涡位置,G auss分布扰动的均方差是最不稳定扰动的350~500倍。可见,提出的最不稳定波入口扰动能够迅速诱发流动本身的K elv in-Helm holtz不稳定,从而提高模拟空间区域的有效性,减小计算量,节约计算时间。     

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

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

摇板式造波机性能实验研究

通过改进后的水池造波系统对工程中常用的规则波和二维非规则波浪的模拟试验，研究了该造波系统的最大造波能力、波浪总谐波变形以及对随机波浪谱的模拟精度等波浪模拟性能．测试数据表明，该造波机控制系统经过技术改进，恢复了二维规则波的模拟功能，提高了造波品质，完全满足船模耐波性试验的要求，可应用于实际工程项目之中。     

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

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

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

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

入口扰动对圆湍射流流动大涡模拟预报的影响

运用大涡模拟方法分别对不加扰动和加扰动的圆湍射流流动进行了数值模拟,比较了两者在大尺度涡结构的演化过程以及时均速度和湍流强度等统计结果方面的预报差异,并与实验结果进行了对照。模拟结果表明,施加了入口扰动的模拟结果在物理真实性方面和统计结果的准确性方面均要明显优于未加扰动的模拟结果。分析了不加扰动和加扰动模拟的本质差异,指出准确预报实际圆湍射流流动应采用施加入口扰动的三维模拟。