自适应网格局部加密方法在可压缩流场模拟中的应用

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/     

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

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

利用NURB作曲线和曲面的插值

给出了非均匀有理三次Ｂ－样条插值的方法，利用曲线（面）方程的矩阵表达式导出了反求顶点问题的方程组，使方程组的系数矩阵呈三对角型，易于求解。同时，分别推出几种情况下的端点（边界）条件。利用该结果可以使得曲线（面）易于调整，运算速度快，便于处理。     

基于旋转坐标系的高阶间断有限元方法非定常湍流数值模拟

基于高阶间断有限元方法(Discontinuous Galerkin method,DGM),对旋转非惯性系下耦合了修正的一方程S-A模型的RANS方程进行了离散求解。为了在稀疏网格上获得更贴近真实的物面形状,使用了多层高阶弯曲网格方法对物面进行拟合。非定常时间推进采用了隐式双时间步方法,每个时间步产生的线性系统采用预处理的方法,即广义最小残差方法(Generalized minimal residual method,GMRES)来求解。计算了旋转圆柱绕流以及经典翼型振荡算例的升力和力矩迟滞曲线,与实验结果以及前人的计算结果对比验证了本文方法的正确性和有效性。     

三角网格模型重建中的误差分析与显示

根据实物模型的坐标测量，数据重建模拟的三角网格曲面表示在实践中有着广泛的应用。在许多对精度有较高要求的应用中，对重建网格模型的误差分析是必不可少的。本文提出并实现了一种对重建三角网格模型进行误差分析的方法，该方法首先对网格模型空间进行划分以提高算法的效率，通过计算测量点与网格模型间的最短距离，得到网格模型的误差分布，最后将计算得到的三角网格模型中每个顶点的误差度量转换为相应的颜色表示，并用彩色云图对分析结果进行直观的显示。     

基于混合网格的直升机旋翼/机身非定常干扰流场数值模拟方法

基于混合网格建立了一个适合于直升机旋翼/机身非定常干扰流场的数值模拟方法。其中,采用混合网格系统生成旋翼/机身干扰流场网格以降低网格生成难度;将二阶MUSCL格式和Roe格式相结合进行空间离散以提高流场的计算精度;同时引入隐式LU-SGS时间格式和OpenMP并行模式以提高干扰流场的计算效率。应用建立的计算方法,对Georgia和Robin旋翼/机身模型进行了数值模拟和分析,与试验结果良好的相关性表明了该方法的有效性。     

网格移动定位服务及关键技术

分析了网格移动定位服务的技术定位。在分析当前移动定位服务应用模式不足的基础上,提出了网格移动定位服务的应用模式,并基于移动定位服务的产业链,构建了网格移动定位服务的体系结构。深入研究了面向网格移动定位服务的元数据服务、资源代理服务以及资源副本服务三个关键技术。实现了基于网格移动定位服务的公交导航服务,并阐述了网格移动定位服务的应用前景,故对未来移动定位服务的应用提供理论和技术参考。     

基于非线性有限元降阶方法的网格加筋筒壳结构屈曲承载特性研究

非线性有限元数值计算方法是结构屈曲性能评估的主要手段,然而该方法需要基于精细网格迭代求解,降低结构非线性计算规模是其关键技术之一。采用一种结构非线性有限元降阶方法,对均匀轴向压载荷作用下的网格加筋筒壳结构进行了细致的承载稳定性特性分析。该方法将Koiter初始后屈曲理论与Newton弧长法相结合,通过多次使用Koiter摄动展开来跟踪结构非线性平衡路径。计算比较了特征值屈曲载荷和计及几何非线性效应的非线性屈曲载荷值;采用侧向小扰动载荷模拟结构初始几何形状缺陷,分析了缺陷幅度对临界屈曲载荷的影响规律,在验证了非线性有限元降阶方法分析精度的同时,细致考核了该方法在结构非线性屈曲分析以及缺陷敏感性分析中的计算效率优势。     

发动机进气流动三维瞬态数值模拟研究

采用任意拉格朗日-欧拉法(ALE法)和所提出的动态网格生成新方法,实现了车用发动机进气流动过程的三维瞬态可压缩湍流数值模拟,湍流模型采用经压缩修正的κ-ε双方程模型.计算结果揭示了发动机进气过程的三维瞬态流动特性,为发动机进气系统的优化设计提供了理论依据.     

非交错网格上不可压流动中的近似投影方法

对于压力增量形式的投影方法,在非交错网格上,仅仅采用近似投影过程,通常不足以抑制压力场的数值振荡.本文通过选取合适的投影变量,构造了一种可以在非交错网格上避免压力振荡的新方法.该方法在不增加计算量的前提下,让前一时刻的压力参与投影过程,有效地割断了前一时刻中可能存在的压力数值振荡向下一时刻传递的途径,避免了数值振荡的积累,提高了数值方法的稳定性.此外,我们在投影变量中考虑了"几何网格雷诺数"的影响,拓宽了本方法的适用范围.数值试验验证了本文提出的近似投影方法是可行的.