梯形空腔热环境数值研究

对梯形空腔内的层流自然对流换热进行了数值模拟，用SIMPLE算法和中心差分格式、二阶迎风格式对该问题（Ra=1×10^3～1×10^7）进行了详细的数值计算，最终得出了梯形空腔在不同Ra数下的等温线与流线，总结分析了其内部换热规律，认为Ra≤1×10^4时空腔内部换热方式以导热为主，Ra〉1×10^4时空腔内部换热方式以对流换热为主，Rn〉1×10^7时空腔内部为湍流，应该用湍流模型进行数值计算研究。     

采用TRIP3.0模拟CHN-T1模型气动特性

采用TRIP3.0软件平台(TRIsonic Platform version 3.0),对第一届航空CFD(Computational Fluid Dynamics)可信度研讨会组委会提供的运输机标模构型(CHiNa-Transport,CHN-T1)进行了流动数值模拟。本文采用粗、中、细三套网格及百亿极细网格进行了网格收敛性研究;采用相应网格分别从气动特性、压力系数分布曲线、表面流态三个方面分析,研究了支撑装置、机翼静弹性变形以及雷诺数效应对CHN-T1构型气动特性的影响。数值模拟结果与试验结果有良好的一致性。数值模拟得到了网格收敛结果;支撑装置对力矩特性影响较大;机翼静弹性变形对气动特性影响较小;雷诺数效应对气动力特性影响较大。     

区域分裂法在湍流流场数值解中的研究

应用区域分裂算法(简称DDM)对湍流流场进行数值研究。对湍流流动控制方程建立一种通用的DDM-1计算格式S_{I,ω,给出重叠区域内变量处理方法。DDM-1的SI,ω计算格式的适用性由突扩圆管分离流流场数值计算所验证,数值计算结果与实验数据比较表明:采用SI,ω计算格式所预测的再附点位置及平均轴向速度、湍动能强度比在整个计算区域内采用单一模型的方法能够更好地与实验值吻合。}      

消失球方法及其对旋翼厚度噪声的分析研究

对旋翼厚度噪声及其数值模拟方法进行了讨论,着重解决了以下两个问题：（1）通过理论和数值分析,讨论了使用消失球公式计算旋翼旋转噪声时计算结果的有效性;（2）通过计算机模拟,深入分析了旋翼厚度噪声的特征。     

基于谐波平衡法的含Iwan模型干摩擦振子非线性振动

通过谐波平衡法,研究了由Iwan模型和1个质量块构成的干摩擦振子系统的自由振动和受迫振动问题.对于自由振动,推导了系统卸载和加载过程中各自对应的非线性运动方程;通过谐波平衡法,求得了各方程的1阶谐波解;为了得到方程正确的数值解,提出了确定零速度时刻的策略;数值算例结果表明:谐波平衡法与数值方法得到结果吻合较好,从而验证了谐波平衡法求解此类问题的有效性.对于系统微滑移下的受迫振动,采用谐波平衡法求得了系统的幅频关系.最后,对Iwan模型阻尼特性的研究发现,模型等效黏性阻尼比与其振幅存在非线性关系.      

前后掠风扇叶片颤振特性对比分析

采用基于能量法的流固耦合数值预测方法对比研究了模态、叶间相位角对前后掠风扇转子叶片气动弹性稳定性的影响.结果表明:对于两种掠形叶片,在所研究的前3阶振动模态下,前掠叶片对应于弯扭耦合振型的气动模态阻尼比最大,其数值为0.801%,而后掠叶片对应于弯扭耦合振型的气动模态阻尼比最小,其数值为0.248%;叶间相位角对两种掠形叶片气动模态阻尼比都有显著影响,在1阶振动模态相同叶间相位角(节径)下,前掠叶片前行波对应的气动模态阻尼比小于逆行波对应的气动模态阻尼比,而后掠叶片与此相反;1阶振动模态所有叶间相位角下,前掠叶片比后掠叶片气弹稳定性更好.      

跨声速级不同转速下静叶的损失特性

为讨论不同转速下静叶的气动及损失特性,在09倍设计转速下对某跨声速单级压气机进行了试验测量,获得了级特性曲线并对数值计算方法进行了校核。利用校核后的数值软件对07~11倍设计转速下压气机静叶内部流场进行了详细的数值研究,探讨了转速变化对静叶内部流场结构的影响机制,分析了流场结构与损失分布的关联性,对比了不同损失评价标准并给出了各自对于叶片优化设计的指导意义。结果表明:随着转速的增加静叶顶部流场持续恶化,而叶根和叶中流场变化较小;流道内的旋涡和压力梯度输运了低能流体团,造成了低能流体在流场中的积聚而形成高损失区。      

CT-1标模大迎角静态气动特性数值模拟

 采用雷诺平均的Navier-Stokes方程,对接网格技术和TRIP（TRI sonic Platform）软件,数值模拟了CT-1标模大迎角静态气动特性。来流马赫数0.5,迎角0°～100°,侧滑角0°和5°,在与1.2 m跨声速风洞试验结果对比的基础上,重点讨论了对流项离散方法、湍流模型、远场距离和不同尾部处理形式对CT-1标模大迎角静态气动特性影响。研究结果表明,不同尾部处理形式对大迎角静态气动特性数值模拟结果影响显著,模拟试验的尾部支杆是对比大迎角静态测力试验结果与计算结果的基本要求。     

