Mass‐lumping or not mass‐lumping for eigenvalue problems

In this article we analyze the effect of mass‐lumping in the linear triangular finite element approximation of second‐order elliptic eigenvalue problems. We prove that the eigenvalue obtained by using mass‐lumping is always below the one obtained with exact integration. For singular eigenfunctions, as those arising in non convex polygons, we prove that the eigenvalue obtained with mass‐lumping is above the exact eigenvalue when the mesh size is small enough. So, we conclude that the use of mass‐lumping is convenient in the singular case. When the eigenfunction is smooth several numerical experiments suggest that the eigenvalue computed with mass‐lumping is below the exact one if the mesh is not too coarse. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 653–664, 2003     

Conformable fractional Sturm‐Liouville equation

In this article, we discuss a conformable fractional Sturm‐Liouville boundary‐value problem. We prove an existence and uniqueness theorem for this equation and formulate a self‐adjoint boundary value problem. We also construct the associated Green function of this problem, and we give the eigenfunction expansions. Finally, we will give some examples.     

Sufficient conditions of non‐uniqueness for the Coulomb friction problem

We consider the Signorini problem with Coulomb friction in elasticity. Sufficient conditions of non‐uniqueness are obtained for the continuous model. These conditions are linked to the existence of real eigenvalues of an operator in a Hilbert space. We prove that, under appropriate conditions, real eigenvalues exist for a non‐local Coulomb friction model. Finite element approximation of the eigenvalue problem is considered and numerical experiments are performed. Copyright © 2003 John Wiley & Sons, Ltd.     

Three‐point bending tests—Part I: Mathematical study and asymptotic analysis

The goal of this work is to study the static behaviour of a three‐dimensional elastic beam when is subjected to a three‐point bending test. In the first part, under suitable compatibility conditions on the applied forces and on the geometry of the beam, we will prove the existence of a unique solution for the associated contact elastic problem; these conditions of compatibility on the data come from the absence of a Dirichlet condition on the beam boundary. In the second part, we will study the asymptotic behaviour of this problem; in particular, we will deduce the one‐dimensional models associated with the displacement components, and we will give the existence and uniqueness of solution for them. Moreover, we will give an expression for the normal axial stress in the beam which is related to the modulus of rupture of brittle materials. In the final part of the work, we will deal with the regularity of the solution for the bending problem and we will prove some properties of the coincidence set. Copyright © 2011 John Wiley & Sons, Ltd.     

Contrôle Optimal de Systèmes Gouvernés par des Problèmes aux Valeurs Propres

In this paper, we study the optimal control of systems governed by an eigenvalue problem. We prove the existence of an optimal control and we obtain first order optimality conditions. We present also two numerical examples and we give the complete solution of both problems, involving approximation by the finite element method and minimization by a gradient method. Finally numerical algorithms are detailed.     

Galerkin method for optimal control of second‐order evolution equations

This paper presents the Galerkin approximation of the optimization problem of a system governed by non‐linear second‐order evolution equation where a non‐linear operator depends on derivative of the state of the system. The control is acting on a non‐linear equation. After giving some results on the existence of optimal control we shall prove the existence of the weak and strong condensation points of a set of solutions of the approximate optimization problems. Each of these points is a solution of the initial optimization problem. Finally we shall give a simple example using the obtained results. Copyright © 2004 John Wiley & Sons, Ltd.     

Optimal Shape Problems for Eigenvalues

ABSTRACT

We consider the first Dirichlet eigenvalue for nonhomogeneous membranes. For given volume we want to find the domain which minimizes this eigenvalue. The problem is formulated as a variational free boundary problem. The optimal domain is characterized as the support of the first eigenfunction. We prove enough regularity for the eigenfunction to conclude that the optimal domain has finite parameter. Finally an overdetermined boundary value problem on the regular part of the free boundary is given.     

Global existence of solutions for non‐small data to non‐linear spherically symmetric thermoviscoelasticity

We consider some initial–boundary value problems for non‐linear equations of thermoviscoelasticity in the three‐dimensional case. Since, we are interested to prove global existence we consider spherically symmetric problem. We examine the Neumann conditions for the temperature and either the Neumann or the Dirichlet boundary conditions for the elasticity equations. Using the energy method, we are able to obtain some energy estimates in appropriate Sobolev spaces enough to prove existence for all time without any restrictions on data. Due to the spherical symmetricity the constants in the above estimates increase with time so the existence for all finite times is proved only. Copyright © 2003 John Wiley & Sons, Ltd.     

Asymptotic expansion and extrapolation for the Eigenvalue approximation of the biharmonic eigenvalue problem by Ciarlet–Raviart scheme

