稀薄流过渡区气固两相喷流的建模与数值模拟

从描述粒子运动的细观层次出发,运用DSMC方法,构造了稀薄流过渡区气固相互作用模型和颗粒碰撞模型,数值模拟了有气固两相横向喷流干扰的二维有限平板高超声速流场,分析了流场结构和物面气动系数分布等特性.结果表明,对于不同的来流Knudsen数和颗粒质量载荷等参数,喷流干扰流场存在差异,同时反映了基于DSMC方法的气固相互作用模型是研究稀薄及过渡条件下的气固两相流动的有效方法.     

气液交叉射流动力学行为的基础研究

为了模拟燃油在高速喷嘴内的快速蒸发混合过程,基于Eulerian—Eulerian方法,用RNGk一8湍流模型建立了高速喷嘴内两相流流动的数学模型,数值模拟了高速喷嘴内气液两相交叉射流的流动规律。改变气液动量比和韦伯数等条件下的两相流场计算结果表明：动量比增大3倍时横向穿透深度略有增大;而较大的韦伯数能起到强化横向射流的破碎效果。计算结果与同条件下的试验结果相符合。结果为研究HiTAC快速速喷嘴内可燃混合气的快速形成奠定了基础。     

乙醇-空气稀液雾值班火焰的大涡模拟研究

为加深对气液两相湍流燃烧现象的认知,检测火焰面模型在液雾燃烧中的适用性,本文在欧拉-拉格朗日架构下使用火焰面/反应进度变量模型(FPV)数值模拟了乙醇-空气稀液雾值班火焰。欧拉坐标系下的气相湍流场使用大涡模拟方法模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑了相间质量、动量、能量交换。燃烧模型中采用碳元素定义混合物分数,在混合物分数方程源项中体现液相对燃烧模型的影响。模拟得到的气相温度分布和液相统计值均和实验数据较好吻合,验证了该燃烧模型对稀液雾扩散类型火焰的适用性。分析瞬时图发现,该稀液雾火焰的最高燃烧温度往往位于当量混合物分数附近,在出口下游20倍直径处火焰完全点着,此处上游FPV模型能给出局部点火熄火现象。蒸发作用在剪切层和点火区域较强,而液滴和火焰作用较弱,单一液滴很少被火焰包围。     

两电极等离子体合成射流性能及出口构型影响仿真研究

通过将火花放电的物理效应等效为气体焦耳加热的过程,在能量方程中引入能量源项,进行了单次能量沉积下两电极等离子体合成射流的唯象模拟。提出将等离子体合成射流对外界流场的动能和热能注入分别作为表征射流"冲击效应"和"热效应"的参数。研究表明在单次放电条件下射流建立的自维持振荡过程中,射流动能和热能主要集中于主射流阶段且射流的"冲击效应"相比"热效应"衰减更快,在一个大气压下两电极激励器总的能量转化效率约为2.3%。分析了出口构型对射流的影响,研究表明收缩孔结构可以有效提高射流速度,但将导致射流动量和饱和频率的降低。     

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

矩方法及其在水凝结研究中的应用与发展

矩方法(Method of Moments,MOM),又可以称作内部坐标矩方法,通过求解液滴尺寸分布函数对内部坐标的各阶矩的输运方程,反映粒子数量及尺寸分布的演化。在各种有关伴随凝结流动的数值方法中,矩方法兼顾了计算效率和液滴分布信息,因此被广泛应用于水凝结的研究。本文首先介绍了经典矩方法,综述了矩方法的最新发展。然后应用经典凝结矩方法,研究了伴随同质凝结的Prandtl-Meyer膨胀流动问题,计算发现随着来流水蒸气饱和度的提高,凝结诱导的激波由定常状态向非定常状态转变,出现自激振荡现象;进一步,将经典凝结矩方法推广到异质凝结领域。通过引入"瞬间活化"假设,简化了异质成核模型,将此方法应用到伴随异质凝结的激波管问题中,进行了参数研究,加深了对异质凝结和流动相互作用的认识;最后针对经典矩方法无法处理相间相互作用的问题,将经典矩方法中的矩方程改写为基于液相平均速度的形式,发展了一套考虑相间质量、动量、能量交换的矩方法,应用此方法研究了伴随同质凝结的旋涡运动问题,分析了旋涡中凝结的产生过程,以及凝结与流动相互作用对液相分布的影响,计算结果显示了凝结生成的液滴在离心力作用下离开旋涡中心的现象。     

