有无化学反应的射流及尾流的稳定性(英文)

射流或尾流可视为无固壁边界而流场速度剖面具有拐点的流动。在无粘不可压流中,有拐点的速度剖面总是不稳定的,因此,射流及尾流总比固壁边界的边界层更不稳定。对于无粘可压缩流,情况有所不同。文献[1,2]曾指出,高M数流动条件下,射流或涡层可能是稳定的。最近,在对超音速燃烧问题的研究中,考察一下化学反应对可压缩射流或尾流稳定性的影响,也许是有意义的。本文评述了有无化学反应的射流或尾流的稳定性问题,发现化学反应对射流或尾流的稳定性有显著的影响。     

流体压缩性对发动机滑油通风系统流动特性的影响

为了研究流体压缩性对滑油通风系统流动特性的影响,通过在元件流动控制方程中引入气流马赫数,建立了考虑压缩性的元件计算模型。模型利用可压损失系数来建立元件进出口流动参数之间的关系,以确定进出口气流的马赫数。以发动机滑油通风系统为对象,进行验证计算,计算结果与试验结果非常接近,最大误差不超过4%。从元件、系统两个方面,通过计算深入分析流体压缩性的影响。结果表明:相比于不可压流动,马赫数较小时,空气流量差异很小;马赫数较大时,空气流量明显小于按照不可压流动计算的流量;对于单支路系统,出口Ma为1时可压流计算的空气流量约为不可压流的78%。应用该模型,进行通风系统的一维仿真分析,有利于较为全面地了解系统内气体流动状态和提高计算精度。     

超声速气流中液滴破碎的大涡模拟

为了研究超声速气流中的液滴破碎过程,采用基于CLSVOF界面追踪方法的大涡模拟方法。该方法中用可压流和不可压流求解器分别计算超声速气流流场和液体流场,并且在分别计算气相或液相时,另一相作为该相计算时的边界条件参与计算,从而实现可压流求解器和不可压求解的融合。通过去散度化液体速度外推解决界面处大液气密度比所带来的动量误差大的问题。利用该方法对来流Ma=1.358的超声速气流中的水滴破碎进行了数值模拟。模拟结果展示了超声速气流中液滴破碎的几何结构变化细节和基本破碎形态,且与实验吻合较好。数值模拟结果揭示了Rayleigh-Taylor不稳定和气动剪切力对液滴破碎的作用过程。破碎过程中,液滴在垂直来流方向上的最大无量纲宽度(Dmax/D0)约为4.5,无量纲破碎时间(tbU∞/D0)约为29.85。     

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

自适应直角坐标网格方法在对流换热数值模拟中的应用

复杂几何边界是对流换热计算面临的主要难点之一。自适应直角坐标网格方法可以比较方便的处理复杂几何边界,但其应用目前主要集中于可压无粘流的Euler方程、可压粘性流的Navier-Stokes方程的外部绕流计算。本文应用自适应直角坐标网格方法计算了两个不可压内流对流换热算例,并把计算结果同商用CFD软件STAR-CD的计算结果进行了对比。自适应直角坐标网格方法在本文的两个算例中得到了同STAR-CD一致的结果。这说明采用该方法数值模拟对流换热问题是适用的。      

翼柱型装药燃烧室三维非定常湍流的数值模拟

采用к－ε双方程湍流模型,以ＳＩＭＰＬＥ计算程式结合ＴＤＭＡ法求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流Ｎ－Ｓ方程。在建立翼柱型装药简化模型的基础上,采用边界标志来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法,求出了多个时间步的非稳态流场结构。计算结果表明：燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡。     

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

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

关于“跨音速绕流的边界积分法”

文献[1]是一篇参加国际会议的文章,由第二作者在会上宣读。当时,我们只准备在会议上透露一信息:我们已可建立可压缩非线性位流方程的边界元解法。对方法的中心思想无意细说,拟在此项研究工作结束后再正式发表专文系统叙述。所以文献[1]写得比较简单显得晦涩。现仅就文献[2]和文献[3]对文献[1]所提的问题作如下补充说明,我们今后发表的正式文章中对此将作更全面的论述。     

一种求解涡轮平面叶栅的跨音松驰法

一、引　　言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法…     

三维水下爆炸气泡的数值模拟

在理想不可压势流假定下,采用高阶有限元离散边界积分方程求解水流场。假定泡内气体遵循等熵关系。用Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方法模拟气泡的演化。计算结果与球泡“精确解”及实验值比较,表明数学模型合理且三维数值方法精度高。计算给出气泡迁移、射流等重要动力学行为和物体受力分析     

一种结合壁面模型模拟湍流流动的间断界面浸入边界法

提出了一种具有鲁棒性的间断界面浸入边界法，用来模拟不可压湍流的物体绕流问题。开发了一种新的OpenFOAM求解器，该求解器结合了压力隐式分裂算法和k-ωSST湍流模型求解Navier-Stokes方程。为了减少对近壁面网格的要求，在实施边界条件时，采用壁面模型计算得到的壁面剪切力将速度和湍流变量修正联系起来。壁面切应力采用远离壁面的镜像点求得，镜像点上流场变量通过反距离线性插值由周围流场信息得到。使用该求解器模拟了平板绕流、NACA0012表面压强和切应力分布、Buice扩散管、圆柱绕流、方柱绕流5个湍流流动问题，本文数值模拟结果与文献结果有较好的一致性，验证了该方法模拟湍流流动的有效性和准确性。     

