非正交曲线坐标湍流族流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴你坐标系统，用坐标变换的方法处理锥形燃烧室曲线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。亲流模型采用工程上常用的ｋ－ε模型。还运用贴体坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提供了理论依据。     

区域法数值模拟短突扩压器流场

在一般曲线坐标下对环形燃烧室突扩压器和其内外环通道流场进行了数值模拟。紊流特征有用工程上常用的κ－ε双方程模型来描述，使用壁面函数法处理曲线壁面边界，采用非交错网格系统和ＳＩＭＰＬＥ算法求解各守恒方程。为了计算方便，利用区域法对复杂的流道进行分块处理，计算所得的速度和压力分布较为合理，表明本方法是可行的。     

气氢亚燃冲压燃烧室三维流场计算

采用数值计算方法模拟气氢冲压模型燃烧室流场。文中采用к－ε模型描述紊流特性，Ｍａｇｎｕｓｓｅｎ的ＥＤＣ模型估算化学反应速率，六通量辐射模型考虑热辐射对燃烧过程的影响，采用非均匀交错网格系统，控制容积法和混合差分格式，离散控制方程，并应用ＳＩＭＰＬＥ算法求解。在冷热态流场计算基础上，研究了不同来流空气和燃料氢各气流参数以及喷氢孔直径，孔数和孔间距对燃烧室出口温度分布，温升效率以及总压恢复系数的影响，     

大涡模拟在加力燃烧室冷态流场中的应用

对带V型稳定器模型加力燃烧室二维冷态流动进行了大涡模拟。采用了代数和k方程两种亚网格尺度模型，并用SIMPLE算法和混合差分格式求解离散方程，用壁面函数处理固壁边界条件，得到了模型加力燃烧室二维瞬态流场的模拟结果。通过对不同亚网格尺度模型和不同入口速度下的瞬态流场的比较，揭示了不同亚网格尺度模型和不同入口速度分布对带V型稳定器燃烧室二维冷态瞬态流场的影响。     

双旋流燃烧室两相喷雾试验和数值研究

采用粒子场脉冲激光全息技术对航空发动机燃烧室中的雾化场进行了测量，得到了燃烧室中燃油液滴直径的空间分布，从而对燃烧室中的雾化过程进行了研究。自主开发完成了适用于航空发动机燃烧室的三维两相数值计算平台，建立了首次雾化模型和二次雾化模型。基于LISA模型和KH-RT模型，对燃烧室中的首次雾化过程和二次雾化过程进行了数值模拟，得到了燃烧室中液雾的空间分布。通过将计算结果与试验结果进行对比，显示开发完成的雾化模型能很好的模拟高温高压，强旋流条件下航空发动机燃烧室的整个喷雾雾化过程。     

波瓣喷管混合流动的数值计算

借助Ｎ－Ｓ方程对波瓣喷管混合装置的流场进行了数值分析。文中采用三维贴体曲线坐标,非交错网格法和代数湍流模型,对波瓣喷管内外两侧的气流通道及波瓣喷管下游的混合段波场采用分块计算方法进行了数值计算。     

煤粉锅炉旋流燃烧器的气动热力性能数值分析

采用数值分析法对旋流燃烧器气动热力性能进行预估，利用代数应用力模型描述紊流特性，紊流燃烧模型为k-ε-g模型，固体壁面采用参数投影法，热辐射对燃烧室内传热过程影响是用修正热通量模型来考虑。     

双旋流燃烧室两相喷雾试验和数值研究（英文）

采用粒子场脉冲激光全息技术对航空发动机燃烧室中的雾化场进行了测量,得到了燃烧室中燃油液滴直径的空间分布,从而对燃烧室中的雾化过程进行了研究。自主开发完成了适用于航空发动机燃烧室的三维两相数值计算平台,建立了首次雾化模型和二次雾化模型。基于LISA模型和KH-RT模型,对燃烧室中的首次雾化过程和二次雾化过程进行了数值模拟,得到了燃烧室中液雾的空间分布。通过将计算结果与试验结果进行对比,显示开发完成的雾化模型能很好的模拟高温高压,强旋流条件下航空发动机燃烧室的整个喷雾雾化过程。     

