非均匀矩形网格的局部网格细化LBM算法研究

传统的格子玻尔兹曼方法 (LBM),特别是基于均匀正方形网格的经典单松弛计算模型(SLBM),其算法鲁棒性和数值稳定性较差,限制了LBM的发展和应用.而网格细化策略可以有效缓解这一窘境,但是传统LBM中网格细化必然会导致计算效率骤降,计算设备要求攀高.为了解决这一问题,文章基于非均匀矩形网格结构,结合插值LBM算法的思路,在保证物面处和流动变化剧烈区域的局部网格细化以及计算精度的前提下,提出了25点拉格朗日插值LBM算法.以经典顶盖驱动方腔内流为算例,开展了包括不同网格分辨率和插值格式的对比分析研究.验证算例既包括了定常流动的数值模拟,也涉及了非定常周期性流动的求解.计算结果表明,相较于其他插值格式,拉格朗日插值格式表现优异;文章局部网格细化工作可以确保物面处及流动变化剧烈区域流动细节的捕捉;数值模拟算法可以为数值仿真提供可信的计算结果;同时大幅降低了总网格数量.因此很大程度上提升了计算效率;数值模拟方法鲁棒性较好,适用于包括定常和非定常流动的数值模拟.     

FVM与LBM分区耦合模拟圆柱绕流

构造了用于模拟远场边界下圆柱绕流的有限容积法(FVM)与格子Boltzmann方法(LBM)的分区耦合模型.模型中,靠近圆柱处采用多块网格的LBM,远离圆柱处采用FVM,并将计算结果同适体网格LBM以及多块网格LBM进行了比较.结果表明,耦合模型能在保证计算精度的前提下,显著提高计算效率.     

某型双立尾战斗机的网格生成及流场解

以某型双立尾战斗机为例，给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程。用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格，用求解椭圆型方程的方法生成块内空间网格。在所生成的双立尾战斗机分块结构网格中分区求解Ｅｕｌｅｒ方程，获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果     

薄膜结构流固耦合的CFD数值模拟研究

基于弱耦合分区求解策略,在CompaqVisualFortran6．5环境下搭建了薄膜结构三维流固耦合效应的CFD数值模拟平台。程序采用模块化编程思想,主要包含几何建模、流体分析、结构分析和数据交换四个模块。其中几何建模模块采用自行编制的膜结构找形分析程序,流体分析模块采用经过二次开发的计算流体力学软件FLUENT6．0,结构分析模块采用自行编制的膜结构动力分析程序MDLFX;在数据交换模块中,编制了基于薄板样条法的插值计算程序,以实现流固交界面上不同区域网格间的数据传递问题,编制了基于代数法和迭代法的动网格变形程序,以实现流固耦合运算中的动网格更新。基于该软件平台,对单向柔性屋盖和鞍形膜结构屋盖进行了流固耦合数值模拟,验证了方法的有效性。     

一种快速稳健的并行多块结构动网格方法

为解决传统网格处理方法不能满足复杂外形在大设计空间内进行优化时对网格质量的要求的问题,提出了一种并行多块动网格方法,该方法基于初始外形的多块结构网格,根据优化过程中个体外形与初始外形拓扑结构相近的特性,利用体样条插值方法来拟合多块结构网格各块顶点的位移,得到几何外形变化后的拓扑结构,再利用无限插值方法并行地移动初始外形多块结构网格的边、面和块内的网格点,进行光顺处理后得到变形后几何外形的空间网格;该方法在保证网格质量的同时,可以极大地提高网格生成效率,本文以某翼身组合体为例结果表明,该方法在大设计空间的复杂外形设计问题中具有很强的实用性。     

流固耦合分析的载荷映射开放式软件架构设计

基于开放式工程与科学计算软件平台SiPESC设计实现了流固耦合分析流场载荷映射软件架构。软件的核心问题是解决计算流体力学(CFD)网格模型与计算结构力学(CSD)网格模型交互界面网格不匹配情况下的流-固载荷映射问题。软件采用插值方法将流场分析得到的物面载荷转换为结构分析的载荷边界条件。软件基于SiPESC平台的微核心+插件的开放式可扩展软件框架进行设计,依托SiPESC.ENGDBS工程数据库管理系统实现大规模数据管理。设计实现的软件框架提供了算法的灵活扩展接口与管理机制,可动态扩展新的插值算法,满足流固耦合分析需要的数据管理与数据转换需求。在该软件框架下,已实现了多种插值算法,并完成验证算例与工程算例的载荷数据转换。算例表明软件功能具备良好的工程适用性,为进一步开发与应用奠定基础。     

