Parallel Computation of Acoustic Eigenvalue Problems Using a Polynomial Jacobi-Davidson Method

以作者查詢圖書館館藏

、以作者查詢臺灣博碩士

、以作者查詢全國書目

、勘誤回報

、線上人數：23

、訪客IP：3.144.91.223

姓名

賴聲泓(Sheng-Hong Lai) 查詢紙本館藏

畢業系所

數學系

論文名稱

(Parallel Computation of Acoustic Eigenvalue Problems Using a Polynomial Jacobi-Davidson Method)

相關論文

★ 非線性塊狀高斯消去牛頓演算法在噴嘴流體的應用	★ 以平行 Newton-Krylov-Schwarz 演算法解 Poisson-Boltzmann 方程式的有限元素解在膠體科學上的應用
★ 最小平方有限元素法求解對流擴散方程以及使用Bubble函數的改良	★ Bifurcation Analysis of Incompressible Sudden Expansion Flows Using Parallel Computing
★ Parallel Jacobi-Davidson Algorithms and Software Developments for Polynomial Eigenvalue Problems in Quantum Dot Simulation	★ An Inexact Newton Method for Drift-DiffusionModel in Semiconductor Device Simulations
★ Numerical Simulation of Three-dimensional Blood Flows in Arteries Using Domain Decomposition Based Scientific Software Packages in Parallel Computers	★ A Parallel Fully Coupled Implicit Domain Decomposition Method for the Stabilized Finite Element Solution of Three-dimensional Unsteady Incompressible Navier-Stokes Equations
★ A Study for Linear Stability Analysis of Incompressible Flows on Parallel Computers	★ Numerical Study of Algebraic Multigrid Methods for Solving Linear/Nonlinear Elliptic Problems on Sequential and Parallel Computers
★ A Parallel Multilevel Semi-implicit Scheme of Fluid Modeling for Numerical Low-Temperature Plasma Simulation	★ Performance Comparison of Two PETSc-based Eigensolvers for Quadratic PDE Problems
★ A Parallel Two-level Polynomial Jacobi-Davidson Algorithm for Large Sparse Dissipative Acoustic Eigenvalue Problems	★ A Full Space Lagrange-Newton-Krylov Algorithm for Minimum Time Trajectory Optimization
★ Parallel Two-level Patient-specific Numerical Simulation of Three-dimensional Rheological Blood Flows in Branching Arteries	★ A Markov Chain Multi-elimination Preconditioner for Elliptic PDE Problems on GPU

檔案

[Endnote RIS 格式]

[Bibtex 格式]

[相關文章]

[文章引用]

[完整記錄]

[館藏目錄]

[檢視]

[下載]

本電子論文使用權限為同意立即開放。
已達開放權限電子全文僅授權使用者為學術研究之目的，進行個人非營利性質之檢索、閱讀、列印。
請遵守中華民國著作權法之相關規定，切勿任意重製、散佈、改作、轉貼、播送，以免觸法。

摘要(中)

在我們日常生活周遭隨處可遇到有關於聲音的問題。例如開車時、搭公車時或搭飛機時。這些問題可以經由數學模型來描述聲音在一些吸振材質加入之下的振動情況，再經由有限元素法離散之後可以得到一個多項式型的特徵值問題。在部份應用中，我們有興趣找出一些低頻率的特徵值，但這些特徵值通常都落在頻譜的內部。並且當網格切的非常細的時後，特徵值問題的維度會變成相當地大。所以我們需要一個可平行化的工具來解這種大型稀疏的多項式特徵值問題。Jacobi-Davidson 演算法提供了一個快速並且有效率的方式來解出這類問題的內部特徵值。我們在 additive Schwarz 的架構下平行化實作 acobi-Davidson 演算法，並且用來解由聲音的問題所產生的多項式型的特徵值問題。經由數值實驗結果，我們列出了一些 additive Schwarz preconditioned Jacobi-Davidson 的平行效能。在解 correction equation 時，經由 Krylov-Schwarz 演算法的幫助之下，Jacobi-Davidson 的效率有著顯著的進步。

