[Getdp] Method of moments
Ruth V. Sabariego
r.sabariego at ulg.ac.be
Thu Jan 17 11:44:53 CET 2008
Dear Lars,
You do have to include the Dof in the definition of the integral
quantity. E is a degree of freedom, so you have to indicate it even if
appears in an integral quantity.
The 'so-called' Dof flag is not necessary if you are going to use the
integral quantity in a postprocessing.
{ Name form_dirichlet; Type FemEquation;
Quantity {
{ Name E; Type Local; NameOfSpace nodal_form0;}
{ Name E2; Type Integral ; NameOfSpace nodal_form0;
[ Helmholtz[]{2D,k0} * Dof{E} ]; In scat; Jacobian Jac ;
Integration I1 ; }
}
Equation {
Galerkin { [ Dof{E},{E} ]; In Domain; Integration I1;
Jacobian Jac;}
Galerkin { [ -inc[],{E} ]; In Domain; Integration I1;
Jacobian Jac;}
Galerkin { [ Dof{E2},{E} ]; In Domain2D; Integration I1;
Jacobian Jac;}
/* g d^2 h - g k^2 h = -dg*dh - k^2 g*h*/
}
}
I've just tried your files with my Getdp and Dof in the integral
quantity definition mand they work here.
Regards,
Ruth
Lars Rindorf wrote:
> Dear Olivier
>
> I tried to do this, but the message consol returns
> GetDP : Operation : Generate[A]
> GetDP : Solver : Loading parameter file 'C:\Documents and Settings\lhr\My Documents\getPD\Greens_sq\solver.par'
> GetDP : Info : Setting System {A,b} to zero
> GetDP : Error : Segmentation violation (invalid memory reference)
>
> In the getpd manual it is describes that the integral type is assembled before the calculation, and thus can't be a degree of freedom (dof). Before the calculation the field is zero, and thus the integral term is effectively ignored.
>
> But I don't really know what to do from here. I have attached the .pro, .geo, and solver paramter files, so you can play with it. Note, that when you postproces the scattered field, i.e. the integral with the green's function, the computed E-field is used.
>
> Thanks
> Lars
>
>
>
> ________________________________
>
> From: Olivier Castany [mailto:castany at quatramaran.ens.fr]
> Sent: Wed 1/16/2008 11:12 PM
> To: Lars Rindorf; getdp at geuz.org
> Subject: Re: [Getdp] Method of moments
>
>
>
>
>> { Name form_neumann; Type FemEquation;
>> Quantity {
>> { Name E; Type Local; NameOfSpace scalar;}
>> { Name E2; Type Integral ; NameOfSpace scalar;
>> [ Helmholtz[]{2D,k0} * {E} ];
>> In feed; Jacobian Jac ; Integration I1 ; }
>> }
>> Equation {
>> Galerkin { [ Dof{E},{E} ]; In Domain; Integration I1; Jacobian Jac;}
>> Galerkin { [ -inc[],{E} ]; In Domain; Integration I1; Jacobian Jac;}
>> Galerkin { [ {E2},{E} ]; In Domain2D; Integration I1; Jacobian Jac;}
>> }
>> }
>> The problem is 2D. "feed" is a 1D boundary. "inc[]" is the incoming field. Domain is the union of the 2D domain (Domain2D) and the 1D domain (feed).
>>
>> The problem is that the Green's function integral E2 is evaluated once during initialization of the problem, and it is not included in the discretization
>> of the matrix. I have tried making a coupled model. This seems to be able to work, but it is not very elegant from a mathematical point of view. I have
>> come rather stuck in the problem.
>>
>
>
> I think you should add the keyword "Dof" in the definition of the integral quantity and in the equations :
>
> { Name form_neumann; Type FemEquation;
> Quantity {
> { Name E; Type Local; NameOfSpace scalar;}
> { Name E2; Type Integral ; NameOfSpace scalar;
> [ Helmholtz[]{2D,k0} * Dof{E} ];
> In feed; Jacobian Jac ; Integration I1 ; }
> }
> Equation {
> Galerkin { [ Dof{E},{E} ]; In Domain; Integration I1; Jacobian Jac;}
> Galerkin { [ -inc[],{E} ]; In Domain; Integration I1; Jacobian Jac;}
> Galerkin { [ Dof{E2},{E} ]; In Domain2D; Integration I1; Jacobian Jac;}
> }
> }
>
> Let me know if it works, I'm interested.
>
> Bye,
>
> O.C.
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp
>
--
Dr. Ir. Ruth V. Sabariego
University of Liege, Institut Montefiore,
Dept. of Electrical Engineering & Computer Science,
Applied & Computational Electromagnetics (ACE),
Sart Tilman Campus, Grande Traverse, 10 (B28), B-4000 LIEGE, Belgium
phone: +32-4-3663737 -- fax: +32-4-3662910 -- http://elap.montefiore.ulg.ac.be/