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

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

用格子Boltzmann方法模拟MKDV方程的行波解

用精确到0(ε^4)的5速格子Boltzmann模型较系统地对MKDV方程：ut θu^2ux βuxxx=0的行波解进行了模拟，并将模拟的数值结果与相应的理论结果进行了比较，模拟值与理论值吻合很好。     

Burgers方程定态激波解的格子Boltzmann方法模拟

用一维５速格子Ｂｏｌｔｚｍａｎｎ方法推导一维Ｂｕｒｇｅｒｓ方程,采用Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ展开的局域平衡分布函数形式并由守恒方程确定其展开系数。最后用格子Ｂｏｌｔｚｍａｎｎ方法模拟了对流系数ｖ＝１,扩散系数Ｄ＝０．３３３３的一维Ｂｕｒｇｅｒｓ方程的激波解。计算机数值模拟结果与理论解析解精确吻合。     

基于格子Boltzmann法的树木流场数值模拟

通过多松弛时间格式的格子Boltzmann法(lattice Boltzmann method,LBM)对树木流场进行数值研究,树木简化为表征元尺度下的多孔介质,为了研究雷诺数Re=10 000下的流场,将大涡模拟结合LBM得到树木周围的流速和风压,并在模型前设置格栅,以此制造扰动并模拟B类场地下的流场。结果表明,树木下游不同位置的平均风速剖面与试验吻合良好;风压从迎风面到背风面显示出明显的变化梯度,而且其最大绝对值主要分布在迎风面和背风面附近; B类地貌下的树木遮蔽区下游风速有所减小,湍流度极大值在树冠顶部附近。     

Burgers方程的格子Boltzmann方法模拟

用一维格子Boltzmann方法构造O(ε^4)的Burgers方程模型．格子Boltzmann方法的数值模拟结果与具有特定边界条件的Burgers方程的解析解精确吻合．     

偏格子中的格子Boltzmann方程

给出偏格子的格子Boltzmann模型. 应用Chapman-Enskog展开和多重尺度技术, 通过 选择平衡态分布函数的高阶矩和偏张量的形式, 得到偏网格上的Euler方程.     

二维反应扩散方程的格子Boltzmann方法模拟

建立了二维反应的散方程五速四方格子模型的格子Boltzann方法，进行计算机模拟演化，分析模拟结果，并与SChen利用七速六方格子模型得到的结果进行比较，得出在模拟求解流速为零情况下的反应扩散方程时，这两种格子模型是一致的结论。     

结合大涡模型的等离子射流Lattice Boltzmann模拟

基于一种结合大涡模型的Lattice Boltzmann方法对等离子体射流湍流流动状态进行模拟，并将模拟结果与实验结果和其他文献的计算结果进行对比。结果表明，与传统的差分法相比，这种混合Lattice Boltzmann方法收敛速度快，且能获得比较满意的结果。     

稳定分层湍流的大涡模拟

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

大涡模拟预报可靠性分析

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

方柱绕流的大涡模拟

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

中心进气式鼓泡反应器内气液流动的大涡模拟

研究鼓泡塔内涡的结构及其演变规律，对揭示反应器内传热、传质性能具有重要意义．用大涡模拟（LES）方法对中心进气式鼓泡塔内瞬态流动性能进行了模拟研究，结果表明，增大气速可使涡结构剧烈变化，增加液体黏度限制了涡尺度的发展，减少气泡尺寸加快了涡的变化速率．模拟结果与实验值吻合良好，大涡模拟方法能有效地揭示鼓泡塔内不同尺度的涡结构及其发展变化规律．     

小尺度新月形沙丘背风侧流场特性的大涡模拟分析

