非结构网格变形方法研究进展

动网格技术中的非结构网格变形是计算流体力学和计算固体力学的关键技术之一.本文在总结现有非结构网格变形方法的基础上, 提出了一种网格变形方法的详细分类, 将现有方法归类为虚拟结构法、偏微分方程法和代数法. 本文综述了各类方法的最新研究进展, 分析并比较了各类方法的特性, 评述了当前网格变形研究的几个主要方向:复杂结构外形在不规则变形下的动网格生成、三维动网格生成、并行动网格生成和动节点技术. 最后简要地探讨了该领域的发展趋势.      

GMRES方法在含动边界流场中的应用

以基于格心的有限体积法为基础，空间二阶精度，采用4阶Runge—Kutta，GMRES隐式方法求解基于ALE形式的Euler方程，网格单元边界处守恒量通量的计算采用了Hanel方法，对NACA0012翼型绕流及运动圆球绕流等问题进行数值模拟，取得了较好的结果．GMRES方法克服了以往隐式方法大量耗费内存的弱点，达到了计算耗时短和占用内存少的统一．     

Two‐dimensional implicit time‐dependent calculations on adaptive unstructured meshes with time evolving boundaries

An implicit multigrid‐driven algorithm for two‐dimensional incompressible laminar viscous flows has been coupled with a solution adaptation method and a mesh movement method for boundary movement. Time‐dependent calculations are performed implicitly by regarding each time step as a steady‐state problem in pseudo‐time. The method of artificial compressibility is used to solve the flow equations. The solution mesh adaptation method performs local mesh refinement using an incremental Delaunay algorithm and mesh coarsening by means of edge collapse. Mesh movement is achieved by modeling the computational domain as an elastic solid and solving the equilibrium equations for the stress field. The solution adaptation method has been validated by comparison with experimental results and other computational results for low Reynolds number flow over a shedding circular cylinder. Preliminary validation of the mesh movement method has been demonstrated by a comparison with experimental results of an oscillating airfoil and with computational results for an oscillating cylinder. Copyright © 2005 John Wiley & Sons, Ltd.     

High-order curvilinear mesh generation technique based on an improved radius basic function approach

A high-order curvilinear hybrid mesh generation technique is developed for high-order numerical method (eg, discontinuous Galerkin method) applications to improve the accuracy for problems with curve boundary. The grid generation technique is based on an improved radius basic function (RBF) approach by which the straight-edge mesh is converted into high-order curve mesh. Firstly, an initial straight-edge mesh is prepared by traditional grid generation software. Then, high-order interpolation points are inserted into the mesh entities such as edges, faces, and cells according to the final demand of mesh order. To preserve the original geometry, the inserted points on solid wall are then projected onto the CAD model using an open source tool “Open Cascade.” Finally, other inserted points in the field near the solid wall are moved to appropriate positions by the improved RBF approach to avoid tangled cells. If we use the original RBF approach, then the inserted points on the edge and face entities normal to the solid boundary in the region of boundary layer will move to improper positions. To overcome this problem, a weighting based on the local grid aspect ratio between normal direction and tangential direction is introduced into the baseline RBF approach. Three typical configurations are tested to validate the mesh generator. Meanwhile, a third-order solution of subsonic flow over an analytical 3D body of revolution in the second International Workshop on High-Order CFD Methods is supplied by a discontinuous Galerkin solver. These numerical tests demonstrate the potential capability of present technique for high-order simulations of complex geometries.     

基于背景网格变形的动态网格移动方法

基于Delaunay背景网格插值技术的动态网格生成方法无需迭代计算,效率较高。但对复杂构形大幅运动的动边界问题,尤其当边界大幅转动时,背景网格极易交叉重叠。重新生成背景网格和重新定位网格节点信息不仅费时而且会导致网格质量的严重下降。本文提出改进的基于背景网格的动态网格变形方法,通过在初始Delaunay背景网格中添加辅助点,生成一层新的背景网格和新的映射关系;采用ball-vertex弹簧法驱动新背景网格的变形,进而牵动目标网格的变形。算例表明,本文提出的动态网格变形方法对所关心区域的网格具有良好保形性,边界可转动更大角度而不会出现网格交叉重叠问题,总体上提高了动态网格更新的效率和质量。     

Three‐dimensional automatic local remeshing for two or more hybrid meshes

In this paper, an automatic local remeshing method is proposed for merging two or more hybrid meshes consisting of prisms, hexahedra, pyramids and tetrahedra to reduce users' efforts and turnaround time for mesh generation as an extension of the local remeshing method for small devices by Ito et al. (AIAA J. 2009; 47 (5):1270–1276. DOI: 10.2514/1.40875 ). Two types of relationships between two volume meshes are considered: (1) a volume mesh containing the other smaller volume mesh inside (parent–child relationship) and (2) a volume mesh connecting to the other volume mesh (sibling relationship). The primary advantage of the proposed method over the existing ones is that it can handle the geometries in the two hybrid meshes intersecting each other. Two examples are shown to demonstrate the local remeshing method. Copyright © 2010 John Wiley & Sons, Ltd.     

