可压流数值模拟中当地DFD方法的改进和应用

阐述了求解守恒型Euler 方程的当地DFD (Domain-Free Discretization) 方法的改进和应用。DFD 离散策略的核心,是解域内点上控制方程的离散形式可与解域外的一些点相关。通过边界附近解的近似形式,外部相关点上的流动变量值得到确定并强加相应的边界条件。与最初的当地DFD方法不同,在解的近似形式构建中,采用了CCST技术 (Curvature-Corrected Symmetry Technique),因此外部相关点上的密度和切向速度分别由等熵和等总焓关系确定。空间离散采用Galerkin 有限体积格式。最后,给出了固定和运动物体可压缩绕流的数值模拟结果,以验证改进的当地DFD方法的可靠性和数值解精度的提高。     

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

研究用于非结构动网格的弹簧近似方法，采用顶点弹簧描述，导出并讨论了弹簧倔强系数的取值。通过引入边界修正和扭转效应修正，对标准弹簧近似方法进行了改进，转动翼型算例的结果表明，改进后的方法大大提高了网格变形能力和网格质量，应用该动网格方法耦合求解基于（Arbitrary Lagrangian-Eulerian,ALE）描述的三维Euler方程，模拟了作俯抑振动的矩形机翼绕流，计算结果与实验数据及文献计算结果十分一致，作为多个自由刚与流体耦合运动问题的简单例证，耦合刚体动力学方程，模拟了激波与双立方体的相互作用，得到了非定常流场结构，研究表明，基于弹簧近似的非结构动网格与有限体积流式流场解算器相结合，是模拟包含运动边界的非定常流动问题的有效方法。     

基于反馈力浸入边界法模拟复杂动边界流动

浸入边界法是模拟流固耦合的重要数值方法之一。本文采用反馈力浸入边界方法,对旋转圆柱和水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟。其中,固体边界采用一系列离散的点近似代替,流体为不可压缩牛顿流体,使用笛卡尔自适应加密网格,利用有限差分法进行求解。固体对流场的作用通过构造适宜的反馈力函数实现。本文首先通过旋转圆柱绕流的计算结果同实验结果进行对比,吻合较好,验证了该计算方法的可靠性。然后针对水电站水力过渡过程中水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟,得到导叶动态绕流后的流场分布特性和涡结构的演化特性。     

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

计算含动边界非定常流动的无网格算法

在无网格算法中考虑了含动边界的流动问题,研究了可以计算处理包含一定位移及扭转动边界非定常流动的算法.创建了无网格算法的动点法则,并引入抗扭方法对弹簧方法进行改进来处理离散点运动,提高了方法的可用度及精度.发展了求解基于无网格的ALE方程组的算法,在点云离散的基础上采用曲面逼近计算空间导数及HLLC格式计算数值通量,运用四步龙格-库塔法进行时间推进.在跨、超音速条件下,计算模拟了典型翼型简谐振动流场,计算结果与实验结果及文献对比吻合,验证了该算法的正确性.     

前后双扑翼水平距离及相位差对其气动性能的影响

通过在动态网格上求解Navier-Stokes方程,对前后双扑翼的非定常粘性流场进行了数值模拟和气动干扰分析,考察了前后翼不同水平距离和不同相位差对其气动力和气动效率的影响。结果表明,扑动前翼和静止后翼间的气动干扰在各种不同水平距离下都有利于气动特性的改善,但气动干扰的作用随着前后翼水平距离的增大而减弱;前后双翼扑动的相位差是影响气动性能的重要参数,两翼间的气动干扰是否有利则与相位差和水平距离有直接联系。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界,避免生成复杂贴体网格及动网格,极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法,将简化附加体积力以源项形式嵌入动量方程迭代中,通过用户自定义函数对CFD软件FLUENT二次开发,实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证,并探讨了插值函数对计算精度的影响。研究表明,通过引入浸入边界模型,能够提高计算效率,并实现结构网格背景下复杂边界和动边界的高效建模。     

非结构重叠网格方法研究及应用

