亚声速风洞非均匀流场的修正计算

本文给出了一种修正风洞非均匀方向场和速度场的方法。对于任意的机翼-尾翼布局和风洞剖面的计算结果和实验数据作了比较,两者吻合得相当好。计算证明在非均匀流场中俯仰力矩曲线可能发生畸变。本方法为修正试验结果提供了一个手段。     

正交各向异性材料构件粘塑性损伤寿命研究

用损伤力学与粘塑性理论相结合的方法,将正交各向异性材料的粘塑性统一本构模型进行了修正和推广,建立了在蠕变与疲劳载荷交互作用下正交各向异性材料的粘塑性损伤统一本构模型。用模拟气膜冷却的空心涡轮叶片的子构件,进行了在蠕变与疲劳载荷作用下构件的损伤试验和寿命试验。将计算结果和试验结果进行了比较。     

涡轮叶片内部通道流动换热算法改进

在使用工程计算方法对涡轮叶片温度场进行计算时,往往将叶片内流通道简化成光滑或带肋的换热管元件,容易忽略 各内流管段之间的影响,造成计算得到的叶片3维温度场与真实温度场存在较大差异。针对上述问题,为了提高对涡轮叶片3维 温度场模拟的准确度,对涡轮叶片内流通道的换热流动算法进行改进。考虑涡轮内部蜿蜒通道中弯转区和弯转效应2种因素对 涡轮内部流动换热的影响,使用试验得到的数据对2种因素影响区域的换热情况进行修正,利用修正后的算法对某工作叶片进行 温度场计算,并对修正前后叶片温度场进行了对比分析。结果表明：采用修正后算法得到的蜿蜒通道内的气体温度相较于修正前 算法得到的沿程升高更多,修正后算法求得的叶片整体平均温度降低,最大温差增大。     

运用SIMPLE方法计算大曲率弯道中的湍流流动

在三维任意贴体非正交曲线坐标系下应用同位网格技术,发展了适于求解复杂流动的压力修正算法,采用速度的协变物理分量为动量方程中速度的求解变量,湍流模型选用标准的κ-ε模型.计算表明,方法具有良好的收敛性、准确性和稳定性,计算值与试验数值吻合良好,而且能预示通道内复杂的流动现象.     

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

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

复杂区域内粘性流动的数值计算

本文应用原始变量法(SIMPLEC)计算了复杂区域内的不可压和可压缩流动问题。计算时,仍使用交错网格,但在求解动量方程时采用协变分量做为原始变量,有效地消除了在任意非正交坐标系下压力和笛卡尔速度分量不耦合而产生的压力波动。在压力修正方程中考虑了密度的影响并构造出一种密度的迎风格式,保证了计算的稳定与收敛,使之预示马赫数变化范围较大的流动成为可能。对由NACA65(12A10)10构成的平面叶栅、由NGTE10C4/30C50构成的平面叶栅的不可压缩流动和平面超声速缩放喷管进行了计算。结果表明,本方法有较好的收敛性和准确性,能够适用于从亚音速到超音速的流动场合,有一定实用价值。      

斜网格修正的低速计算方法及应用

采用无坐标变换的一种有限体积法对RANS方程进行离散求解,离散方法考虑了对网格非正交情况的二阶修正,通过构造一个动量插值实现了同位网格下的压力修正算法,并发展了相应的分区算法.利用以上方法对轿车绕流和三峡坝区主流摆动问题进行了数值计算,计算结果与实验吻合较好.     

地火转移轨道误差分析与第1次中途修正计算

针对地火转移过程中出现的各种误差,基于地火转移轨道的误差传递矩阵分析误差发散的性质,在此基础上讨论如何选取轨道中途修正的时机,并基于该矩阵对地火转移轨道第一次中途轨道修正的速度增量进行估算。与微分修正方法的严格计算结果的比较表明,基于该方法定性研究地火转移轨道第1次中途修正速度增量变化和选取合适的轨道机动时机是可行的。使用蒙特卡洛数值模拟对上述方法和微分修正方法进行计算和比较,结果表明,第1次中途修正速度增量大小差异不超过1.2m/s,相对误差不超过6%。在轨道控制精度大约为1m/s的情况下使用该方法代替微分修正方法进行计算,可以节省大量的计算时间。     

状态方程对叶片外物损伤数值模拟结果的影响

分别采用Bammann粘塑性本构模型和修正Zerilli-Armstrong本构模型在LS-DYNA软件中对钢珠冲击TC4钛合金平板试件的冲击损伤过程进行了数值模拟。通过对修正Zerilli-Armstrong本构模型选择线性多项式状态方程和Gruneisen状态方程,研究了状态方程在发动机叶片FOD数值模拟中的作用。通过对比分析钢珠在不同速度下冲击TC4试件产生的弹坑外貌、尺寸以及残余应力分布,结果表明,是否使用状态方程来描述叶片材料在冲击载荷下的体积压缩或膨胀行为,对FOD数值模拟结果的影响很小。在发动机叶片的FOD数值模拟中,材料本构模型的选择和评估更为重要。     

基于本征正交分解的叶片碰摩系统降阶方法

