低展弦比跨音速轴流风扇转子流场三维数值模拟

应用基于时间推进的有限差分法求解跨音速压气机风扇转子内部三维粘性流场．该方法以新型ＬＵ隐式格式和改良型 ＭＵＳＣＬ ＴＶＤ格式为基础,对三维可压缩雷诺平均 Ｎａｖｉｅ－Ｓｔｏｋｅｓ方程和低雷诺数ｑ－ω双方程湍流模型进行求解．计算得到了ＮＡＳＡ Ｌｅｗｉｓ ６７跨音速、低展弦比轴流风扇转子的性能曲线并重点分析了近最高效率点工况和近失速工况下的内部流场．计算与实验结果的对比表明此方法能够得到三维粘性流场的流动特性且计算精度较高,可用来模拟跨音速风扇转子内部流动．     

机翼跨声速非常绕流IA型反命题变域变分有限元解

以文[1]提出的二维振荡机翼含激波跨声速非定常绕流IA型反命题变分原理为基础，构建求解IA型反命题的有限元解法。构造了三维时空可变节点有限元来捕获自由尾涡面和翼面几何形状，跨声速流中的激波用人工密度法捕获。在远场边界上采用简化的无反射边界条件，新型非定常Kutta条件被用于处理尾缘条件。用该方法，根据翼型跨声速非定常绕流翼面压力分布求解IA型反命题，得到了NACA64A010翼型的几何形态，计算结果令人满意。     

先进双发民用飞机跨音速绕流数值模拟

采用求解椭圆型和抛物型偏微分方程的方法及多块风格技术，生成机翼／机身／短舱／挂架复杂组合体外形的贴体及与边界正交的计算网格。应用有限体积法高效求解欧拉方程，模拟先进双发民用飞机的跨音速绕流，并计入机翼边界的粘性效应，获得存在短舱、挂架干扰的翼身组合体气动特性，计算表明，短舱，挂架对超临界机翼的气动特性有显著的影响，计算与实验模型的机翼表面压力分布（测量值）符合良好。     

轴流叶轮机三维粘性非定常流动计算的双时间方法

采用双时间方法研究了单级跨音速压气机的三维非定常粘性流场计算问题,以ＬＵ－ＳＧＳ隐式解法快速获得虚拟时间域的定常解。对某跨音速单级压气机同时进行了双时间方法和显式方法的非定常流数值计算,双时间法与ＬＵ－ＳＧＳ隐式解法相结合的非定常粘流计算模式可以显著地节省计算时间。     

三维机翼跨音速颤振的分布式并行计算

以多台微机组成分布式计算系统,利用无限插值方法(TFI)生成三维多块贴体运动网格,以Navier Stokes方程为控制方程,求解三维机翼的跨音速非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义位移响应的时间历程,根据广义位移的时间历程的衰减、等幅和发散振荡等情况确定机翼跨音速颤振临界条件.通过算例验证,计算结果与实验结果和理论分析相吻合.     

机翼跨声速非定常绕流IA型反命题变域变分有限元解

以文[1]提出的二维振荡机翼含激波跨声速非定常绕流I_A型反命题变分原理为基础,构建求解I_A型反命题的有限元解法。构造了三维时空可变节点有限元来捕获自由尾涡面和翼面几何形状,跨声速流中的激波用人工密度法捕获。在远场边界上采用简化的无反射边界条件,新型非定常Kutta条件被用于处理尾缘条件。用该方法,根据翼型跨声速非定常绕流翼面压力分布求解I_A型反命题,得到了NACA64A010翼型的几何形状,计算结果令人满意。     

密实炸药床跨音速两相反应流的数值模拟

本文建立了密实炸药床跨音速两相反应流的二维非定常连续介质数学模型。将SIMPLE型的数值计算方法引入到密实炸药床燃烧转爆轰跨音速流动过程的数值计算中,其中对跨音速流动的处理和可压缩两相流压力校正方程的建立提出了改进方法。以无起爆药雷管作为算例。结果表明,本法较好地克服了跨音速两相流数值解中通常出现的数值振荡现象。     

