Development of a meshless hybrid boundary node method for Stokes flows

The meshless hybrid boundary node method (HBNM) is a promising method for solving boundary value problems, and is further developed and numerically implemented for incompressible 2D and 3D Stokes flows in this paper. In this approach, a new modified variational formulation using a hybrid functional is presented. The formulation is expressed in terms of domain and boundary variables. The moving least-squares (MLS) method is employed to approximate the boundary variables whereas the domain variables are interpolated by the fundamental solutions of Stokes equation, i.e. Stokeslets. The present method only requires scatter nodes on the surface, and is a truly boundary type meshless method as it does not require the ‘boundary element mesh’, either for the purpose of interpolation of the variables or the integration of ‘energy’. Moreover, since the primitive variables, i.e., velocity vector and pressure, are employed in this approach, the problem of finding the velocity is separated from that of finding pressure. Numerical examples are given to illustrate the implementation and performance of the present method. It is shown that the high convergence rates and accuracy can be achieved with a small number of nodes.     

Fast HdBNM for large-scale thermal analysis of CNT-reinforced composites

Because of their high thermal conductivities, carbon nanotubes (CNT) have promising potential in development of fundamentally new composites. To study the influence of CNTs distribution on the overall properties of a composite, the modeling of a Representative Volume Element (RVE) including a large number of CNTs that are randomly distributed and oriented is necessary. However, analysis of such a RVE using standard numerical methods faces two severe difficulties, namely the discretization of the geometry and a very large computational scale. In this paper, the first difficulty is alleviated by employing the Hybrid Boundary Node Method (HdBNM), which is a form of the boundary type meshless methods. To overcome the second difficulty, the Fast Multipole Method (FMM) is combined with the HdBNM to solve a simplified mathematical model. RVEs containing various numbers of CNTs with different lengths, shapes and alignments have been analyzed, resulting in valuable insights gained into the thermal behavior of the composite material.     

Meshless analysis of two-dimensional Stokes flows with the Galerkin boundary node method

In this paper, the Galerkin boundary node method (GBNM) is developed for the solution of stationary Stokes problems in two dimensions. The GBNM is a boundary only meshless method that combines a variational form of boundary integral formulations for governing equations with the moving least-squares (MLS) approximations for construction of the trial and test functions. Boundary conditions in this approach are included into the variational form, thus they can be applied directly and easily despite the MLS shape functions lack the property of a delta function. Besides, the GBNM keeps the symmetry and positive definiteness of the variational problems. Convergence analysis results of both the velocity and the pressure are given. Some selected numerical tests are also presented to demonstrate the efficiency of the method.     

A hybrid radial boundary node method based on radial basis point interpolation

The hybrid boundary node method (HBNM) retains the meshless attribute of the moving least squares (MLS) approximation and the reduced dimensionality advantages of the boundary element method. However, the HBNM inherits the deficiency of the MLS approximation, in which shape functions lack the delta function property. Thus in the HBNM, boundary conditions are implemented after they are transformed into their approximations on the boundary nodes with the MLS scheme.This paper combines the hybrid displacement variational formulation and the radial basis point interpolation to develop a direct boundary-type meshless method, the hybrid radial boundary node method (HRBNM) for two-dimensional potential problems. The HRBNM is truly meshless, i.e. absolutely no elements are required either for interpolation or for integration. The radial basis point interpolation is used to construct shape functions with delta function property. So unlike the HBNM, the HRBNM is a direct numerical method in which the basic unknown quantity is the real solution of nodal variables, and boundary conditions can be applied directly and easily, which leads to greater computational precision. Some selected numerical tests illustrate the efficiency of the method proposed.     

Reliability-based design optimization of knuckle component using conservative method of moving least squares meta-models

