高速旋转冲压增程弹用进气道复杂流场数值模拟

采用块结构网格与二阶精度流场分区求解技术,对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟,得到了高速旋转工况下对应于不同来流马赫数和攻角,以及临界工况时超音速进气道内外流场复杂的波系结构。在旋转工况下,进气道流场结构与不旋转时相似,但旋转速度以及攻角的存在对冲压发动机进气道的总体性能产生了负面影响,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大;冲压发动机进气道在较低马赫数工作时的综合性能优于在较高马赫数工作时的综合性能。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

二维三轴编织复合材料压缩失效行为的细观有限元模拟

为研究典型二维三轴编织复合材料（2DTBC）的压缩破坏机理,建立了细观有限元模拟方法体系。提出了反映编织复合材料真实几何特性的单胞模型建模策略,根据Murakami-Ohno损伤理论建立了各向异性损伤模型来模拟纤维束中的损伤起始和扩展行为,通过引入波动系数描述了纤维束的起伏状态,并采用内聚力单元来模拟界面分层。在此基础上,分析得到了二维三轴编织复合材料在压缩载荷下的破坏过程,研究了压缩载荷下纤维束和界面层的损伤演化,探讨了纤维束波动对压缩性能的影响规律。通过与相关试验结果对比,该模型能够准确预测二维三轴编织复合材料在面内压缩载荷下的力学响应和主要失效行为,以及自由边效应。细观失效过程分析结果表明,二维三轴编织复合材料轴向压缩的破坏是由轴向纤维束的纤维压缩失效主导的;横向压缩破坏则是由偏轴纤维束的纤维压缩失效引起的。     

二元曲面可调进气道流量系数精确预测方法

为了满足二元曲面可调进气道模态转换马赫数范围(来流马赫数为2.2~3.2)的流量要求,针对唇口平移、转动和转动+平移三种调节方案,基于理论分析和基准进气道的流场,提出了一种流量系数精确预测方法,并通过数值计算进行验证且获得了进气道的总体性能。结果表明:调节后的进气道流量系数与预测值完全相等,而且无需多次试算,符合设计预期,可拓展应用于轴对称进气道。相对基准进气道,唇口前移时流量系数和压缩效率同时增加,来流马赫数为2.5时出口总压恢复系数相等而增压比增加了14.6%;在降低相同流量系数条件下,后移唇口使得增压比和压缩效率均降低,来流马赫数为2.5时出口总压恢复系数基本相等而增压比减小了12.9%,转动唇口使增压比进一步减小了9.1%,唇口后移方案性能更优。      

螺旋桨滑流对发动机进气道气动性能的影响

针对螺旋桨/短舱/进气道一体化构型,采用多重参考系(MRF)法和滑移网格非定常法进行了数值模拟研究,实现了螺旋桨/短舱/进气道一体化内外流耦合流场仿真,研究了不同工况下螺旋桨滑流对进气道气动性能的影响,并与无桨构型的进气道性能进行了对比,最后通过数值分析解释了螺旋桨桨盘后总压上升的原因。结果表明:多重参考系法和滑移网格非定常法对进气道出口的时均参数的计算结果吻合较好,多重参考系法可用于螺旋桨及进气道的气动性能参数计算;螺旋桨做功使进气道总压恢复系数增高,螺旋桨转速越大,总压恢复系数的增量就越大,且螺旋桨滑流使进气道出口畸变指数明显增大;通过数值分析推导的桨盘前后总压关系公式与仿真结果的误差均在3%以内。      

进口水平投影可控的流线追踪内收缩进气道设计

为了满足两侧进气布局飞行器的乘波前体与进气道一体化设计要求,提出了一种进口水平投影可控的流线追踪内收缩进气道设计方法。基于马赫数分布可控的轴对称基准流场,在指定进口水平投影为椭圆的条件下,采用该方法设计了内收缩进气道并在设计点(Ma=5.4)和接力点(Ma=4.0)对其进行数值研究。结果表明,设计点时进气道都能保持基准流场的波系结构和沿程压力分布,无粘时可以全捕获自由来流,喉道性能与基准流场几乎相等。有粘条件下,设计点和接力点时进气道具有较高的压缩效率和良好的流量捕获能力,接力点的流量系数高达0.85。该设计方法为内收缩进气道与乘波前体的一体化设计提供了新途径。     

