具有无穷时滞的非线性Volterra反应扩散系统解的先验界

详细讨论了具有无穷时滞的非线性Volterra反应扩散系统解的先验界与吸引性,运用上下解方法和单调迭代方法,得到了在某些假定下当t→∞时解趋于先验界及在某些特殊情况下t→∞时解趋于正常数的充分条件。     

跨声压气机内确定应力耦合问题研究

基于三维定常Denton程序发展了二维定常、二维非定常、三维非定常计算程序,以二维跨声压气机叶栅、三维跨声单级压气机近失速点为算例,进行确定应力耦合问题的研究.结果表明:在确定应力应用到时均方程的耦合处理中存在一些因素导致时均方程的解和非定常解的时均值有一定偏差,它们分别是空间差分离散引入的误差、黏性相关项引入的误差、掺混面处理引入的误差.在多排耦合问题中掺混面的处理方式至关重要,对耦合计算结果有较大的影响,简单的掺混面处理并不能获得理想的效果.      

三维跨声速多级压气机内确定应力耦合问题研究

为开展确定应力耦合效果的研究,基于三维定常Denton程序发展了三维非定常计算程序,以多级跨声速压气机近失速点为算例,进行确定应力耦合精度的研究。除定常计算和非定常计算外,还采用了三种计算方式,分别是多级环境下非定常时均方程耦合确定应力的解(COU2),提取非定常精确解时均的值作为S2进口边界条件进行单独S2计算域的定常计算(ST1);在ST1基础上耦合确定应力的解(COU1)。结果表明:三维压气机算例中出口总压和马赫数的径向曲线对比显示COU1的误差水平分别为ST1的53%和34%;验证了确定应力项的耦合可以使定常计算的结果更接近非定常计算的时均解;对比发现ST1计算精度较COU2高,误差水平分别降低至37%和32%。说明在该算例中精确的边界条件比精确的确定应力发挥更大的作用。非定常时均结果中近压力面的熵值比近吸力面的熵值要高,反映出非定常相互作用中上游尾迹向压力面迁移的效应,Adamczyk时均方程耦合确定应力后得到了类似的结果。     

旋流畸变影响低速轴流压气机数值模拟研究

为研究进气旋流畸变对轴流压气机性能和稳定性影响的机理,设计了一种叶片式旋流畸变发生器,分别采用定常和非定常计算方法对旋流畸变发生器和两级低速轴流压气机进行了联合数值模拟研究。定常计算结果表明:旋流畸变发生器的叶片弯角越大、对涡强度越高;对涡旋流导致压气机总压比和效率降低,稳定边界流量明显增加;对涡强度越大,对压气机性能和稳定性影响越严重;对涡的反向涡对压气机影响较大,同向涡影响较小;反向涡造成部分叶片的叶根区域吸力面附面层提前分离,是导致压气机失稳的主因;非定常解的总压比和效率均略高于定常解,两种计算方法得到的稳定边界差异明显,非定常解的稳定边界流量比定常解小19.2%,非定常计算结果更接近试验值。     

基于N-S方程的非定常地面效应研究

研究了二维翼型在移动地面和固定地面条件下的地面效应及流动特性问题.流场求解采用雷诺平均N-S方程,湍流模型采用k-ω和SST模型.对NACA0015对称翼型在离地速度为1 m/s情况下的非定常绕流进行了计算,选取移动地面条件下不同高度的定常解与非定常计算结果进行了比较.结果表明,考虑非定常效应时升力和阻力的变化明显小于定常解的结果,相对而言力矩的差别较小,固定地面条件下的升力系数和阻力系数略大于移动地面条件下的计算结果.     

变攻角下端壁非定常脉动抽吸对高负荷压气机叶栅性能的影响

