二维及三维贴体网格的一种生成方法

本文讨论了二维及三维情况下应用Ｌａｐｌａｃｅ方程生成贴体网格的方法，并应用ＡＦ１格式求解方程以大大加速迭代收敛的速度，对于二维及三维情况，分别讨论了相应用Ｄｉｒｉｃｈｌｅｔ边界问题解法。对于二维情况，还讨论了翼型的Ｏ－、Ｃ－和Ｈ－型三种典型的网格生成方法。着重讨论了有关割缝上的网格点的处理方法，以保证网格的光滑性。三维情况下，以简化航天飞机外形为例，讨论了网格的生成方法，尤其着重讨论了边界点的处理     

壁面冷却的三股混合气流的数值计算

对双层壁扩压器与波瓣喷管组合的三股气流的引射混合进行了流场温度场的数值计算。计算域划分为各子域 ,通过求解拉普拉斯方程 ,分别生成三维贴体曲线坐标网格 ,边界网格加密且波瓣处正交 ,各子域组合成整个网格。为避免因波瓣造成的网格强烈的非正交而引起解的发散 ,本文采用了同位网格 SIMPLEC计算方法和Chen-Kim修正的 k-ε湍流模型。固体区域采用大粘性解决流固耦合 ,还结合了线性欠松弛和拟瞬态欠松弛。结果证明 :双层壁间有外界冷气流被引射进入 ,形成壁面的冷却气流 ,相对单层壁扩压器 ,双层壁扩压器的壁面温度明显降低。对于小扩压角 ( 0～ 5°) ,随着扩压角的增大 ,壁面平均温度降低 ;但当扩压角较大时 ( 5～ 1 0°) ,壁面平均温度则增加     

三角网格模型的数据分块算法

对测量获得的三角网格模型进行分块。首先根据三角片中心点的高斯曲率和平均曲率的符号将三角片标示成8种类型。然后任取一个三角片作为生长点，沿着边界向外扩展“生长”，将具有相同类型的三角片组成一个连通的块。最后采用两种方法对块实施优化：根据平均曲率差值最小的原则，将小块融合到邻近块中．用两点间的最短路径来光滑边界。本方法中的三角网格和块的拓扑均可以是任意的，分块结果符合零件的特征构造。     

双流混合边界的正交贴体网格生成

通过数值实验表明，用泊松方程变换形式生成的贴体坐标网格，其边界点的β值与源项呈连续单调的函数关系。根据这种关系，可以用方程迭代求根方法中的增值寻根法和对分法，使网格边界正交化。由薄壁结构相隔的两种流体混合边界正交网格生成的实例表明，这种方法对于生成适应于红外抑制器全流场计算的网格是很有效的。     

多尺度湍流模式在压缩拐角绕流中的应用

 应用4 个多尺度涡粘性湍流模式对一个二维压缩拐角流动进行了数值模拟,通过与实验结果以及两个基准模式(一个单尺度k-ε两方程模式和一个零方程模式) 的计算结果的比较,结合压缩性对湍能输运的作用、湍流的非平衡性对多尺度湍流模式在拐角非平衡湍流中的应用进行了评估。     

二维有限元网格全自动生成方法

主要研究二维有限元网格全自动生成方法。提出一种新的节点生成方法 ,通过实例验证该方法可靠。采用三角形形状因子评定三角形网格质量 ,可减少三角化的时间 ,优化三角网格质量。     

含内嵌激波流场的自适应网格求解

为了精确求解含内嵌激波的流场，应采用自适应网格技术。本文用Ｂｒａｃｋｂｉｌｌ和Ｓａｌｔｚｍａｎ的基于变分原理的自适应网格生成方法得到了在“人造斜激波”附近自动加密的网格，然后将该方法用于圆柱体底部流场的ＮＳ方程的求解。计算结果表明，用自适应网格求解比用固定网格求解捕捉的激波宽度较窄，质量较高。     

基于IPDG方法的嵌套网格技术

在内罚间断伽辽金（IPDG）方法框架内应用嵌套网格技术求解了二维可压缩Navier-Stokes方程。通过把黏性通量作为辅助变量实现了Navier-Stokes方程的降阶,继而用间断有限元方法进行离散得到空间半离散方程。时间推进采用Newton-Krylov隐式方法。利用库埃特流动精确解验证了该方法的精度。在此基础上,对包括有黏NACA0012翼型绕流、定常和非定常的圆柱绕流在内的若干典型测试算例开展数值模拟,进一步验证了该方法的鲁棒性和可靠性。在嵌套网格中应用了h型网格自适应技术用于提高激波分辨率,对NACA0012翼型无黏跨声速绕流开展数值模拟,结果表明：自适应之后的单元数相比无自适应时只增加了8.4%,但是激波分辨率得到了显著提高,激波附近的计算结果也与实验值符合得更好,从而验证了该方法的有效性。     

四面体网格边界恢复改进算法

Delaunay三角化是生成四面体网格的主流方法,而通过该方法生成的初始网格无法保留所有边界约束,因此边界恢复是Delaunay四面体网格生成的必要步骤。边界恢复的难点在于确定辅助点（Steiner点）的数量和位置,以及降低Steiner点引发的负面效应。为此,本文提出了一种改进的四面体网格边界恢复方法,旨在提高约束边界恢复成功率,并减少加入的Steiner点数量。该方法总体流程由“两层+两轮”迭代构成,流程充分利用了约束边界恢复方法和体内Steiner点插入方法的优势,降低了边界恢复对边界约束的破坏。方法首先针对递归壳变换算法容易陷入局部最优解的问题,通过引入模拟退火算法加以解决,提出了一种改进的约束边界恢复方法;然后对于方法中的“网格实体—边界实体求交”关键环节,引入或实现了AABB树、几何精确的“线面相交”函数、交点数查询哈希算法三个关键方法加以优化;最后对于拓扑变换无法恢复的边界,实现了体内、边界两种Steiner点插入方法。实验阶段取用了Thingi10k数据集的4 000余例面网格样本进行测试,结果表明该方法有效提高了约束边界恢复的成功率,并能普遍性地减少边界恢复引入的St...     

基于舵面控制的外挂物分离仿真方法研究

基于舵面控制的外挂物分离仿真是一个典型的多学科耦合问题,需要对其气动性能、飞行性能和控制性能综合进行分析。本文主要围绕数值仿真中计算流体力学（CFD）和控制方程与刚体动力学（RBD）方程的耦合求解、高效嵌套网格技术、耦合飞行控制系统（FCS）的气动/运动/控制的多学科耦合仿真的三个关键技术进行研究。针对以上技术难点,发展了一套基于舵面控制的外挂物分离仿真方法。仿真结果表明,通过舵面控制,有效改善了外挂物在分离过程中的姿态变化剧烈的现象,提升了分离安全性与分离品质;发展的基于舵面控制的外挂物分离仿真方法能够处理复杂运动边界的非定常问题,检验飞行控制律,具有一定的工程价值。