二元喷管无粘流的有限体积法分析

本介绍了一种计算收敛－扩散喷管内流特性的有限体积法，它将求解域分成许多有限体积单元，将欧拉方程在小体积单元上积分，从而得到半离散方程，然后采用Mac Cormack格式进行时间推进，求得定常解。计算结果与试验值吻合良好，是一种有效的计算具有喉道小曲率半径不连续的收敛－扩散喷管特性的方法。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

高阶辐射换热有限元算法研究

随着航空航天技术的发展,飞行器热结构所需承受的温度越来越高,辐射换热变得非常重要;当前大部分商业软件对于辐射问题的有限元计算方法还是基于一致表面温度和辐射热流假设,使得计算精度和网格密度的矛盾越来越严重.进行了选用高阶单元、采用高斯积分精确计算单元表面变辐射热流方法的研究,从而摆脱了一致表面温度和辐射热流的假设,使得在相同网格密度的情况下计算精度大大提高;同时,从包含辐射换热问题的有限元计算方程出发,采用与有限元数值计算时相同的积分方案,只在独立的积分点处计算辐射热流,克服了积分方法计算效率低的缺点.经与ANSYS的计算结果对比,应用辐射热流积分方法于高阶单元能大大提高计算精度;并且在相同计算精度条件下,此方法的计算效率更高,具有一定的实用价值.     

一种无条件稳定的时域推进法与有限元混合分析瞬态热传导问题

本文将一种新的无条件稳定的时域推进法与有限元结合,用于分析瞬态热传导问题。根据微-积分型热传导控制方程,对时间变量在小区闻内插值积分,变初-边值问题为一系列离散时刻的边界值问题,再应用有限元法求解之。由于推进中的每一时刻的解,都严格满足原给定的初始条件,就消除了累积误差的影响。此外,这种时域推进法对时间变量的数值积分也优于数值微分。因而,与现行的基于微分型控制方程,采用时-空有限元模型同时离散,再逐步求解的直接积分算法相比,计算精度可大大提高。不难预料,当求解较长时间后的瞬态值时,本算法的优越性会更加明显。     

极坐标算法对飞机CCIP弹道解算的优化

采用极坐标算法对某型飞机CCIP攻击模式进行弹道解算优化,该算法进行坐标变换和积分变量转换后,使弹道的积分变为一个具有确定积分上限的定积分,有效地简化了计算过程,节省了计算时间,提高了计算精度,能较好满足火控系统的精度和实时性要求。     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

Hoff型夹层板振动分析的边界积分方程法

 本文利用静力问题的基本解建立了Hoff型夹层板振动分析的边界积分方程。在数值实现过程中,边界上和区域内的未知函数均釆用线性插值以提高计算精度。数值结果表明本文方法具有未知量少和计算精度高等优点。对于低阶频率,本文结果与精确解的误差不大于1％。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

翼型跨音速小扰动势流-边界层计算

翼型绕流的外部无粘流,采用跨音速小扰动速势方程和改进小扰动速势方程对三种差分格式进行对比计算:(1)Murman-Cole非守恒式;(2)守恒式;(3)Engquist-Osher式。在相同的初值条件下,它们的计算稳定性、收敛性和收敛解基本相同,但激波位置(1)在前,(2)在后,(3)在中间。对人工时间阻尼项εφ_(xt)的作用进行了线化分析,所得结论与数值实验结果相符。二级精度格式计算表明:采用这种方案能在几乎不增加计算机时的前提下,大大提高差分计算精度。本文最后采用边界层积分方程法,计算翼型边界层,并与势流进行迭代计算。在弱激波条件下,迭代收敛解与实验值接近。     

PSO节点寻优的样条逼近微分

样条节点分布对逼近精度的影响很大,样条节点寻优模型的雅可比矩阵求解复杂。利用经典参数的PSO（粒子群算法）解算样条最优节点分布模型,以逼近残差的平方和为目标函数,每步对节点序列排序后再计算目标函数,可以提高样条逼近精度,为样条节点寻优提供了一种较好的实现方法。仿真计算表明,在一定的节点个数情况下,PSO节点寻优的逼近效果比Carl de Boor的NEWNOT程序中的方法要好。     

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

高精度的机翼静气动弹性数值模拟方法研究

