湍流模型对喷流计算精度的影响研究(英文)

采用三维RANS代码MBNS3D对包含了内流、外流及后体干扰的三种简化尾喷流,进行了数值模拟研究。通过壁面压力分布及喷流截面速度型与实验数据的对比,研究了SA和SST两种湍流模型及其修正(SACC,SADES,SSTCC)对喷流计算精度的影响。本文计算表明,SSTCC模拟效果最为精确;利用SA模型进行脱体涡模拟,在同样的网格下优势不明显。湍流模型可压缩修正改善了模拟喷流的膨胀趋势,但对附面层分离位置的预测,在本文的工作中,还不如原始模型。两类模型相比,SST模型的结果好于SA模型。     

电脉冲除冰系统电磁脉冲力仿真分析

为了获得更精确的电脉冲除冰(EIDI)系统中电磁力在飞机蒙皮上的分布及其随时间变化的情况,针对系统的工作原理与工作过程进行了仿真研究。计算分析了EIDI系统电路中电流随时间变化的关系;通过麦克斯韦方程组的求解完成了EIDI系统电磁场的仿真;计算分析了电磁力随时间变化的规律,获得了蒙皮上的磁感应强度分布、涡电流密度分布和瞬态电磁力密度分布。考虑线圈工作时蒙皮对其电感的影响,建立了更接近真实电路的激励加载方法,得到了准确的电磁力变化趋势;考虑趋肤效应对蒙皮厚度方向上涡电流密度、电磁力密度分布的影响,用有限元分析方法使蒙皮上电磁力分布的计算更为全面合理。      

纵向通风公路隧道空气污染浓度适时预报算法

公路隧道空气污染模型的计算通常采用有限体积法,耗时费力,在通风系统的动态调控中不易使用.根据隧道内污染物浓度沿隧道长度近似线性分布的特点,从二阶偏微分方程表示的隧道通风污染模型中推导出了动态预测隧道出口空气污染物浓度的线性迭代方程,对方程中新出现的变量提出了可行的计算方法,同时对方程所描述的物理过程进行了解释．用西汉高速秦岭隧道内的CO浓度的实测数据对方程进行了验证,发现此方程能够在离线的情况下粗略预报隧道出口CO浓度值．为进一步改善预报精度,利用实测数据对上述方程进行修正,提出了一种带有观测误差修正的预报算法,并用现场数据对修正后的预报算法进行验证,结果表明修正后的预报算法能够更准确地在线预测隧道出口实际CO浓度值．     

跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程,选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果,采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据,对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

基于雷达散射截面约束的气动优化设计

采用拉格朗日乘子法优化设计了雷达散射截面约束条件下的锥形融合气动外形.拉格朗日乘子法中的极小化问题采用动态演化的优化设计方法求解.该方法是一种基于非定常演化的优化设计方法,即在求解非定常流动支配方程的时候同时履行优化过程,较其它基于定常解的优化方法具有高得多的计算效率.其中的雷达散射截面通过求解非结构的笛卡儿网格上的时域麦克斯韦方程来得到,而升阻比则通过求解锥形流方程来计算.通过优化设计,得到了M∞=8.0时,升阻比为4.98,雷达散射截面只有1.66m2的锥形融合气动外形.     

应用自由变形法的翼型几何外形参数化研究

对一种新的基于NURBS方法的自由变形参数化方法进行了研究,该方法继承了NURBS方法的优点并避免了其缺点。利用该方法对两副翼型分别进行了参数化和优化设计、反设计,流动控制方程为N-S方程,目标函数梯度计算方法采用了离散共轭方法。算例表明,自由变形参数化方法能够满足设计需要。     

充液挠性航天器的建模与动力学分析

借助于Kane方法建立了充液中心刚体上带有挠性附件的航天器的开链多体动力学模型.采用了混合坐标法和弹簧-质量块液体晃动等效力学模型,所得方程形式简洁.由于引入了2级附件,该模型具有更强的通用性.以跟踪与数据中继卫星为例,给出了数值仿真结果.      

一种TVD方法在三维非结构网格中的应用

以三维非结构网格的显式有限体积法为基础,采用了一种TVD方法求解三维Euler方程,使用Roe通量来计算网格单元边界处的守恒量通量.为了验证方法的可行性,用该方法模拟三维爆炸问题,得出的结论是我们的方法可行,稳定且有效.     

基于粒子群算法的翼型优化设计

采用粒子群算法（PSO）对层流翼型进行了以提高升阻比为目标的优化设计。翼型的设计达到了设计要求,优化设计后的翼型气动特性也有显著地改善,这表明了粒子群算法应用于翼型气动优化设计的可行性。在优化设计的过程中,粒子采用递减惯性权重,以加强粒子初期的全局搜索能力与后期的局部搜索能力。翼型由解析函数线性叠加法表示,目标函数和粒子的适应度由基于二维欧拉方程的流场数值解来提供。