Maxwell气固相互作用模型对稀薄高超声速凹腔绕流流动特征和热环境的影响

针对高空航天飞机等再入飞行器表面缺陷或防热瓦缝隙导致的局部压力过高和气动加热问题,采用直接模拟Monte Carlo （DSMC）方法研究了飞行高度为80 km的稀薄流区高超声速凹腔绕流问题,考虑气固相互作用（GSI）模型对凹腔流场特征和表面压力、热流的影响。结果表明：稀薄流条件（80 km）下,GSI为完全漫反射时,在凹腔前缘分离的剪切层再次附着在后缘,在凹腔内部形成一个充满腔体的单涡结构;随着GSI从完全漫反射向镜面反射变化,气体与凹腔表面之间的切向动量交换减弱,即黏性剪切作用减弱,外部气流被卷入凹腔的程度减弱,导致涡结构不断减小直至消失,凹腔底部逐渐出现所谓的"死水区"。与完全漫反射相比,镜面反射或近镜面反射会导致凹腔上游侧面的峰值压力和峰值热流以及下游侧面的峰值压力剧烈增大,在飞行器设计中,应特别留意上述表面的压力载荷和热载荷。     

进口边界层对采用弯叶片的平面扩压静叶栅流场性能的影响

给定不同型式的进口边界层,在两种不同亚音速条件下对一平面扩压静叶栅的弯叶片流场进行了数值模拟。结果表明弯叶片对扩压叶栅的改善的能力受进口边界层的特征影响。这种影响分为两个方面:(1)边界层厚度的影响和(2)边界层动量损失厚度的影响。边界层越厚或动量损失厚度越大,在低马赫数条件下弯叶片对吸力面角区密流增加越明显,从而更大程度地提高了端区的流动性能,降低了叶栅损失。在高马赫数条件下,若边界层越厚或动量损失厚度越大,角区密流虽变化不大,但因端区损失较大,其性能的提高会给叶栅总性能的改善带来较大的收益。      

舰船燃气轮机高压涡轮颗粒沉积特性研究

为了探究典型舰船燃气轮机高压涡轮叶栅通道内部的颗粒沉积分布规律,论文以分析随燃气进入涡轮叶栅通道内的颗粒的动力学及颗粒与壁面相互作用的特性为出发点,探究颗粒与涡轮叶片表面撞击后发生黏附与剥离以及沉积与反弹的相互作用准则,在此基础上构建相应的颗粒壁面沉积模型,采用数值模拟方法,建立气-固两相流无量纲方程组,嵌入UDF用户自定义函数,并比较分析临界速度模型和临界黏度模型对于模拟颗粒沉积分布的异同。结果表明,在本文工况条件下,对于临界速度模型而言,叶片表面和下端壁处的颗粒沉积效率都在动量Stokes数增大到约0.1时开始从100%减小,在动量Stokes数增加到约10后减小到0;上端壁处的颗粒沉积效率在动量Stokes数增加到约1时开始从100%减小,在动量Stokes数增大到约100后减小到0。而对于临界黏度模型而言,叶片表面和下端壁处的颗粒沉积效率随着动量Stokes数的增大而不断增大,二者都在动量Stokes数约为1时达到100%,而上端壁处的沉积效率却先增大后减小,最后再增大,在动量Stokes数约为0.1颗粒沉积效率降低到最小值约为80%,后开始增加到100%。     

有蒸发稀液雾冷态射流的大涡模拟研究

为研究气相湍流场中液滴弥散和蒸发过程的特点和影响因素,加深对气液耦合作用的认知,用欧拉-拉格朗日相结合的办法数值模拟了稀液雾冷态受限射流和自由射流。欧拉坐标系下的气相湍流场使用大涡模拟技术模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑相间质量、动量、能量交换,采用相平衡蒸发模型模拟液相的蒸发效应。模拟得到的液相速度场、粒径分布和质量流率等均和实验数据较好吻合。有回流的受限射流中初始惯性对Stokes数大的液滴弥散影响较大,而Stokes数在1附近的小液滴弥散更受气相流场结构影响。液雾单位体积蒸发速率和液滴弥散分布密切相关。自由射流中液滴受初始惯性和受气相流场作用大体一致,大小液滴的弥散分布基本满足自相似性。     

谱方法的基本问题与有限谱法