利用N-S方程模拟机翼气动弹性的一种计算方法

利用一种双时间方法求解三维非定常N-S方程,得到与任意非定常运动对应的气动力,在求解非定常气动力的同时,在时间域内用二阶龙格 库塔方法求解机翼弹性运动方程,从而模拟粘性流动中的气动弹性全过程.为保证网格生成效率,采用无限插值理论生成O-H型代数网格,考虑了机翼变形时的网格生成问题,并得到计算结果.     

二阶STC格式及其对跨音速湍流流动的数值模拟

ＳＴＣ格式是一种求解积分形式物理守恒律的计算方法,具有守恒性好,精度高的特点。目前求解Ｅｕｌｅｒ方程的ＳＴＣ格式已经建立起来,但求解Ｎ－Ｓ方程的ＳＴＣ格式尚在发展之中。为了使ＳＴＣ格式能够有效地求解粘性流动的问题,本文构造出了求解曲线坐标系下的二维Ｎ－Ｓ方程的二阶ＳＴＣ格式。用该格式对跨音速湍流流动问题进行了计算。     

应用隐式SMAC格式计算射流放水阀内部非定常不可压粘性流场

为了分析射流放水阀内部的非定常不可压粘性流场,在一般曲线坐标系下求解了以逆变速度为参数的不可压Ｎ－Ｓ方程和压力椭圆方程。计算采用了基于隐式ＳＭＡＣ格式的有限差分格式。在这种格式中,非定常流场的速度是通过使用Ｃｒａｎｋ－Ｎｉｃｈｏｌｓｏｎ格式,并在每一时间步处以牛顿迭代的方法得到的;而压力椭圆方程是通过使用Ｔｓｃｈｅｂｙｓｃｈｅｆｆ ＳＬＯＲ格式并交替改变计算方向而求解。通过使用交错网格和迎风差分,诸如改进的ＱＵＩＣＫ格式;计算过程中的误差和数值不稳定能够得到有效地控制。计算结果和实验观察得到的流动图谱对比表明,二者十分吻合。     

一先进跨音速离心压气机叶轮内部粘性流动数值模拟

随着计算机和计算技术的发展，三元粘性数值计算被广泛地应用于叶轮机械内部流场分析，本文使用三元Ｎ－Ｓ方程对这一先进跨音速离心压气机叶轮内部流场进行数值模拟，一是对这一成功设计进行分析验证。二是使用数值方法进行内部流场细节的模拟掌握其内部粘性流动特点，为今后设计和改进高压比跨音速离心压气机提供经验。     

跨音速全位势方程的机翼气动性能计算

采用守恒型的跨音速全位势方程计算了机翼的气动性能,包括表面压力分布,力和力矩系数。计算中未计及粘性效应。坐标变换采用Joukowsky变挟。计算采用Jameson超松弛逐次线迭代格式(SLOR)和逐次加密网格技术。以ONERA M6机翼作为算例。计算结果表明与其他采用全位势方程的计算结果以及实验结果符合很好。     

NACA-65-0012机翼的非定常反问题解

利用有限差分法解非定常欧拉方程是振荡机翼非定常气动反问题的一种求解方法。本文重点讨论了非定常气动反问题求解的计算方案与反问题固壁边界条件的处理,编写了相应的计算程序。本文以绕前缘点做小幅俯仰周期振动的NACA-65-0012翼型为初始机翼,在振幅分别为0°、2．51°的条件下进行了气动反问题计算,经过多次迭代后计算翼型的压力分布基本与目标压力吻合,这表明本文的反问题求解方案、反问题边界处理及计算程序是可行的。     

数值分析带二次流跨音速尾喷管流场

