Application of the pseudo-spectral method to 2D eigenvalue problems in elasticity

A pseudo-spectral approach to 2D vibrational problems arising in linear elasticity is considerede using differentiation matrices. The governing partial differential equations and associated boundary conditions on regular domains can be translated into matrix eigenvalue problems. Accurate results are obtained to the precision expected in spectral-type methods. However, we show that it is necessary to apply an additional “pole” condition to deal with ther=0 coordinate singularity arising in the case of a 2D disc.     

Mixed impedance transmission problems for vibrating layered elastic bodies

Direct scattering problems for partially coated piecewise homogenous and inhomogeneous layered obstacles in linear elasticity lead to mixed impedance transmission problems for the steady‐state elastic oscillation equations. For a piecewise homogenous isotropic composite body, we employ the potential method and reduce the mixed impedance transmission problem to an equivalent system of boundary pseudodifferential equations. We give a detailed analysis of the corresponding pseudodifferential operators, which live on the interface between the layers and on a proper submanifold of the boundary of the composite elastic body, and establish uniqueness and existence results for the original mixed impedance transmission problem for arbitrary values of the oscillation frequency parameter; this is crucial in the study of inverse elastic scattering problems for partially coated layered obstacles. We also investigate regularity properties of solutions near the collision curves, where the different boundary conditions collide, and establish almost best Hölder smoothness results. Further, we analyze the asymptotic behavior of the stress vector near the collision curve and derive explicit formulas for the stress singularity exponents. The case of Lipschitz surfaces is briefly treated separately. In the case of a composite body containing homogeneous or inhomogeneous finite anisotropic inclusions, we develop an alternative hybrid method based on the so‐called nonlocal approach and reduce the mixed transmission problem to an equivalent functional‐variational equation with a sesquilinear form that ‘lives’ on a bounded part of the layered composite body and its boundary. We show that this sesquilinear form is coercive and that the corresponding variational equation is uniquely solvable. Copyright © 2015 John Wiley & Sons, Ltd.     

Some remarks on the constant in the strengthened CBS inequality: Estimate for hierarchical finite element discretizations of elasticity problems

For a class of two‐dimensional boundary value problems including diffusion and elasticity problems, it is proved that the constants in the corresponding strengthened Cauchy‐Buniakowski‐Schwarz (CBS) inequality in the cases of two‐level hierarchical piecewise‐linear/piecewise‐linear and piecewise‐linear/piecewise‐quadratic finite element discretizations with triangular meshes differ by the factor 0.75. For plane linear elasticity problems and triangulations with right isosceles triangles, formulas are presented that show the dependence of the constant in the CBS inequality on the Poisson's ratio. Furthermore, numerically determined bounds of the constant in the CBS inequality are given for plane linear elasticity problems discretized by means of arbitrary triangles and for three‐dimensional elasticity problems discretized by means of tetrahedral elements. Finally, the robustness of iterative solvers for elasticity problems is discussed briefly. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 469–487, 1999     

Application of the Pseudo-Spectral Method to 2D Eigenvalue Problems in Elasticity

A pseudo-spectral approach to 2D vibrational problems arising in linear elasticity is considered using differentiation matrices. The governing partial differential equations and associated boundary conditions on regular domains can be translated into matrix eigenvalue problems. Accurate results are obtained to the precision expected in spectral-type methods. However, we show that it is necessary to apply an additional pole condition to deal with the r=0 coordinate singularity arising in the case of a 2D disc.     

TV‐like regularization for backward parabolic problems

This paper investigates an inverse problem for parabolic equations backward in time, which is solved by total‐variation‐like (TV‐like, in abbreviation) regularization method with cost function ∥u_x∥². The existence, uniqueness and stability estimate for the regularization problem are deduced in the linear case. For numerical illustration, the variational adjoint method, which presents a simple method to derive the gradient of the optimization functional, is introduced to reconstruct the unknown initial condition for both linear and nonlinear parabolic equations. The conjugate gradient method is used to iteratively search for the optimal approximation. Numerical results validate the feasibility and effectiveness of the proposed algorithm. Copyright © 2016 John Wiley & Sons, Ltd.     

Sturm–Liouville properties for Atkinson's semi‐definite p‐Laplacian eigenvalue problems

We define Atkinson's semi‐definite p‐Laplacian eigenvalue problems, which include the regular p‐Laplacian eigenvalue problems with L ¹ coefficient functions. Then we show that the Sturm oscillation theorem also holds for this eigenvalue problem.     

Discontinuous Galerkin methods for periodic boundary value problems

