Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. Copyright © 2004 John Wiley & Sons, Ltd.     

Application of the finite point method to high‐Reynolds number compressible flow problems

In this work, the finite point method is applied to the solution of high‐Reynolds compressible viscous flows. The aim is to explore this important field of applications focusing on two main aspects: the easiness and automation of the meshless discretization of viscous layers and the construction of a robust numerical approximation in the highly stretched clouds of points resulting in such domain areas. The flow solution scheme adopts an upwind‐biased scheme to solve the averaged Navier–Stokes equations in conjunction with an algebraic turbulence model. The numerical applications presented involve different attached boundary layer flows and are intended to show the performance of the numerical technique. The results obtained are satisfactory and indicative of the possibilities to extend the present meshless technique to more complex flow problems. Copyright © 2014 John Wiley & Sons, Ltd.     

Evaluation of RBFs for volume data interpolation in CFD

A greedy method for choosing an optimum reduced set of control points is integrated with RBF interpolation and evaluated for the purpose of interpolating large‐volume data sets in CFD. Given a function defined at a set of points, the greedy method selects a small subset of these points that is sufficient to keep the interpolation error at all the remaining points below a chosen bound. This is equivalent to a type of data compression and would have useful storage, post‐processing, and computational applications in CFD. To test the method in terms of both the point selection scheme and the suitability of reduced control point volume interpolation, a trial application of the interpolation to velocity fields in CFD volume meshes is considered. To optimise the point selection process, and attempt to be able to capture multiple length scales, a variable support radius formulation has also been included. Structured and unstructured mesh cases are considered for aerofoils, a wing case and a wing‐body case. For smooth volume functions, the method is shown to work well, producing accurate velocity interpolations using a very small number of the cells in the mesh. For general complex fields including large gradients, the method is still shown to be effective, although large gradients require more interpolation points to be used.Copyright © 2013 John Wiley & Sons, Ltd.     

An efficient meshfree point collocation moving least squares method to solve the interface problems with nonhomogeneous jump conditions

We are going to study a simple and effective method for the numerical solution of the closed interface boundary value problem with both discontinuities in the solution and its derivatives. It uses a strong‐form meshfree method based on the moving least squares (MLS) approximation. In this method, for the solution of elliptic equation, the second‐order derivatives of the shape functions are needed in constructing the global stiffness matrix. It is well‐known that the calculation of full derivatives of the MLS approximation, especially in high dimensions, is quite costly. In the current work, we apply the diffuse derivatives using an efficient technique. In this technique, we calculate the higher‐order derivatives using the approximation of lower‐order derivatives, instead of calculating directly derivatives. This technique can improve the accuracy of meshfree point collocation method for interface problems with nonhomogeneous jump conditions and can efficiently estimate diffuse derivatives of second‐ and higher‐orders using only linear basis functions. To introduce the appropriate discontinuous shape functions in the vicinity of interface, we choose the visibility criterion method that modifies the support of weight function in MLS approximation and leads to an efficient computational procedure for the solution of closed interface problems. The proposed method is applied for elliptic and biharmonic interface problems. For the biharmonic equation, we use a mixed scheme, which replaces this equation by a coupled elliptic system. Also the application of the present method to elasticity equation with discontinuities in the coefficients across a closed interface has been provided. Representative numerical examples demonstrate the accuracy and robustness of the proposed methodology for the closed interface problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1031–1053, 2015     

薄板弯曲分析的节点积分高阶无网格法

采用高阶无网格法求解薄板弯曲问题,在已发展的线性曲率光顺方案的基础上,通过引入泰勒展开技术,建立了能够精确再现纯弯曲和线性弯曲模式的节点积分方法。与之相比,目前无网格薄板分析主要采用的节点积分方法仅能精确再现纯弯曲模式。数值结果表明,本文方法可精确通过纯弯曲和线性弯曲试验,且能得到光滑、无振荡的弯矩场。与标准的高斯积分方法和目前已存在的节点积分方法相比,本文方法在计算精度、效率以及弯矩分布等方面均展现出显著优势。     