铁（Ⅲ）—苯芴酮—Tritonx—100体系测定铬若美中的铁

本文研究了铁 (Ⅲ ) -苯芴酮—Tritonx -1 0 0体系测定铬若美中的微量铁。铁与苯芴酮的二元配合物的摩尔吸光系数值可达 6 4× 1 0 4 ·L·mol- 1·cm- 1,加入表面活性剂Tritonx -1 0 0后 ,铁与苯芴酮、Tritonx -1 0 0生成三元配合物 ,摩尔吸光系数值提高到 2 4× 1 0 5·L·mol- 1·cm- 1,Fe含量在 1 4 5～ 2 8μg/2 0ml范围内服从耳定律 ,回收率 97%～ 99% ,反应选择性好 ,用拟定的方法测定铬若美 (铬和镍合金 )中的微量铁 ,测试结果令人满意     

体积成形数值模拟技术的研究现状及发展趋势

 简要回顾了金属体积成形三维数值模拟技术的发展历程 ,对国内外在此领域的研究现状及发展趋势作了较全面的分析与展望 ,指出研究的重点由宏观转向微观、单一分散转向耦合集成、由共性通用转向特性专用型的模拟。介绍了 4种数值模拟新方法——无网格数值模拟、自适应有限元数值模拟、刚塑性有限边界单元法、基于特征的有限元数值模拟。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

计及边界层汇流效应的多段翼型失速特性计算

 本文给出计及边界层汇流效应的多段翼型失速特性的解法。用高阶奇点分布面元法（Panel Method）求位流解,然后进行各翼段粘性尾迹形状迭代,并解出正常边界层和汇流边界层特性。当翼段上有后缘分离时,还要确定分离尾迹的形状,用位移厚度当量源（汇）模拟粘性效应。进行粘/位流迭代直至收敛。上述方法在超过多段翼型失速迎角时仍然有效。计算结果与实验数据比较,符合良好。     

基于笛卡尔网格的非定常气动力计算

在非结构自适应笛卡尔网格基础上构建了一种基于笛卡尔动网格的非定常气动力计算方法。为了适应非定常流场对网格品质要求较高的特点,网格生成采用了基于几何外形、模型表面曲率为基础的自适应生成及加密方法。采用显格式龙格库塔推进方法和中心格式空间离散的有限体积法对NACA0012翼型,三维机翼的跨音速刚性振动进行了数值模拟。对伪时间收敛精度控制,真实时间步长对结果的影响进行了讨论。将计算结果与试验数据进行了比较,得到了良好的效果。     

潜入喷管背壁流场模拟显示

为了研究固体火箭发动机潜入喷管背壁流场的特点,设计了一种冷模拟实验装置.利用引入外加烟雾以及片光照明的显示技术进行背壁区域的流动显示,并用照相机和摄象机记录显示结果.结果表明:在雷诺数较高、潜入度较大和喷管潜入部分型面为凹弧面条件下,背壁区域容易形成旋涡;同时,针对实验条件进行内流场的数值计算,实验和计算结果相比较,二者在定性上具有良好的一致性.     

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

二维喷管两相流动实验理论研究

 一、实验研究 （1）实验装置 图1是实验装置系统示意图。整个系统分为气源,颗粒输送器,过渡段,实验喷管和测量仪器五个部分。实验中气相用压缩空气,颗粒相用铝粉。 气源由两个贮气罐构成,输入端接高压管路,输出端接过渡段。颗粒经由粒子输送轴上的螺旋形通道,被送到过渡段与气相混合。颗粒质量流量通过辅的转速来控制。过     

Identification of vortices in a transonic compressor flow and the stall process

IntroductionVortex is often the essential element re-sponsible for triggering stall and surge in com-pressors.So the identification and the analysisof the vortices in compressor are usually quitehelpful for the understanding of the compressorbehavior.Howe…     

轴向前倾叶片流场的数值分析

对一种新式的风机叶片即轴向前倾叶片流场进行数值计算, 采用时间推进法计算无粘流场, 采用非正交曲线坐标系附面层积分方程求解附面层参数, 得到了一些有意义的结果。      

An inverse code based on a residual-correction method for wing design

A fast, efficient and reliable code for wing design has been developed at INTA coupling a residual-correction method by Bauer, McFadden and Garabedian to an inviscid solver for wing analysis in transonic flow. Smoothing procedures, including Bézier cubic splines, are used to avoid irregularities of the wing surface, as well as the twist distribution. A non-conservative full potential finite differences code is used as the flow solver. The redesign of two wings in transonic flow are presented with promising results. The computing time required to achieve convergence up to a good level of accuracy is less than 1 hour in a Workstation of 300 Mflops. The code is actually a powerful tool for preliminary wing design. Future work will be addressed to the addition of viscous effects and to studies on the optimization of target pressure distributions.