基于三维非结构网格的连续伴随优化方法

基于非结构网格,采用Euler方程、连续伴随方法以及自由形面变形技术(Free-form deformation,FFD)结合序列二次规划(Sequential quadratic programming,SQP)优化算法,开发出了三维飞行器的气动外形优化代码。应用Fortran 90语言,动态存储数组和指针,节省了内存。采用Roe形式的二阶迎风格式和LU-SGS隐式时间离散格式求解连续伴随方程,同时采用一种新的参数化坐标求解方法,使得FFD技术更加鲁棒。和有限差分方法求得的梯度进行比较,验证了连续伴随方法所求梯度的正确性,实现了ONERA M6机翼的减阻优化,证明本文开发的优化代码能够实现三维气动外形的减阻优化。     

基于风洞试验的非定常气动力微分方程建模方法

为准确描述大攻角非定常运动中非定常气动力的气动特性,研究了以微分方程为基础的非定常气动建模方法。在某战斗机大振幅偏航、滚转单自由度及偏航-滚转耦合运动风洞试验的基础上,对偏航力矩、滚转力矩建立微分方程模型,分析研究不同频率、同一侧滑角情况下,气动力的迟滞时间;然后对高频运动下的气动力进行拟合,进而建立高频运动下的气动力微分方程模型。研究表明:模型可精确预测不同机动下非定常气动特性;建模方法具有较强的工程可行性。     

多冲量近程交会最优化模型

文章基于Lawden方程对椭圆参考轨道的近程最优交会问题进行了研究,并提出了一种混合遗传算法求解最优近程交会问题。首先在一定假设条件下给出了目标在椭圆参考轨道的近距离相对运动模型——Lawden方程,构建了多脉冲最优交会问题模型并进行了理论分析。性能指标选为轨道交会过程中燃料消耗和时间消耗加权最小的多目标优化指标,优化参数为脉冲大小和脉冲施加时刻,终端状态受到相对位置和相对速度的约束。然后介绍了具有较强全局和局部寻优能力的混合遗传算法。最后以四脉冲为例进行仿真计算。仿真结果表明,是否考虑第一次脉冲位置,总燃料消耗变化不明显。因此,追踪航天器一旦捕获到目标信息即可施加第一次脉冲。仿真结果还证明了混合遗传算法在求解最优交会问题时的有效性。因此,混合遗传算法对基于Lawden方程的椭圆参考轨道近程最优交会问题的求解可行。     

不同反压下椭圆形隔离段流场特征与气动性能

在反压条件下,对椭圆形等直隔离段进行了全三维数值仿真研究.结果表明:在对称的几何条件和边界条件以及均匀来流条件下,当来流马赫数较小时,流场与激波串结构基本对称;当来流马赫数较大(大于2)时,流场与激波串结构出现明显不对称现象,并且随着来流马赫数增大,这种不对称现象加剧.椭圆形隔离段抗反压能力随着长短半轴比减小逐渐增强,出口马赫数和温升比基本不受长短半轴比影响.不同来流马赫数下,当隔离段承受最大反压时,出口马赫数变化不大.      

Oscillatory thermocapillary flows in simulated floating-zones with time-dependent temperature boundary conditions

This study deals with numerical simulations of the Maxus sounding rocket experiment on oscillatory Marangoni convection in liquid bridges. The problem is investigated through direct numerical solution of the non-linear, time-dependent, three-dimensional Navier-Stokes equations. In particular, a liquid bridge of silicon oil 2[cs] with a length L = 20 [mm] and a diameter D = 20 (mm) is considered. A temperature difference ΔT = 30 [K] is imposed between the supporting disks, by heating the top disk and cooling the bottom one with different rates of ramping. The results show that the oscillatory flow starts as an ‘axially running wave’, but after a transient time the instability is described by the dynamic model of a ‘standing wave’, with an azimuthal spatial distribution corresponding to m = 1 (where m is the critical wave number). After the transition, the disturbances become larger and the azimuthal velocity plays a more important role and the oscillatory field is characterized by a travelling wave. The characteristic times for the onset of the different flow regimes are computed for different rates of ramping.     

HIERARCHICAL ITERATION METHOD TO SOLVE DYNAMIC DIFFERENTIAL EQUATION

1CauseoftheProblemTherearemultitalysystemsindynamicanalysswherethemassandinertiaOfthemaintheyaremuchbiggerthanthoseoftheappendagds,buttheactionofasmallapPendagecanchangethecharacterofmOionofthemaintalydramaticallylforexample,amultitheysystemlike8pin-stdssateiliteanditspassivenutatindamPer.ThesmalldamPerhastheabili-tytochangethesateilitenutationmotbofromanunstablestatustoastableone.Generally,theenergydissipationmethodcanbeusedtoanalyzethiskindOfdampersys-tem.Unfortunately,ifthesystemdynamicr…     

CYCLIC VISCOPLASTIC DESCRIPTION FOR THE MEAN STRAIN EFFECT

ＣＹＣＬＩＣＶＩＳＣＯＰＬＡＳＴＩＣＤＥＳＣＲＩＰＴＩＯＮＦＯＲＴＨＥＭＥＡＮＳＴＲＡＩＮＥＦＦＥＣＴＣＹＣＬＩＣＶＩＳＣＯＰＬＡＳＴＩＣＤＥＳＣＲＩＰＴＩＯＮＦＯＲＴＨＥＭＥＡＮＳＴＲＡＩＮＥＦＦＥＣＴ￥ＹａｎｇＸｉａｎｊｉｅ（ＳｏｕｔｈｗｅｓｔＪ...     

物面突变引起的二维分离流动的数值模拟

本文用一个单步二阶精度格式,对由物体截面突变引起的二维超声速分离流动进行了数值模拟。计算了近尾迹流场,并详细研究了来流M数和Re数对绕后台阶分离流场结构的影响,得到了与理论和前人计算一致的结果。另外,对超声速绕竖舨流动进行了数值计算,给出了细致的流场结构图画。     

超声速有翼飞行器的无粘分区数值模拟

本文采用ＭａｃＣｏｒｍａｃｋ两步显式格式，用空间推进法求解Ｅｕｌｅｒ方程来对有翼飞行器的超声速绕流流场进行数值模拟，应用特征边条件处理场面边界及装配激波的计算，在翼缘及表面斜率不连续处，采用局部解析处理，提高了算法的可靠性，复杂外形的计算采用分区处理技术，网格由代数形式生成，在大攻角的计算中，通过对横流速度型的修正改善了背风面压力分布特性，解决了推进中断问题，同时模拟出背风面涡。     

旋转坐标系N-S方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明：所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。