基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

Chebyshev finite spectral method with extended moving grids

A Chebyshev finite spectral method on non-uniform meshes is proposed. An equidistribution scheme for two types of extended moving grids is used to generate grids. One type is designed to provide better resolution for the wave surface, and the other type is for highly variable gradients. The method has high-order accuracy because of the use of the Chebyshev polynomial as the basis function. The polynomial is used to interpolate the values between the two non-uniform meshes from a previous time step to the current time step. To attain high accuracy in the time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme is used. To avoid numerical oscillations caused by the dispersion term in the Korteweg-de Vries (KdV) equation, a numerical technique on non-uniform meshes is introduced. The proposed numerical scheme is validated by the applications to the Burgers equation (nonlinear convectiondiffusion problems) and the KdV equation (single solitary and 2-solitary wave problems), where analytical solutions are available for comparisons. Numerical results agree very well with the corresponding analytical solutions in all cases.     

Towards consistent hybrid overset mesh methods for rotorcraft CFD

The overset mesh method chimera is popular within the rotorcraft research community, because the use of multiple, non‐matching grids make the CFD simulations of bodies in relative motion much simpler. Consequently, the relative motion between the helicopter blades and fuselage can be accurately accounted for. In this paper, the method for treating overset grids within CFD codes is presented. It is compatible with multi‐block, structured‐grid solvers. The proposed method is based on hierarchy of overset, non‐matching grids, whose cells are automatically identified as computational or non‐computational and localised with respect to all grids they overlap with. The efficiency of the method relies on the hierarchical, multi‐step approach, for the overset mesh localisation and the use of a tree search. Because of the high efficiency of the algorithm, the search for overlapping cells can be carried out on‐the‐fly, during time‐marching of the unsteady, implicit CFD solver. In addition, the algorithm is suitable for parallel execution. The method has been demonstrated for several flows, ranging from simple aerofoils to rotor‐body interaction. The paper presents and demonstrates the method and shows that it has a low CPU overhead. It also highlights the limitations of the method and suggests remedies for improvement. Copyright © 2013 John Wiley & Sons, Ltd.     

Development of circular function‐based gas‐kinetic scheme (CGKS) on moving grids for unsteady flows through oscillating cascades

In this paper, the circular function‐based gas‐kinetic scheme (CGKS), which was originally developed for simulation of flows on stationary grids, is extended to solve moving boundary problems on moving grids. Particularly, the unsteady flows through oscillating cascades are our major interests. The main idea of the CGKS is to discretize the macroscopic equations by the finite volume method while the fluxes at the cell interface are evaluated by locally reconstructing the solution of the continuous Boltzmann Bhatnagar–Gross–Krook equation. The present solver is based on the fact that the modified Boltzmann equation, which is expressed in a moving frame of reference, can recover the corresponding macroscopic equations with Chapman–Enskog expansion analysis. Different from the original Maxwellian function‐based gas‐kinetic scheme, in improving the computational efficiency, a simple circular function is used to describe the equilibrium state of distribution function. Considering that the concerned cascade oscillating problems belong to cases that the motion of surface boundary is known a priori, the dynamic mesh method is suitable and is adopted in the present work. In achieving the mesh deformation with high quality and efficiency, a hybrid dynamic mesh method named radial basic functions‐transfinite interpolation is presented and applied for cascade geometries. For validation, several numerical test cases involving a wide range are investigated. Numerical results show that the developed CGKS on moving grids is well applied for cascade oscillating flows. And for some cases where nonlinear effects are strong, the solution accuracy could be effectively improved by using the present method. Copyright © 2017 John Wiley & Sons, Ltd.     

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

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

A local grid refinement method for three‐dimensional turbulent recirculating flows

A local grid refinement method is presented and applied to a three‐dimensional turbulent recirculating flow. It is based on the staggered grid arrangement. The computational domain is covered by block‐structured subgrids of different refinement levels. The exchange of information between the subgrids is fully conservative and all grids are treated implicitly. This allows for a simultaneous solution of one variable in all grids. All variables are stored in one‐dimensional arrays. The solver selected for the solution of the discretised finite difference equations is the preconditioned bi‐conjugate gradient (Bi‐CG) method. For the case examined (turbulent flow around a surface‐mounted cube), it was found that the latter method converges faster than the line solver. The locally refined mesh improved the accuracy of the pressure distribution on cube faces compared with a coarse mesh and yielded the same results as a fine single mesh, with a 62% gain in computer time. Copyright © 1999 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.     

不可压N-S方程高效算法及二维槽道湍流分析

构造了基于非等距网格的迎风紧致格式,并将其与三阶精度的Adams半隐方法相结合,构造了求解不可压N－S方程高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项,解决了边界处的残余散度问题;同时还利用快速Fourier变换将方程的隐式部分解耦,离散后的代数方程组利用追赶法求解,大大减少了计算量。通过对二维槽道流动的数值模拟,证实了所构造的数值方法具有精度高,稳定性好,能抑制混淆误差等优点,同时具有很高的计算效率,是进行壁湍流直接数值模拟的有效方法。在数值模拟的基础上对二维槽道流动进行了分析,得到了Reynolds数从6000到15000的二维流动饱和态解（所谓“二维槽道湍流”）;定性及定量结果均与他人的数值计算结果吻合十分理想。对流场进行了分析,指出了“二维湍流”与三维湍流统计特性的区别。     