高超声速马赫反射

通过实验和数值计算模拟了激波（马赫数７．６）绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项、粘性项和化学源项分别采用了迎风ＴＶＤ格式、中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

进口速度分布对短突扩压器性能的影响

在任意曲线坐标系下，对环型燃烧室的三维流场进行了数值计算，研究了四种不同的进口速度分布对扩压器性能的影响。计算结果表明，进口速度大小和方向都对扩压器性能有重要影响，进口速度分布均匀以及方向贴近扩压器壁面，有利于提高扩压器的性能。计算中采用标准ｋ－ε双方程紊流模型，采用控制容积法进行离散，在非交错网格体系下用ＳＩＭＰＬＥ法求解     

主燃孔对双旋流燃烧室流场的影响

通过试验和数值模拟的方法研究了主燃孔几何参数与布局对双旋流燃烧室内三维冷热态流场的影响.用粒子图像测速仪(PIV)对不同工况的冷热态流场进行测量,用多点温度耙对燃烧室出口温度分布进行测量,用多点取样管和气体分析仪对燃烧室出口气体成分进行分析.用商业软件Fluent对不同主燃孔时的燃烧室冷热态流场进行数值模拟.计算结果与试验数据相当吻合,表明在相同的进口条件下,燃烧室主燃孔的几何参数和布局变化对燃烧室的冷热态流场、出口温度及燃气成分分布有较大的影响.     

采用CARS试验技术与UFPV数值方法研究航空发动机燃烧室

在自主开发的软件平台上，采用基于URANS的方法计算航空发动机燃烧室的三维两相燃烧流动，考虑了液态燃油从液膜-液滴-燃气-燃烧的完整物理化学过程。其中，颗粒相采用LISA一次破碎模型，KH-RT二次破碎模型和标准的蒸发模型，湍流燃烧模型采用可以考虑非稳态燃烧特性的非稳态火焰面/反应进度变量方法，得到了航空发动机燃烧室中温度、组分浓度和燃油液滴的颗粒直径分布规律。同时，采用CARS光学手段测量燃烧室主燃区的温度分布，并将数值计算结果与光学试验测量值进行比较，数值计算结果和试验值吻合较好，数值计算误差小于7.3%。说明了本文的数值计算方法和UFPV方法在计算航空发动机燃烧室的两相燃烧流动时具有较高的精度。     

跨洞庭湖大桥主梁节段定常气动力研究

采用风洞实验与数值模拟相结合的手段，对跨洞庭湖大桥主梁节段的定常气动力特性进行了研究；流场数值模拟是采用ｋ－ε湍流模型封闭的雷诺时均Ｎ－Ｓ方程，求解方法是建立在交错网格基础上的压力－速度场迭代求解的ＳＩＭＰＬＥ算法；实验是在国防科技大学ＫＤ－０３ 低速风洞进行；研究表明实验与计算结果较为一致     

采用CARS试验技术与UFPV数值方法研究航空发动机燃烧室（英文）

在自主开发的软件平台上,采用基于URANS的方法计算航空发动机燃烧室的三维两相燃烧流动,考虑了液态燃油从液膜-液滴-燃气-燃烧的完整物理化学过程。其中,颗粒相采用LISA一次破碎模型,KH-RT二次破碎模型和标准的蒸发模型,湍流燃烧模型采用可以考虑非稳态燃烧特性的非稳态火焰面/反应进度变量方法,得到了航空发动机燃烧室中温度、组分浓度和燃油液滴的颗粒直径分布规律。同时,采用CARS光学手段测量燃烧室主燃区的温度分布,并将数值计算结果与光学试验测量值进行比较,数值计算结果和试验值吻合较好,数值计算误差小于7.3%。说明了本文的数值计算方法和UFPV方法在计算航空发动机燃烧室的两相燃烧流动时具有较高的精度。     

