一个解析三维非结构网格生成方法

给出了一个三维非结构网格生成方法。提出弱区域指示函数描述三维封闭流场形状,并证明了该函数的存在不唯一性。利用该函数基于Delaunay方法依据流场尺度生成网格结构。其中,边界引力有效解决了流场边界恢复和内嵌边界问题,弹簧振子和正四面体趋进技术有效调整网格结构并提高网格质量。最后,弹簧振子和边界引力随着弱区域指示函数的演变与监测函数的指导成功应用到自适应网格生成。     

双边界代数法在喷管内流场网格生成中的应用

用一种双边界代数方法生成二维网格，通过在超限插值函数中引入边界的一阶导数实现网格与边界是的正交，采用几何一致性条件保证物理区域不变形，发展一系列拉伸方程控制网格在边界或给定位置的间距。应用此方法生成了满意的喷管内流场网格，网格生成效率明显优于微分方程网格生成方法。     

新型重叠网格洞面优化方法及其应用

高度自动化的洞面优化方法是重叠网格关键技术之一。通过对割补法的分析,针对其洞面切割后洞边界位置不确定性及不可预知性问题,提出将其与物面距离优化准则相结合的新型混合洞面优化方法。该方法针对各物体间体网格,采用物面距离优化准则,将洞边界位置控制在物体中间的位置,使重叠区域更合理、可预测;针对背景网格,将其物面距离设置为较大值,采用物面距离优化准则将其落入物体网格内部的单元一并挖去;针对物面处重叠的网格,采用割补法对其进行切割和填补,保证物面及近壁面网格得以有效的重叠。该方法在保证高适用性的前提下提高了网格重叠质量,自动化程度高,无需人工干预。3个典型复杂流动算例计算结果与实验结果吻合良好,网格重叠区域流场变量传递正确,等值线过渡光滑,流场刻画准确,证明该方法准确可靠。     

圆柱管多分支联接三维分块贴体网格的耦合生成

利用分块耦合技术生成了两圆柱管成任意角度,管径不相等的三分支联接的三维贴体网格,所介绍的网格生成方法在分块生成网格时能保证穿过分界面上的网格线是光滑连续的。此外,用方法生成网格时间通过边界上网格的疏密或者稳定边界上网格线之间的夹角来控制内部区域的网格疏密。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

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

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

有限元三维网格自动生成

讨论了为提高有限元法求解流动问题的精度和效率而寻求自动生成三维网格的方法。给出了利用代数变换、输入极少量区域边界数据即可插值生成三维网格的方法。最后举出了自动生成三维网格的例子。     

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

用代数法生成飞机翼型进气道等复杂形状物体的网格

本文论述了研制网格生成软件对于促进计算流体力学研究和应用的重要意义。着重介绍了网格生成的多面方法,编制了多面法网格生成程序,生成了飞机翼型、进气道等复杂形状物体的网格.该方法可快速生成二维或三维边界正交网格,同时还可有效地控制区域内网线的分布和边界上网格点的分布.     

并行化非结构重叠网格隐式装配技术

重叠网格技术广泛应用于复杂外形及多体分离问题的数值模拟,但是超大规模重叠网格装配的鲁棒性和效率依然是实际应用中的瓶颈。基于格心型的有限体积格式,发展了一种并行化的非结构重叠网格隐式装配技术。该技术采用节点的壁面距离作为插值边界的判断准则,并通过物理边界推进确定出活跃区域。此外,基于网格分区策略实现了隐式装配技术的并行化,为了能够适应超大规模的计算网格,在每一个网格子区内单独建立用于查询的交替数字树（ADT）数据结构,整体网格点在各个子区内进行查询,并通过信息通讯接口（MPI）传递查询结果并进行逻辑判断。该并行化的非结构重叠网格隐式装配技术具有较高的自动化程度,通过多体装配算例以及机翼外挂物分离算例,证明了该方法的实用性。     

通道中有稳定器的湍流流场分区计算方法

采用任意曲线坐标系对含有Ｖ型稳定器的通道内湍流流场进行数值计算。对含有Ｖ型槽的复杂流场中网格分区方法作了详细的分析,提出一个新的网格划分方法。数值计算采用曲线坐标系下非交错网格的ＳＩＭＰＬＥ方法,计算时对整个流场进行分区迭代,直到得到收敛结果。文中对一含有Ｖ型稳定器的直通道湍流流场进行了数值计算,结果合理。     