针对叶片碰摩响应求解问题,提出了采用本征正交分解进行动力学降阶的方法。通过对快照矩阵进行本征正交分解生成投影子空间,将系统动力学方程投影到子空间进行模型降阶,并结合数值积分方法进行了碰摩响应求解。基于降阶模型分析了不同转速及侵入量参数下的叶片碰摩响应,并与全阶模型进行了对比。结果表明:降阶模型的时域响应幅值偏差小于5%,计算效率了提升98.4%;通过改变叶片转速、侵入量参数验证了降阶模型的鲁棒性,并且发现随着转速、侵入量的增加,本征正交模态能量在低阶与高阶之间发生转移,并呈现出不同的传递规律,由转速引起的模态能量转移与结构的固有频率存在一定关系。该方法及结论可为叶片碰摩分析及故障诊断提供依据。      

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

非结构混合网格高雷诺数流动计算中一种改进的LU-SGS隐式格式(英文)

The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology,low numerical complexity and modest memory require...     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

基于RANS/LES方法的超声速底部流场数值模拟

分别采用基于两方程k-ω剪切应力输运(SST)湍流模型的延迟DES(DDES)、更改的DDES(MDDES)和改进的DDES(IDDES)方法,并引入可压缩修正,结合三阶MUSCL-Roe和五阶WENO-Roe两种空间离散格式,针对超声速底部的复杂流动现象,开展了数值模拟研究。计算结果表明本文方法能够捕捉到超声速底部流动中丰富的湍流结构,通过分析计算结果对超声速底部的流动机理有了进一步的认识,为下一步的超声速底部流动减阻改进和雷诺平均NavierStokes/大涡模拟(RANS/LES)方法在非定常高可压缩性流动中的应用提供了参考。通过对比分析不同空间离散格式的计算结果研究了数值耗散对计算的影响,五阶WENO-Roe格式的计算结果与实验结果吻合良好;对不同RANS/LES混合方法的计算结果进行了对比分析,结果表明IDDES方法在近壁区的表现优于DDES和MDDES方法。     

双管爆轰单引射增推机理的数值研究

基于二维Euler方程,结合高精度Roe/HLL(Harten-Lax-Van Leer)混合格式与自适应网格加密技术,对引射模态下的双管爆轰非定常单次引射复杂流场进行了数值模拟.结果分别描述了方管爆轰波喷出及引射管内二次爆轰波的传播过程.双爆轰波喷出爆轰管后,以欠膨胀射流形式向外传播,同时,它们相互作用延长了引射管尾...     

矢通量分裂-TVD杂交新格式及其应用

在贴体曲线坐标系中,本文提出了一个高分辨率的新格式,它是Steger-Warming的矢通量分裂与Harten的TVD(Total Variation Diminishing)格式的杂交,用于计算跨音速定常流动和捕获激波。典型的跨音速内流和叶栅流的算例表明:在60×15网格下捕获的激波过渡区不超过两个网格,证明其分辨率较高;分析激波附近的数值结果,没出现“低亏、过跳”、伪振荡现象;与实验结果的比较表明,对定常流在无人为附加耗散项的条件下它能给出较高质量无数值波动的激波流场解。     

基于动量源方法的直升机旋翼/机身流场数值模拟

建立了一套基于动量源方法的直升机旋翼/机身流场数值模拟方法.将旋翼对流场的作用以动量源项的形式代表,建立了一个包含动量源项的CFD方法.针对直升机旋翼/机身流场的特点以及动量源方法对网格系统的要求,采用了一种混合网格生成方法.通过结合旋翼桨叶的运动方式、几何特征及气动特性,建立一个包含动量源项的N-S方程的旋翼流场计算方法和迭代流程,并采用该方法进行直升机旋翼、机身和旋翼/机身下洗流场的数值模拟,得到了关于旋翼/机身干扰流场的一些有意义的结论.     

二维进气道不启动流场非定常特性的混合LES/RANS模拟

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)方程模拟方法结合5阶WENO(weighted essentially non-oscillatory)格式对马赫数为3的来流中、内收缩比为1.5的不启动状态下的二维进气道进行了计算,再现了不启动进气道中的非定常流场.计算结果表明:所采用的模拟方法对入口处的平均绝热壁温、摩擦速度和雷诺应力的计算精度较好,进气道不启动流场中激波波系和分离区存在大时空尺度的低频运动,其占主导的特征频率和典型的激波/湍流边界层干扰问题中激波和分离区的低频频率接近,且进气道出现了间歇性的启动状态.      

旋翼VPM/CFD耦合模型的建立及其在小前进比气动分析中的应用

为了提高旋翼计算流体力学（CFD）流场计算效率及克服其在尾迹捕捉上的不足,将旋翼黏性涡粒子方法（VPM）与CFD分析相结合,建立了一个新的旋翼VPM/CFD耦合气动分析模型。在该模型中,采用VPM分析以实现对旋翼尾迹中黏性涡的高效捕捉而不引入数值耗散,而采用CFD分析用于精确地模拟旋翼桨叶近体区域内复杂的流动现象,同时也为VPM分析提供一个较高精度的涡源模型。至于两者信息交换,则使用集中涡源法将CFD信息传递至VPM分析,而VPM计算得到的尾迹信息则通过边界修正施加至CFD域的远场边界上,从而可以鲁棒地实现CFD域与VPM域的耦合计算。在此基础上,对“Helishape 7A旋翼”小前进比前飞桨-涡干扰（BVI）状态进行了较为深入的气动分析,计算结果表明：与全CFD计算比较,建立的VPM/CFD耦合分析模型可以有效地避免旋翼尾迹区桨尖涡的数值耗散,从而更加可靠地捕捉桨-涡干扰状态下的桨叶非定常气动载荷脉动,同时对于本文算例,计算效率可以提高30%以上。