跨声速非等熵势函数方程的数值求解

基于动量方程的特点,本文将只与温度有关的“等熵密度”与熵的变化分离,而由散度型动量方程直接计算出气体经过激波的熵增。与经典势函数方程迭代,可以方便地求解非等熵跨声速流动。从对任意回转面叶栅的跨声速流动计算来看,这种非等熵势函数得到的激波位置可比等熵势函数前移1～2个网格,强度有所减弱,和试验结果较接近。由于各流线的熵不相等,对Kutta条件和出口边界条件作了修正。本文建议的方法也可用于计算流函数的熵增。     

叶轮机械全三维气动杂交问题的数值求解方法

本文介绍了一种求解叶轮机械全三维气动杂交问题新的数值方法。本方法使用了非正交曲线座标系和非正交速度分量的势函数方程, 并构造了守恒形式的近似分解因式格式, 通过不断修正叶片通道宽度, 使叶面上速度达到预先要求的速度分布。根据本方法使用 FORTRANTT编写了计算程序, 并对单转子压气机叶片进行了修正计算, 其结果十分令人满意。      

逼近无控旋转目标航天器的混合势函数安全制导

研究了追踪航天器逼近无控旋转目标航天器的安全制导问题，逼近过程中，追踪航天器需要躲避空间中的障碍物，同时需要避免与目标航天器的太阳能电池帆板和天线等附件发生碰撞。建立了视线坐标系下的两航天器间的相对运动方程，采用四元数描述目标航天器的姿态运动。将参考位置设为引力源，设计了吸引势函数。针对安全逼近问题，建立了球面安全区和锥面安全走廊，设计了安全势函数。将障碍物假设为具有一定半径的球体，设计了障碍物势函数。吸引势函数、安全势函数和障碍物势函数一起组成了混合势函数。为了解决整个势场中除参考位置外还可能存在其他局部极小点问题，对混合势函数进行了修正，保证参考位置位于混合势函数的最低点。利用Lyapunov稳定性理论对混合势函数进行了稳定性分析，推得符合要求的控制加速度，使追踪航天器沿着混合势函数的负梯度方向逼近无控旋转目标航天器。最后通过数值仿真验证了该方法的有效性。     

再入气动环境类电池帆板材料微观响应变形行为分子动力学模拟研究

针对服役期满大型航天器无控飞行轨道衰降再入大气层解体过程及落区难以提前预测,再入解体后生成的碎片可能造成地面危害等问题,采用分子动力学模拟方法,选取MEAM势函数,构建了碳元素质量分数0.215%的含碳钢分子动力学模型,计算了不同温度下材料的平衡态晶格常数,并通过晶格常数与温度的关系,计算了模型的线膨胀系数,验证了MEAM势函数在所建立的仿真模拟系统合理性;使用经过验证的分子动力学模型与MEAM势函数,模拟了钢制平板在Ma_∞=8.37,Kn_∞=0.01,γ=1.4的高超声速再入气动环境中,通过结构动态热力响应变形行为有限元算法计算出的部分状态下材料微观演化行为,初步证明了分子动力学模拟方法在服役期满大型航天器再入大气层解体过程的分析仿真计算中的作用,为实现分子动力学方法同动态热力响应有限元算法的耦合奠定了基础。     

快速CFD计算工具在民机概念优化设计中的应用

高亚声速民机概念设计和优化阶段需进行大量的反复迭代,其中气动力计算的准确性与时效性是一对很难处理的矛盾体。Navier-Stokes(N-S)方程可进行较为精确的计算,但需要付出的时间和资源代价很大;工程估算方法可以快速得到估算结果,但精度较差,尤其表现在跨声速阶段的波阻估算。本文采用基于全速势方程的快速CFD计算工具与工程估算相结合的气动计算方法,避免了N-S方程耗时长和工程估算精度差的问题,实现了在民机概念设计和优化程序中的集成和应用。     

非结构网格生成技术方法研究

利用点源或线源所影响的势函数分布规律,通过建立势函数与网格背景信息之间的关系,提供非结构网格生成过程中所需的网格背景参数;在物体边界节点已知的情况下,适当布置点源和线源个数,用最小二乘法确定各点源和线源的强度;通过对部分点源或线源的影响势函数施加方向控制开关,使所生成的网格在期望的区域得到加密;算例表明本文方法是一种有效、快速的处理方法,与传统用背景网格提供背景信息的方法相比,在计算量、程序量、光滑性和灵活性方面都能表现出良好的性质。     

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

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