计算流体力学中的网格生成方法及其应用

讨论了多种网格生成方法及其流场的数值模拟。给出了一种生成现代双立尾战斗机外形的分区对接网格和Ｅｕｌｅｒ方程分区解。给出了一种带副翼偏转的三维机翼分区网格及绕流的Ｅｕｌｅｒ方程解。以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法。     

转速对高速旋转冲压增程弹用进气道影响数值模拟

结合分区对接网格技术以及二阶精度区域分解算法，对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟。得到了高速旋转工况下对应于不同来流攻角和旋转角速度，临界工况时，超声速进气道内外流场复杂的波系结构。随着旋转角速度的提高，进气道总压恢复系数和动能效率均有所降低，而流场畸变指数则显著增大。特别是当转速达到20kr／min，进气道总压恢复系数和动能效率下降趋势以及流场畸变指数增大趋势更明显。攻角的存在对冲压发动机进气道的总体性能产生了负面影响。     

带有双层壁扩压器的波瓣喷管混合流场的数值计算

采用三维贴体曲线坐标网格,边界网格加密且正交,在整个计算区域进行全场计算。为避免因波瓣造成的网格强烈的非正交引起解的发散,采用了Chen和Kim修正的k-ε湍流模型及同位网格SIMPLC计算方法,对带有双层壁扩压器的波瓣喷管混合流动进行了数值计算和分析。对波瓣及双层壁,采用大粘性的方法解决流固耦合。计算结果表明:在双层壁间有外界冷气流被引射进入,形成了壁面的冷却气流;自波瓣出口截面沿流向产生的环流速度场,强化了主次流的掺混,速度分布渐趋均匀。计算结果与实验数据符合良好,二者在离开波瓣52 mm的混合管内相对主流速度的最大误差为21.65 %。     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

尾喷管内部及其射流流场的数值模拟

通过求解三维Ｎ－Ｓ方程,对某型发动机尾喷管在亚临界及超临界状态下的性能进行了数值模拟,尤其对喷管射流流场及其对喷管内部流动的影响进行了研究。差分格式采用具有ＴＶＤ性质的三阶精度的Ｇｏｄｕｎｏｖ格式,射流场及内部流场简采用分区算法。数值试验表明,所采用的算法可以很好地处理来自于任意方向上的任意间断。     

三维复杂无粘流场的自适应八叉树结构直角网格算法模拟

运用自适应八叉树结构直角网格对三维复杂流场进行了划分 ,并采用点面判别法进行外法向量的判断 ,使与之关联的与物面相交网格的体积计算更为简便准确。在此基础上 ,采用以中心差分为基础的 Jameson的有限体积法进行 Euler方程求解 ,在通量计算过程中 ,采用以面为基础的算法 ,减少了计算工作量。对多个外挂、翼身组合体及全机等复杂外形绕流流场进行了模拟 ,均取得了较好的效果     

欧拉方程多重网格方法及三维跨声速大迎角流动数值计算

本文给出一种粗、细网格时间步长可调整匹配至最优组合的多重网格方法以获得在最细网格上最大的加速收敛效益，二重网格中节省机时可达３５％。文中的数值方法可用于飞行器跨声速大迎角无粘有旋的涡干扰流场的模拟研究。     

用于复杂流场数值模拟的重叠／嵌套网格方法

介绍了一类用于复杂流场数值模拟的重叠／嵌套网格方法。它是一类多块结构化网格,可用于有或无多重网格加速收敛技术的有限体积、有限差分等计算。其基本思想是把复杂流场的计算区域化分成相对简单的子区域,每个子区域上的网格可以独立生成;网格界面之间的流场信息交换通过三线性插值来实现,通过修改流场解算方法来排除洞中的点。     

多体全自由飞过程的数值模拟

以三维非结构粘性直角网格重构为构架,数值求解了固定空间网格瞬间运动边界的N-S方程,逼近了流场的非定常过程,结合六自由度的弹道方程,对单体或多体运动过程进行了全自由飞数值模拟,数值模拟结果和级间分离的风洞自由飞试验进行了对比,初步的计算结果说明,全自由飞数值模拟在一定程度上可以描述风洞试验的飞行轨迹,并且适应能力更强,整个全自由飞数值模拟系统可以为更复杂的单体机动飞行或多体运动过程提供快速实用的数值模拟方法.