数值模拟在钱塘江涌潮分析中的应用──Ⅱ.计算结果和分析

关于钱塘江涌潮分析中用数值模拟的计算结果和分析,对于一维情况,采用有精确解问题的计算、水跃的数值模拟和溃坝的数值模拟来检验数值方法的正确性,并将该方法应用于钱塘江涌潮的计算．对于二维情况,对存在理论解的直道水激波与有实验数据的湾道溃坝和非线性水波爬坡问题进行了计算,并将该方法应用于钱塘江涌潮的计算,数值计算结果表明了方法的有效性．     

数值模拟在钱塘江涌潮分析中的应用──Ⅱ.计算结果和分析

关于钱塘江涌潮分析中用数值模拟的计算结果和分析,对于一维情况,采用有精确解问题的计算、水跃的数值模拟和溃坝的数值模拟来检验数值方法的正确性,并将该方法应用于钱塘江涌潮的计算．对于二维情况,对存在理论解的直道水激波与有实验数据的湾道溃坝和非线性水波爬坡问题进行了计算,并将该方法应用于钱塘江涌潮的计算,数值计算结果表明了方法的有效性．     

A method for computing curved meshes via the linear elasticity analogy,application to fluid dynamics problems

We propose and analyze an algorithm for the robust construction of curved meshes in two and three dimensions. The meshes are made of curved simplexes. The algorithm starts from a mesh made of straight simplexes, and using a linear elasticity analogy applied on well‐chosen data, one can generate a curved mesh. Note that if the initial mesh has a boundary layer, this method allows to conserve it on the final mesh. This algorithm is used on several airfoils in two and three dimensions, including a turbulent M6 wing. Copyright © 2014 John Wiley & Sons, Ltd.     

基于二重网格的动网格松弛法改进

相比传统的弹簧法等方法,基于球松弛算法的动网格松弛法在复杂边界大变形条件下可以得到质量更高的边界网格以及更大的极限变形量,但该方法在时间效率上还有提升的空间.引入二重网格,采用动网格松弛法进行稀疏网格的网格变形,将边界位移传递到整个网格计算域;再利用二重网格映射,将稀疏网格位移映射到原有计算网格的节点上.算例表明,改进...     

三维非结构聚合多重网格法数值模拟研究

在三维非结构网格上应用聚合式多重网格技术来加速Euler方程的收敛过程．自行设计了一种高效率的网格聚合方法．采用四重三维非结构网格，在每一层网格上采用有限体积法进行计算．通过对M6翼型的数值求解验证了多重网格加速收敛的高效性．     

Numerical diffusion for flow‐aligned unstructured grids with application to estuarine modeling

The benefits of unstructured grids in hydrodynamic models are well understood but in many cases lead to greater numerical diffusion compared with methods available on structured grids. The flexible nature of unstructured grids, however, allows for the orientation of the grid to align locally with the dominant flow direction and thus decrease numerical diffusion. We investigate the relationship between grid alignment and diffusive errors in the context of scalar transport in a triangular, unstructured, 3‐D hydrodynamic code. Analytical results are presented for the 2‐D anisotropic numerical diffusion tensor and verified against idealized simulations. Results from two physically realistic estuarine simulations, differing only in grid alignment, show significant changes in gradients of salinity. Changes in scalar gradients are reflective of reduced numerical diffusion interacting with the complex 3‐D structure of the transporting flow. We also describe a method for utilizing flow fields from an unaligned grid to generate a flow‐aligned grid with minimal supervision. Copyright © 2013 John Wiley & Sons, Ltd.     

基于Voronoi cells的二维不规则自适应网格的生成及其应用

基于Voronoi cells的数据结构和算法，给出了一种二维不规则自适应网格的生成方法，用VisualC＋＋语言在微机上开发了Windows环境下网格自动生成的可视化软件。既可以得到Voronoi cells网格，也可得到相应的Delaunay triangles网格，网格生成的实例表明，本文方法所得到的网格非常适合于多尺度系统的流动问题的计算，具有较好的应用前景。     

Comparison of adjoint‐based and feature‐based grid adaptation for functional outputs

Adjoint‐based and feature‐based grid adaptive strategies are compared for their robustness and effectiveness in improving the accuracy of functional outputs such as lift and drag coefficients. The output‐based adjoint approach strives to improve the adjoint error estimates that relate the local residual errors to the global error in an output function via adjoint variables as weight functions. A conservative adaptive indicator that takes into account the residual errors in both the primal (flow) and dual (adjoint) solutions is implemented for the adjoint approach. The physics‐based feature approach strives to identify and resolve significant features of the flow to improve functional accuracy. Adaptive indicators that represent expansions and compressions in the flow direction and gradients normal to the flow direction are implemented for the feature approach. The adaptive approaches are compared for functional outputs of three‐dimensional arbitrary Mach number flow simulations on mixed‐element unstructured meshes. Grid adaptation is performed with h‐refinement and results are presented for inviscid, laminar and turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.     

