MacCormack格式在计算非定常激波时的特点

前言 MacCormack格式(以下简称为MC格式)近年来在国内外用得较多。在非线性气体力学问题的计算中,它与著名的L-W格式或L-W二步格式相比有许多不同。其主要差别是MC格式的数值解中激波后的过头(Overshoot)较小。对于多维问题,MC格式对时间变量分步(Time Splitting)是十分有效的。 为了更好地了解MC格式,本文讨论MC格式的一些性质,并计算了两个典型问题。     

一种求解双曲守恒律方程的新方法

近年来,在计算流体力学中,WENO(WeightedEssentiallyNon-oscillatory)格式得到了很大的发展,这也促进了基于WENO格式基础上的CWENO(CentralWENO)格式的发展。WENO格式的主要思想是通过凸组合的方式进行插值重构,以得到高阶逼近,并且在间断附近防止伪震荡的产生。本文构造了一种新的五阶紧凑CWENO格式,该格式所选模板节点少,并且保留了原有中心格式的优点。最后通过一些数值算例来验证该格式的高精度,本性无震荡性质。     

ISOMS——计算粘性流场的一种隐式差分格式

用稳定化方法求解可压缩完全Navier-Stokes方程(CCNS方程)是计算复杂流场的一种有效的手段。我们在文献[1]、[2]中利用新的开关函数,构造了一种新的混合格式,称为SOMS格式,用于数值求解CCNS方程。本文是将SOMS格式推广到隐式混合格式,称为ISOMS格式(The implicit second order monotone scheme)。它具有如下几个特点:ⅰ)它是隐式格式,稳定性良好;ⅱ)穿过激波是单调的;ⅲ)对于     

一种MWENO格式的构造和应用

在结构网格上提出了一种新的在混合模板集上守恒量插值的非线性权系数计算方法,构造了一种多重加权实质无波动激波捕捉格式,即所谓的MWENO(Multi-Weighted Essentially Non-Oscillatory)格式.与DWENO(Double Weighted Essentially Non-Oscillatory)格式相比非线性权系数的计算不用人为调节.随后文中给出了几个典型算例.通过MWENO格式和WENO格式的计算结果比较,表明MWENO格式具有较高的间断分辨率.     

基于GPU和隐式格式的CFD并行计算方法

从图形处理器(GPU)架构特点出发,提出了基于数据并行的隐式计算流体力学(CFD)求解方法,空间离散格式采用迎风Roe格式,计算网格适用于结构和非结构网格。采用统一计算设备架构(CUDA)技术实现了GPU上的隐式CFD并行计算。分别在Intel Core2 Quad3.0GHzCPU和NVIDIAGTX280 GPU上进行了计算,结果表明隐式格式计算速度是显式格式6倍以上,采用显式格式的计算加速比达到28倍,采用隐式格式计算加速比达到了28.7倍,同时计算加速比随计算规模的增加而增加。计算结果和实验结果较为吻合。     

非均匀笛卡尔坐标系网格优化频散相关保持格式研究

为克服传统频散相关保持格式 (即 DRP格式 )在数值求解非均匀网格问题时由于寄生波的产生而导致的数值不稳定的问题 ,Cheong &Lee提出了一种网格优化频散相关保持格式 (即 GODRP格式 )。根据其优化思想 ,在非均匀笛卡尔坐标系网格条件下 ,推导得到了更为通用的 GODRP格式优化系数公式。针对 Cheong&Lee文中算例比较不是在同一基准上进行的不足之处 ,设计了一个极端恶劣条件下的非均匀笛卡尔网格 ,采用 GODRP格式和 DRP格式分别求解了二维初始扰动波传播的算例。计算结果表明 ,由于 GODRP格式具有与网格相关的自带粘性 ,不用再另外加入人工粘性便可获得精度较高的数值结果 ,避免了传统 DRP格式在给定粘性时人为经验的影响因素。对于非均匀笛卡尔网格 ,GODRP格式确实比传统 DRP格式更加通用 ,并可以得到更为准确的数值模拟效果      

一种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离散方法对半离散格式进行时间离散,得到时空一致高精度全离散格式。最后,在文中给出几个经典数值算例用以验证本格式的计算能力。     

无波动,无自由参数的耗散差分格式（NND）在喷流计算中的应用

本文用文献[1]建立的无波动、无自由参数的耗散差分格式(NND格式),对无粘喷流流动进行了数值模拟,并就一典型的自由射流流动,与国外高精度、高分辨率的PLM差分格式(Piecewise-Linear Method)的计算结果和实验结果作了比较,得到了很一致的结果,且所用的计算时间少于PLM格式。     

双重加权实质无波动激波捕捉格式的改进和应用

