显式滤波大涡模拟算法

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

湍流大涡模拟的时空关联方法

本文介绍了湍流大涡模拟的时空关联方法,它用于检验和发展大涡模拟的亚格子尺度模型,其目的是提高大涡模拟方法对于欧拉和拉格朗日时空关联的预测准确度.现在广泛使用的湍流涡黏亚格子模型的理论基础是能量平衡方程,它们能够正确预测湍流的空间能量谱,但不能正确预测湍流的时间能量谱.本文针对大涡模拟的时间尺度问题,从时空耦合的观点出发,综述了作者近年来在时空关联方法方面的研究工作.理论分析和数值模拟结果表明,大涡模拟总是高估湍流的时间尺度,其物理机制是涡黏亚格子模型仅能耗散小尺度能量,但不能反映亚格子涡对大尺度涡的随机反馈作用.小尺度湍流会影响颗粒的聚集程度和相对扩散.建立亚格子湍流的拉格朗日时间尺度模型,为携带颗粒湍流大涡模拟中基于随机微分方程的颗粒相亚格子模型提供封闭参数,并用于颗粒湍流扩散的大涡模拟.     

群体互扰下的风场绕流特性和风压分布分析

研究群体互扰下的风场绕流和风压分布问题.改进了流场模拟大涡模型的亚格子应力的湍动黏度,以改善大涡模型对壁面湍流的描述,在壁面处得到与直接数值模拟(DNS)相近的结果.基于此种改进的大涡湍流模型,在单物体和3个物体绕流情形下,分别选择结构刚性模型和弹性模型开展了流场和风压的数值模拟,得出相应的绕流特性和风压分布并与运用其他湍流模型的计算结果进行对比分析,以验证本文提出的改进技术.同时,计算表明,单物体和3个物体情形下的流场特性和风压分布显著不同.     

旋转槽道中Taylor-G6rtler涡对雷诺应力输运的影响

我们通过直接数值模拟研究了旋转槽道湍流中Taylor-Gortler涡对雷诺应力输运的作用．采用流向平均的方法检测Taylor-Gortler涡，并定义了TG脉动．我们发现，和无旋转湍流槽道不同，旋转槽道湍流的吸力面湍流扩散对雷诺应力的输运是主要的；其中Taylor—Gortler涡在湍流扩散中起重要作用．我们关注中低旋转数的旋转槽道流动，并阐明旋转数对雷诺应力输运的影响．     

单流道泵内部流动研究方法探讨

该文将单流道泵叶轮和涡壳作为一个整体，利用势流理论计算了单流道泵内部流动，对流场进行初步分析，得到了叶轮转至不同角度时水对叶轮的径向力以及水功率．同时，作者还利用湍流模型对单流道泵内三维湍流流动进行了数值模拟，并通过两种方法的结果比较和与实验对比以及对应用前景分析等验证了研究方法的实用性．     

考虑旋转和曲率修正的螺旋桨梢涡流场模拟

应用显式代数雷诺应力湍流模型对螺旋桨尾流中梢涡流场分布进行了数值研究,为了避免过高地预报梢涡涡核内湍流黏性耗散,对湍流模型进行了旋转和曲率修正．应用全六面体网格对螺旋桨计算域进行网格划分,为了避免数值离散误差,对梢涡区域进行了网格加密处理．计算结果表明：提出的尾流中梢涡流场分布数值模拟方法能够准确预报螺旋桨梢涡流场的分布及涡核位置,并准确反映了梢涡形成和发展过程中梢涡内主涡和次涡的关系,与实验测量结果基本一致．     

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

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

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

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

大涡模拟在水利中的研究进展

大涡模拟作为一种研究湍流流动的有效手段，在国际上已经得到了广泛的应用．本文在介绍了大涡模拟的发展历程及其理论方法的基础上，就国内外大涡模拟在水利中的最新研究成果作了较为详细的阐述和系统的分析，指出了其中应该注意的问题和今后的发展方向。     

高雷诺数湍流的熵格子玻尔兹曼方法大涡模拟计算

