圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

湍流通用物理方程的推导及讨论

湍流的通用方程已经获得。利用变形率张量与摄动变形位移的点积分解湍流脉动速度,并使用密度加权的侧偏系综平均方法取代Reynolds平均方法,对Navier-Stokes方程进行平均。在二阶精度上没有出现未知相关量,从而获得了具有二阶精度的封闭的湍流通用物理方程。方程中无相关量,不引入任何经验系数,无需构造任何湍流模型,即可描述任意复杂湍流场并解释湍流物理本质。方程在长时标上描述湍流平均行为,在短时标上描述相应尺度层次的湍流结构现象,从而统一了湍流的统计及结构二大学派。方程简化后得到湍流的规范方程,是混合的Burgers-Korteweg-de Vries方程,对该方程的分析给出了湍流转捩准则。文中对湍流的间歇性、拟序结构、各向异性、湍流能量逆转等现象的物理本质进行了初步讨论。      

地月系低能转移轨道设计

分析了不同转移轨道间的性能指标,以黄道面与白道面交线作为庞加莱截面实现不同系统轨道的拼接,利用其中耗能最小的日地系稳定流形和地月系的不稳定流形设计了一条低能耗转移轨道.最后对不同轨 道进行了能量分析和对比,优化了国内提出的低能探月轨道.仿真结果表明,文中方法比直接通过霍曼转移节约21.1%的能量,对探月工程的轨道设计具...     

使用GAO-YONG方程组对湍流边界层拟序结构特性的研究

与以往的湍流模型不同,GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。其级数形式的能量方程与非线性现象多尺度层次相对应,具备了描述湍流统计平均和拟序结构的双重功能。本文采用交错网格和SIMPLE方法求解GAO-YONG不可压湍流方程组,考察了GAO-YONG方程组中能量方程阶数以及计算网格对拟序结构的影响,证明了该方程在描述拟序结构方面具有一定的适定性。      

二元翼型结冰数值模拟研究

对二元翼型前缘结冰数值模拟进行了研究。采用中心有限体积法求解N—S方程;四阶龙格一库塔法求解水滴运动轨迹方程;逐次二分法寻找水滴撞击极限,然后计算局部水滴收集率;分析翼型表面的传质、传热项,建立相应质量、能量方程并迭代求解,得到翼型表面结冰质量;根据角度二分法生成结冰冰形。预测的NACA0012结冰冰形与实验及Lewice模拟结果吻合很好,结冰过程中的特征参数与冰型及相应冰形特征吻合。     

跨声速三维Euler方程的显式／隐式格式有限元解法

本文采用有限体积Ｇａｌｅｒｋｉｎ法和非结构网格数值求解了跨声速三维Ｅｕｌｅｒ方程,通量的求解采用ＶａｎＬｅｅｒ－阶逆风矢能量发裂方式,并通过外插使其上升为二阶精度,分别采用了显式及隐式两种格式求解。     

直升机发动机燃油泵自吸油地面模拟试验中用正、负压模拟燃油过载的研究

本文推导了在势函数π=fz的有势力场中一维绝热定常流的能量方程,并对该方程进行分析,提出了在直升机燃油系统地面模拟试验中用正、负压模拟从油箱到发动机吸油口管路中燃油的过载来进行在各种飞行高度、各种供油状态下的发动机燃油泵自吸油地面模拟试验。     

Gao-Yong理性湍流方程

对近十年来Gao-Yong理性湍流方程的研究进行综述,给出矢量形式的可压湍流封闭方程组的推导过程。侧偏统计平均保留了湍流脉动量的一阶统计信息,在引入加权漂移速度对称性及正交各向异性后,导出了漂移流的连续方程、动量方程及机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组。方程组不含任何经验系数,不使用壁面函数,保留了NS方程的均化的非线性特性。其级数形式的能量方程与非线性现象多尺度层次现象相对应,具备了描述湍流统计平均流动及拟序结构流动的双重功能。已进行了大量算例验证,验证结果证明了Gao-Yong湍流方程对广泛范围的复杂湍流问题的适定性。     

等谱AKNS方程的约化

通过对等谱AKNS方程的约化,构造出非线性Schrodinger和mKdV方程的新双Wronskian解;分别推导出这2个方程的双Wronskian形式的有理解。     

基于Lagrange方法的单旋翼飞行器动力学建模

为了进行微小型单旋翼飞行器的动力学建模,通过建立多个坐标系来反映各部分间的相对运动.首先,利用坐标变换得到位置、速度及加速度等向量,并代入拉格朗日方程得到运动学模型;然后,对模型进行数值求解,得到飞行器的姿态响应.仿真结果表明,飞行器定点盘旋时合外力为零,能量保持不变;爬升或前飞时有非保守力做正功,能量增大.     

某型涡轴发动机等压差活门建模分析

考虑某型涡轴发动机燃油流量调节器通道截面的局部能量损失,推导出流经计量油针的燃油流量方程,在此基础上详细分析等压差活门的功用、结构组成及其工作原理,以连续方程和力平衡方程为基础,采用线性化处理方法,建立等压差活门的数学模型。通过分析等压差活门系统结构,得到等压差活门系统的传递函数,进而对其稳定性进行分析,计算单位阶跃输入的稳态误差。结果表明:以流量方程和力平衡方程为基础建立等压差活门数学模型的方法可行,分析出等压差活门稳定工作条件和稳态误差的主要影响因素,并且为燃油调节器数学模型的建立以及仿真计算奠定了基础。     