航空、航天和兵器技术等领域的研究中存在大量包含运动边界的流场。非结构重叠网格方法是一种高效的处理动边界问题的新方法。围绕相对运动的每个物体单独生成非结构网格,在网格重叠区域通过搜索和插值完成网格系之间的信息传递,提出了动态八叉树搜索算法,发展了绝对坐标系和相对坐标系相结合的流场求解方式,采用二阶精度Van Leer/Hanel格式和四阶Runge-Kutta法分别进行空间和时间离散,形成了一种新的非结构重叠网格算法。对三维Riemann问题的求解结果与精确解能很好吻合,证明了本文的重叠网格算法具有较好的时空离散精度和插值精度。对7.62mm步枪射击过程进行了数值模拟,描述了弹丸离开膛口后膛口流场的发展过程,与实验结果体现的发展过程较为吻合,验证了本文提出的非结构网格算法体系具有较好的计算性能,是研究含动边界复杂流场的一种有效手段。     

时刻追踪多介质界面运动的动网格方法

在对可压缩多介质流动的数值模拟中,定义介质界面为一种内部边界,由网格的边组成,界面边两侧对应两种不同介质中的网格。通过求解Riemann问题追踪介质界面上网格节点的运动,同时采用局部重构的动网格技术处理介质界面的大变形问题,并将介质界面定义为网格变形边界,以防止该边界上网格体积为负。运用HLLC格式求解ALE方程组得到整个多介质流场的数值解。最后从几个多介质流模型的计算结果可以看出,本文的动网格方法是可行的,而且可以时刻追踪介质界面的运动状态。     

一种改进的分块隐式数值方法

提出了一种改进的分块隐式数值方法，在贴体坐标和交错网格下以逆变速度分量和压力Ｕ，Ｖ，Ｗ，ｐ为基本求解变量，由此克服了原分块隐式数值方法求解复杂边界流动时的困难．９０°弯管流动数值计算初步表明，本文提出的方法合理、可行     

Improvement of mass conservation and reduction of spurious oscillations for the local DFD method

By introducing a mass source/sink term into the continuity equation, the mass-conservation property of the local domain-free discretisation (DFD) method is improved to reduce the spurious oscillations in the simulation of moving-boundary problems. The mass source/sink term is constructed by evaluating the mass flux through the solid part of the control volume split by the immersed boundary. Additionally, some special treatments are given for the multi-valued nodes associated with thin bodies. The introduction of source/sink term has also been extended to three-dimensional problems. Unlike the ghost-cell immersed boundary method, coupling the mass source/sink algorithm with the local DFD method is implemented on triangular and tetrahedral meshes. Compared to the hybrid reconstruction formulation, the properties of the original local DFD method can be preserved better. Numerical experiments for two- and three-dimensional moving-boundary problems show that the present mass source/sink treatment can reduce spurious oscillations effectively.     

A local domain‐free discretization method to simulate three‐dimensional compressible inviscid flows

In this paper, the domain‐free discretization method (DFD) is extended to simulate the three‐dimensional compressible inviscid flows governed by Euler equations. The discretization strategy of DFD is that the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior‐dependent points are updated at each time step by extrapolation along the wall normal direction in conjunction with the wall boundary conditions and the simplified momentum equation in the vicinity of the wall. Spatial discretization is achieved with the help of the finite element Galerkin approximation. The concept of ‘osculating plane’ is adopted, with which the local DFD can be easily implemented for the three‐dimensional case. Geometry‐adaptive tetrahedral mesh is employed for three‐dimensional calculations. Finally, we validate the DFD method for three‐dimensional compressible inviscid flow simulations by computing transonic flows over the ONERA M6 wing. Comparison with the reference experimental data and numerical results on boundary‐conforming grid was displayed and the results show that the present DFD results compare very well with the reference data. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of local DFD method to simulate unsteady flows around an oscillating circular cylinder

In this paper, the recently proposed local domain‐free discretization (DFD) method is applied to simulate incompressible flows around an oscillating circular cylinder. It is found that it is very easy for the local DFD method to handle such moving boundary flow problems. This is because it does not need to move the mesh, which is indeed needed in traditional methods. Numerical experiments show that the present numerical results agree very well with the available data in the literature, and that the local DFD method is an effective tool for the computation of moving boundary flow problems. Copyright © 2008 John Wiley & Sons, Ltd.     

A local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.     

Adaptive mesh refinement‐enhanced local DFD method and its application to solve Navier–Stokes equations

