Finite integration method for partial differential equations

A finite integration method is proposed in this paper to deal with partial differential equations in which the finite integration matrices of the first order are constructed by using both standard integral algorithm and radial basis functions interpolation respectively. These matrices of first order can directly be used to obtain finite integration matrices of higher order. Combining with the Laplace transform technique, the finite integration method is extended to solve time dependent partial differential equations. The accuracy of both the finite integration method and finite difference method are demonstrated with several examples. It has been observed that the finite integration method using either radial basis function or simple linear approximation gives a much higher degree of accuracy than the traditional finite difference method.     

Analysis of a non-standard mixed finite element method with applications to superconvergence

We show that a non-standard mixed finite element method proposed by Barrios and Gatica in 2007, is a higher order perturbation of the least-squares mixed finite element method. Therefore, it is also superconvergent whenever the least-squares mixed finite element method is superconvergent. Superconvergence of the latter was earlier investigated by Brandts, Chen and Yang between 2004 and 2006. Since the new method leads to a non-symmetric system matrix, its application seems however more expensive than applying the least-squares mixed finite element method. Dedicated to Ivan Hlaváček on the occasion of his 75th birthday     

对流扩散方程的有限体积-有限元方法的误差估计

本文结合有限体积方法和有限元方法处理非线性对流扩散问题,非线性对流项利用有限体积方法处理,扩散项利用有限元方法离散,并给近似解的误差估计。     

A finite element method for granular flow through a frictional boundary

A finite element method for the flow of dry granular solids through a domain involving a frictional contact boundary is formulated. The granular material is assumed as a compressible viscous-elastic–plastic continuum. Based on the principles of continuum mechanics, a complete set of equations is developed. The resulting boundary value problem is solved by the finite element method in space and by the finite difference method in time. The derivation of the finite element equations and the mathematical framework of the numerical technique are presented, together with two illustrative examples to demonstrate the validity of the technique.     

Finite bases for finitely generated,relatively congruence distributive quasivarieties

In [21], D. Pigozzi has proved in a non-constructive way that every relatively congruence distributive quasivariety of finite type generated by a finite set of finite algebras is finitely axiomatizable. In this paper we show that the non-constructive parts of Pigozzi's argument can be replaced by constructive ones. As a result we obtain a method of constructing a finite set of quasi-equational axioms for each relatively congruence distributive quasivariety generated by a given finite set of finite algebras of finite type. The method can also be applied to finitely generated congruence distributive varieties.Presented by Joel Berman.     

An operator splitting method with FDM and FEM for the convection- diffusion equation using boundary-fitted coordinate system

An operator splitting method combining finite difference method and finite element method is proposed in this paper by using boundary-fitted coordinate system. The governing equation is split into advection and diffusion equations and solved by finite difference method using boundary-fitted coordinate system and finite element method respectively. An example for which analytic solution is available is used to verified the proposed methods and the agreement is very good. Numerical results show that it is very efficient.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

Analysis of a finite volume element method for the Stokes problem

In this paper we propose a stabilized conforming finite volume element method for the Stokes equations. On stating the convergence of the method, optimal a priori error estimates in different norms are obtained by establishing the adequate connection between the finite volume and stabilized finite element formulations. A superconvergence result is also derived by using a postprocessing projection method. In particular, the stabilization of the continuous lowest equal order pair finite volume element discretization is achieved by enriching the velocity space with local functions that do not necessarily vanish on the element boundaries. Finally, some numerical experiments that confirm the predicted behavior of the method are provided.     

A wavelet-based stochastic finite element method of thin plate bending

A wavelet-based stochastic finite element method is presented for the bending analysis of thin plates. The wavelet scaling functions of spline wavelets are selected to construct the displacement interpolation functions of a rectangular thin plate element and the displacement shape functions are expressed by the spline wavelets. A new wavelet-based finite element formulation of thin plate bending is developed by using the virtual work principle. A wavelet-based stochastic finite element method that combines the proposed wavelet-based finite element method with Monte Carlo method is further formulated. With the aid of the wavelet-based stochastic finite element method, the present paper can deal with the problem of thin plate response variability resulting from the spatial variability of the material properties when it is subjected to static loads of uncertain nature. Numerical examples of thin plate bending have demonstrated that the proposed wavelet-based stochastic finite element method can achieve a high numerical accuracy and converges fast.     

三维热传导型半导体问题的交替方向特征有限元方法

提出交替方向特征有限元方法,对电场位势方程采用混合元格式,对电子,空穴浓度方程采用交替方向特征有限元格式,对温度方程提出交替方向格式.应用向量积计算及先验估计理论和技巧,得到最佳的L2误差估计.     

Determination of the equilibrium boundaries of a two-layer system convection stability : PMM vol. 43, no. 6, 1979, pp. 1008–1013

The method which permits the determination of a large number of critical Rayleigh numbers and respective critical motions is applied to the determination of convection stability boundaries of a two-layer system. The method consists essentially of reducing the problem to the algebraic problem of eigenvalues by discretizing the equations by the method of finite elements or finite differences. Only the scheme of the method based on discretization by finite elements is then presented.     

An Iterative FEM for Solving Wave Problems of a Halfspace

