水润滑轴承微观界面流场特性的数值模拟

为了研究水润滑高分子材料轴承的微观界面的流场特性,本文以单个微凹体中的流场作为研究对象,建立了单个微凹体流动模型。采用多松弛格子玻尔兹曼方法MRT-LBM (multi-relaxation time lattice Boltzmann method)对单个微凹体的流体流动问题进行数值模拟。研究了流场中速度的变化规律,以及流场中流线分布随雷诺数的变化规律。结果表明当雷诺数较低时,流场处于稳定性流动状态;当雷诺数达到一定程度时,流场处于周期流动状态;当雷诺数很高时,流场处于湍流运动状态。探究微观流场中复杂的变化规律,为研究水润滑高分子材料轴承的润滑和摩擦性能提供一些理论支撑。     

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

可压缩剪切层与边界层流动相互作用的稳定性研究

为了探究可压缩剪切层与边界层组合流动稳定性规律,采用线性稳定性理论对这种组合流动的线性稳定性问题进行研究。平均流动模型取双曲正切函数描述的可压缩剪切层与可压缩相似边界层速度剖面的直接叠加。采用整体数值方法 (Global)求解组合流动线性稳定性方程,获得全部不稳定模态。结果表明,在这样的组合流动中,存在两个不稳定模态,一个与剪切层相关的不稳定模态,另一个是与边界层相关的不稳定模态。在此基础上,研究了剪切层与壁面之间的距离、扰动频率、雷诺数、马赫数、速度梯度、温度梯度对这种组合流动稳定性的影响。马赫数增大对不稳定模态起抑制作用,剪切层对边界层的不稳定模态的促进、抑制作用与马赫数、扰动频率有关,边界层对剪切层不稳定模态基本上起抑制作用。     

不可压缩流动的高精度非标准格子玻尔兹曼方法

首先回顾了高精度非标准格子玻尔兹曼方法的发展历程,基于高精度通量重构格式,发展了一种通量重构格子玻尔兹曼方法（FRLBM）。采用两种求解方法：一种是将碰撞项隐式处理,直接求解离散速度玻尔兹曼方程（直接法）;另一种是先执行碰撞步,再求解纯对流方程（分步法）。通过收敛性研究,比较了这两种方法的精度和稳定性。研究结果表明,在小时间步长下两种方法误差近似,都能取得高阶精度;然而当时间步长增大,直接法误差几乎不变,分步法误差出现明显上升。由此表明,当取得近似误差时,直接法可以采用较大时间步长,计算效率更高,且直接法的稳定性略占优。接着通过模拟顶盖驱动方腔流验证了FRLBM捕捉流场细节的能力,并且比较了基于半隐格式的显式方法和一阶、二阶及三阶隐式-显式Runge-Kutta格式的时间离散,在不同雷诺数下的最大允许Courant-Friedrichs-Lewy数,数值结果表明二阶隐式-显式Runge-Kutta格式效果最优。最后数值模拟了圆柱绕流,验证了FRLBM计算复杂外形绕流的可靠性。     

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

重型燃气轮机高雷诺数CDA叶型转捩特性数值计算

为研究重型燃气轮机的压气机叶片在高雷诺数工况下的气动性能,基于Gamma-Theta转捩模型的雷诺时均方程对某可 控扩散叶型进行了数值计算。通过对比不控制马赫数与控制马赫数,分析高雷诺数对可控扩散叶型气动性能及转捩特性的影响。 结果表明：在不控制马赫数条件下,在零攻角时,雷诺数从7×10 5 增大为9×10 5 ,总压损失增加了约391.95%;在高雷诺数工况下随 着雷诺数的增大,叶片流动损失不断增大,叶片可用攻角范围减小,同时在叶片吸力面出现激波,干扰转捩的产生。在控制马赫数 条件下,当Ma=0.6时,在零攻角工况下,雷诺数从8.2×10 5 增大为1×10 7 ,总压损失减小了约38.98%,吸力面转捩起始点从4.78%弦 长处前移至1.11%弦长处;在高雷诺数工况下,叶片流动损失随着雷诺数的增大不断减小,吸力面转捩位置前移。     

低速流数值模拟的低耗散预处理格式