为进一步探究非定常脉冲抽吸控制高负荷压气机叶栅流动分离的机理,考察非定常脉动抽吸在变攻角下的适应性和可行性,采用非定常数值方法,系统研究了变攻角下,非定常脉动抽吸对流场性能的影响,并将其与传统的定常抽吸进行了比较分析。结果表明,在设计攻角下,保证相同的时均抽吸量,非定常脉动抽吸控制效果明显优于定常抽吸;在时均抽吸量ms=0.4%时,在给定的激励频率范围内,非定常脉动抽吸都展现出更好的性能,在最优频率时,损失减小了9.4%,静压升提高了12.9%,相比于定常抽吸损失减小了4.2%,静压升提高了4.7%。在变攻角下,在给定的激励频率范围内,非定常脉动抽吸控制效果相比于定常抽吸仍具有较大优势;但大攻角下,非定常脉动抽吸和定常抽吸控制效果均有所下降。     

超声速进气道加速/减速过程起动/不起动现象研究

为了研究加速/减速过程中超声速进气道的起动/不起动现象,对二元超声速进气道的二维流场进行了定常/非定常数值模拟.分析了在来流马赫数为1.75~2.05范围内4种不同加速度绝对值对进气道气动性能的影响,并对定常/非定常模拟方法的结果进行了对比.结果表明:与定常模拟时相比,非定常数值模拟时进气道的流场特征在加速起动过程中和减速不起动过程中均存在明显的滞后现象,且加速度绝对值越大,滞后现象越明显;非定常模拟时的进气道处于起动状态时,其气动性能存在迟滞现象;当加速度绝对值一定时,进气道起动过程耗费的时间大于不起动过程耗费的时间.      

旋翼干扰流场数值模拟中的新型Blend低速预处理方法

针对定常低速预处理在双时间步长内收敛性差、旋翼干扰流场数值模拟中非定常低速预处理自动退化为无预处理的问题,在非定常低速预处理方法的基础上,采用量级分析方法研究了对流项和黏性项量级,通过组合定常和非定常预处理特征值矩阵,建立了一种新型Blend低速预处理方法。通过圆柱非定常扰流、前飞直升机旋翼/机身干扰非定常流场和直升机/舰船干扰非定常流场的数值模拟,对各低速预处理方法的收敛性和预测精度进行了细致研究。结果表明,本文方法解决了在旋翼干扰流场数值模拟时非定常低速预处理自动退化为无预处理的问题,相比定常低速预处理大幅提高了收敛效率。     

离心叶轮与扩压器相互干扰数值模拟

采用定常与非定常数值模拟方法,对某小型离心压气机进行数值计算,首先比较了定常结果与非定常时均结果的特性;然后结合离心叶轮出口和径向扩压器进口典型位置非定常压力脉动的频域特性,研究了离心叶轮与扩压器之间的转静干涉作用;最后着重分析了离心叶轮带与不带分流叶片对径向扩压器进口压力脉动频率特性的影响。结果表明:在相同流量下,非定常时均结果的压比偏高约3.7%,效率偏高0.1%;在带有分流叶片的离心压气机中,离心叶轮叶片出口与扩压器叶片进口的脉动除了叶片通过频率的非定常脉动外,分流叶片尾缘会产生与之强度相当的2倍叶片通过频率的非定常压力脉动。同时,分流叶片会对径向扩压器进口产生更强的0.5倍叶片通过频率非定常压力脉动,这些频率的非定常压力脉动很可能会对离心压气机叶片的疲劳强度造成影响。     

振动叶栅中尾迹/激波相干的数值研究

研究了尾迹/激波相干对振动叶栅非定常流动的影响.数值方法利用了先进的谐波方法.该方法基于傅立叶变换,把叶栅非定常流变量分解为一时间平均量和若干谐波分量,继而,时域非定常N-S方程被转化为频域定常傅立叶系数方程.定常傅立叶系数方程的求解采用了伪时间推进方法.为加速数值解收敛,采用了本地时间步长和多重网格技术.通过NASA67转子算例结果表明:上游尾迹对振动叶栅非定常流动具有显著影响.当存在尾迹相干时,叶片通道内的激波在轴向方向上表现出较大的漂移运动,叶片所受的非定常力幅值变大,同时失去了时间周期性.      

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

航空燃气轮机挤压油膜阻尼器油膜温度场的理论分析