摘要(英)

The acoustic problems usually happens around us in our daily life when we drive a car, take a bus or take a plane. From the problems of acoustic vibrations with damping, a polynomial eigenvalue problem is obtained by applying the Galerkin finite element method. For particular applications, we are interested in finding some selected low frequency eigenvalues which are located within the interior of the spectrum. The size of the resulting eigenproblem is typically large especially for very fined mesh case so that the parallel polynomial eigensolver is need to deal with such problem. The Jacobi-Davidson method provides a fast and efficient manner for solving the interior eigenvalues for the large sparse polynomial eigenvalue problems. We proposed an Jacobi-Davidson method based on an additive Schwarz framework in parallel implementation and used it to solve the polynomial eigenvalue problem arising from the acoustic. And we showed some parallel performance of the additive Schwarz preconditioned Jacobi-Davidson method by numerical experiments. With help of Krylov-Schwarz algorithm for the correction equation, the efficiency of JD algorithm is greatly improved.

關鍵字(中)

★ 平行計算
★ 多項式型
★ 特徵值
★ 有限元素法

關鍵字(英)

★ finite element method
★ acoustic
★ eigenvalue
★ precondition
★ additive Schwarz
★ Jacobi-Davidson method

論文目次

List of Figures . . . . . . . . . . . . . . . . . . . . . . . . vi
List of Tables . . . . . . . . . . . . . . . . . . . . . . . . vii
1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . 1
2 Polynomial acoustic eigenvalue problems . . . . . . . . . . . . 4
2.1 Mathematical acoustic model . . . . . . . . . . . . . . . . 4
2.1.1 The continuity equation and Euler’s equation . . . . . 4
2.1.2 Acoustic wave equation . . . . . . . . . . . . . . . . 6
2.1.3 Boundary conditions and time-separation assumption . . 7
2.2 Galerkin ﬁnite element discretization . . . . . . . . . . . 9
2.3 Some acoustic polynomial eigenvalue problems . . . . . . . 10
2.4 A special case with semi-analytical solution available . . 12
3 Parallel solution algorithms for polynomial eigenvalue problem 17
3.1 Additive Schwarz preconditioned Jacobi-Davidson algorithm 17
3.1.1 Preconditioned correction equation . . . . . . . . . . 20
3.1.2 Restricted additive Schwarz . . . . . . . . . . . . . 21
4 Numerical experiment of ASPJD . . . . . . . . . . . . . . . . 22
4.1 A benchmark problem and setup for numerical experiments . 22
4.2 Parallel code validation . . . . . . . . . . . . . . . . . 24
4.3 ASPJD tunning . . . . . . . . . . . . . . . . . . . . . . 28
4.4 Parallel scalability analysis . . . . . . . . . . . . . . 31
4.5 Comments on other possible correction equation solvers . . 33
4.6 Linearized approaches . . . . . . . . . . . . . . . . . . 33
5 Concluding remarks . . . . . . . . . . . . . . . . . . . . . . 37
Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . 38
Appendix . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
A. Bilinear basis and the local stiffess and mass matrices . . 41

參考文獻