We propose and analyze the Ciarlet–Raviart mixed scheme for solving the biharmonic eigenvalue problem with bilinear finite elements. We derive a higher order convergence rate for eigenvalue and eigenfunction approximations. Furthermore, we give an asymptotic expansion of the eigenvalue error, from which an efficient extrapolation and an a posteriori error estimate for the eigenvalue are given. Finally, numerical experiments illustrating the theoretical results are reported. This author was supported by China Postdoctoral Sciences Foundation.     

复Monge－Ampere方程的特征值问题

讨论了复空间中强拟凸域上的复Ｍｏｎｇｅ－Ａｍｐｅｒｅ方程的特征值问题，证明了特征值问题解的存在唯一性，并给出了这个特征值与一类复空间中复Ｌａｐｌａｃｅ算子的第一特征值的关系，最后利用特征值及特征函数的存在性讨论了一类复Ｍｏｎｇｅ－Ａｍｐｅｒｅ方程的解的存在性及其分歧．     

Existence and approximation results for shape optimization problems in rotordynamics

We consider a shape optimization problem in rotordynamics where the mass of a rotor is minimized subject to constraints on the natural frequencies. Our analysis is based on a class of rotors described by a Rayleigh beam model including effects of rotary inertia and gyroscopic moments. The solution of the equation of motion leads to a generalized eigenvalue problem. The governing operators are non-symmetric due to the gyroscopic terms. We prove the existence of solutions for the optimization problem by using the theory of compact operators. For the numerical treatment of the problem a finite element discretization based on a variational formulation is considered. Applying results on spectral approximation of linear operators we prove that the solution of the discretized optimization problem converges towards the solution of the continuous problem if the discretization parameter tends to zero. Finally, a priori estimates for the convergence order of the eigenvalues are presented and illustrated by a numerical example.     

Blowup of solutions for a class of non‐linear evolution equations with non‐linear damping and source terms

We consider the blowup of solutions of the initial boundary value problem for a class of non‐linear evolution equations with non‐linear damping and source terms. By using the energy compensation method, we prove that when p>max{m, α}, where m, α and p are non‐negative real numbers and m+1, α+1, p+1 are, respectively, the growth orders of the non‐linear strain terms, damping term and source term, under the appropriate conditions, any weak solution of the above‐mentioned problem blows up in finite time. Comparison of the results with the previous ones shows that there exist some clear condition boundaries similar to thresholds among the growth orders of the non‐linear terms, the states of the initial energy and the existence and non‐existence of global weak solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

Multiparameter spectral problem for some weakly coupled systems of Hamiltonian equations

We consider a multiparameter spectral problem for a weakly coupled system of ordinary differential equations in which every equation is Hamiltonian and contains two unknown functions. Using the notion of the number of an eigenvalue for a problem with one such equation, we give a statement of the problem of finding the desired eigentuple of values for the problem with several equations. We prove the existence and uniqueness of a solution of this problem and suggest and study a numerical solution method.     

A non‐linear and non‐local boundary condition for a diffusion equation in petroleum engineering

This article deals with a boundary value problem for Laplace equation with a non‐linear and non‐local boundary condition. This problem comes from petroleum engineering and is used to obtain an estimation of well productivity. The non‐linear and non‐local boundary condition is written on the well boundary. On the outer reservoir boundaries, we have both Dirichlet and Neumann conditions. In this paper, we prove the existence and uniqueness of a solution to this problem. The existence is proved by Schauder theorem and the uniqueness is obtained under more restricted conditions, when the involved operator is a contraction. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical analysis of a non‐singular boundary integral method: Part I. The circular case

In order to numerically solve the interior and the exterior Dirichlet problems for the Laplacian operator, we present here a method which consists in inverting, on a finite element space, a non‐singular integral operator. This operator is a geometrical perturbation of the Steklov operator, and we precisely define the relation between the geometrical perturbation and the dimension of the finite element space, in order to obtain a stable and convergent scheme. Furthermore, this numerical scheme does not give rise to any singular integral. The scheme can also be considered as a special quadrature formula method for the standard piecewise linear Galerkin approximation of the weakly singular single layer potential, the special quadrature formula being defined by the introduction of a neighbouring curve. In the present paper, we prove stability and we give error estimates of our numerical scheme when the Laplace problem is set on a disk. We will extend our results to any domains by using compact perturbation arguments, in a second paper. Copyright © 2001 John Wiley & Sons, Ltd.     

An efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations

In this paper, we propose an efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations. Firstly, we carry out a truncation from a three‐dimensional unbounded domain to a bounded spherical domain. By using spherical coordinate transformation and spherical harmonic expansion, we transform the original problem into a series of one‐dimensional eigenvalue problem that can be solved effectively. Secondly, we introduce a weighted Sobolev space to treat the singularity in the effective potential. Using the property of orthogonal polynomials in weighted Sobolev space, the error estimate for the approximate eigenvalues and corresponding eigenfunctions are proved. Error estimates show that our numerical method can achieve spectral accuracy for approximate eigenvalues and eigenfunctions. Finally, we give some numerical examples to demonstrate the efficiency of our algorithms and the correctness of the theoretical results.     

