高阶精度DG/FV混合方法在二维粘性流动模拟中的推广

DG/FV混合方法因其具有紧致性、易于推广至高阶及相比同阶DGM计算量、存储量小等优点,已成功应用于一维/二维标量方程和Euler方程的求解。在此基础上,将该方法推广于二维三角形/矩形混合网格上的NavierStokes方程数值模拟,将格式形式精度提高至4~5阶。物理量的空间重构及离散使用DG/FV混合重构方法;无粘通量计算采用Roe格式;粘性通量计算采用BR2格式;时间方向离散采用高阶显式R-K方法或隐式方法。利用该方法计算了有解析解的Couette流动问题以验证几种格式的数值精度阶,并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明DG/FV混合方法达到了设计的精度阶,在较粗的网格上亦能得到高精度的计算结果;定性分析和数值结果表明相比同阶DG方法单步计算量减少约40%。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法．在背景网格建立和三维初始推进面生成的前提下，利用推进阵面法生成了翼身组合体的三维空间网格，并应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

超声速底部喷流干扰流场数值模拟

数值模拟了不同马赫数,不同喷流压比下的轴对称超声速底部喷流干扰流场,采用LU隐式算法进行数值求解并引入了Baldwin-Lomax代数湍流模型.采用分区网格将弹身与底部区域合为一个整体进行计算,得到了清晰的流场结构和弹体表面及底部的压力分布,并与试验结果进行了比较.数值模拟结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同, 喷流对底压分布有明显影响,进而对轴向力系数影响显著.     

基于自适应混合网格的脱体涡模拟

基于混合网格和CGNS（CFD General Notation System）数据结构,建立了一种各向同性加密/稀疏的网格自适应方法。在悬空点的后处理中,让含有悬空点的单元转化为任意多面体,从而简化了自适应单元剖分模版,同时自适应网格单元之间可完全相容,自适应生成的网格能够直接用于可处理任意多面体的流场求解器。将该自适应方法与脱体涡模拟（DES）算法相结合,开展了65°后掠三角翼大迎角流动的数值模拟应用,并与初始网格的模拟结果进行了详细比较。对比表明：采用网格自适应方法适当增加局部网格量,能够以较小的成本迅速提高三角翼背风区的空间分辨率,增强数值模拟对小尺度涡系结构的解析能力,从而弥补了基于混合网格的脱体涡模拟中常用二阶格式计算的空间分辨率相对偏低、不利于湍流多尺度结构精细模拟的不足。     

固体火箭发动机喉衬用轴编C/C复合材料的细观热结构特性分析

为研究固体火箭发动机喉衬用轴编C/C复合材料的细观热结构特性,以组分材料之间界面分析为基础,完成热结构参数的实验测定与基于代表性体积单元的等效热物理参数的预测,得出轴编C/C复合材料的等效热膨胀系数与等效热导率系数。分析表明:基于代表性体积单元程序温度周期性边界条件的应用,可以较精确预测复合材料等效热膨胀系数和等效热传导率;得出了组分材料界面相对等效热结构参数的影响,组分材料界面相采用一定厚度的单元模拟更加接近实验数据;讨论了等效热膨胀系数和等效热导率随编织参数的变化规律。     

氧浓度对同轴乙烯层流扩散火焰中碳烟形貌及粒径演变的影响研究

为探究氧体积分数对乙烯同轴射流扩散火焰中碳烟颗粒生成及演变过程的影响,采用SiC纤维沉积法和热泳探针采样法,对层流扩散火焰中不同径向和轴向位置处的碳烟生成特性进行了研究。研究表明,氧体积分数增加,使火焰同一位置温度升高。在氧体积分数低于31%、火焰高度低于30mm时,SiC纤维上碳烟沉积物形态由火焰中心位置处的表面光滑的类液态演变为粗糙的凹凸块状,随温度升高及氧化作用加强,渐渐被氧化为疏松的团簇状、絮状,最终过渡为致密的团簇状和纤维网状。相同氧体积分数下,碳烟颗粒平均粒径随火焰高度增加呈现先增大后减小的趋势。在氧体积分数21%、火焰高度30mm处初生粒子直径达到最大值,为41.8nm。在火焰根部,碳烟颗粒平均粒径随氧浓度升高而增大,而火焰较高位置处则呈现了相反趋势。     

基于HLLC格式的理想MHD方程数值模拟

针对理想MHD(magnetohydrodynamics)方程数值求解困难的问题,基于原始的HLLC (Harten-Lax-Van Leer Contact wave)近似黎曼解方法,发展出一种新的适用于理想MHD问题的通量计算格式。控制方程采用有限体积法离散,时间推进采用隐式的LU-SGS（lower-upper symmetric Gauss-Seidel）格式,并且采用双曲型散度清除技术来抑制磁场散度的累积。通过一维激波管问题的数值模拟表明,HLLC-MHD格式能准确的分辨并捕捉复杂的磁流体力学波系结构,保证数值计算的精度,最大的数值计算误差不超过10%;通过二维的Rotor问题的数值模拟研究表明,HLLC-MHD格式能够应用于多维理想MHD问题的数值模拟,并且能够准确捕捉多维情况下磁流场中的阿尔文波;通过对比有、无散度清除的Rotor问题计算结果表明,双曲型散度清除技术可以将磁场散度峰值从50降低到2,有效抑制高磁场散度区域的散度累积,同时也会将误差传播到低散度区域,并且引起边界处散度的累积,影响计算的稳定性。      

二元喷管跨音速流场的数值计算

本采用有限体积法进行空间离散，Runge-Kutta法进行时间推进的办法，通过求解Euler方程来模拟二元收敛－扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡，在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解，计算结果与试验数据和其他数值方法吻合良好。     

氰酸酯/环氧树脂的湿热老化机理

采用热失重、热机械分析、正电子淹没技术(PAL)对氰酸酯/环氧树脂体系湿热老化前后的性能进行了测试,对体系的吸水机理进行了研究。结果表明:该体系的吸水率为1.8%,湿热老化后对体系的热失重温度没有影响,体系的储存模量降低了95MPa,Tg降低了约7℃,损耗峰变宽并向低温方向移动。湿热老化后自由体积尺寸减小,自由体积浓度增大,体系中的水分起着增塑剂的作用,导致基体中形成更多的微裂纹。     

准一维可压缩瞬变管流的有限体积模型(Ⅰ)流场的有限体积模型

针对现有流体系统一维动力学模型缺乏统一的理论推导过程和能够涵盖众多模型的基本方程组的现象,从适用于连续介质控制体的欧拉型积分形式的流动守恒方程出发,经过合理模化推导得到了可压缩变截面管流准一维瞬变流动积分形式和微分形式的守恒方程.通过保留能量方程中的变体积项考虑了由于控制体体积的改变所输出的膨胀功.通过对守恒方程在空间交错网格上的有限体积离散建立了一种计算准一维可压缩变截面管流瞬变流场的有限体积模型,该模型可考虑变体积、变物性、轴向热传导、重力场影响,一方面,在深化和阐明了有限元状态变量模型体系的理论基础的同时,拓宽了其对瞬变流动的适用范围,另一方面,结合阀芯节流模型,从此模型出发可推导出管路系统常见元件的流场模型.