ON SMOOTH LU DECOMPOSITIONS WITH APPLICATIONS TO SOLUTIONS OF NONLINEAR EIGENVALUE PROBLEMS

We study the smooth LU decomposition of a given analytic functional A-matrix A（A） and its block-analogue. Sufficient conditions for the existence of such matrix decompositions are given, some differentiability about certain elements arising from them are proved, and several explicit expressions for derivatives of the specified elements are provided. By using these smooth LU decompositions, we propose two numerical methods for computing multiple nonlinear eigenvalues of A（A）, and establish their locally quadratic convergence properties. Several numerical examples are provided to show the feasibility and effectiveness of these new methods.     

Convergence theory for the exact interpolation scheme with approximation vector as the first column of the prolongator and Rayleigh quotient iteration nonlinear smoother

We extend the analysis of the recently proposed nonlinear EIS scheme applied to the partial eigenvalue problem. We address the case where the Rayleigh quotient iteration is used as the smoother on the fine-level. Unlike in our previous theoretical results, where the smoother given by the linear inverse power method is assumed, we prove nonlinear speed-up when the approximation becomes close to the exact solution. The speed-up is cubic. Unlike existent convergence estimates for the Rayleigh quotient iteration, our estimates take advantage of the powerful effect of the coarse-space.     

AN ANISOTROPIC NONCONFORMING FINITE ELEMENT METHOD FOR APPROXIMATING A CLASS OF NONLINEAR SOBOLEV EQUATIONS

An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.     

SUPERCONVERGENCE ANALYSIS OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS OF THE STATIONARY BENARD TYPE

In this paper, we consider the finite element approximation of the distributed optimal control problems of the stationary Benard type under the pointwise control constraint. The states and the co-states are approximated by polynomial functions of lowest-order mixed finite element space or piecewise linear functions and the control is approximated by piecewise constant functions. We give the superconvergence analysis for the control; it is proved that the approximation has a second-order rate of convergence. We further give the superconvergence analysis for the states and the co-states. Then we derive error estimates in L^∞-norm and optimal error estimates in L^2-norm.     

The Rayleigh quotient iteration for quadratic eigenvalue problems