An iterative finite element method for solving wave problems of a halfspace is presented in this paper. The halfspace is first truncated by introducing a fictive finite boundary on which some fictive boundary conditions must be imposed. A finite computational domain is in each iteration subjected to actual boundary conditions on real boundary and to fictive Dirichlet or Neumann boundary conditions on the fictive boundary. The radiation condition is satisfied by using DtN operator. The DtN operator is not introduce in the finite element formulation on the fictive boundary so any finite elements can be used. The method is simple and specially useful for computing higher harmonics.     

三维热传导型半导体问题的特征混合元方法和分析

本文研究三维热传导型半导体态问题的特征混合元方法及其理论分析,其数学模型是一类非线性偏微分方程的初边值问题,对电子位势方程提出混合元逼近,对电子,空穴浓度方程笔挺表限元逼近;对热传导方程采用对时间向后差分的Ｇａｌｅｒｋｉｎ逼近,应用微分方程先验估计理论和技巧得到了最优阶Ｌ＾２误差估计。     

Non-overlapping Domain Decomposition Method for a Nodal Finite Element Method

A new approach is proposed for constructing nonoverlapping domain decomposition procedures for solving a linear system related to a nodal finite element method. It applies to problems involving either positive semi-definite or complex indefinite local matrices. The main feature of the method is to preserve the continuity requirements on the unknowns and the finite element equations at the nodes shared by more than two subdomains and to suitably augment the local matrices. We prove that the corresponding algorithm can be seen as a converging iterative method for solving the finite element system and that it cannot break down. Each iteration is obtained by solving uncoupled local finite element systems posed in each subdomain and, in contrast to a strict domain decomposition method, is completed by solving a linear system whose unknowns are the degrees of freedom attached to the above special nodes.     

Numerical treatment for solving fractional Riccati differential equation

This paper presents an accurate numerical method for solving fractional Riccati differential equation (FRDE). The proposed method so called fractional Chebyshev finite difference method (FCheb-FDM). In this technique, we approximate FRDE with a finite dimensional problem. The method is based on the combination of the useful properties of Chebyshev polynomials approximation and finite difference method. The Caputo fractional derivative is replaced by a difference quotient and the integral by a finite sum. By this method the given problem is reduced to a problem for solving a system of algebraic equations, and by solving this system, we obtain the solution of FRDE. Special attention is given to study the convergence analysis and estimate an error upper bound of the obtained approximate formula. Illustrative examples are included to demonstrate the validity and applicability of the proposed technique.     

不可压缩流问题低次元稳定有限体积数值方法研究

分析了R^d,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈Hd,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈H¹情况下,本文方法得到的解与稳定有限元方法解之间具有O(h1情况下,本文方法得到的解与稳定有限元方法解之间具有O(h²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

On finite difference approximation of a matrix-vector product in the Jacobian-free Newton-Krylov method

The Jacobian-free Newton-Krylov (JFNK) method is a special kind of Newton-Krylov algorithm, in which the matrix-vector product is approximated by a finite difference scheme. Consequently, it is not necessary to form and store the Jacobian matrix. This can greatly improve the efficiency and enlarge the application area of the Newton-Krylov method. The finite difference scheme has a strong influence on the accuracy and robustness of the JFNK method. In this paper, several methods for approximating the Jacobian-vector product, including the finite difference scheme and the finite difference step size, are analyzed and compared. Numerical results are given to verify the effectiveness of different finite difference methods.     

A finite element method for interface problems in domains with smooth boundaries and interfaces

This paper is concerned with the analysis of a finite element method for nonhomogeneous second order elliptic interface problems on smooth domains. The method consists in approximating the domains by polygonal domains, transferring the boundary data in a natural way, and then applying a finite element method to the perturbed problem on the approximate polygonal domains. It is shown that the error in the finite element approximation is of optimal order for linear elements on a quasiuniform triangulation. As such the method is robust in the regularity of the data in the original problem.     

Compact difference scheme for two-dimensional fourth-order nonlinear hyperbolic equation

High-order compact finite difference method for solving the two-dimensional fourth-order nonlinear hyperbolic equation is considered in this article. In order to design an implicit compact finite difference scheme, the fourth-order equation is written as a system of two second-order equations by introducing the second-order spatial derivative as a new variable. The second-order spatial derivatives are approximated by the compact finite difference operators to obtain a fourth-order convergence. As well as, the second-order time derivative is approximated by the central difference method. Then, existence and uniqueness of numerical solution is given. The stability and convergence of the compact finite difference scheme are proved by the energy method. Numerical results are provided to verify the accuracy and efficiency of this scheme.     

On the semi-discrete stabilized finite volume method for the transient Navier–Stokes equations

A stabilized finite volume method for solving the transient Navier–Stokes equations is developed and studied in this paper. This method maintains conservation property associated with the Navier–Stokes equations. An error analysis based on the variational formulation of the corresponding finite volume method is first introduced to obtain optimal error estimates for velocity and pressure. This error analysis shows that the present stabilized finite volume method provides an approximate solution with the same convergence rate as that provided by the stabilized linear finite element method for the Navier–Stokes equations under the same regularity assumption on the exact solution and a slightly additional regularity on the source term. The stability and convergence results of the proposed method are also demonstrated by the numerical experiments presented.