[Getdp] 2D Helmholtz eigenvalue equation

[Bernhard Kubicek]

Hello!
That's quite a very nice .pro file. But your error happened in the .geo
file.
The physical line describing the boundary was set to the Line Loop(3).
Instead you have to use the individual circle arcs:
Physical Line(1001) = {1,2};
After remeshing, everything works nicely.
I am just fooling around to simulate a guitar top using your script :)

Very nice greetings,
Bernhard

On Mon, 2007-10-22 at 16:06 +0200, Lars Rindorf wrote:
> 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
> 
>  
> 
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp