HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

Implicit time accurate simulation of unsteady flow

Implicit time integration was studied in the context of unsteady shock‐boundary layer interaction flow. With an explicit second‐order Runge–Kutta scheme, a reference solution to compare with the implicit second‐order Crank–Nicolson scheme was determined. The time step in the explicit scheme is restricted by both temporal accuracy as well as stability requirements, whereas in the A‐stable implicit scheme, the time step has to obey temporal resolution requirements and numerical convergence conditions. The non‐linear discrete equations for each time step are solved iteratively by adding a pseudo‐time derivative. The quasi‐Newton approach is adopted and the linear systems that arise are approximately solved with a symmetric block Gauss–Seidel solver. As a guiding principle for properly setting numerical time integration parameters that yield an efficient time accurate capturing of the solution, the global error caused by the temporal integration is compared with the error resulting from the spatial discretization. Focus is on the sensitivity of properties of the solution in relation to the time step. Numerical simulations show that the time step needed for acceptable accuracy can be considerably larger than the explicit stability time step; typical ratios range from 20 to 80. At large time steps, convergence problems that are closely related to a highly complex structure of the basins of attraction of the iterative method may occur. Copyright © 2001 John Wiley & Sons, Ltd.     

On the Realizability of Nonlinear Stress–Strain Relationships for Reynolds Stress Closures

A class of recently developed explicit algebraic stress models based on tensorially quadratic stress--strain relations [7] is subjected to a systematical realizability analysis. It is found that these models, which are of particular interest for their rigorous derivation from linear second-moment closure models, tend to produce inappropriate unrealizable results like negative turbulence energy components, even in simple shear flows. The cause of the defect is identified in conjunction with a set of realizability-furnishing constraints on the model coefficients. With the exception of the silent normal stress component in accelerated flow, the nature and rationale of the explicit algebraic stress model suggested by Gatski and Speziale [7] can be extended to maintain the realizability principle. Results obtained from the corresponding quasi-realizable quadratic eddy-viscosity model are reported in comparison with other nonlinear modelling practices.     

Hybrid grid generation for viscous flow analysis

Cartesian grid with cut‐cell method has drawn attention of CFD researchers owing to its simplicity. However, it suffers from the accuracy near the boundary of objects especially when applied to viscous flow analysis. Hybrid grid consisting of Cartesian grid in the background, body‐fitted layer near the object, and transition layer connecting the two is an interesting alternative. In this paper, we propose a robust method to generate hybrid grid in two‐dimensional (2D) and three‐dimensional (3D) space for viscous flow analysis. In the first step, body‐fitted layer made of quadrangles (in 2D) or prisms (in 3D) is created near the object's boundary by extruding front nodes with a speed function depending on the minimum normal curvature obtained by quadric surface fitting. To solve global interferences effectively, a level set method is used to find candidates of colliding cells. Then, axis‐aligned Cartesian grid (quadtree in 2D or octree in 3D) is filled in the rest of the domain. Finally, the gap between body‐fitted layer and Cartesian grid is connected by transition layer composed of triangles (in 2D) or tetrahedrons (in 3D). Mesh in transition layer is initially generated by constrained Delaunay triangulation from sampled points based on size function and is further optimized to provide smooth connection. Our approach to automatic hybrid grid generation has been tested with many models including complex geometry and multi‐body cases, showing robust results in reasonable time. Copyright © 2012 John Wiley & Sons, Ltd.     

The implementation of an aeronautical CFD flow code onto distributed memory parallel systems

The parallelization of an industrially important in‐house computational fluid dynamics (CFD) code for calculating the airflow over complex aircraft configurations using the Euler or Navier–Stokes equations is presented. The code discussed is the flow solver module of the SAUNA CFD suite. This suite uses a novel grid system that may include block‐structured hexahedral or pyramidal grids, unstructured tetrahedral grids or a hybrid combination of both. To assist in the rapid convergence to a solution, a number of convergence acceleration techniques are employed including implicit residual smoothing and a multigrid full approximation storage scheme (FAS). Key features of the parallelization approach are the use of domain decomposition and encapsulated message passing to enable the execution in parallel using a single programme multiple data (SPMD) paradigm. In the case where a hybrid grid is used, a unified grid partitioning scheme is employed to define the decomposition of the mesh. The parallel code has been tested using both structured and hybrid grids on a number of different distributed memory parallel systems and is now routinely used to perform industrial scale aeronautical simulations. Copyright © 2000 John Wiley & Sons, Ltd.     

An aspect ratio agglomeration multigrid for unstructured grids

A robust aspect ratio‐based agglomeration algorithm to generate high quality of coarse grids for unstructured and hybrid grids is proposed in this paper. The algorithm focuses on multigrid techniques for the numerical solution of Euler and Navier–Stokes equations, which conform to cell‐centered finite volume special discretization scheme, combines vertex‐based isotropic agglomeration and cell‐based directional agglomeration to yield large increases in convergence rates. Aspect ratio is used as fusing weight to capture the degree of cell convexity and give an indication of cell stretching. Agglomeration front queue is established to propagate inward from the boundaries, which stores isotropic vertex and also high‐stretched cell marked with different flag according to aspect ratio. We conduct the present method to solve Euler and Navier–Stokes equations on unstructured and hybrid grids and compare the results with single grid as well as MGridGen, which shows that the present method is efficient in reducing computational time for large‐scale system equations. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical investigation of a water flow in a rib-roughened channel by using Reynolds stress model

In this paper, water flow in a rib-roughened channel is investigated numerically by using Reynolds stress turbulence models (RSM) on a three-dimensional (3-D) domain. Computational results for mean streamwise velocity component and turbulent kinetic energy show good agreements with available experimental data. Five rib pitch-to-height ratios (p/h) of 1, 5, 10, 15 and 20 are analysed for six different Reynolds numbers (Re) of 3000, 7000, 12,000, 20,000, 40,000 and 65,000. Velocity vectors, streamlines and Reynolds stresses are showed for these ratios and Re numbers. Streamlines revealed that Reynolds numbers do not affect flowfield but play an important role in the Reynolds stresses.     

A data-driven memory model for solving turbulent flows with the pseudo-direct numerical simulation method

It is well known that the inherent three-dimensional and unsteady nature of turbulent flows is a stumbling block for all approaches aimed at resolving their spatial and temporal variability. The pseudo-direct numerical simulation (P-DNS) method for turbulent flows, proposed by the authors in a previous publication, focused on resolving the spatial variability, leaving the task of solving the temporal evolution to a highly simplified, parameter dependent model, to be adjusted in a case by case basis. Although some auspicious results were obtained, the applicability of P-DNS for problems of industrial interest required a more sophisticated method to deal with the temporal variability. In this sense, the present work proposes a new, parameter free, data-driven memory model for P-DNS. The model is based on the study of off-line DNS solutions of turbulent flows transitioning between statistically steady states in simple domains. The new P-DNS model is tested and successfully compared against existing methods in selected three-dimensional turbulent flows.     

A spectral quadrilateral multidomain penalty method model for high Reynolds number incompressible stratified flows

A two‐dimensional quadrilateral spectral multidomain penalty method (SMPM) model has been developed for the simulation of high Reynolds number incompressible stratified flows. The implementation of higher‐order quadrilateral subdomains renders this model a nontrivial extension of a one‐dimensional subdomain SMPM model built for the simulation of the same type of flows in vertically nonperiodic domains (Diamessis et al., J. Comp. Phys, 202 :298‐322, 2005). The nontrivial aspects of this extension consist of the implementation of subdomain corners, the penalty formulation of the pressure Poisson equation (PPE), and, most importantly, the treatment of specific challenges that arise in the iterative solution of the SMPM‐discretized PPE. The two primary challenges within the framework of the iterative solution of the PPE are its regularization to ensure the consistency of the associated linear system of equations and the design of an appropriate two‐level preconditioner. A qualitative and quantitative assessment of the accuracy, efficiency, and stability of the quadrilateral SMPM solver is provided through its application to the standard benchmarks of the Taylor vortex, lid‐driven cavity, and double shear layer. The capacity of the flow solver for the study of environmental stratified flow processes is shown through the simulation of long‐distance propagation of an internal solitary wave of depression in a manner that is free of numerical dispersion and dissipation. The methods and results presented in this paper make it a point of reference for future studies oriented toward the reliable application of the quadrilateral SMPM model to more complex environmental stratified flow process studies. Copyright © 2014 John Wiley & Sons, Ltd.     

A novel implicit numerical treatment for non‐linear turbulence models using high and low Reynolds number formulations

