关于超声速剪切流动的数值模拟

本文通过采用混合型的四阶精度频谱关系保持格式和基于信号特征传播的边界处理方法,对对流马赫数分别为0.5、0.9和1.4的时间发展和超声速空间发展的剪切流动进行了数值模拟,目的是对剪切层不稳定结构的演化方式、相同对流马赫数下时间发展和超声速空间发展的剪切流动的异同作初步研究.研究表明,对时间发展的剪切流动来说(1)当Mc=1.4时,剪切层内流动演化的形态为首先压力等值线在其极值点处产生中心型结构,且极大值点和极小值点交替分布.然后中心型结构进一步发展并产生横向分裂,从而在剪切层内等压力线出现鞍点结构.随后剪切层内重新产生中心型的压力极值点结构,进而又出现鞍点结构.如此发展下去,剪切层内的压力等值线呈现中心和鞍点的组合.剪切层外的压缩波沿特征线方向传播并加强,形成类激波的结构,最后在流场中包含许多压缩和膨胀波.在计算时间内流向的周期性未受破坏.(2)当Mc=0.5时,初始的流向的周期结构经过一定时间后将丧失其稳定性.在计算域为两个扰动周期的情况下,扰动形成的旋涡的后期演化表现为涡的对并.(3)对Mc=0.9的流动,演化图象介于前两者之间.对具有相同对流马赫数的空间问题来说,由于粘性和非线性发展对周期性的破坏,随着空间距离的增长,流动演化与时同问题的差异增大.相比较而言,高对流马赫数的时间问题与空间问题的差异明显大于低Mc下的差异.     

超声速高度欠膨胀冲击射流的大涡模拟

冲击射流广泛应用于短距起飞垂直降落飞行器(SVTOL)等航空航天领域.本文采用大涡模拟方法对高度欠膨胀的超声速冲击射流的流场进行了数值模拟.本文数值模拟得到了高度欠膨胀冲击射流流场中的激波结构和内外剪切层中不同尺度的涡结构.数值结果观察到了马赫盘的振荡,以及在斜激波、马赫盘及大尺度涡结构的共同作用下,射流内外剪切层之间的环形激波的生成与消失的周期过程.并对流场内剪切层的涡结构的演化进行了研究,数值结果显示内剪切层的大尺度涡结构的形成与马赫盘的振荡相关,在内外剪切层的作用下形成了壁射流区内外交错的涡结构.     

颗粒两相流自由剪切层大涡模拟

采用大涡模拟(LES)方法数值仿真了可压缩流自由剪切层内颗粒两相流流场结构.气相控制方程采用二维盒式滤波法滤除小于网格尺度的小涡以得到LES控制方程,而流场中小涡对大涡的影响采用K-ωSST湍流模型模拟.在Lagrange坐标系下用颗粒轨道模型仿真颗粒运动,用双线性插值方法实现Euler坐标系到肌Lagrange坐标系的转换,并考虑了颗粒速度和温度的改变对气相流场的反作用.分析了仿真结果中气相自由剪切层的涡结构特性,不同St数颗粒在流场中的运动特性,颗粒对流场的反作用效果,得到了合理的结论.     

边界层厚度对腔体气动声学特性影响数值模拟

为了研究来流边界层厚度对开式腔体气动声学特性的影响,基于分离涡模拟方法,计算了来流马赫数为2.0条件下,不同来流边界层厚度与腔体深度比时,长深比为5.88的腔体流动特性,得到了该腔体声压级的频谱特性.计算结果表明:随着来流边界层厚度增加,形成的剪切层稳定性增强,失稳后上下摆动幅度减少,失稳生成的大尺度涡与超声速主流的相互作用减弱,使得大尺度涡发展到腔体后缘时所具有的平动动能和转动动能降低.大尺度涡撞击腔体后缘在腔体内形成的气动噪声的声压级降低,最大减小幅度达7.5dB.同时各阶模态的频率也发生偏移,偏移值在100Hz左右.基于新的假设重新推导了Rossiter公式,明确了经验常数的物理意义,并以此解释了频率偏移现象.      

高速可压缩来流与凹腔驻涡相互作用及其物质输运特性研究

