光滑及粗糙明槽湍流流动大涡模拟

采用大涡模拟（Ｌａｒｇｅ Ｅｄｄｙ Ｓｉｍｕｌａｔｉｏｎ）方法和Ｓｍａｇｏｒｉｎｓｋｙ亚格子模型数值计算了光滑及粗糙明槽湍流运动。以槽高为参考长度的雷诺数分别为３０００，２２８００和１０＾５，数值离散方法为有限体积法，自由面的处理使用了简单了刚盖假定，固壁面采用了壁函数，与低雷诺数直接模拟（Ｄｉｒｅｃｔ Ｓｉｍｕｌａｔｔｉｏｎ）和物理实验结果相比，本文数模呈现了良好的预测效果，估计的湍流宏尺度（Ｍ     

子格（SGS）模型在内流湍流中的大涡模拟

本文将子格模型应用于大涡模拟，对直方管，弯方管和弯曲槽道内充分展湍流，进行了数值模拟，并将所得的结果与实验数据作比较，计算表明，应用于子格模型的大涡模拟可以用来模拟比较复杂的内流，研究湍流运动的结构，检验湍流模型。     

槽道中方形障碍物绕流的大涡模拟

本文应用二阶全展开ETG有限元离散格式与大涡模拟相结合的方法对固定在槽道一侧的方形障碍物的绕流进行了数值模拟，计算了雷诺数为40000的湍流情况下的绕流流场，将模拟结果与实验结果及他人的数值结果进行了对比，符合较好，表明大涡模拟与所论的具有较高精度的数值算法相结合，适合于具有大尺度涡的绕流运动流场的分析。     

双圆柱绕流的大涡模拟

基于FLUENT软件平台上的三维大涡模拟方法,对定常来流流场中的单圆柱及不同排列形式的双圆柱进行了水动力计算,并对圆柱周围的瞬时流场、时均流场以及圆柱受力进行了分析.分析结果表明,当Re=6×104时,串列双圆柱中上游圆柱受力明显大于下游圆柱,下游圆柱Gd、Cl的变化频率满足fd=2fl,两圆柱的脱涡频率近似相等,涡街的相位相同;并列双圆柱的Cd基本相等,且变化趋势有一定相似性,脱涡频率均接近单个圆柱,涡街的相位相反;交错45.排列的两个圆柱Cd基本相等,脱涡频率接近,下游圆柱的涡街相位明显滞后于上游圆柱.     

大涡模拟法的二维形式

本文给出了大涡模拟（ＬＥＳ０法的二维简化形式，选择了标准的后台阶湍流进行考题。计算结果与有关实验十分吻合。为大涡模拟法的工程应用打下了基础。     

绕方柱流脉动压力场的大涡模拟

采用有限差分法，对雷诺数为22000的绕方柱流的脉动压力场进行了三维大涡模拟。采用时间分解算子法，将N－S方程分为对流分步、扩散分步和传播分步，再分别进行计算，在大涡模拟中，用两层模型处理Smagorinsky模型在近壁区模拟结果的发散问题。计算的压力场以及阻力系数、升力系数和Strouhal数均与试验结果较为吻合。     

叶顶间隙对轴流泵轮缘泄漏流动影响的大涡模拟

为揭示不同尺寸叶顶间隙下轴流泵轮缘区湍流特征,采用大涡模拟方法对轴流泵在设计流量下的内部非定常流动进行了数值计算,分析了5种叶顶间隙下泵轮缘间隙区的流场结构。通过泵外特性参数预测值与实验结果的对比,证实本文所提出的方法可较准确反映泵内流动特征。结果表明:随轮缘间隙增加,泵扬程、功率和效率均呈下降趋势;叶顶间隙内泄漏流速度沿径向逐渐增大;随间隙尺寸增加,轮缘间隙主泄漏涡强度及与叶片夹角均增大,泄漏涡的非稳定性增强;当间隙δ/D2大于1.0‰,间隙主泄漏涡发展至相邻叶片正面;当δ/D2大于1.5‰,叶顶间隙内出现间隙分离涡,间隙区形成多个次泄漏涡,其影响范围随间隙尺寸的增加而增大。     

三方柱绕流的大涡模拟及频谱分析

高雷诺数条件下,无论来流是均匀还是非均匀,多方柱绕流情况下方柱受力及尾流的相互干扰都是相当复杂的。为了研究方柱受力及下游尾流的相互干扰,基于大涡模拟紊流模型对后品字等边布置的三方柱绕流进行了数值模拟,数值模拟结果表明,上游两个方柱的阻力系数要明显小于下游方柱,而上游两方柱升力系数远大于下游方柱,且下游方柱的升力系数基本上在零左右振荡,通过频谱分析发现了它们之间的区别。方柱的尾流存在着强烈的非对称相互干扰,且涡有明显的三维性。     

三方柱绕流的大涡模拟及频谱分析

高雷诺数条件下,无论来流是均匀还是非均匀,多方柱绕流情况下方柱受力及尾流的相互干扰都是相当复杂的.为了研究方柱受力及下游尾流的相互干扰,基于大涡模拟紊流模型对后品字等边布置的三方柱绕流进行了数值模拟,数值模拟结果表明,上游两个方柱的阻力系数要明显小于下游方柱,而上游两方柱升力系数远大于下游方柱,且下游方柱的升力系数基本...     

Large eddy simulation of water flow over series of dunes

Large eddy simulation was used to investigate the spatial development of open channel flow over a series of dunes. The three-dimensional filtered Navier-Stokes (N-S) equations were numerically solved with the fractional-step method in sigma coordinates. The subgrid-scale turbulent stress was modeled with a dynamic coherent eddy viscosity model proposed by the authors. The computed velocity profiles are in good agreement with the available experimental results. The mean velocity and the turbulent Reynolds st...     