Non‐linear turbulence models can be seen as an improvement of the classical eddy‐viscosity concept due to their better capacity to simulate characteristics of important flows. However, application of non‐linear models demand robustness of the numerical method applied, requiring a stable discretization scheme for convergence of all variables involved. Usually, non‐linear terms are handled in an explicit manner leading to possible numerical instabilities. Thus, the present work shows the steps taken to adapt a general non‐linear constitutive equation using a new semi‐implicit numerical treatment for the non‐linear diffusion terms. The objective is to increase the degree of implicitness of the solution algorithm to enhance convergence characteristics. Flow over a backward‐facing step was computed using the control volume method applied to a boundary‐fitted coordinate system. The SIMPLE algorithm was used to relax the algebraic equations. Classical wall function and a low Reynolds number model were employed to describe the flow near the wall. The results showed that for certain combination of relaxation parameters, the semi‐implicit treatment proposed here was the sole successful treatment in order to achieve solution convergence. Also, application of the implicit method described here shows that the stability of the solution either increases (high Reynolds with non‐orthogonal mesh) or preserves the same (low Reynolds number applications). Additional advantages of the procedure proposed here lie in the possibility of testing different non‐linear expressions if one considers the enhanced robustness and stability obtained for the entire numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

CFD模拟方法的发展成就与展望

本文综述了计算流体力学(computational fluid dynamics, CFD),尤其是计算空气动力学的发展概况.从计算方法、网格技术、湍流模型、大涡模拟等方面分别总结了CFD所取得的成就,分析了当前存在的问题、困惑,展望了其发展趋势.在CFD计算方法中,主要介绍了中心格式、迎风格式、TVD格式、WENO格式、紧致格式以及间断Galerkin有限元方法,对不同方法的原理和特性进行了系统阐述.网格技术包括结构网格、非结构网格、混合网格以及重叠网格,重点讨论了重叠网格的若干关键技术.在湍流模型中,对目前的模型进行分类介绍,包括线性涡黏性模型、二阶矩模型、非线性模型等,还介绍了转捩模型、DES方法以及SAS方法等.在大涡模拟方法中,就其中若干相关的研究方向进行了探讨,包括滤波方法、亚格子模型、收敛标准、数值格式等. 文中还包含了作者在相关领域的若干研究成果.     

多层建筑屋顶风场的数值模拟

基于计算流体动力学软件Fluent17.2,以浙江温州大学某栋带女儿墙的多层试验教学楼为研究对象,采用RNG k-ε湍流模型对其进行数值模拟.通过对不同风向角下的数值模拟结果进行对比分析,探究屋顶风场与未受扰来流风场的区别,屋顶不同高度风场的差异,以及不同风向角来流情况下近屋面风场的分布特点、变化规律.结果表明:在离屋...     

A novel finite point method for flow simulation

A novel finite point method is developed to simulate flow problems. The mashes in the traditional numerical methods are supplanted by the distribution of points in the calculation domain. A local interpolation based on the properties of Taylor series expansion is used to construct an approximation for unknown functions and their derivatives. An upwind‐dominated scheme is proposed to efficiently handle the non‐linear convection. Comparison with the finite difference solutions for the two‐dimensional driven cavity flow and the experimental results for flow around a cylinder shows that the present method is capable of satisfactorily predicting the flow separation characteristic. The present algorithm is simple and flexible for complex geometric boundary. The influence of the point distribution on computation time and accuracy of results is included. Copyright © 2002 John Wiley & Sons, Ltd.     

Assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number

This paper focuses on the assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number. The Taylor–Green vortex at Re = 1600 is considered. The results are compared with those obtained using a pseudo‐spectral solver, converged on a 512³ grid and taken as the reference. The temporal evolution of the dissipation rate, visualisations of the vortical structures and the kinetic energy spectrum at the instant of maximal dissipation are compared to assess the results. At an effective resolution of 288³, the fourth‐order accurate discontinuous Galerkin method (DGM) solution (p = 3) is already very close to the pseudo‐spectral reference; the error on the dissipation rate is then essentially less than a percent, and the vorticity contours at times around the dissipation peak overlap everywhere. At a resolution of 384³, the solutions are indistinguishable. Then, an order convergence study is performed on the slightly under‐resolved grid (resolution of 192³). From the fourth order, the decrease of the error is no longer significant when going to a higher order. The fourth‐order DGM is also compared with an energy conserving fourth‐order finite difference method (FD4). The results show that, for the same number of DOF and the same order of accuracy, the errors of the DGM computation are significantly smaller. In particular, it takes 768³ DOF to converge the FD4 solution. Finally, the method is also successfully applied on unstructured high quality meshes. It is found that the dissipation rate captured is not significantly impacted by the element type. However, the element type impacts the energy spectrum in the large wavenumber range and thus the small vortical structures. In particular, at the same resolution, the results obtained using a tetrahedral mesh are much noisier than those obtained using a hexahedral mesh. Those obtained using a prismatic mesh are already much better, yet still slightly noisier. Copyright © 2013 John Wiley & Sons, Ltd.     

