A RECONSTRUCTION METHOD FOR FINITE VOLUME FLOW FIELD SOLVING BASED ON INCREMENTAL RADIAL BASIS FUNCTIONS

由于流场参数重构中, 用于重构的基网格单元的物理参数波动量相对于均值较小, 径向基函数（RBF） 直接插值方法重构会产生较大的数值振荡, 论文提出了一种增量RBF 插值方法, 并用于有限体积的流场重构步, 明显改善了插值格式的收敛性和稳定性. 算例首先通过简单的一维模型说明该方法的有效性, 当目标函数波动量相对于均值为小量时, 增量RBF 插值能够抑制数值振荡; 进一步通过二维亚音速、跨音速定常无黏算例、静止圆柱绕流非定常算例以及超音速前台阶算例来说明该方法在典型流场数值求解中的通用性和有效性. 研究表明增量RBF 重构方法可陡峭地捕捉激波间断, 可有效改善流场求解的收敛性和稳定性, 数值耗散小, 计算效率高.     

基于拉氏变换的弹性动力学插值型重构核粒子法

插值型重构核粒子法的形函数具有离散点插值特性和不低于核函数的高阶光滑性,因而不仅可以直接施加本质边界条件,同时也保证了较高的计算精度．本文将弹性动力学方程作拉氏变换后,在变换域内用插值型重构核粒子法求解,最后再借助Durbin数值反演方法求得时间域的解．针对典型的弹性动力学问题,给出了插值型重构核粒子法的数值算例,并验证了本文方法的有效性．     

无网格法求解一类分段连续型延迟偏微分方程

无网格法是基于散点信息求解偏微分方程问题的数值方法,无网格法可减少或完全消除对网格的依赖,数值实施更加灵活.因此,考虑采用基于径向基函数的无网格插值法求解一类分段连续型延迟偏微分方程.首先,利用θ-加权有限差分法得到方程时间上的离散格式,利用基于径向基函数的无网格插值法近似空间导数,得到了全离散数值格式.采用的基函数是Multiquadric (MQ)径向基函数,MQ径向基函数在精度及稳定性等方面都优于其他径向基函数.其次,采用傅里叶分析方法对该方法进行稳定性分析,得到了该方法稳定的条件,且该条件只与时间步长有关.最后,通过数值算例验证了方法的收敛性和稳定性,从而说明了方法的有效性和适用性.     

基于拉格朗日插值的无网格直接配点法和稳定配点法

由于无网格法中大多数近似函数均为有理式,不具有Kronecker delta性质,因此难以精确地施加本质边界条件.边界误差较大容易导致整个求解域求解结果精度低,甚至引起数值不稳定现象.文章在无网格直接配点法和稳定配点法中引入拉格朗日插值函数作为形函数,构建了拉格朗日插值配点法(LICM)和拉格朗日插值稳定配点法(SLICM).由于拉格朗日插值具有Kronecker delta性质,可以像有限元法一样简单而精确地施加本质边界条件,提高这两种方法的数值求解精度.稳定配点法基于子域对强形式方程进行积分,可以满足高阶积分约束,即可以保证形函数在积分形式下也满足高阶一致性条件,实现精确积分.同时,进行子域积分还可以减少离散矩阵的条件数,从而提高算法的稳定性.进一步提高拉格朗日插值稳定配点法的精度和稳定性.通过数值算例验证这两种方法的精度、收敛性和稳定性,结果表明基于拉格朗日插值的配点法的精度优于基于重构核近似的配点法,拉格朗日插值稳定配点法的精度和稳定性均优于拉格朗日插值配点法.     

两点边值问题的小波配点法

根据多分辨分析,提出用任意连续的尺度函数构造区间上的插值基函数,形成以尺度函数为基础的求解两点边值问题的小波配点法.该方法中,尺度函数不受紧支撑、插值等性质的限制,计算复杂度小,数值解收敛性由多分辨分析理论保证.同时,给出边值条件的积分处理方法,能够方便地处理任意边界条件,当尺度函数不具有高阶导数时,该方法也能有效使用.数值算例表明,该方法是一个高效、高精度的算法.     

势问题的复变量重构核粒子法

在重构核粒子法的基础上,引入复变量,提出复变量重构核粒子法,在构造形函数时采用一维基函数建立二维问题的修正函数.应用于势问题,具有计算量小、精度高的优点.数值算例证明了方法的有效性.     

应用变域变分有限元法解平面叶栅反命题

反命题作为一种可变(未知)边界问题近年来得到了广泛的研究。本文给出了亚声速平面叶栅反命题计算的势函数变域变分有限元解法。变域变分通过把可变边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而使可变边界求解和流场分析可以完全耦合进行。本文针对亚声速平面叶栅的反命题，根据泛函的驻值必要条件，介绍了变域变分有限元方法的求解过程，最后给出了两个数值算例。这两个算例均采用四节点矩形单元的插值基函数，第一个算例用于检验程序的可靠性，第二个算例设计了一个给定叶面马赫数分布的叶型，并与试验结果进行对照。     