一类多线程并行四面体网格优化算法

提出了一类多线程并行、组合了光滑化与局部重连操作的四面体网格优化算法。采用传统的基于拓扑的数据分解策略实现了并行光滑化算法,利用图染色算法将待光滑化的点分解成多个独立点集。同时提出将一类基于几何的数据分解策略应用于局部重连操作的并行化,在每个局部重连操作涉及的几何区域中定义一个特征点,然后沿希尔伯特(Hilbert)曲线对特征点进行排序,曲线的均匀分解对应局部重连操作在各线程的分配。这一分配策略的优点是使并行执行局部重连操作时重连区域相互干涉的情形极少出现。因此,当干涉情形出现时,可选择放弃产生干涉的操作,并行优化效率和效果并无明显的负面影响。最后,数值实验验证了本文算法的效率和有效性。     

A staggered control volume scheme for unstructured triangular grids

The purpose of this work is to introduce and validate a new staggered control volume method for the simulation of 2D/axisymmetric incompressible flows. The present study introduces a numerical procedure for solving the Navier–Stokes equations using the primitive variable formulation. The proposed method is an extension of the staggered grid methodology to unstructured triangular meshes for a control volume approach which features ease of handling of irregularly shaped domains. Two alternative elements are studied: transported scalars are stored either at the sides of an element or at its vertices, while the pressure is always stored at the centre of an element. Two interpolation functions were investigated for the integration of the momentum equations: a skewed mass-weighted upwind function and a flow-oriented exponential shape function. The momentum equations are solved over the covolume of a side or of a vertex and the pressure–velocity coupling makes use of a localized linear reconstruction of the discontinuous pressure field surrounding an element in order to obtain the pressure gradient terms. The pressure equation is obtained through a discretization of the continuity equation which uses the triangular element itself as the control volume. The method is applied to the simulation of the following test cases: backward-facing step flow, flow over a two-dimensional obstacle and flow in a pipe with sudden contraction of cross-sectional area. All numerical investigations are compared with experimental data from the literature. A grid convergence and error analysis study is also carried out for flow in a driven cavity. Results compared favourably with experimental data and so the new control volume scheme is deemed well suited for the prediction of incompressible flows in complex geometries. © 1997 John Wiley & Sons, Ltd.     

Adaptive meshing techniques for viscous flow calculations on mixed element unstructured meshes

An adaptive refinement strategy based on hierarchical element subdivision is formulated and implemented for meshes containing arbitrary mixtures of tetrahedra, hexahedra, prisms, and pyramids. Special attention is given to keeping memory overheads as low as possible. This procedure is coupled with an algebraic multigrid flow solver, which operates on mixed element meshes. Inviscid flows, as well as viscous flows, are computed on adaptively refined tetrahedral, hexahedral, and hybrid meshes. The efficiency of the method is demonstrated by generating an adapted hexahedral mesh containing 3 million vertices on a relatively inexpensive workstation. Copyright © 2000 John Wiley & Sons, Ltd.     

An adaptive Cartesian grid generation method for ‘Dirty’ geometry

Traditional structured and unstructured grid generation methods need a ‘water‐tight’ boundary surface grid to start. Therefore, these methods are named boundary to interior (B2I) approaches. Although these methods have achieved great success in fluid flow simulations, the grid generation process can still be very time consuming if ‘non‐water‐tight’ geometries are given. Significant user time can be taken to repair or clean a ‘dirty’ geometry with cracks, overlaps or invalid manifolds before grid generation can take place. In this paper, we advocate a different approach in grid generation, namely the interior to boundary (I2B) approach. With an I2B approach, the computational grid is first generated inside the computational domain. Then this grid is intelligently ‘connected’ to the boundary, and the boundary grid is a result of this ‘connection’. A significant advantage of the I2B approach is that ‘dirty’ geometries can be handled without cleaning or repairing, dramatically reducing grid generation time. An I2B adaptive Cartesian grid generation method is developed in this paper to handle ‘dirty’ geometries without geometry repair. Comparing with a B2I approach, the grid generation time with the I2B approach for a complex automotive engine can be reduced by three orders of magnitude. Copyright © 2002 John Wiley & Sons, Ltd.     

非结构动网格在可动边界问题中的应用研究

本文对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法并通过耦合求解基于(ArbitraryLagrangian-Eulerian ALE)描述的Euler方程,模拟了谐和振动NACA0012翼型及M6机翼的跨音速绕流,计算结果与参考文献提供的结果及实验结果吻合良好。