Interpolated characteristic interface conditions for zonal grid refinement of high-order multi-block computations

When numerically calculating fluid flows around complex geometries using the high-order finite-difference method on structured grids, grid singularities are frequently observed, even if the grids are carefully generated. In multi-block computations with the generalised characteristic interface conditions (GCIC), decomposed subdomains (blocks) do not overlap but are connected by the inviscid characteristic relations. In the original theory of the GCIC, discontinuity of the metrics on the interface can be accurately treated; however, discontinuity of the grid lines on the interface is not allowed. This article proposes a theoretical extension to the GCIC by incorporating high-order interpolation methods; this extension is called GCIC with interpolation (GCIC + I). The basic concept and solution procedure of the multi-block computation with the GCIC + I are presented in detail, and two benchmark tests are conducted to validate the proposed theory.     

基于非结构/混合网格的高阶精度格式研究进展

尽管以二阶精度格式为基础的计算流体力学(CFD) 方法和软件已经在航空航天飞行器设计中发挥了重要的作用, 但是由于二阶精度格式的耗散和色散较大, 对于湍流、分离等多尺度流动现象的模拟, 现有成熟的CFD 软件仍难以给出满意的结果, 为此CFD 工作者发展了众多的高阶精度计算格式. 如果以适应的计算网格来分类, 一般可以分为基于结构网格的有限差分格式、基于非结构/混合网格的有限体积法和有限元方法,以及各种类型的混合方法. 由于非结构/混合网格具有良好的几何适应性, 基于非结构/混合网格的高阶精度格式近年来备受关注. 本文综述了近年来基于非结构/混合网格的高阶精度格式研究进展, 重点介绍了空间离散方法, 主要包括k-Exact 和ENO/WENO 等有限体积方法, 间断伽辽金(DG) 有限元方法, 有限谱体积(SV) 和有限谱差分(SD) 方法, 以及近来发展的各种DG/FV 混合算法和将各种方法统一在一个框架内的CPR (correctionprocedure via reconstruction) 方法等. 随后简要介绍了高阶精度格式应用于复杂外形流动数值模拟的一些需要关注的问题, 包括曲边界的处理方法、间断侦测和限制器、各种加速收敛技术等. 在综述过程中, 介绍了各种方法的优势与不足, 其间介绍了作者发展的基于"静动态混合重构" 的DG/FV 混合算法. 最后展望了基于非结构/混合网格的高阶精度格式的未来发展趋势及应用前景.     

Adaptive thermo-fluid moving boundary computations for interfacial dynamics

In this study,we present adaptive moving boundary computation technique with parallel implementation on a distributed memory multi-processor system for large scale thermo-fluid and interfacial flow computations.The solver utilizes Eulerian-Lagrangian method to track moving(Lagrangian) interfaces explicitly on the stationary(Eulerian) Cartesian grid where the flow fields are computed.We address the domain decomposition strategies of Eulerian-Lagrangian method by illustrating its intricate complexity of the computation involved on two different spaces interactively and consequently,and then propose a trade-off approach aiming for parallel scalability.Spatial domain decomposition is adopted for both Eulerian and Lagrangian domain due to easy load balancing and data locality for minimum communication between processors.In addition,parallel cell-based unstructured adaptive mesh refinement(AMR) technique is implemented for the flexible local refinement and even-distributed computational workload among processors.Selected cases are presented to highlight the computational capabilities,including Faraday type interfacial waves with capillary and gravitational forcing,flows around varied geometric configurations and induced by boundary conditions and/or body forces,and thermo-fluid dynamics with phase change.With the aid of the present techniques,large scale challenging moving boundary problems can be effectively addressed.     

Efficient high-order immersed interface methods for heat equations with interfaces

An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.     

双曲型守恒律的高阶、高分辨有限体积法

有限体积法是一种离散积分形式守恒律的数值方法。它可以吸收有限元法和有限差分法的一些重要思想与技巧。由于它可方便地利用多种类型的网格(结构网格和非结构网格),从而非常适用于处理复杂计算区域,目前已成为一种在计算流体力学中十分重要的方法。本文将针对二维双曲守衡律,对高精度、高分辨的有限体积法及其近年来的进展做一简要介绍。      

Pressure‐based unified solver for gas‐liquid two‐phase flows where compressible and incompressible flows coexist

