Sharif Digital Repository

This paper presents three schemes of 2D meshless finite volume (MFV) method, referred to as MFV with overlapping control volumes (MFV1), MFV with irregular non-overlapping control volumes (MFV2) and MFV with regular non-overlapping control volumes (MFV3). The methods utilize the local symmetric weak form of system equation and the interpolation functions constructed using the weighted multi-triangles method (WMTM) which is recently developed by the present authors. The proposed formulation involves only integrals over the boundaries of control volumes. The performance of the proposed schemes is studied in three benchmark problems. A comparative study between the predictions of the above MFV... 

The rupturing process of a membrane, located between two fluids at the center of a three-dimensional channel, is numerically investigated. The smoothed particle hydrodynamics (SPH) and the finite element method (FEM) are used, respectively, for modeling the fluid and solid phases. A range of pressure differences and membrane thicknesses are studied and two different rupturing processes are identified. These processes differ in the time scale of the rupture, the location of the rupture initiation, the level of destruction and the driving mechanism. © 2019 by the authors  

In this paper, two adaptive refinement techniques are presented for enhancing the capability of the new kind of finite volume method called meshless finite volume (MFV) method in which the moving least squares (MLS) approximation technique is employed. The proposed algorithms perform by inserting new nodes in large error regions identified using the Zienkiewicz-Zhu (Z-Z) error estimator and the T-Belytschko (TB) stress recovery scheme. In the first method referred as CV-based adaptive refinement, the new nodes are inserted at the vertices of control volumes with higher errors. The second method, referred as GA-based adaptive refinement, contains two schemes where an adaptive refinement... 

A major disadvantage of conventional meshless methods as compared to finite element method (FEM) is their weak performance in dealing with constraints. To overcome this difficulty, the penalty and Lagrange multiplier methods have been proposed in the literature. In the penalty method, constraints cannot be enforced exactly. On the other hand, the method of Lagrange multiplier leads to an ill-conditioned matrix which is not positive definite. The aim of this paper is to boost the effectiveness of the conventional reproducing kernel particle method (RKPM) in handling those types of constraints which specify the field variable and its gradient(s) conveniently. Insertion of the gradient term(s),... 

It is well-known that the conventional reproducing kernel particle method (RKPM) is unfavorable when dealing with the derivative type essential boundary conditions [1-3]. To remedy this issue a group of meshless methods in which the derivatives of a function can be incorporated in the formulation of the corresponding interpolation operator will be discussed. Formulation of generalized moving least squares (GMLS) on a domain and GMLS on a finite set of points will be presented. The generalized RKPM will be introduced as the discretized form of GMLS on a domain. Another method that helps to deal with derivative type essential boundary conditions is the gradient RKPM which incorporates the... 

In this paper, the truly Meshless Local Petrov-Galerkin (MLPG) method is extended for computation of unsteady incompressible flows, governed by the Navier-Stokes equations (NSE), in vorticity-stream function formulation. The present method is a truly meshless method based only on a number of randomly located nodes. The formulation is based on two equations including stream function Poisson equation and vorticity advection-dispersion-reaction eq. (ADRE). The meshless method is based on a local weighted residual method with the Heaviside step function and quartic spline as the test functions respectively over a local subdomain. Moving Least Square approximation (MLS) is employed in shape... 

In this paper, an anisotropic weight function in the elliptic form is introduced for the Element Free Galerkin Method (EFGM). In the circular (isotropic) weight function, each node has one characteristic parameter that determines its domain of influence. In the elliptic weight function, each node has three characteristic parameters that are major influence radius, minor influence radius and the direction of the major influence. Using the elliptic weight function each point of the domain may be affected by a less number of nodes in certain conditions. Thus, the computational cost of the method is decreased. In addition, the dependency of the solution on the method that is used for the... 

The potential application of SWCNTs as mass nanosensors is examined for a wide range of boundary conditions. The SWCNT is modeled via nonlocal Rayleigh, Timoshenko, and higher-order beam theories. The added nano-objects are considered as rigid solids, which are attached to the SWCNT. The mass weight and rotary inertial effects of such nanoparticles are appropriately incorporated into the nonlocal equations of motion of each model. The discrete governing equation pertinent to each model is obtained using an effective meshless technique. The key factor in design of a mass nanosensor is to determine the amount of frequency shift due to the added nanoparticles. Through an inclusive parametric... 

A new meshless method called gradient reproducing kernel particle method (GRKPM) is proposed for numerical solutions of one-dimensional Burgers' equation with various values of viscosity and different initial and boundary conditions. Discretization is first done in the space via GRKPM, and subsequently, the reduced system of nonlinear ordinary differential equations is discretized in time by the Gear's method. Comparison with the exact solutions, which are only available for restricted initial conditions and values of viscosity, approves the efficacy of the proposed method. For challenging cases involving small viscosities, comparison with the results obtained using other numerical schemes... 

The meshless element-free Galerkin method was generalized and an algorithm was developed for 3D dynamic modeling of deformable bodies in real time. The efficacy of the algorithm was investigated in a 3D linear viscoelastic model of human spleen subjected to a time-varying compressive force exerted by a surgical grasper. The model remained stable in spite of the considerably large deformations occurred. There was a good agreement between the results and those of an equivalent finite element model. The computational cost, however, was much lower, enabling the proposed algorithm to be effectively used in real-time applications