为了开发基于熵格子玻尔兹曼方法的通用计算流体动力学求解器,建立了熵格子玻尔兹曼方法的大涡模拟标准Smagorinsky亚格子尺度湍流模型的计算模型,探究了其进行高雷诺数湍流模拟的有效性,并且尝试展开了基于熵格子玻尔兹曼方法大涡模拟亚格子尺度模型的近壁面研究.得到结论:熵格子玻尔兹曼方法的大涡模拟能够有效进行高雷诺数湍流的数值模拟计算;近壁面处理的Catalano方法对壁面处流场速度的计算结果优于Van Driest方法,而Van Driest方法的优越性表现为模拟计算方法的一致性.     

基于RANS和LES湍流模型的二甲醚超临界喷雾数值模拟

分别采用RANS和LES湍流模型耦合KH-RT破碎模型对二甲醚超临界非燃烧喷雾实验进行数值模拟,研究两种湍流模型的适用性,并对比分析二甲醚亚临界喷雾和异辛烷超临界喷雾的发展过程。同时,基于LES湍流模型探讨了环境温度和压力以及喷射温度和压力的变化对二甲醚超临界工况下喷雾贯穿距的影响。结果表明,RANS和LES湍流模型均能较好地模拟二甲醚超临界喷雾,LES湍流模型的模拟结果与喷雾实验观察结果更吻合;与二甲醚亚临界喷雾相比,其超临界喷雾的雾化效果更好;在相同初始条件下,二甲醚超临界喷雾的雾化效果要优于异辛烷超临界喷雾;二甲醚超临界喷雾贯穿距随环境温度、燃油喷射温度和喷射压力的升高而增大,随环境压力的升高而减小。     

不稳定近地面层湍流的大涡模拟

采用加密网格的大涡模式获取边界层风、温场的高分辨率模拟结果,并据以分析近地面层大气的湍流特性。结果表明,较小的网格尺度使次网格湍流贡献率大为降低,模式计算结果对次网格参数的依赖性减小,边界层整体特征得到更好的反映。同时,模拟出的近地面层通量 廓线关系及湍流速度特征与实际观测结果吻合甚好,表明模式具有反映近地面层平均运动和湍流特性的能力。     

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

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

基于谐波合成法的大涡模拟脉动风场生成方法研究

为准确模拟大涡模拟入口处的脉动信息,在充分考虑脉动风场的功率谱、相关性、风剖面等参数前提下,运用谐波合成方法生成了满足目标风场湍流特性的随机序列数,通过对FLUENT软件平台进行二次开发,将生成的随机序列数赋给大涡模拟的入口边界,从而实现了大涡模拟的脉动输入.基于风洞试验数据,分别建立了两种模拟脉动风场的数值模型,一为没有任何障碍物的空风洞,运用谐波合成方法生成的随机序列数作为入口边界来生成脉动信息;二为与真实风洞一致的尖劈粗糙元风洞,采用平均风作为入口边界,利用尖劈粗糙元对风场的扰动来产生脉动信息.通过对比两种数值模型发现:基于谐波合成方法生成的脉动风场可作为大涡模拟的入口边界,可为大涡模拟脉动入口研究提供参考.     

基于大涡模拟的湍流非预混燃烧混合分数概率密度函数

混合分数概率密度函数(probability density function,PDF)反映了湍流对燃料和氧化剂混合过程的影响,在湍流非预混燃烧的理论研究和数值模拟中有非常重要的作用。该文基于大涡模拟(large eddy simulation,LES)对非预混火焰中的混合分数PDF进行了研究。利用LES预测的SandiaFlame D的速度和温度的均值和均方根分布与实验结果符合很好,瞬态温度场显示了合理的湍流火焰形态。混合分数PDF在反应区为钟形分布,在贫燃侧和富燃侧为钟形分布或单调形分布,取决于当地流场状态。对简化PDF模型的研究表明:β函数模型对钟形PDF和单调形PDF的预测效果都很好;截尾Gauss函数模型只能较好地预测钟形分布PDF;多点δ函数模型的预测能力与截尾Gauss函数模型的预测能力类似;双δ函数模型的预测结果偏差较大。     

