[Getdp] Method of moments
Olivier Castany
castany at quatramaran.ens.fr
Wed Jan 16 23:12:11 CET 2008
> { 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.