This article considers the extension of well‐known discontinuous Galerkin (DG) finite element formulations to elliptic problems with periodic boundary conditions. Such problems routinely appear in a number of applications, particularly in homogenization of composite materials. We propose an approach in which the periodicity constraint is incorporated weakly in the variational formulation of the problem. Both H¹ and L² error estimates are presented. A numerical example confirming theoretical estimates is shown. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Analysis of Two Linearized Problems Modeling Viscous Two‐Layer Flows

Two problems that appear in the linearization of certain free boundary value problems of the hydrodynamics of two viscous fluids are studied in the strip‐like domain Π = {x = (x₁, x₂) ∈ ℝ² : x₁ ∈ ℝ¹, (0 < x₂ < h_*) ∨ (h_* < x₂ < 1)}. The first problem arises in the linearization of a two‐layer flow down a geometrically perturbed inclined plane. The second one appears after the linearization of a two‐layer flow in a geometrically perturbed inclined channel with one moving (smooth) wall. For this purpose the unknown flow domain was mapped onto the double strip Π. The arising linear elliptic problems contain additional unknown functions in the boundary conditions. The paper is devoted to the investigation of these boundary problems by studying the asymptotics of the eigenvalues of corresponding operator pencils. It can be proved that the boundary value problems are uniquely solvable in weighted Sobolev spaces with exponential weight. The study of the full (nonlinear) free boundary value problems will be the topic of a forthcoming paper.     

Efficient algorithms for three‐dimensional axial and planar random assignment problems

Beautiful formulas are known for the expected cost of random two‐dimensional assignment problems, but in higher dimensions even the scaling is not known. In three dimensions and above, the problem has natural “Axial” and “Planar” versions, both of which are NP‐hard. For 3‐dimensional Axial random assignment instances of size n, the cost scales as Ω(1/ n), and a main result of the present paper is a linear‐time algorithm that, with high probability, finds a solution of cost O(n^–1+o(1)). For 3‐dimensional Planar assignment, the lower bound is Ω(n), and we give a new efficient matching‐based algorithm that with high probability returns a solution with cost O(n log n). © 2013 Wiley Periodicals, Inc. Random Struct. Alg., 46, 160–196, 2015     

A posteriori error estimates for discontinuous Galerkin method to the elasticity problem

This work concerns with the discontinuous Galerkin (DG) method for the time‐dependent linear elasticity problem. We derive the a posteriori error bounds for semidiscrete and fully discrete problems, by making use of the stationary elasticity reconstruction technique which allows to estimate the error for time‐dependent problem through the error estimation of the associated stationary elasticity problem. For fully discrete scheme, we make use of the backward‐Euler scheme and an appropriate space‐time reconstruction. The technique here can be applicable for a variety of DG methods as well.     

Correction of finite element eigenvalues for problems with natural or periodic boundary conditions

Recent results of Andrew and Paine for a regular Sturm-Liouville problem with essential boundary conditions are extended to problems with natural or periodic boundary conditions. These results show that a simple asymptotic correction technique of Paine, de Hoog and Anderssen reduces the error in the estimate of thekth eigenvalue obtained by the finite element method, with linear hat functions and mesh lengthh, fromO(k ⁴ h ²) toO(kh ²). Numerical results show the correction to be useful even for low values ofk.     

Long time approximations for solutions of wave equations associated with the Steklov spectral homogenization problems

The interest in the use of quasimodes, or almost frequencies and almost eigenfunctions, to describe asymptotics for low‐frequency and high‐frequency vibrations in certain singularly perturbed spectral problems, which depend on a small parameter ε, has been recently highlighted in many papers. In this paper we deal with the low frequencies for a Steklov‐type eigenvalue homogenization problem: we consider harmonic functions in a bounded domain of ?², and strongly alternating boundary conditions of the Dirichlet and Steklov type on a part of the boundary. The spectral parameter appears in the boundary condition on small segments T^ε of size O(ε) periodically distributed along the boundary; ε also measures the periodicity of the structure. We consider associated second‐order evolution problems on spaces of traces that depend on ε, and we provide estimates for the time t in which standing waves, constructed from quasimodes, approach their solutions u^ε(t) as ε→0. Copyright © 2009 John Wiley & Sons, Ltd.     

Memory‐efficient Arnoldi algorithms for linearizations of matrix polynomials in Chebyshev basis

