A meshless solution of two-dimensional unsteady flow

The boundary element method (BEM) is a very useful numerical method for groundwater flow models. Particularly, this method was used to solve problems in homogeneous domains. However, it presents even greater difficulties than the other numerical methods when coping with non-homogeneities which are so characteristic in the groundwater hydraulics. Recently, meshless method which is based on a local boundary integral approach is introduced. It uses distributed nodal points, covering the domain. These points can be randomly spread over the domain. Every node is surrounded by a simple surface centered at the collocation point and the boundary integral equation is written on this local boundary. The unknown variables, in the local sub-domains, are approximated by some of the interpolation method. In this paper the combination of radial basis functions and the dual reciprocity method is used to solve the time-dependent groundwater flow.     

Meshless boundary integral equations with equilibrium satisfaction

A known feature of any mixed interpolation boundary integral equations (BIE)-based methods is that equilibrium is not generally guaranteed in the numerical solution. Here, a complete meshless technique, based on the boundary element-free method (BEFM) with complete equilibrium satisfaction for 2D elastostatic analysis is proposed. The BEFM adopted is a meshless method based on boundary integral equations such as local boundary integral equation (LBIE) method and boundary node method (BNM), differing from them with respect to the integration domain and the approximation scheme.     

A Trefftz function approximation in local boundary integral equations

 In the present paper the Trefftz function as a test function is used to derive the local boundary integral equations (LBIE) for linear elasticity. Since Trefftz functions are regular, much less requirements are put on numerical integration than in the conventional boundary integral method. The moving least square (MLS) approximation is applied to the displacement field. Then, the traction vectors on the local boundaries are obtained from the gradients of the approximated displacements by using Hooke's law. Nodal points are randomly spread on the domain of the analysed body. The present method is a truly meshless method, as it does not need a finite element mesh, either for purposes of interpolation of the solution variables, or for the integration of the energy. Two ways are presented to formulate the solution of boundary value problems. In the first one the local boundary integral equations are written in all nodes (interior and boundary nodes). In the second way the LBIE are written only at the interior nodes and at the nodes on the global boundary the prescribed values of displacements and/or tractions are identified with their MLS approximations. Numerical examples for a square patch test and a cantilever beam are presented to illustrate the implementation and performance of the present method. Received 6 November 2000     

A meshless local boundary integral equation method for dynamic anti-plane shear crack problem in functionally graded materials

This paper presents a meshless local boundary integral equation method (LBIEM) for dynamic analysis of an anti-plane crack in functionally graded materials (FGMs). Local boundary integral equations (LBIEs) are formulated in the Laplace-transform domain. The static fundamental solution for homogeneous elastic solids is used to derive the local boundary–domain integral equations, which are applied to small sub-domains covering the analyzed domain. For the sub-domains a circular shape is chosen, and their centers, the nodal points, correspond to the collocation points. The local boundary–domain integral equations are solved numerically in the Laplace-transform domain by a meshless method based on the moving least–squares (MLS) scheme. Time-domain solutions are obtained by using the Stehfest's inversion algorithm. Numerical examples are given to show the accuracy of the proposed meshless LBIEM.     

An advanced meshless LBIE/RBF method for solving two-dimensional incompressible fluid flows

The present work presents a meshless local boundary integral equation (LBIE) method for the solution of two-dimensional incompressible fluid flow problems governed by the Navier–Stokes equations. The method uses, for its meshless implementation, nodal points spread over the analyzed domain and employs in an efficient way the radial basis functions (RBF) for the interpolation of the interior and boundary variables. The inverse matrix of the RBF is computed only once for every nodal point and the interpolation functions are evaluated by the inner product of the inverse matrix with the weight vector associated to the integration point. This technique leads to a fast and efficient meshless approach, the locality of the method is maintained and the system matrices are banded with small bandwidth. The velocity–vorticity approach of the Navier–Stokes equations is adopted and the LBIEs are derived for the velocity and the vorticity field, resulting in a very stable and accurate implementation. The evaluation of the volume integrals is accomplished via a very efficient and accurate technique by triangularizing the local area of the nodal point to the minimum number of well formed triangles. Numerical examples illustrate the proposed methodology and demonstrate its accuracy.     

A meshless local boundary integral equation method for two-dimensional steady elliptic problems