本文概述了谱方法的特点，讨论了谱方法的四个基本问题及解决方法，由此得出谱方法局域化的重要性，重点介绍了作者提出的有限谱法两个发展阶段，特别是最近的基于截断傅立叶积分的有限谱法，列举了有限谱法的一系列计算格式，给出了计算高超声速流，各向同性湍流，地震波的传播，分区算法，复杂形状物体，激波管问题，KdV方程，湍流和激波相干扰等问题的计算结果。     

连续体结构拓扑优化方法评述

连续体结构拓扑优化在优化中能产生新的构型,对实现自动化智能结构设计具有重要意义。目前,连续体结构拓扑优化方法主要有：均匀化方法、变厚度法、变密度法、渐进结构优化方法、水平集法、独立连续映射方法。本文首先系统回顾了以上方法的发展历程,介绍了它们的研究现状。其次,通过对比以上拓扑优化方法对若干典型算例的优化结果,表明以上方法都有较好的减重效果。最后,对以上方法进行了总结,列出了它们的优缺点和发展方向。     

一种改进型洞映射法

针对传统洞映射法存在挖洞曲面必须封闭的局限性和消耗大量内存的缺点，提出了改进型洞映射法。主要做了两点改进：1）采用“反向标记法”判别洞映射笛卡尔网格单元属性，使得挖洞曲面不再局限于封闭的情况；2）运用“矢量射线法”进一步判断洞映射笛卡尔网格边缘单元内的计算网格点属性，使得笛卡尔网格分辨率低时仍能精确挖洞。研究表明：相比于传统洞映射法，改进型洞映射法保持了很高的计算效率和精准挖洞的能力，还极大地降低了内存消耗，并且具有较强的可靠性。     

柔性决策

提出了柔性决策的思想,并对决策科学的发展进行了回顾。为解决柔性决策问题,提出了极大－极小方法、人机交互方法和集包含方法。     

含裂纹铆接加筋板应力强度因子的计算

本文提出了选择复连通域边界(多环路)计算J积分的一种新方法,成功地解决了单连通域内有铆接点时计算J积分的问题。文中还提出了一种“虚设铆钉法”,采用这种方法,可利用同一网格有限元模型计算不同裂纹长度下裂纹尖端的应力强度因子,这大大简化了模型准备工作并减少了计算结构总刚度矩阵的工作量。本文采用刚度导数法、J积分法及能量释放率法,计算了含裂纹铆接加筋板的应力强度因子,通过计算结果对比,对各种方法作了粗浅的评价。     

基于水平集方法的连续体结构拓扑优化方法研究

考虑到现有各种连续体结构拓扑优化方法的优缺点,提出了一种新的连续体结构拓扑优化新策略。该方法将无网格法和水平集方法结合起来,充分发挥各自的优点,使得该拓扑优化新策略具有良好的收敛性和计算精度。     

典型翼型参数化方法的翼型外形控制能力评估

现代飞行器的精细化设计需求要求翼型参数化方法具有精准的翼型外形控制能力。依据标准的风洞模型公差要求,以五个典型翼型作为测试翼型,分别对Hicks-Henne参数化方法、样条参数化方法、PARSEC参数化方法和CST翼型参数化方法的翼型外形控制能力进行评估,依据每个翼型参数化方法的评估结果并结合其固有特点进行适应性分析。结果表明:四种翼型参数化方法对翼型外形的控制能力存在较大差别,具有不同的适用范围。评估典型翼型参数化方法的所得结果可指导设计人员依据具体的设计需求选择适当的翼型参数化方法。     

基于混合 MOM-PO的复合电磁散射问题分析

海面(或地面)环境中目标电磁散射特性,对实际战场环境下军用目标的雷达探测、雷达特性形成以及雷达目标识别等技术具有基础支撑作用。因此,快速准确地分析目标和粗糙面复合电磁散射特性具有重要理论意义和应用价值。文章采用基于矩量法(MOM)和物理光学法(PO)的混合方法MOM-PO分析了目标和粗糙面复合电磁散射问题。与传统的矩量法相比,该方法在保证精确度的前提下大大加快了计算速度。     

关于电磁场的自然单元法

有限元法和无网格法在电磁场数值计算中已经得到了广泛的应用,然而有限元法存在前处理网格剖分问题,无网格法存在计算时间长、边界条件和不连续面难处理等问题。针对以上问题,本文提出了利用自然单元法求解电磁场的方法,该方法根据场域中的离散点的信息构造自然邻点插值函数,最后算出场值。通过算例证明了该方法在电磁场计算中的可行性。