本文重新研究了双重加权实质无波动(DWENO)激波捕捉格式的构造过程，对模板光滑系数的定义，以及两重权值的计算方法进行了改进，并通过一维Riemann问题测试了这些改进，确定了参数的选择范围。最后，本文将改进后的DWENO(r=3)格式进一步应用于一维、二维典型问题的数值模拟，并比较了DWENO格式和WENO格式的精度和分辨率。     

一种混合格式的构建及其在多尺度流动中的应用

面向非定常湍流数值模拟方法精度提升问题,构建了实用的矢通量分裂(FVS)形式高阶精度混合格式。该格式平滑区具有五阶精度线性紧致格式的高分辨率特证,而强间断和强梯度区呈二阶耗散型格式特征,进而保证全流场计算的稳定性。基于所建立的混合格式开展有限体积大涡模拟(LES)计算,将其应用于三维高Reynolds数非定常湍流预测中。结果显示,所建立的数值格式能够可靠预测湍流多尺度涡结构特征,且计算量小,具有较高分辨率和鲁棒性。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。     

超音速流中球头反向喷流流场的数值模拟

采用CSCM(Conservative Supra-Characteristics Method)格式结合激波捕捉法求解了Navier-Stokes方程。通过算例对喷口总压等于2.08倍到6.54倍波后总压等4种情况进行了计算。得到壁面压力分布、等密度线和典型流场速度矢最图。与实验结果和其它数值计算值进行了比较和分析,得到本文计算出的弓形激波脱体距离在喷口压力低时与实验值符合得很好;压力高时差些。     

用CSCM格式模拟二维收缩管道内的流动

 采用Lombard等人提出的CSCM(Conservative Supra-Characteristics Method)差分方法数值模拟二维收缩管道内波系干扰及激波边界层干扰等问题。所得结果同已往计算结果进行了比较。通过压力场和速度场的计算结果可看到Ma数对流场的影响。在Ma数较小时激波角较大,可看到激波经两次相交一次反射逐渐变弱的情况。Re数改变对流场也有明显影响,且给出由于激波干扰引起边界层分离的计算结果。最后还给出收缩管道内激波由正激波逐步形成斜激波的非定常过程。     

Coupling fluid (CFD) and structural (FE) models using finite interpolation elements

In the field of aeroelasticity, interactions between elastic structures and fluid flow are investigated. Recently, numerical aeroelastic models have been built composing those of the combining fluid dynamics (CFD) and the computational structural dynamics (CSD) domains. Since the fluid and the structural models differ in their formulation and discretization, an interface model has to be introduced that represents the connectivity and physical interaction between the two single domain models. In the following, a scheme for coupling fluid (CFD) and structural models (FE) in space is presented which is based on finite interpolation elements. It is applied to static aeroelastic problems, in order to predict the equilibrium of elastic wing models in transonic fluid flow. The structure is represented by finite elements and the related equations are solved using commercial FE analysis codes. The transonic fluid flow is described by the three-dimensional Euler equations, solved by an upwind scheme procedure. The resulting coupled field problem containing the fluid and the structural state equations, is solved by applying a partitioned solution procedure. In each solution step the interface and boundary conditions are exchanged and updated. Here, a fixed-point iteration scheme is used for the coupled aeroelastic equations.     

Multivariate interpolation for fluid-structure-interaction problems using radial basis functions

A multivariate interpolation scheme for coupling fluid (CFD) and structural models (FE) in three-dimensional space is presented using radial basis functions. For the purpose of numerical aeroelastic computations, a selection of applicable functions is chosen: a classical without compact support, and some recently presented smooth compactly supported radial basis functions. The scheme is applied to a typical static aeroelastic problem, the prediction of the equilibrium of an elastic wing model in transonic fluid flow. The resulting coupled field problem containing the fluid and the structural state equations is solved by applying a partitioned solution procedure. The structure is represented by finite elements and the related equations are solved using a commercial FE analysis code. The transonic fluid flow is described by the three-dimensional Euler equations, solved by an upwind scheme procedure.     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。     

跨声速风扇的流固耦合的颤振分析(英文)

发展了时间推进的叶片颤振的数值方法,采用了时间推进求解流体和固体相互作用的过程。气动模型是基于求解三维雷诺平均的N-S方程,采用了多块结构化网格的有限体积格式,对流和耗散通量的计算使用了二阶迎风格式和中心格式。在变形的动网格上流体运动的守恒型方程的求解采用双时间步,隐式格式和多重网格方法。叶片振动采用了振型叠加的线性气动弹性模型。该方法在气动弹性标准算例4进行了验证,并用于求解跨音速风扇的颤振问题。     

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。