基于GPU和隐式格式的CFD并行计算方法

从图形处理器(GPU)架构特点出发,提出了基于数据并行的隐式计算流体力学(CFD)求解方法,空间离散格式采用迎风Roe格式,计算网格适用于结构和非结构网格。采用统一计算设备架构(CUDA)技术实现了GPU上的隐式CFD并行计算。分别在Intel Core2 Quad3.0GHzCPU和NVIDIAGTX280 GPU上进行了计算,结果表明隐式格式计算速度是显式格式6倍以上,采用显式格式的计算加速比达到28倍,采用隐式格式计算加速比达到了28.7倍,同时计算加速比随计算规模的增加而增加。计算结果和实验结果较为吻合。     

重叠网格洞面优化技术的改进与应用

针对传统割补法在壁面狭缝处网格重叠失败的问题,提出了以单元顶点状态进行洞面优化的vertex方法及以单元大小进行洞面优化的volume方法的两种改进方法。这两种方法将传统割补法中切割和填补两步搜索方法改为一步搜索方法进行重叠网格的洞面优化。通过数值试验证明:两种洞面优化方法都能很好地解决壁面狭缝网格重叠问题,网格重叠区域流场数值传递正确,计算结果与试验值吻合,优化效率高,收敛速度快。volume优化方法与vertex优化方法相比可以减少由网格大小不对等引起的数值误差,计算结果和计算精度均优于vertex优化方法。     

火星大气环境模拟装置设计及仿真分析研究

对火星表面大气环境特性进行了研究,通过选取合适的计算方法并结合FLUENT流体有限元计算软件对火星表面稀薄气体内部环流进行了模拟仿真分析,提出了以动量源模拟风扇段内流的仿真方法,并进行了可行性讨论。进一步实现了针对圆柱形模拟装置多工况下的内部气体流场稳态和非稳态计算仿真,并对计算结果进行了分析讨论,为火星大气环境模拟装置的设计提供了技术支持和参考。     

弹体结构平方米级扇形区域脉冲外压加载技术

针对弹体外表面平方米级区域的高峰值、短脉宽压力加载难题,依据一维可压缩流设计了高压腔、膜片腔、试验腔的气腔加载系统,附加了与试验件随形、扩展加压面积、均化压力分布的水腔,实现了平方米级区域1 MPa峰值、10 ms上升沿的脉冲压力加载要求,建立了弹性水腔的两自由度压力分析模型,获得了水腔系统的动特性参数,计算拟合曲线与...     

关于涡方法的讨论

本文讨论了如何提高二维涡方法的计算精度和减小涡方法计算量的问题。文中指出了改进的途径,同时列举了当前这些方面的研究进展情况。     

空间飞行体与压缩区内等离子体非稳态相互作用的控制方程

本文详细研究了压缩区内飞行体与等离子体之间的非稳态相互作用，得到了这种相互作用的控制方程。作为辐射源的飞行体天线系统，可以激发出可探测到的电磁弧波。带电粒子密度扰动对电磁弧波的影响表现为一个振幅的衰减因子。这表明，可通过对这种电磁弧波的探测而发现隐身飞行体的踪迹。     

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

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

基于增广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,无量纲化拟合速度平均偏差明显小于速度平均不确定度。实验结果证明了无量纲化对数校准数学模型应用于可压缩流体热线探针校准的可行性。     

高滑翔伞伞绳叉联结构优化设计

为获得翼伞伞绳最优叉联结构设计方案,文章提出了 1种高滑翔翼伞伞绳叉联结构的优化设计方法。首先,通过计算流体力学方法进行翼伞气动分析,确定最优气动性能时对应的伞绳安装角及翼伞伞衣气动载荷;然后,通过几何分析和受力分析,建立伞绳结点处的力系平衡方程及伞绳长度几何控制方程;最后,在满足伞绳结构安全的前提下,以伞绳总长度最小为优化目标,采用基于外罚函数和 Davidon-Fletcher-Powell(DFP)算法原理的优化算法获得不同约束力条件对应的翼伞伞绳叉联点位置方案。优化结果表明:与优化前叉联方案相比,该方法伞绳长度可减小 9.0%,阻力系数降低 8.8%。此外,通过增加伞绳受力约束,可实现翼伞滑翔性能的进一步提高。