arrow
Return

Eigenvalue problem derivatives computation for a complex matrix using the adjoint method

delete2023-02-01
delete10
PRE
AI
S
Sicheng He *
Y
Yayun Shi
E
Eiríkur Jónsson
J
Joaquim R. R. A. Martins
DOI:10.1016/j.ymssp.2022.109717delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
Derivatives of eigenvalues and eigenvectors with respect to design variables are required for gradient-based optimization in many engineering design problems. However, for the general-ized and standard eigenvalue problems with general complex and non-Hermitian coefficient matrices, no method can accurately compute the eigenvalue and eigenvector derivatives while remaining efficient for large numbers of design variables. In this paper, we develop an adjoint method to compute complex eigenvalue and eigenvector derivatives with machine precision. For the special case when only the eigenvalue derivative is required, we propose a reverse algorithmic differentiation (RAD) formula using a newly developed dot product identity for complex functions. We verify the proposed method against the finite differences (FD) for a simple algebraic example with a 3-by-3 complex non-Hermitian matrix and a plane Poiseulle flow stability problem that is modeled as a generalized eigenvalue problem. The adjoint method is demonstrated to scale well with the number of design variables, matching the FD reference to about 5 to 7 digits.
Keywords:
Adjoint method
Eigenvalue problem
Sensitivity analysis
Algorithmic differentiation

Journal

Mechanical Systems and Signal Processing cover
Mechanical Systems and Signal Processing
IF:
8.9
Papers:
1.3W
Citations:
6.6W

Organization

U
university of michigan system
Scholars:
9.1W
Papers: 8.6W
Citations: 133
Cited Papers

Cited Papers

3-D high-fidelity hydrostructural optimization of cavitation-free composite lifting surfaces
err2021-07-01
err17
errOAAI
errLiao, Yingqian; Martins, Joaquim R. R. A.; Young, Yin Lu
errShare
errSave
Editorial
err2000-06-01
err0
errOAAI
errDavid Thompson
errShare
errSave
researcher View more