处理动边界问题的无网格/直角网格混合算法

把无网格/直角网格混合算法发展用于涉及动边界的非定常流动计算问题.采用无网格/直角网格混合算法,计算区域整体采用直角网格,只在物面附近嵌入局部无网格.对于动边界问题,本文让主体直角网格保持固定,用原混合算法中的局部无网格区跟随运动边界,使得动边界物面附近有不变的点云结构,避免了点云重构;无网格区与直角网格交界面随边界运动跟踪确定,由于主体基于直角网格,方便了交界面的跟踪.用发展的混合算法,结合非定常Euler方程双时间推进方法,先对翼型作俯仰振荡的非定常跨声速流动进行了数值模拟,计算结果与现有文献和实验数据进行了比较与验证;结合刚体运动方程,成功模拟了二维外挂物投放过程,算例展示出发展的混合算法适合处理大位移运动边界问题.     

一种基于三阶WENO重构的无网格算法研究

提出了基于三阶WENO重构的无网格算法,以期替代传统的无网格线性逼近重构,提高无网格算法的求解精度。基于无网格点云卫星点分布特征,通过沿点云中心点与卫星点连线方向引入局部一维坐标系,并结合虚拟点的设置,构成了在点云中实施三阶WENO重构所要求的模板。算法涉及的卫星点连线中点处的左右状态值,则利用模板上各点的流场值,采用WENO重构方法计算给出。对于模板中设置的虚拟点上的流场值,本文则利用已有的无网格点云信息,基于最近节点流场值插值得到。文中给出了采用上述WENO重构方法,结合Roe的近似Riemann求解器确定通量,并采用三阶TVD Runge-Kutta方法进行时间推进求解Euler方程的具体实施过程。对模型问题的典型流动进行了计算分析,通过比较数值解和精确解,验证了算法获得的数值解能逼近三阶精度。给出的激波管流动和超声速半圆柱绕流算例展示出所发展的算法捕捉激波等间断问题具有更高分辨率,能体现出WENO重构"基本无振荡"的特性。文中最后给出了非定常激波过弯道绕双圆柱流场的数值模拟算例,展示了所发展的算法处理含非定常激波演化的复杂流动问题的效果。     

模拟存在运动激波流场的无网格自适应算法

将无网格自适应算法发展用于求解存在运动激波的非定常流动问题。为了降低初始布点要求,同时提高对运动激波的分辨率,提出了一种结合点云加密和点云粗化的动态点云技术。采用该技术,在非定常计算的每一个时间步,首先通过基于压强梯度的流场探针识别出激波所在位置,然后对激波附近点云进行加密,从而提高对激波的分辨率,同时对远离激波的点云进行粗化,尽可能减少冗余节点,提高计算效率。采用所提出的无网格自适应算法求解了N-S方程,首先对一维运动激波问题进行了数值模拟,通过将数值模拟结果与精确解比较,验证了所提算法求解存在运动激波的非定常流动的可行性,在此基础上,将方法推广用于二维非定常流动问题的求解,分别对瞬时激波碰撞圆柱以及双弧翼型激波振荡流动进行了数值模拟,计算结果表明,方法可以显著提高对激波的分辨率,同时相对于直接采用较密布点的非自适应算法,自适应算法在计算效率上具有优势。     

用于非定常粘性流动计算的无网格算法(英文)

把无网格算法发展用于求解涉及动边界的非定常粘性流动问题。在处理粘性流动的布点问题时,通过在远离物面的区域采用无粘流动计算时所用的各向同性点云,而在物面附近引入各向异性点云,从而在保证物面法线方向布点较密以准确模拟边界层的同时,有效控制了布点总量,减少了计算时间。对于动边界问题,本文在上述方法获得的初始布点的基础上,采用基于扰动衰减规律的点云移动技术,快速得到物面运动到其他位置时流场求解所需点云。在所获得的点云上,采用一种带权系数的二次极小曲面逼近方法来离散Navier-Stokes方程的空间导数,权系数的引入使得流场求解时对各向同性点云和各向异性点云无需分开考虑,而是采用统一的求解方法,从而简化了编程求解。用发展的无网格算法,结合Navier-Stokes方程求解双时间推进方法,并耦合Spalart-Allmaras湍流模型,先成功地模拟出NLR7301翼型后缘摆动诱发的非定常流,并通过与实验比较,验证了本方法.接着进行了NACA0012失速模拟,给出了动态失速过程,展示出用本文方法处理复杂动边界问题的效果。     

HYBRID CARTESIAN GRID/GRIDLESS METHOD FOR CALCULATING VISCOUS FLOWS OVER MULTI-ELEMENT AIRFOILS

A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils. The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leaves only small region adjacent to the aerodynamic bodies to be filled with the cloud of points used in the gridless methods, which results in a better combination of the computational efficiency of the Cartesian grid and the flexi- bility of the gridless method in handling complex geometries. The clouds of points in the local gridless region are implemented in an anisotropic way according to the features of the thin boundary layer of the viscous flows over the airfoils, and the clouds of points at the vicinity of the interface between the grid and the gridless regions are also controlled by using an adaptive refinement technique during the generation of the unstructured Cartesian grid. An implementation of the resulting hybrid method is presented for solving two-dimensional compressible Navier-Stokes （NS） equations. The simulations of the viscous flows over a RAE2822 airfoil or a two-element airfoil are success- fully carried out, and the obtained results agree well with the available experimental data.     

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

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