基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

基于线性回归近似替代模型的分布式卫星系统不确定性分析方法

分布式卫星系统在空间科学领域具有广泛的应用前景，同时也对概念设计阶段的分析工作提出了更高的要求.在分布式卫星系统概念设计阶段，评价不确定性参数对于最终探测效能的影响有着重要的工程应用价值.传统的不确定性分析方法存在解析困难、数值模拟计算效率低且耦合关系表征不明确的缺点.本文结合概念设计阶段不确定性分析的特点和需求，提出一种线性回归近似替代模型，将不确定性参数对于探测效能影响的结果计算降维为一个线性组合系数的求解问题.利用一个空间多点探测物理场分界面任务作为典型事例进行仿真验证，结果表明，相比传统数值模拟方法，本文方法在计算效率上具有明显优势且对关键不确定性参数有一定识别作用.      

风雷软件LES开发设计与验证

基于国家数值风洞风雷软件开源框架，设计开发了LES(Large Eddy Simulation)湍流模型，主要包括Fourier谱/有限差分方法 LES求解器和有限体积/有限差分方法 LES求解器。简要介绍了采用的不可压缩流动求解的投影法、Fourier谱/有限差分混合方法、亚格子模型等理论方法，给出设计的软件框架和计算流程，尤其说明采用的松/紧2种耦合模式。通过数值模拟不可压缩槽道湍流、亚临界雷诺数圆柱绕流、NACA0012临界攻角的低频振荡算例，验证求解器的计算精度和复杂湍流模拟能力。基于风雷开源框架设计的开源LES模型，具备高精度数值格式、亚格子模型、湍流统计等通用模块，可为国内学者提供一个LES湍流模拟研究的开放平台。     

气冷涡轮导叶流热耦合计算及机理

针对气冷涡轮叶片的多场耦合特性，利用流热耦合（CHT）方法，对采用不同气冷结构的高压涡轮导叶进行数值模拟。在内冷涡轮导叶算例中，对比实验数据选取精度较高的流热耦合计算方案，分析该内冷涡轮导叶的多场特性及耦合机理。在此基础上，以带有气膜冷却孔及内冷通道的气冷涡轮导叶为研究对象，重点围绕冷却射流与主流的相互作用，讨论近壁边界层中流热耦合关系及气冷效率影响因素等相关问题。结果表明:采用流热耦合计算方法及合适的湍流转捩模型有利于提高数值精度；气冷涡轮导叶的流场温度场密切耦合，流动换热特性互相影响；冷气射速低时，增加冷气流量可提高气膜冷却效率，冷气量达到一定值时，冷气流量增加将导致气膜冷却孔后上游冷却效果变差，下游冷却效果变好；冷气射速较高时，将与主流相互作用产生复杂流动结构（如肾形涡、马蹄涡等），对温度分布存在一定影响。      

左行运动激波主导的管内流动特性研究

为了分析左行运动激波主导的管内流动特征，本文采用非定常数值仿真方法，对亚声速进口条件下等直管道内左行运动激波传播与演化特性、左行运动激波/边界层干扰特征开展研究。研究结果表明：在出口周期性强压力脉动干扰下管内存在连续的左行运动激波，该左行运动激波传播特征具有相似性，激波强度、传播速度按幂函数规律衰减。气流经过左行运动激波后总压、总温、静压阶跃式升高，随后受膨胀波影响气流总压、总温、静压下降；左行运动激波/边界层干扰诱发形成翼型回流区，该回流区随运动激波强度衰减逐渐减小。理论与数值分析表明存在左行运动激波后速度为零和运动激波两侧总压相等的两个临界状态。波前马赫数低于临界值或左行运动激波强度高于临界值时，左行运动激波后为倒流、波后总压高于波前。     

SABRE预冷器计算模型及其在整机模型中的应用

预冷器的性能对SABRE等预冷组合循环发动机具有重要影响,为实现发动机方案设计阶段预冷器的快速设计与评估,建立了预冷器准二维快速评估模型。将SABRE预冷器的几何结构简化为一个扇环形区域,沿径向和周向将该区域划分为二维节点。应用守恒方程及传热关联式完成单个节点计算,再求解节点矩阵的平衡方程组,计算内外流体特定节点上参数的二维分布,得到预冷器出口参数。将模型嵌入发动机总体性能程序中,实现了发动机设计及非设计点的预冷器性能计算的功能。与文献数据对比结果表明,预冷器模型传热计算误差小于5%,摩擦阻力误差小于10%。整机计算结果显示,Ma0~5范围内,预冷器空气侧温降范围为143K至932K,温降随飞行马赫数升高单调上升。预冷器传热有效度范围为0.896-0.945,空气侧总压恢复系数范围为0.852-0.904。     

双射流环量控制翼型的控制力矩特性研究

为了探究环量控制技术在飞行控制性能方面的优势，在定常流场中对定常射流环量控制翼型的控制力矩作用机理展开了研究，采用数值仿真的方法，对比分析了单射流、双射流产生的虚拟舵面与传统舵面作用下的气动力系数的变化规律，并基于无舵面飞行器CCSCAOON对其气动力矩的控制特性进行了验证。验证结果表明:单射流作用下的虚拟舵面能够提供用于飞行器所需的滚转和俯仰力矩，且作用机理相似，控制性能优于传统舵面；无论是单射流还是双射流，在大迎角下虚拟舵面的气动控制特性较差，限制了环量控制的使用迎角；双射流较单射流而言，升阻比特性和控制力矩特性较好；双射流下的虚拟舵面通过调节下射流口动量系数，能够有效降低偏航力矩与滚转、俯仰力矩之间的耦合效应。      

冷压过程对球堆积结构渗透率的影响

空间轴承的供油系统主要目标是实现对飞轮轴承的长时间、稳定、微量供油，利用微孔介质的复杂结构可能有效实现这一目标，但目前的渗流机理和规律尚不明晰，主要问题在于多孔介质结构难以调控、多孔结构与渗透率之间的关系复杂、渗透率受装配压力影响等.芯阀材料一般经过两步制成：先对颗粒物质进行冷压，再对冷压得到的结构进行烧结得到成品.目前的生产过程主要以经验为主，机理不清，难以精准调控.本文在数值上结合离散单元法与格子玻尔兹曼方法对颗粒堆积结构在冷压过程中的渗透率变化进行了研究.采用离散单元法模拟颗粒堆积结构冷压过程，再用格子玻尔兹曼方法进行渗流模拟，得到不同摩擦角下得到的冷压结构的渗透率变化.本研究为全面揭示多孔介质冷压结构与渗透率的关系及其微观机理奠定了基础.     

四旋翼倾转飞行器操纵冗余设计

建立了四旋翼倾转飞行器旋翼、机翼、机身的非线性气动模型和飞行动力学模型,并在Matlab/Simulink环境下搭建了仿真模型,根据仿真模型计算各个操纵面的操纵效率,确定了不同状态下各个通道需要用到的操纵方式,结合飞行动力学模型和操纵方式,在纵向通道上进行了全包线的配平计算,表明了操纵策略的可行性。