We propose a pressure‐based unified solver for gas‐liquid two‐phase flows where compressible and incompressible flows coexist. Unlike the original thermo–Cubic Interpolated Propagation Combined Unified Procedure (CIP‐CUP) method proposed by Himeno et al (Transactions of the Japan Society of Mechanical Engineers, Series B, 2003), we split the advection term of the governing equations into a conservation part and into the rest. The splitting of advection term has two advantages. One is the high degree of freedom in choosing discretization schemes such as central‐difference schemes, upwind schemes, and Total Variation Diminishing (TVD) schemes. The other is the ease of implementation on unstructured grids. The advantages enable the analyses of various flows such as turbulent and supersonic ones in actual complicated boundaries. Therefore, the solver is useful for practical analyses. The solver was validated on the following test cases: subsonic single‐phase flows, incompressible single‐phase turbulent flows, and incompressible gas‐liquid two‐phase flows. With unstructured grids, we obtained the equivalent results as the ones with structured grids. After the validations, subsonic jet impinging on a water pool was calculated and compared with experimental results. It was confirmed that the calculated results were consistent with the experimental ones.     

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

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

DNS of premixed turbulent V-flame: coupling spectral and finite difference methods

To allow for a reliable examination of the interaction between velocity fluctuations, acoustics and combustion, a novel numerical procedure is discussed in which a spectral solution of the Navier–Stokes equations is directly associated to a high-order finite difference fully compressible DNS solver (sixth order PADE). Using this combination of high-order solvers with accurate boundary conditions, simulations have been performed where a turbulent premixed V-shape flame develops in grid turbulence. In the light of the DNS results, a sub-model for premixed turbulent combustion is analyzed. To cite this article: R. Hauguel et al., C. R. Mecanique 333 (2005).     

High-order ghost-point immersed boundary method for viscous compressible flows based on summation-by-parts operators

A high-order immersed boundary method is devised for the compressible Navier-Stokes equations by employing high-order summation-by-parts difference operators. The immersed boundaries are treated as sharp interfaces by enforcing the solid wall boundary conditions via flow variables at ghost points. Two different interpolation schemes are tested to compute values at the ghost points and a hybrid treatment is used. The first method provides the bilinearly interpolated flow variables at the image points of the corresponding ghost points and the second method applies the boundary condition at the immersed boundary by using the weighted least squares method with high-order polynomials. The approach is verified and validated for compressible flow past a circular cylinder at moderate Reynolds numbers. The tonal sound generated by vortex shedding from a circular cylinder is also investigated. In order to demonstrate the capability of the solver to handle complex geometries in practical cases, flow in a cross-section of a human upper airway is simulated.     

Simulation of three-dimensional incompressible flows in generalized curvilinear coordinates using a high-order compact finite-difference lattice Boltzmann method

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to provide an accurate 3-D incompressible flow solver. A high-order spectral-type low-pass compact filtering technique is applied to have a stable solution. All boundary conditions are implemented based on the solution of the governing equations in the 3-D generalized curvilinear coordinates. Numerical solutions of different 3-D benchmark and practical incompressible flow problems are performed to demonstrate the accuracy and performance of the solution methodology presented. Herein, the 2-D cylindrical Couette flow, the decay of a 3-D double shear wave, the cubic lid-driven cavity flow with nonuniform grids, the flow through a square duct with 90° bend and the flow past a sphere at different flow conditions are considered for validating the present computations. Numerical results obtained show the accuracy and robustness of the present solution methodology based on the implementation of the high-order compact finite-difference lattice Boltzman method in the generalized curvilinear coordinates for solving 3-D incompressible flows over practical and realistic geometries.     

A high-order Runge-Kutta discontinuous Galerkin method with a subcell limiter on adaptive unstructured grids for two-dimensional compressible inviscid flows

A robust, adaptive unstructured mesh refinement strategy for high-order Runge-Kutta discontinuous Galerkin method is proposed. The present work mainly focuses on accurate capturing of sharp gradient flow features like strong shocks in the simulations of two-dimensional inviscid compressible flows. A posteriori finite volume subcell limiter is employed in the shock-affected cells to control numerical spurious oscillations. An efficient cell-by-cell adaptive mesh refinement is implemented to increase the resolution of our simulations. This strategy enables to capture strong shocks without much numerical dissipation. A wide range of challenging test cases is considered to demonstrate the efficiency of the present adaptive numerical strategy for solving inviscid compressible flow problems having strong shocks.     

A PARALLEL BLOCK-STRUCTURED MULTIGRID METHOD FOR THE PREDICTION OF INCOMPRESSIBLE FLOWS

In this paper a parallel multigrid finite volume solver for the prediction of steady and unsteady flows in complex geometries is presented. For the handling of the complexity of the geometry and for the parallelization a unified approach connected with the concept of block-structured grids is employed. The parallel implementation is based on grid partitioning with automatic load balancing and follows the message-passing concept, ensuring a high degree of portability. A high numerical efficiency is obtained by a non-linear multigrid method with a pressure correction scheme as smoother. By a number of numerical experiments on various parallel computers the method is investigated with respect to its numerical and parallel efficiency. The results illustrate that the high performance of the underlying sequential multigrid algorithm can largely be retained in the parallel implementation and that the proposed method is well suited for solving complex flow problems on parallel computers with high efficiency.