改进格网型中的跨音速分析和设计计算

 <正> 近年来采用不同的主控方程计算跨音速绕流的方法有了不断的发展和改进,越来越多的算例表明了计算流体力学的迅速发展,其中解全速势方程的跨音速计算方法也有了很大进展,并已为国外航空工业界广泛地应用于设计和分析计算中。其中有限体积方法提高了处理复杂外形的能力,若采用多重网格技术则可改进计算效率。但是为了发展一种普遍适用和可靠的格网系统以耦合有限体积算法仍有许多工作要做     

时间相关法计算喷管跨音速流动

本文用时间相关法完成了定常、无粘跨音速喷管流场的计算,内部点参数用Mac Cormack差分格式计算、壁边界点参数用简化了的特征边界条件计算,轴上点参数用反射原理计算。计算表明:计算是收敛的,计算结果与实验数据符合良好。 定常、无粘、跨音速喷管流动的控制方程是混合型的(即在亚音速区域是椭圆型,而在超音速区是双曲型的),给数值求解带来很大的困难。为了克服这一困难,广泛采用时间相关法,即认为定常流动方程的解是相应的非定常流方程的解在时间趋于无穷时的渐近解。因为非定常流的控制方程不论在亚音速区,还是在超音速区域,都是双曲型的,可以用统一的方法来求解,而且易于求解;另一方面,它可适用于形状比较复杂的喷管。 本文首先阐述了时间相关法的计算方法,然后列出算例喷管的计算结果。通过与实验数据的比较,证明计算是符合实际的。     

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

IV: GEODESY: Remarks on the Role of Height Datum in Altimetry-Gravimetry Boundary-Value Problems

The problem of global geoid determination is usually solved using satellite altimetry data on the oceans, together with an oceanographic model of sea surface topography, and gravity anomaly data on the continents. Such data, however, enable to obtain only potential differences with respect to a reference surface whose absolute potential is unknown. This situation suggests to modify the classical mixed boundary-value problem of physical geodesy by inserting into the boundary conditions an unknown additive constant, that must be determined by imposing a suitable additional constraint. Yet, such formulation of the boundary-value problem, from the point of view of its mathematical properties, is not unconditionally well-posed, and, furthermore, does not reflect faithfully the available physical model, as the present knowledge of ocean circulation does not allow to connect along coastlines the reference surfaces defined on the oceans and on the continents. The introduction of two different unknown additive constants, one for the oceans and one for the earth, to be determined by imposing two additional constraints, gives rise to a more faithful picture of the present physical knowledge, and, at the same time, to a new well-posed formulation of the boundary-value problem. This revised version was published online in August 2006 with corrections to the Cover Date.     

一种求解涡轮平面叶栅的跨音松驰法

一、引　　言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法…     

改进跨音速位势方法的讨论

 简要地回顾并讨论了传统位势方法之不足。对改进传统位势方法的两种必要的修正(熵修正和涡修正)的量级在小扰动条件下作了分析和估值。结果表明涡修正是更高阶小量,因而作为第一步修正可采用非等熵位势流来近似。采用熵修正激波算子来计及非等熵的激波突跃条件,并考虑了由熵修正带来的尾迹条件的变化。算例表明非等熵位势方法能得到预期的改进位势方法的效果。     

高阶面元法计算多段翼型之位流

本文给出了一种求解多段翼型位流的高阶面元法。面元和奇点的分布分别采用了二阶和三阶的形式。由于高阶面元法较为精确地反映了绕流体的几何形状和奇点表面分布,本文给出了比低阶面元法更为精确的解,并且消除了低阶面元法计算中经常出现的翼型后缘处解的波动。     

计算机翼跨音速绕流的加罚Galerkin有限元法