关于“跨音速流的边界积分法”一文的某些问题

 最近杨岞生教授和J.C.Wu教授提出了一种跨音速流的边界积分法。该方法以脉冲压强为基础建立非线性可压缩位流的边界积分表示式,利用边界元法给出了二维、三维机翼跨音速位流的计算方法。我们认为文献[1]在推导该方法的基本方程（即主管脉冲压强的方程）和求此方程的基本解这两个方面存在一些问题。     

运动物体与涡量场相互作用的不可压理论 涡量场对运动物体的作用力

本文与文献[1]一起构成运动物体与涡量场相互作用的不可压理论。本文证明:不可压流体对物体的作用方式同物面产生涡量的方式直接对应。特别是,匀速运动的三维物体所受的合力等于物面涡量流矢矩的面积分,纯运动学效应造成的表面力变化完全对消。这不仅加深了关于流体对物体作用机理的理解,而且为分析和改善物体受力状况提供了新的原理。     

湿压缩对压气机转子失稳边界工况性能的影响

针对某型跨声速压气机转子的失速边界工况,采用数值模拟方法系统研究了不同加湿条件对其性能的影响。获得了水滴粒径、速度、温度和喷水量与转子进气流量、压比、效率、比压缩功等性能参数之间的关系曲线。结果表明:在失速边界上对跨声速压气机进气加湿,可以增加空气流量,提高总压比和效率,降低比压缩功。通过对流场的分析,发现进气加湿后,向低能区注入了高动量的水滴颗粒,削弱叶顶阻塞团,将主流与泄漏流交界面和激波位置向下游推进,减小激波强度,降低叶顶载荷,增大主流区载荷,改善了叶顶流场。提出了以水滴的面积流量定量评价不同进气加湿条件对叶顶流场的作用效果。      

中心锥式进气道跨音速内外流场的混合差分计算

一、引言 文献[1]曾用跨音速定常势流的混合差分法成功地计算了皮托管式轴对称进气道内外跨音速流场,本文进一步用上述方法计算了中心锥式进气道的跨音速流场。中心锥式进气道的通道形状大多为收缩-扩张型,内通道壁面扰动较强。同时,由于中心体的存在,使传递进气道出口反压扰动的横向松弛线被隔断,因此,在皮托管式进气道计算中行之     

适于旋翼CFD模拟的高效预定边界运动嵌套网格方法

通过引入新的预定边界嵌套策略,并结合提出的"逆向边界"和"Local Direct-Map"(LDP)技术,建立了一套预定边界嵌套网格方法。"逆向边界"通过预估计算对边界进行动态调整克服了透视图方法很难明确避开流动非线性区的缺点;LDP方法解决了传统的基于Inverse-map透视图嵌套方法中分辨率与计算效率矛盾的问题。对不同操纵特性下相同嵌套网格的洞边界分布特性进行检验,结果表明预定边界嵌套方法保持了高鲁棒性并且嵌套效率提高了16.7倍。为进一步验证对非线性流场模拟的有效性以及对旋翼复杂运动的适应性,建立了一套适用于旋翼非定常流场计算的可压雷诺平均N-S方程数值模拟方法,分别对悬停和前飞状态下的C-T旋翼和UH-60A旋翼进行了数值分析。模拟结果与试验值吻合较好,表明了预定边界嵌套网格方法能够有效地用于旋翼非定常流场和气动特性的数值模拟分析。     

超声速主流中逆向喷流流场的数值模拟

本文在差分格式NND-2基础上^[1]，提出了一个修正格式，形式简单，应用方便。从薄层近似的NS方程出发，用该格式计算了超声速主流中存在逆向喷流的钝体绕流的粘性流场，成功地捕捉到流场内的各种波系和涡系结构。计算的结果与文献[2]中的实验结果比较表明二者符合很好。     

基于深度神经网络的含运动边界非定常流场预测方法研究

针对飞行器设计过程中对流固耦合系统快速预测的需求,探索基于数据驱动的非定常流场建模策略,缩短流场演化求解耗时,从而加快流固耦合系统模拟速度。流固耦合系统中流场演化求解部分等价于含运动边界的非定常流场演化。本文提出了一种基于神经网络的流场预测模型,来学习并预测运动边界附近的非定常流场演化过程。此神经网络可以基于当前流场状态及边界运动信息预测出下一个时刻的流场状态。通过不同振动频率及振幅下的运动圆柱绕流问题,测试了本文提出的神经网络模型的预测精度及泛化能力。神经网络的预测结果和计算流体力学仿真结果中流场结构吻合度较高,且通过对预测的流场数据中边界上的压力积分得到的气动力也具有较高的精度。测试结果展示了该神经网络模型具有良好的预测性能,因此该方法可以用于快速、准确获得运动边界周围非定常流场状态。     

双旋流模型燃烧室的燃烧特性

同轴双旋流进气能建立起有利于燃烧的流场,文献[1]的结果表明,同轴的内外两股旋流反向旋转时,边界层内的紊流强度和雷诺应力很高,对强化混合和燃烧有利。据文献[2]报导,用这种双旋流改装了JT-8D火焰筒,取消了火焰筒上的主燃孔与掺混孔,试验实测的慢车状态燃烧效率高达99%。表明这种燃烧室方案在燃气轮机和工业燃烧装置中有着广阔的应用前景。本文重点研究两个径向式叶片旋流器在同轴反向旋转条件下的燃烧特性,包括燃烧效率、熄火边界和出口温度场等。      

二维可压粘流N—S方程的边界元解法

本文给出了一种求解二维可压Navier-Stokes(简称N-S)方程的边界元法。使用线化技术把控制方程变成变系数的线性偏微分方程,其基本解由组合代数的方法构造,从而给出了N-S方程的积分表达式。本文计算了几例低R_e数的跨声速翼型绕流,结果表明边界元法可以应用到可压粘流中。