This paper discusses reliability-based design optimization (RBDO) of an automotive knuckle component under bump and brake loading conditions. The probabilistic design problem is to minimize the weight of a knuckle component subject to stresses, deformations, and frequency constraints in order to meet the given target reliability. The initial design is generated based on an actual vehicle specification. The finite element analysis is conducted using ABAQUS, and the probabilistic optimal solutions are obtained via the moving least squares method (MLSM) in the context of approximate optimization. For the meta-modeling of inequality constraint functions, a constraint-feasible moving least squares method (CF-MLSM) is used in the present study. The method of CF-MLSM based RBDO has been shown to not only ensure constraint feasibility in a case where the meta-model-based RBDO process is employed, but also to require low expense, as compared with both conventional MLSM and non-approximate RBDO methods.     

Fast multipole DBEM analysis of fatigue crack growth

A fast multipole method (FMM) based on complex Taylor series expansions is applied to the dual boundary element method (DBEM) for large-scale crack analysis in linear elastic fracture mechanics. Combining multipole expansions with local expansions, both the computational complexity and memory requirement are reduced to O(N), where N is the number of DOF. An incremental crack-extension analysis based on the maximum principal stress criterion and the Paris law is used to simulate the fatigue growth of numerous cracks in a 2D solid. Some examples are presented to validate the numerical scheme.     

A node enrichment adaptive refinement in Discrete Least Squares Meshless method for solution of elasticity problems

In this paper, an adaptive refinement procedure is proposed to be used with Discrete Least Squares Meshless (DLSM) method for accurate solution of planar elasticity problems. DLSM method is a newly introduced meshless method based on the least squares concept. The method is based on the minimization of a least squares functional defined as the weighted summation of the squared residual of the governing differential equation and its boundary conditions at nodal points used to discretize the domain and its boundaries. A Moving Least Square (MLS) method is used to construct the shape function making the approach a fully least squares based approach. An error estimate and adaptive refinement strategy is proposed in this paper to increase the efficiency of the DLSM method. For this, a residual based error estimator is introduced and used to discover the region of higher errors. The proposed error estimator has the advantages of being available at the end of each analysis contributing to the efficiency of the proposed method. An enrichment method is then used by adding more nodes to the area of higher errors as indicated by the error estimator. A Voronoi diagram is used to locate the position of the nodes to be added to the current nodal configuration. Efficiency and effectiveness of the proposed procedure is examined by adaptively solving two benchmark problems. The results show the ability of the proposed strategy for accurate simulation of elasticity problems.     

A posteriori error estimates and adaptive procedures for the meshless Galerkin boundary node method for 3D potential problems

The Galerkin boundary node method (GBNM) is a boundary only meshless method that combines variational formulations of boundary integral equations with the moving least-squares approximations. This paper presents the mathematical derivation of a posteriori error estimates and adaptive refinement procedures for the GBNM for 3D potential problems. Two types of error estimators are developed in detail. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two successive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the GBNM solution itself and its L²-orthogonal projection. The reliability and efficiency of both types of error estimators is established. That is, these error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization technique is introduced to accommodate the non-local property of integral operators for the needed local and computable a posteriori error indicators. Convergence analysis results of corresponding adaptive meshless procedures are also given. Numerical examples with high singularities illustrate the theoretical results and show that the proposed adaptive procedures are simple, effective and efficient.     

New implementation of MLBIE method for heat conduction analysis in functionally graded materials

The meshless local boundary integral equation (MLBIE) method with an efficient technique to deal with the time variable are presented in this article to analyze the transient heat conduction in continuously nonhomogeneous functionally graded materials (FGMs). In space, the method is based on the local boundary integral equations and the moving least squares (MLS) approximation of the temperature and heat flux. In time, again the MLS approximates the equivalent Volterra integral equation derived from the heat conduction problem. It means that, the MLS is used for approximation in both time and space domains, and we avoid using the finite difference discretization or Laplace transform methods to overcome the time variable. Finally the method leads to a single generalized Sylvester equation rather than some (many) linear systems of equations. The method is computationally attractive, which is shown in couple of numerical examples for a finite strip and a hollow cylinder with an exponential spatial variation of material parameters.     

