COMSOL 6.2 - Eigenvalue

) study and study step are used to compute the eigenvalues and eigenmodes of a linear or linearized model in a generic eigenvalue formulation where the eigenvalues are not necessarily frequencies. The Eigenvalue study gives you full control of the eigenvalue formulation, in contrast to the eigenfrequency study that is adapted for specific physics interfaces. The Eigenvalue study is typically used for equation-based modeling.

Selecting an Eigenvalue study gives a solver configuration with an Eigenvalue Solver.


	The Results While Solving, Physics and Variables Selection, Values of Dependent Variables, Store in Output, Mesh Selection, Adaptation and Error Estimates, and Geometric Entity Selection for Adaptation sections and the Include geometric nonlinearity check box are described in Common Study Step Settings. There is also detailed information in the Physics and Variables Selection and Values of Dependent Variables sections.

Study Settings

From the list, choose (the default), , or :

The ARPACK algorithm is based on an algorithmic variant of an Arnoldi process When the matrix A is symmetric it reduces to a variant of the Lanczos process. For its settings, see Study Settings for ARPACK below.

The FEAST algorithm uses an inverse residual iteration algorithm and seeks to accelerate the convergence of the subspace eigenvalue problem. It can be effective for finding clustered eigenvalues by using an ellipse to search within. For its settings, see Study Settings for FEAST.

The LAPACK algorithm is useful for finding all eigenvalues for a filled matrix. This option is only applicable for small eigenvalue problems. For its settings, see Study Settings for LAPACK (Filled Matrix).

For more information about these eigenvalue solvers, see The Eigenvalue Solver Algorithms.

Study Settings for ARPACK

From the list, select a search method:

•

(the default), to specify some search criteria manually. See Around Shift Eigenvalue Search Settings below.

•

, to define an eigenvalue search region in a complex plane. See Rectangular Eigenvalue Search Region Settings below and The Eigenvalue Solver Algorithms.

Around Shift Eigenvalue Search Settings

By default, the physics interfaces suggest a suitable number of eigenvalues to search for. To specify the number of eigenvalues, select the check box in front of the field to specify the number of eigenvalues you want the solver to return (default: 6).


	In a solid 3D model, the first six eigenvalues are typically zero and correspond to the rigid-body modes of a 3D geometry. You may therefore need to specify a larger number of eigenvalues. The largest number of eigenvalues that the solver can compute is equal to the number of unconstrained degrees of freedom minus two.

From the list, choose a suitable unit (default rad/s).

By default, the physics interfaces suggest a suitable value around which to search for eigenvalues. To specify the value to search for eigenvalues around (shift), select the check box in front of the field; you can then specify a value or expression around which the eigenvalue solver should look for solutions to the eigenvalue equation (default: 0).

Use the list to control how the eigenvalue solver searches for eigenvalues around the specified shift value:

•

Select (the default value) to search for eigenvalues that are closest to the shift value when measuring the distance as an absolute value.

•

Select to search for eigenvalues with a larger real part than the shift value.

•

Select to search for eigenvalues with a smaller real part than the shift value.

•

Select to search for eigenvalues with a larger imaginary part than the shift value.

•

Select to search for eigenvalues with a smaller imaginary part than the shift value.

Rectangular Eigenvalue Search Region Settings

Use the field to specify the approximate number of eigenvalues you want the solver to return (default: 20). The value of the will affect the used by the algorithm; see the section of the Eigenvalue Solver. It means that increasing the value of the will increase the memory requirement and the computational time. If the solver indicates that the value of the is smaller than the actual number of eigenvalues in the given region, it will perform a search for more eigenvalues, which increases the computational time; see The Eigenvalue Region Algorithm. Within limits it is often more efficient to provide a too large value of than a too small.

In the field, you can specify a maximum number of eigenvalues to limit the eigenvalue solver’s search for additional eigenvalues (default: 200).

Under , you define a (defaults: rad/s) and the size of the search region for eigenvalues as a rectangle in the complex plane by specifying the , , , and in the respective text fields. The search region also works as an interval method if the and are equal; the eigenvalue solver then only considers the real axis and vice versa.

Symmetry and Consistency Settings

If you have chosen from the list, the check box is available and selected by default to increase confidence that the solver finds all eigenvalues in the search region. The work required for performing the consistency check constitutes a significant part of the total work of the eigenvalue computation.

Eigenvalue computations can be performed with a nonsymmetric solver or, if applicable, a real symmetric solver. From the list, choose (the default) or . For the option there is the option to select the check box. This consistency check increases the computational time and memory requirements but provides a rigorous check of the applicability of the real symmetric solver.


	When using a search for eigenvalues, the nonsymmetric eigenvalue solver is always used, even when the problem to solve is symmetric.


	The eigenvalue solver can in some cases return more than the desired number of eigenvalues (up to twice the desired number). These are eigenvalues that the eigenvalue solver finds without additional computational effort.

Study Settings for FEAST

