Convergence in L space for the homogenization problems of elliptic and parabolic equations in the plane

We study the convergence rate of an asymptotic expansion for the elliptic and parabolic operators with rapidly oscillating coefficients. First we propose homogenized expansions which are convolution forms of Green function and given force term of elliptic equation. Then, using local L^p-theory, the growth rate of the perturbation of Green function is found. From the representation of elliptic solution by Green function, we estimate the convergence rate in L^p space of the homogenized expansions to the exact solution. Finally, we consider L2(0,T:H1(Ω)) or L∞(Ω×(0,T)) convergence rate of the first order approximation for parabolic homogenization problems.     

Absolutely stable explicit schemes for reaction systems.

We introduce two numerical schemes for solving a system of ordinary differential equations which characterizes several kinds of linear reactions and diffusion from biochemistry, physiology, etc. The methods consist of sequential applications of the simple exact solver for a reversible reaction. We provide absolute stability and convergence of the proposed explicit methods. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A dual iterative substructuring method with a penalty term

An iterative substructuring method with Lagrange multipliers is considered for second order elliptic problems, which is a variant of the FETI-DP method. The standard FETI-DP formulation is associated with the saddle-point problem which is induced from the minimization problem with a constraint for imposing the continuity across the interface. Starting from the slightly changed saddle-point problem by addition of a penalty term with a positive penalization parameter η, we propose a dual substructuring method which is implemented iteratively by the conjugate gradient method. In spite of the absence of any preconditioners, it is shown that the proposed method is numerically scalable in the sense that for a large value of η, the condition number of the resultant dual problem is bounded by a constant independent of both the subdomain size H and the mesh size h. Computational issues and numerical results are presented. This work was partially supported by the SRC/ERC program of MOST/KOSEF(R11-2002-103).     

A Locking-Free Nonconforming Finite Element Method for Planar Linear Elasticity

We propose a locking-free nonconforming finite element method based on quadrilaterals to solve for the displacement variable in the pure displacement boundary value problem of planar linear elasticity. The method proposed in this paper is optimal and robust in the sense that the convergence estimates in the energy and L ²-norms are independent of the Lamé parameter .