小波插值方法自适应数值求解时间进化微分方程

应用小波自相关函数的插值性质,得到任意给定函数的插值小波表达式,然后对其直接求导,可以得到函数导数的表达式。导数运算不再应用差分算法,扩展了小波方法在数值求解微分方程中的应用。由于小波基函数的有限支撑特点,小波方法可以有效地处理微分方程中解的局部突变问题。通过设定小波系数阀值,实现了求解过程的自适应。本文给出了两个算例,结果表明了算法的自适应特点及其向二维空间问题推广的有效性。     

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

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

双材料界面裂纹的改进广义有限差分法

采用数值方法进行断裂力学分析时,裂纹尖端奇异区域处理的好坏直接关系到最终断裂力学参数的求解精度。与传统均匀介质不同,复合材料界面裂纹渐近位移和应力场表现出剧烈的振荡特性,许多用于表征经典的平方根和负平方根物理场渐近性的传统方法也因此失效。论文提出了一种改进的广义有限差分法,该方法基于多元函数泰勒级数展开和移动最小二乘法的思想,将节点变量的各阶导数由相邻点集函数的加权线性累加来近似,具有无网格、无数值积分、数据准备简单、稀疏矩阵快速求解等优点。为提高该方法求解断裂力学问题的计算精度和数值稳定性,论文引入了裂尖奇异区域局部点簇的自动创建技术和一种基于局部点簇几何尺寸的矩阵正则化算法。数值算例表明,所提算法稳定,效率高,在不增加计算量的前提下,显著提高了裂尖近场力学参量和断裂力学参数的求解精度和数值稳定性。     

A Control Volume Radial Basis Function Techniques for the Numerical Simulation of Saturated Flows in Semi-confined Aquifer

In this work, a Control Volume Radial Basis Function technique (CV-RBF) is adapted to solve ground water flow in the saturated zone of the semi-confined aquifer. The CV-RBF method differs from classical CV methods in the way that the flux at the cell surfaces is computed. A local RBF interpolation of the field variable is performed at the centres of the cell being integrated and its neighbours. This interpolation is then used to reconstruct the solution and its gradient in the integration points which support the flux computation. In addition, it is required that such interpolation satisfies the governing equation in a certain number of points placed around the cell centres. In this way, the local interpolations become equivalent to local boundary-value problems. The CV-RBF method is combined with a local remeshing technique in order to track the phreatic surface, where the gradients required to satisfy the kinematic condition are computed by the same local RBF interpolations used for the flux computation. The proposed numerical approach is validated in a series of three-dimensional groundwater flow problems where the operations of recharging and extracting water from a semi-confined aquifer are modelled.     

Walsh函数有限体积法的多重网格特征研究

Walsh函数有限体积法(FVM-WBF)是一种能够在网格内部捕捉间断的新型数值方法. 持续增加Walsh基函数数目能够稳步提高FVM-WBF方法的求解分辨率, 但计算量暴发式增长和收敛速度下降的问题也会同步出现. 针对Walsh基函数数目增加而引起的计算效率问题, 本文分析了Walsh基函数及其系数所能影响的网格单元局部均值区域尺度, 发现其中隐含类似多重网格的尺度特征, 据此提出一种结合多重网格策略的FVM-WBF方法. 在定常流场计算中根据各级Walsh基函数影响尺度的不同, 对每级Walsh基函数设置满足其稳定性约束的时间步长, 在时间推进求解的过程中快速消除不同波长的数值误差, 实现多重网格的加速收敛效果. 选取NACA0012翼型和二维圆柱的定常无黏绕流问题作为算例, 对引入多重网格策略的FVM-WBF方法和不考虑多重网格策略的FVM-WBF方法进行对比测试. 数值结果证实: 新发展的FVM-WBF方法具备多重网格的关键特征, 在不增加任何特殊处理和计算量的情况下, 只需通过时间步长的调整, 就能够达到多重网格的加速效果, 显著提升计算效率.      

建筑风压场数值模拟中压力波动的平抑

在以同位网格为基础的简单流场压力计算中,通常采用动量插值方法来平抑流场中的压力波动现象;但是对于建筑风场等复杂的钝体绕流问题,由该平抑方法得到的收敛风压场仍可能存在小幅波动。为彻底解决同位网格格式下的压力波动,除采用动量插值方法外,本文提出了在压力校正方程的界面流速中添加压力梯度差值项的方法。算例分析表明,该方法计算得到的建筑风压场完全避免了压力波动现象,风压解与试验结果吻合良好。     

