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

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

MRT-LBM三重网格局部加密算法研究

应用MRT-LBM计算流场中运动粒子上的力和力矩时,可使用三重网格局部加密方法来提高计算结果的精度.将流场计算区域划分为互不重合的粗糙区、过渡区和加密区,相邻区域界面重合节点上的物理参数如密度、速度和应力等保持连续,不同区域的分布函数通过过渡区边界点进行传递,运算程序采用碰撞-迁移算法,给出了数据的初始化过程和加密的网格结构.通过对含非对称放置的粒子Couette流动,分别采用标准Boltzmann、全场加密、局部二重加密和三重加密的4种网格进行计算.数值计算结果表明,三重网格加密算法减少了计算结果的波动.对方腔流顶盖下方左右奇异角落处应用三重局部网格加密,结果显示,对雷诺数1 000的方腔流动,沿腔中心线的速度分布与经典文献结果对比效果良好,压力轮廓图的噪声明显降低,应力振荡明显减少.模拟结果证实了所建立的局部加密方法的有效性.      

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

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

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

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

Cavitation flow simulation for a centrifugal pump at a low flow rate

A computational model on the basis of the combination of the full cavitation model and the modified re-normalization group (RNG) k-ε turbulence model are used to simulate the cavitating flow in a centrifugal pump at a low flow rate. A filter function is introduced for the RNG k-ε turbulence model in order to consider the compressibility in the cavitating flow. The calculated values of the net positive suction head available (H NPSHa ) agree well with the experimental measurements. The cavities in the impeller are obtained and their evolvement with H NPSHa are revealed.     

Entropy and its application in turbulence modeling

The entropy concept was introduced into the turbulence modeling strategy in the present work. First, the turbulent boundary-layer was described from the point of energy dissipation. Based on the theoretical analysis and direct numerical simulations, the relationship between the entropy increment and viscosity dissipation was systematically investigated. Then, an entropy function fs was proposed to distinguish the turbulent boundary-layer from the external flow. This function is universal, independent of the inflow conditions or any specific turbulence model. With this function, a new version of delayed-detached-eddy simulation method SDES was constructed and verified with the supersonic boundary-layer flow and the cavity-ramped flow. Initial results showed that this method could successfully avoid the modeled stress depletion problem inherited from the original DES method.     

复杂铝型材挤压成形有限体积仿真

根据型材挤压成形属于特大变形、热-力耦合、坯料头部存在变形死区等特点,采用欧拉描述有限体积法(FVM)模拟了3级复杂型材的挤压成形过程,避免了传统刚塑性有限元法无法模拟变形死区,拉格朗日描述有限变形弹塑性有限元法难以重划特大变形网格等难题．模拟过程中,采用欧拉描述FVM仿真坏料流动,采用拉格朗日描述有限元法仿真模具变形和应力,采用显式动态接触算法仿真坏料与模具之间的相互作用．确定了兼顾仿真精度和速度的有限体积网格的最小尺寸和单元最大数量,提高了计算效率;建立了6063-T5铝合金在不同变形温度下的系列热-本构-摩擦模型,获得了金属流动、力能、温度、密度变化分布．为工艺和模具参数优化提供了依据．     

Mach8的平板可压缩湍流边界层直接数值模拟及分析

对来流马赫数等于8,壁温等于10.03倍来流参考温度的平板可压缩湍流边界层做了直接数值模拟,计算涵盖了从层流到转捩以及最终充分发展湍流的全空间演化过程.对湍流的统计特征做了详细的分析,结果表明,在当前的计算工况下,湍流边界层核心区平均速度剖面仍然满足对数率,且卡门常数基本不变;可压缩效应明显增强,由于采用近似恢复温度的等温壁条件,使得近壁区温度较高,导致当地声速增大,使得湍流马赫数绝对值较低,造成内在压缩性效应不强,与经典强雷诺比拟相比,除在数量上产生一些偏差外,强雷诺比拟关系近似成立,且Morkovin假设依然有效;对扩展自相似性和标度率分析表明,对于平板可压缩湍流边界层而言,高马赫数流动使得其适用范围减小;压缩性效应对近壁湍动能,条带结构,涡等值面分布的影响得到分析.     

理查德森数对多孔介质T型管道近壁面处流动与传热的影响

