飞机构型Euler绕流结构网格方法计算的可信度

采用三种方法验证飞机构型Ｅｕｌｅｒ绕流结构网格方法计算。第一,同一网格,采用不同的方法或程序进行计算,采用Ｊａｍｅｓｏｎ中心差分有限体积方法和Ｓｔｅｇｅｒ－Ｗａｒｍｉｎｇ流量分裂迎风格式有限体积方法,为计算提供容错手段;第二,网格细化和网格点重新分布,并通过网格收敛指标进行评估;第三,提出新的多块拼接网格拓朴,采用不同的网格拓扑进行计算。     

二维非结构网格生成及Euler方程计算的方法研究

用离翼型表面最小作为阵面推进法中的参数选择依据。生成二维问题的非结构网格。这种方法戏了传统的背景网格观念,直接提供网格生成过程中所需的背景信息。在求解Ｅｌｕｅｒ方程时,用格心格式的有限体积法作空间离散,用四步Ｒｕｎｇｅ－Ｋｕｔｔａ作时间推进,采用不同的加速收敛措施获得定常流动。提出了两这界条件的构造办法,并 不同边界条件对结果的影响。     

二维非结构自适应多重网格的Euler方程解

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中。利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡。为提高计算效率,采用了自适应多重网格法。在生成非结构网格时,采用了一种新颖的数据结构——数组链接表。最后给出了NACA0012翼型和带10%鼓包的管道流动的几个Euler方程的算例,并进行了分析比较。     

非结构混合网格上的NS方程求解方法

提出了一套较为通用的,完全自动化的非结构混合网格生成方法.在物面粘性作用区,采用一种改进的推进层方法生成三棱柱形和金字塔形网格;在其他流动区域采用阵面推进方法生成四面体网格.采用一种改进精度的格心有限体积法对三维NS方程进行了求解,在加速收敛措施方面,提出了一种新的当地时间步长取定方法来减小质量较差的网格单元对流场计算稳定性和收敛速度的不利影响.以M6机翼和DLR/F4翼身组合体外形的粘性流场作为数值算例,验证了上述网格生成和流场求解方法的正确性和实用性.     

一种非结构/结构多层混合网格方法及其应用

对于粘性绕流的数值模拟,在自适应直角网格基础上,结合三角形非结构网格和结构化网格,利用其各自的优势和特点,提出一种生成混合杂交网格的思路和方法.在物面附近生成适合粘性流计算的大长宽比结构化网格,在远场分布自适应直角网格,快速离散计算空间.对于复杂的多体问题,采用三角形网格来连接各体网格,并运用网格合并的方法,保证各网格之间的光滑过渡与连接,提高网格质量.针对一些二维、三维外形的绕流问题,在上述网格基础上,采用B-L代数湍流模型和中心有限体积法,完成Navier-Stokes和Euler方程数值模拟的对比计算,结果表明网格生成和流场计算是正确的.     

多分量流计算的高分辨KFVS有限体积方法

论及高分辨分子动力学通向量分裂（ＫＦＶＳ）有限体积方法的推广。在方法中提出了适当修改Ｍａｘｗｅｌｌ平衡分布用以修复Ｅｕｌｅｒ方程。基于熟知的Ｅｕｌｅｒ方程与Ｂｏｌｔｚｍａｎｎ方程的关系提出了一类求解多分量Ｅｕｌｅｒ方程的高分辨分子动力爱向量分裂（ＫＦＶＳ）有限体积方法,应用该方法不需要求解任何Ｒｉｅｍａｎｎ问题或求附加的非守恒压力方程也需要任何非守恒修正。数值计算表明,数值解在物质界面附近无振荡,     

用非结构网格与欧拉方程计算复杂区域的二维流动

提出用Delaunay三角化方法生成非结构网格的一种过程。所生成的网格可用于复杂多连通域内的可压流计算。采用Euler方程和格心有限体积法,研制出程序,给出了算例。     

飞机空速校准用静压测量装置轮廓结构绕流影响分析

介绍了针对一种面向飞机空速校准的外置式静压测量方法,为了减小测量装置轮廓结构绕流对静压测量的影响,从而提高静压测量精度,应用FLUENT对装置轮廓结构进行绕流影响分析,得到绕流影响最小的轮廓结构类型,并采用正交试验方法对收缩角度β、过渡长度L、最大直径ΦD、最小直径Φd共4个参数进行优化;仿真结果表明,测量装置轮廓采用圆弧型收缩过渡结构绕流影响最小,轮廓最大直径和过渡长度对静压测量绕流影响较为显著;最优参数组合为:β=33°,L=40 mm,ΦD=50 mm,Φd=22 mm,优化后的结构对静压管的绕流影响减小了37.5%,提高了静压测量精度。     

无结构网格有限体积方法及在热对流中的应用研究

首次将无结构三角网格的有限体积方法和压强连接半隐式算法相结合,用于求解非平行壁管道中的热对流问题。并由此分析了化学汽相淀积薄膜生长的均匀性问题。计算结果对于分析一类管道中热和动量输运现象均有普遍指导意义。     

三维多介质欧拉方法混合网格的有限体积计算格式

在采用Youngs界面重构技术的基础上，对三维欧拉方法混合网格的计算格式进行了研究。运用Youngs技术确定界面后，混合网格内每一部分物质一般不再是正规的六面体结构，可能是非规则的四面体、五面体、六面体或七面体。本文采用对非规则网格适应性很强的有限体积法对每一部分分别进行计算，给出了混合网格内每种物质的压力、速度、能量等的计算公式，比较有效地解决了混合网格的计算问题。     

非结构网格上温度扩散方程的能流计算方法