一种新型高鲁棒性动网格技术及其应用

首先对四元数进行李代数空间指数映射,解决了多个四元数插值问题,并结合距离倒数插值方法实现网格边界扰动向空间网格的传播,建立了新型高鲁棒性的四元数变形网格技术. 针对该型动网格技术中由于大型矩阵运算量引起的运算效率低问题,同时利用四元数方法在动网格变形中具备与物面边界高阶一致性的特点,提出了分层次变形策略,避免了面向全流场网格节点的大型矩阵运算;进一步基于无限插值技术较强的逻辑保持能力,建立了面向结构网格分层混合变形方法. 充分利用多区域重叠、对接网格变形技术中隐含的并行性,基于对等式编程思想及MPI 库函数对动网格程序进行并行化编程,建立了高效高鲁棒性的变形网格技术. 以某型客机翼身组合体气动弹性分析为范例,研究了不同方法之间的计算效率以及鲁棒性,进一步将分层混合变形网格技术应用于某型支线客机全机型架外形设计与修正,验证了所建立的新型动网格技术的高效性与鲁棒性.      

基于非重叠Mortar方法的比例边界等几何分析

在比例边界等几何分析的建模剖分过程中会出现交界面网格非匹配现象,给计算分析带来困难。为了能够处理此类问题,提出了基于非重叠Mortar方法的比例边界等几何分析。该方法能在将全域分解为若干子域时,针对每个子域分别建模和剖分网格,交界面网格无需逐点匹配;采用交叉点修正的非均匀有理B-样条基函数构造Lagrange乘子空间,子域交界面连续条件可通过Mortar条件满足;并根据交界面连续条件进行自由度凝聚,得到对称正定的系数矩阵,可直接求解。分片试验、U形结构和半无限空间上的柔性基础等数值算例验证了本文方法的有效性,计算精度满足要求。     

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

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

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

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

Numerical study of radial basis function interpolation for data transfer across discontinuous mesh interfaces

Allowing discontinuous or non‐matching mesh spacing across zonal interfaces within a computational domain offers many advantages, particularly in terms of easing the mesh generation process, reduction of required mesh densities, and relative motion between mesh zones. This paper presents a numerical study of a universal method for interpolating solution data across such interfaces. The method utilises radial basis functions (RBFs) for n‐dimensional volume interpolation, and treats the available solution data points simply as arbitrary clouds of points, eliminating all connectivity requirements and making it applicable to a wide range of computational problems. Properties of the developed meshless interface interpolation are investigated using analytic functions, and three issues are considered: the achievable order of spatial accuracy of the RBF interpolation alone and comparison with a variable order polynomial; the effect of a combined RBF and polynomial interpolation; and the ability of the method to recover frequency content. RBF interpolation alone is shown to achieve fourth‐order to sixth‐order spatial accuracy in one and two dimensions, and in three dimensions, using a small number of data points, third‐order and above is achievable even for a 3 : 1 discontinuous cell spacing ratio, that is a 27 : 1 volume ratio, across the interface. Hence, it is inefficient to include polynomial terms, since improving on the RBF spatial accuracy results in a significant increase in the system size and deterioration in conditioning. It is also shown that only five points per wavelength are required to capture both frequency and amplitude content of periodic solutions to less than 0.01% error.Copyright © 2013 John Wiley & Sons, Ltd.     

Mesh adaptation for simulation of unsteady flow with moving immersed boundaries

In this work, an approach for performing mesh adaptation in the numerical simulation of two‐dimensional unsteady flow with moving immersed boundaries is presented. In each adaptation period, the mesh is refined in the regions where the solution evolves or the moving bodies pass and is unrefined in the regions where the phenomena or the bodies deviate. The flow field and the fluid–solid interface are recomputed on the adapted mesh. The adaptation indicator is defined according to the magnitude of the vorticity in the flow field. There is no lag between the adapted mesh and the computed solution, and the adaptation frequency can be controlled to reduce the errors due to the solution transferring between the old mesh and the new one. The preservation of conservation property is mandatory in long‐time scale simulations, so a P1‐conservative interpolation is used in the solution transferring. A nonboundary‐conforming method is employed to solve the flow equations. Therefore, the moving‐boundary flows can be simulated on a fixed mesh, and there is no need to update the mesh at each time step to follow the motion or the deformation of the solid boundary. To validate the present mesh adaptation method, we have simulated several unsteady flows over a circular cylinder stationary or with forced oscillation, a single self‐propelled swimming fish, and two fish swimming in the same or different directions. Copyright © 2012 John Wiley & Sons, Ltd.     