Recently, the domain‐free discretization (DFD) method was presented to efficiently solve problems with complex geometries without introducing the coordinate transformation. In order to exploit the high performance of the DFD method, in this paper, the local DFD method with the use of Cartesian mesh is presented, where the physical domain is covered by a Cartesian mesh and the local DFD method is applied for numerical discretization. In order to further improve the efficiency of the solver, the newly developed solution‐based adaptive mesh refinement (AMR) technique is also introduced. The proposed methods are then applied to the simulation of natural convection in concentric annuli between a square outer cylinder and a circular inner cylinder. Numerical experiments show that the present numerical results agree very well with available data in the literature, and AMR‐enhanced local DFD method is an effective tool for the computation of flow problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Extension of domain‐free discretization method to simulate compressible flows over fixed and moving bodies

This paper is the first endeavour to present the local domain‐free discretization (DFD) method for the solution of compressible Navier–Stokes/Euler equations in conservative form. The discretization strategy of DFD is that for any complex geometry, there is no need to introduce coordinate transformation and the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior dependent points are updated at each time step to impose the wall boundary condition by the approximate form of solution near the boundary. Some points inside the solution domain are constructed for the approximate form of solution, and the flow variables at constructed points are evaluated by the linear interpolation on triangles. The numerical schemes used in DFD are the finite element Galerkin method for spatial discretization and the dual‐time scheme for temporal discretization. Some numerical results of compressible flows over fixed and moving bodies are presented to validate the local DFD method. Copyright © 2006 John Wiley & Sons, Ltd.     

一种简单的局部离散速度统一气体动理论格式

统一气体动理论格式UGKS(Unified Gas-Kinetic Scheme)是一种适用于从连续流到自由分子流的全流域计算格式。在该格式中一般使用统一的离散速度空间。而在高速流动中,不同节点的分布函数往往差异很大。为了保证计算的精度,离散速度空间必须满足所有节点的需要,占用了大量的内存。采用局部的均匀离散速度空间,离散速度的范围随节点状态的变化而变化,从而降低了内存的需要,并通过引入背景网格避免了不同节点离散速度的插值。最后,通过两个一维算例对该方法进行了测试。测试结果显示,采用局部离散速度空间能够得到可靠的结果,并且在模拟高速流动时计算效率明显提高。     

An immersed boundary method for simulation of inviscid compressible flows

In this paper, an immersed boundary method for simulating inviscid compressible flows governed by Euler equations is presented. All the mesh points are classified as interior computed points, immersed boundary points (interior points closest to the solid boundary), and exterior points that are blanked out of computation. The flow variables at an immersed boundary point are determined via the approximate form of solution in the direction normal to the wall boundary. The normal velocity is evaluated by applying the no‐penetration boundary condition, and therefore, the influence of solid wall in the inviscid flow is taken into account. The pressure is computed with the local simplified momentum equation, and the density and the tangential velocity are evaluated by using the constant‐entropy relation and the constant‐total‐enthalpy relation, respectively. With a local coordinate system, the present method has been extended easily to the three‐dimensional case. The present work is the first endeavor to extend the idea of hybrid Cartesian/immersed boundary approach to compressible inviscid flows. The tedious task of handling multi‐valued points can be eliminated, and the overshoot resulting from the extrapolation for the evaluation of flow variables at exterior points can also be avoided. In order to validate the present method, inviscid compressible flows over fixed and moving bodies have been simulated. All the obtained numerical results show good agreement with available data in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

A mesh-free radial basis function–based semi-Lagrangian lattice Boltzmann method for incompressible flows

In this paper, a local radial basis function–based semi-Lagrangian lattice Boltzmann method (RBF-SL-LBM) is proposed. This is a mesh-free method that can be used for the simulation of incompressible flows. In this method, the collision step is performed locally, which is the same as in the standard LBM. In the meanwhile, the steaming step is solved in a semi-Lagrangian framework. The distribution functions at the departure points, which may be not the grid points in general, are computed by the local radial basis function interpolation. Several numerical tests are conducted to validate the present method, including the lid-driven cavity flow, the steady and unsteady flow past a circular cylinder, and the flow past an NACA0012 airfoil. The present results are in good agreement with those published in the previous literature, which demonstrates the capability of RBF-SL-LBM for the simulation of incompressible flows.