Controllability method for acoustic scattering with spectral elements

被引:9
作者
Heikkola, Erkki
Monkola, Sanna
Pennanen, Anssi
Rossi, Tuomo
机构
[1] Numerola Oy, FI-40101 Jyvaskyla, Finland
[2] Univ Jyvaskyla, Dept Math Informat Technol, FI-40014 Jyvaskyla, Finland
关键词
exact controllability; Helmholtz equation; spectral element method; mass lumping;
D O I
10.1016/j.cam.2006.02.046
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We formulate the Helmholtz equation as an exact controllability problem for the time-dependent wave equation. The problem is then discretized in time domain with central finite difference scheme and in space domain with spectral elements. This approach leads to high accuracy in spatial discretization. Moreover, the spectral element method results in diagonal mass matrices, which makes the time integration of the wave equation highly efficient. After discretization, the exact controllability problem is reformulated as a least-squares problem, which is solved by the conjugate gradient method. We illustrate the method with some numerical experiments, which demonstrate the significant improvements in efficiency due to the higher order spectral elements. For a given accuracy, the controllability technique with spectral element method requires fewer computational operations than with conventional finite element method. In addition, by using higher order polynomial basis the influence of the pollution effect is reduced. (c) 2006 Elsevier B.V. All rights reserved.
引用
收藏
页码:344 / 355
页数:12
相关论文
共 15 条
[1]   Controllability methods for the computation of time-periodic solutions;: application to scattering [J].
Bristeau, MO ;
Glowinski, R ;
Périaux, J .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 147 (02) :265-292
[2]  
Cohen G C., 2002, HIGHER ORDER NUMERIC, DOI DOI 10.1007/978-3-662-04823-8
[3]   RADIATION BOUNDARY-CONDITIONS FOR ACOUSTIC AND ELASTIC WAVE CALCULATIONS [J].
ENGQUIST, B ;
MAJDA, A .
COMMUNICATIONS ON PURE AND APPLIED MATHEMATICS, 1979, 32 (03) :313-357
[4]   A novel multigrid based preconditioner for heterogeneous Helmholtz problems [J].
Erlangga, Y ;
Oosterlee, C ;
Vuik, C .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2006, 27 (04) :1471-1492
[5]   On a class of preconditioners for solving the Helmholtz equation [J].
Erlangga, YA ;
Vuik, C ;
Oosterlee, CW .
APPLIED NUMERICAL MATHEMATICS, 2004, 50 (3-4) :409-425
[6]   A discontinuous Galerkin method with plane waves and Lagrange multipliers for the solution of short wave exterior Helmholtz problems on unstructured meshes [J].
Farhat, C ;
Wiedemann-Goiran, P ;
Tezaur, R .
WAVE MOTION, 2004, 39 (04) :307-317
[7]  
Glowinski R, 2003, MATHEMATICAL AND NUMERICAL ASPECTS OF WAVE PROPAGATION, WAVES 2003, P805
[8]  
HEIKKOLA E, 2006, SCI COMPUTING B, V3
[9]   Algebraic multigrid for higher-order finite elements [J].
Heys, JJ ;
Manteuffel, TA ;
McCormick, SF ;
Olson, LN .
JOURNAL OF COMPUTATIONAL PHYSICS, 2005, 204 (02) :520-532
[10]   The perfectly matched layer for the ultra weak variational formulation of the 3D Helmholtz equation [J].
Huttunen, T ;
Kaipio, JR ;
Monk, P .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2004, 61 (07) :1072-1092