A novel meshless local boundary integral equation (LBIE) method is proposed for the numerical solution of two-dimensional steady elliptic problems, such as heat conduction, electrostatics or linear elasticity. The domain is discretized by a distribution of boundary and internal nodes. From this nodal points’ cloud a “background” mesh is created by a triangulation algorithm. A local form of the singular boundary integral equation of the conventional boundary elements method is adopted. Its local form is derived by considering a local domain of each node, comprising by the union of neighboring “background” triangles. Therefore, the boundary shape of this local domain is a polygonal closed line. A combination of interpolation schemes is taken into account. Interpolation of boundary unknown field variables is accomplished through boundary elements’ shape functions. On the other hand, the Radial Basis Point Interpolation Functions method is employed for interpolating the unknown interior fields. Essential boundary conditions are imposed directly due to the Kronecker delta-function property of the boundary elements’ interpolation functions. After the numerical evaluation of all boundary integrals, a banded stiffness matrix is constructed, as in the finite elements method. Several potential and elastostatic benchmark problems in two dimensions are solved numerically. The proposed meshless LBIE method is also compared with other numerical methods, in order to demonstrate its efficiency, accuracy and convergence.     

Micromechanics of fibrous composites subjected to combined shear and thermal loading using a truly meshless method

In this study, a micromechanical model is presented to study the combined normal, shear and thermal loading of unidirectional (UD) fiber reinforced composites. An appropriate truly meshless method based on the integral form of equilibrium equations is also developed. This meshless method formulated for the generalized plane strain assumption and employed for solution of the governing partial differential equations of the problem. The solution domain includes a representative volume element (RVE) consists of a fiber surrounded by corresponding matrix in a square array arrangement. A direct interpolation method is employed to enforce the appropriate periodic boundary conditions for the combined thermal, transverse shear and normal loading. The fully bonded fiber–matrix interface condition is considered and the displacement continuity and traction reciprocity are imposed to the fiber–matrix interface. Predictions show excellent agreement with the available experimental, analytical and finite element studies. Comparison of the CPU time between presented method and the conventional meshless local Petrov–Galerkin (MLPG) shows significant reduction of the computational time. The results of this study also revealed that the presented model could provide highly accurate predictions with relatively small number of nodes and less computational time without the complexity of mesh generation.     

Inverse heat conduction problems by meshless local Petrov–Galerkin method

The meshless local Petrov–Galerkin (MLPG) method is used to solve stationary and transient heat conduction inverse problems in 2-D and 3-D axisymmetric bodies. A 3-D axisymmetric body is generated by rotating a cross section around an axis of symmetry. Axial symmetry of geometry and boundary conditions reduce the original 3-D boundary value problem to a 2-D problem. The analyzed domain is covered by small circular subdomains surrounding nodes randomly spread over the analyzed domain. A unit step function is chosen as test function in deriving the local integral equations (LIEs) on the boundaries of the chosen subdomains. The time integration schemes are formulated based on the Laplace transform technique and the time difference approach, respectively. The local integral equations are non-singular and take a very simple form. Spatial variation of the temperature and heat flux (or of their Laplace transforms) at discrete time instants are approximated on the local boundary and in the interior of the subdomain by means of the moving least-squares (MLS) method. Singular value decomposition (SVD) is applied to solve the ill-conditioned linear system of algebraic equations obtained from the LIE after MLS approximation. The Stehfest algorithm is applied for the numerical Laplace inversion, in order to retrieve the time-dependent solutions.     

Heat Conduction Analysis of 3-D Axisymmetric and Anisotropic FGM Bodies by Meshless Local Petrov–Galerkin Method

The meshless local Petrov–Galerkin method is used to analyze transient heat conduction in 3-D axisymmetric solids with continuously inhomogeneous and anisotropic material properties. A 3-D axisymmetric body is created by rotation of a cross section around an axis of symmetry. Axial symmetry of geometry and boundary conditions reduces the original 3-D boundary value problem into a 2-D problem. The cross section is covered by small circular subdomains surrounding nodes randomly spread over the analyzed domain. A unit step function is chosen as test function, in order to derive local integral equations on the boundaries of the chosen subdomains, called local boundary integral equations. These integral formulations are either based on the Laplace transform technique or the time difference approach. The local integral equations are nonsingular and take a very simple form, despite of inhomogeneous and anisotropic material behavior across the analyzed structure. Spatial variation of the temperature and heat flux (or of their Laplace transforms) at discrete time instants are approximated on the local boundary and in the interior of the subdomain by means of the moving least-squares method. The Stehfest algorithm is applied for the numerical Laplace inversion, in order to retrieve the time-dependent solutions.     

A meshless solution to two-dimensional convection–diffusion problems

An integral equation domain decomposition method has been implemented in a meshless fashion. The method exploits the advantage of placing the source point always in the centre of circular sub-domains in order to avoid singular or near-singular integrals. Three equations for two-dimensional (2D) or four for three-dimensional (3D) potential problems are required at each node. The first equation is the integral equation arising from the application of the Green's identities and the remaining equations are the derivatives of the first equation in respect to space coordinates. Radial basis function interpolation is applied in order to obtain the values of the field variable and partial derivatives at the boundary of the circular sub-domains, providing this way the boundary conditions for solution of the integral equations at the nodes (centres of circles). Dual reciprocity method (DRM) has been applied to convert the domain integrals into boundary integrals, though the approach is general and can be applied without the DRM. The accuracy and robustness of the method has been tested on a convection–diffusion problem. The results obtained using the current approach have been compared with previously reported results obtained using the finite element method (FEM), and the DRM multi-domain approach (DRM-MD) showing similar level of accuracy.     

