三维槽道两相湍流的大涡模拟

利用基于动态粘性亚网格模型的大涡模拟方法,分别模拟计算了Ret=180的三维槽道湍流内的三种颗粒的运动状况．统计得到的流体的平均速度、脉动速度等结果与有关文献中的DNS结果一致,不同的颗粒在流场中(特别是近壁区)表现出不同的行为,并再现了颗粒在壁面附近局部富集的现象．     

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

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

从不稳定到近中性大气边界层的大涡模拟

用大涡模拟方法对不稳定到近中性层结的大气边界层情况进行模拟。模拟条件为：地面运动学热通量0.0001～0.2 K·m·s^-1、地转驱动风速0～15 m·s^-1,对应的无因次稳定度参数 - z/L 为0.0062～368。模拟结果显示,大气边界层直到 - z/L＝0.668的近中性情况,其热通量和动量通量等仍保留着典型对流边界层的基本特征,可由模式较好地求解,且在近地面与观测结果相符较好。对更接近中性的情况,模式分辨率限制了对边界层湍流特性的求解。在大风、弱不稳定边界层条件下模拟发现有滚动涡现象出现,其稳定度范围大致对应- z/L 的取值0.668～24.6。     

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

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

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

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

微槽道内气体流动的数值模拟

基于微尺度流动的滑移模型,建立了二维微槽道内气体流动的计算方法,并开发了计算程序.对相关文献中微槽道内的气体流动进行了计算,得到了压力降和速度分布剖面图,计算结果与试验结果吻合较好.在此基础上,进一步研究了克努森数Kn、切向动量调节系数σv及雷诺数Re对速度分布的影响.当Kn>0.001时,边界滑移速度随着Kn的增大而增大;σv的变化对滑移速度的影响十分显著,σv减小时,边界滑移速度增大;在滑流区范围内与低Re时,Re的变化对边界滑移速度的影响很小.     

方柱绕流的大涡模拟

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

离心泵叶轮内部固液两相流动的大涡模拟

为提高离心泵的抗磨蚀能力 ,对泵内固液两相流的流态进行深入研究。在多相紊流运动双流体模型的基础之上 ,建立了两相紊流运动的大涡模拟。利用大涡模拟 ,对离心泵叶轮内部的固液两相流动进行数值计算。计算结果得出液相压力分布和相对速度场 ,及固相的颗粒数分布和相对速度场。其中压力和颗粒的分布与同一离心泵在日本所做合作研究量测的实验结果进行比较 ,两者相互吻合     

稳定分层湍流的大涡模拟

稳定分层湍流是大气和海洋流动中常见的自然现象, 是预测大气和海洋流动的重要物理机制. 我们用大涡模拟方法研究稳定分层流体湍流, 以期获得与实际大气和海洋中稳定分层湍流相符的性质. 我们在不同初始条件的流场和不同的流动特征参数（雷诺数Re和弗罗德数Fr）下用大涡模拟方法获得充分发展湍流场. 通过分析湍动能和湍流特征尺度随时间的演化考察稳定分层湍流演化过程; 并通过三维能谱、水平和垂直能谱来检验数值模拟是否和真实大气、海洋中实测结果一致. 研究结果表明, 大涡模拟可以用较少网格数获得和实测一致的稳定分层湍流性质. 不同的初始流场都能演化为湍流, 并且发现演化过程的后期总是内波能量占优势. 研究证实要获得和实测一致的稳定分层湍流的关键是特征参数必须满足一定条件： 特征雷诺数和弗罗德数不一定要和实际大气和海洋中的参数相等, 但是要求雷诺数足够大, 弗罗德数足够小（例如Re〉102, Fr〈0.1）; 特别重要的是组合参数ReFr2必须大于10.     

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

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

四阶紧致格式有限体积法湍流大涡模拟

为准确预测不可压复杂湍流,提出了一种应用于大涡模拟的高精度有限体积法。该方法在非交错网格上数值求解大涡模拟方程,空间离散采用有限体四阶紧致格式,时间推进用四阶Runge-Kutta法,压力速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动和振荡平板上方的流动,证实了该方法具有近四阶精度;并在此基础上,采用动力Smagorinsky亚格子应力模式,成功地完成了充分发展槽道湍流的大涡模拟计算,所得结果与直接数值模拟结果吻合良好。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

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

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

气固两相圆湍射流流动的大涡模拟

针对气固圆湍射流流动采用基于双向耦合的大涡模拟方法进行了模拟,研究了直径为105μm的颗粒的存在对圆湍射流的影响。单相流动的预报结果在射流形态以及时均速度和湍流强度等统计结果方面均与实验结果符合良好。两相流动的时均速度和湍流强度等统计结果与实验结果符合较好,与单相流动的模拟结果相比,105μm的颗粒加入流场后,由于其较快的响应特征起到了阻尼作用,从而降低了气相场的湍流强度。     

大涡模拟中入流边界脉动风速模拟相关参数取值分析

大涡模拟中入流边界条件脉动风速的处理对计算流体动力学的数值模拟结果的正确性和精度有着至关重要的影响,文章对在运用DSRFG法生成湍流脉动风的算法中,研究了几个相关参数(包括脉动风速能量谱离散化的区间样本数Kmax,能量谱密度模拟所选取三角级数向量组的个数N,空间尺度参数因子θ1和时间尺度参数因子θ2)的数值选取,对入流脉动风速功率谱和脉动风速竖向空间相关性模拟结果的影响.模拟计算结果表明:相对于其他3个参数,空间尺度参数因子θ1的选取相对比较离散.结合DSRFG算法和文章提出的最佳参数取值,模拟得到与实际大气风场特性吻合较好的冯卡门谱及其竖向相关特征.     

大涡模拟预报可靠性分析

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

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

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

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

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