高精度ENO格式的有效实现

基本无振荡（ＥＮＯ）有限差分格式是近几年出现的一类高精度差分格式，重构多项式的形成是ＥＮＯ格式构造的核心部分，一般均以网格平均值作为插值点。文中由通量值作为插值点，构造了一个简单的高精度格式，在求解高维欧拉方程时，较节省计算时间。引入三 阶ＴＶＤＲｎｇｅＫｕｔｔａ时间离散，就二维欧拉方程的激波反射与双马赫反射问题进行了数值试验，试验结果表明格式所具有的高精度和基本无振荡性质，能够满足实际计算的需要     

Euler方程无网格算法及布点技术

研究了Bitina的显式无网格算法。在该算法的基础上,给出了Euler方程无网格离散形式,运用RungeKutta显式时间推进格式推进求解。文中的耗散项采用了与非结构网格上类似的守恒型耗散算子,边界条件的处理借鉴了结构化网格处理技术。此外,本文还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了二维翼型的典型绕流。     

用任意拉格朗日--欧拉方法求解流体力学方程组

从积分形式的二维Lagrange流体力学方程组出发，用有限体积格式进行计算，针对不规则的四边形网格，在用次级网格的方法确定网眼内物理量的梯度的基础上，提出一种体积加权的方法来构造网眼内的线性插值多项式。将新网格进行了细划，先确定新网格小网眼中心点的密度，将小网眼的质量直接计算出来，再计算新网格的质量，然后确定相应的密度，从而实现高精度重映。最后用ALE方法进行数值模拟，得到了预期的效果。     

基于结构网格点的Euler方程无网格自适应求解算法

提出了一种自适应无网格Euler方程求解算法。构造了基于无网格点云结构的布点加密技术,实现了借助压力梯度变化准则的流场局部加密;借助流场结构化剖分技术,给出初始无网格节点分布,并采用四步显式Runge—Kutta时间格式推进计算,求解了绕翼型的典型流动问题。数值算例表明,本文提出的方法能有效地提高流场的分辨率,如捕捉激波等流动特征清晰明了。     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

基于扫掠法的六面体网格生成算法及实现

提出了一种基于扫掠法的六面体网格生成算法,该算法通过源和目标曲面参数域间的近似最小二乘仿射映射将源曲面的网格投影到目标曲面,并利用背景网格插值的方法计算内结点,从而将2.5维实体生成六面体网格。该算法避免了空间网格投影过程中耗时的“寻根”过程,并且不需要导致效率降低的网格优化,能够处理各种2.5维实体,生成的网格质量较好。     

非结构网格生成技术

描述了非结构网格生成的过程,网格的生成方法基于Delaunay 三角形方法,并对物面附近网格进行拉伸以便获得在边界层内高拉伸比的三角形单元。网格的生成过程分两部分:(1) 按Delaunay 三角形方法生成适于无粘情况下的非结构网格。(2) 对近物面的部分三角形单元进行拉伸,使生成的网格能适用于粘流。使用上述方法本文构造了二维弯管内绕叶栅区域的非结构网格     

基于基因算法与博奕论的翼型跨音速Euler方程气动优化

将基因算法与博奕论中的 Nash平衡相结合 ,构建了多目标基因优化算法—— Nash基因算法 (NashGAs) ,并对初始翼型为 NACA0 0 1 2二维翼型进行给定跨音速流动下的形状增升优化。计算中应用 Bézier曲线对翼型形状进行参数化 ,避免了非流线型的不合理形状产生 ;采用有限元非结构网格 AUSM+通量分裂格式Euler方程数值解进行个体适应度值评估 ;应用动态网格技术调整计算网络 ,节约了 CPU机时 ;最后给出了优化结果     

二维网格局部动态细化的关键技术及其在切削仿真中的应用

为了解决仿真过程中网格畸变导致的计算终止以及计算精度与效率低等问题,本文利用Python对ABAQUS有限元仿真软件进行了二次开发,提出一种二维局部网格动态细化算法。该算法通过研究网格细化准则和细化方法以及物理场的传递过程,完成了整个算法程序的开发,并应用该算法实现了在二维切削仿真中局部网格动态细化。与全局加密网格的模型相比,在保证仿真精度误差为5%以内的情况下,采用局部网格动态细化时计算效率提高了210%。最后,通过试验验证了采用该算法所建立的二维切削仿真模型的准确性,仿真结果与试验结果基本吻合。     