Using the FEAST eigenvalue solver, you can define an ellipse to explicitly exclude the eigenvalues outside the ellipse. Choosing a suitable search region determines which eigenvalues that will be returned for results analysis.

From the list, select a search method:

•

(the default), to define an eigenvalue search contour as an ellipse in a complex plane. See Ellipse Search Region settings below.

•

, to define an eigenvalue search contour by half of the whole ellipse. See Half Ellipse Search Region below.

Ellipse Search Region

You define the by specifying the (default: rad/s), , , , and in the respective text fields. It is important that the real semiaxis that you specify in the field is large enough to enclose the eigenvalues of interest. In the field, you specify the imaginary semiaxis in a similar way. In the field, specify the rotation angle in degrees from the vertical axis in the range of −180 degrees to 180 degrees.

From the list, select the method for evaluating the number of eigenvalues inside the eigenvalue search ellipse:

•

(the default) to use stochastic estimation to evaluate the number of eigenvalues. After the stochastic estimation finishes, the eigenvalue solver automatically calculates the eigenvalues inside the eigenvalue search contour, using the number of eigenvalues calculated from stochastic estimation as the (default: 6).

•

to specify the number of eigenvalues inside the eigenvalue search contour manually in the field (default: 6).

For an eigenvalue problem for which you know how many eigenvalues there are, the solver can compute the eigenvalues directly. For that case, choose from the list and enter the desired number of eigenvalues in the field. If you do not know the number of eigenvalues in the defined region, there are two cases: You can just estimate the number of eigenvalues without computing them, or you can let the software estimate the number of eigenvalues and then compute them automatically afterward. The former is achieved by clicking the button (

) at the top right of the section. The latter is achieved by choosing from the list, and then click . Both cases require doing stochastic estimation, which needs the setting of a field (default: 6).


	The Stochastic Estimation works as a Compute to Selected step. Because the Stochastic Estimation action does not solve the eigenvalue problem, there will not be any valid solution from this action. So, if the Stochastic Estimation action is done for a single step or the first step in a multistep study, there will not be any solution that can be used for postprocessing and results analysis. If the action is done for a step after the first step in a multistep study, the main solution of the sequence corresponds to the study step right before the eigenstep for which the estimation is done.

There is an option to select the check box. If selected, linear system factorizations are stored from the first FEAST iteration and reused in later iterations.

If the check box is selected in the dialog box, select the check box to run the FEAST eigenvalue solver in parallel. See Running FEAST in a Parallel MPI Mode for more information.

Half Ellipse Search Region

You define the by specifying the (default: rad/s), , , and in the respective text fields.

From the list, select the method for evaluating the number of eigenvalues inside the eigenvalue search contour:

•

to specify the number of eigenvalues inside the eigenvalue search contour manually in the field (default: 6).

In the field (default: 6), specify the initial guess of the search subspace dimension, which can be interpreted as an initial guess for the number of eigenvalues inside the half contour. This setting is only available when using a stochastic estimation.

There is an option to select the check box. If selected, linear system factorizations are stored from the first FEAST iteration and reused in later iterations.

If the check box is selected in the dialog box, select the check box to run the FEAST eigenvalue solver in parallel. See Running FEAST in a Parallel MPI Mode for more information.

Symmetry and Consistency Settings

For the eigenvalue search method, eigenvalue computations can be performed with a nonsymmetric solver or, if applicable, a real symmetric solver. From the list, choose (the default) or . For the option, and for the eigenvalue search method, there is the option to select the check box. This consistency check increases the computational time and memory requirements but provides a rigorous check of the applicability of the real symmetric solver.

Study Settings for LAPACK (Filled Matrix)

To specify a shift to use in the modes computation, select the check box and then enter a shift (in rad/s) in the associated text field.

In the field, enter an upper limit on the matrix size (default: 2000).

From the list, choose a suitable unit (default: rad/s).

Values of Linearization Point

Use the settings in this section to specify a linearization point.

From the list, choose (the default) to use linearization point settings controlled by the physics interfaces. Choose to specify the linearization point using the list. Select:

•

to use the expressions specified on the nodes under a specific physics interface as a linearization point.

•

to use a solution as a linearization point.

Use the list to specify which solution to use from the available studies. Select:

•

to use a linearization point that is identically equal to zero.

•

Any other available solution to use it as a linearization point. It can be the current solution in the sequence, or a solution from another sequence, or a solution that was stored with the Solution Store node. You select a stored solution by changing to the name of the stored solution. Choose a solution using the list (see under Common Study Step Settings).

Study Extensions

This section contains some optional extensions of the study, such as auxiliary sweeps (see Common Study Step Settings). Adding an auxiliary parametric sweep adds an Eigenvalue Parametric attribute node to the Eigenvalue Solver.

Distribute Parametric Solver

If you are running an auxiliary sweep and want to distribute it by sending one parameter value to each compute node, select the check box. To enable this option, click the button (

) and select in the dialog box.