An Adaptive Fast Multipole Boundary Element Method for Three-dimensional Potential Problems

An adaptive fast multipole boundary element method (FMBEM) for general three-dimensional (3-D) potential problems is presented in this paper. This adaptive FMBEM uses an adaptive tree structure that can balance the multipole to local translations (M2L) and the direct evaluations of the near-field integrals, and thus can reduce the number of the more costly direct evaluations. Furthermore, the coefficients used in the preconditioner for the iterative solver (GMRES) are stored and used repeatedly in the direct evaluations of the near-field contributions. In this way, the computational efficiency of the adaptive FMBEM is improved significantly. The adaptive FMBEM can be applied to both the original FMBEM formulation and the new FMBEM with diagonal translations. Several numerical examples are presented to demonstrate the efficiency and accuracy of the adaptive FMBEM for studying large-scale 3-D potential problems. The adaptive FMBEM is found to be about 50% faster than the non-adaptive version of the new FMBEM in solving the model (with 558,000 elements) for porous materials studied in this paper. The computational efficiencies and accuracies of the FMBEM as compared with the finite element method (FEM) are also studied using a heat-sink model. It is found that the adaptive FMBEM is especially advantageous in modeling problems with complicated domains for which free meshes with much more finite elements would be needed with the FEM.     

An Exact Asymptotically Efficient Confidence Bound for Reliability in the Case of the Weibull Distribution

This paper presents a simple method for obtaining exact lower confidence bounds for reliabilities (tail probabilities) for items whose life times follow a Weibull distribution where both the “shape” and “scale” parameters are unknown. These confidence bounds are obtained both for the censored and non-censored cases and are asymptotically efficient. They are exact even for small sample sizes in that they attain the desired confidence level precisely. The case of an additional unknown “location” or “shift” parameter is also discussed in the large sample case. Tables are given of exact and asymptotic lower confidence bounds for the reliability for sample sizes of 10, 15, 20, 30, 50 and 100 for various censoring fractions.     

Application of the multi-level time-harmonic fast multipole BEM to 3-D visco-elastodynamics

This article extends previous work by the authors on the single- and multi-domain time-harmonic elastodynamic multi-level fast multipole BEM formulations to the case of weakly dissipative viscoelastic media. The underlying boundary integral equation and fast multipole formulations are formally identical to that of elastodynamics, except that the wavenumbers are complex-valued due to attenuation. Attention is focused on evaluating the multipole decomposition of the viscoelastodynamic fundamental solution. A damping-dependent modification of the selection rule for the multipole truncation parameter, required by the presence of complex wavenumbers, is proposed. It is empirically adjusted so as to maintain a constant accuracy over the damping range of interest in the approximation of the fundamental solution, and validated on numerical tests focusing on the evaluation of the latter. The proposed modification is then assessed on 3D single-region and multi-region visco-elastodynamic examples for which exact solutions are known. Finally, the multi-region formulation is applied to the problem of a wave propagating in a semi-infinite medium with a lossy semi-spherical inclusion (seismic wave in alluvial basin). These examples involve problem sizes of up to about 3×10⁵ boundary unknowns.     

A fast multipole boundary element method for solving the thin plate bending problem

A fast multipole boundary element method (BEM) for solving large-scale thin plate bending problems is presented in this paper. The method is based on the Kirchhoff thin plate bending theory and the biharmonic equation governing the deflection of the plate. First, the direct boundary integral equations and the conventional BEM for thin plate bending problems are reviewed. Second, the complex notation of the kernel functions, expansions and translations in the fast multipole BEM are presented. Finally, a few numerical examples are presented to show the accuracy and efficiency of the fast multipole BEM in solving thin plate bending problems. The bending rigidity of a perforated plate is evaluated using the developed code. It is shown that the fast multipole BEM can be applied to solve plate bending problems with good accuracy. Possible improvements in the efficiency of the method are discussed.