多重网格法求解轴流压气机内部流场

提出了一种能够与显式时间推进法、有限体积差分格式、粘性体积力法有效地配合使用的有效地配合使用的多重网格法来解轴流压气机内部流场。作者对美国NASA37^#跨音速轴流压气机转子内部流场的数值计算表明,多重网格法的计算结果与试验结果吻合, 同时该多重网格法的收敛速度比不用多重网格法加快了1．5倍。     

A correction for Reynolds-averaged-Navier–Stokes turbulence model under the effect of shock waves in hypersonic flows

Studies on the unphysical increase of turbulent quantities for RANS simulation induced by shock waves in hypersonic flows are carried out. Numerical experiments on the hypersonic flow over a blunt body reveal that the phenomenon of unphysical increase of turbulent quantities across the detached shock wave is induced by the strain-rate-based production terms of the k-$\omega$ and k-$\omega$ SST turbulence models, which leads to the over-prediction of aerothermal prediction. While this phenomenon does not occur for Spalart–Allmaras (S–A) turbulence model because of its vorticity-based production term. In order to eliminate this unphysical phenomenon, and to maintain the accuracy of the original models for boundary layer and separation flows, a new correction method for the k-$\omega$ and k-$\omega$ SST models is proposed: by comparing the orders of magnitude between the strain-rate-based and vorticity-based production terms, the vorticity-based production term is used near the shock waves, while the original strain-rate-based production term is still used in other regions. Finally, the correction method is applied to turbulence and transition flows over blunt bodies, and the numerical results show that the correction method effectively eliminates the unphysical increase of turbulent quantities across shock waves and improves the accuracy of aerothermal and transition onset location prediction.     

An efficient and accurate algebraic interface capturing method for unstructured grids in 2 and 3 dimensions: The THINC method with quadratic surface representation

We present in this paper an efficient and accurate volume of fluid (VOF) type scheme to compute moving interfaces on unstructured grids with arbitrary quadrilateral mesh elements in 2D and hexahedral elements in 3D. Being an extension of the multi‐dimensional tangent of hyperbola interface capturing (THINC) reconstruction proposed by the authors in Cartesian grid, an algebraic VOF scheme is devised for arbitrary quadrilateral and hexahedral elements. The interface is cell‐wisely approximated by a quadratic surface, which substantially improves the numerical accuracy. The same as the other THINC type schemes, the present method does not require the explicit geometric representation of the interface when computing numerical fluxes and thus is very computationally efficient and straightforward in implementation. The proposed scheme has been verified by benchmark tests, which reveal that this scheme is able to produce high‐quality numerical solutions of moving interfaces in unstructured grids and thus a practical method for interfacial multi‐phase flow simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

Deconvolution‐based nonlinear filtering for incompressible flows at moderately large Reynolds numbers

We consider a Leray model with a deconvolution‐based indicator function for the simulation of incompressible fluid flow at moderately large Reynolds number (in the range of a few thousands) with under‐resolved meshes. For the implementation of the model, we adopt a three‐step algorithm called evolve–filter–relax that requires (i) the solution of a Navier–Stokes problem, (ii) the solution of a Stokes‐like problem to filter the Navier–Stokes velocity field, and (iii) a final relaxation step. We take advantage of a reformulation of the evolve–filter–relax algorithm as an operator‐splitting method to analyze the impact of the filter on the final solution versus a direct simulation of the Navier–Stokes equations. In addition, we provide some direction for tuning the parameters involved in the model based on physical and numerical arguments. Our approach is validated against experimental data for fluid flow in an idealized medical device (consisting of a conical convergent, a narrow throat, and a sudden expansion, as recommended by the U.S. Food and Drug Administration). Numerical results are in good quantitative agreement with the measured axial components of the velocity and pressures for two different flow rates corresponding to turbulent regimes, even for meshes with a mesh size more than 40 times larger than the smallest turbulent scale. After several numerical experiments, we perform a preliminary sensitivity analysis of the computed solution to the parameters involved in the model. Copyright © 2015 John Wiley & Sons, Ltd.     

Adaptive mesh refinement and coarsening for aerodynamic flow simulations

A new mesh refinement technique for unstructured grids is discussed. The new technique presents the important advantage of maintaining the original grid skewness, thanks to the capability of handling hanging nodes. The paper also presents an interpretation of MacCormack's method in an unstructured grid context. Results for a transonic convergent–divergent nozzle, for a convergent nozzle with a supersonic entrance and for transonic flow over a NACA 0012 airfoil are presented and discussed. Copyright © 2004 John Wiley & Sons, Ltd.