为了准确获得新月形沙丘背风侧湍流流场流动特性,采用基于Smagorinsky的亚格子尺度涡黏模型的大涡模拟(LES)方法,对简化了的实验缩比小尺度沙丘模型绕流气相流场进行了数值研究,模拟了相同来流风速、不同沙丘迎风坡度和高度下背风侧风场的湍流流动模式,获得了背风侧区域地表面摩擦系数分布,比较了回流区长度和湍流强度分布.模拟结果表明:LES方法能较好地揭示背风侧湍流流场特性;流场重附距离随沙丘迎风坡度和高度的增加而增大;背风侧回流区内湍流强度总体上比回流区外大,回流区内沙丘坡脚位置及重附点位置处湍流强度最大值出现在贴近地表附近,随沙丘高度的增加而增大;回流区中部湍流强度最大值出现在回流区顶部,随沙丘高度的增加先增大后基本不变.     

类车体气动性能的大涡模拟

以MIRA车体气动性能的风洞试验数据为基础,对采用大涡模拟方法解算非定常特征显著且具有大分离流动结构的近地钝体外部绕流场所涉及的迭代步数、时间步长、网格方案等影响因素开展研究.采用对比分析方法对3种亚格子湍流模型的计算准确性进行研究.提出适用于三厢车型的大涡模拟数值仿真策略.     

存在热浮力的剪切湍流中颗粒近壁运动的大涡模拟

采用基于动力学亚格子模型和热通量模型的大涡模拟方法,结合Lagrange粒子追踪技术,数值研究了存在热浮力作用的剪切湍流中细颗粒群的输运特性.一方面探讨在槽道剪切湍流中不同热分层效应与粒子输运的动态关联性,以及不同Stokes数的细颗粒在同一分层湍流环境下的浓度分布;另一方面探讨动量和热量耦合作用下影响粒子迁移、扩散和沉积的不同因素.此外,基于一些典型计算结果,探讨了稳定和不稳定分层对颗粒运动各向异性程度的影响,研究了不同颗粒的平均速度、脉动强度的统计特征量以及与湍流流场之间的相关性,并对颗粒在槽道不同区域的体积浓度分布和局部积聚现象进行了分析.     

恒定有压扩散流的大涡模拟

采用大涡模拟（LES）方法对恒定有压扩散流进行了数值模拟，对大尺度漩涡直接求解Navier-Stokes方程，小尺度漩涡采用标准Smagorinsky亚格子模型模拟．针对不同雷诺数下扩散段内的水流进行数值模拟，得到了各雷诺数流动时的时均速度场，分析了不同流态下扩散段内恒定有压扩散流时均流动特性．通过对不同流态下扩散段内连续瞬态流场的分析，探讨了不同雷诺数流动时扩散段内瞬态流场结构特点以及分离区内漩涡的变化规律，指出在恒定流状态下扩散段内存在流动的局部不稳态现象（主流的非周期性摆动），并对此现象进行分析．计算结果和试验结果吻合较好．     

大跨度屋盖风荷载的大涡模拟研究

应用一种新的湍流脉动流场产生方法DSRFG(discretizing and synthesizing random flow generation)模拟风场实际的湍流边界条件.采用一种新的大涡模拟LES(Large Eddy Simulation)的亚格子模型,基于linux系统下软件平台Fluent6.3的并行计算技术,计算了长沙机场扩建航站楼屋盖结构在5个风向下的风荷载.根据对机场屋盖的平均风压和脉动风压系数分布规律的分析,给出了各风向角下屋面的表面风压分布特性和最不利风向角,为长沙机场扩建航站楼的抗风设计提供了依据.并为进一步发展此类复杂建筑结构在复杂湍流环境下的风荷载数值风洞技术提供参考.     

应用大涡模拟对水平管道进行流场分析

为了准确描述湍流三维非稳态随机波动这一高度非线性的特性，简述湍流物理现象的特点，引出基于不同角度建模的湍流模型，在理论上对大涡模拟进行分析，并运用并行计算机、FLUENT6．2软件作为计算平台，针对水平管道中两股不同入射速率的流股相互作用下的流场进行大涡模拟数值分析，并与标准模型进行比较．