Solving Navier–Stokes equation for flow past cylinders using single-block structured and overset grids

Many complex fluid motions are driven by physical processes of instability, transition and turbulence dependent upon nonlinear mechanisms. Here, we solve the flow past cylinder(s) using single-block structured and overset grids by computing Navier–Stokes equation in two-dimensions. The suitability of a compact scheme in discretizing convection and diffusion terms are investigated first by looking at relevant numerical properties. Also, for the overset grid method, one of the methods is identified that shows the best results in minimizing interpolation error at sub-domain boundaries for an analytical test function. We provide extensive comparisons with experimental and other computational results for flow past a single cylinder, utilizing both single-block structured and Chimera or overset grids. Apart from showing instability of this flow calculated by these methods, we also compare the computed vorticity and velocity data using these two grids by employing the proper orthogonal decomposition (POD). We have analyzed and developed an overset grid method with compact scheme that does not need any filtering to control error. This has been ascertained by performing POD analysis. To show that the developed method is capable of handling complex geometries, we have computed flow past two cylinders in side-by-side arrangement. Results obtained capture the known flow characteristics for this arrangement well using relatively fewer number of grid points.     

并行重叠/变形混合网格生成技术及其应用

为了适用于柔性变形、相对运动等复杂动边界问题,建立了并行环境下重叠和变形相结合的动态混合网格生成技术.通过计算区域分解以及分布式并行实现了重叠和变形技术的结合,其中重叠网格采用了并行化的隐式装配方法,并发展了两种并行化查询策略.变形网格则采用了并行化的径向基函数（RBF）插值方法.并行化动态网格生成方法大幅提高了动态网格生成效率,有利于处理大规模的动边界问题.在此基础上,发展了基于变形/重叠动态混合网格的流动/运动/控制一体化数值模拟方法,进一步改进了耦合模拟软件平台——HyperFLOW.典型应用算例证明了该动态混合网格技术及一体化算法的实用性.      

包含动边界的非定常流场动网格数值模拟

发展建立了用于多体分离等包含动边界的非定常流场数值模拟方法,建立了笛卡尔坐标系下边界以任意速度运动的控制体上的流动控制方程,结合几何守恒律确定网格速度,发展了基于结构网格的动网格方法,网格移动量采用加权插值方法得到.通过数值模拟二维翼型及三维机翼的强迫振荡非定常流场,表明该方法可以数值模拟包含动边界的非定常流场.     

Provably stable overset grid methods for computational aeroacoustics

The simulation of sound generating flows in complex geometries requires accurate numerical methods that are non-dissipative and stable, and well-posed boundary conditions. A structured mesh approach is often desired for a higher-order discretization that better uses the provided grids, but at the expense of complex geometry capabilities relative to techniques for unstructured grids. One solution is to use an overset mesh-based discretization where locally structured meshes are globally assembled in an unstructured manner. This article discusses recent advancements in overset methods, also called Chimera methods, concerning boundary conditions, parallel methods for overset grid management, and stable and accurate interpolation between the grids. Several examples are given, some of which include moving grids.     

应用笛卡尔非结构切割网格进行外挂物投放的数值模拟

描述了一种新的网格生成技术,即笛卡尔非结构切割网格技术,采用叉树数据结构,完成了几种单段和多段翼型以及三维机翼的网格生成.应用中心有限体积法,对其绕流问题进行Euler方程数值模拟,并将计算结果与实验数据进行对比.在机翼绕流数值模拟的基础上,求解出机翼带外挂物的分离投放的流场计算问题.     

子母弹多体分离过程的非定常CFD/RBD数值仿真

