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

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

弹性力学的无单元Galerkin方法的误差估计

基于移动最小二乘法在Sobolev空间W^k,p(Ω)中的误差估计以及弹性力学问题的变分弱形式中出现的双线性形式的连续性和强制性,研究了弹性力学问题的无单元Galerkin方法的误差分析以及数值解的误差和影响域半径之间的关系,给出了弹性力学问题的无单元Galerkin方法在Sobolev空间中的误差估计定理,并证明了当节点和形函数满足一定条件时该误差估计是最优阶的.从误差分析中可以看出,数值解的误差与权函数的影响域半径密切相关.最后,通过算例验证了结论的正确性. 关键词： 无网格方法 无单元Galerkin方法 弹性力学 误差估计     

积分节点迎风偏移的节点积分无单元Galerkin方法

建立求解稳态对流-扩散方程的-种稳定、高效的无单元Galerkin方法.该方法计算积分时采用基于局部Taylor展开的节点积分,并根据对流占优的程度对积分节点进行自适应迎风偏移.与传统的使用稳定化的无单元Galerkin方法相比,该方法是-种不依赖于背景网格积分的纯无网格方法,具有更好的稳定性和较高的计算效率,其程序实施更为简便.     

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

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

运用改进的无单元Galerkin方法分析机场道面断裂力学问题

在改进的无单元Galerkin方法的基础上,将能反映裂纹尖端附近应力奇异性的特征项r~(1/2)引入改进的移动最小二乘法的基函数中,将断裂力学和改进的无单元Galerkin方法结合,研究了线弹性断裂力学的改进的无单元Galerkin方法,并对含反射裂缝的机场复合道面层状体系结构进行了数值分析.本文的理论为机场复合道面断裂力学分析提供了一种新方法.     

耦合径向基函数与多项式基函数的无网格方法

耦合径向基函数和多项式基函数,形成一种新的近似函数.该近似函数对散乱分布的离散数据点进行逼近时,只需节点信息,不需要划分网格.详细描述了耦合近似函数的建立、属性、插值行为及其形函数和形函数导数的性质.最后引入修正变分原理和单位分解积分技术求解边值问题,并给出了计算实例,表明耦合径向基函数和多项式基函数是一种非常有效的方法.     

Runge-Kutta间断Galerkin法在求解Navier-Stokes方程中的应用

Cockburn & Shu^[1] 在1988年提出了一种TVB Runge-Kuta局部投影的间断Galerkin有限元方法应用于Euler方程的求解,并取得了成功。文章将该方法进一步应用到Navier-Stokes方程的求解。     

基于广义交替数值通量的局部间断Galerkin方法求解二维波动方程

波的传播往往在复杂的地质结构中进行,如何有效地求解非均匀介质中的波动方程一直是研究的热点.本文将局部间断Galekin(local discontinuous Galerkin, LDG)方法引入到数值求解波动方程中.首先引入辅助变量,将二阶波动方程写成一阶偏微分方程组,然后对相应的线性化波动方程和伴随方程构造间断Galerkin格式;为了保证离散格式满足能量守恒,在单元边界上选取广义交替数值通量,理论证明该方法满足能量守恒性.在时间离散上,采用指数积分因子方法,为了提高计算效率,应用Krylov子空间方法近似指数矩阵与向量的乘积.数值实验中给出了带有精确解的算例,验证了LDG方法的数值精度和能量守恒性;此外,也考虑了非均匀介质和复杂计算区域的计算,结果表明LDG方法适合模拟具有复杂结构和多尺度结构介质中的传播.     

非线性Schrdinger方程的直接间断Galerkin方法

讨论一维和二维非线性Schrdinger(NLS)方程的数值求解.基于扩散广义黎曼问题的数值流量,构造一种直接间断Galerkin方法(DDG)求解非线性Schrdinger方程.证明该方法L2稳定性,并说明DDG格式是一种守恒的数值格式.对一维NLS方程的计算表明,DDG格式能够模拟各种孤立子形态,而且可以保持长时间的高精度.二维NLS方程的数值结果显示该方法的高精度和捕捉大梯度的能力.     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

An improved element-free Galerkin method for solving generalized fifth-order Korteweg-de Vries equation

In this paper, an improved element-free Galerkin (IEFG) method is proposed to solve the generalized fifth-order Korteweg-de Vries (gfKdV) equation. When the traditional element-free Galerkin (EFG) method is used to solve such an equation, unstable or even wrong numerical solutions may be obtained due to the violation of the consistency conditions of the moving least-squares (MLS) shape functions. To solve this problem, the EFG method is improved by employing the improved moving least-squares (IMLS) approximation based on the shifted polynomial basis functions. The effectiveness of the IEFG method for the gfKdV equation is investigated by using some numerical examples. Meanwhile, the motion of single solitary wave and the interaction of two solitons are simulated using the IEFG method.     

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.     

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.     

Solving unsteady Schr?dinger equation using the improved element-free Galerkin method

By employing the improved moving least-square (IMLS) approximation, the improved element-free Galerkin (IEFG) method is presented for the unsteady Schrödinger equation. In the IEFG method, the two-dimensional (2D) trial function is approximated by the IMLS approximation, the variation method is used to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. Because the number of coefficients in the IMLS approximation is less than in the moving least-square (MLS) approximation, fewer nodes are needed in the entire domain when the IMLS approximation is used than when the MLS approximation is adopted. Then the IEFG method has high computational efficiency and accuracy. Several numerical examples are given to verify the accuracy and efficiency of the IEFG method in this paper.     

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.     

Analysis of the generalized Camassa and Holm equation with the improved element-free Galerkin method

In this paper, we analyze the generalized Camassa and Holm （CH） equation by the improved element-free Galerkin （IEFG） method. By employing the improved moving least-square （IMLS） approximation, we derive the formulas for the generalized CH equation with the IEFG method. A variational method is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Because there are fewer coefficients in the IMLS approximation than in the MLS approximation, and in the IEFG method, fewer nodes are selected in the entire domain than in the conventional EFG method, the IEFG method should result in a higher computing speed. The effectiveness of the IEFG method for the generalized CH equation is investigated by 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.     

Topology optimization using the improved element-free Galerkin method for elasticity

The improved element-free Galerkin(IEFG) method of elasticity is used to solve the topology optimization problems.In this method, the improved moving least-squares approximation is used to form the shape function. In a topology optimization process, the entire structure volume is considered as the constraint. From the solid isotropic microstructures with penalization, we select relative node density as a design variable. Then we choose the minimization of compliance to be an objective function, and compute its sensitivity with the adjoint method. The IEFG method in this paper can overcome the disadvantages of the singular matrices that sometimes appear in conventional element-free Galerkin(EFG) method. The central processing unit(CPU) time of each example is given to show that the IEFG method is more efficient than the EFG method under the same precision, and the advantage that the IEFG method does not form singular matrices is also shown.