Complex variable element-free Galerkin method for viscoelasticity problems

Based on the complex variable moving least-square (CVMLS) approximation, the complex variable element-free Galerkin (CVEFG) method for two-dimensional viscoelasticity problems under the creep condition is presented in this paper. The Galerkin weak form is employed to obtain the equation system, and the penalty method is used to apply the essential boundary conditions, then the corresponding formulae of the CVEFG method for two-dimensional viscoelasticity problems under the creep condition are obtained. Compared with the element-free Galerkin (EFG) method, with the same node distribution, the CVEFG method has higher precision, and to obtain the similar precision, the CVEFG method has greater computational efficiency. Some numerical examples are given to demonstrate the validity and the efficiency of the method.     

势问题的无单元Galerkin方法的误差估计

在高维情况下,首先研究了无单元Galerkin方法的形函数构造方法——移动最小二乘法在Sobolev空间W^k,p(Ω)中的误差估计.然后,在势问题的无单元Galerkin方法的基础上,研究了势问题的通过罚函数法施加本质边界条件的无单元Galerkin方法在Sobolev空间中的误差估计.当节点和形函数满足一定条件时,证明了该误差估计是最优阶的.从误差分析中可以看出,数值解的误差与权函数的影响半径密切相关.最后,通过算例验证了结论的正确性. 关键词： 无网格方法 无单元Galerkin方法 势问题 误差估计     

Analysis of variable coefficient advection-diffusion problems via complex variable reproducing kernel particle method

The complex variable reproducing kernel particle method （CVRKPM） of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method （RKPM） and the element- free Galerkin （EFG） method.     

黏弹性问题的改进的复变量无单元Galerkin方法

基于改进的复变量移动最小二乘法,提出了二维黏弹性问题的改进的复变量无单元Galerkin方法.采用改进的复变量移动最小二乘法建立形函数,根据Galerkin积分弱形式建立求解方程,并用罚函数法施加本质边界条件,推导了二维黏弹性问题的改进的复变量无单元Galerkin方法的计算公式.最后,通过实际算例,将计算结果与复变量无单元Galerkin方法及有限元法的结果进行了对比,说明了本文方法具有更高的计算精度和计算效率.     

The element-free Galerkin method of numerically solving a regularized long-wave equation

The element-free Galerkin (EFG) method is used in this paper to find the numerical solution to a regularized long-wave (RLW) equation. The Galerkin weak form is adopted to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. The effectiveness of the EFG method of solving the RLW equation is investigated by two numerical examples in this paper.     

An improved interpolating element-free Galerkin method with nonsingular weight function for two-dimensional potential problems

In this paper, an improved interpolating moving least-square (IIMLS) method is presented. The shape function of the IIMLS method satisfies the property of Kronecker δ function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. And the number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square (MLS) approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin (IIEFG) method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin (EFG) method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has a higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.     

An improved interpolating element-free Galerkin method with a nonsingular weight function for two-dimensional potential problems

In this paper, an improved interpolating moving least-square (IIMLS) method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker δ function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square (MLS) approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin (IIEFG) method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin (EFG) method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.     

Element-free Galerkin （EFG） method for a kind of two-dimensional linear hyperbolic equation

The present paper deals with the numerical solution of a two-dimensional linear hyperbolic equation by using the element-free Galerkin (EFG) method which is based on the moving least-square approximation for the test and trial functions. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for hyperbolic problems needs only the scattered nodes instead of meshing the domain of the problem. It neither requires any element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for two-dimensional hyperbolic problems is investigated by two numerical examples in this paper.     

A new complex variable element-free Galerkin method for two-dimensional potential problems

In this paper, based on the element-free Galerkin (EFG) method and the improved complex variable moving least- square (ICVMLS) approximation, a new meshless method, which is the improved complex variable element-free Galerkin (ICVEFG) method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square (CVMLS) approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.     

一类各项异性半线性椭圆方程自然边界元与有限元耦合法

将冯康和余德浩提出的自然边界归化方法用于研究一类半线性椭圆方程外区域问题，提出一种自然边界元与有限元的耦合算法、针对某一类半线性椭圆方程，应用变分原理，研究其弱解性及Galerkin逼近，得到有限元解的误差估计及收敛阶O(h^n)，最后给出相应数值例子。     

A meshless method for the compound KdV-Burgers equation