超额定状态下二元超声速进气道的流动特性

为了揭示超额定工作状态下超声速进气道内的复杂流动机理,对一设计马赫数为2.0的二元超声速进气道开展了数值模拟研究,获得了其在不同来流马赫数状态、不同节流状态下的流场结构.结果表明:当进气道工作在额定状态时,随着节流程度的增加,其激波串的核心区由偏向下壁面摆至偏向上壁面;而在超额定状态下,由于依次受到唇罩内侧分离包、唇罩激波等的影响,激波串核心区则由偏向上壁面转而摆至偏向下壁面.在来流马赫数为2.5的节流状态下,其唇罩激波与前体斜激波相交形成了马赫杆等复杂波系结构,而来流马赫数为3.0状态却并未形成此类现象.在上述两种超额定工作状态下,前体斜激波的上透射激波均在高反压条件下演化为正激波形态,而唇罩激波的下透射激波形态也发生了明显的改变.      

一种正向乘波前体设计方法初步研究

以定楔角乘波体设计方法为基础,研究了影响高超/超声速乘波体"乘波"的主要因素,给出了前体前缘实际气流压缩角的确定方法及影响因素,可知在相同的来流马赫数和压缩角δ下,随着前缘角θ和气流与前缘夹角α的增加,实际气流偏转角γ减小。据此,基于幂函数进气道前体构形,给出了前缘激波不脱体的限制条件及具体的判定方法,分析了乘波体典型几何特征参数对前缘激波不脱体的影响规律,结果显示在相同的来流马赫数和压缩角度下,增大前缘形状因子n,减小前体的长宽比L/W及增大前缘角均有利于激波不脱体。根据给出的前体几何参数对前缘激波脱体的影响规律曲线,对一种"前体几何外形构造+前缘激波附体条件限制"的正向前体乘波器工程设计方法进行了研究,给出了具体设计流程,并进行了初步的数值仿真验证,表明通过该方法设计的乘波前体流动特征与预期的结果吻合,说明文中所给出的激波附体条件及影响规律是可信的,乘波前体设计方法是可行的。     

进气道堵盖打开过程数值模拟

为获得整体式固冲发动机转级过程中进气道堵盖(包括入口堵盖与出口堵盖)打开过程的流动形态,建立了进气道二维模型,利用Fluent动网格技术和UDF方法,开展了进气道堵盖打开过程非稳态流场研究。结果表明,在入口堵盖打开前,进气道前端形成强烈的弓形激波;在入口打开、出口未开的过程中,沿进气道轴向各监测点压强呈现周期性变化,振荡频率为100 Hz左右,出口堵盖位置压强振荡幅值为0.53 MPa;在出口打开后,补燃室残余热量形成的压强峰导致进气道在短时间内无法起动,随着背压降低至小于进气道再起动反压,进气道完全起动。     

旋转光滑直通道湍流流动二维热线实验

采用二维热线测量了旋转光滑直通道内不同流向位置上的平均速度和雷诺应力。实验结果表明:较高的当地旋转参数使得旋转对平均速度的影响区域扩大并且导致了前缘面附近湍流流动的再层流化。后缘面附近无量纲主流平均速度型严格按照旋转数顺序依次排列,且在对数律区呈现对数律分布。与此同时,所有无量纲雷诺应力分量在后缘面附近基本不受旋转影响。再层流化导致了前缘面附近无量纲主流平均速度型无法在对数律区维持对数律分布,且所有雷诺应力分量都随着旋转数和流向位置半径增大不断衰减;经过u-v象限分析,再层流化现象的直接原因被归于湍流脉动生成减少。