讨论非结构网格上温度扩散方程的能流计算方法.应用有限点方法(Finite Point Method,简称FPM)导出基于有限点两点公式和三点公式的能流计算公式,该公式适用于任意多边形及非匹配网格等非结构网格;给出网格角点温度新的计算公式.数值试验表明:基于两点公式的离散解和基于三点公式的离散解均具有平方阶的收敛速度;基于三点公式的离散解的精度总优于基于两点公式的离散解.     

求解可压缩流的高精度非结构网格WENO有限体积法

提出-种基于最小二乘重构和WENO限制器的非结构网格高精度有限体积方法.用中心网格的某些邻居网格建立重构多项式,给出-定的原则搜索和存储足够多的邻居网格以建立重构多项式,采用最小二乘法求解重构多项式的系数.用-种通用的方法控制重构邻居个数,以减少存储和计算,采用WENO限制器和旋转Riemann求解器以达到统-的高精度并且抑制守恒律方程求解中的非物理振荡.为检验上述算法,以基于节点的梯度重构,Bath and Jesperson限制器的二阶算法为基准,给出三阶和四阶格式与二阶格式以及高阶格式若干经典算例计算结果的对比和分析.     

基于非结构四边形网格的WENO格式

基于非结构四边形网格发展求解双曲守恒律的三阶加权基本无振荡（WENO）格式.针对任意非结构四边形网格选取重构模板,并给出基于线性多项式的三阶线性重构.但对于一般的非结构四边形网格,会出现非常大的线性权和负权,使得非线性重构的WENO格式对光滑问题也不稳定.本文给出一个处理非常大的线性权的优化重构方法,对优化后得到的负线性权采用分裂方法进行处理.对于非线性权,提出一种考虑局部网格和物理量间断的新光滑度量因子.采用优化重构方法和新的非线性权,当前的三阶WENO格式在质量很差的网格上也具有很好的稳定性.理论的三阶精度在数值精度测试算例中得到验证,同时一范数和无穷范数的误差绝对值不依赖于网格质量;具有强间断的数值结果证明了当前格式的有效性.     

一类自适应运动网格的ALE方法

从积分形式的二维Lagrange流体力学方程组出发,使用ENO高阶插值多项式,推广了四边形结构网格下的一阶有限体积格式,构造一类结构网格下的高精度有限体积格式.结合有效的守恒重映方法,发展一类高精度的ALE方法,并结合自适应运动网格技术,进行ALE方法的数值模拟,得到预期的效果.     

计算圆柱绕流拉力系数的自适应迎风有限元方法

研究圆柱绕流中拉力系数的自适应计算问题.对二维定常粘性不可压Navier-Stokes方程使用迎风有限元格式,通过引入对偶问题,得到圆柱拉力系数的加权后验误差估计器,实现了自适应网格加密.数值实验表明格式的可行性以及估计器的稳健性.     

在冷气干扰下叶栅绕流的控制方程及求解方法

1控制方程高温涡轮的冷却空气量随着燃气轮机性能的提高而在不断增加，因而在数值模拟气冷涡轮流场时则须考虑喷射冷气对燃气主流的影响。本文采用文献［1］的“质量源”模型，推导出具有统一形式的适用于叶栅三维和周向平均混合流动的控制方程。在绝热和忽略粘性应力作功率条件下，将文献出中的通用形式的控制方程在柱坐标系中展开，然后再经坐标变换，将其转换到任意曲线坐标系《，刀J中，得到了在冷气干扰下的叶栅三维绕流的双曲守恒型控制方程（用矩阵形式表示）：aU／＆＋aE／a＋OF／的十OG／0（一H（1）式中仅三项含有冷气参数，…     

基于树网格的格子Boltzmann方法以及曲线边界的处理

将格子玻尔兹曼方法(LBM)和BFECC技术相结合,在四叉树网格上实现,并用于二维不可压缩流动的求解.二维四叉树(三维八叉树)网格是一种非均匀网格,利用树结构进行数据存储,具有非常高的检索效率.BFECC是一种数值技巧,可以把求解输运方程的奇数阶格式通过误差修正提高一阶精度.基于四叉树网格的LBM,引入BFECC后,能有效减小输运过程中非均匀网格插值产生的误差,在总体上实现二阶精度,和经典的LBM方法相当.最后给出二维顶盖驱动方腔流动算例和二维圆柱绕流算例,从中可以看出格式的有效性.     

考虑密度与温度耦合关系的非结构网格SIMPLEC算法

提出一种新的密度压力速度和温度耦合的算法,使用静焓方程构造温度修正量与压力修正量之间的关系,类似Rhie-Chow方法用动量方程插值,网格面上的温度压力修正量插值采用能量方程的形式.CGNS API(CFD General Notation System Application Programming Interface)作为非结构化网格求解器的前处理和后处理,基于网格的FVM(Finite Volume Method)作为偏微分方程求解方法.燃烧模拟采用Arrinius-Eddy Dissipation模型.最后给出两个算例作为说明.     

二维不可压流函数N-S方程的多重网格方法

通过对二维不可压缩N-S方程的涡量-流函数方程组消去涡量而得到仅以流函数为求解变量的控制方程,从而 使不可压N-S方程的求解个数减到最少。求解方法采用本文提出的二阶精度的九节点紧致差分格式,因此无须对靠近边 界的网格点作特殊处理。为了加快迭代收敛速度,采用多重网格加速技术。数值实验结果验证了方法的精确性和可靠性。     

含运动物体流场计算的自适应非结构二维网格生成方法

采用自动插点的Delaunay方法和局部网格重新生成方法,对含运动边界自适应非结构二维网格生成方法进行了系统研究,实现了带有非定常运动物体流场的数值计算,并进行了数值实验.