城市冠层植被大气环境特性大涡模拟

 植被通过生物物理过程影响着周围微气候环境,对于城市植被环境流动的研究需要建立合理的植被模型。本研究发展了一套适于模拟植被边界层流动的大涡模拟数值方法,并通过水平均匀模型植被算例对其进行了验证。在此基础上,对空气流过三维非均匀分布钝体元阵列模型植被进行了大涡模拟,并将不同植被数值模型的计算结果与相应的风洞实验数据进行了比较。结果表明,树干部分采用刚性圆柱模拟、树冠用阻力模型的“阻力冠”方法相比于传统的“阻力元”方法能更好反映植被流动的影响。     

基于虚拟叶片方法的复杂地形风电场气动性能分析

基于计算机流体力学(computational fluid dynamics,CFD)方法和虚拟叶片(virtul blade method,VBM)方法研究了风力机尾迹与复杂地形的相互作用,以实际风电场为例,分析了造成复杂地形风电场中风力机振动的主要原因,提出了一种基于CFD的快速分析风电场气动性能的方法,比较了雷诺平均(Reynolds-averaged,Navier-Stokes,RANS)和大涡横拟法(large eddy simulation,LES)求解方法在捕捉地形效应和尾迹及预测风电场功率性能方面的差异。研究分析发现,机组尾流和复杂地形效应造成的湍流强度过大是造成该风电场风力机运行时剧烈震动的原因,LES方法预测尾流的延迟衰减和更明显的尾流干扰作用导致风电场的发电量低于RANS方法。RANS方法过度预测湍流,导致湍流动量更快的扩散和尾流快速恢复。因此,CFD结合VBM分析复杂地形条件下风电场气动性能的方法,可实现复杂地形大型风电场性能的快速分析,可应用于风电场相应气动控制策略的制定和复杂地形风电场的微观选址。     

亚临界雷诺数下圆柱体尾流结构的数值模拟

在典型亚临界雷诺数下,运用基于Smagorinsky亚格子模型的大涡模拟(LES)方法,对均匀来流作用下的圆柱体绕流问题进行三维数值模拟.通过比较基于时间积分的平均积分分量,验证了数值计算结果的准确性,然后分析了尾流区的流场结构,并给出了平均流场以及湍流流场的流场特征.此外,解释了数值结果与早期试验结果间存在微小差异的原因,同时指出回转长度是尾流结构特性中最重要的特征参数.     

自然通风条件下炼钢厂煤气泄漏扩散研究

为评估炼钢厂中突发煤气泄漏的后果,通过建立二维简化平面模型,采用LES湍流模型对煤气在厂区的泄露扩散过程进行了数值模拟,得到了CO浓度随时空的分布规律.计算结果表明,在泄露过程中,泄漏口近的CO浓度分布比较集中;泄漏停止后,在自然风的作用下CO浓度并不能很快降低,而是呈现一个平缓的发展过程.本结果对突发煤气泄漏的报警和预案的制定有一定参考价值.     

Hydrodynamics numerical investigation of hoistable masts for underwater vehicles

Using the unsteady incompressible Navier-Stokes equation as the governing equation, the large eddy simulation (LES) model is implemented to investigate the shedding of vortices, the flow pattern of turbulence, the unsteady pressure fluctuation and the time history of the lift coefficient and drag coefficient of hoistable masts with various mast shapes and various arrangements in this paper. Combining the FFT, combined time-frequency transform and wavelet power spectrum analysis, the characteristics of unsteady pressure can be obtained in both time and frequency domain. It shows that the main frequency of pressure fluctuation is near the frequency of vortex shedding in time domain using the FFT method. It can be inferred from the combined time-frequency transform that the unsteady pressure fluctuation has obviously the peak value and the second peak value in time domain. It could indicate that the fluctuation power varies from the fluctuation frequency through the power spectrum analysis. By the data analysis, it shows that the vortex shedding is the dominant cause of the periodically pressure fluctuation. And the interaction pattern of wake and interplay between wake and the walls of masts under different arrangements are also discussed in this paper.