半无限深理想势流中圆柱形结构固有振动的求解

本文利用加权残数法，求解了浸沉半无限深理想势流中圆柱形结构的固有振动。数值算例表明，这种方法，计算量适中，试函数的选取合理，适用于工程中对水下下航行体固有振动的求解。     

后机身收缩方式对通航飞机气动特性的影响研究

通航飞机是民用飞机的重要组成部分,其外形特征对气动阻力具有不可忽视的影响。采用基于 N-S方程的数值模拟方法,分析典型通航飞机阻力极曲线的桶形特征以及机身上表面轮廓线曲率对阻力的影响;研究机身上表面、下表面和侧面 3 种收缩方式下,通航飞机气动特性的差异,并对基础外形和改进外形气动特性进行对比。结果表明：收缩率增大时,...     

水下爆轰燃气泡形态与激波传播过程研究

针对脉冲爆轰发动机在水下工作过程中形成的燃气射流问题,搭建了水下爆轰燃气实验系统,研究了第一个爆轰循环中燃气泡发展变化过程。建立了基于气液两相双流体模型的脉冲爆轰发动机水下喷射模型,采用时-空守恒元和求解元方法,模拟了爆轰波与水相互作用形成的激波的传播及衰减过程。研究结果表明：燃气射流冲击水面时,燃气泡形态在膨胀阶段受到水介质的阻滞作用呈现“豌豆”状,其轴向与径向尺寸有不同的发展规律;同时燃气泡内由于气水冲击作用和燃气扩散受限一直保持较高压力;前导激波传播速度远大于燃气泡发展速度,脱离燃气泡后激波压力值迅速衰减至常压量级,且在中心轴线上的衰减最为剧烈,导致其强度指向性发生改变;前导激波在气液交界面处产生反射,回传的反射激波与后续气流形成的拦截激波碰撞在燃气泡内出现回击现象,使管口附近形成较高的压力峰值。     

基于增广Burgers方程的音爆远场计算及应用

音爆的精确模拟对于超声速客机的低音爆研究与设计意义重大。由于计算能力的限制,客机巡航高度的音爆全场直接模拟目前还难以实现。现有的音爆预测方法一般分为两步,先通过超声速线化理论或计算流体力学的方法得到音爆近场的过压值（Over-pressure）分布,再通过声学理论将近场声压信号推进至远场,以获得飞行器的地面音爆信号。在远场计算中,传统的波形参数法没有考虑音爆传播过程中的经典吸收和分子驰豫效应所造成的声能损失,得到的激波没有厚度,导致计算得到的远场声压级不准确。基于算子分裂法,开展了非线性声学中的增广Burgers方程的数值解法研究。通过计算第二届音爆预测研讨会（SBPW-2）发布的两个标准算例,验证了该方法可以实现地面音爆波形的精确预测。发现在近场声压信号前加入一段无幅值的缓冲信号可以有效提升"N"波上升时间的模拟精度。网格收敛性研究表明适当加密计算网格是有必要的。在此基础上研究了大气声吸收对地面波形的影响,发现分子驰豫效应的影响要强于经典吸收。最后,研究了不同湿度、温度对地面音爆波形的影响,发现干燥、低温的环境对音爆信号的过压值有抑制作用。     

热线探针对数校准方法研究及改进

开展了可压缩流体中热线探针校准方法的研究,以满足其在各种速度测量场合的使用需求。研究了对数校准数学模型,发现校准系数求解过程中存在矩阵奇异性过强的问题,导致在速度小扰动条件下方程求解稳定性差。对对数校准数学模型进行了参数无量纲化及添加正向偏置的改进,建立了无量纲化对数校准数学模型。在马赫数为0.3~0.5,引射压力为150~300 k Pa范围内进行了校准实验,利用对数校准数学模型对实验数据进行拟合,拟合优度为0.997 61,拟合速度平均偏差为1.378 m/s,校准系数求解过程中系数矩阵条件数为1.595×108,矩阵奇异性过强,加入速度小扰动(1 m/s)后,拟合优度为0.379 74,拟合速度平均偏差为43.81 m/s,方程求解稳定性差。利用无量纲化对数校准数学模型对实验数据进行拟合,拟合优度为0.998 95,拟合速度平均偏差为1.203 m/s,校准系数求解过程中系数矩阵条件数为3.655×102,且无量纲化方法不受速度小扰动影响。对流体速度进行不确定度分析,速度平均不确定度为3.168 m/s,无量纲化拟合速度平均偏差明显小于速度平均不确定度。实验结果证明了无量纲化对数校准数学模型应用于可压缩流体热线探针校准的可行性。     

固体燃料冲压发动机内流物性参数自动计算

通常固体燃料冲压发动机内流流动过程的模拟计算中,热力参数及输运系数都视作常数, 实际上, 它们随当地温度、组分等的改变而在发生变化。为确定热力参数和输运系数在各个截面上的变化规律, 研究了热力组分的平衡算法, 提供了简单实用的计算模型。计算结果显示, 把物性参数视作常数会带来较大误差, 说明了考虑物性参数在各个截面上的变化是非常必要的     

盘型磁流变流体阻尼器悬臂转子系统响应时间测试

在不同的磁流变流体、激励方式、转子转速以及磁场强度情况下系统地测量了盘型磁流变流体阻尼器悬臂转子系统的响应时间。结果发现转子系统的响应时间由初始响应时间、快速响应时间和慢变响应时间3部分组成。转子系统的初始响应时间一般小于100ms,快速响应时间一般在50～400ms的范围内,两者均不可能为几个毫秒。初始响应时间及快速响应时间不仅与磁场强度、磁流变流体的组分和转子转速有关,而且还与激励方式有关。突加电流过程中的响应时间比去掉电流过程中的响应时间短。无论是在突加电流还是在去掉电流的过程中,都存在着一个由磁化或去磁过程引起的初始滞后时间。      

计算流体力学的时空观:模型的时空关联性及算法的时空耦合性

流体力学中波的有限传播、粒子的碰撞、各种力之间相互作用,无不体现时空关联效应.本文从计算方法的视角探讨计算流体力学的时空观,即流体力学模型的时空关联性和计算方法的时空耦合性.从流体力学微团法建模出发,明确模型时空关联性的涵义,建立有限体积格式的基本原理,阐述算法时空耦合的必要性,实现流体力学基本控制方程物理建模与有限体...     

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

跨音速翼型的气动最优化设计方法

将计算流体动力学(CFD)与最优化技术相结合,通过数值求解欧拉方程,对翼型绕流流场作出了数值模拟.再利用几何、流动和最优化控制方程,反复迭代求得在一定约束条件下气动性能最优的翼型.本文以NACA0012为原始翼型,选取两种设计工况,都取得了满意的结果.