Bounds and extremal domains for Robin eigenvalues with negative boundary parameter

We present some new bounds for the first Robin eigenvalue with a negative boundary parameter. These include the constant volume problem, where the bounds are based on the shrinking coordinate method, and a proof that in the fixed perimeter case the disk maximises the first eigenvalue for all values of the parameter. This is in contrast with what happens in the constant area problem, where the disk is the maximiser only for small values of the boundary parameter. We also present sharp upper and lower bounds for the first eigenvalue of the ball and spherical shells. These results are complemented by the numerical optimisation of the first four and two eigenvalues in 2 and 3 dimensions, respectively, and an evaluation of the quality of the upper bounds obtained. We also study the bifurcations from the ball as the boundary parameter becomes large (negative).Comment: 26 pages, 20 figure

Numerical minimization of dirichlet laplacian eigenvalues of four-dimensional geometries

We develop the first numerical study in four dimensions of optimal eigenmodes associated with the Dirichlet Laplacian. We describe an extension of the method of fundamental solutions adapted to the four-dimensional context. Based on our numerical simulation and a postprocessing adapted to the identification of relevant symmetries, we provide and discuss the numerical description of the eighth first optimal domains.The work of the first author was partially supported by FCT, Portugal, through the program “Investigador FCT” with reference IF/00177/2013 and the scientific project PTDC/MATCAL/4334/2014. The work of the second author was supported by the ANR, through the projects COMEDIC, PGMO, and OPTIFORMinfo:eu-repo/semantics/publishedVersio

Asymptotic behaviour and numerical approximation of optimal eigenvalues of the Robin Laplacian

We consider the problem of minimising the $n^{th}-$eigenvalue of the Robin Laplacian in $\mathbb{R}^{N}$. Although for $n=1,2$ and a positive boundary parameter $\alpha$ it is known that the minimisers do not depend on $\alpha$, we demonstrate numerically that this will not always be the case and illustrate how the optimiser will depend on $\alpha$. We derive a Wolf-Keller type result for this problem and show that optimal eigenvalues grow at most with $n^{1/N}$, which is in sharp contrast with the Weyl asymptotics for a fixed domain. We further show that the gap between consecutive eigenvalues does go to zero as $n$ goes to infinity. Numerical results then support the conjecture that for each $n$ there exists a positive value of $\alpha_{n}$ such that the $n^{\rm th}$ eigenvalue is minimised by $n$ disks for all $0<\alpha<\alpha_{n}$ and, combined with analytic estimates, that this value is expected to grow with $n^{1/N}$

Analysis of a class of boundary value problems depending on left and right Caputo fractional derivatives

In this work we study boundary value problems associated to a nonlinear fractional ordinary differential equation involving left and right Caputo derivatives. We discuss the regularity of the solutions of such problems and, in particular, give precise necessary conditions so that the solutions are C1([0, 1]). Taking into account our analytical results, we address the numerical solution of those problems by the augmented-RBF method. Several examples illustrate the good performance of the numerical method.P.A. is partially supported by FCT, Portugal, through the program “Investigador FCT” with reference IF/00177/2013 and the scientific projects PEstOE/MAT/UI0208/2013 and PTDC/MAT-CAL/4334/2014. R.F. was supported by the “Fundação para a Ciência e a Tecnologia (FCT)” through the program “Investigador FCT” with reference IF/01345/2014.info:eu-repo/semantics/publishedVersio

On the behavior of clamped plates under large compression

We determine the asymptotic behavior of eigenvalues of clamped plates under large compression by relating this problem to eigenvalues of the Laplacian with Robin boundary conditions. Using the method of fundamental solutions, we then carry out a numerical study of the extremal domains for the first eigenvalue, from which we see that these depend on the value of the compression, and start developing a boundary structure as this parameter is increased. The corresponding number of nodal domains of the first eigenfunction of the extremal domain also increases with the compression.This work was partially supported by the Funda ̧c ̃ao para a Ciˆencia e a Tecnologia(Portugal) through the program “Investigador FCT” with reference IF/00177/2013 and the projectExtremal spectral quantities and related problems(PTDC/MAT-CAL/4334/2014).info:eu-repo/semantics/publishedVersio

Detection of holes in an elastic body based on eigenvalues and traces of eigenmodes

We consider the numerical solution of an inverse problem of finding the shape and location of holes in an elastic body. The problem is solved by minimizing a functional depending on the eigenvalues and traces of corresponding eigenmodes. We use the adjoint method to calculate the shape derivative of this functional. The optimization is performed by BFGS, using a genetic algorithm as a preprocessor and the Method of Fundamental Solutions as a solver for the direct problem. We address several numerical simulations that illustrate the good performance of the method.info:eu-repo/semantics/publishedVersio

Determination of elastic resonance frequencies and eigenmodes using the method of fundamental solutions

In this paper, we present the method of fundamental solutions applied to the determination of elastic resonance frequencies and associated eigenmodes. The method uses the fundamental solution tensor of the Navier equations of elastodynamics in an isotropic material. The applicability of the the method is justified in terms of density results. The accuracy of the method is illustrated through 2D numerical examples for the disk and non trivial shapes.info:eu-repo/semantics/publishedVersio

Numerical optimisation of Dirac eigenvalues

Motivated by relativistic materials, we develop a numerical scheme to support existing or state new conjectures in the spectral optimisation of eigenvalues of the Dirac operator, subject to infinite-mass boundary conditions. We study the optimality of the regular polygon (respectively, disk) among all polygons of a given number of sides (respectively, arbitrary sets), subject to area or perimeter constraints. We consider the three lowest positive eigenvalues and their ratios. Roughly, we find results analogous to known or expected for the Dirichlet Laplacian, except for the third eigenvalue which does not need to be minimised by the regular polygon (respectively, the disk) for all masses. In addition to the numerical results, a new, mass-dependent upper bound to the lowest eigenvalue in rectangles is proved and its extension to arbitrary quadrilaterals is conjectured.Comment: 19 pages, 26 figure