[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>