油膜温度是挤压油膜阻尼器的一项重要工作参数。本文运用有限元法对挤压油膜阻尼器油膜温度场进行分析,即根据油膜边界条件,联立求解适用于阻尼器的雷诺方程、能量方程、滑油粘温特性方程得到油膜温度分布。文中进一步分析航空发动机阻尼器的工作环境,提出启动阶段,滑油温升由于阻尼减振发热所致;稳态阶段主要由于油膜环向滑油的强迫对流换热,并提出用试凑法来求解稳态滑油工作温度。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

Godunov方法在复杂流场数值模拟中的应用

本文根据Ｖａｎ Ｌｅｅｒ提出的方法给出了求解三维非定常Ｅｕｌｅｒ方程和三维定常Ｅｕｌｅｒ方程的二阶Ｇｏｄｕｎｏｖ方法。方法要点为：在各体积离散单元中,用流动变量的分段线性分布代替一阶方法中的分段常量分布,在确定线性分布的斜率中引入单调性限制条件;在离散单元界面构造Ｒｉｅｍａｎｎ问题,由Ｒｉｅｍａｎｎ间断解确定单元侧面的通量;用预测－修正两步法对离散方程推进求解,用上述方法结合单区网格或多区重迭网格     

一个n阶微分方程的两点边值问题

有一些n阶非线性微分方程的两点边值问题已被人们讨论过,大多是在微分方程的右端函数有界、满足李普希兹条件或在相应边值问题存在上下解的情况下讨论的。[3]、[4]在这篇文章里,讨论了一个n阶非线性微分方程的两点线性边值问题,是首先通过n-2次累次积分将原方程化成了一个二阶微积分方程边值问题,并用拓扑度理论讨论了解的存在性,同时给出了解的唯一性条件。推广了文献[1]的结果。     

用单调相关格式计算跨音速定常和非定常翼型绕流

在跨音速非定常势流计算方面,人们一直在努力寻求处理包括激波运动的非线性问题的能力和计算效率的提高。Ballhaus和Goorjian（1977年）用交替方向隐式差分法（ADI）成功地直接求解了时间域的非定常跨音速小扰动方程（此法称为“LTRAN2”）。文献[5]引用坐标变换也做了类似LTRAN2的工作。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

气冷叶片表面局部对流换热系数的数值计算

本文提出了联立求解二维可压缩边界层方程组和稳定导热方程的迭代法,从而修正了文献[1]假定物体表面温度等于来流总温的局限;同时提供了用这种方法计算气冷叶片表面局部对流换热系数的完整的FORTRAN程序。文中用近似积分法求解边界层方程组,用有限元素法计算叶片的稳定温度场。算例的结果与有关的理论及实验结果相接近。     

时间依赖边界条件在亚声速自由射流数值模拟中的应用

对亚声速自由射流进行数值模拟时，无法预先给定远场边界和无反射入流边界上的速度值。以往采用的黎曼不变量在自由边界上给定无反射条件的处理方法，只给出了这些边界法线方向上的一维特征关系。时间依赖边界条件是一种高精度的边界处理方法，可以很方便地处理计算中遇到的各种边界，广泛地应用于DNS中。我们将这种边界条件用于定常流动的计算中，计算了出口马赫数为0．4和0．7的轴对称自由射流，获得了与实验结果吻合的数值模拟结果。     

采用近似边界求解壁板颤振的CBS有限元算法

为了简化Euler方程求解壁板颤振问题,将近似边界条件成功地推广到CBS(characteristic-based split)有限元方法求解壁板颤振领域,把运动壁面的影响转化为Euler方程的边界条件,避免了对动网格的处理.通过对近似边界方法、动网格方法以及3阶活塞理论的数值模拟结果进行对比发现:首先,加载近似边界的CBS有限元方法操作非常简便,且模拟结果与动网格方法完全吻合;其次,Euler方程可以准确捕捉到壁板前后两端不光滑尖点所引起的压力突跳,而3阶活塞理论则无法做到;最后,采用近似边界的CBS有限元方法在求解无黏可压缩流动中细长体小变形的流-固耦合问题上具有潜力和优势.