用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

基于组合优化策略的月球软着陆最优轨道设计

基于Pontryagin极大值原理,把求解月球软着陆燃料最优化问题归结为终端自由型两点边值问题.采用粒子群算法和单纯形算法接力优化的组合优化策略,在初始猜测值的邻域内进行搜索,充分利用粒子群算法的全局搜索能力迅速缩小搜索范围,然后利用单纯形算法的局部搜索优势快速获得优化结果.该优化策略最大的优势是使粒子群算法的全局搜索能力和单纯形算法的局部搜索能力同时得到最大化的发挥.仿真证明该优化方法在考虑一些实际工程约束的情况下,能较快速而准确的获得月球软着陆优化轨迹,具有一定的优越性.     

低速风洞均匀吸气地板下阻塞对地板边界层影响的研究

为了研究低速风洞地板下部阻塞对地板上表面边界层的影响以及控制方法,在国防科技大学KD-03低速风洞利用均匀吸气地板系统进行了实验研究.地板下部阻塞对地板上表面主流区的流动产生干扰,并影响边界层的分布和发展,随着地板下部阻塞度的增加,地板上表面边界层的厚度有增加的趋势;在某一阻塞度和吸气系数下,来流速度越大,地板下部的阻塞对主流区流动的影响越小;地板的均匀吸气使边界层的厚度显著降低,也有效减小因地板下部阻塞引起的主流区流动的不均匀性.该实验的研究结论为8m×6m风洞均匀吸气地板系统研制提供了参考.     

弹性支撑条件下分段轴压阶梯梁自由振动及稳定性分析

分析了在各种经典及其派生弹性支撑边界条件下任意多段分段常轴压阶梯梁的自由振动和稳定问题,得到了形式统一的封闭解析解。首先采用分布传递函数方法,得到了轴压等截面梁在各种边界条件下自由振动及稳定性分析的解析解,然后根据分段轴压阶梯梁横截面弯曲刚度和轴向压力沿轴线的变化情况,将梁分成多段子梁,再利用各子梁的解析解以及各子梁间的位移连续和力平衡条件,得到粱的各阶自由振动频率和失稳载荷及其相应的模态形状。最后通过三阶梯梁的三个算例验证了本文方法的正确性,并针对各种边界条件,计算了各种边界条件下分段轴压四阶梯梁自由振动的前三阶固有频率和失稳载荷参数,分析了固有频率和失稳载荷与其模态形状的关系,探讨了弹性支撑条件下轴压阶梯梁固有频率随弹簧常数的变化趋势。     

一种渐变无反射边界条件的建立与验证

建立并验证了一种适用于非线性流动的新无反射边界条件,以简单的迁移方程模型为基础,提出了一种通量渐变消失无反射边界条件。通过乘以合适的函数,在渐变层内偏微分方程中的空间偏导项逐渐消失,从而在边界点上得到只含有时间变量的常微分方程,来解决自由区域处边界造成的定解问题。将渐变边界条件推广到了一维和二维非线性Euler方程,其方程形式简单,便于实现。在计算域边界处采用空间滤波和二维情形时角区内引入牛顿阻尼区提高了渐变边界条件的数值稳定性。在非线性情况下对二维Euler方程渐变边界条件进行了一系列的数值验证,得到了令人满意的数值结果。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

不同波系配置的鼓包压缩面流动特性实验研究

为探究波系配置对鼓包压缩面流动特性的影响,通过纹影和基于纳米粒子示踪的平面激光散射系统(NPLS)观测以及压力测量等实验手段,研究了三种不同波系配置的鼓包压缩面诱导的复杂流动,获得了鼓包压缩面上的压力分布、边界层分布特性以及鼓包压缩面上的精细流场结构。结果表明:三种不同波系配置的鼓包压缩面均具有边界层排移能力,在鼓包对称面上,沿流向边界层不断降低,边界层不断向两侧排移,在鼓包侧面形成堆积。尽管如此,压缩过程的不同,边界层的排移过程存在一定差异,对于单锥鼓包,其最先将低能边界层排移掉,而双锥和等熵锥鼓包略为迟缓,但三种鼓包末端的边界层厚度大体一致。此外,三种鼓包压缩面对称面上前中部压力梯度存在较大差异,但鼓包末端具有相同的压力分布。因此,在保证总偏转角一致的情况下,具体选用哪种鼓包配波方式对于边界层排移、增压能力而言没有显著的差异。     

局部凹槽对高超声速边界层中第二模态扰动演化的影响

采用直接数值模拟方法研究了局部矩形凹槽对来流马赫数为6.0的高超声速平板边界层中第二模态扰动演化的影响。引入了透射系数(定义为凹槽下游与上游扰动渐近幅值的比率)来量化凹槽的局部散射效应。数值结果表明:对于较浅的凹槽,频率较低的第二模态扰动被促进,而高频扰动的规律相反。对于大多数情况,凹槽深度的增加使得透射系数减小,这表明凹槽对扰动的促进作用减弱或抑制作用增强。当凹槽的深度超过某一临界值时,透射系数与凹槽深度的关系变为正相关,这表明了另一机制的出现。临界深度受扰动的频率影响:频率越低,临界深度越大。      

微结构的拓扑表达模型及自由排布阵列加工方法

制造过程中需要微结构阵列的数学模型从而进行图纸设计、刀具补偿计算、切深数据生成等.本文从拓扑原理出发,结合方向关系与距离关系,得到微结构阵列空间关系的统一表达模型.搭建快刀伺服加工实验平台,利用统一表达模型,并以具体样例加工了微透镜阵列,实现了自由排布的微结构加工.     

模拟弹塑性碰撞的修正SPH方法

从基本的无网格光滑粒子法SPH(Smoothed Particle Hydrodynamics)近似出发,修正了模拟固体力学中大变形弹塑性碰撞的SPH方法.在边界处采用修正的边界条件,弹塑性分析过程中采用增量理论计算应力,迭代过程中用守恒光滑法进行滤波修正消除拉力不稳定.对SPH方法进行了程序实现,给出了杆弹塑性碰撞的算例.计算分析表明,SPH方法节点的影响域较大、精度较相同节点间距有限元法的结果有一定差距,但是通过增加粒子数量可以提高SPH的精度,保持了其简单性和计算大变形的特性.