Novel memory‐efficient Arnoldi algorithms for solving matrix polynomial eigenvalue problems are presented. More specifically, we consider the case of matrix polynomials expressed in the Chebyshev basis, which is often numerically more appropriate than the standard monomial basis for a larger degree d. The standard way of solving polynomial eigenvalue problems proceeds by linearization, which increases the problem size by a factor d. Consequently, the memory requirements of Krylov subspace methods applied to the linearization grow by this factor. In this paper, we develop two variants of the Arnoldi method that build the Krylov subspace basis implicitly, in a way that only vectors of length equal to the size of the original problem need to be stored. The proposed variants are generalizations of the so‐called quadratic Arnoldi method and two‐level orthogonal Arnoldi procedure methods, which have been developed for the monomial case. We also show how the typical ingredients of a full implementation of the Arnoldi method, including shift‐and‐invert and restarting, can be incorporated. Numerical experiments are presented for matrix polynomials up to degree 30 arising from the interpolation of nonlinear eigenvalue problems, which stem from boundary element discretizations of PDE eigenvalue problems. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical solution of singular boundary value problems via Chebyshev polynomial and B-spline

In this paper, we present B-spline method for numerically solving singular two-point boundary value problems for certain ordinary differential equation having singular coefficients.These problems arise when reducing partial differential equation to ordinary differential equation by physical symmetry. To remove the singularity, we first use Chebyshev economizition in the vicinity of the singular point and the boundary condition at a point x=δ (in the vicinity of the singularity) is derived. The resulting regular BVP is then efficiently treated by employing B-spline for finding the numerical solution. Some examples have been included and comparison of the numerical results made with other methods.     

Boundary value problems for higher order linear impulsive differential equations

In this paper higher order linear impulsive differential equations with fixed moments of impulses subject to linear boundary conditions are studied. Green's formula is defined for piecewise differentiable functions. Properties of Green's functions for higher order impulsive boundary value problems are introduced. An appropriate example of the Green's function for a boundary value problem is provided. Furthermore, eigenvalue problems and basic properties of eigensolutions are considered.     

Eigenvalue problems to compute almost optimal points for rational interpolation with prescribed poles

Explicit formulas exist for the (n,m) rational function with monic numerator and prescribed poles that has the smallest possible Chebyshev norm. In this paper we derive two different eigenvalue problems to obtain the zeros of this extremal function. The first one is an ordinary tridiagonal eigenvalue problem based on a representation in terms of Chebyshev polynomials. The second is a generalised tridiagonal eigenvalue problem which we derive using a connection with orthogonal rational functions. In the polynomial case (m = 0) both problems reduce to the tridiagonal eigenvalue problem associated with the Chebyshev polynomials of the first kind. Postdoctoral researcher FWO-Flanders.     

Boundary value problems with eigenvalue depending boundary conditions

We investigate some classes of eigenvalue dependent boundary value problems of the form where A ? A⁺ is a symmetric operator or relation in a Krein space K, τ is a matrix function and Γ₀, Γ₁ are abstract boundary mappings. It is assumed that A admits a self‐adjoint extension in K which locally has the same spectral properties as a definitizable relation, and that τ is a matrix function which locally can be represented with the resolvent of a self‐adjoint definitizable relation. The strict part of τ is realized as the Weyl function of a symmetric operator T in a Krein space H, a self‐adjoint extension à of A × T in K × H with the property that the compressed resolvent P_K (à – λ)^–1|_K k yields the unique solution of the boundary value problem is constructed, and the local spectral properties of this so‐called linearization à are studied. The general results are applied to indefinite Sturm–Liouville operators with eigenvalue dependent boundary conditions (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Local mesh refinement for the discretization of Neumann boundary control problems on polyhedra

This paper deals with optimal control problems constrained by linear elliptic partial differential equations. The case where the right‐hand side of the Neumann boundary is controlled, is studied. The variational discretization concept for these problems is applied, and discretization error estimates are derived. On polyhedral domains, one has to deal with edge and corner singularities, which reduce the convergence rate of the discrete solutions, that is, one cannot expect convergence order two for linear finite elements on quasi‐uniform meshes in general. As a remedy, a local mesh refinement strategy is presented, and a priori bounds for the refinement parameters are derived such that convergence with optimal rate is guaranteed. As a by‐product, finite element error estimates in the H¹(Ω)‐norm, L²(Ω)‐norm and L²(Γ)‐norm for the boundary value problem are obtained, where the latter one turned out to be the main challenge. Copyright © 2015 John Wiley & Sons, Ltd.     

Higher order convergence results for the Rayleigh-Ritz method applied to eigenvalue problems: 2. Improved error bounds for eigenfunctions

The application of the Rayleigh-Ritz method for approximating the eigenvalues and eigenfunctions of linear eigenvalue problems in several dimensions is investigated. The object is to improve upon known error estimates for the approximate eigenfunctions. Results for the Galerkin approximation of the eigenfunctions are developed under varying assumptions on the boundary conditions and domain of definition of the eigenvalue problem. These results, coupled with a previous result relating Galerkin and Rayleigh-Ritz approximation of the eigenfunctions, are then used to obtain improved error estimates for the approximate eigenfunctions in theL ² and uniform norms.This research was supported in part by AEC Grant (11-1)-2075.