在预处理Roe格式中加入可调参数,给出了一种低耗散预处理Roe格式,并分析了可调参数的取值规律.在极低速区,低耗散预处理Roe格式反应了流动的实际耗散;在与来流速度相近的区域,低耗散预处理R0e格式可以保证流场的稳定性.对低马赫数、低雷诺数定常圆柱绕流和低马赫数、高雷诺数NACA0012翼型绕流所做的数值实验,得到了比...     

雷诺数对增升装置流动特性影响的计算研究Ⅰ——气动力特性和汇流边界层

在1×10~630×10~6的雷诺数范围内,马赫数为0.197的情况下,使用数值计算方法研究了雷诺数对NHLP-2D翼型的气动力特性和流动特性的影响。建立的数值模型考虑了汇流边界层的网格处理,与已有试验和计算结果对比分析表明本数值模型可信。计算结果表明,当雷诺数大于1.5×10~7时,雷诺数对气动力系数的影响明显减小,且小迎角下气动力随雷诺数呈线性变化趋势。汇流边界层高度随雷诺数增大而降低,缝翼和主翼产生的尾迹强度随雷诺数的增大而减弱,同时尾迹宽度逐渐减小。在高雷诺数下,襟翼尾缘处仍存在较强的缝翼尾迹,说明尾迹/边界层的相互融合作用随雷诺数增大而减小。本文为后续雷诺数对缝道流动特性的影响研究提供了基础。     

高转速转静盘腔流动换热相似性数值研究

为了得到高转速转静盘腔流动换热实验中的相似准则,在定几何、变物性、可压缩和有耗散的情况下,对转静盘腔中的控制方程进行无量纲化,得到了除无量纲位置外的8个无量纲准则数。并以高转速转静盘腔模型为基础,采用数值方法对影响发动机工况（高温）和实验工况（常温）换热相似性的主要无量纲准则数进行了研究。研究表明：流动雷诺数、旋转雷诺数、流动马赫数、旋转马赫数和转子表面温度与来流气流温度之比是影响换热相似性的主要因素。计算结果显示：在热边界条件相似时,仅需保证发动机状态与实验状态的雷诺数相等就可保证局部努塞尔数相同;仅需保证发动机状态与实验状态的马赫数相等就可保证绝热壁温分布相同;但只有保证实验状态和发动机状态的雷诺数和马赫数都相等时,无量纲换热量才保持一致。     

激波结构内流动问题的气体运动论描述

从介观Boltzmann速度分布函数理论出发,发展计及分子粘性碰撞截面与扩散碰撞截面,可描述各流域一维气体流动问题的Boltzmann简化速度分布函数方程及其气体运动论数值计算方法。通过对不同Knudsen数下非定常激波管流动及不同马赫数定常正激波结构问题数值模拟,研究分析不同流区的激波突跃变化过程以及近连续流、稀薄过渡流特有的分子输运现象,揭示不同马赫数、不同分子模型的激波内流动与传热变化规律,证实基于Bo-ltzmann模型方程的气体运动论数值计算方法用于激波结构内流动研究的准确可靠性。     

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

风雷软件LES开发设计与验证

基于国家数值风洞风雷软件开源框架，设计开发了LES(Large Eddy Simulation)湍流模型，主要包括Fourier谱/有限差分方法 LES求解器和有限体积/有限差分方法 LES求解器。简要介绍了采用的不可压缩流动求解的投影法、Fourier谱/有限差分混合方法、亚格子模型等理论方法，给出设计的软件框架和计算流程，尤其说明采用的松/紧2种耦合模式。通过数值模拟不可压缩槽道湍流、亚临界雷诺数圆柱绕流、NACA0012临界攻角的低频振荡算例，验证求解器的计算精度和复杂湍流模拟能力。基于风雷开源框架设计的开源LES模型，具备高精度数值格式、亚格子模型、湍流统计等通用模块，可为国内学者提供一个LES湍流模拟研究的开放平台。     

基于二维紧凑型四阶格子模型的激波管流动模拟

