求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

用虚拟区域法结合虚拟边界条件求解对称及管壁流动问题

分析了用基于Lagrange乘子的虚拟区域法数值求解时,翼型后缘存在的奇异问题.为了保持直角坐标网格结构简单及可利用快速算法的特点,提出引入虚拟边界条件控制虚拟流动的方法处理奇异问题.用虚拟区域法数值模拟了翼型不可压对称绕流;并把该方法推广应用到管壁流动的数值模拟,给出了喷管内流场的算例.数值结果表明,虚拟边界条件的引入有助于控制Lagrange乘子的变化,从而改善数值模拟的精度.     

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

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

基于局部无网格的混合算法研究

研究了一种求解任意外形绕流的局部无网格混合算法.该算法计算区域整体采用了直角网格,只在物面附近嵌入局部无网格处理.整体直角网格的规则性与正交性,不仅使得网格生成快,而且计算格式实施简单,避免了与非正交性相关的计算项,从而引入误差相对较小.局部无网格处理只要求物面附近布点离散,具有灵活性,使得发展的算法适合处理任意外形.先以理想不可压流动问题为例,给出了其控制方程Laplace方程的具体求解实施过程,并给出了直槽道内柱体绕流的计算结果;接着把方法推广用于求解可压缩流动的Euler方程,对喷管等典型流动问题进行了数值模拟,并与传统有限体积法结果进行了比较,在捕捉的波系结构、表面压力系数分布等方面具有良好的一致性,表明发展的新方法是可行的.     

动边界问题无网格算法及其动态点云

提出基于Delaunay图映射的点云运动处理技术.该技术使点云随物体运动再生成问题变成简单的映射关系,无需迭代,且能保持点云结构不变,故有效地避免了因点云结构变动而产生的附加计算量.运用该技术,结合非定常Euler方程双时间步无网格算法,成功地模拟了翼型的跨声速非定常流动.其中物理时间迭代采用二阶隐式格式,伪时间遮代采用四步龙格一库塔显式格式.文中给出了不同翼型的计算结果,并与实验结果进行了比较,取得了令人满意的结果.