基于RNG k-ε湍流模型钝体绕流的数值模拟

采用RNG k-ε湍流模型对钝体绕流流场进行了数值模拟．计算采用非交错网格的有限体积法求解二维不可压N-S方程．计算结果表明，RNG k-ε湍流模型可以成功地模拟绕钝体的不稳定、非定常和剧烈分离流动。     

定常二维层流运动流函数方程数值解法

在分析各种常见数值解法的基础上,采用中心差分和迎风格式两种方法对二维层流运动流函数方程进行离散,然后给出该方程的求解方法和边界条件。     

二维微槽道内的流动模拟

直接模拟蒙特卡罗方法(DSMC)通过计算机跟踪仿真分子的运动,求解稀薄气体流动规律.选用硬球模型模拟仿真分子间的碰撞,编写DSMC的计算程序,模拟处于滑流和过渡流动领域的二维微槽道内气体流动.结果表明,气体流动在边界附近出现速度滑移和温度跳跃.计算结果和理论值相吻合.     

边界元——有限元耦合方法分析

本文对边界元－有限元耦合方法分析二维弹性力学问题进行了探讨。采用非连续元离散边界积分方程，解决了耦合分析中的自由度约束问题。给出了利用非连续元实现边界元－有限元耦合的具体步骤。数值算例表明了本文方法的有效性。     

流—固耦合问题边界元—有限元耦合方法分析

利用边界元－有限元耦合方法对流－固耦合振动问题进行了分析，假设流体控制方程为Ｌａｐｌａｃｅ方程，利用非连续边界元对流体域进行离散，从而有铲地解决了边界元分析中的“角点效应”问题。固体以平面梁为模型，采用有限元进行离散，对非连续边界元和有限元的耦合问题进行了分析，通过对悬臂梁在一侧受液体作用时的瞬态响应分析的数值解同解析解的比较，表明了本文所给方法的有效性，同时为利用边界元－有限元耦合方法对流－固耦     

对流扩散方程的半拉格朗日有限差分法

针对二维对流扩散方程提出了半拉格朗日Crank-Nicolson有限差分格式.化对流扩散方程成拉格朗日形式,然后进行时间和空间离散.时间的离散采用Crank-Nicolson格式,空间离散采用有限差分格式.数值实验表明,新建格式是稳定和收敛的.     

二维非连续边界元分析积分计算的精确表达式

以二维弹性力学问题为研究对象，采用线性非连续元离散边界积分方程，给出了系数矩阵计算的精确表达式，对二维弹性力学问题进行了数值计算，对非连续边界元配位点对计算结果精度的影响进行了讨论，结果表明准奇异积分计算是配位点影响计算结构精度的主要因素。     

考虑黏性及压缩效应的旋涡运动特性

采用扩散速度法处理黏性作用项的离散涡法，数值模拟了旋涡在黏性不可压流场中的运动；推导了可压缩流中的旋涡运动方程，对不可压缩流中的Biot-Savart方程进行了可压缩修正，进而用离散涡方法求解出非定常、不稳定、可压缩流场中的旋涡运动。分析了胀量及黏性扩散作用对旋涡运动的影响，通过大量长时间行为的数值试验得出了一系列结果。结果表明：为了真实地反映旋涡在实际流体中的运动规律，必须考虑黏性及压缩性的影响。     

方管中颗粒随黏弹性流体迁移的力学特性

为研究黏弹性流体中颗粒迁移的力学成因,对方管Poiseuille流动中悬浮颗粒的迁移进行数值研究。选用Giesekus模型的黏弹性流体和球形刚性颗粒作为研究对象,忽略流体的惯性效应以研究流体弹性效应和剪切变稀效应对颗粒受力的影响。控制方程采用有限元法和伽辽金 最小二乘法(GLS)方法进行求解,并取得了较好的收敛性。颗粒运动模型采用基于“相对运动模型的准定常算法”,对通道中颗粒的横向升力的分布特征进行研究。研究结果表明：颗粒在不同横向位置上所受的横向升力决定了颗粒在实际流动中的迁移方向,且受到流体弹性和剪切变稀的影响。     

基于有限体积法的离散裂缝模型两相流动模拟

正确认识裂缝性油藏流体流动规律是其高效开发的基础，离散裂缝模型因其对裂缝予以显式降维处理，能够大大提高计算效率而逐渐成为当前研究的热点。有限体积法在控制体单元上具有明确的物理意义，能够严格保证局部质量守恒，适用于流场性质不连续的情况。基于Delaunay 三角网的对偶网格构建有限体积控制单元，推导了离散裂缝模型油水两相流动控制方程的有限体积数值计算格式，实现了离散裂缝模型油水两相数值模拟，通过算例验证了该方法的正确性与高效性。研究表明：该方法能够严格保证局部质量守恒，计算过程中收敛性较好，适用于离散裂缝油藏的数值模拟，与单孔隙介质模型计算结果具有很好的一致性，计算效率提高2.5∼3.0 倍。     

多级离心压缩机级间静止部件气动特性与旋涡结构的数值研究

