The Method of Polynomial Particular Solutions for Solving Nonlinear Poisson-Type Equations

In this paper,the method of polynomial particular solutions is used to solve nonlinear Poisson-type partial differential equations in one,two,and three dimensions.The condition number of the coefficient matrix is reduced through the implementation of multiple scale technique,ultimately yielding a stable numerical solution.The methodological process can be divided into two main parts:first,identifying the corresponding polynomial particular solutions for the linear differential operator terms in the governing equations,and second,employing these polynomial particular solutions as basis function to iteratively solve the remaining nonlinear terms within the governing equations.Additionally,we investigate the potential improvement in numerical accuracy for equations with singularities in the analytical solution by shifting the computational domain a certain distance.Numerical experiments are conducted to assess both the accuracy and stability of the proposed method.A comparison of the obtained results with those produced by other numerical methods demonstrates the accuracy,stability,and efficiency of the proposed method in handling nonlinear Poisson-type partial differential equations.

Numerical Simulation of Non-Linear Schrodinger Equations in Arbitrary Domain by the Localized Method of Approximate Particular Solution

The aim of this paper is to propose a fast meshless numerical scheme for the simulation of non-linear Schrodinger equations.In the proposed scheme,the implicit-Euler scheme is used for the temporal discretization and the localized method of approximate particular solution(LMAPS)is utilized for the spatial discretization.The multiple-scale technique is introduced to obtain the shape parameters of the multiquadric radial basis function for 2D problems and the Gaussian radial basis function for 3D problems.Six numerical examples are carried out to verify the accuracy and efficiency of the proposed scheme.Compared with well-known techniques,numerical results illustrate that the proposed scheme is of merits being easy-to-program,high accuracy,and rapid convergence even for long-term problems.These results also indicate that the proposed scheme has great potential in large scale problems and real-world applications.

A Comparative Study on Polynomial Expansion Method and Polynomial Method of Particular Solutions

In this study,the polynomial expansion method(PEM)and the polynomial method of particular solutions(PMPS)are applied to solve a class of linear elliptic partial differential equations(PDEs)in two dimensions with constant coefficients.In the solution procedure,the sought solution is approximated by the Pascal polynomials and their particular solutions for the PEM and PMPS,respectively.The multiple-scale technique is applied to improve the conditioning of the resulted linear equations and the accuracy of numerical results for both of the PEM and PMPS.Some mathematical statements are provided to demonstrate the equivalence of the PEM and PMPS bases as they are both bases of a certain polynomial vector space.Then,some numerical experiments were conducted to validate the implementation of the PEM and PMPS.Numerical results demonstrated that the PEM is more accurate and well-conditioned than the PMPS and the multiple-scale technique is essential in these polynomial methods.

A Revisit on the Derivation of the Particular Solution for the Differential Operator ∆^(2) ± λ^(2)

In this paper,we applied the polyharmonic splines as the basis functions to derive particular solutions for the differential operator ∆^(2) ± λ^(2).Similar to the derivation of fundamental solutions,it is non-trivial to derive particular solutions for higher order differential operators.In this paper,we provide a simple algebraic factorization approach to derive particular solutions for these types of differential operators in 2D and 3D.The main focus of this paper is its simplicity in the sense that minimal mathematical background is required for numerically solving higher order partial differential equations such as thin plate vibration.Three numerical examples in both 2D and 3D are given to validate particular solutions we derived.

The Recursive Formulation of Particular Solutions for Some Elliptic PDEs with Polynomial Source Functions

In this paper we develop an efficient meshless method for solving inhomogeneous elliptic partial differential equations.We first approximate the source function by Chebyshev polynomials.We then focus on how to find a polynomial particular solution when the source function is a polynomial.Through the principle of the method of undetermined coefficients and a proper arrangement of the terms for the polynomial particular solution to be determined,the coefficients of the particular solution satisfy a triangular system of linear algebraic equations.Explicit recursive formulas for the coefficients of the particular solutions are derived for different types of elliptic PDEs.The method is further incorporated into the method of fundamental solutions for solving inhomogeneous elliptic PDEs.Numerical results show that our approach is efficient and accurate.

Efficient Algorithms for Approximating Particular Solutions of Elliptic Equations Using Chebyshev Polynomials

In this paper,we propose efficient algorithms for approximating particular solutions of second and fourth order elliptic equations.The approximation of the particular solution by a truncated series of Chebyshev polynomials and the satisfaction of the differential equation lead to upper triangular block systems,each block being an upper triangular system.These systems can be solved efficiently by standard techniques.Several numerical examples are presented for each case.

Position optimization of particular solution sources for distributed source boundary point method by volume velocity-matching

Choosing particular solution source and its position have great influence on accu- racy of sound field prediction in distributed source boundary point method. An optimization method for determining the position of particular solution sources is proposed to get high accu- racy prediction result. In this method, tripole is chosen as the particular solution. The upper limit frequency of calculation is predicted by setting 1% volume velocity relative error limit using vibration velocity of structure surface. Then, the optimal position of particular solution sources, in which the relative error of volume velocity is minimum, is determined within the range of upper limit frequency by searching algorithm using volume velocity matching. The transfer matrix between pressure and surface volume velocity is constructed in the optimal position. After that, the sound radiation of structure is calculated by the matrix. The results of numerical simulation show that the calculation error is significantly reduced by the proposed method. When there are vibration velocity measurement errors, the calculation errors can be controlled within 5% by the method.