[1] Satish Balay, Kris Buschelman, William D. Gropp, Dinesh Kaushik, Matthew G. Knepley, Lois Curfman McInnes, Barry F. Smith, and Hong Zhang. PETSc Web page, 2001. http://www.mcs.anl.gov/petsc.
[2] A. Bermudez, R. G. Duran, and J. Rodriguez, R. Solomin. Finite element analysis of a quadratic eigenvalue problem arising in dissipative acoustics. SIAM J. Numer. Anal., 38:267–291, 2000.
[3] J.G.L. Booten, H.A. van der Vorst, P.M. Meijer, and H.J.J. te Riele. A Preconditioned Jacobi-Davidson Method for Solving Large Generalized Eigenvalue Problems. Centrum voor Wiskunde en Informatica, 1994.
[4] T.-M. Huang W. Wang F.-N. Hwang, Z.-H. Wei. A parallel additive Schwarz preconditioned Jacobi-Davidson algorithm for polynomial eigenvalue problems in quantum dot simulation. Journal of Computational Physics, 229:2932 – 2947, 2010.
[5] D.R. Fokkema, G.L.G. Sleijpen, and H.A. van der Vorst. Jacobi-Davidson style QR and QZ algorithms for the reduction of matrix pencils. SIAM J. Sci. Comput., 20:94–125, 1998.
[6] M.B. Van Gijzen. The parallel computation of the smallest eigenpair of an acoustic problem with damping. Int. J. Numer. Meth. Engng., 45:765–777, 1999.
[7] D.R. Fokkema H.A. van der Vorst G.L.G. Sleijpen, A.G.L. Booten. Jacobi-Davidson type methods for generalized eigenproblems and polynomial eigenproblems. BIT, 36:595–633, 1996.
[8] J.-S. Guo, W.-W. Lin, and C.-S. Wang. Numerical solutions for large sparse quadratic eigenvalue problems. Linear Algebra Appl., 225:57–89, 1995.
[9] V. Hernandez, J.E. Roman, A. Tomas, and V. Vidal. SLEPc home page. http://www.grycap.upv.es/slepc, 2006.
[10] M. Hochbruck and D. Lochel. A multilevel Jacobi-Davidson method for polynomial pde eigenvalue problems arising in plasma physics. SIAM J. Sci. Comput., 06 2010.
[11] T.J.R. Hughes. The Finite Element Method: Linear Static and Dynamic Finite Element Analysis. Dover Publications, 2000.
[12] M. Nool and A. van der Ploeg. A parallel Jacobi-Davidson-type method for solving large generalized eigenvalue problems in magnetohydrodynamics. SIAM J. Sci. Comput., 22:95–112, 2000.
[13] A. Bermudez R. Rodriuez. Modelling and numerical solution of elastoacoustic vibrations with interface damping. Int. J. Numer. Meth. Engng., 46:1763–1779, 1999.
[14] Y. Saad. Numerical methods for large eigenvalue problems. Manchester University Press, 1992.
[15] G.L.G. Sleijpen and H.A. van der Vorst. A Jacobi-Davidson iteration method for linear eigenvalue problems. SIAM J. Matrix Anal. Appl., 17:401–425, 1996.
[16] G.L.G. Sleijpen and H.A. van der Vorst. A Jacobi-Davidson iteration method for linear eigenvalue problems. SIAM Rev., 42(2):267–293, 2000.
[17] T. Betcke, N.J. Higham, V. Mehrmann, C. Schroder, F. Tisseur. NLEVP: A collection of nonlinear eigenvalue problems. April 2008.
[18] J.-L. Liu W. Wang T.-M. Hwang, W.-W. Lin. Jacobi-Davidson methods for cubic eigenvalue problems. Numer. Linear Algebra Appl., 12:605–624, 2005.
[19] W.-C. Wang W. Wang T.-M. Hwang, W.-W. Lin. Numerical simulation of three dimensional pyramid quantum dot. J. Comput. Phys., 196:208–232, 2004.
[20] F. Tisseur. Backward error and condition of polynomial eigenvalue problems. Linear Algebra Appl., 309:339–361, 2000.
[21] Meerbergen K. Tisseur, F. The quadratic eigenvalue problem. SIAM Rev., 43(2):235–286, 2001.
[22] H. Voss. A Jacobi-Davidson method for nonlinear and nonsymmetric eigenproblems. Comput. & Structures, 85(17-18):1284–1292, 2007.
[23] W.-W. Lin J.-L. Liu W. Wang, T.-M. Hwang. Numerical methods for semiconductor heterostructures with band nonparabolicity. J. Comput. Phys., 190:141–158, 2003.
[24] M. Sarkis X.-C. Cai. A restricted additive Schwarz preconditioner for general sparse linear systems. SIAM J. Sci. Comput, 21:792–797, 1999.

指導教授

黃楓南(Feng-Nan Hwang)

審核日期

2010-7-28

推文