对多孔介质T型方管内冷热流体混合过程的流动和热传递进行大涡模拟,获得了近壁面处湍流混合的瞬时速度和温度,分析了3种理查德森数（Ri=-0.517,-0.093,0.517）对流动与传热的影响。数值结果表明：无论理查德森数Ri为正或负,在热流体区,多孔介质骨架温度较低,而在冷流体区,多孔介质骨架温度较高,说明热量从高温流体通过多孔介质向低温流体传递;最大温度波动发生在近壁面处的上部;由于多孔介质的湍动增强作用,减弱了浮升力对流动的影响,不同理查德森数对流动和传热的影响不够明显。     

修正的PANS模型在云状空化流动计算中的应用评价

将考虑空化流动密度变化的密度函数引入PANS模型,基于实验结果评价了这种修正的PANS模型在云状空化流动计算中的应用.分别采用不同的模型控制参数f_k值,计算了绕Clark-Y型水翼云状空化流动,获得了空穴形态的发展历程、升力系数曲线波动情况,并与实验结果进行了对比.结果表明,修正PANS模型继承了标准PANS模型的特点,即随着f_k值的减小,可以获得流场中越来越丰富的多相湍流流动细节,很好地预测绕翼型云状空化流动的非定常特性;修正PANS模型考虑了空化区域气液混合介质对湍流黏性的影响,可以更全面地描述反向射流作用下空穴的断裂及多尺度空泡团的脱落细节,所得翼型升力系数曲线的波动情况更接近于实验值.      

基于方格网状等离子体激励器的翼型湍流减阻实验

等离子体流动控制技术具有结构简单、响应迅速等特点，已成为流动控制领域的研究热点。为减小飞机的湍流摩擦阻力， 提出了一种基于方格网状等离子体激励器的新型湍流减阻方法，研究了其放电特性与诱导流动特性，并在风洞中获得该激励器减小NACA0012翼型湍流摩擦阻力的参数规律。结果表明，静止条件下，方格网状激励诱导的射流速度与占空比成正比，而随脉冲频率的增大先增加后减小，诱导射流的最大瞬时速度为1.75 m/s。来流速度为15 m/s时，激励能使翼型湍流摩擦阻力减小3.5%。方格网状激励诱导产生的射流使近壁面流体整体抬升，破坏近壁面涡结构，进而抑制湍流生成，实现摩擦减阻。     

典型标准水面船型阻力和黏性流场的计算

本文通过采用PISO(Pressure Implicit with Splitting of Operators)算法求解不可压缩RANS方程,对两类典型标准水面船型阻力和黏性流场问题进行数值计算验证.其中自由面的模拟采用单相Level set方法,湍流模式选取SST k-ω模型.首先以标准数学船型Wigley船为例,进行网格收敛性验证,讨论分析时间步长的选取和壁面网格参数的确定对计算结果的影响,为复杂问题的求解提供经验参考;然后在合适的网格及时间步长下计算不同Froude数(Fr)时Wigley船体受到的总阻力系数、摩擦阻力系数、船体表面波高图等,并分别与模型试验结果和经验公式估算结果对比,结果吻合很好,从而验证了该算法的可靠性.此外,还模拟计算了第二类标准船型DTMB5415的阻力系数和黏性流场,总阻力计算结果与实验结果吻合很好,船体波高图与实验结果相比稍有误差.数值计算结果和分析充分证明本文算法和计算程序的可靠性.     

离心通风机内部流场三维瞬态计算

采用三维时均ＮａｖｉｅｒＳｔｏｋｅｓ方程和ｋε模型分别模拟了离心通风机内部流场的流动;几何拓扑结构使用多块结构以减少计算机内存的浪费;运用贴体坐标系统描述了复杂的几何边界;使用二阶中心差分以离散扩散项和源项;对于控制方程的对流项的离散,为了克服中心差分在网格Ｐｅ数大于２后解的不稳定和迎风差分时不考虑Ｐｅ影响的缺陷,应用混合差分（ＨＤＳ）格式,并对时间参数采用二阶向后差分;求解方程使用Ｓｉｍｐｌｅ算法．通过这些方法和手段,得出了不同时刻叶轮内流道速度、湍动能及耗散能的分布．这为了解旋转对流道内速度分布、壁面摩擦力和紊流结构等影响提供了理论依据,为离心风机设计、节能改造提供了详细的数据和有效的手段     

同位网格有限体积法在变密度浆液流场计算中的应用

