可压缩流动模拟的多重网格-扰动域推进方法

针对可压缩流动数值模拟的加速需求，在作者前期提出的扰动域推进计算框架下，融入多重网格加速技术，提出了多重网格-扰动域推进方法。建立了粗、细网格上动态计算域的更新逻辑与高效算法，并针对多重网格方法的特点，对粗网格初始化、增大对流动态域等步骤进行了改进。设计了新型粗网格生成策略，可消除结构网格粗网格生成对单元数目的限制。典型算例验证表明，多重网格-扰动域推进方法可同时降低粗、细网格上的计算量并减少总迭代步数；通过在粗、细网格上同时采用随数值扰动传播而增大、随求解收敛而缩小的动态计算域，多重网格-扰动域推进方法可在传统多重网格方法加速的基础上再将计算效率提升57.9%。     

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

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

局部网格加密方法计算不可压N—S方程

本文发展了二维非定常不可压N-S方程的局部网格加密方法(LMR):需加密区域预先指定;仅对空间方向加密,分粗细二层网格,细网格覆盖在粗网格上;粗和细网格上分别用显式和隐式差分格式。结合压力修正法类的余量型差分格式,恰当地处理了粗细网格之间的信息传递,使得粗细网格交接面上质量守恒,即满足连续性条件。粗网格通过插值给细网格边界值影响细网格;细网格通过粗细网格压力Poisson方程耦合求解影响粗网格,并且由于压力Poisson方程从动量方程形成,在细网格覆盖下的粗网格上,压力Poisson方程是细网格方程的组合,从而耦合求解时粗细网格压力Poisson方程不需迭代。 本文中计算了二维方腔内的自然对流模型问题,Rayleigh数是10~6。粗网格用显式ULWC格式,细网格用隐式余量型近似因式分解格式。Poisson方程用快速直接算法PO1STG(在FISHPACK中)。粗网格距△=1/16,细网格距△=1/64,加密1/4区域,计算到定常局部网格加密方法所需CPU时间比同等均匀网格(△=1/64)节约一倍以上。     

非结构多重网格流场数值模拟研究

本文在非结构网格上应用多重网格技术来加速Euler方程的收敛过程。此方法在一系列粗网格和细网格之间进行值的传递，并且粗网格和细网格之间结构上无关。在每一重网格上用有限体积法进行计算。采用本文的方法数值求解了NACA0012翼型的绕流。     

二维Euler方程的自适应非结构直角网格算法

本文用四叉树型非结构直角网格和有限体积法求解了二维Ｅｕｌｅｒ方程,在基本网格生成过程中,采用了以几何外形为基础的网格加细方法,包括物面边界网格加细和曲率加细;并提出了一个用自适应的方法在激波处及其附近局部加细基本网格过程,提高了激波分辨率,在计算中将结构网格经常采用的以中心差分为基础的Ｊａｍｅｓｏｎ的有限体积法做了推广,减少了计算工作量,节省了机时,为今后推广到三维创造条件,本文对单个翼型,多段翼     

网格自适应策略在高升力构型计算中的应用

提出了一种基于多尺度网格误差分析的新型混合网格自适应局部加密策略。在初始网格基础上,通过多重网格计算,比较次细网格和细网格的计算结果,揭示出网格敏感区域,然后对这些网格敏感区域加密,计算自适应加密后的网格得到最终计算结果。相比于常规网格自适应策略,提出的新策略特别适合多种流动现象并存的高升力构型绕流流场计算中的敏感区域判断。NASA高升力全展构型的数值试验充分表明这种新策略的有效性。     

一种求解S形进气道三维粘性流场的多重网格算法

提出了一种基于多重网格技术,可快速求解Ｓ形进气道粘性内流场的实用工程算法：从Ｂｒａｎｄｔ的多重网格理论出发,采用隐式Ｂｅａｍ－Ｗａｒｍｉｎｇ时间推进格式,结合“锯齿式”迭代技术交替在粗、细两种网格上对雷诺平均三维非定常Ｎ－Ｓ方程进行求解。用此算法数值模拟了某型号Ｓ形进气道的内流场,与试验结果对比显示：计算得到的流场符合一般Ｓ形管道流动特性,所捕获的二次流图谱清晰合理,进气道出口截面处总压恢复系数与试验值较为接近。与普通单层网格计算方法相比,计算结果精度高,收敛速度快。     

HyperFLOW 软件数值模拟TrapWing 高升力外形北大核心CSCD

采用High-Lift研讨会提供的梯形翼外形(TrapWing),利用自主研发的基于结构、非结构网格的通用CFD软件HyperFLOW进行了数值计算,以评估其对复杂外形低速流的模拟能力。分别采用了三套不同拓扑结构的计算网格,包括两套非结构/混合网格和一套多块结构网格,每套网格又分为粗、中、细三种密度不同的网格数量以考察其网格收敛性。利用Richardson插值法,对计算结果开展了可信度分析。结果表明,不管是结构网格还是非结构网格,HyperFLOW均建立了可接受的可信度;对于高升力外形数值模拟,SA湍流模型要比SST湍流模型模拟的更准确;在失速迎角附近,现有的二阶精度解算器仍需持续改进。     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。     

压气机叶片排中三维流场计算的多网格法

以三维有限差分形式的MacCormack格式为基础, 对相对圆柱坐标系下弱守恒型控制方程, 导出了残差在粗、细网格上的传播与回插公式, 形成了一种多网格算法。计算表明, 这种算法大幅度地提高了计算收敛速度, 而且方法本身简便易行。      