一种基于块结构的非结构网格快速生成方法

以阵面推进法为基础,提出了一种用预制模块的方法实现在平面域内生成非结构网格的新方法,即模块方法。模块是由一组三角形构成,文中以三角形模块为例,介绍了网格模块的预制方法。网格的生成是根据空间步长的背景信息,从预存的网格模块中选取当地合适的模块,直接复制或变形到计算区域实现的。由于生成过程中每一步引入的是预制的网格模块而不是传统的单个三角形单元,因此可期望新方法要比传统方法快得多。算例显示了新方法的快速网格生成过程,该方法具有推广应用前景。     

三维非结构网格的生成及其应用

采用阵面推进法生成了高质量的三维四面体非结构网格，该方法的特点是在生成单元的同时引入结点，与其他预布结点的方法相比，其好处是能更好地控制网格单元尺度在整个计算域的合理２并提高网格质量。用Ｄｅｌａｕｎａｙ方法连接存储有控制参数的结点生成了背景网格。网格生成过程中采用了合理的数据结构，有效地提高了程序运行效果，最后给出了Ｍ６机翼的数值算例。     

GENERATION AND APPLICATION OF UNSTRUCTURED ADAPTIVE MESHES WITH MOVING BOUNDARIES

INTRODUCTIONMany engineering applications require thenumerical prediction of strongly unsteady flowsinteracting with moving bodies.A case in pointis that the detailed calculation of sphere motionsin an unsteady field induced by blast waves is akey issue to determine the drag force exerting up-on the sphere by the surrounding gas and,inturn,to reveal the mechanism of unsteady effectson the drag coefficient.To cope with such pro-blems,one should simulate unsteady flows withmoving boundary co…     

NEW APPROACHES OF MUSCL SCHEMES FOR COMPRESSIBLE FLUIDS

The MUSCL scheme for compressible gas dynamicsis studied in this paper. We propose a new type of Euleri-an MUSCL scheme, which evaluates the intercell flux athalf time step. The intercell flux is comptuted from char-acteristic equations with the data which traced backthrough the approximatecharacteristic from the edge athalf time step. The data is classified in several differenttypes, depending on the characteristic directions of data.We also present a general procedure of oscillation-fre…     

双曲型守恒律的无振荡中心差分格式

讨论了双曲型守恒律的一类无振荡中心差分格式。Ｈ．Ｎｅｓｙａｈｕ和Ｅ．Ｔａｄｍｏｒ研究了以交错型ＬａｘＦｒｉｅｄｒｉｃｈ格式（ＬｘＦ）为模块的无振荡中心格式的构造与熵不等式。此类格式利用高阶的ＭＵＳＣＬ型插值替代一阶分片常数逼近，减少了ＬａｘＦｒｉｅｄｒｉｃｈ格式的过多数值粘性，建立了一维标量非线性双曲型守恒律的一类高分辨格式。讨论以非交错ＬａｘＦｒｉｅｄｒｉｃｈ格式为模块建立起的差分格式。证明了此格式具有二阶精度、ＴＶＤ性质并在一定条件下满足熵条件。     

旋转坐标系下三维MUSCL格式的构成

MUSCL格式是近年来新发展的许多高分辨率格式之一，但直接将该格式应用到旋转坐标系时，会使计算发散或得到非物理解，不能计算旋转机械的内部流场。本文在原有MUSCL格式的基础上，根据旋转坐标系下雷诺平均N－S方程，重新构成了MUSCL格式，并且与原有MUSCL格式进行了比较，分析了两者相同与不同之处，利用该格式对两个典型算例（旋转直管和跨音速压气机）进行了计算，通过与试验结果对比表明，重新构成的MUSCL格式能够计算包括跨音速气机在内的复杂旋转机械流场，从而扩展了MUSCL格式的应用范围。     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

真实火箭燃气流场的测量与数值仿真

对用激光－莫尔偏折法测出的真实火箭燃气射流流场进行了数值仿真，在计算中采用了二阶精度的ＭＵＳＣＬ格式，将三维问题转化为二维问题，结果表明，理论计算与实验测量结果符合较好。