基于非结构动网格的非定常激波装配法

从ALE方程出发,通过修正体积的新思路构建离散几何守恒律,基于非结构动网格的信息传递方法得到网格重构以后新网格的流场参数,采用一阶格式捕捉法得到的流场确定装配法的初始激波位置,成功地把激波装配法应用到非定常流动的模拟。模拟超声速飞行环境下整流罩分离的非定常流固耦合过程,在马赫数6条件下装配法和捕捉法的计算结果完全一致,随着马赫数增加,二阶精度格式的捕捉法需要进行调整限制器参数、时间步长等人工干预,但是装配法直到马赫数20计算过程没有出现异常。该基于非结构动网格技术的激波装配方法,与传统的激波装配方法相比,能够用于复杂外形产生的非规则形状激波和非定常流动产生的运动激波,解决二阶精度有限体积法用于高马赫数飞行条件下多体分离过程的流动模拟时遇到的稳定性问题。     

基于离散等价方程的非结构网格有限差分法

激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

特定条件下高阶WENO格式计算结果误差

通过比较一阶迎风格式和五阶WENO格式模拟激波、接触间断、曲线坐标系下的均匀流和激波正规反射等4个简单流场得到的数值结果,发现WENO格式模拟的激波和接触间断在从初始间断变化成数值过渡区的过程中出现的非物理波动比一阶迎风格式的结果更加明显,流场结构也更加复杂;同时,由坐标变换而产生的几何诱导误差和边界近似模型误差也明显...     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

一种改进的非定常激波装配算法

基于非结构动网格技术和格心型有限体积方法,提出一种改进的非定常激波装配算法,进一步拓展了其在包含有运动激波的非定常流场的应用范围。首先,针对激波在直/曲壁面传播这类问题,分别建立了壁面间断节点的运动模型;其次,为保证激波在大范围运动时装配阵面不产生失真,基于Bézier曲线拟合方法实现了间断节点分布的自动重构;接着,通过嵌入局部网格自动重构模块,提高了算法的计算效率和自动化程度;最后,对于激波相交点的运动,设计了一种根据位移推算速度的方法进行装配。数值算例表明,所提算法能够有效地处理激波传播问题,相比激波捕捉方法可以提取更多的流场信息,同时可以获得流场间断更加直观清晰的图谱。     

捕捉间断的高精度数值方法

为发展适用于捕捉超声速流场中各种间断的高精度算法,将通量限制的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的通量限制型差分格式.通过在时间方向上利用一阶精度格式计算的一维定常激波,以及在时间方向采用多步Runge-Kutta方法计算的一维非定常激波管问题上的数值试验与二阶精度的TVD格式所计算的结果比较,表明新方法比二阶精度方法在间断的捕捉上具有明显的优势.通过新方法的计算结果与精确解的比较,表明新方法的准度也是非常令人满意的.     

用于气动声学计算的非均匀网格紧致差分格式

为克服传统紧致差分格式在数值求解非均匀网格问题时产生的寄生波,构造了一种新的高精度紧致差分格式.通过泰勒展开分析方法,详细给出了格式系数的通用形式;利用傅里叶分析方法,分析了数值耗散、色散误差.以3对角6阶精度紧致差分格式求解均匀扰动网格问题为例,计算表明:色散值和耗散值随扰动因子的增加而更加趋近于精确值;当扰动因子大于0.213时,格式不稳定,当扰动因子小于等于0.213时,格式渐近稳定;对一维对流波和二维波传播的模拟计算所得数值解与精确解吻合,体现了该格式在求解非均匀网格问题时的优越性.      

NS方程激波计算的摄动有限差分方法

摄动有限差分(PFD)方法从一阶迎风差分格式出发,将差分系数展开为网格步长的幂级数,通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式,因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导,将摄动有限差分格式引入时间相关法的计算,并构造了守恒形式的摄动有限差分格式,然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明:本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率,而且保持了一阶迎风格式的无振荡性质。     

一种用于间断装配的流场结构辨识算法

使用激波装配法时,初始激波是否准确将会对计算过程产生影响。为了确定初始激波的位置,提出了一种新的流场结构辨识算法。该算法以捕捉法计算得到的流场作为系统观测数据,根据密度、压力等参数从该数据中获取激波和接触间断等流动特征周围的网格节点作为离散点集。通过将该离散点集分割成若干子区域,在各子区域内进行分片拟合,最终将离散点集拟合成连续光滑的实体模型,并将此作为初始激波面。在二维方法的基础上,通过引入单位球模型成功将该辨识算法拓展到三维应用。结果表明,采用该方法获得的间断曲面(激波和接触间断)与捕捉法流场中的间断分布吻合较好,作为初始间断面用于装配法可快速得到收敛解。该方法解决了应用激波装配法时确定初始间断面的难题。此外,该方法还可用于网格自适应方法。选择不同流动参数,可以获得相应流场特征结构的空间曲面,在此曲面的基础上可进行网格局部加密或重剖分。该流场结构辨识算法用于网格自适应具有网格尺度自由设置的优势。