Shape morphing and topology optimization of fluid channels by explicit boundary tracking

An integrated shape morphing and topology optimization approach based on the deformable simplicial complex methodology is developed to address Stokes and Navier‐Stokes flow problems. The optimized geometry is interpreted by a set of piecewise linear curves embedded in a well‐formed triangular mesh, resulting in a physically well‐defined interface between fluid and impermeable regions. The shape evolution is realized by deforming the curves while maintaining a high‐quality mesh through adaption of the mesh near the structural boundary, rather than performing global remeshing. Topological changes are allowed through hole merging or splitting of islands. The finite element discretization used provides smooth and stable optimized boundaries for simple energy dissipation objectives. However, for more advanced problems, boundary oscillations are observed due to conflicts between the objective function and the minimum length scale imposed by the meshing algorithm. A surface regularization scheme is introduced to circumvent this issue, which is specifically tailored for the deformable simplicial complex approach. In contrast to other filter‐based regularization techniques, the scheme does not introduce additional control variables, and at the same time, it is based on a rigorous sensitivity analysis. Several numerical examples are presented to demonstrate the applicability of the approach.     

An immersed boundary solver for inviscid compressible flows

In this paper, a simple and efficient immersed boundary (IB) method is developed for the numerical simulation of inviscid compressible Euler equations. We propose a method based on coordinate transformation to calculate the unknowns of ghost points. In the present study, the body‐grid intercept points are used to build a complete bilinear (2‐D)/trilinear (3‐D) interpolation. A third‐order weighted essentially nonoscillation scheme with a new reference smoothness indicator is proposed to improve the accuracy at the extrema and discontinuity region. The dynamic blocked structured adaptive mesh is used to enhance the computational efficiency. The parallel computation with loading balance is applied to save the computational cost for 3‐D problems. Numerical tests show that the present method has second‐order overall spatial accuracy. The double Mach reflection test indicates that the present IB method gives almost identical solution as that of the boundary‐fitted method. The accuracy of the solver is further validated by subsonic and transonic flow past NACA2012 airfoil. Finally, the present IB method with adaptive mesh is validated by simulation of transonic flow past 3‐D ONERA M6 Wing. Global agreement with experimental and other numerical results are obtained.     

An improved Rhie–Chow interpolation scheme for the smoothed‐interface immersed boundary method

Rhie–Chow interpolation is a commonly used method in CFD calculations on a co‐located mesh in order to suppress non‐physical pressure oscillations arising from chequerboard effects. A fully parallelized smoothed‐interface immersed boundary method on a co‐located grid is described in this paper. We discuss the necessity of modifications to the original Rhie–Chow interpolation in order to deal with a locally refined mesh. Numerical simulation with the modified scheme of Choi shows that numerical dissipation due to Rhie–Chow interpolation introduces significant errors at the immersed boundary. To address this issue, we develop an improved Rhie–Chow interpolation scheme that is shown to increase the accuracy in resolving the flow near the immersed boundary. We compare our improved scheme with the modified scheme of Choi by parallel simulations of benchmark flows: (i) flow past a stationary cylinder; (ii) flow past an oscillating cylinder; and (iii) flow past a stationary elliptical cylinder, where Reynolds numbers are tested in the range 10–200. Our improved scheme is significantly more accurate and compares favourably with a staggered grid algorithm. We also develop a scheme to compute the boundary force for the direct‐forcing immersed boundary method efficiently. Copyright © 2016 John Wiley & Sons, Ltd.     

A local mesh refinement approach for large‐eddy simulations of turbulent flows