BEM FOR MODAL ANALYSIS OF 3-D ANISOTROPIC STRUCTURES

The boundary element method for the modal analysis of freevibration for 3-D anisotropic structures using particular solutionshas been developed. The complete polynomials of order two are used toconstruct the particular solutions for general anisotropic materials.The numerical results for 3-D free vibra- tion analysis of anisotropic cantilever beam by the method presented is in goodagreement with the results us- ing the Ritz technique. Foranisotropic materials, the numerical results calculated form theproposed method are in good agreement with the results from MSC.NASTRAN.

Bending analysis of magnetoelectroelastic nanoplates resting on Pasternak elastic foundation based on nonlocal theory

Based on the nonlocal theory and Mindlin plate theory,the governing equations(i.e.,a system of partial differential equations(PDEs)for bending problem)of magnetoelectroelastic(MEE)nanoplates resting on the Pasternak elastic foundation are first derived by the variational principle.The polynomial particular solutions corresponding to the established model are then obtained and further employed as basis functions with the method of particular solutions(MPS)to solve the governing equations numerically.It is confirmed that for the present bending model,the new solution strategy possesses more general applicability and superior flexibility in the selection of collocation points.The effects of different boundary conditions,applied loads,and geometrical shapes on the bending properties of MEE nanoplates are evaluated by using the developed method.Some important conclusions are drawn,which should be helpful for the design and applications of electromagnetic nanoplate structures.

STRESS RATE INTEGRAL EQUATIONS OF ELASTOPLASTICITY

The stress rate integral equations of elastoplasticity are deduced based on Ref. [1] by consistent methods. The point at which the stresses and/or displacements are calculated can be in the body or on the boundary, and in the plastic region or elastic one. The existence of the principal value integral in the plastic region is demonstrated strictly, and the theoretical basis is presented for the paticular solution method by unit initial stress fields. In the present method, programming is easy and general, and the numerical results are excellent.

BOUNDARY ELEMENT METHOD FOR MODEL ANALYSIS OF 2-D COMPOSITE STRUCTURE

The boundary element method is used for he modal analysis of freevibration of 2-D composite structures in this paper. Since theparticular solution method is used to treat the terms of body forces(inertial forces) in the equation of motion, only static fundamentalsolutions are needed in solving the problem. For an isotropiccantilever beam, the numerical results obtained by using the BEMpresented in this paper are in good agreement, with those of usingFEM or other BEM, but this BEM can also be used to analyze problemsfor anisotropic materials.

The Method of Fundamental Solutions for Solving Convection-Diffusion Equations with Variable Coefficients

A meshless method based on the method of fundamental solutions(MFS)is proposed to solve the time-dependent partial differential equations with variable coefficients.The proposed method combines the time discretization and the onestage MFS for spatial discretization.In contrast to the traditional two-stage process,the one-stage MFS approach is capable of solving a broad spectrum of partial differential equations.The numerical implementation is simple since both closed-form approximate particular solution and fundamental solution are easy to find than the traditional approach.The numerical results show that the one-stage approach is robust and stable.

Using Gaussian Eigenfunctions to Solve Boundary Value Problems

Kernel-basedmethods are popular in computer graphics,machine learning,and statistics,among other fields;because they do not require meshing of the domain under consideration,higher dimensions and complicated domains can be managed with reasonable effort.Traditionally,the high order of accuracy associated with these methods has been tempered by ill-conditioning,which ariseswhen highly smooth kernels are used to conduct the approximation.Recent advances in representing Gaussians using eigenfunctions have proven successful at avoiding this destabilization in scattered data approximation problems.This paper will extend these techniques to the solution of boundary value problems using collocation.The method of particular solutions will also be considered for elliptic problems,using Gaussian eigenfunctions to stably produce an approximate particular solution.

A Contour Integral Method for Linear Differential Equations in Complex Plane

This paper presents the contour integral method for solving the linear constant coefficient ordinary differential equations in complex plane,and obtains the uniform expressions of the general solutions.Firstly,by using Residue Theorem,the general form of the contour integral representation for the homogeneous complex differential equation is obtained,which can be degenerated to classical results in real line.As for inhomogeneous complex differential equations with constant coefficients,we construct the integral expression of the particular solution for any continuous forcing term,and give rigorous proof via Residue Theorem.Thus the general solutions of inhomogeneous complex differential equations are also given.The main purpose of this paper is to give a foundation for a complete theory of linear complex differential equations with constant coefficients by a contour integral method.The results can not only solve the inhomogeneous complex differential equation well,but also explain the forms that are difficult to be understood in the classical solutions.

A Boundary Meshless Method for Solving Heat Transfer Problems Using the Fourier Transform

Fourier transform is applied to remove the time-dependent variable in the diffusion equation.Under non-harmonic initial conditions this gives rise to a non-homogeneous Helmholtz equation,which is solved by the method of fundamental solutions and the method of particular solutions.The particular solution of Helmholtz equation is available as shown in[4,15].The approximate solution in frequency domain is then inverted numerically using the inverse Fourier transform algorithm.Complex frequencies are used in order to avoid aliasing phenomena and to allow the computation of the static response.Two numerical examples are given to illustrate the effectiveness of the proposed approach for solving 2-D diffusion equations.