为研究子母弹在抛撒时的干扰流场特性,选取多舱段子母弹为计算模型,基于课题组自主开发的非结构混合网格Reynolds平均Navier-Stokes方程求解程序HUNS3D,结合非结构嵌套网格技术,耦合六自由度刚体运动学方程,使用了改进的4阶Adams预估-校正法求解六自由度刚体运动方程.利用跨声速下典型外挂物分离作为验证算例,仿真结果与实验结果高度拟合,验证了求解器的精度.对锁定母弹自由度和释放母弹的自由度两种计算状态进行数值模拟.仿真结果表明:由于激波干涉作用,子弹与母弹之间有较强的耦合作用;释放母弹自由度后,子弹的气动力参数发生了较大变化.      

A multi-block ADI finite-volume method for incompressible Navier–Stokes equations in complex geometries

An efficient second-order accurate finite-volume method is developed for a solution of the incompressible Navier–Stokes equations on complex multi-block structured curvilinear grids. Unlike in the finite-volume or finite-difference-based alternating-direction-implicit (ADI) methods, where factorization of the coordinate transformed governing equations is performed along generalized coordinate directions, in the proposed method, the discretized Cartesian form Navier–Stokes equations are factored along curvilinear grid lines. The new ADI finite-volume method is also extended for simulations on multi-block structured curvilinear grids with which complex geometries can be efficiently resolved. The numerical method is first developed for an unsteady convection–diffusion equation, then is extended for the incompressible Navier–Stokes equations. The order of accuracy and stability characteristics of the present method are analyzed in simulations of an unsteady convection–diffusion problem, decaying vortices, flow in a lid-driven cavity, flow over a circular cylinder, and turbulent flow through a planar channel. Numerical solutions predicted by the proposed ADI finite-volume method are found to be in good agreement with experimental and other numerical data, while the solutions are obtained at much lower computational cost than those required by other iterative methods without factorization. For a simulation on a grid with O(10⁵) cells, the computational time required by the present ADI-based method for a solution of momentum equations is found to be less than 20% of that required by a method employing a biconjugate-gradient-stabilized scheme.     

导弹超声速尾退分离干扰特性试验

超/高超声速尾退分离在防热、保形、隐身、多次投放、回收等方面具有明显优势,有望成为高超声速飞行器载荷投放的优选方案。由此面临一类新的多体分离问题:超/高超声速尾退分离问题（aft super/hypersonic ejection separation,ASES）。超/高超声速尾退分离问题本质上是带空腔底部流动与多体分离构成的耦合问题,具有流场结构复杂、气动非定常非线性非对称效应显著的特点。针对超声速尾退分离问题,采用网格测力和轨迹捕获（captive trajectory system,CTS）风洞试验方法探索了尾退分离干扰流场的结构,发现可根据流场结构和舵效变化分为低速-亚声速无激波、高亚声速-跨声速弱激波、超声速激波和准自由流弱干扰4种典型干扰特征,揭示了尾流场影响后不同区域的全弹气动特性和舵效特性以及控制律、攻角、高度和Mach数对分离位移和姿态的影响规律。相关结论将有助于增强对尾退分离问题的认识,对尾退分离技术的工程实践具有参考价值。      

无人机投放分离特性仿真与优化

用于远程投放作战的无人机在与母机分离的过程中, 其姿态变化与运动轨迹直接影响母机的安全性和自身稳定性。基于重叠网格技术, 提出了一种流体-刚体运动耦合数值模拟方法, 对无人机动态投放分离过程进行仿真。采用WPFS标模进行计算, 发现仿真与试验结果吻合度高, 数值模拟方法得到有效验证。采用上述数值模拟方法实现无人机投放过程的仿真, 得到无人机在投放分离过程中位置与姿态随时间的变化规律。进一步基于仿真结果, 完成了无人机弹射方案的优化设计。      

Nested Cartesian grid method in incompressible viscous fluid flow

In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow–structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier–Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure.     

Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a synthetic porous medium using the LB method and on the quadtree grid structure. To this end, the LB method was used on nonuniform grids coupled with a technique for image reconstruction which resulted in the quadtree grids for simulation of fluid flow through porous media. Accuracy and efficiency of this algorithm is compared against the conventional LB method based on uniform grids. While the decrease in computational time in the proposed LB method on nonuniform grids is found to be significant regarding the size of the initial and reconstructed images, the same level of accuracy is obtained when compared with the conventional LB method on uniform grids.     