任意回转面亚、跨音流场计算方法

本文所设计的任意回转面亚、跨音流场的计算软件 ,能根据叶型和进出口条件分析流动特性 ,并自动选择精确和节省机时的合适的数值解法。用户在给出叶型坐标和定解条件后 ,可以得到叶面马赫数分布或压力分布 ,全流场的等马赫数或等压图     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

四边简支压电热弹性体层合开口壳的精确解

根据压电热弹性体的本构关系和热传导关系的相似性,通过增维方式将热传导关系写进压电热弹性体的本构关系中,直接得到了一个新的增维的压电热弹性体本构关系。依据此增维的本构关系,结合压电热弹性体平衡方程和热平衡方程,成功地导出了压电热弹性体层合开口壳机、电、热耦合问题的齐次状态方程。这类方法大大简化了人们在分析压电热弹性体耦合问题时通常要求解的非齐次微分方程和关于温度的平衡方程和导热方程的二阶微分方程的繁琐方法,同时,齐次状态方程的提出大大减少了数值计算的工作量。     

迷宫密封的理论与实验研究

 根据相似原理。通过实验回归分析法,得出压力损失系数ζ与齿腔几何特性和流动特性的关系式。然后应用能量定律、连续方程、状态方程和Bernoulli方程与其联立,求解后可获得迷宫密封的全部特性。并可归纳成流量系数α、穿越系数β和迷宫函数ψ对全压比λ的特性曲线。应用这类曲线,可方便地计算出迷宫密封的泄漏率。     

弹流润滑螺旋锥齿轮热摩擦行为分析

针对直升机主减速器中的高速重载螺旋锥齿轮,建立了点接触热弹流润滑分析数学模型,包括Reynolds方程、能量方程和膜厚方程等,采用数值方法求解弹流润滑状态下的齿面摩擦因数.模型中考虑了润滑油黏度和密度随压力、温度的变化,并通过轮齿承载接触分析(LTCA)获得齿面真实载荷和卷吸速度、相对滑动速度等运动学变量.基于真实齿面点建立了螺旋锥齿轮单齿模型,考虑滑动摩擦生热和不同表面上的热边界条件,通过有限元稳态热分析和瞬态热分析得到了轮齿本体温度场和接触点瞬时闪温,并与现有文献和算例齿轮台架试验结果进行对比.      

喷液式固体火箭发动机的稳定性

提出了一种喷液式推力大小可控的固体火箭发动机方案。对反馈控制系统的各个环节建立了运动方程,进而求得时变非线性系统在平衡点的常系数拟线性控制方程,并确定了系统的平衡稳定边界方程和稳定准则。分析指出,喷液式推力可控固体发动机和普通固体发动机在稳定性上有本质差别。小推力工况下发动机工作初期稳定性最差。喷注比或燃烧室自由容积越大,稳定性越好,压强指数在中n和高n值时有利于稳定。     

Numerical investigation of transcritical liquid film cooling in a methane / oxygen rocket engine

Transcritical film cooling was investigated by numerical study in a methane cooled methane/oxygen rocket engine.The respective time-averaged Navier-Stokes equations have been solved for the compressible steady three-dimensional(3-D) flow.The flow field computations were performed using the semi-implicit method for pressure linked equation(SIMPLE) algorithm on several blocks of nonuniform collocated grid.The calculation was conducted over a pressure range of 202 650.0 Pa to 1.2×107 Pa and a temperature range of 120.0 K to 3 568.0 K.Twenty-nine different cases were simulated to calculate the impact of different factors.The results show that mass flow rate,length,diameter,number and diffused or convergence of film jet channel,injection angle and jet array arrangements have great impact on transcritical film cooling effectiveness.Furthermore,shape of the jet holes and jet and crossflow turbulence also affect the wall temperature distribution.Two rows of film arranged in different axial angles and staggered arrangement were proposed as new liquid film arrangement.Different radial angles have impact on the film cooling effectiveness in two row-jets cooled cases.The case of in-line and staggered arrangement are almost the same in the region before the second row of jets,but a staggered arrangement has a higher film cooling effectiveness from the second row of jets.     

电子设备冷却用液体冷板表面温度场的一种工程算法

本文提出折合厚度概念,把冷板简化成一个带内热源、导热系统各异的三维稳态导热模型实体,根据能量平衡方程,建立一组适于工程应用的三维非均匀热流冷板数学模型,进行数值求解,计算结果与试验结果比较表明本文数学模型和计算方法是可行的。     

二次空气系统一维非定常计算方法

研究了燃气轮机二次空气系统的非定常计算方法.通过将整个二次空气系统划分为两部分,特征"长"元件以及局部损失元件.特征线法用于单个特征"长"元件的求解,通过将连续方程、动量方程、能量方程构成的偏微分方程组转换为常微分方程组并联立求解,可以充分考虑内部流动与换热的相互影响;基于压力修正的流体网络法则用于由多个局部损失元件组成的局部流体网络,首先对动量方程组进行求解,进而对各自的能量方程进行求解.如何在采用这两种方法进行求解的上述两类元件之间进行数据交互将被讨论.此外,所提出的时间步长定义方法以及改进的元件计算模型使得二次空气系统的非定常分析得以实现.计算程序经与某系统定常实验数据及单个元件的非定常计算数据进行对比,在10s的计算时间内,基于特征线法得到结果与CFD数据偏差达到10.2%,结果表明该计算程序可以满足二次空气系统非定常计算的要求.