自适应与双目标优化机翼的气动特性比较

在二维翼型自适应的研究基础上,用Powell法优化计算了对前缘后掠角为35°,展弦比为3.9,梢根比为0.17,机翼剖面为NACA65006翼型的梯形翼的前后缘舵面偏转角,从而获得了在亚跨声速时升阻比大而在超声速时阻力系数小的自适应机翼的最优气动外形.采用了并行遗传算法,计算了要求亚跨声速升阻比大同时超声速阻力小的气动双目标优化机翼的外形.讨论了优化机翼相对于原始机翼的气动增益.与二维一样,三维数值算例也证明了自适应机翼可获得明显的气动增益.      

探针式光纤探头的有效检测深度

研究了探针式光纤探头的有效检测深度。采用扩散方程描述光在强散射生物组织中的传播规律,运用纽曼边界条件,使用灵敏度矩阵对有效检测深度进行了研究。灵敏度矩阵和光源与检测器的空间位置分布密切相关。光源与检测器对的灵敏度矩阵定义为：光源位置和检测器位置分别作为光源的光场分布的乘积。设计了6组针对不同参数的模型,参数主要包括光源、检测器光纤芯径、光源、检测器光纤中心距、光纤探头所处的深度、生物组织的吸收系数和散射系数。对探针式光纤探头在各组模型中的有效探测深度进行了仿真,重点研究了六种参数对检测深度的影响。所得结论可以推广到分析光纤检测器的空间分辨率和多光纤检测器的研究中。     

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。     

弱撞击对接机构动力学特性建模

针对空间飞行器捕获对接过程中复杂的多体动力学问题,鉴于弱撞击对接机构(Low impact docking mechanism,LIDM)不同于以往常见的刚性对接机构,其在对接过程中产生碰撞力远小于空间飞行器自身的重量,且自身对接环处安装有六维力传感器,因此可以将捕获过程中LIDM处产生的碰撞力等效为作用于对接环质心处...     

混合式CRP吊舱推进器水动力性能数值模拟及试验

建立混合式CRP吊舱推进器数值模型并置于数值空泡水洞内,结合RANS方程和SST k-ω湍流模型,运用滑移网格方法对混合式CRP吊舱推进器在均匀流场中水动力性能进行了数值模拟,得到了敞水性能预报结果,并分析了混合式CRP推进器的非定常力;同时,利用空泡水洞、吊舱动力仪及长轴动力仪对混合式CRP吊舱推进器敞水性能进行了试验研究,通过3组重复性试验表明了试验方法的稳定性;并将3组试验平均值与将数值模拟所得的敞水性能结果进行了比较,彼此吻合较好。结果表明,本文所用试验方法和数值计算方法对吊舱推进器模型敞水特征的预报具有较高的可信度,能达到工程应用要求的要求;和国外学者在空泡水洞内得到的不附带吊舱相同对转桨敞水试验数据进行比较,分析了吊舱对前后桨敞水性能的影响。结果表明,吊舱对前桨的敞水性能基本没有影响,但对后桨的敞水性能影响显著。     

离散速度方法及其在跨流域问题中的应用研究进展

往返大气层跨流域飞行器的气动力和气动热计算一直是当前流体力学研究的难点和热点之一。由于该类飞行器面临的流动场景不再是单一的连续流或稀薄流,采用Navier-Stokes方程求解器或DSMC方法均不能获得全流域的准确结果。近年来,以不依赖于连续性假设的Boltzmann模型方程为基础,通过在位置空间和速度空间同时离散求解该方程,实现了跨流域问题的统一求解。本文对该类算法的研究进展进行回顾和分析,着重介绍气体动理论统一算法（Gas kinetic unified algorithm, GKUA）、统一气体动理学格式（Unified gas kinetic scheme, UGKS）和改进离散速度方法（Improved discrete velocity method, IDVM）3种数值途径,分析它们的基本假设和实现方式,关注它们目前取得的进展和应用情况。同时,本文还将IDVM进一步扩展到非定常情形,以便用于非定常跨流域问题求解。最后,本文对该类算法存在的一些问题进行讨论,期望在后续的研究中能予以解决。     

求解代数Riccati方程的一种新算法

本文讨论了Hamiltonian矩阵在辛相似变换下的标准形,由此给出代数Riccati方程存在非负定解的一个充分条件,提出了求解代数Riccati方程的一种新的算法。该方法节省运算量,且精度较高,尤其适合于解阶数不太高、系数矩阵为满阵的Riccati方程。最后给出了一个数值例子,并将该方法与其他方法作了比较.     

李雅普诺夫方程AX-XB=C的误差界

本文对李雅普诺夫方程AX-XB=C及(A+δA)(X+δX)-(X+δX)(B+δB)=C+δC的解X及X+δX给出了相对误差‖δx‖/‖X‖的一个上界。     

双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

串列叶栅粘性紊流流场有限分析数值研究

本文将近年发展起来的有限分析方法(Finite Analytic Method)用于曲线坐标系上紊流N-S方程的数值计算。计算了来流雷诺数为2.0×10~5,气流攻角分别为0°,10°,-10°三种情况的二维粘性紊流串列叶栅流场。文中用k—ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。数值计算结果与实验结果相比较,吻合程度令人满意。