In this paper, a local mesh refinement (LMR) scheme on Cartesian grids for large‐eddy simulations is presented. The approach improves the calculation of ghost cell pressures and velocities and combines LMR with high‐order interpolation schemes at the LMR interface and throughout the rest of the computational domain to ensure smooth and accurate transition of variables between grids of different resolution. The approach is validated for turbulent channel flow and flow over a matrix of wall‐mounted cubes for which reliable numerical and experimental data are available. Comparisons of predicted first‐order and second‐order turbulence statistics with the validation data demonstrated a convincing agreement. Importantly, it is shown that mean streamwise velocities and fluctuating turbulence quantities transition smoothly across coarse‐to‐fine and fine‐to‐coarse interfaces. © 2016 The Authors International Journal for Numerical Methods in Fluids Published by John Wiley & Sons Ltd     

不可压缩N-S方程的稳定分步算法及耦合离散

在有限增量微积分(finite increment calculus, FIC)的理论框架下，通过引入一个附加变量，发展了压力稳定型分步算法，有效改善了经典 分步算法的压力稳定性，同时还避免了标准FIC方法中存在的空间高阶导数的计算. 为保证 数值方法同时具有较快的计算速度和较好的健壮性，发展了有限元与无网格的耦合空间离散 方法. 该方案可在网格发生扭曲的区域采用无网格法空间离散以保证求解的精度和稳定性， 而在网格质量较好的区域以及本质边界上保留使用有限元法空间离散以提高计算效率和便于 施加本质边界条件. 方腔流考题的数值模拟结果突出地显示了所发展的压力稳定型分步算 法比经典分步算法具有更好的压力稳定性，能够有效消除速度-压力插值空间违反LBB条件而 导致的压力场的虚假数值振荡. 平面Poisseuille流动和一个典型型腔充填过程的数值模拟 结果, 表明了发展的耦合离散方案相对于单一的有限元法和单一的无网格法在综合考虑计 算效率和算法健壮性方面的突出优点.     

一类面向高阶精度自适应流动计算的流场插值方法

网格自适应技术和高阶精度数值方法是提升计算流体力学复杂问题适应能力的有效技术途径. 将这两项技术结合需要解决一系列技术难题, 其中之一是高阶精度流场插值. 针对高阶精度自适应流动计算, 提出一类高精度流场插值方法, 实现将前一迭代步网格中流场数值解插值到当前迭代步网格中, 以延续前一迭代步中的计算状态. 为实现流场插值过程中物理量守恒, 该方法先计算新旧网格的重叠区域, 然后将物理量从重叠区域的旧网格中转移到新网格中. 为满足高阶精度要求, 先采用k-exact最小二乘方法对旧网格上的数值解进行重构, 获得描述物理量分布的高阶多项式, 随后采用高阶精度高斯数值积分实现物理量精确地转移到新网格单元上. 最后, 通过一个具有精确解的数值算例和一个高阶精度自适应流动计算算例验证了本文算法的有效性. 第一个算例结果表明当网格规模固定不变时, 插值精度阶数越高, 插值误差越小; 第二个算例显示本文方法可以有效缩短高精度自适应流动计算的迭代收敛时间.      

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

Influence of mesh deformation on the quality of large eddy simulations

The influence of mesh motion on the quality of large eddy simulation (LES) was studied in the present article. A three‐dimensional, turbulent pipe flow (Re_τ=360) was considered as a test case. Simulations with both stretching and static meshes were carried out in order to understand how mesh motion affects the turbulence statistics. The spatial filtering of static and moving mesh direct numerical simulation (DNS) data showed how an ideal LES would perform, while the comparison of DNS cases with static and moving meshes revealed that no significant numerical errors arise from the mesh motion when the simulation is fully resolved. The comparison of the filtered fields of the DNS with a moving mesh with the corresponding LES fields revealed different responses to mesh motion from different numerical approaches. A straightforward test was applied in order to verify that the moving mesh works consistently in LES: when the mesh is stretched in the streamwise direction, the moving mesh results should be in between the two extremal resolutions between which the mesh is stretched. Numerical investigations using four different LES approaches were carried out. In addition to the Smagorinsky model, three implicit LES approaches were used: linear interpolation (non‐dissipative), the Gamma limiter (dissipative), and the scale‐selective discretisation (slightly dissipative). The results indicate that while the Smagorinsky and the scale‐selective discretisation approaches produce results consistent with the resolution of the non‐static mesh, the implicit LES with linear interpolation or the Gamma scheme do not. Copyright © 2016 John Wiley & Sons, Ltd.