[Getdp] coupled electrical field - particle densities
Kubicek Bernhard
Bernhard.Kubicek at arsenal.ac.at
Wed Nov 15 18:08:48 CET 2006
Hello,
once again.
As promised, there is now a wiki-example for something similar: http://www.geuz.org/getdp/wiki/Avalanche2d
nice greetings,
Bernhard
-----Ursprüngliche Nachricht-----
Von: getdp-bounces at geuz.org [mailto:getdp-bounces at geuz.org] Im Auftrag von Reinhold Ingo
Gesendet: Dienstag, 14. November 2006 18:28
An: getdp at geuz.org
Betreff: [Getdp] coupled electrical field - particle densities
Dear all,
I'm not too familiar with GetDP and hope to get some answers by
posting my problem here.
I'm trying to compute a problem with coupled electrical field and
particle densities. From the literature it is suggested that for
example electron densities in partial discharges are computed by
something like
dt(ne)=div(D grad(ne))-div( ve*ne)+C*ne
The electrical field is then
div( epsr grad(phi))-e/eps0*ne=0
However what confuses me now is, that when I'm trying to compute the
problem in GetDP I run into a lot of problems. Maybe I'm doing
something wrong in the formulation.
If I'm right the second equation gives a weak formulation of
Galerkin { [ epsr[] * Dof{d v} , {d v} ]; In Volume; Jacobian Vol;
Integration Int; }
GlobalTerm { [ -elemc*Dof{NE}/(eps0) , {V} ]; In Boundary;}
Since NE is defined as a GlobalQuantity my problem of understanding
is now how to arrange the first equation in such a way, that both
work together in a non-linear algorithm, even though the density is
in the plane, while NE in the above given formula is only at the
boundary and set up in a function space like this
FunctionSpace {
{ Name Elec; Type Form0;
BasisFunction {
{ Name sn ; NameOfCoef vn ; Function BF_Node ; Support Region
[{Volume}] ; Entity NodesOf[ All,Not Boundary] ; }
{ Name sf ; NameOfCoef vfu ; Function BF_GroupOfNodes ; Support
Volume ; Entity GroupsOfNodesOf[Boundary] ; }
}
GlobalQuantity {
{ Name GlobalElectricPotential ; Type AliasOf ; NameOfCoef
vfu ; }
{ Name GlobalElectronDensity ; Type AssociatedWith ; NameOfCoef
vfu; }
}
Constraint {
{ NameOfCoef GlobalElectronDensity ; EntityType
GroupsOfNodesOf ; NameOfConstraint ElectronDensity; }
{ NameOfCoef GlobalElectricPotential ; EntityType
GroupsOfNodesOf ; NameOfConstraint ElectricScalarPotential ; } }
}
}
My guess was that the electron density is then calculated using
Galerkin { Dt[ Dof{NE} , {NE} ]; In Volume; Jacobian Vol;
Integration Int; }
Galerkin { [ Dof{NE}*we , {d ne} ]; In Volume; Jacobian Vol;
Integration Int; }
Galerkin { [ -De*Dof{d ne} , {d ne} ]; In Volume; Jacobian Vol;
Integration Int; }
Galerkin { [ -C*Dof{NE} , {NE} ]; In Volume; Jacobian Vol;
Integration Int; }
but didn't succeed.
It would be really great if someone could give me hint how to
formulate it.
Thanks in advance and nice greetings from Chemnitz,
Ingo
_______________________________________________
getdp mailing list
getdp at geuz.org
http://www.geuz.org/mailman/listinfo/getdp
More information about the getdp
mailing list