非结构化网格生成中背景网格对源的响应

背景网格对源的响应方法对非结构网格的生成具有重要意义,如果背景网格参数能恰当地反映出源的变化,则所生成的非结构网格就能满足网格设计的要求。将背景网格的三角形单元分为7个子域,根据源在单元中的位置对网格进行局部加密或变形,使每一个源点都与一个背景网格结点相对应,再通过求解热传导方程,使背景网格参数分布光滑,基本反映网格设计的要求,从而提高了非结构化网格的生成质量。     

自适应笛卡尔网格超声速黏性流动数值模拟

复杂外形/流场的高质量网格的生成往往需要占用大量人力资源,而自适应笛卡尔网格方法能够自动化生成高质量网格,具有很好的工程实用价值和应用前景。基于笛卡尔网格方法,采用叉树数据结构进行数据的存储和访问,分别从几何特征和流场解特征出发进行网格的自适应加密和粗化,发展了一种二维情况下自动、高效的自适应笛卡尔网格生成方法。从浸入边界方法出发,结合虚拟镜像对称方法和曲率修正技术进行黏性物面边界条件的处理,同时建立了多值点问题的处理技术,发展了一种在笛卡尔网格下可有效模拟黏性物面边界条件的方法。针对自适应笛卡尔网格非均匀的特点,发展了悬挂网格的处理方法,并构建了适用于自适应笛卡尔网格的黏性数值求解器。通过典型算例的考核,验证了所发展的自适应笛卡尔网格生成技术和构建的数值求解器具有较高的精度和可靠性。     

高效确定重叠网格对应关系的距离减缩法及其应用

提出了“距离减缩法”,可以通过搜索有限几个数目很少的点来寻找到目标单元,从而可以快速确定重叠网格之间的对应关系,应用于重叠网格洞点识别与插值单元寻址两个方面,并且还可以结合Inverse Map方法高效建立网格之间的映射关系。该方法将三维搜索化为准一维搜索,大幅度提高搜索效率,且算法通用性强、易于实现。应用表明：利用该方法可以迅速识别洞点,与等参变换结合可以有方向地快速锁定插值单元,找点效率高,适用于计算流体力学(CFD)三维非定常计算中有相对运动的重叠网格系统在每一个时间站位都需要重新进行洞点识别与插值单元寻址的计算,可以有效节省计算时间,具有工程实用价值。     

一种适合于旋翼前飞非定常流场计算的新型运动嵌套网格方法

针对旋翼流场计算中关键的运动嵌套网格技术,提出了一个新的网格单元关系快速判断方法。该方法通过计算特殊网格点在背景网格上的位置来建立桨叶网格在背景网格中所占据范围的关系,并用于嵌套网格洞单元的识别,又进一步将该方法推广用于Inverse Map的建立,在此基础上,建立了一套用于前飞流场计算的运动嵌套方法。通过多种前飞旋翼测试算例,分别评估了背景网格、桨叶网格尺寸改变对计算时间的影响,同时验证了该方法的可靠性,结果表明构建网格嵌套关系所用时间随桨叶网格和背景网格尺寸的变化幅度较小,且具有良好的鲁棒性,能够满足旋翼前飞状态计算的需要。最后,在基于前飞非定常N-S/Euler流场求解器中,采用该方法模拟了UH-60A和7A旋翼悬停/前飞状态的气动特性,通过计算时间的对比体现出该方法的高效性,并且流场计算结果与试验数据对比表明该方法可以有效地用于旋翼前飞非定常流场的计算。     

分层嵌套重叠网格自适应树结构动态组装方法

采用重叠网格可以有效地进行复杂流动的大规模数值模拟,特别是包含运动部件(如旋翼、投弹)的动态模拟。本文将树结构的自适应直角网格用于重叠网格组装过程中的切割和贡献单元的搜索,大大加快重叠网格的组装速度。通过二叉树自适应直角网格对物体外形进行离散,实现切割过程的快速定位;采用八叉树自适应直角网格对流场区域进行离散,高效地搜索贡献单元。使用基于壁面距离准则的重叠区域最小化方法和分层嵌套重叠策略,能提高重叠网格组装的效率和质量。对于具有运动部件的动态重叠网格问题,采用多个二/八叉树减少组装过程中信息更新的冗余计算,从而大幅度减少重叠网格组装的时间消耗。实际算例的重叠网格组装结果说明本文发展的重叠网格组装方法具有很高的计算效率,可以满足运动边界复杂流动问题的动态计算要求。     

X—38飞行器无粘绕流的非结构网格模拟

本文采用非结构网格，利用DLR所发展的计算软件TAU对X－38飞行器高超声速无粘绕流进行数值计算。为了生成有效的计算网格及提高对流场细节的分辨率，在求解过程中应用了网格自适应,并研究了网格密度对气动力的影响程度；同时与现有的实验结果和用结构网格得到的结果进行了比较，符合程度令人满意。     

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

PVM平台下通量守恒的嵌入自适应网格解激波和底部流

运用嵌入式自适应网格技术来提高结构网格求解复杂流场的分辨率。该方法根据流场的形成,将密集网格嵌入到物理量剧烈变化的区域,使求解精度得以提高的同时,避免较大的计算开销,从而达到求解效率和求解精度之间的平衡。 Ｐ Ｖ Ｍ 平台使得被嵌入的网格可以在网络中不同的机器上完成,大幅度降低了墙上时间。原始网格与被嵌入网格之间数值通量守恒的实现增强了该方法在捕捉间断时的有效性和可靠性。     

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

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

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

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