Multiscale finite element discretizations based on local defect correction for the biharmonic eigenvalue problem of plate buckling

In this paper, we study the multiscale finite element discretizations about the biharmonic eigenvalue problem of plate buckling. On the basis of the work of Dai and Zhou (SIAM J. Numer. Anal. 46[1] [2008] 295‐324), we establish a three‐scale scheme, a multiscale discretization scheme, and the associated parallel version based on local defect correction. We first prove a local priori error estimate of finite element approximations, then give the error estimates of multiscale discretization schemes. Theoretical analysis and numerical experiments indicate that our schemes are suitable and efficient for eigenfunctions with local low smoothness.     

Postprocessing and higher order convergence for the mixed finite element approximations of the eigenvalue problem

In this paper, we propose a method to improve the convergence rate of the lowest order Raviart-Thomas mixed finite element approximations for the second order elliptic eigenvalue problem. Here, we prove a supercloseness result for the eigenfunction approximations and use a type of finite element postprocessing operator to construct an auxiliary source problem. Then solving the auxiliary additional source problem on an augmented mixed finite element space constructed by refining the mesh or by using the same mesh but increasing the order of corresponding mixed finite element space, we can increase the convergence order of the eigenpair approximation. This postprocessing method costs less computation than solving the eigenvalue problem on the finer mesh directly. Some numerical results are used to confirm the theoretical analysis.     

Quasi‐optimal preconditioners for finite element approximations of diffusion dominated convection–diffusion equations on (nearly) equilateral triangle meshes

The paper is devoted to the spectral analysis of effective preconditioners for linear systems obtained via a finite element approximation to diffusion‐dominated convection–diffusion equations. We consider a model setting in which the structured finite element partition is made by equilateral triangles. Under such assumptions, if the problem is coercive and the diffusive and convective coefficients are regular enough, then the proposed preconditioned matrix sequences exhibit a strong eigenvalue clustering at unity, the preconditioning matrix sequence and the original matrix sequence are spectrally equivalent, and under the constant coefficients assumption, the eigenvector matrices have a mild conditioning. The obtained results allow to prove the conjugate gradient optimality and the generalized minimal residual quasi‐optimality in the case of structured uniform meshes. The interest of such a study relies on the observation that automatic grid generators tend to construct equilateral triangles when the mesh is fine enough. Numerical tests, both on the model setting and in the non‐structured case, show the effectiveness of the proposal and the correctness of the theoretical findings. Copyright © 2014 John Wiley & Sons, Ltd.     

Riesz basis property of root vectors of non‐self‐adjoint operators generated by aircraft wing model in subsonic airflow

This paper is the third in a series of several works devoted to the asymptotic and spectral analysis of a model of an aircraft wing in a subsonic air flow. This model has been developed in the Flight Systems Research Center of UCLA and is presented in the works by Balakrishnan. The model is governed by a system of two coupled integro‐differential equations and a two‐parameter family of boundary conditions modeling the action of the self‐straining actuators. The differential parts of the above equations form a coupled linear hyperbolic system; the integral parts are of the convolution type. The system of equations of motion is equivalent to a single operator evolution–convolution equation in the energy space. The Laplace transform of the solution of this equation can be represented in terms of the so‐called generalized resolvent operator, which is an operator‐valued function of the spectral parameter. This generalized resolvent operator is a finite‐meromorphic function on the complex plane having the branch cut along the negative real semi‐axis. Its poles are precisely the aeroelastic modes and the residues at these poles are the projectors on the generalized eigenspaces. In the first two papers (see [33, 34]) and in the present one, our main object of interest is the dynamics generator of the differential parts of the system. This generator is a non‐self‐adjoint operator in the energy space with a purely discrete spectrum. In the first paper, we have shown that the spectrum consists of two branches, and have derived their precise spectral asymptotics with respect to the eigenvalue number. In the second paper, we have derived the asymptotical approximations for the mode shapes. Based on the asymptotical results of the first two papers, in the present paper, we (a) prove that the set of the generalized eigenvectors of the aforementioned differential operator is complete in the energy space; (b) construct the set of vectors which is biorthogonal to the set of the generalized eigenvectors in the case when there might be not only eigenvectors but associate vectors as well; and (c) prove that the set of the generalized eigenvectors forms a Riesz basis in the energy space. To prove the main result of the paper, we made use of the Nagy–Foias functional model for non‐self‐adjoint operators. The results of all three papers will be important for the reconstruction of the solution of the original initial‐boundary‐value problem from its Laplace transform in the forthcoming papers. Copyright © 2000 John Wiley & Sons, Ltd.