Eigenvalue
Solve a PDE eigenvalue problem.
Syntax
model.sol(sname).create(fname,"Eigenvalue")
model.sol(sname).feature(fname).set(pname,value)
 
model.sol(sname).feature(fname).create(fname2,LinearType)
 
model.sol(sname).feature(fname).create(fname2,"Advanced")
 
Here LinearType is any of the allowed linear solver feature types.
Description
Operation feature.
For both linear and nonlinear problems, the eigenvalue problem is that of the linearization about a solution U0. If the eigenvalue appears nonlinearly, COMSOL Multiphysics reduces the problem to a quadratic approximation around a value λ0 specified by the property eigref. The discretized form of the problem reads
where K, D, E, N, and NF are evaluated for U = U0 and λ = λ0. Λ is the Lagrange multiplier vector, and λ is the eigenvalue. The eigenvalue name can be given by the property eigname. The linearization point U0 can be given with the property U. The shift, described below, is compensated according to the linearization point for the eigenvalue. Therefore, changing the linearization point has no effect at all for linear or quadratic eigenvalue problems. The eigenvalue search method can be manual or a region in the complex plane (controlled by the property eigmethod) for the ARPACK eigenvalue solver or by stochastic estimation or a manually defined half or full contour in the complex plane for the FEAST eigenvalue solver.
The eigenvalue feature accepts the following properties:
true | false
Perform eigenvalue region consistency check when eigsolver is arpack and eigmethod is region.
auto | on | off
Rotation angle of the ellipse contour (180 – 180 degrees), when eigsolver is feast and searchcon is wholecon.
Lower bound of search interval, when eigsolver is feast and searchcon is halfcon.
Center of the ellipse contour, when eigsolver is feast and searchcon is wholecon.
Horizontal radius of the ellipse contour, when eigsolver is feast and searchcon is wholecon.
Imaginary semiaxis, when eigsolver is feast and searchcon is halfcon.
Name of the controlling study step or user if the feature is controlled manually.
Upper bound of search interval, when eigsolver is feast and searchcon is halfcon.
true | false
true | false
manual | region
Eigenvalue search method when eigsolver is arpack: around a shift or using a rectangular region.
true | false
arpack | feast | lapack
true | false
average | maximum | mass
auto | gaussf | trapez | zolot
Integration type for estimation: automatic, Gauss, or trapezoidal, when eigsolver is feast. Also, you can choose Zolotarev (zolot) when searchcon is halfcon.
on | off
init | solution
zero | solution object
current | solution store
Maximum matrix size, when eigsolver is set to lapack.
Approximate number of eigenvalues, when neigsmethodnhm is set to manual.
stochest | manual
auto | manual
auto | manual
Number of integration points, when nintgptsmethodhm is manual.
halfcon | wholecon
Specify the eigenvalue search contour when eigsolver is feast: a whole contour or a half contour (Hermitian problem).
Eigenvalue search location when eigsolver is arpack or lapack.
on | off
Eigenvalue shift (when shiftregselect = manual).
auto | manual
auto | all | positive integer
true | false
on | off
relerrtrace | maxrelres
true | false
Enter transformed values, when eigsolver is feast and eigmethod is manual.
true | false
true | false
auto | false
Specify where to look for the desired eigenvalues with the property shift. Enter a real or complex scalar; the default value is 0, meaning that the solver tries to find eigenvalues close to 0.
For more information about the eigenvalue solver, see Eigenvalue Solver in the COMSOL Multiphysics Reference Manual.