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

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

非结构混合网格高超声速绕流与磁场干扰数值模拟

对均匀磁场干扰下的二维钝头体无粘高超声速流场进行了基于非结构混合网格的数值模拟.受磁流体力学方程组高度非线性的影响及考虑到数值模拟格式的精度,目前在此类流场的数值模拟中大多使用结构网格及有限差分方法,因而在三维复杂外形及复杂流场方面的研究受到限制.本文主要探索使用非结构网格(含混合网格)技术时的数值模拟方法.控制方程为耦合了Maxwell方程及无粘流体力学方程的磁流体力学方程组,数值离散格式采用Jameson有限体积格心格式,5步Runge-Kutta显式时间推进.计算模型为二维钝头体,初始磁场均匀分布.对不同磁感应强度影响下的高超声速流场进行了数值模拟,并与有限的资料进行了对比,得到了较符合的结果.     

THE INFLUENCE OF GLOBAL-DIRECTION STENCIL ON GRADIENT AND HIGH-ORDER DERIVATIVES RECONSTRUCTION OF UNSTRUCTURED FINITE VOLUME METHODS 1)

模板选择方式对非结构有限体积方法的计算准确性会产生显著影响. 在之前的工作中, 基于局部方向模板存在的问题, 我们探索了一种更加简单有效的全局方向模板选择方法, 并将其应用于二阶精度非结构有限体积求解器. 基于该方法找到的模板单元均沿着壁面法向与流向, 可有效捕捉流场变化, 反映流动的各向异性, 并且模板选择过程脱离了对网格拓扑的依赖, 避免了局部方向模板选择方法中复杂的阵面推进与方向判断过程, 克服了在大压缩比三角形网格上模板单元偏离壁面法向的现象, 同时在二阶精度求解器上得到了较高的计算精度与计算准确性. 为了进一步验证全局方向模板在高阶精度非结构有限体积方法中应用的可行性, 本文初步测试了该模板对变量梯度及高阶导数重构的影响. 经检验, 在不同类型的网格上, 采用全局方向模板得到的变量梯度与高阶导数误差明显低于局部方向模板, 同时也低于共点模板的计算误差. 此外, 在高斯积分点处由全局方向模板得到的变量点值与导数误差同样在三种模板中最低. 因此该模板选择方法在非结构有限体积梯度与高阶导数重构方面具有较好的数值表现, 具备在高阶精度非结构有限体积求解器中应用并推广的可行性.     

全局方向模板对非结构有限体积梯度与高阶导数重构的影响

模板选择方式对非结构有限体积方法的计算准确性会产生显著影响.在之前的工作中,基于局部方向模板存在的问题,我们探索了一种更加简单有效的全局方向模板选择方法,并将其应用于二阶精度非结构有限体积求解器.基于该方法找到的模板单元均沿着壁面法向与流向,可有效捕捉流场变化,反映流动的各向异性,并且模板选择过程脱离了对网格拓扑的依赖,避免了局部方向模板选择方法中复杂的阵面推进与方向判断过程,克服了在大压缩比三角形网格上模板单元偏离壁面法向的现象,同时在二阶精度求解器上得到了较高的计算精度与计算准确性.为了进一步验证全局方向模板在高阶精度非结构有限体积方法中应用的可行性,本文初步测试了该模板对变量梯度及高阶导数重构的影响.经检验,在不同类型的网格上,采用全局方向模板得到的变量梯度与高阶导数误差明显低于局部方向模板,同时也低于共点模板的计算误差.此外,在高斯积分点处由全局方向模板得到的变量点值与导数误差同样在三种模板中最低.因此该模板选择方法在非结构有限体积梯度与高阶导数重构方面具有较好的数值表现,具备在高阶精度非结构有限体积求解器中应用并推广的可行性.     

结合并行策略与距离减缩法的搜索算法及其在CAA中的应用

为提高流场与声场信息传递效率,建立了一种耦合MPI并行策略与改进距离减缩法的搜索算法。在完成点搜索后,采用一种适合于结构和非结构网格的形函数插值算法进行流场插值,实现了流场信息从流场网格到声场网格的快速传递。针对网格点搜索算法效率的验证,选用二维30P30N三段翼为研究对象,在不同的子进程数下进行对比分析。结果表明,相对于传统的距离减缩算法,本文提出的改进算法能有效提高CAA网格点的搜索效率,并且效率随着子进程数的增加而提高。在此基础上,将建立的流场与声场信息传递技术模块应用于CAA方法中,并对二维NACA0012翼型的后缘噪声问题进行计算分析,计算结果反映基于流场与声场信息快速传递算法的CAA方法能有效模拟宽频气动噪声问题。     