气动结构耦合是研究静气动弹性问题的关键所在。采用三维Euler方程为控制方程,以有限元方法计算结构弹性变形,用细平板样条法(TPS)实现气动结构数据间的转换,进行机翼静气动弹性数值的模拟。算例结果表明,本方法是一种具有很高计算精度和工程实用性的气动弹性模拟方法。     

中心进气转静系转盘风阻扭矩数值模拟

采用计算流体力学(CFD)商业软件CFX从流动参数及几何结构两个方面对中心进气转静系盘腔内转盘风阻扭矩特性进行了数值模拟研究,同时搭建一个转盘风阻扭矩实验台测量自由盘风阻扭矩.中心进气转静系转盘风阻扭矩与经典经验公式进行了对比,结果表明:利用数值模拟方法,得到的无量纲扭矩系数与经典经验关系式计算得到的数据总体相差小于10%,利用湍流参数来表征转盘转速以及叠加流流量,对于分析转盘风阻扭矩非常便利.利用湍流参数定量分析时需注意叠加流以及转速的量级同样影响着风阻扭矩的数值.几何结构对于风阻扭矩也有影响,但相比流动参数影响较小.      

A novel method for solving shortest tool length based on compressing 3D check surfaces relative to tool postures

Solving the shortest tool length quickly under a known tool trajectory in multi-axis machining of complex channel parts is an urgent problem in industrial production. To solve this problem, a novel and efficient method is proposed which is featured by extracting only a few necessary curves from the check surface instead of sampling the entire surface. By rotating and compressing the 3D check surface relative to all tool postures, the boundaries of the area occupied by the 2D compressed surfaces are the essential elements for determining the shortest tool length. A tracking-based numerical algorithm is introduced to efficiently solve the silhouette curves which are formed in compressing. To define the multi-taper shaped tool holding system (THS) which is commonly used in production, a characterization model for THS profile is established. A model for solving the shortest tool length is finally constructed based on the critical interference relationship between the THS profile and all compressed boundary curves. For acceleration, the boundary splines are segmented according to their knot vectors. Then a new concept called the axis-aligned tool length box (AATB) is introduced, which can provide a conservative range of tool length for a spline segment. By scanning the AATBs of all spline segments, the very few effective spline segments that may ultimately determine the shortest tool length are filtered out. This acceleration method makes the solution for the shortest tool length more focused and efficient. The results of experimental examples are also reported to validate the efficiency and accuracy of the proposed algorithm.     

轴向通流旋转盘腔内换热的数值模拟

以数值模拟的方法,采用旋转坐标系稳态方程,研究了轴向通流旋转盘腔内的换热。主要讨论了流动对换热的作用以及旋转系下各力对换热的影响,给出了盘腔内的换热随各无量纲参数的变化规律。研究结果表明:哥氏力的增大削弱了轴向通流旋转盘腔内的换热,惯性力和浮升力的增大增强了换热;反映在无量纲参数上,随着进口雷诺数的提高,盘腔内的换热增强;随着瑞利数的提高(提高转速),盘腔内的换热经历一个缓慢变化-突增-缓慢变化的过程,换热的突然增强是冷气流穿透盘腔所致。      

求解叶轮机粘性流场的再生解析数值方法

应用再生核方法求解叶轮机粘性流场是一种新颖的解析数值方法, 它可以解决计算中的存贮、计算速度和精度问题。本文给出了用再生核方法求解叶轮机粘性流场的思路和方法, 以及如何在曲线坐标系中应用再生核方法, 并给出了流函数积分方程的求解过程。最后对平面扩压管道进行了数值试验, 结果与有关文献是吻合的。      

压气机和涡轮转子三维温度场计算研究

为了对压气机盘和涡轮盘低循环疲劳寿命分析提供温度场数据,采用ANSYS计算软件的热分析模块对某型航空发动机高压压气机转子瞬态温度场和低压涡轮转子三维稳态温度场进行计算研究。重点分析旋转盘腔、旋转轴系、封严篦齿、榫头装配间隙等部位的换热规律;计算结果与相关文献进行对比,验证转子温度场计算方法的可行性。计算结果表明：从慢车到最大状态过程中,高压压气机盘最大径向温差先增加后减小直至稳定,中心孔附近较厚区域温度梯度最大;低压涡轮盘中下部沿轴向盘面温差很小,整个轮盘高温区域集中在轮盘盘缘。     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。