A p-version finite element method for nonlinear elliptic variational inequalities in 2D

This article introduces and analyzes a p-version FEM for variational inequalities resulting from obstacle problems for some quasi-linear elliptic partial differential operators. We approximate the solution by controlling the obstacle condition in images of the Gauss–Lobatto points. We show existence and uniqueness for the discrete solution u _p from the p-version for the obstacle problem. We prove the convergence of u _p towards the solution with respect to the energy norm, and assuming some additional regularity for the solution we derive an a priori error estimate. In numerical experiments the p-version turns out to be superior to the h-version concerning the convergence rate and the number of unknowns needed to achieve a certain exactness of the approximation.     

Asymptotic and boundary integral formulae for the computational reconstruction of small conductivity imperfections

In this paper we investigate the electrostatic problem of determining conductivity profiles from the knowledge of boundary currents and voltages. We obtain an improved estimate for the voltage potential of a two-dimensional conductor having finitely many circular inclusions and piecewise constant conductivity profile. We derive an asymptotic expansion for the voltage potential in terms of the reference voltage potential and the location, size, and conductivity of the inhomogeneities. This representation is used to formulate the nonlinear least squares problem for estimating the location and size of the inhomogeneities. Required boundary data for the voltage potential are generated numerically by solving a system of integral equations. Computational experiments are presented to demonstrate the effectiveness of our identification procedure.     

Riccati-based preconditioner for computing invariant subspaces of large matrices

Summary. This paper introduces and analyzes the convergence properties of a method that computes an approximation to the invariant subspace associated with a group of eigenvalues of a large not necessarily diagonalizable matrix. The method belongs to the family of projection type methods. At each step, it refines the approximate invariant subspace using a linearized Riccati's equation which turns out to be the block analogue of the correction used in the Jacobi-Davidson method. The analysis conducted in this paper shows that the method converges at a rate quasi-quadratic provided that the approximate invariant subspace is close to the exact one. The implementation of the method based on multigrid techniques is also discussed and numerical experiments are reported. Received June 15, 2000 / Revised version received January 22, 2001 / Published online October 17, 2001     

On a generalized $L_2$ projection and some related stability estimates in Sobolev spaces

Summary. In this paper we investigate a stability estimate needed in hybrid finite and boundary element methods, especially in hybrid coupled domain decomposition methods including mortar finite elements. This stability estimate is equivalent to the stability of a generalized projection in certain Sobolev spaces. Using piecewise linear trial spaces and appropriate piecewise constant test spaces, the stability of the generalized projection is proved assuming some mesh conditions locally. Received April 11, 2000 / Revised version received February 15, 2001 / Published online July 25, 2001     

A supraconvergent scheme for the Korteweg-de Vries equation

Summary A finite-difference method for the integration of the Korteweg-de Vries equation on irregular grids is analyzed. Under periodic boundary conditions, the method is shown to be supraconvergent in the sense that, though being inconsistent, it is second order convergent. However, such a convergence only takes place on grids with an odd number of points per period. When a grid with an even number of points is used, the inconsistency of the method leads to divergence. Numerical results backing the analysis are presented.     

On the convergence of Newton's method when estimating higher dimensional parameters

In this paper, we consider the estimation of a parameter of interest where the estimator is one of the possibly several solutions of a set of nonlinear empirical equations. Since Newton's method is often used in such a setting to obtain a solution, it is important to know whether the so obtained iteration converges to the locally unique consistent root to the aforementioned parameter of interest. Under some conditions, we show that this is eventually the case when starting the iteration from within a ball about the true parameter whose size does not depend on n. Any preliminary almost surely consistent estimate will eventually lie in such a ball and therefore provides a suitable starting point for large enough n. As examples, we will apply our results in the context of M-estimates, kernel density estimates, as well as minimum distance estimates.     

On the fast reduction of a quasiseparable matrix to Hessenberg and tridiagonal forms

In this paper we design a fast new algorithm for reducing an N × N quasiseparable matrix to upper Hessenberg form via a sequence of N − 2 unitary transformations. The new reduction is especially useful when it is followed by the QR algorithm to obtain a complete set of eigenvalues of the original matrix. In particular, it is shown that in a number of cases some recently devised fast adaptations of the QR method for quasiseparable matrices can benefit from using the proposed reduction as a preprocessing step, yielding lower cost and a simplification of implementation.     

The tetrahedron algebra and its finite-dimensional irreducible modules

Recently Terwilliger and the present author found a presentation for the three-point sl₂ loop algebra via generators and relations. To obtain this presentation we defined a Lie algebra ? by generators and relations and displayed an isomorphism from ? to the three-point sl₂ loop algebra. In this paper we classify the finite-dimensional irreducible ?-modules.     

Fabrication of Long-Period Fibre Gratings Using 800 nm Femtosecond Laser Pulses

Long-period fibre gratings inside standard single-mode optical communication fibres are successfully fabricated with infrared femtosecond laser pulses. The refractive index perturbations are well confined within the fibre core by choosing the proper laser focusing parameters and translation speed of the fibre during the direct laser writing process. With the self-focusing effect considered and at a constant average irradiation dose of 1.62 × 10^3 J/（cm^2μm）, the threshold intensity for fabricating long-period gratings with infrared femtosecond laser pulses is determined to be 5.13 × 10^13 W/cm2.     

Finite element error estimates for non-linear elliptic equations of monotone type

Summary In this paper we shall consider the application of the finite element method to a class of second order elliptic boundary value problems of divergence form and with gradient nonlinearity in the principal coefficient, and the derivation of error estimates for the finite element approximations. Such problems arise in many practical situations — for example, in shock-free airfoil design, seepage through coarse grained porous media, and in some glaciological problems. By making use of certain properties of the nonlinear coefficients, we shall demonstrate that the variational formulations associated with these boundary value problems are well-posed. We shall also prove that the abstract operators accompanying such problems satisfy certain continuity and monotonicity inequalities. With the aid of these inequalities and some standard results from approximation theory, we show how one may derive error estimates for the finite element approximations in the energy norm.