紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

差分格式的优化组合

用偏微分方程组解实际的三维紊流流动问题时,如何选择各方程中对流项和扩散项的差分格式将对计算的稳定性、准确性与经济性影响很大.对一个实际的流动问题,计算机模拟验证了迎风差分、中心差分、混和格式、PLDS(Power Law Differencing Scheme)、QUICK(二次迎风插值)和OPTIMAL(Oscillation Preventing Technical algorithm based on the second-order Interpolation Method for Advection Linked convection)的各种组合对计算的稳定性、准确性与经济性的影响.对中心差分应用在雷诺时均方程中出现的数值解的振荡问题提出了在局部网格使用迎风差分的具体对策.对QUICK应用在κ,ε输运方程上出现的数值解的发散问题提出了QUICK和迎风差分格式交替使用的具体方法.      

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

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

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

激波捕捉法中信号传输的监控和检测

依据本文引入的模拟流场中信息传输的监测因子θ,能十分容易地判定流场的激波。为了提高捕捉激波的质量,文中提出了通过选择合适的限控函数R(θ)构造高精度TVD差分格式的原理及利用由R(θ)提供的信息,用最小二乘法确定激波位置的设想。     

计算气动声学中的高精度紧致差分格式研究

优化了五对角紧致差分格式.通过Fourier分析,将优化目标转化为一个求带有约束的多元非线性函数的最小值问题,利用序列二次规划(SQP)方法获得最佳系数.通过三种措施保证优化的格式具有高精度和分辨率:①直接对量化波数误差积分求其最小值;②采用绝对误差准则,使各种波长的波具有同一误差限;③优化的波数空间和精确求解区间一致.通过调整内部和边界格式的Taylor精度及误差限,保证整个格式的稳定性,并从理论上证明了优化格式具有渐进稳定性.一维和二维基准算例体现了优化格式的性能改进.      

三阶中心无振荡格式性能分析

三阶中心无振荡格式Cn3(Centered nonoscillatory scheme of third order)使用对称模板构造具有三阶精度的插值公式.利用单调区域与精确区域修正原始插值参数,该格式能够获得间断附近无振荡、光滑区域高精度的计算结果.通过一维和二维典型算例,将Cn3格式与三阶和五阶WENO格式(Weighted Essentially Non-Oscillatory schemes)进行比较,重点分析了3种格式的间断分辨率、计算稳定性和数值耗散性.分析可见,Cn3格式能够精确、稳定地捕捉激波和接触间断,同时对光滑流动区域的小尺度流动结构保持较低的耗散,值得进一步研究及推广应用.     

二维Euler方程组的高分辨率隐式差分格式的并行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二维Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合于多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。     

非线性双曲型守恒律的一类应用通量限制器的二步二阶格式

本文提出了计算非线性双曲型守恒律弱解的一类新的二步二阶精度有限差分格式,这类高精度非线性差分格式的构造应用了三点单调差分格式的数值通量和通量限制器。文章证明,在拟CFL条件限制下,格式为总变差衰减差分格式(TVD格式),因此不发生常用的诸如二步二阶Lax-Wendroff格式、Maccomark格式和其他古典二阶格式计算不连续解时在其周围出现的剧烈振荡现象。 为了减少Roe推广技术引起的计算复杂性,文中我们基于Davis,S.F.的思想,使用了Rusanov数值通量,把标量差分格式推广到方程组,构造了求解双曲型守恒律方程组的新的预测-修正格式,较之其他一般的预测-修正格式而言,仅增加了很少的计算复杂性,初步的数值试验表明,本文提出的格式无振荡和具有极好的激波分辨率。