A generalized element-free Galerkin method for Stokes problem

A generalized element-free Galerkin (GEFG) method is developed in this paper for solving Stokes problem in primitive variable form. To obtain stable numerical results for both velocity and pressure, extended terms are only introduced into the approximate space of velocity in a special way as that in the generalized finite element method. Theoretical analysis shows that the GEFG method implies a stabilized formulation similar to that in the variational multiscale element-free Galerkin (VMEFG) method. Numerical results show the efficiency of the present method and reveal that both computational errors and CPU times of the present method are less than those of the VMEFG and the finite element methods.     

The element-free Galerkin method for the nonlinear p-Laplacian equation

The element-free Galerkin (EFG) method is developed in this paper for solving the nonlinear p-Laplacian equation. The moving least squares approximation is used to generate meshless shape functions, the penalty approach is adopted to enforce the Dirichlet boundary condition, the Galerkin weak form is employed to obtain the system of discrete equations, and two iterative procedures are developed to deal with the strong nonlinearity. Then, the computational formulas of the EFG method for the p-Laplacian equation are established. Numerical results are finally given to verify the convergence and high computational precision of the method.     

Convergence of optimal control problems governed by second kind parabolic variational inequalities

We consider a family of optimal control problems where the control variable is given by a boundary condition of Neumann type. This family is governed by parabolic variational inequalities of the second kind. We prove the strong convergence of the optimal control and state systems associated to this family to a similar optimal control problem. This work solves the open problem left by the authors in IFIP TC7 CSMO2011.     

势问题的复变量无单元Galerkin方法

为提高无单元Galerkin(Element-Free Galerkin, EFG)方法的计算效率,将复变量移动最小二乘法与EFG方法结合,利用控制方程的积分弱形式并采用Lagrange乘子法引入边界条件,提出势问题的复变量无单元Galerkin(Complex Variable EFG,CVEFG)方法,并推导相关公式.与传统的EFG方法相比,该方法采用复变量移动最小二乘法可以减少试函数中的待定系数,从而减少计算量、提高计算效率. 最后,给出数值算例验证该方法的有效性.     

A meshless Galerkin method for Stokes problems using boundary integral equations

A meshless Galerkin scheme for the simulation of two-dimensional incompressible viscous fluid flows in primitive variables is described in this paper. This method combines a boundary integral formulation for the Stokes equation with the moving least-squares (MLS) approximations for construction of trial and test functions for Galerkin approximations. Unlike the domain-type method, this scheme requires only a nodal structure on the bounding surface of a body for approximation of boundary unknowns, thus it is especially suitable for the exterior problems. Compared to other meshless methods such as the boundary node method and the element free Galerkin method, in which the MLS is also introduced, boundary conditions do not present any difficulty in using this meshless method. The convergence and error estimates of this approach are presented. Numerical examples are also given to show the efficiency of the method.     

Penalty method for variational inequalities with multivalued mappings. Part I

The regularization method of variational inequalities is generalized by means of penalty operators to a class of variational inequalities with multivalued mappings. Translated from Kibernetika i Sistemnyi Analiz, No. 4, pp. 57–69, July–August, 2000.     

Error analysis of a two-grid discontinuous Galerkin method for non-linear parabolic equations

Discontinuous Galerkin (DG) approximations for non-linear parabolic problems are investigated. To linearize the discretized equations, we use a two-grid method involving a small non-linear system on a coarse gird of size H and a linear system on a fine grid of size h. Error estimates in H¹-norm are obtained, O(h^r+H^r+1) where r is the order of the DG space. The analysis shows that our two-grid DG algorithm will achieve asymptotically optimal approximation as long as the mesh sizes satisfy h=O(H^(r+1)/r). The numerical experiments verify the efficiency of our algorithm.     

On a Galerkin boundary node method for potential problems

A Galerkin boundary node method (GBNM), for boundary only analysis of partial differential equations, is discussed in this paper. The GBNM combines an equivalent variational form of a boundary integral equation with the moving least-squares (MLS) approximations for generating the trial and test functions of the variational formulation. In this approach, only a nodal data structure on the boundary of a domain is required, and boundary conditions can be implemented directly and easily despite of the fact that the MLS shape functions lack the delta function property. Formulations of the GBNM using boundary singular integral equations of the second kind for potential problems are developed. The theoretical analysis and numerical results indicate that it is an efficient and accurate numerical method.     

A three-step wavelet Galerkin method for parabolic and hyperbolic partial differential equations

A three-step wavelet Galerkin method based on Taylor series expansion in time is proposed. The scheme is third-order accurate in time and O(2^?jp ) accurate in space. Unlike Taylor–Galerkin methods, the present scheme does not contain any new higher-order derivatives which makes it suitable for solving non-linear problems. The compactly supported orthogonal wavelet bases D6 developed by Daubechies are used in the Galerkin scheme. The proposed scheme is tested with both parabolic and hyperbolic partial differential equations. The numerical results indicate the versatility and effectiveness of the proposed scheme.     

