可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

基于PIV速度场测量的压强梯度计算

介绍了基于PIV速度场测量结果计算压强梯度的基本原理和3种不同的计算方法,并通过对流高斯涡模型系统研究了离散格式选取、PIV参数设置和流场特征等对压强梯度计算的影响.结果表明:PIV速度测量的精度是压强梯度计算的关键性影响因素,1%噪声对结果的影响比其他过程高约2个数量级;时间分辨率和空间分辨率的设置并非越高越好,而是有个合适的中间值;拉格朗日方法在大多数情况下优于欧拉方法,但当流动结构较复杂时性能较差,应根据流动类型和流场特征合理选取压强梯度计算方法.      

Revisiting the space-time gradient method: A time-clocking perspective,high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternativ...     

多重网格技术在复杂粘性流场计算中的应用及研究

流体力学控制方程的数值求解过程中,当网格加密、粘性效应强或来流马赫数低的时候,流场收敛非常缓慢。本文在多块网格、LUSGS隐式格式基础上采用多重网格技术,计算了NLR7301case#1两段翼型的低速绕流及DLRF4跨声速流场,在复杂二维、三维粘性流场求解中获得了满意的加速收敛效果。文中还比较了V循环、W循环的加速效率,为多重网格的工程应用奠下基础。     

采用改进的CE/SE方法模拟方管中氢氧爆轰波的稳定传播结构

采用改进的时空守恒元/解元（CE/SE）格式和新型二步化学反应模型，对小扰动情况下方管通道中的氢氧气相爆轰传播过程进行了数值模拟，对不同截面尺寸和两种扰动模式下的爆轰波三波线结构进行了详细计算和分析。数值模拟花费了相对较少的计算代价，清晰地得到了方形截面通道中3种相对稳定的爆轰波传播结构，即直角模式、对角模式和螺旋模式。随着截面尺寸的减小，爆轰波胞格图样与波阵面结构随之产生适应性的变形，乃至出现胞格数目的减少和部分三波线的消失；最终，直角模式和对角模式都会转化为顺时针或逆时针的螺旋模式；临界尺寸下，得到了转换过程中新的过渡结构和压力脉冲振幅的变化，并以此初步讨论了直角和对角传播模式的稳定性。     

基于带波动算子非线性Schr(o)dinger方程数值分析的守恒差分算法

对一类带波动算子的非线性Schr(o)dinger方程进行了数值分析,提出了一个含参数的二阶守恒差分格式,根据参数选取的差异,该格式既可隐式计算也可显式计算.对初值条件进行了中心差分离散,使其具有二阶精度,从而与守恒格式的精度一致.利用矩阵理论证明了差分解的存在惟一性,并利用一个重要的不等式在先验估计的基础上,运用能量估计的方法证明了该格式按无穷范数以二阶精度收敛到真实解.数值实验表明该格式具有较高的计算效率.     

一种改进的MUSCL格式

给出了一种改进的二阶精度的ＭＵＳＣＬ格式，并用来求解Ｅｕｌｅｒ方程，改进后的ＭＵＳＣＬ格式在线化插值步中，用Ｖａｎ Ｌｅｅｒ的插值斜率对特片变量进行插值，并对与熵波、剪切波对应的特片变量使用了小耗散的插值斜率，在迎风步采用Ｒｏｅ的近似黎曼解求得迎风通量。比较一维和二维的计算结果均可以看出，改进后的ＭＵＳＣＬ格式在消除数值振荡的同时，明显提高了流场间断的分辨率。     

基于数据库的柔性MIS系统初探

本文针对传统MIS系统的局限,介绍了柔性MIS系统的概念和定义,着重讨论了柔性MIS系统的总体结构和各组成模块的功能及设计要点。主要阐述了库文件结构的定义及自动维护的原理和方法。给出了通过功能模块语句中条件成分的表达及处理方法、相关性的描述方法、屏幕记录界面和格式报表的生成原理、字段的属性分类以及关键字段、公共字段、代码字段的描述及处理。还简要涉及了各模块间的信息传输。