1.引言 全速位Galerkin有限元法是跨音速有限元法中重要方法之一。处理超临界流时,为了自动捕捉到激波,通常采用上风技术,但因此导致了质量守恒条件得不到完全满足。文献[1]从Bateman变分原理出发,应用加罚函数方法,找出了强制满足质量守恒条件在Galerkin有限元方程中所应加入的附加项。本文将文献[1]的方法应用到超临界跨音速绕机翼流动,克服了一般Galerkin有限元法中由于采用上风技术而出现的激波不够陡和位     

电力作动器中永磁容错电机的电感和谐波分析

 针对永磁容错电机大漏感的特点,利用磁矢量有限元法对六相十极永磁容错电机进行分析,得出其绕组自感的组成成分,尤其是电机内漏感分布;同时建立六相十极永磁容错电机的磁路模型,利用磁路法推导出较高精度的绕组自感及其组成部分的解析式,通过有限元验证,该解析式的计算精度在5.6%以内。此解析式可推广到一般的永磁容错电机中去,进而为永磁容错电机的初始设计和性能分析提供了理论依据。最后,对永磁容错电机内谐波漏感进行了量化分析,为后续的电机控制及其热分析打下基础。     

On the Treatment of Inertial Measurement Unit Data in Optimal Linear Navigation Policies

The inclusion of the integrated nongravitational acceleration outputs from an inertial measurement unit in an optimal linear navigation policy is treated. This formulation appears to be considerably different nt from the treatments of IMU data for space vehicle navigation that have previously been discussed and promises to yield considerably better results in many instances. It has the advantage that the data are processed in an optimal fashion that makes the best use of knowledge of the general system characteristics. Further, the IMU enters the formulation in the same manner as any other sensor and, therefore, provides a more uniform method of treatment. Two methods of dealing ng with the IMU are discussed and then applied to a simple example. The advantages of the second formulation are observed to be quite significant.     

颗粒冲刷条件下三元乙丙绝热材料烧蚀特性实验

为了研究三元乙丙基础配方、抗过载配方、无SiO2配方和无纤维配方等不同配方绝热材料在颗粒冲刷条件下的烧蚀特性和机理,基于高过载地面模拟试验发动机,开展了烧蚀试验研究。分析了实验后绝热材料的表面以及断面形貌特征和宏观烧蚀率结果,探讨了三元乙丙基础配方和抗过载配方绝热材料炭化层中致密结构形成的机理,进行了不同配方材料的烧蚀模式和机理的分析。研究结果表明:相同颗粒冲刷条件下,不同配方三元乙丙材料烧蚀模式和破坏机理不同,烧蚀率从小到大依次为抗过载配方基础配方无SiO2配方无纤维配方;其中基础配方和抗过载配方三元乙丙材料在烧蚀过程中炭化层形成"致密/疏松"结构,表面形成的致密结构可以提高材料的抗烧蚀性能,无SiO2和无纤维三元乙丙材料以基体的粒状和块状剥落破坏为主;不同配方绝热层炭化层的形貌和结构与纤维、SiO2含量以及烧蚀环境密切相关。     

丁羟推进剂燃速计算图象法

本文用PEM模型编制的燃速计算程序,计算了一定压力,一定过氯酸铵(AP)粒径级配条件下不同配方的燃速,求出配方组分与燃速关系方程,绘制出AP/Al/HTPB系统的等燃速三角图,在图中由配方组分可迅速查找出燃速,也可由给定的燃速选择合适的配方.     

Micro-mechanical model for predicting the elasto-plastic behavior of composites based on secant formulation method

To predict the Elasto-Plastic Behaviors (EPBs) of aligned inclusions reinforced composites, this paper develops an interpolative Mori-Tanaka/Double-Inclusion (MT-DI) homogenization model with the secant formulation, and gives the numerical implementation algorithms of the developed MT-DI model with the secant formulation. The Finite Element (FE) homogenization method is implemented to provide the “exact” EPBs of the composites and thus validate the MT-DI model with the secant formulation. The MT-DI model with the 2nd-order secant formulation is validated to provide the more accurate predictions, while the MT-DI model with the 1st-order secant formulation always gives the stiffer predictions. The results show that using the macro-stress and macro-strain as the inputs, the MT-DI model with the secant formulation gives the identical predictions. The predictions of the MT-DI model with the secant formulation vary between those of the MT and DI models with the secant formulation. Meanwhile, the MT-DI model with the secant formulation does not predict the accurate EPBs for the phases of the composites.