Analyses of Circular Magnetoelectroelastic Plates with Functionally Graded Material Properties

A meshless method based on the local Petrov–Galerkin approach is proposed for plate bending analysis with material containing functionally graded magnetoelectroelastic properties. Material properties are considered to be continuously varying along the plate thickness. Axial symmetry of geometry and boundary conditions for a circular plate reduces the original 3D boundary value problem into a 2D problem in axial cross section. Both stationary and transient dynamic conditions for a pure mechanical load are considered in this article. The local weak formulation is employed on circular subdomains in the axial cross section. Subdomains surrounding nodes are randomly spread over the analyzed domain. The test functions are taken as unit step functions in derivation of the local integral equations (LIEs). The moving least-squares (MLS) method is adopted for the approximation of the physical quantities in the LIEs. After performing the spatial integrations, one obtains a system of ordinary differential equations for certain nodal unknowns. That system is solved numerically by the Houbolt finite-difference scheme as a time-stepping method.     

Three dimensional static and dynamic analysis of thick functionally graded plates by the meshless local Petrov–Galerkin (MLPG) method

In this paper, three dimensional (3D) static and dynamic analysis of thick functionally graded plates based on the Meshless Local Petrov–Galerkin (MLPG) is presented. Using the kinematics of a three-dimensional continuum, the local weak form of the equilibrium equations is derived. A weak formulation for the set of governing equations is transformed into local integral equations on local sub-domains using a Heaviside step function as test function. In this case, governing equations corresponding to the stiffness matrix do not involve any domain integration or singular integrals. Nodal points are distributed in the 3D analyzed domain and each node is surrounded by a cubic sub-domain to which a local integral equation is applied. The meshless approximation based on the three-dimensional Moving Least-Square (MLS) is employed as shape function to approximate the field variable of scattered nodes in the problem domain. The Newmark time integration method is used to solve the system of coupled second-order ODEs. Effective material properties of the plate, made of two isotropic constituents with volume fractions varying only in the thickness direction, are computed using the Mori–Tanaka homogenization technique. Numerical examples for solving the static and dynamic response of elastic thick functionally graded plates are demonstrated. As a result, the distributions of the deflection and stresses through the plate thickness are presented for different material gradients and boundary conditions. The effects of the volume fractions of the constituents on the centroidal deflection are also investigated. The numerical efficiency of the proposed meshless method is illustrated by the comparison of results obtained from previous literatures.     

A local Petrov-Galerkin approach with moving Kriging interpolation for solving transient heat conduction problems

A meshless Local Petrov-Galerkin approach based on the moving Kriging interpolation (Local Kriging method; LoKriging hereafter) is employed for solving partial different equations that govern the heat flow in two- and three-dimensional spaces. The method is developed based on the moving Kriging interpolation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for domain integral in symmetric weak form. As the shape functions possess the Kronecker delta function property, essential boundary conditions can be implemented without any difficulties. The traditional two-point difference method is selected for the time discretization scheme. For computation of 3D problems, a novel local sub-domain from the polyhedrons is used for evaluating the integrals. Several selected numerical examples are presented to illustrate the performance of the LoKriging method.     

A meshless local boundary integral equation method for solving transient elastodynamic problems

A new local boundary integral equation (LBIE) method for solving two dimensional transient elastodynamic problems is proposed. The method utilizes, for its meshless implementation, nodal points spread over the analyzed domain and employs the moving least squares (MLS) approximation for the interpolation of the interior and boundary variables. On the global boundary, displacements and tractions are treated as independent variables. The local integral representation of displacements at each nodal point contains both surface and volume integrals, since it employs the simple elastostatic fundamental solution and considers the acceleration term as a body force. On the local boundaries, tractions are avoided with the aid of the elastostatic companion solution. The collocation of the local boundary/volume integral equations at all the interior and boundary nodes leads to a final system of ordinary differential equations, which is solved stepwise by the -Wilson finite difference scheme. Direct numerical techniques for the accurate evaluation of both surface and volume integrals are employed and presented in detail. All the strongly singular integrals are computed directly through highly accurate integration techniques. Three representative numerical examples that demonstrate the accuracy of the proposed methodology are provided.     

Inverse fracture problems in piezoelectric solids by local integral equation method