本文用有限元法数值分析带二次流可调收—扩喷管跨音速流动,紊流粘性采用ｋ—ε双方程模型来描述,并利用加权余量法中的Ｇａｌｅｒｋｉｎ法建立有限元方程,求解可压缩粘性流体流动Ｎ－Ｓ方程,喷管壁面采用壁面函数。通过四种不同喷管、不同工况下尾喷管内部流场计算,并用实验数据进行对比表明：计算结果合理,程序可靠,可利用有限元法对喷管超音速流动进行数值研究。     

复杂流动的LBGK模拟

ＬＢＧＫ方法在保留了格子气模型的一些优点的同时，克服了格子气模型的不足之处．应用ＬＢＧＫ方法，可以对更复杂的流体进行模拟．本文讨论了一种九点ＬＢＧＫ模型，并用Ｃｈａｐｍａｎ－Ｅｎｓｋｏｇ展开方法和多尺度技术证明其在二阶精度上表现为标准的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程．用该模型对圆柱绕流和管排流动进行模拟的结果显示，该模型能较好地模拟复杂流动现象，并具有一定工程应用背景．     

基于扩散速度方法粘性流场的离散涡数值模拟

基于扩散速度方法粘性流场的离散涡数值模拟吴文权，李春，卜千里（上海机械学院上海２０００９３）关键词：扩散速度模型，粘性流场，离散涡方法。ＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＯＦＶＩＳＣＯＵＳＦＬＯＷＦＩＥＬＤＳＵＳＩＮＧＤＩＦＦＵＳＩＯＮＶＥＬＯＣ...     

推求等效组态的原子谱项波函数

本文提出了推求复杂组态Ｌ－Ｓ耦合谱项波函数的系统方法，着重解决了推求多重谱项的谱项波函数这一难题，运用ＰＡＳＣＡＬ语言编制了计算程序，利用本程序可快速、准确地求出等效组态的ｓ＾Ｎ，ｐ＾Ｎ，ｄ＾Ｎ，ｇ＾Ｎ，ｉ＾Ｎ，ｋ＾Ｎ的Ｌ－Ｓ耦合谱项和谱项波函数     

数值模拟圆柱绕流旋涡运动及尾流不稳定性分析

1引言流体绕过圆柱所产生的非定常旋涡运动以及由此引起的流动不稳定性在理论和实践上都具有重要的意义。数值模拟圆柱绕流旋涡产生及演化过程，探讨圆柱尾流涡街产生的机制，控制尾迹不同速度型以抑制涡街的产生，避免涡激振动在工程上造成破坏作用具有重大实际意义。为使问题简化，本文以二维圆柱绕流作为研究对象。该流动涉及到非定常分离，旋涡的形成、运动及发展，流动不稳定性质改变等许多未完全解决的问题。B。had等山对圆柱突然起动问题作了一系列实验研究。文献门对二维圆柱绕流问题作了系统数值研究。本文采用文献门提出的差分格…     

振动叶栅中三维振荡粘性流场的压力密度消去法

1前言三维非定常流场的求解是目前国内外的一个热点研究课题山。文献［2］完成了三维可压非定常欧拉流场的求解,这一方法在求解三维非定常欧拉流场时,应用了四阶fringe－Kutta方法对控制方程进行积分,用中心差分进行空间离散,采用了四阶人工粘性项来保证计算格式稳定,计算稳定性要求严格,时间步长不能大,计算时间长。本文从非定常三维粘性N－S方程组出发,通过合理的数学方法,消去压力及密度项,得到只包含振荡速度矢量项对空间的偏微分方程组,在已知定常速度场后,这一方程组很容易求解。2基本方程在以角速度为n作旋转的相对坐…     

有限元方法在二维六角形组件堆芯临界问题计算中的应用

应用有限元方法对二维六角组件维芯的少群中子扩散临界计算进行了研究，并编制了二维少群中子扩散方程有限元计算程序２ＤＦＭＥ。以快堆３００ＭＷ－ＬＭＦＢＲ临界问题为例，作了数值计算。且与国内外有限元程序和有限差分程序ＣＩＴＡＴＩＯＮ的计算结果进行了分析比较。