浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动.该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递.分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法.流场...     

任意多面体网格上的欧拉方程数值算法

发展并验证了一种新的支持多面体网格的欧拉方程离散算法,采用Fortran 95编写了支持任意网格拓扑及三维复杂外形的通用求解器.在空间离散上,基于径向基函数理论提出了一种新的梯度计算方法,并采用动能通量分裂格式来得到准确并且稳定的二阶精度重构.该方法不依赖于网格几何形状,因此对网格质量不敏感.由于在时间方向采用了点隐松...     

三维贴体坐标系下燃烧室热态流场的大涡模拟

大涡模拟三维贴体坐标系下环形燃烧室火焰简热态紊流瞬态流场。利用椭圆方程方法生成三维贴体网格，计算中采用k方程亚网格尺度模型估算亚网格紊流粘性；亚网格EBU燃烧模型估算化学反应速率；热通量辐射模型估算辐射通量。并在非交错网格系下采用sIMPLE算法和混合差分格式求解离散方程，利用壁面函数处理固壁边界条件。计算结果与实验结果的比较表明，采用大涡模拟方法能更真实反映环形燃烧室火焰简内紊流化学反应流气流结构和燃烧过程。     

带隔热屏加力燃烧室热态流场计算

在任意曲线坐标系下对带有横向波纹隔热屏、外冷却通道和尾喷口的涡喷加力燃烧室的三维热态流场进行数值模拟.利用椭圆型微分方程和区域法生成三维贴体网格, k-ε模型预估紊流特性, EBU-Arre-henius紊流燃烧模型估算化学反应速率, 六通量模型预估辐射通量以考虑其对壁面和气流温度的影响.在非交错网格系统下应用混合差分格式离散控制方程, SIMPLE算法求解离散方程, 数值分析两种进口气流参数分布对加力室热态流场和壁温分布的影响.所得计算值与试验数据基本相符, 表明计算方法合理.      

结构/直角切割网格下时域有限体积法在计算电磁中的应用研究

 采用有限体积法直接求解时域麦克斯韦方程组，研究了二维圆柱、翼型和三维球体、立方体、双球体等典型完全导电目标的电磁散射问题。空间离散使用近似黎曼解构建网格通量，并针对二维结构网格、三维自适应直角切割网格分别给出重构方法，时间方向采用二步龙格库塔法。计算与理论、实验或文献结果符合良好。     

环形燃烧室两相喷雾燃烧的大涡模拟

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对环形燃烧室火焰筒气液两相紊流瞬态反应流进行数值模拟.采用椭圆偏微分方程生成三维贴体网格,计算中所采用的数学度模型有:k方程亚网格尺度模型估算亚网格紊流粘性;亚网格EBU燃烧模型估算化学反应速率;热通量辐射模型估算辐射换热.并在非交错网格体系下气相采用SIMPLE算法和混合差分格式求解,液相采用随机离散模型(Stachasttc Separated Flow,简称SSF),在拉格朗日坐标系下追踪各油珠群沿各自轨道运动、质量损失及能量变化.通过计算结果与实验数据相比较,表明在三维贴体坐标系下对燃烧室火焰筒两相紊流油雾燃烧流场进行大涡模拟,采用欧拉-拉格朗日两相大涡模拟方法能反映两相紊流化学反应流流动及实际燃烧过程.     

隐式多重网格法求解叶轮机械三维跨声速湍流流场

本文提出了一种高效率求解叶轮机械三维跨声速湍流流场的隐工多重网格算法。它以矢通量分裂格式为基础，采用ＦＡＳ型多重网格加速技术。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切.使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格.应用特征恢复技术恢复凹面特征并改善凸面处的网格质量.流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量.加入Spalart-Allmaras一方程湍流模型以获得湍流解.上下对称高斯-赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性.算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算.     

Godunov方法在复杂流场数值模拟中的应用

本文根据Ｖａｎ Ｌｅｅｒ提出的方法给出了求解三维非定常Ｅｕｌｅｒ方程和三维定常Ｅｕｌｅｒ方程的二阶Ｇｏｄｕｎｏｖ方法。方法要点为：在各体积离散单元中,用流动变量的分段线性分布代替一阶方法中的分段常量分布,在确定线性分布的斜率中引入单调性限制条件;在离散单元界面构造Ｒｉｅｍａｎｎ问题,由Ｒｉｅｍａｎｎ间断解确定单元侧面的通量;用预测－修正两步法对离散方程推进求解,用上述方法结合单区网格或多区重迭网格     

带副翼三维机翼绕流的N-S方程解

提出了一种求解带副翼偏转三维机翼绕流的N-S方程计算方法。采用对接分区网格与分区求解算法的结合,有效地求解绕此外形的复杂流动。提出了一种满足通量守恒的内边界耦合条件。数值方法中选用vanLeer分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD性质,采用中心差分格式来离散粘性通量项。数值算例表明该方法是求解带操纵面偏转的机翼绕流的有效方法。     

