An efficient large-scale mesh deformation method based on MPI/OpenMP hybrid parallel radial basis function interpolation

An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption. Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation, the present method employs different surface points sets and influence radius for each volume point movement, accompanied with efficient geometry searching strategy. The deformed surface points, also called Control Points (CPs), are stored in each processor. The displacement of spatial points is interpolated by using only 20–50 nearest control points, and the local influence radius is set to 5–20 times the maximum displacement of control points. To shorten the searching time for the nearest control point clouds, an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique. Besides, an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications. Three 3D cases, including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements, are used to test the present mesh deformation method. The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 45° and more than 80% parallel efficiency with 1024 MPI processors. In addition, the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.     

High-order compact finite volume methods on unstructured grids with adaptive mesh refinement for solving inviscid and viscous flows

In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The high-order variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4est based on the forest of octrees is adopted. Several two- and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.     

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.     

太阳帆板地面展开的气动阻力数值仿真

为研究太阳帆板地面展开试验中气动阻力的影响,基于动网格技术提出了一种时变空气阻力计算方法。该方法采用CFD商业软件进行仿真计算,分析了三帆板地面低速展开的气动阻力影响。分析结果表明,该帆板在运动过程中所受的阻力较小,与工程设计的预估结果相近。     

用动网格方法模拟导弹发射过程中的燃气射流流场

采用局部非结构网格弹性变形方法和网格再生方法相结合的动网格技术，计算导弹与发射简具有相对运动的燃气射流非定常流场。计算中将安置在发射简内的导弹作为运动实体，随着导弹运动，相应流场计算域边界发生变化。计算结果表明，采用动网格技术模拟导弹发射过程中燃气射流的非定常流场，具有较高的精度。     

三角网格模型的数据分块算法

对测量获得的三角网格模型进行分块。首先根据三角片中心点的高斯曲率和平均曲率的符号将三角片标示成8种类型。然后任取一个三角片作为生长点，沿着边界向外扩展“生长”，将具有相同类型的三角片组成一个连通的块。最后采用两种方法对块实施优化：根据平均曲率差值最小的原则，将小块融合到邻近块中．用两点间的最短路径来光滑边界。本方法中的三角网格和块的拓扑均可以是任意的，分块结果符合零件的特征构造。     

二维有限元网格全自动生成方法

主要研究二维有限元网格全自动生成方法。提出一种新的节点生成方法 ,通过实例验证该方法可靠。采用三角形形状因子评定三角形网格质量 ,可减少三角化的时间 ,优化三角网格质量。     

Simple and robust h-adaptive shock-capturing method for flux reconstruction framework

In this paper, a simple and robust shock-capturing method is developed for the Flux Reconstruction(FR) framework by combining the Adaptive Mesh Refinement(AMR) technique with the positivity-preserving property. The adaptive technique avoids the use of redundant meshes in smooth regions, while the positivity-preserving property makes the solver capable of providing numerical solutions with physical meaning. The compatibility of these two significant features relies on a novel limiter designed for...     

一种新型颤振激励系统的特性研究

在非结构网格的基础上，以欧拉方程为流场的主控方程，采用格心格式空间离散的有限体积法和四步Runge—Kutta显式时间推进的方法．对一种新型颤振激励系统的流场进行了数值研究。研究了该系统的工作原理和激励力特性。研究结果表明，该系统能很好地激励飞机结构。     

关于 CFD 高精度保真的数值模拟研究

本文由四部分组成。第一部分研究了利用高阶格式、网格、物理模型和高性能计算机求解 NS 方程时,计算数据的高精度保真性的要求以及高精度保真的指标。第二部分研究了利用大规模计算开展大涡模拟和直接数值模拟存在的问题。第三部分分析了在高雷诺数下 NS 方程计算方法和网格的关系。最后提出了建立计算数值验证、确认的新方法。