[Getdp] 2D Helmholtz eigenvalue equation
Lars Rindorf
lhr at com.dtu.dk
Mon Oct 22 16:06:47 CEST 2007
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>