[Getdp] 2D Helmholtz eigenvalue equation

[Lars Rindorf]

Dear all
 
I'm trying to model a Helmholtz equation on a circular domain in 2D. The script is adapted from the 1eigen script from the getdp webpage where the same equation was solved in 1D. The model works fine, but when I try to impose Dirichlet boundary conditions, then the bc is not enforced as the field is non-zero at the boundary. 
 
What is wrong in my script?
 
Kind regards
Lars Rindorf

_________________________________________________

Lars Rindorf, M. Sc., PhD student

Fibers and Nonlinear Optics Group

COM*DTU Department of Communication, Optics and Materials

Technical University of Denmark

Building 345v/267

DK-2800 Lyngby

Denmark

Phone: +45 45256365

Fax: +45 45936581

Email: LHR at com.dtu.dk

 

-------------- next part --------------
A non-text attachment was scrubbed...
Name: eigen2D.pro
Type: application/octet-stream
Size: 2931 bytes
Desc: eigen2D.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20071022/87f44b52/attachment.pro>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: solver.par
URL: <http://www.geuz.org/pipermail/getdp/attachments/20071022/87f44b52/attachment.par>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: eigen2D_2.geo
Type: application/octet-stream
Size: 639 bytes
Desc: eigen2D_2.geo
URL: <http://www.geuz.org/pipermail/getdp/attachments/20071022/87f44b52/attachment.geo>