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.     

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…     

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     

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

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

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.     

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.     

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.     

一种基于颗粒接触的二维无网格方法及其在高速冲击模拟中的应用

为解决基于连续介质力学的离散元方法(CDEM)在高速冲击模拟中因网格畸变导致的系统能量发散问题,提出了一种基于颗粒接触的二维无网格方法(PCMM)。该方法基于颗粒间复杂丰富的接触信息构建三角形单元,通过接触对的演化更新实现旧单元(满足删除条件的单元)的删除及新单元(满足创建条件的单元)的重建,通过在单元内引入流体弹塑性模型实现高速冲击问题的模拟。给出了三角形单元创建的3个必备条件:组成单元的3个颗粒必须彼此接触,任意一个内角必须在30°~150°之间,任意一条边长必须大于平均半径的0.5倍。弹性杆撞击、泰勒杆、碎片云、子弹入射靶板等算例的结果表明了PCMM方法在模拟高速冲击问题方面的正确性及合理性。