结合并行策略与距离减缩法的搜索算法及其在CAA中的应用

为提高流场与声场信息传递效率,建立了一种耦合MPI并行策略与改进距离减缩法的搜索算法。在完成点搜索后,采用一种适合于结构和非结构网格的形函数插值算法进行流场插值,实现了流场信息从流场网格到声场网格的快速传递。针对网格点搜索算法效率的验证,选用二维30P30N三段翼为研究对象,在不同的子进程数下进行对比分析。结果表明,相对于传统的距离减缩算法,本文提出的改进算法能有效提高CAA网格点的搜索效率,并且效率随着子进程数的增加而提高。在此基础上,将建立的流场与声场信息传递技术模块应用于CAA方法中,并对二维NACA0012翼型的后缘噪声问题进行计算分析,计算结果反映基于流场与声场信息快速传递算法的CAA方法能有效模拟宽频气动噪声问题。     

基于Navier-Stokes／Kirchhoff方程的悬停旋翼跨声速气动声学计算

研究了Kirchhoff积分面是否有盖有底,以及是否计及旋翼网格上的流场值,这两个因素对噪声预测结果的影响.发展了一种基于重叠网格的计算悬停旋翼远场噪声的数值方法.数值计算过程分为流场模拟和声场模拟两部分.悬停旋翼流场的数值模拟是在两个相互重叠的网格上进行的:在高质量的旋翼网格上求解Navier-Stokes方程,用于模拟旋翼附近的粘性流动和近场尾涡的捕捉;在远离粘性区域处布置符合悬停流场物理特征的圆柱形背景网格,控制方程为Euler方程,用于远场尾涡的捕捉.计算得到的流场信息插值到用于声场计算的Kirchhoff积分面上.观测点处的噪声可以认为是由这个完全包含桨叶的Kirchhoff积分面上的面元(声源)发声得到.远场声波的传播由Kirchhoff积分公式描述.计算结果表明:采用有盖有底的Kirchhoff积分面并且同时计及旋翼网格流场值时,计算得到的HSI噪声与实验值吻合最好.     

求解对流扩散方程的一种高效的有限体积法

考虑无结构三角网格上求解对流扩散方程的有限体积法.引入一种梯度函数的计算方法,将现有方法中计算解变量在网格单元中心和网格单元边界的梯度的两个独立过程改造成一个过程来完成,发展了一种求解对流扩散方程的高效的有限体积法.数值实验结果表明,该方法完全达到了已有方法同样的精度,而在计算速度上有明显的提高.     

非结构网格热化学非平衡辐射流场数值计算

针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法. 根据原子分子光谱理论逐条计算了100$\sim$1\,500\,nm间N, O, N$^{ + }$, O$^{ + }$的谱线以及N计算流体力学; 辐射; 热化学非平衡; 非结构网格; 有限体积法针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法.根据原子分子光谱理论逐条计算了100～1 500nm间N,O,N+,O+的谱线以及N2,O2,NO,N+2等分子的10个谱带,特别分析了NO的β'带,γ'带,δ带和ε带的辐射特性.采用耦合辐射的双温模型计算热化学非平衡流场,辐射源项通过直接求解辐射输运方程RTE(radiative transport equation)获得.在空间和方向上分别离散后,利用有限体积法求解不同方向上的辐射输运方程.计算得出了再入飞行器前驻点的辐射强度分布.采用该数值方法计算了MUSES-C模型在速度为11.6km/s时的绕流流场及前驻点处的辐射热流密度.并通过对比分析了热辐射对流场的影响.     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

A local adaptive grid procedure for incompressible flows with multigridding and equidistribution concepts