为了深入认知凹腔驻涡与高速可压缩来流之间的相互作用过程以及来流马赫数对上述动力学过程的内在影响规律,对长深比为1.0的单凹腔基础构型在常温、常压、Ma=0.3~0.7的高速可压缩来流下进行冷态非定常数值模拟,并基于拉格朗日视角研究来流与凹腔驻涡相互作用及其物质输运特性。结果表明：凹腔上方剪切层内旋涡结构仍然具有类似自由剪切层的周期性演化过程,主流与凹腔通过三次掺混进行物质输运;随着主流马赫数的增加,旋涡运动方式及凹腔-主流物质掺混机制都发生了本质上的改变,剪切层内旋涡和凹腔主涡涡量均呈现数量级式增长,剪切层旋涡结构脱落频率近似线性增长并逐渐卷吸至凹腔中部,进入凹腔的主流流体占比降低,但扩散程度和速度加快;此外,源于凹腔的流体流出唇口线后经剪切层的夹带作用回流占比增大,降低了凹腔-主流的整体质量交换水平;驻留参数整体衰减,高马赫数时仅为低速工况的1/2,显著降低了回流区对于流体驻留作用的贡献。     

基于间断有限元方法的可压缩混合层中大尺度结构和压缩性效应的数值模拟研究

基于间断有限元方法(DGM),对二维和三维、对流马赫数Mc=0.4和0.8的空间发展的可压缩混合层进行了数值模拟,并对混合层中的大尺度结构及压缩性效应进行了研究.重点观察了大尺度结构随压缩性变化的全景图像,文中采用λ2方法定义三维流动中的大尺度结构,流场中能够观察到大尺度的涡结构森林.当前数值模拟表明混合层的湍动能和雷诺应力随对流马赫数的增加而降低,这对RANS和LES模型的建立具有一定的参考意义.本研究也表明,间断有限元方法能够成功地应用于可压缩湍流的数值模拟.     

时间发展的二维混合层的物理分析和数值模拟

本文由四个部份组成.第一部份研究时间发展的二维剪切层流动的拓朴结构,提出了以下拓朴规律(1)沿零u线,流线的结点和鞍点是交替分布的;流线的鞍点和结点分别近似对应于压力值最大和最小的位置点;(2)压力值最大和最小的位置点是等压力线的中心点或鞍点;(3)如果在零u线上某点"0”,出现     

用于可压缩自由剪切流动的湍流混合长度

抓住可压缩流动变密度特性,构造出基于有效涡量的三维von Karman混合长度。湍流模型采用仅依赖湍动能k的单方程KDO(Kinetic Dependent Only)模型,引入新构造的混合长度替换旧尺度得到CKDO模型。为了验证其描述可压缩自由剪切湍流的能力,选择无壁面束缚、密度梯度大和可压缩效应强的自由剪切混合层为算例,其对流马赫数Ma_c=0.8。计算结果表明,KDO模型对混合层的速度分布有着良好的控制和模拟,而经可压缩修正后的CKDO模型与原模型及其他可压缩修正模型相比,所计算的速度分布、主雷诺剪切力和混合层厚度与试验结果更加接近,说明了该混合长度对可压缩混合层这种自由剪切湍流有着良好的刻画能力。     

三维超声速混合层内标量混合的大涡模拟

对三维对流马赫数0.62的超声速混合层流动的标量混合进行大涡模拟, 控制方程对流项采用五阶精度的WENO格式求解, 小尺度涡的作用采用一方程LDKM亚格子模型处理, 过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了混合层流场大尺度拟序结构以及标量场的演化过程, 研究表明标量混合过程受混合层内涡系演化所控制, 标量场具有明显的三维特征.模拟得到的速度、组分及其脉动的统计时均结果和实验结果相符较好.      

圆柱跨声速绕流中的激波/湍流相互作用大涡模拟研究

采用大涡模拟方法研究了圆柱跨声速绕流中的激波/湍流相互作用问题,来流马赫数M∞取为0.75,基于圆柱直径D的雷诺数为2×105。计算结果表明,圆柱分离点处出现一道斜激波,并且以与涡脱泻Strouhal数一致的特征频率向上游传播。激波运动导致流场中出现反对称的流动模态,剪切层中压力信号的功率谱曲线中存在0.4、-1和-5次方的斜率关系,剪切层中的剪切应力角约为0°,脉动速度以流向脉动速度为主,并且沿剪切层的大尺度结构组织性减小。     

弯曲混合管引射系统气动特性实验研究

对波瓣喷管和圆喷管与弯曲混合管匹配的引射系统内流气动特性进行了实验研究, 结果表明:(1)与圆喷管相比, 波瓣喷管可以将系统的引射能力提高约100%, 混合管出口的峰值温度降低近10%, 而两者的总压损失几乎是一致的;(2)因弯曲引起的二次流及主流冲击弯曲壁面的扩散影响, 混合管椭圆形排气口的长轴方向, 尾缘温度明显地高于前缘温度;而在短轴方向, 温度分布呈马鞍形状。      