八叉树网格到非结构混合网格的转换

在数值模拟中, 非结构网格的优势是可以采用相同的数值格式统一处理任意复杂的计算区域, 但在网格生成过程中难度大, 并且不容易控制网格质量。树结构网格可以认为是介于结构网格和非结构网格之间的一种网格, 目前已经有相对成熟的方法快速在复杂区域内生成二维四叉树网格和三维八叉树网格。在实际应用中, 数值方法往往需要在连接协调的非结构网格上做离散, 树结构网格中不同尺寸的网格之间连接不是协调的, 在应用上会受到很多限制。文章实现了树结构网格到非结构混合网格的转换, 这种转换在二维情况下就是将四叉树网格转换为非结构三角形和四边形的混合网格, 三维情况下则将八叉树网格转换为非结构混合网格。这一转换过程的难点在于需要考虑数千种不同的八叉树单元, 并给出能实现连接协调的非结构混合网格划分。可以出现的网格单元包括六面体、三棱柱、金字塔和四面体这4种不同情况。通过特别的分类, 实现了程序的自动生成, 这种程序自动生成技术一方面可以避免人工编写大量程序时的失误, 另一方面也使得对数以千计的不同情况的处理成为可能。通过对几个简单网格的测试, 对网格数据转换方法做了初步的验证。      

物体穿越激波非定常流场的统一分区加强解

用数值方法求解非定常可压缩Navier Stokes方程,模拟了物体以超音速穿越与之同向运动激波及逆向运动激波的全过程.为确保数值模拟中的激波和物体的相对运动有足够的网格分辨率,使用了统一分区加强解算法,详细给出了追击穿越和碰撞穿越的流场结构和气动力,并且通过和单区结果的对比,说明统一分区加强解算法是一种处理复杂几何外形和提高网格分辨率的简单有效方法.     

非设计工况下动静叶相互干扰的非定常流动特性

本文对非设计工况下动静叶栅相互干扰的非定常流动特性进行数值研究。非设计工况主要考虑大攻角下的分离 流动,动静叶干扰的非定常流动的数值求解在N—S方程的基础上采用分区计算的方法来完成。数值分析结果揭示了在分 离流动本身的非定常特性与动静叶栅相互干扰的非定常特性的双重影响下的流场情况。     

A high-order cut-cell method for numerical simulation of hypersonic boundary-layer instability with surface roughness

Laminar-turbulent transition of hypersonic boundary layers can be affected significantly by the existence of surface roughness. Currently many important mechanisms of roughness-induced transition are not well understood. In recent years, direct numerical simulation (DNS) has been extensively applied for investigating instability and transition mechanisms of hypersonic boundary layers. Most of the past DNS studies, however, have been based on body-fitted grids for smooth surfaces without roughness. Due to complex geometry of embedded roughness, the use of body-fitted grids can be very difficult for flow with arbitrary surface roughness. In this paper, we present a new high-order cut-cell method to overcome the natural complexities in grid generation around arbitrary surface of roughness. The new method combines a non-uniform-grid finite-difference method for discrete grid points near the solid boundary and a shock-fitting method for the treatment of the bow shock. The non-uniform-grid finite-difference formulas are expressed in a general explicit form so that they can be applied to different multi-dimensional problems without any modification. The computational accuracy of new algorithms of up to O(h⁴) are tested on several one- and two-dimensional elliptic equations in irregular domains. In addition, the new method is applied to the simulation of the receptivity process of a Mach 5.92 flow over a flat plate under the combined effect of an isolated surface roughness element and surface blow and suction. A good agreement is found between the unsteady flow results and those obtained by a Linear Stability Theory (LST).     