An adaptive grid solution procedure is developed for incompressible flow problems in which grid refinement based on an equidistribution law is performed in high-error-estimate regions that are flagged from a preliminary coarse grid solution. Solutions on the locally refined and equidistributed meshes are obtained using boundary conditions interpolated from the preliminary coarse grid solution, and solutions on both the refined and coarse grid regions are successively improved using a multigrid approach. For this purpose, suitable correction terms for the coarse grid equations are derived for all variables in the flagged regions. This procedure with Local Adaptation, Multigridding and Equidistribution (LAME) concepts is applied to various flow problems to demonstrate the accuracy improvements obtained using this method.     

Convergence acceleration of segregated algorithms using dynamic tuning additive correction multigrid strategy

A convergence acceleration method based on an additive correction multigrid–SIMPLEC (ACM‐S) algorithm with dynamic tuning of the relaxation factors is presented. In the ACM‐S method, the coarse grid velocity correction components obtained from the mass conservation (velocity potential) correction equation are included into the fine grid momentum equations before the coarse grid momentum correction equations are formed using the additive correction methodology. Therefore, the coupling between the momentum and mass conservation equations is obtained on the coarse grid, while maintaining the segregated structure of the single grid algorithm. This allows the use of the same solver (smoother) on the coarse grid. For turbulent flows with heat transfer, additional scalar equations are solved outside of the momentum–mass conservation equations loop. The convergence of the additional scalar equations is accelerated using a dynamic tuning of the relaxation factors. Both a relative error (RE) scheme and a local Reynolds/Peclet (ER/P) relaxation scheme methods are used. These methodologies are tested for laminar isothermal flows and turbulent flows with heat transfer over geometrically complex two‐ and three‐dimensional configurations. Savings up to 57% in CPU time are obtained for complex geometric domains representative of practical engineering problems. Copyright © 1999 John Wiley & Sons, Ltd.     

THE PARALLEL BLOCK ADAPTIVE MULTIGRID METHOD FOR THE IMPLICIT SOLUTION OF THE EULER EQUATIONS

A method capable of solving very fast and robust complex non-linear systems of equations is presented. The block adaptive multigrid (BAM) method combines mesh adaptive techniques with multigrid and domain decomposition methods. The overall method is based on the FAS multigrid, but instead of using global grids, locally enriched subgrids are also employed in regions where excessive solution errors are encountered. The final mesh is a composite grid with uniform rectangular subgrids of various mesh densities. The regions where finer grid resolution is necessary are detected using an estimation of the solution error by comparing solutions between grid levels. Furthermore, an alternative domain decomposition strategy has been developed to take advantage of parallel computing machines. The proposed method has been applied to an implicit upwind Euler code (EuFlex) for the solution of complex transonic flows around aerofoils. The efficiency and robustness of the BAM method are demonstrated for two popular inviscid test cases. Up to 19-fold acceleration with respect to the single-grid solution has been achieved, but a further twofold speed-up is possible on four-processor parallel computers.     

A robust multi‐grid pressure‐based algorithm for multi‐fluid flow at all speeds

This paper reports on the implementation and testing, within a full non‐linear multi‐grid environment, of a new pressure‐based algorithm for the prediction of multi‐fluid flow at all speeds. The algorithm is part of the mass conservation‐based algorithms (MCBA) group in which the pressure correction equation is derived from overall mass conservation. The performance of the new method is assessed by solving a series of two‐dimensional two‐fluid flow test problems varying from turbulent low Mach number to supersonic flows, and from very low to high fluid density ratios. Solutions are generated for several grid sizes using the single grid (SG), the prolongation grid (PG), and the full non‐linear multi‐grid (FMG) methods. The main outcomes of this study are: (i) a clear demonstration of the ability of the FMG method to tackle the added non‐linearity of multi‐fluid flows, which is manifested through the performance jump observed when using the non‐linear multi‐grid approach as compared to the SG and PG methods; (ii) the extension of the FMG method to predict turbulent multi‐fluid flows at all speeds. The convergence history plots and CPU‐times presented indicate that the FMG method is far more efficient than the PG method and accelerates the convergence rate over the SG method, for the problems solved and the grids used, by a factor reaching a value as high as 15. Copyright © 2003 John Wiley & Sons, Ltd.     

非均匀介质弹性波动方程的不规则网格有限差分方法

