Fast Method of Particular Solutions for Solving Partial Differential Equations

Method of particular solutions (MPS) has been implemented in many science and engineering problems but obtaining the closed-form particular solutions, the selection of the good shape parameter for various radial basis functions (RBFs) and simulation of the large-scale problems are some of the challenges which need to overcome. In this dissertation, we have used several techniques to overcome such challenges. The closed-form particular solutions for the Matérn and Gaussian RBFs were not known yet. With the help of the symbolic computational tools, we have derived the closed-form particular solutions of the Matérn and Gaussian RBFs for the Laplace and biharmonic operators in 2D and 3D. These derived particular solutions play an important role in solving inhomogeneous problems using MPS and boundary methods such as boundary element methods or boundary meshless methods. In this dissertation, to select the good shape parameter, various existing variable shape parameter strategies and some well-known global optimization algorithms have also been applied. These good shape parameters provide high accurate solutions in many RBFs collocation methods. Fast method of particular solutions (FMPS) has been developed for the simulation of the large-scale problems. FMPS is based on the global version of the MPS. In this method, partial differential equations are discretized by the usual MPS and the determination of the unknown coefficients is accelerated using a fast technique. Numerical results confirm the efficiency of the proposed technique for the PDEs with a large number of computational points in both two and three dimensions. We have also solved the time fractional diffusion equations by using MPS and FMPS

Numerical solution of second order linear hyperbolic telegraph equation

This paper is of about a numerical solution of the second order linear hyperbolic telegraph equation. To solve numerically the second order linear hyperbolic telegraph equation, the cubic B-spline collocation method is used in space discretization and the fourth order one-step method is used in time discretization. By using the fourth order one-step method, it is aimed to obtain a numerical algorithm whose accuracy is higher than the current studies. The efficiency and accuracy of the proposed method is studied by two examples. The obtained results show that the proposed method has higher accuracy as intended.This work has been supported by the Scientific Research Council of Eskisehir Osmangazi University under project No. 2018-2090.Publisher's Versio

A high-performance boundary element method and its applications in engineering

As a semi-numerical and semi-analytical method, owing to the inherent advantage, of boundary-only discretisation, the boundary element method (BEM) has been widely applied to problems with complicated geometries, stress concentration problems, infinite domain problems, and many others. However, domain integrals and non-symmetrical and dense matrix systems are two obstacles for BEM which have hindered the its further development and application. This thesis is aimed at proposing a high-performance BEM to tackle the above two drawbacks and broaden the application scope of BEM. In this thesis, a detailed introduction to the traditional BEM is given and several popular algorithms are introduced or proposed to enhance the performance of BEM. Numerical examples in heat conduction analysis, thermoelastic analysis and thermoelastic fracture problems are performed to assess the efficiency and correction of the algorithms. In addition, necessary theoretical derivations are embraced for establishing novel boundary integral equations (BIEs) for specific engineering problems. The following three parts are the main content of this thesis. (1) The first part (Part II consisting of two chapters) is aimed at heat conduction analysis by BEM. The coefficient matrix of equations formed by BEM in solving problems is fully-populated which occupy large computer memory. To deal with that, the fast multipole method (FMM) is introduced to energize the line integration boundary element method (LIBEM) to performs better in efficiency. In addition, to compute domain integrals with known or unknown integrand functions which are caused by heat sources or heterogeneity, a novel BEM, the adaptive orthogonal interpolation moving least squares (AOIMLS) method enhanced LIBEM, which also inherits the advantage of boundary-only discretisation, is proposed. Unlike LIBEM, which is an accurate and stable method for computing domain integrals, but only works when the mathematical expression of integral function in domain integrals is known, the AOIMLS enhanced LIBEM can compute domain integrals with known or unknown integral functions, which ensures all the nonlinear and nonhomogeneous problems can be solved without domain discretisation. In addition, the AOIMLS can adaptively avoid singular or ill-conditioned moment matrices, thus ensuring the stability of the calculation results. (2) In the second part (Part III consisting of four chapters), the thermoelastic problems and fracture problems are the main objectives. Due to considering thermal loads, domain integrals appear in the BIEs of the thermoelastic problems, and the expression of integrand functions is known or not depending on the temperature distribution given or not, the AOIMLS enhanced LIBEM is introduced to conduct thermoelasticity analysis thereby. Besides, a series of novel unified boundary integral equations based on BEM and DDM are derived for solving fracture problems and thermoelastic fracture problems in finite and infinite domains. Two sets of unified BIEs are derived for fracture problems in finite and infinite domains based on the direct BEM and DDM respectively, which can provide accurate and stable results. Another two sets of BIEs are addressed by employing indirect BEM and DDM, which cannot ensure a stable result, thereby a modified indirect BEM is proposed which performs much more stable. Moreover, a set of novel BIEs based on the direct BEM and DDM for cracked domains under thermal stress is proposed. (3) In the third part (Part IV consisting of one chapter), a high-efficiency combined BEM and discrete element method (DEM) is proposed to compute the inner stress distribution and particle breakage of particle assemblies based on the solution mapping scheme. For the stress field computation of particles with similar geometry, a template particle is used as the representative particle, so that only the related coefficient matrices of one template particle in the local coordinate system are needed to be calculated, while the coefficient matrices of the other particles, can be obtained by mapping between the local and global coordinate systems. Thus, the combined BEM and DEM is much more effective when modelling a large-scale particle system with a small number of distinct possible particle shapes. Furthermore, with the help of the Hoek-Brown criterion, the possible cracks or breakage paths of a particle can be obtained

Transient bioheat transfer analysis in biological tissues by fundamental-solution-based numerical methods

Taylor's expansion approach was applied to linearize the nonlinear term in the original nonlinear bioheat transfer governing equation. Then the DRM and the MFS was established to obtain the particular and homogeneous solutions. The influence of blood perfusion rate on temperature distribution in the skin tissue was analysed by changing the coefficients in the three expressions of blood perfusion rate. Numerical results showed that the variation of blood perfusion rate plays a significant role in the temperature distribution within the skin tissue. Finally, a meshless numerical scheme combining the operator splitting method (OSM), the RBF interpolation and the MFS was developed for solving transient nonlinear bioheat problems in two-dimensional skin tissue. In the numerical scheme, the nonlinearity caused by the temperature-dependent blood perfusion rate (TDBPR) is taken into consideration. In the procedure, the OSM is used to separate the Laplacian operator and the nonlinear source term, and then second-order time-stepping schemes are employed for approximating two splitting operators in order to convert the original governing equation into a linear nonhomogeneous Helmholtz-type governing equation (NHGE) at each time step. The full fields consisting of the particular and homogeneous solutions are enforced to fit the NHGE at interpolation points and the boundary conditions at boundary collocations to determine unknowns at each time step. The proposed method was verified by comparison with other methods. Furthermore, the sensitivity of the coefficients in cases of a linear and an exponential relationship of TDBPR was investigated to reveal their bioheat effect on the skin tissue

Impedance Transformers

Non