模型燃烧室湍流亚网格尺度模型

分别采用Smagorinsky-Lilly和k方程亚网格尺度模型,在任意曲线坐标系下对模型燃烧室湍流流场进行大涡模拟.分别用QUICK和二阶隐式差分格式离散控制方程及其时间项,并用SIMPLE算法进行求解.计算结果与PIV测量数据比较表明:大涡模拟方法能较好模拟旋流器后面回流区的瞬态变化,两种亚网格尺度模型都能较好模拟湍流流动,但k方程亚网格尺度模型稍优于Smagorinsky-Lilly亚网格尺度模型,因此采用k方程亚网格模型更适用于大涡模拟实际燃烧室复杂湍流流场.     

基于Greenwood叶间燃烧室实验模型的数值计算

为分析新型叶间燃烧室的燃烧性能，以叶间燃烧室试验模型为研究对象，用Realizable k-ε湍流模型、非预混PDF 燃烧模型，数值模拟了燃烧流场特性。燃油颗粒大部分在环形腔顶部区域蒸发混合，火焰在环形腔底部连通，形成旋流燃烧，涡轮径向斜槽内出现温度较高的火焰区。涡轮径向斜槽结构以及主流的卷吸作用对燃烧产物与一次气流的掺混存在影响。贫油燃烧方式工况燃烧性能较优，富油燃烧方式工况的径向平均出口温度线型较好。计算结果与相同操作条件下的试验结果符合较好。     

一个非常规前体机身的流动显示研究

描述了有类似Ｅｒｉｃｋｓｏｎ前体的非常规机身的初步流动显示研究结果，包括测力、表面油流和激光片光流动显示结果和用层流Ｎ－Ｓ方程进行数值模拟的结果。计算和试验的参数范围：α＝０°～５０°，β＝０°～２０°，虽然计算与试验所用的外形在后部有一些不同，但是两者在涡的位置方面显示了良好的一致性。同时研究也表明，大攻角的流动特性可以通过改变机身前体形状进行控制。通过研究还表明，这类前体在改善大攻角横侧方向安定性方面具有很大的潜力。     

矩形喷管内外喷流流场的数值模拟

基于 Ｆｒａｖｅ 平均的 Ｎ Ｓ方程和 Ｂ／ Ｌ 湍流模型，采用了 Ｊａｍ ｅｓｏｎ 格式和 Ｍ Ｕ Ｓ Ｃ Ｌ 格式，对外流在亚声、超声和零 Ｍ ａｃｈ 数三种状态时矩形喷管的喷流流场进行了数值模拟，探讨了内外喷流及喷管内的流动特征，并与实验数据进行了比较，计算结果与实验结果吻合较好，验证了计算所使用的 Ｎ Ａ Ｐ Ａ 软件的可靠性，从而为用计算流体力学的方法模拟具有两股来流的三维湍流流场提供了一种有效的手段。     

三维曲线坐标系N—S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错风格设置，该算法双ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的了出现，在压力修正方程中对压力梯度的离散附加一项１－δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

双重互易边界元法在求解N—S方程中的应用

提出了用双互易边界元法求解Ｎ－Ｓ方程的新思路。用Ｐｅａｃｅｍａｎ－Ｒａｃｈｆｏｒｄ算子分理解将时间相依的Ｎ－Ｓ方程分解为线性和非线性子问题，线性问题用共轭梯度法解除压强－速度的耦合；非线性问题进行局部线化。对所得的Ｐａｓｓｉｏｎ方程及相关类型方程采用双重互易的边界元法求解，消除了传统Ｎ－Ｓ方程边界元解法的区域积分问题。