通过CFD技术对某型多级离心压缩机内的流动进行了数值模拟,重点研究了多级离心压缩机级间装置如回流器、两类弯道内气动特点及其旋涡结构的特征,旨在揭示气流经弯道与回流器后气动参数的变化规律.结果表明,不同位置的弯道对流动的影响是截然不同的.由于这些静止部件对流体的扰动,使次级动叶的进口来流工况点偏离设计值,易于造成整机气动性能的下降.因此,多级压缩机的动叶设计应对来流的不均匀性加以充分考虑.     

用SMAC方法对二维非定常流动的数值模拟

使用ＳＭＡＣ算法,从非定常的Ｎ－Ｓ方程出发,对二维非定常流场进行数值模拟。并在非交错网格上,针对其压力算法的不足,通过求解压力Ｐｉｏｉｓｓｏｎ方程来求解压力。计算与实验的结果比较,显示了本文的算法较好地模拟了二维低Ｒｅ圆柱与方形钝体的绕流。     

不可压粘性流动N-S耦合方程的有限元解法—I应用

应用有限元方法求解不可压粘性流动Ｎ－Ｓ耦合方程，将此方法用于驱动方腔和突扩管道内的层流分析，结果表明此方法具有良好的椭圆型性质，能预测分离流动．用该方法结合两层湍流代数模型，对压气机叶栅内不可压流动进行模拟分析，与实验结果的比较证明方法是有效的．     

圆柱绕流远场涡对的数值模拟

基于离散涡方法求得非定常、不稳定流场,数值模拟了三种不同时刻高雷诺数下圆柱绕流结构的发展,从流谱图、等涡量线图和涡谱图可以清晰地看出从近场的初生卡门涡街,过渡到远场的二次涡街的过程,计算结果发现：远场离散涡有形成二个涡的涡对及三个涡的涡对的趋势,计算结果说明了流体运动中涡对结构的本质：由于来流是均匀的,没有加入任何拔动,当流体流过钝体时产生具有剧烈分离的不稳定流动,因此在远场形成的二次涡对及卡门涡     

不可压粘性流动N－S耦合方程的有限元解法──I理论

本文对不可压粘性流动非定常Ｎ－Ｓ方程组在耦合条件下应用有限元方法直接求解，并对在高雷诺数下如何提高计算的精度、稳定性以及收敛速度等方面进行了讨论，从而建立起不可压粘性流动Ｎ－Ｓ耦合方程的有限元解法．     

河道中钝体周围的床面冲刷试验研究

对泄洪闸闸墩、淹没薄壁长立方体以及圆柱形桥墩周围的河床床面冲刷进行了较为详细的试验。以试验所获得的大量测试数据、录像资料及照片为基础,分析了钝体周围的水流速度场、水流旋涡场、并重点分析了钝体周围的床面冲刷与钝体周围水流旋涡结构的关系,得出了关于河道中钝体周围床面冲刷的一些规律性结论。     

布置错位通道的微混合器内流体热混合特性

基于常规平面微混合器内流体混合特性的研究,对布置错位通道的微混合器进行了三维数值模拟,分析了指定Reynolds数范围内微通道内流体的热混合特性。计算模型采用有限体积法离散、SIMPLEC算法进行层流计算。结果表明,温度在指定范围变化时,其对混合性能的扰动作用较小。此外,由于错位通道结构所产生的漩涡在促进混合的同时降低了流体工质的温度。     

轴流式叶栅内部流动变工况性能的计算分析

应用有限元方法求解不可压黏性流动Ｎ－Ｓ耦合方程，对ＮＡＣＡ叶型的叶栅在不同冲角的条件下，叶栅内部的流动作了较为完整的计算分析。结果包括叶片表面的压力分布、叶片表面法向的速度分布以及表面摩擦系数。在计算中，采用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型和壁面函数以提高计算的速度。结果表明：冲角不仅对叶片表面的压力分布，而且对叶片附近的流动状态，特别是在叶栅流道的后半部，有很大的影响。     

向心透平蜗壳轴向对称性影响的数值模拟研究

基于三维Naiver-Stokes方程组，采用分块结构化网格，用数值计算的方法模拟了矩形截面轴向对称和非对称蜗壳内气体的流动．湍流模型采用Baldwin—Lomax代数模型．模拟的结果表明，对称蜗壳和非对称蜗壳在出口气流特征上基本一致，非对称蜗壳气流的周向均匀性要优于对称蜗壳；非对称涡壳的后1／4流程出现了横向压力梯度，这个梯度容易引起二次流动损失．     

新型涡流空化器的实验研究

目前利用空化释放出能量是充分利用能量的一种新方法,为此设计了一种新型的涡流空化器,定量分析了空化 引发?OH（羟自由基）产量与空化强度的关系。文中以亚甲基蓝作为羟自由基的捕捉剂,使用紫外一可见分光光度计间 接测定了水力空化产生的?OH产量,并利用该方法对涡流空化器空化效果进行了实验研究。分析了旋流腔直径、空化 器出液管直径和空化器进液管流速等工艺参数对自由基产量的影响。力求寻出有利于空化发生以及提高空化作用强度 的适宜空化器的设计工艺。结果表明：进液管流速越大,出液管直径越小,旋流腔直径越小,?OH摩尔浓度越大。与孔 板水力空化器对比,涡流空化器单位压降的?OH的摩尔涨度是孔板水力空化器的13. 99 - 68. 19倍,说明涡流空化器能 量利用率更高,空化效果更好。