[Getdp] Circular cavity
Jakub Tesař
krlesch at seznam.cz
Mon Apr 12 10:48:17 CEST 2010
Dear GetDP developers,
I´d like to solve resonance frequency and electric field in loaded circular cavity with an EigenSolve procedure. The .pro file is bellow. I became this request:
Frequency= 0Hz
and
No system available for EigenSolve: check `DtDt` and `GenerateSeparate`
Can somebody help me?
Thank you for your answer
Jakub Tesar, Masaryk univerzity in Brno, Czech republic
_______________________________________________________
Group {
Wall = Region[29];
Vol = Region[30];
Sur = Region[{Wall}] ;
Tot = Region[{Vol,Wall}] ;
Vol_S = Region[{Wall}] ;
Vol_C = Region[{Vol}] ;
}
Function {
nu0 = 1 / (4.e-7 * Pi) ;
ep0 = 8.854187817e-12 ;
nu [ Vol ] = nu0 ;
epsilon [ Vol ] = ep0 * 1.00058;
sigma [ Vol ] = 0 ;
}
Constraint {
{ Name PEC ;
Case {
{ Region Wall; Value 0. ; }
}
}
}
Include "IntJac.pro"
Function{
DefineFunction[epsilon, sigma, nu, J_s, K_s, NxH];
DefineVariable[ Val_Rint, Val_Rext, Freq, DFreq ] ;
}
Group{
DefineGroup[ Vol, Vol_C, Vol_S, Sur, Pml] ;
}
FunctionSpace {
{ Name Hcurl_e_3D; Type Form1;
BasisFunction {
{ Name se; NameOfCoef ee; Function BF_Edge; Support Tot; Entity EdgesOf[All]; }
}
Constraint {
{ NameOfCoef ee; EntityType EdgesOf ; NameOfConstraint PEC; }
}
}
}
Formulation {
{ Name MW_e_3D; Type FemEquation;
Quantity {
{ Name e; Type Local; NameOfSpace Hcurl_e_3D; }
}
Equation {
Galerkin { DtDt [ epsilon[] * Dof{e} , {e} ];
In Vol; Integration I1; Jacobian JVol; }
Galerkin { [ nu[] * Dof{d e} , {d e} ];
In Vol; Integration I1; Jacobian JVol; }
Galerkin { Dt [ J_s[] , {e} ];
In Vol_S; Integration I1; Jacobian JVol; }
}
}
}
Resolution {
{ Name MW_e_3D;
System {
{ Name A; NameOfFormulation MW_e_3D; Type ComplexValue; Frequency Freq;}
}
Operation {
GenerateSeparate [A]; EigenSolve[A,10,1,0]; SaveSolutions[A];
}
}
}
PostProcessing {
{ Name MW_e_3D; NameOfFormulation MW_e_3D; NameOfSystem A;
Quantity {
{ Name e; Value{ Local{ [ {e} ] ; In Vol; } } }
{ Name e2; Value{ Local{ [ SquNorm[{e}] ] ; In Vol; } } }
}
}
}
PostOperation {
{ Name MW_e_3D; NameOfPostProcessing MW_e_3D;
Operation {
Print[ e, OnElementsOf Vol , File "e.pos" ] ;
Print[ e2, OnElementsOf Vol , File "e2.pos" ] ;
}
}
}
-------------- next part --------------
A non-text attachment was scrubbed...
Name: volny.geo
Type: text/plain charset=us-ascii
Size: 1137 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20100412/57c9eeb8/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: volny.msh
Type: text/plain charset=us-ascii
Size: 4317625 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20100412/57c9eeb8/attachment.msh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: volny.pro
Type: text/plain charset=us-ascii
Size: 2097 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20100412/57c9eeb8/attachment.pro>