为了提高格子的稳定性，使用Hermite展开方法，构建了新的二维四阶紧凑型格子模型，即D2Q37A。比较了D2Q37A和与Philippi给出的紧凑型格子模型(D2Q37B)的稳定性。在相同的碰撞频率下，与D2Q37B相比，D2Q37A可以模拟初始密度比更高的一维激波管流动。这表明D2Q37A与现有格子模型相比，具有更好的稳定性。详细给出了适用于高阶格子模型的边界条件实现方式。此边界条件实现方式保留了体现LBM（lattice Boltmann method）粒子特性的迁移 碰撞机制。用以上给出的格子模型和边界条件处理方式模拟激波管流动，得到的模拟结果和解析解吻合得很好。这表明所给出的边界处理方式是可行的。此边界格式同样可以用于其他类型的流动和边界。      

微尺度离心叶轮间隙泄漏流动堵塞模型

通过理论分析结合数值模拟,开展了针对微尺度离心叶轮间隙泄漏流动堵塞模型的研究.结合微尺度旋转机械低雷诺数和高马赫数的工作条件,建立了考虑可压缩性和黏性损失的间隙泄漏流动堵塞模型.通过数值模拟分析了微尺度间隙泄漏流动特征及黏性在其中的作用,判定黏性作用对间隙泄漏流动的驱动作用不可忽略,黏性作用与压差作用对间隙泄漏流动的驱动作用效果相当.同时,黏性作用还会导致间隙泄漏流动的初始堵塞增大.采用某微尺度离心叶轮对间隙泄漏流动堵塞模型进行了验证,从数值计算结果中提取气动堵塞结果并与模型预测结果进行对比,表明改进的间隙泄漏流动堵塞模型能更准确地对间隙泄漏流动堵塞进行评估,在微尺度条件下具备更好的适用性.      

低马赫数流动数值模拟方法的研究

在文献[1，2]提出的一种低马赫数流动数值模拟新方法的基础上，本文作者进一步优化了算法，确定了该方法适用的马赫数范围，讨论了参数σ对计算的影响。大量的定常／非定常、低雷诺数流动的计算结果表明这种方法具有计算精度高、收敛速度快以及计算花销小的特点。     

微型飞行器相关技术研究进展

微型飞行器是航空领域一个新的研究热点，涉及多个技术研究领域。围绕与微型飞行器相关的低雷诺数空气动力学问题，进行了绕圆柱体低雷诺数非定常流动的数值模拟方法研究、低雷诺数机翼粘性非定常流动机理研究、扑翼的气动力估算方法研究、微流动控制的数值模拟研究和微型飞机的气动实验研究，指出了微型飞行器面临的问题，并展望了其发展趋势。     

基于双eN方法的短舱层流转捩影响因素

发展自然层流短舱对提升现代民机的经济性和环保性具有重要意义，而对影响短舱层流转捩的因素进行研究有助于更好地开展短舱的层流设计。本文基于线性稳定性分析方法，将双e^N方法同RANS方程求解器耦合，建立了一套可同时计算T-S（Tollmien-Schlichting）波和横流（CF）驻波诱导转捩的流动转捩预测方法，通过标准椭球算例验证了该方法的正确性，进而研究了来流马赫数、雷诺数、湍流度以及迎角对短舱转捩的影响。结果表明：马赫数和迎角会带来压力梯度的明显改变从而引起转捩位置发生变化；而在此构型的高雷诺数工况下，雷诺数和湍流度对转捩位置影响相对较小，转捩位置移动的区域不超过短舱长度的5%。因此在设计阶段，在高雷诺数条件下保持层流设计要尽量避免较大的逆压梯度，保持顺压梯度。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。     

跨声速高雷诺数圆球绕流的实验研究

本文采用无支撑干扰、洞壁影响很小的带加压抽空系统的弹道靶设备,进行了雷诺数在10~5相似文献   

高雷诺数下的翼型绕流LBM数值模拟

采用格子Boltzmann方法(LBM)的二维9速度(D2Q9)模型和贴体网格,通过引入非均匀网格插值方法和非平衡态外推边界处理,分别结合Baldwin-Lomax(B-L)湍流模型和Spalart-Allmaras(S-A)湍流模型,对高雷诺数Re≥5×105下的NACA0012翼型绕流进行了数值模拟和对比研究,两者的结果与CFL3D的结果和实验结果均吻合的很好,相比之下,采用S-A模型能更好地预测失速迎角,其处理分离流动的能力要强于B-L模型。改进后的LBM适用于非均匀贴体网格,曲边边界,计算简单,并可应用于更复杂的高雷诺数流动中。