For the direct solution of quadratic eigenvalue problems of the form (λ² M + P + Q ) x = 0 , a generalization of the Rayleigh quotient iteration is presented. Numerical simulations show good convergence for problems where the eigenvalues have nonzero imaginary part. The method is used to calculate eigenvalue paths of parameter dependent problems in structural dynamics. Bifurcations with double eigenvalues, which can occur in the path, are passed by using a perturbation of the velocity dependent matrix P . (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

RESTRICTED ADDITIVE SCHWARZ METHOD FOR A KIND OF NONLINEAR COMPLEMENTARITY PROBLEM

In this paper, a new Schwarz method called restricted additive Schwarz method （RAS） is presented and analyzed for a kind of nonlinear complementarity problem （NCP）. The method is proved to be convergent by using weighted maximum norm. Besides, the effect of overlap on RAS is also considered. Some preliminary numerical results are reported to compare the performance of RAS and other known methods for NCP.     

基于辛本征空间的线性阻尼振动系统动力学分析

对于考虑阻尼项和陀螺项的一般线性动力学振动系统,建立基于辛本征空间展开求解的一般方法.基于Rayleigh商本征值的模态展开方法被广泛应用于复杂结构动力系统振动分析,但对于很多机械系统,由于其不能有效考虑陀螺效应的影响,其适用性却受到很大限制.该文首先讨论了无阻尼系统Rayleigh商本征值问题与辛本征值问题的对应关系,表明前者实际可由后者的一种退化形式给出（也即忽略陀螺效应）,而后者更具有一般性.在此基础上,进一步基于辛本征空间本征向量展开,推导了同时考虑阻尼和陀螺系统的一般线性动力学系统的有效求解方法.数值算例选取不考虑陀螺效应及考虑陀螺效应的两种线性阻尼振动系统对所提出的方法进行了验证,分析结果表明了该文所建立方法的正确性和有效性.     

ERROR ESTIMATES FOR THE TIME DISCRETIZATION FOR NONLINEAR MAXWELL＇S EQUATIONS

This paper is devoted to the study of a nonlinear evolution eddy current model of the type δtB（H） ＋ △↓× （ △↓ × H） =0 subject to homogeneous Dirichlet boundary conditions H × v = 0 and a given initial datum. Here, the magnetic properties of a soft ferromagnet are linked by a nonlinear material law described by B（H）. We apply the backward Euler method for the time discretization and we derive the error estimates in suitable function spaces. The results depend on the nonlinearity of B（H）.     

LOCALLY STABILIZED FINITE ELEMENT METHOD FOR STOKES PROBLEM WITH NONLINEAR SLIP BOUNDARY CONDITIONS

Based on the low-order conforming finite element subspace （Vh, Mh） such as the P1-P0 triangle element or the Q1-P0 quadrilateral element, the locally stabilized finite element method for the Stokes problem with nonlinear slip boundary conditions is investigated in this paper. For this class of nonlinear slip boundary conditions including the subdifferential property, the weak variational formulation associated with the Stokes problem is an variational inequality. Since （Vh, Mh） does not satisfy the discrete inf-sup conditions, a macroelement condition is introduced for constructing the locally stabilized formulation such that the stability of （Vh, Mh） is established. Under these conditions, we obtain the H1 and L2 error estimates for the numerical solutions.     

Homogenization of a nonlinear degenerate parabolic differential equation

This paper is devoted to the homogenization of a nonlinear degenerate parabolic problem ɑtu∈-div(D(x/∈, u∈,▽u∈)+ K(x/∈, u∈))= f(x) with Dirichlet boundary condition. Here the operator D(y, s,s) is periodic in y and degenerated in ▽s. In the paper, under the two-scale convergence theory, we obtain the limit equation as ∈→ 0 and also prove the corrector results of ▽u∈ to strong convergence.     

A geometric theory for preconditioned inverse iteration I: Extrema of the Rayleigh quotient

The discretization of eigenvalue problems for partial differential operators is a major source of matrix eigenvalue problems having very large dimensions, but only some of the smallest eigenvalues together with the eigenvectors are to be determined. Preconditioned inverse iteration (a “matrix-free” method) derives from the well-known inverse iteration procedure in such a way that the associated system of linear equations is solved approximately by using a (multigrid) preconditioner. A new convergence analysis for preconditioned inverse iteration is presented. The preconditioner is assumed to satisfy some bound for the spectral radius of the error propagation matrix resulting in a simple geometric setup. In this first part the case of poorest convergence depending on the choice of the preconditioner is analyzed. In the second part the dependence on all initial vectors having a fixed Rayleigh quotient is considered. The given theory provides sharp convergence estimates for the eigenvalue approximations showing that multigrid eigenvalue/vector computations can be done with comparable efficiency as known from multigrid methods for boundary value problems.     

On the Determination of a Density Function by Its Autoconvolution Coefficient

We investigate eigenvalues and eigenvectors of certain linear variational eigenvalue inequalities where the constraints are defined by a convex cone as in [4], [7], [8], [10]-[12], [17]. The eigenvalues of those eigenvalue problems are of interest in connection with bifurcation from the trivial solution of nonlinear variational inequalities. A rather far reaching theory is presented for the case that the cone is given by a finite number of linear inequalities, where an eigensolution corresponds to a (+)-Kuhn-Tucker point of the Rayleigh quotient. Application to an unlaterally supported beam are discussed and numerical results are given.     

Local convergence analysis of several inexact Newton-type algorithms for general nonlinear eigenvalue problems

We study the local convergence of several inexact numerical algorithms closely related to Newton’s method for the solution of a simple eigenpair of the general nonlinear eigenvalue problem $T(\lambda )v=0$ . We investigate inverse iteration, Rayleigh quotient iteration, residual inverse iteration, and the single-vector Jacobi–Davidson method, analyzing the impact of the tolerances chosen for the approximate solution of the linear systems arising in these algorithms on the order of the local convergence rates. We show that the inexact algorithms can achieve the same order of convergence as the exact methods if appropriate sequences of tolerances are applied to the inner solves. We discuss the connections and emphasize the differences between the standard inexact Newton’s method and these inexact algorithms. When the local symmetry of $T(\lambda )$ is present, the use of a nonlinear Rayleigh functional is shown to be fundamental in achieving higher order of convergence rates. The convergence results are illustrated by numerical experiments.     

Interaction of faults under slip‐dependent friction. Non‐linear eigenvalue analysis

We analyse the evolution of a system of finite faults by considering the non‐linear eigenvalue problems associated to static and dynamic solutions on unbounded domains. We restrict our investigation to the first eigenvalue (Rayleigh quotient). We point out its physical significance through a stability analysis and we give an efficient numerical algorithm able to compute it together with the corresponding eigenfunction. We consider the anti‐plane shearing on a system of finite faults under a slip‐dependent friction in a linear elastic domain, not necessarily bounded. The static problem is formulated in terms of local minima of the energy functional. We introduce the non‐linear (static) eigenvalue problem and we prove the existence of a first eigenvalue/eigenfunction characterizing the isolated local minima. For the dynamic problem, we discuss the existence of solutions with an exponential growth, to deduce a (dynamic) non‐linear eigenvalue problem. We prove the existence of a first dynamic eigenvalue and we analyse its behaviour with respect to the friction parameter. We deduce a mixed finite element discretization of the non‐linear spectral problem and we give a numerical algorithm to approach the first eigenvalue/eigenfunction. Finally we give some numerical results which include convergence tests, on a single fault and a two‐faults system, and a comparison between the non‐linear spectral results and the time evolution results. Copyright © 2004 John Wiley & Sons, Ltd.     

Equations of second order in time with quasilinear damping: existence in Orlicz spaces via convergence of a full discretisation

A nonlinear evolution equation of second order with damping is studied. The quasilinear damping term is monotone and coercive but exhibits anisotropic and nonpolynomial growth. The appropriate setting for such equations is that of monotone operators in Orlicz spaces. Global existence of solutions in the sense of distributions is shown via convergence of the backward Euler scheme combined with an internal approximation. Copyright © 2015 John Wiley & Sons, Ltd.     

复合材料特征值的高阶多尺度Rayleigh商校正

本文讨论了周期结构复合材料特征值的多尺度计算,提出了高阶多尺度Rayleigh商校正算法,并给出了收敛性分析. 最后,通过大量数值实验结果表明,新算法是有效且必要的.     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP --In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

A MONOTONE COMPACT IMPLICIT SCHEME FOR NONLINEAR REACTION-DIFFUSION EQUATIONS

A monotone compact implicit finite difference scheme with fourth-order accuracy in space and second-order in time is proposed for solving nonlinear reaction-diffusion equations. An accelerated monotone iterative method for the resulting discrete problem is presented. The sequence of iteration converges monotonically to the unique solution of the discrete problem, and the convergence rate is either quadratic or nearly quadratic, depending on the property of the nonlinear reaction. The numerical results illustrate the high accuracy of the proposed scheme and the rapid convergence rate of.the iteration.     

A geometric theory for preconditioned inverse iteration IV: On the fastest convergence cases

In order to compute the smallest eigenvalue together with an eigenfunction of a self-adjoint elliptic partial differential operator one can use the preconditioned inverse iteration scheme, also called the preconditioned gradient iteration. For this iterative eigensolver estimates on the poorest convergence have been published by several authors. In this paper estimates on the fastest possible convergence are derived. To this end the convergence problem is reformulated as a two-level constrained optimization problem for the Rayleigh quotient. The new convergence estimates reveal a wide range between the fastest possible and the slowest convergence.     

Computation of eigenvalues and solutions of regular Sturm–Liouville problems using Haar wavelets

The paper presents a novel method for the computation of eigenvalues and solutions of Sturm–Liouville eigenvalue problems (SLEPs) using truncated Haar wavelet series. This is an extension of the technique proposed by Hsiao to solve discretized version of variational problems via Haar wavelets. The proposed method aims to cover a wider class of problems, by applying it to historically important and a very useful class of boundary value problems, thereby enhancing its applicability. To demonstrate the effectiveness and efficiency of the method various celebrated Sturm–Liouville problems are analyzed for their eigenvalues and solutions. Also, eigensystems are investigated for their asymptotic and oscillatory behavior. The proposed scheme, unlike the conventional numerical schemes, such as Rayleigh quotient and Rayleigh–Ritz approximation, gives eigenpairs simultaneously and provides upper and lower estimates of the smallest eigenvalue, and it is found to have quadratic convergence with increase in resolution.