基于伴随方程的翼型数值优化设计方法研究

对Ｅｕｌｅｒ方程为流动模型的翼型绕流,在微分方程生成的网格中,以Ｂ样条控制顶点为设计变量,用连续伴随方程法对翼型气动外形进行数值优化设计。给出了两个亚声速的反设计算例,验证了该方法的可行性。     

偏微分方程组的一种化简方法

以李群为工具,给出了一种化简两个自变量,两个未知函数的一阶偏微分方程组的方法,在一定条件下,可将偏微分方程组化简为一个阶偏微分方程,若能求出此一阶偏微分方程的解,则只需再解一个一阶常微分方程,即可得出原偏微分方程线的解,该方法可用于某些一阶偏微分方程组的求解,其中包括非定常完全气体一维均熵方程和波动方程。     

固体火箭发动机燃烧室三维流动数值计算

采用k-ε双方程湍流模型,以SIMPLE计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流N-S方程.在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构.计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡.      

低中等雷诺数超声速轴对称射流气动声场

采用线性抛物化稳定性方程(LPSE)方法对低、中等雷诺数条件下的超声速轴对称射流进行了研究。结果表明:射流中不稳定波的增长特性决定了声辐射的特性。在低雷诺数条件下,不稳定波的频率范围较窄,且由于不稳定波的增长区间较长,因此声辐射的范围较宽,声源不集中,声辐射的方向角不明确。在中雷诺数条件下,不稳定波的频率范围较宽,且由于不稳定波的增长区间较短,因此声辐射的范围较窄,声源集中,声辐射的方向角较明确;且随着斯特劳哈尔数增加,声辐射的方向角增加。在低雷诺数条件下,射流基本流的法向速度对射流气动声场的影响较弱,而在中等雷诺数条件下,其对射流的气动声场影响较强。      

基于CE/SE法的翼型流场激波捕捉计算

在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波.     

光纤陀螺分形噪声的预白化滤波方法

光纤陀螺分形噪声往往具有长相关特性,导致陀螺输出信号慢漂移.为了滤除这类分形噪声,提出了一种分数阶预白化差分方法,通过该方法可以将分形噪声转化为高斯白噪声;利用小波变换,在小波变换域中采用Bayes软阈值去噪方法去除该白化噪声,从而形成了一种光纤陀螺分形噪声的预白化滤波方法.对VG949p型光纤陀螺实测数据处理结果表明,该方法能够有效地去除其中的分形噪声,抑制了陀螺低频漂移,同时也滤掉了其高频噪声项.      

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

机翼非平行边界层稳定性研究

从Navier-Stocks(N-S)方程导出曲线坐标系下的抛物化稳定性方程(Parabolicstabilityequation,PSE),研究机翼非平行的可压缩边界层稳定性问题。发展了求解PSE的高效数值方法:引进法向变换,使得在临界层与壁面之间的扰动量变化最快的区域有更多法向网格点;采用包含边界邻域在内的完全四阶精度的法向差分格式,这对方程精确离散至关重要;以及全局法和局部法相结合的数值方法及其新的迭代公式,能大大加速收敛并得到更精确的特征值。算例分析研究了扰动增长因子和形状函数等演化曲线。     

基于结构网格点的Euler方程无网格自适应求解算法

提出了一种自适应无网格Euler方程求解算法。构造了基于无网格点云结构的布点加密技术,实现了借助压力梯度变化准则的流场局部加密;借助流场结构化剖分技术,给出初始无网格节点分布,并采用四步显式Runge—Kutta时间格式推进计算,求解了绕翼型的典型流动问题。数值算例表明,本文提出的方法能有效地提高流场的分辨率,如捕捉激波等流动特征清晰明了。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。