大涡模拟水环境中污染物团的运动规律

本文采用大涡模拟方法封闭引入浮力特性并经空间平均后的Navier-Stokes紊流方程，对水环境中污染物的运动规律进行研究。在用有限差分法对控制方程进行求解时，用破开算子法将控制方程分为对流、扩散和源项三部分，其中对流项用QUICKET程式求解，扩散项用中心差分求解，源项用Gauss-Seidel迭代法求解。数值计算结合物理模型结果表明，污染物团在水环境中运动除有时均定量关系之外，还有三大定性规律：第一，污染物团在时间和空间上具有不规则性；第二，污染物团后部是空的；第三，在同一条件下进行多次试验（数学和物理），污染物团的运动是不一样的。计算表明大涡模拟要优于通常的 κ-ε 紊流模型，它不但能够象 κ-ε 紊流模型那样计算出污染物在水环境中的时均运动规律，而且能够计算出其在空间上的不对称发展过程，能够准确地再现其在物理模型中的不对称运动特性。计算结果无论从定性还是从定量关系上都与试验结果吻合一致。     

特征线算子分裂有限元的圆柱绕流大涡模拟

为研究圆柱绕流流场特性,将大涡模拟与特征线算子分裂有限元相结合,建立了大涡模拟特征线算子分裂有限元模型,对单圆柱和串列双圆柱绕流问题进行了数值模拟,所得结果与现有研究结果吻合较好。模拟结果表明:对单圆柱绕流,随着雷诺数的增大,圆柱近尾流区上下交替的涡旋逐渐靠近通过圆柱几何中心的水平线,且涡脱落位置逐渐靠近圆柱。对Re=1 000的串列双圆柱绕流,临界间距在圆柱直径的2.25~2.5倍之间;当两圆柱间距小于临界间距时,上游圆柱后方无明显涡旋脱落,间隙处压力较稳定;大于临界间距时,有涡旋脱落,上游圆柱尾流区上下表面交替出现强负压区。     

不等直径串列圆柱绕流大涡模拟

为研究背负式海底管线中增设的小直径附属管线对主管线的水动力影响,将大涡模拟中经典Smagorinsky亚格子模型与特征线算子分裂有限元法结合,并引入出口对流边界条件,完善了基于特征线算子分裂有限元的大涡模拟方法。通过自编程序数值模拟Re=1 000的单圆柱绕流,计算结果与相关文献吻合较好,验证了该算法计算圆柱绕流的有效性,并分析了Re=1 000时不同直径比、间距比情况下的串列双圆柱绕流,根据流场的不同涡脱落形态及两圆柱平均阻力系数、升力系数随直径比、间距比变化的规律得到了不同直径比条件下的临界间距范围。达到临界间距后,流场由单一涡脱落状态转变为双涡旋脱落状态。最后分析了两圆柱平均阻力系数及升力系数在临界间距后急剧增加的原因,为背负式海底管线的布局优化提供了理论依据。     

不同长细比圆柱绕流的大涡模拟

自从人们对层流圆柱绕流现象有了系列研究及清楚的认识后,逐渐把目光投向湍流的圆柱绕流问题,但高雷诺数下的圆柱绕流具有很强的三维特性,因此在数值计算时模型轴向长度的选取是大家关心的问题。该文基于开源软件OpenFOAM运用LES计算了4种不同长细比圆柱的静止圆柱绕流(包括一个二维算例),重点分析和对比了结果中的一些基础参数如St数、平均阻力系数等,与实验结果吻合良好。文中还进一步分析了时均顺流向速度U在流场中沿流向和横向的分布,发现长细比为π的圆柱算例中,计算域的轴向尺度已经能满足圆柱下游靠近圆柱处的流场中的三维结构的发展。     

用大涡模拟方法计算尾水管内非定常周期性湍流

本文运用大涡模拟-双方程模型进行尾水管内部非定常周期性湍流计算。而且本文对尾水管复杂计算域的网格划分，提出一种更加合理的多块结构网格。计算结果用波形显示了尾水涡带和压力脉动的主要性质和运动特征。     

非对称复式断面明渠中污染物输移扩散的大涡模拟

运用大涡模拟计算了两个非对称复式断面明渠内的水流运动和污染物扩散，其中由亚格子尺度涡运动引起的亚格子应力与物质扩散通过Smagorinsky-Lilly亚网格模型模拟。预测结果与Shiono等使用激光多普勒流速计(LDA)和激光诱导荧光技术(LIF)测量的实验数据进行了对比，计算得到的流场结构和污染物浓度分布均较好地吻合于测量结果。成功地再现了在其中一个复式断面明渠中观测到的横向浓度分布呈现显著双峰值的现象，根据计算的流场得出复式断面明渠中流向交替变化的强烈横向流动是产生这种现象的主要原因。     

横流中射流及其温度标量输运大涡模拟

将动力拟序涡黏性亚格子应力模型拓展到温度标量亚格子模型中,数值模拟了横流条件下有、无温度标量场作用的射流,得到的横流条件下浮力射流的温度和速度分布与Anwar的试验值吻合一致。在此基础上,分析了有、无温度标量场作用下射流回流区域大小和射流轨迹线特性,对比分析了回流区域内涡心和分离点处湍动能和耗散率、拟涡能以及边界层处拟序结构等湍流特性。计算结果表明:温度场的作用使射流的回流区域增大,射流速度轨迹线高度增加,回流区域内湍流的湍动能增加,边界处拟序结构的周期性不如无温度场时明显。