PEEK卡箍载荷及耐环境性能研究

为支撑型号设计选用,针对P型和马鞍型两种结构的聚醚醚酮卡箍进行了拉伸载荷试验和多种环境试验,分析了振动与冲击、温度老化、液体浸渍、湿度、盐雾等环境试验条件对产品性能的影响,试验得到了卡箍最易失效的受力方向以及不同安装方式、不同温度条件下的拉伸载荷数据。研究表明:PEEK卡箍在试验环境温度下拉伸性能比较稳定;在湿热环境下易出现开胶现象;在振动与冲击、温度老化等环境载荷下能保证结构完整和性能稳定,可以在一定条件下代替金属箍带卡箍。     

An online data driven actor-critic-disturbance guidance law for missile-target interception with input constraints

In this article, we develop an online robust actor-critic-disturbance guidance law for a missile-target interception system with limited normal acceleration capability. Firstly, the missiletarget engagement is formulated as a zero-sum pursuit-evasion game problem. The key is to seek the saddle point solution of the Hamilton Jacobi Isaacs(HJI) equation, which is generally intractable due to the nonlinearity of the problem. Then, based on the universal approximation capability of Neural Networks(N...     

定常三维分离准则探讨

本文从奇点理论拓扑规律出发,说明在垂直于分离线及物面的截面上的截面流线有半鞍点存在,进而说明Lighthill等人认为的极限流线要在分离线以前抬起来是不合适的。合理的说法应该是极限流线总是同剪切应力线一致。只有在进入奇点以后,极限流线才抬起来离开物面。本文最后表明,单纯把二维分离准则用到三维分离的情况是不够的,并给出三维分离准则的表达式。     

混合初边值问题解的渐近性质及其差分解

 本文讨论了当t→∞时拟线性方程组初边值解的一些渐近性质。在假设非定常问题的解唯一、定常问题的解存在时,对同一边界条件,证明了当t→∞时非定常方程组的解趋于定常问题的解,讨论了用差分方法求解这种问题时,边界点上一些量的精确算法;为了消除在物体尖角附近出现的局部非物理解,得到正确解,我们还给出了一个简单的耗散型物面边界格式;最后给出了一些尖角凹形体的算例。     

带收益的车辆路径问题研究综述

车辆路径问题是企业实现物流配送的关键环节,对带收益的车辆路径问题的研究进行了综述。根据目前该问题的研究进展,对相关的研究进行了分类;分析了该类问题的特点,探讨了相应的0-1整数规划模型及集划分模型,总结了求解该问题的精确算法与启发式算法。介绍了该问题在实际中的应用,展望了其研究前景,为相关研究指出了方向。     

旋涡运动研究的某些进展

 本文简要介绍研究旋涡运动在以下问题上的某些结果:低速不同后掠角三角翼在各个迎角下的九种分离流类型及其边界;应用微分方程定性论与拓扑学对三维分离流与旋涡流的分析;旋涡破裂形态,对三角翼前缘涡破裂的实验研究与理论分析;受控分离与旋涡的干扰,二旋涡的位移、绕转与合并等。     

叶片型面刀位计算方法分析

为实现汽轮机７００ｍｍ扭叶片型面的数控砂带磨削精加工,通过分析研究有关现有文献的刀位计算方法,提出采用鞍形曲面法,对由离散型值点给出的汽轮机扭叶片,可以确定其宽行线接触加工的刀位轨迹,以实现叶片型面的数控砂带磨削精加工。     

磁悬浮控制系统混沌预测

 根据谐波平衡原理,理论上对非线性磁悬浮控制系统的倍周期分叉及混沌行为进行了分析,给出了系统产生倍周期分叉的近似条件,确定了系统混沌存在近似的位置,这为保证此系统实现稳定控制及优化控制器设计提出了一种有利的理论依据。     

考虑使用因素的涡扇发动机排气温度换算方法

针对某型涡扇发动机厂内排气温度换算值验收合格而外场地面检查实测值偶有不合格的问题,分析了只考虑大气温度单一因素的排气温度换算方法的缺陷。提出了一种综合考虑非标准大气和调节规律使用因素的换算方法,建立了数学模型,计算得到了非标准大气的修正系数和调节规律的修正系数,并经试车试验验证。结果表明:提出的考虑使用因素的换算方法符合发动机的实际使用条件,所获得的修正系数与试车试验数据的相对误差小于1.3%,有效解决了排气温度厂内验收合格而外场地面开车不合格的问题。