非均匀网格湍流大涡模拟高精度有限体积解法

为准确预测不可压复杂湍流,提出了一种可用于大涡模拟均匀或非均匀网格上的高精度有限体积法。该方法空间离散采用有限体四阶紧致格式,时间推进采用四阶Runge-Kutta法,压力-速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动证实了该方法具有近四阶的空间精度;并在此基础上,采用动态Smagor-insky亚格子应力模式,成功地实现了充分发展槽道湍流和后台阶湍流流动的大涡模拟计算,所得结果与直接数值模拟结果吻合良好,且采用非均匀网格可在比均匀网格数少的离散系统上得到同样满意的结果。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

高精度数值方法与非线性扰动研究

通过数值模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算的精确性。本文还研究了波-波与涡-涡两种典型的非线性扰动模型,结果表明:在自由剪切层中,波-波模型能加快扰动的增长速度,并激发出了λ结构,但仍然经过了线性失稳阶段;在钝体头部的超声速边界层中,涡-涡模型通过涡-涡相互作用,产生“发卡”结构,进而出现湍流斑。在此过程中,没有线性失稳的迹象。     

一种MWENO谱体积(SVMWENO5)格式

构造一种新的谱体积(spectral volume)格式来求解双曲型守恒律,记为SVMWENO5(Spectral Volume method byMulti-Weighted Essentially Non-Oscillatory)格式。其主要思想是:第一步,将空间计算区域划分为一系列单元,称为谱体积,等分每个谱体积为一些子单元,称为控制体积(control volume)。第二步,在谱体积内部采用类似MWENO5(Multi-Weighted Essentially Non-Oscillatory)的格式进行重构,而在谱体积的边界处采用传统WENO5格式进行处理。第三步,利用Runge-Kutta TVD离散方法对半离散格式进行时间离散,得到时空一致高精度全离散格式。最后,在文中给出几个经典数值算例用以验证本格式的计算能力。     

广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况，开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质，应用合适的高阶精度格式对各项分别进行离散求解，并采用经典几何声学射线法计算声爆传播路径，实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟，验证了本文高阶离散方法的可靠性。进一步的结果分析表明，采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性，在同等网格数目上可获得更高分辨率的计算结果，且其网格收敛性远优于传统二阶精度格式，同时在计算效率上也有一定优势；在Burgers方程的各效应项中，非线性项的影响更为明显，因而采用高阶离散格式时其优势更为突出；热黏性吸收项对数值结果影响较小，采用高阶离散格式对计算精度的提升效果也不明显，实际计算中仍可以采用传统离散格式，甚至忽略该项的贡献。     

一类新型自适应反扩散近似Riemann求解器及其应用

对于包含激波、剪切层等复杂结构的流动问题，为了精确模拟剪切层等精细结构，且保证激波计算的稳定性，必须采用低耗散且强鲁棒的数值通量方法。传统的HLL近似Riemann求解器的耗散性较大，Roe、HLLEM和HLLC等近似Riemann求解器在计算某些含有强激波的物理问题时会出现非物理解，容易导致不稳定。针对这一问题，本文在Riemann求解器中通过合理设计反扩散矩阵，发展了一类具有自适应反扩散的新型Riemann求解器，并将其应用到高阶加权紧致格式，实现了高阶精度求解。通过典型数值算例验证了新型方法的计算精度和稳定性，结果表明本文提出的新型自适应反扩散Riemann求解器克服了传统Riemann求解器的缺陷，既能准确识别剪切层等精细结构，又能保证激波解的稳定性。     

WENO格式的一种新型光滑因子及其应用

光滑因子是构造WENO格式的关键，决定了WENO格式能否达到最优收敛精度以及在间断附近能否保持本质无振荡特性。针对广泛应用的五阶WENO格式，采用算子函数近似导数关系的方法，设计了一种新型光滑因子。在间断区域，与传统的光滑因子相比，新型光滑因子对间断的识别更准确。在光滑区域，新型光滑因子使得新型WENO格式的非线性权更接近线性权。理论分析和数值验证表明新型WENO格式即使在一阶临界点也能保持一致五阶精度。一维典型激波问题的数值结果表明，与现有WENO格式相比，新型WENO格式提高了短波的分辨率，降低了格式的耗散，同时能够准确识别间断。对于典型包含激波和剪切层的流动问题，新型WENO格式不仅能够准确地捕捉强激波，而且能够精细模拟剪切层和声波等流场结构。     

悬停共轴刚性旋翼直升机涡尾迹高精度数值模拟

针对共轴刚性旋翼直升机流场涡尾迹的高精度数值模拟需求，结合结构化运动嵌套网格技术、高精度数值格式和非定常RANS方法，建立了共轴刚性旋翼直升机的非定常涡流场高精度数值模拟方法及高效并行化自主软件平台。运用所发展的方法和软件平台，进行了类“X-2”构型共轴刚性旋翼直升机悬停状态涡尾迹结构的高精度数值模拟。结果表明，所发展的数值模拟方法具有较高的计算效率和计算精度，能够捕捉到精细的涡尾迹结构及其非定常生成、发展和演化过程。     

一种分裂形式CPR格式在欠解析流动中的稳定性研究

基于重构的修正过程（correction procedure via reconstruction, CPR）方法是一种适用于非结构网格、紧致高效的高阶方法，但在离散非线性对流项时容易因为混淆误差的累积而出现数值不稳定。本文研究了基于LG (Legendre-Gauss)点分裂形式的CPR方法在欠解析流动中的稳定性，并将其与子单元限制技术结合，求解含激波的欠解析流动问题。结果表明：经过边界通量修正的LG点分裂形式CPR格式能够满足离散守恒律，且与子单元限制相结合后守恒律依然成立；在无激波的欠解析流动问题中，相比于散度形式CPR格式，分裂形式明显提高了计算的稳定性，而且比使用LGL(Legendre-Gauss-Lobatto)点的分裂形式CPR格式具有更小的数值误差；在含激波的欠解析流动求解中，相比基于LGL点的间断伽辽金谱元法子单元限制策略，基于LG点的分裂形式CPR格式子单元限制策略具有更高的分辨率和更小的振荡。     

标准化参数定义方法及其工具设计

不同卫星任务的数据特点需要设计不同的遥测数据格式，造成不同卫星任务之间遥测数据格式难以复用，需要重复进行遥测数据格式设计、协调、测试、验证以及地面相应设备的开发和解析处理工作。为解决以上问题，基于包应用标准（PUS）、可扩展标记语言（XML）遥测遥控信息交换（XTCE）的描述方法和XML格式化数据单元（XFDU）的封装方法，提出了参数定义方法，实现了工程参数从类型、参数到数据格式的标准化定义，并设计相应的参数定义工具进行验证。工具可用于航天器设计阶段，设计标准化遥测数据格式及其描述和封装，便于异构系统之间的数据转换。同时也可以用于仿真、验证和任务执行阶段的工作。该方法有助于实现星上端按需、灵活定义数据格式，提升星载应对紧急、未知、非预期等情况时的自主能力。     

抛物积分微分方程的Crank-Nicolson全离散格式下的超收敛分析

文章基于低阶协调的双线性元在矩形网格下的高精度积分恒等式，在时间方向使用具有二阶精度的Crank-Nicolson离散格式，再利用插值与投影相结合的技巧，给出了抛物积分微分方程的全离散格式下的超逼近和超收敛的误差估计。最后，通过数值试验验证了理论分析的正确性。