The meshless local Petrov–Galerkin (MLPG) method is used to solve the inverse fracture problems in two-dimensional (2D) piezoelectric body. Electrical boundary conditions on the crack surfaces are not specified due to unknown dielectric permittivity of the medium inside the crack. Both stationary and transient dynamic boundary conditions are considered here. The analyzed domain is covered by small circular subdomains surrounding nodes spread randomly over the analyzed domain. A unit step function is chosen as test function in deriving the local integral equations (LIE) on the boundaries of the chosen subdomains. The Laplace-transform technique is applied to eliminate the time variation in the governing equation. The local integral equations are nonsingular and take a very simple form. The spatial variation of the Laplace transforms of displacements and electrical potential are approximated on the local boundary and in the interior of the subdomain by means of the moving least-squares (MLS) method. The singular value decomposition (SVD) is applied to solve the ill-conditioned linear system of algebraic equations obtained from the LIE after MLS approximation. The Stehfest algorithm is applied for the numerical Laplace inversion to retrieve the time-dependent solutions.     

Inverse heat conduction problems in three-dimensional anisotropic functionally graded solids

A meshless method based on the local Petrov–Galerkin approach is applied to inverse transient heat conduction problems in three-dimensional solids with continuously inhomogeneous and anisotropic material properties. The Heaviside step function is used as a test function in the local weak form, leading to the derivation of local integral equations. Nodal points are randomly distributed in the domain analyzed, and each node is surrounded by a spherical subdomain in which a local integral equation is applied. A meshless approximation based on the moving least-squares method is employed in the implementation. After performing spatial integrations, we obtain a system of ordinary differential equations for certain nodal unknowns. A backward finite-difference method is used for the approximation of the diffusive term in the heat conduction equation. A truncated singular-value decomposition is used to solve the ill-conditioned linear system of algebraic equations at each time step. The effectiveness of the meshless local Petrov–Galerkin (MLPG) method for this inverse problem is demonstrated by numerical examples.     

Application of local boundary integral equation method into micropolar elasticity

A new meshless method for solving boundary value problems in micropolar elasticity is presented. The method is based on the local boundary integral equation (LBIE) method with the moving least squares approximation of physical quantities. Randomly scattered nodes are utilized for interpolation of field data. Every node is surrounded by a simple surface centered at the collocation point in the LBIE method. On the surface of subdomains the LBIEs are written. Fundamental solutions corresponding to uncoupled governing equations are derived. To eliminate the traction vector in the LBIE, the modified fundamental solution is introduced.     

Elastodynamic problems by meshless local integral method: Analytical formulation

In this paper, analytical forms of integrals in the meshless local integral equation method in the Laplace space are derived and implemented for elastodynamic problems. The meshless approximation based on the radial basis function (RBF) is employed for implementation of displacements. A weak form of governing equations with a unit test function is transformed into local integral equations. A completed set of the local boundary integrals are obtained in closed form. As the closed form of the local boundary integrals are obtained, there are no domain or boundary integrals to be calculated numerically. Several examples including dynamic fracture mechanics problems are presented to demonstrate the accuracy of the proposed method in comparison with analytical solutions and the boundary element method.     

A meshless method for two-dimensional diffusion equation with an integral condition

This paper presents a new approach based on the meshless local Petrov–Galerkin (MLPG) and collocation methods to treat the parabolic partial differential equations with non-classical boundary conditions. In the presented method, the MLPG method is applied to the interior nodes while the meshless collocation method is applied to the nodes on the boundaries, and so the Dirichlet boundary condition is imposed directly. To treat the complicated integral boundary condition appearing in the problem, Simpson's composite numerical integration rule is applied. A time stepping scheme is employed to approximate the time derivative. Finally, two numerical examples are presented showing the behavior of the solution and the efficiency of the proposed method.     

Analysis of Microelectromechanical Systems (MEMS) Devices by the Meshless Point Weighted Least-squares Method

In this paper, the characteristics of microelectromechanical systems (MEMS) devices are analyzed by a meshless method—point weighted least-squares (PWLS) method. In the present meshless method, field nodes and collocation points are adopted. The field nodes are used to construct the trial functions based on locally supported interpolation domains. The collocation points that can be independent of the field nodes are adopted to evaluate the total residuals of the problem domain and its boundaries. The least-squares technique is used to obtain the solution of the problem by minimizing the functional of the summation of weighted residuals. The present meshless method possesses some advantages compared with the conventional collocation methods, e.g., it is very stable for both regularly or irregularly nodal distributions; the displacement and derivative boundary conditions can be easily enforced; and the final coefficient matrix is symmetric. Several one-dimensional and two- dimensional MEMS devices that are governed by the nonlinear equations are studied by the present PWLS method. The simulated results are compared with those obtained by other simulation approaches and experimental results. It is shown that the PWLS method is very efficient and accurate for the analysis of MEMS devices.