从弹性波动方程出发，提出了一种新的空间不规则网格有限差分方法，并用于求解非均匀各向异性介质中的弹性波正演问题。这种方法简单易行，对于复杂几何结构，例如低速层、套管井和非平面界面等，在较细的不规则网格上进行离散，计算时间和占用内存更少。与多重网格差分方法相比，该方法不需要粗、细网格之间的插值，所有网格差分计算在同一次空间迭代中完成。具有复杂几何交界面的模型计算，包括地下透镜体、套管井眼等，在确定弹性常数和密度后，用不规则网格的差分方法更易实现。该方法使用了Higdon吸收边界条件解决人工边界反射问题，引入了新的稳定性条件和网格频散条件，很好地消除了非物理散射波。理论模型的效值计算表明，该方法具有良好的稳定性和计算精度，在模拟非均匀介质弹性波传播时，比相同精度的规则网格有限差分方法计算速度更快。该方法易于推广到非结构网格和三维问题中。     

Surface triangulation for polygonal models based on CAD data

This paper presents an approach to the generation of unstructured surface meshes for Computer‐Aided Design (CAD) surface models represented as lists of polygons with minimum user interventions. Stereolithography (STL) data are adopted as an interface between a CAD system and the surface grid generator. STL files often include problems such as overlapping surfaces, gaps, and intersections. They have to be revised quickly and automatically before the surface models are used for the background grid of the surface grid generation. In this paper, we describe an automatic revision method for use with STL‐defined surface models. The advancing front method using geometric features is adopted directly on the modified STL surfaces. The capability of the method is demonstrated for several 3D surface models. Copyright © 2002 John Wiley & Sons, Ltd.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

Large eddy simulation with direct resolution of subgrid motion using a grid free vortex particle method

The paper presents a novel Large Eddy Simulation approach with a direct resolution of the subgrid motion of fine concentrated vortices. The method, proposed first by Kornev (2018), is based on combination of a grid based and the grid free computational vortex particle (VPM) methods. The large scale flow structures are simulated on the grid whereas the concentrated structures are modeled using VPM. Due to this combination the advantages of both methods are strengthened whereas the disadvantages are diminished. The procedure of the separation of small concentrated vortices from the large scale ones is based on LES filtering idea. The flow dynamics is governed by two coupled transport equations taking two-way interaction between large and fine structures into account. The fine structures are mapped back to the grid if their size grows due to diffusion. Algorithmic aspects specific for three dimensional flow simulations are discussed. Validity and advantages of the new approach are illustrated for a well tried benchmark test of the decaying homogeneous isotropic turbulence using the experimental data of Comte-Bellot and Corrsin (1971).     

Seismic propagation simulation in complex media with non-rectangular irregular-grid finite-difference

This paper presents a finite-difference (FD) method with spatially non-rectangular irregular grids to simulate the elastic wave propagation. Staggered irregular grid finite difference operators with a second-order time and spatial accuracy are used to approximate the velocity-stress elastic wave equations. This method is very simple and the cost of computing time is not much. Complicated geometries like curved thin layers, cased borehole and nonplanar interfaces may be treated with nonrectangular irregular grids in a more flexible way. Unlike the multi-grid scheme, this method requires no interpolation between the fine and coarse grids and all grids are computed at the same spatial iteration. Compared with the rectangular irregular grid FD, the spurious diffractions from “staircase” interfaces can easily be eliminated without using finer grids. Dispersion and stability conditions of the proposed method can be established in a similar form as for the rectangular irregular grid scheme. The Higdon‘s absorbing boundary condition is adopted to eliminate boundary reflections. Numerical simulations show that this method has satisfactory stability and accuracy in simulating wave propagation near rough solid-fluid interfaces. The computation costs are less than those using a regular grid and rectangular grid FD method.     

Computational properties of vertical grids for a nonhydrostatic anelastic model

Computational dispersion properties of all vertically staggered grids, which are presently available, are analysed in terms of frequency and group velocity components using the second-order centre difference scheme for a nonhydrostatic anelastic approximation system with a general method. The inertial-gravitational waves with a horizontal scale of a hundred-, ten- and one-kilometres are considered. The comparison analysis shows that the Charny-Phillips (CP) and Lorenz grids are suitable for waves at all abovementioned horizontal scales, while the Lorenz time staggered and Charny-Phillips time staggered grids are applicable only to waves with a horizontal scale less than 10 km. The unstaggered (N) grid is not suitable for simulating waves at any horizontal scale. In an idealised flow numerical test, the result on the CP grid has much less error than that on the N grid.