为研究自膨胀浆液的扩散机理,基于现代计算流体动力学理论,建立了一种变密度浆液二维流场计算方法.该方法采用结构化同位网格剖分计算区域,利用有限体积法离散浆液流动控制方程,借助动量插值技术构建压力修正方程;采用压力耦合方程组的半隐式算法(SIMPLE)迭代求解动量离散方程和压力、速度修正方程.通过变密度浆液在二维矩形狭槽中自由膨胀算例对该算法进行检验.结果表明,该方法求解正确,并具有较高的计算精度,为建立膨胀性高聚物注浆材料在二维裂隙中流动扩散仿真分析方法奠定了基础.     

封闭方腔自然对流的涡结构和传热特性

为研究封闭方腔自然对流的涡结构和传热特性,对普朗特数,Pr=0.71,方腔高长比A=H/L=1,瑞利数Ra=1.58×109的二维封闭方腔自然对流进行了直接数值模拟. 给出了水平边界层发展阶段和垂直边界层转捩阶段的流场结构以及Nu数分布. 结果表明,压力梯度对方腔的水平速度变化起着决定性的作用;高Ra数下的自然对流在逆压梯度作用下水平方向形成一系列的涡,这些涡使得水平边界层流动形成分层结构;壁面Nu数与速度梯度?v/?y有着密不可分的关系.      

三维水翼初始梢涡空泡数的尺度效应

针对梢涡流场和初始梢涡空泡数的尺度效应问题,利用大涡模拟(LES)湍流模型对三维水翼的梢涡流场流速进行模拟计算;为减少误差,对梢涡流域网格进行了局部加密处理,对未发生空化时梢涡内的轴向速度和切向速度进行计算.结果表明,LES湍流模型的流场流速计算结果与实验值吻合较好.同时,介绍了经典初始梢涡空泡数尺度效应公式的推导过程,并利用数值计算的速度环量和涡核半径修正尺度效应公式.     

侧向支柱对冲压增程弹用进气道性能影响研究

结合分区对接网格技术和二阶精度区域分解算法,对某冲压增程弹丸进气道在不同来流攻角和不同侧向支柱形状工况下的内外复杂流场进行了数值模拟.得到了临界工况下超声速进气道内外粘性流场复杂的波系结构,分析了流场结构特性和激波波系结构.当进气道扩压段有侧向支柱存在时,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大.在三种形状侧向支柱中,采用两端削尖支柱形状的进气道的性能最优.攻角的存在也在一定程度上降低了进气道性能.     

船体绕流计算网格自动生成研究

应用保角变换函数和多项式拟合函数，建立了覆盖船体表面及其外场的数值计 算网格自动生成的数学模型和计算机程序。此项研究的主要贡献为：应用边界层薄层 理论计算结果控制网格的总体扩展，应用边界层速度剖面的指数分布律控制网格的局 部扩展；从而保证了计算网格的有效性。可用于船体势流与粘性流的绕流计算，也可 用于船型的计算机辅助设计。     

展向和流向间距不同的粗糙元壁湍流ADV实验研究

为了揭示粗糙元对壁湍流的影响,采用声学Doppler流速仪(ADV)对布有展向或流向圆柱形粗糙元的壁湍流进行了实验测量。研究了不同Reynolds数下不同间距粗糙元对壁湍流流向平均速度和湍流度的影响。结果表明:在同一Reynolds数下,粗糙元展向间距为7倍直径时的壁湍流流向平均速度最小,湍流度最大。粗糙元流向布置的壁湍流流向平均速度随粗糙元间距的减小而减小,湍流度随粗糙元间距的减小而增大。相比流向粗糙元,展向粗糙元对湍流流动特性影响更大。研究结果对壁湍流流动控制和强化换热的工程应用具有指导意义。     

基于密度分域的混合湍流模型的应用与评价

结合数值计算和实验结果评价了一种基于密度分域的混合湍流模型.实验采用高速录像和LDV测量技术观察了绕Clark-y型水翼云状空化随时间的流场结构变化和速度分布.数值模拟中,基于非定常空化的流动现象,建立了一种基于密度分域的混合湍流模型,并分别应用不同的湍流模型计算了绕水翼非定常云状空化流动.通过与实验结果的对比表明,分域湍流模型能够更好地调整流场内的湍流黏性,精确捕捉到空穴形态和空泡脱落的非定常细节,时均流场速度分布与实验值基本吻合.