Mesh Deformation Method Based on Mean Value Coordinates Interpolation

Mesh deformation technique is widely used in many application fields, and has received a lot of attentions in recent years. This paper focuses on the methodology and algorithm of algebraic type mesh deformation for unstructured mesh in numerical discretization. To preserve mesh quality effectively, an algebraic approach for two and three dimensional unstructured mesh is developed based on mean value coordinates interpolation combined with node visibility analysis.The proposed approach firstly performs node visibility analysis to find out the visible boundary for each grid point to be moved, then evaluates the mean value coordinates of each grid point with respect to all vertices on its visible boundary. Thus the displacements of grid points can be calculated by interpolating the boundary movement by the mean value coordinates. Compared with other methods, the proposed method has good deformation capability and predictable computational cost, with no need to select parameters or functions. Applications of mesh deformation in different fields are presented to demonstrate the effectiveness of the proposed approach. The results of numerical experiments exhibit not only superior deformation capability of the method in traditional applications of fluid dynamic grid, but also great potential in modeling for large deformation analysis and inverse design problems.     

Developing a physical influence upwind scheme for pressure-based cell-centered finite volume methods

In the framework of a cell-centered finite volume method (FVM), the advection scheme plays the most important role in developing FVMs to solve complicated fluid flow problems for a wide range of Reynolds numbers. Advection schemes have been widely developed for FVMs employing pressure-velocity coupling methodology in the incompressible flow limit. In this regard, the physical influence upwind scheme (PIS) is developed for a cell-centered finite volume coupled solver (FVCS) using a pressure-weighted interpolation method for linking the pressure and velocity fields. The well-known exponential differencing scheme and skew upwind differencing scheme are also deployed in the current FVCS and their numerical results are presented. The accuracy and convergence of the present PIS are evaluated solving flow in a lid-driven square cavity, a lid-driven skewed cavity, and over a backward-facing step (BFS). The flow within the lid-driven square cavity is numerically solved at Reynolds numbers from 400 to 10 000 on a relatively coarse mesh with respect to other reported solutions. The lid-driven skewed cavity test case at Reynolds number of 1000 demonstrates the numerical performance of the present PIS on nonorthogonal grids. The flow over a BFS at Reynolds number of 800 is numerically solved to examine capabilities of current FVCS employing the current PIS in inlet-outlet flow computations. The numerical results obtained by the current PIS are in excellent agreement with those of benchmark solutions of corresponding test cases. Incorporating implicit role of pressure terms in a pressure-weighted interpolation method and development of PIS provides satisfactory solution convergence alongside the numerical accuracy for the current FVCS. A particular numerical verification is performed for the V velocity calculation within the BFS flow field, which confirms the reliability of present PIS.     

Non‐intrusive reduced‐order modelling of the Navier–Stokes equations based on RBF interpolation

We present a new non‐intrusive model reduction method for the Navier–Stokes equations. The method replaces the traditional approach of projecting the equations onto the reduced space with a radial basis function (RBF) multi‐dimensional interpolation. The main point of this method is to construct a number of multi‐dimensional interpolation functions using the RBF scatter multi‐dimensional interpolation method. The interpolation functions are used to calculate POD coefficients at each time step from POD coefficients at earlier time steps. The advantage of this method is that it does not require modifications to the source code (which would otherwise be very cumbersome), as it is independent of the governing equations of the system. Another advantage of this method is that it avoids the stability problem of POD/Galerkin. The novelty of this work lies in the application of RBF interpolation and POD to construct the reduced‐order model for the Navier–Stokes equations. Another novelty is the verification and validation of numerical examples (a lock exchange problem and a flow past a cylinder problem) using unstructured adaptive finite element ocean model. The results obtained show that CPU times are reduced by several orders of magnitude whilst the accuracy is maintained in comparison with the corresponding high‐fidelity models. Copyright © 2015 John Wiley & Sons, Ltd.     

A NUMERICAL STUDY OF THE FLOW OVER ELLIPSOIDAL OBJECTS INSIDE A CYLINDRICAL TUBE

A numerical scheme is developed to obtain the flow field around one, two and five ellipsoidal objects inside a cylindrical tube. The scheme uses the Galerkin finite element technique and the primitive variable(u–v–p) formulation. The two-dimensional incompressible Navier–Stokes equations are solved numerically by using the direct mixed interpolation method. A Picard iteration scheme is used for the solution of the resulting system of non-linear algebraic equations. The computer code is verified by checking with known analytical solutions for the flow past a sphere. Results for the shear stress distributions along the ellipsoids, forces and drag coefficients are obtained for different geometric ratios and Reynolds numbers. Some of the intermediate computational results on the velocity fields developed are also reported.     

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

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