[Getdp] Helmholtz 2D green function
Guillaume Demésy
guillaume.demesy at fresnel.fr
Wed Mar 5 16:18:38 CET 2014
Dear GetDP users & developers,
I was expecting to get the (Helmholtz) 2D Green function (i.e. exactly i/4 H_0^1(k_0 r)) by solving \laplacian E + k0^2 E=\delta in GetDP.
It seems that my vertex delta condition does not work properly.
My formulation (also attached .geo and .pro) looks like:
FunctionSpace {
{ Name Hgrad; Type Form1P;
BasisFunction {
{ Name sn; NameOfCoef un; Function BF_PerpendicularEdge_1N; Support Omega; Entity NodesOf[Omega]; }
{ Name sn2; NameOfCoef un2; Function BF_PerpendicularEdge_2E; Support Omega; Entity EdgesOf[Omega]; }
}
Constraint {
}
}
}
Formulation {
{
Name helmoltz_scalar; Type FemEquation;
Quantity {
Name u; Type Local; NameOfSpace Hgrad;}
}
Equation {
Galerkin { [k0^2*epsilonr[] * Dof{u},{u} ];
In Omega; Jacobian JVol; Integration Int_1; }
Galerkin {[-1/mur[] *Dof{Curl u} , {Curl u}];
In Omega; Jacobian JVol; Integration Int_1; }
Galerkin { [Vector[0., 0., -1] , {u}];
In Source_point; Jacobian JLin; Integration Int_1; }
}
}
}
Would you have any insight?
Many thanks!
Guillaume Demésy
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140305/70e59e09/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: green2D.geo
Type: application/octet-stream
Size: 2712 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140305/70e59e09/attachment.geo>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140305/70e59e09/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: green2D.pro
Type: application/octet-stream
Size: 5689 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140305/70e59e09/attachment.pro>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140305/70e59e09/attachment-0002.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: parameters_gmsh_getdp.dat
Type: application/octet-stream
Size: 218 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140305/70e59e09/attachment.dat>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20140305/70e59e09/attachment-0003.html>