吸气式飞行器进气道唇口三维激波/激波干扰

对吸气式飞行器进气道唇口处三维曲面激波/弓形激波干扰流场进行数值模拟,利用典型三维气动干扰试验对采用的数值计算方法进行验证。利用拼接网格技术及逆距离加权插值方法获得入口处流场的非守恒变量,作为激波干扰研究的入口边界条件。数值模拟表明,唇口处激波干扰流动的三维效应十分显著,曲面激波与弓形激波产生斜交,尽管唇口前缘半径很小,但Edney提出的6类激波干扰类型可能沿唇口展向方向同时存在;第Ⅲ和Ⅳ类激波/激波干扰的诱导使得唇口热流分布异常严酷;激波相交处形成斜向“伤疤”状局部高热流条带,峰值热流达到参考热流的4~6倍,可能引起唇口结构的局部烧蚀或破坏。      

合排气二元收敛喷管与发动机流量匹配设计方法

发展了一种基于喷管三维黏性数值分析的内外涵气流质量流量匹配设计方法.该方法主要是通过建立以喷管待定几何参数为未知量的非线性无量纲内外涵气流质量流量匹配方程组,然后采用逐步线性迭代法结合喷管三维黏性数值分析对该非线性方程组进行求解,得到满足发动机内外涵流量匹配的喷管设计方案.通过数值算例验证分析,结果表明：①该方法可以快...     

壁压信息法低速大攻角洞壁干扰修正

 本文将壁压信息法用于实壁低速八角形风洞,获得了一个三维模型大攻角实验洞壁干扰修正结果;给出了处理八角形风洞切边效应的方法;研究了等效奇点配置诸因素（数目、展宽、后掠角和攻角）的影响;讨论了不同等效奇点布置对修正结果的唯一性问题;研究了壁压条形状及测压点布置对修正结果的影响。     

超燃冲压发动机三维喷流干扰流场的数值研究

用隐式有限体积法求解三维N－S方程数值研究超燃冲压发动机三维平行、垂直喷流干扰流场。计算中采用了具有当地极值递减（LED）性质的JST和SLIP中心差分格式。为了减缓计算网格数很多时对计算机内存的压力以及提高计算效率，运用了通量守恒的分区计算方法。     

轴对称矢量喷管有效喉道调节方法

基于三维雷诺平均Navier-Stokes方程对轴对称矢量喷管内流场进行了数值模拟,分析了不同落压比下轴对称矢量喷管有效喉道及流量系数随矢量偏转角的变化规律.研究发现:矢量偏转角超过一定值时,轴对称矢量喷管有效喉道位置发生倾斜,有效喉道面积减小,流量系数降低,矢量偏转角越大落压比越低,流量系数降低幅度越大.根据研究结果提出了一种针对矢量偏转状态的轴对称矢量喷管有效喉道调节方法,方法以落压比和矢量偏转角为输入参数,考虑了轴对称矢量喷管几何喉道面积调节前后流量系数的变化.该调节方法能够为发动机控制系统提供更精准的输入,提高控制精度,矢量偏转前后流量相差不超过0.4%,调节时间缩短至少10%,可为推力矢量发动机工作状态调节提供参考.      

圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

基于STL文件和有限元法的在轨航天器关键部件热计算

针对STL(stereo lithography)文件在传递复杂几何实体模型信息方面具有精度高的特点,提出了一种基于STL文件计算复杂结构角系数和外热流的方法.根据单元和节点的拓扑关系识别六面体网格边界单元并进行辐射换热计算.详细阐述了基于STL文件和有限元法进行复杂结构的角系数和外热流的计算方法.研究了导热-辐射耦合计算的有限元方法.最后利用有限元法计算航天器关键部件在轨的三维瞬态温度场.      

基于轴对称比拟的高超声速复杂外形气动热预测方法研究

轴对称比拟法是高超声速飞行器气动热计算的一种有效方法。针对轴对称比拟法应用时公式推导繁琐、变量迭代复杂、计算量大的缺陷,直接在笛卡尔坐标系下采用三维线性方程拟合物面方程及流场变量,推导出相应的流线尺度因子计算公式,扩展了驻点区热流密度的计算方法,提出了一种驻点区下游流线推进格式。以改进的轴对称比拟法为基础,将边界层外无粘流场数值方法与边界层内气动热工程算法进行耦合,发展了一套适用于三维复杂外形飞行器的气动热计算方法。通过对球头钝锥和双椭球算例进行验证,结果表明:方法计算效率较高,适用范围较广,热流计算结果和实验数据吻合良好。     

挤压油膜阻尼器的热平衡分析

利用变分有限元法对挤压油膜阻尼器（简称ＳＦＤ）的热平衡问题进行了分析计算。文中根据雷诺边界条件，通过联立求解ＳＦＤ的雷诺方程、能量方程得到了ＳＦＤ的压力分布、滑油流量等参数。为ＳＦＤ的设计与计算提供了更为可靠的理论依据。