E(FG)<Superscript>2</Superscript>: a new fixed-grid shape optimization method based on the element-free galerkin mesh-free analysis

We propose a shape optimization method over a fixed grid. Nodes at the intersection with the fixed grid lines track the domain’s boundary. These “floating” boundary nodes are the only ones that can move/appear/disappear in the optimization process. The element-free Galerkin (EFG) method, used for the analysis problem, provides a simple way to create these nodes. The fixed grid (FG) defines integration cells for EFG method. We project the physical domain onto the FG and numerical integration is performed over partially cut cells. The integration procedure converges quadratically. The performance of the method is shown with examples from shape optimization of thermal systems involving large shape changes between iterations. The method is applicable, without change, to shape optimization problems in elasticity, etc. and appears to eliminate non-differentiability of the objective noticed in finite element method (FEM)-based fictitious domain shape optimization methods. We give arguments to support this statement. A mathematical proof is needed.     

Differential quadrature domain decomposition method for a class of parabolic equations

In this paper the differential quadrature method (DQM) and the domain decomposition method (DDM) are combined to form the differential quadrature domain decomposition method (DQDDM), in which the boundary reduction technique (BRM) is adopted. The DQDDM is applied to a class of parabolic equations, which have discontinuity in the coefficients of the equation, or weak discontinuity in the initial value condition. Two numerical examples belonging to this class are computed. It is found that the application of this method to the above mentioned problems is seen to lead to accurate results with relatively small computational effort.     

Jackson-type inequalities for h-p clouds and error estimates

The interest in meshfree methods for solving boundary-value problems has grown rapidly in recent years. A meshless method that has attracted much interest in the community of computational mechanics is the h-p clouds method. For this kind of applications it is fundamental to analyze the orders of approximation. In this paper we prove Jackson-type inequalities for h-p cloud functions. These inequalities set up a general framework for the theoretical analysis of high order error estimates of the h-p clouds method, with the same remarkable features of finite element theory.     

Nonzero solutions for a class of set-valued variational inequalities in reflexive Banach spaces

In this paper, we study the existence of nonzero solutions for a class of set-valued variational inequalities involving set-contractive mappings by using the fixed point index approach in reflexive Banach spaces. Some new existence theorems of nonzero solutions for this class of set-valued variational inequalities are established.     

Finite time blow-up for a class of parabolic or pseudo-parabolic equations

In this paper, we study the initial boundary value problem for a class of parabolic or pseudo-parabolic equations:

$u_t?a\Delta u_t?\Delta u+bu=k\left(t\right)|u{|}^{p?2}u,(x,t)\in \Omega \times (0,T),$

where $a\ge 0$, $b>?{?}_1$ with ${?}_1$ being the principal eigenvalue for $?\Delta$ on $H_0^1\left(\Omega \right)$ and $k\left(t\right)>0$. By using the potential well method, Levine’s concavity method and some differential inequality techniques, we obtain the finite time blow-up results provided that the initial energy satisfies three conditions: (i) $J(u_0;0)<0$; (ii) $J(u_0;0)\le d\left(\infty \right)$, where $d\left(\infty \right)$ is a nonnegative constant; (iii) $0<J(u_0;0)\le C\rho \left(0\right)$, where $\rho \left(0\right)$ involves the $L^2$-norm or $H_0^1$-norm of the initial data. We also establish the lower and upper bounds for the blow-up time. In particular, we obtain the existence of certain solutions blowing up in finite time with initial data at the Nehari manifold or at arbitrary energy level.     

A priori and a posteriori error estimates for the method of lumped masses for parabolic optimal control problems

In this paper, we examine the method of lumped masses for the approximation of convex optimal control problems governed by linear parabolic equations, where the lumped mass method is used for the discretization of the state equation. We derive some a priori and a posteriori error estimates for both the state and control approximations with control constraints of obstacle type. Numerical experiments are given to show the efficiency and reliability of the lumped mass method.     

A class of projection and contraction methods for asymmetric linear variational inequalities and their relations to Fukushima's descent method

For solving asymmetric linear variational inequalities, we present a class of projection and contraction methods under the general G-norm. The search direction of our methods is just a convex combination of two descent directions of Fukushima's merit function. However, we use the direction to reduce the distance function (1/2)u − u*²_G, where μ* is a solution point of the problem. Finally, we report some numerical results for spatial price equilibrium problems by using the presented methods.     

An iterative approximation scheme for solving a split generalized equilibrium,variational inequalities and fixed point problems

In this paper, we consider a common solution of three problems in Hilbert spaces including the split generalized equilibrium problem, the variational inequality problem and fixed point problem. For finding the solution, we present a new iterative method and prove the strongly convergence theorem under mild conditions. Moreover, some numerical examples are given in the last section.