The element-free Galerkin (EFG) method for numerically solving the compound Korteweg-de Vries-Burgers (KdVB) equation is discussed in this paper.The Galerkin weak form is used to obtain the discrete equation and the essential boundary conditions are enforced by the penalty method.The effectiveness of the EFG method of solving the compound Korteweg-de Vries-Burgers (KdVB) equation is illustrated by three numerical examples.     

Element-free Galerkin (EFG) method for analysis of the time-fractional partial differential equations

The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin (EFG) method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α (0<α ≤1). The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.     

The dimension splitting element-free Galerkin method for 3D transient heat conduction problems

By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IEFG) method is used for 2D transient heat conduction problems, and the finite difference method is applied in the splitting direction. The discretized system equation is obtained based on the Galerkin weak form of 2D problem; the essential boundary conditions are imposed with the penalty method; and the finite difference method is employed in the time domain. Four exemplary problems are chosen to verify the efficiency of the DSEFG method. The numerical solutions show that the efficiency and precision of the DSEFG method are greater than ones of the IEFG method for 3D problems.     

Element-free Galerkin method for a kind of KdV equation

The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the element-free Galerkin (EFG) method which is based on the moving least-squares approximation. A variational method is used to obtain discrete equations, and the essential boundary conditions are enforced by the penalty method. Compared with numerical methods based on mesh, the EFG method for KdV equations needs only scattered nodes instead of meshing the domain of the problem. It does not require any element connectivity and does not suffer much degradation in accuracy when nodal arrangements are very irregular. The effectiveness of the EFG method for the KdV equation is investigated by two numerical examples in this paper.     

An element-free Galerkin (EFG) method for generalized Fisher equations (GFE)

A generalized Fisher equation(GFE) relates the time derivative of the average of the intrinsic rate of growth to its variance.The exact mathematical result of the GFE has been widely used in population dynamics and genetics,where it originated.Many researchers have studied the numerical solutions of the GFE,up to now.In this paper,we introduce an element-free Galerkin(EFG) method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics.Compared with other numerical methods,the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem.The Galerkin weak form is used to obtain the discrete equations,and the essential boundary conditions are enforced by the penalty method.In comparison with the traditional method,numerical solutions show that the new method has higher accuracy and better convergence.Several numerical examples are presented to demonstrate the effectiveness of the method.     

Meshless analysis of an improved element-free Galerkin method for linear and nonlinear elliptic problems

We first give a stabilized improved moving least squares(IMLS) approximation, which has better computational stability and precision than the IMLS approximation. Then, analysis of the improved element-free Galerkin method is provided theoretically for both linear and nonlinear elliptic boundary value problems. Finally, numerical examples are given to verify the theoretical analysis.     

光滑粒子流体动力学方法固壁处理的一种新型排斥力模型

与传统网格法相比, 光滑粒子流体动力学方法不能直接施加壁面边界条件, 这就限制了该方法在工程中的应用.为此, 本文基于Galerkin加权余量法并结合传统排斥力方法, 推导出一种新的排斥力公式来施加壁面边界条件.该方法不含未知参数, 能在不减小边界粒子尺寸的情形下有效地防止流体粒子穿透壁面, 同时可避免邻近边界的流体粒子的速度及压力振荡. 分别通过静止液柱算例、液柱坍塌算例、容器中液体静止算例及溃坝算 例来验证本文方法的有效性, 并与传统边界处理方法进行对比, 结果表明: 本文方法克服了传统方法存在的缺陷, 是一种有效的固壁边界处理方法. 关键词： 光滑粒子流体动力学法 固壁边界 排斥力 加权余量法     

An element-free Galerkin(EFG) method for numerical solution of the coupled Schrdinger-KdV equations

The present paper deals with the numerical solution of the coupled Schrdinger-KdV equations using the elementfree Galerkin(EFG) method which is based on the moving least-square approximation.Instead of traditional mesh oriented methods such as the finite difference method(FDM) and the finite element method(FEM),this method needs only scattered nodes in the domain.For this scheme,a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method.In numerical experiments,the results are presented and compared with the findings of the finite element method,the radial basis functions method,and an analytical solution to confirm the good accuracy of the presented scheme.     

Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and the other based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.     

Analysis of elastoplasticity problems using an improved complex variable element-free Galerkin method

In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin(ICVEFG) method is presented for two-dimensional(2D) elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.