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

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

跨/超临界环境下液氮射流的机理研究

以开源软件OpenFOAM为平台,开发出适于跨/超临界射流的模型,以研究跨/超临界条件下的喷雾混合特性.模型基于真实流体状态方程和扩展的压力隐式分割算法(PISO),并结合雷诺平均(RANS)湍流模型,对跨/超临界条件下的液氮射流进行模拟,重点分析射流的伪沸腾特性.研究结果表明:在跨临界射流伪沸腾发生之前,射流表面形成...     

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

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

用格子Boltzmann大涡模拟方法计算电磁力作用下的翼型绕流

利用格子Boltzmann大涡模拟(LBM-LES)方法,对较大雷诺数Re =2.4×105下翼型绕流的电磁控制进行数值研究.结果表明,LBM-LES方法计算过程简单,容易并行,适合处理该问题.     

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

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

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

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

大涡模拟预报可靠性分析

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

一种改进的旋转湍流亚格子尺度模型及其应用

基于数学物理分析, 发展了一种新型的满足旋转湍流建模规则(即模型的渐近物性坐标不变性条件)的动力学亚格子尺度应力模型, 并针对展向旋转槽道湍流问题, 计算验证了该模型的正确性. 文中采用大涡模拟方法和该动力学亚格子尺度模型, 进一步研究了展向旋转槽道湍流的统计量和大涡Reynolds应力输运方程主要项的变化特性, 分析了旋转槽道近壁区的湍流拟序结构.     

超临界流体跨临界循环最优压力研究

介绍了超临界流体跨临界循环与传统亚临界循环的不同特点,对跨临界循环的最优压力进行了理论分析和研究,对超临界流体的跨临界循环的最优压力影响因素进行了分析和比较,探讨了超临界流体跨临界循环在热泵和干燥领域的应用前景.     

显式滤波大涡模拟算法

在传统大涡模拟算法的基础上引入显式滤波技术，有效地控制了截断误差．对对流项和湍流应力项进行二次显式滤波，并对湍流应力的分解形式作了相应改进．以三维平板通道为计算模型编写了二阶精度计算程序，对传统、显式滤波大涡模拟2种算法的计算结果进行了对比，并就亚滤波模型的重建程度对计算精度的影响进行了探讨．计算结果表明：2种算法对流线方向平均速度分布模拟差异不大，但显式滤波大涡模拟算法可明显改进脉动速度分布；随着对亚滤波模型重建程度的提高，平均速度和脉动速度分布均有明显改进．     

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

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

涡轮桨搅拌槽内混合过程的大涡模拟

在FLUENT 6.1软件平台和网络并行计算硬件平台上,采用大涡模拟(LES)的方法对涡轮桨搅拌槽内的混合过程进行了数值模拟。利用滤波函数对N av ier-Stokes方程进行空间滤波,对大尺寸的涡直接进行求解,而被滤掉的比网格小的旋涡通过Sam agorinsky-L illy亚格子模型求解,对搅拌桨区域采用滑移网格技术。结果表明:大涡模拟对尾涡的预报优于雷诺平均(RAN S)方法,混合时间以及示踪剂响应曲线模拟结果和实验结果吻合较好,且优于RAN S方法。大涡模拟方法为准确预测搅拌槽内湍流流动的非稳态及周期性脉动特性提供了一种有效的工具。     

大涡模拟理论及其应用综述

对大涡模拟(Large Eddy Simulation-LES)理论的发展过程、基本理论及其应用现状和前景进行了系统的论述,并概要地介绍了运用LES理论,数值模拟方柱绕流的流动分离现象和动量射流与线性波相互作用的数值模拟成果,表明了LES理论对紊流细部结构具有较强的数值模拟能力。可望用以解决大量的工程实际问题。     

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

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

超临界流体萃取过程混合模型的建立

采用混合建模技术,通过集成人工神经网络模型与单元过程机理模型建立了超临界流体萃取过程的混合模型,较好地解决了超临界过程的模拟计算、经济评价、设计和运行优化问题.通过软件集成的方法实现了商业模拟软件Aspenplus与ANN的集成,大大降低了模型的实现难度和实施时间.     

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

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

方柱绕流的大涡模拟

采用物理分步法对方柱绕流的流动分离现象进行了大涡模拟(LES).物理分步法将N-S方程剖分为对流、扩散及波动过程,对不同的物理过程分别采用适宜的方法进行求解.对流过程采用逆风——中心差分格式求解;扩散过程采用中心差分格式求解;压力方程采用SOR迭代法求解.计算的方柱四周的平均压力系数的分布与物理模型试验结果吻合良好;模拟的流动现象也得到了物理模型试验的证实.     

热力羽流的大涡模拟

采用大涡模拟的方法对两无穷大平板间受热形成的羽流进行了数值模拟.控制方程为三维不可压N-S方程组,数值格式采用了高精度的伪谱方法,亚格子模型则采用了TNS模型.首先对槽道流进行了模拟,通过比较大涡模拟和直接数值模拟的结果对此方法进行了验证.随后数值模拟了热力羽流的输运过程,分析讨论了相关的温度场和速度场,并比较了不同Rayleigh数对流动特性的影响.结果表明当Rayleigh数比较大时,热浮力对流动的影响比较大,形成了比较复杂的流动结构.     

方柱绕流的大涡模拟

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

气固流动大涡模拟和两相湍流模型的评价

用基于气体Smagorinsky亚网格应力模型和颗粒动力论模型的双流体大涡模拟(LES)以及统一二阶矩两相湍流模型(USM-RANS),对轴对称突扩气固流动进行了数值模拟。结果表明:LES模拟的颗粒轴向平均速度和均方根脉动速度,以及USM-RANS计算得到的模拟结果都与实验结果吻合较好。LES得到的颗粒-气体轴向脉动速度关联值比USM-RANS模拟值更接近实验值,这表明USM-RANS模拟还有改进余地,并基本上得到了LES的验证。LES的瞬态模拟能显示气固湍流流动的各向异性两相湍流结构和颗粒弥散的发展过程,而USM-RANS无法实现。