High-order non-uniform grid schemes for numerical simulation of hypersonic boundary-layer stability and transition

The direct numerical simulation of receptivity, instability and transition of hypersonic boundary layers requires high-order accurate schemes because lower-order schemes do not have an adequate accuracy level to compute the large range of time and length scales in such flow fields. The main limiting factor in the application of high-order schemes to practical boundary-layer flow problems is the numerical instability of high-order boundary closure schemes on the wall. This paper presents a family of high-order non-uniform grid finite difference schemes with stable boundary closures for the direct numerical simulation of hypersonic boundary-layer transition. By using an appropriate grid stretching, and clustering grid points near the boundary, high-order schemes with stable boundary closures can be obtained. The order of the schemes ranges from first-order at the lowest, to the global spectral collocation method at the highest. The accuracy and stability of the new high-order numerical schemes is tested by numerical simulations of the linear wave equation and two-dimensional incompressible flat plate boundary layer flows. The high-order non-uniform-grid schemes (up to the 11th-order) are subsequently applied for the simulation of the receptivity of a hypersonic boundary layer to free stream disturbances over a blunt leading edge. The steady and unsteady results show that the new high-order schemes are stable and are able to produce high accuracy for computations of the nonlinear two-dimensional Navier–Stokes equations for the wall bounded supersonic flow.     

解三维扩散方程的积分内插法格式

研究三维扩散方程的数值模拟.在非正规六面体网格上,使用积分内插法建立扩散方程差分格式,涉及到27个相邻网格,适用于大变形网格上带间断系数的拟线性扩散方程的计算.叙述差分格式的建立,推导通量流和网格顶点温度的计算公式,给出了数值试验结果.     

Euler及N-S方程的二维四边形/叉树混合网格自适应算法

描述了一种基于直角叉树网格的Euler和N-S方程自适应算法.由于考虑了粘性的作用,提出并使用了四边形叉树混合网格的方法,在几何表面附近生成贴体的四边形网格,外流场使用直角叉树网格.采用中心有限体积法,对Euler及N-S方程进行数值求解,对N-S方程的计算中加入了B-L代数湍流模型.在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离等特性的捕捉和分辨能力.采用上述方法,数值分析了单段和多段翼型的绕流问题.     

圆柱绕流问题的初期流动数值模拟

本文继Zheng-huan Teng[1]中介绍的解Navier-stokes方程的椭圆涡团法,研究了一种新的变形涡团法,用以模拟不可压粘性流体绕圆柱的不定常流动。圆柱在静止流体中突然起动并做匀速直线运动。对整个流动区域构造完全Navier-Stokes方程的解并不容易,近十年出现很多数值研究,本文对算法有所推进。把本文的方法称作变形涡团法是因为圆柱边界附近的流体中用椭圆涡团,远离边界时用圆形涡团。计算圆柱绕流比平板绕流在满足附着条件上更为困难,本文分析了怎样在圆柱边界上给出适当的附着条件的数值方法。在算例中雷诺数分别取200、550、3000,得到了不定常边界层分离,二次涡等复杂的物理现象,这些数值结果与近年实验结果[2]是一致的。     

Algebraically explicit analytical solutions for the unsteady non-Newtonian swirling flow in an annular pipe

This paper presents two algebraically explicit analytical solutions for the incompressible unsteady rotational flow of Oldroyd-B type in an annular pipe. The first solution is derived with the common method of separation of variables. The second one is deduced with the method of separation of variables with addition developed in recent years. The first analytical solution is of clear physical meaning and both of them are fairly simple and valuable for the newly developing computational fluid dynamics. They can be used as the benchmark solutions to verify the applicability of the existing numerical computational methods and to inspire new differencing schemes, grid generation ways, etc. Moreover, a steady solution for the generalized second grade rheologic fluid flow is also presented. The correctness of these solutions can be easily proven by substituting them into the original governing equation.