Unified semi‐analytical wall boundary conditions for inviscid,laminar or turbulent flows in the meshless SPH method

Wall boundary conditions in smoothed particle hydrodynamics (SPH) is a key issue to perform accurate simulations. We propose here a new approach based on a renormalising factor for writing all boundary terms. This factor depends on the local shape of a wall and on the position of a particle relative to the wall, which is described by segments (in two‐dimensions), instead of the cumbersome fictitious or ghost particles used in most existing SPH models. By solving a dynamic equation for the renormalising factor, we significantly improve traditional wall treatment in SPH, for pressure forces, wall friction and turbulent conditions. The new model is demonstrated for cases including hydrostatic conditions for still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. The latter case is also compared with a finite volume and volume‐of‐fluid scheme. The performance of the model for a two‐dimensional laminar flow in a channel is demonstrated where the profiles of velocity are in agreement with the theoretical ones, demonstrating that the derived wall shear stress balances the pressure gradient. Finally, the performance of the model is demonstrated for flow in a schematic fish pass where both the velocity field and turbulent viscosity fields are satisfactorily reproduced compared with mesh‐based codes. Copyright © 2012 John Wiley & Sons, Ltd.     

On the use of characteristic‐based split meshfree method for solving flow problems

This study presents characteristic‐based split (CBS) algorithm in the meshfree context. This algorithm is the extension of general CBS method which was initially introduced in finite element framework. In this work, the general equations of flow have been represented in the meshfree context. A new finite element and MFree code is developed for solving flow problems. This computational code is capable of solving both time‐dependent and steady‐state flow problems. Numerical simulation of some known benchmark flow problems has been studied. Computational results of MFree method have been compared to those of finite element method. The results obtained have been verified by known numerical, analytical and experimental data in the literature. A number of shape functions are used for field variable interpolation. The performance of each interpolation method is discussed. It is concluded that the MFree method is more accurate than FEM if the same numbers of nodes are used for each solver. Meshfree CBS algorithm is completely stable even at high Reynolds numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

The characteristic‐based split (CBS) meshfree method for free surface flow problems in ALE formulation

A semi‐implicit characteristic‐based split (CBS) meshfree algorithm in the arbitrary Lagrangian Eulerian (ALE) framework is proposed for the numerical solution of incompressible free surface flow problem in the paper. The algorithm is the extension of general CBS method which was initially introduced in finite element framework, this is due to the fact that CBS method not only can enhance the stability, but also avoid LBB condition when equal order basis function is used to approximate velocity and pressure variables. Meanwhile, a simple way for node update and node speed calculation is developed which is used to capture the free surface exactly. The numerical solutions are compared with available analytical and numerical solutions, which shows that the proposed method has better ability to simulate the free surface incompressible flow problem. Copyright © 2009 John Wiley & Sons, Ltd.     

h-ADAPTIVITY ANALYSIS BASED ON MULTIPLE SCALE REPRODUCING KERNEL PARTICLE METHOD

IntroductionAs one of the meshfree methods, multiple scale reproducing kernel particle method(RKPM)[1,2]bears, besides the common features of all meshfree methods, a uniqueinherent feature of multiresolution analysis by which the structural response can b…     

基于多尺度再生核质点法的h型自适应分析

基于多尺度再生核质点法(RKPM)的多分辨分析特性，给出了适于求解二维问题的h型自适应分析方法，讨论了不同高梯度指标对检测高梯度区的影响，同时还提出了基于邻节点搜索(SNN)及局部Delaunay三角分解(LDT)技术的适于任意节点分布的无网格自适应分析节点加密技术．通过对二维线弹性平面应力及平板弯曲问题的h型自适应分析，说明了不合适的高梯度指标将对自适应分析的收敛造成不良影响，且LDT节点加密技术较SNN技术的效率更高．数值算例的结果说明自适应分析的有效性、稳定性及良好的收敛特性．