[Getdp] 2d DC Conduction

Laird Bolt lbolt1 at irf.com
Tue Feb 28 22:12:58 CET 2006

I apologise in advance  for sending such a newbie question to the list.  In Payment if someone could
set me on my way I will add my simple example(s) to the main example list.

In an effort to walk (more of a crawl) before I run I am trying to simply
calculate the voltage and current density distribution in a flat bar of copper (no temperature effects assumed)

Fragments for this exist in the mail archives as well as documentation.  However,
I must be doing something wrong as my results are always zero when using Dirichlet BC at both ends
and constant (but non-zero) when I use dirichlet at one end and Neumann BC at the other.

I assume that the null Neumann BC along the long sides is implicit.

I have attached my .geo, .pro and solver.par files for the Dirichlet BC (1V at one end 0v at the other)

THe Formulation was :

Formulation {
  { Name Voltage_Steady ; Type FemEquation ;
    Quantity { 
      { Name Volt ; Type Local ; NameOfSpace Vgrad ; }
    Equation {
      Galerkin { [ invrho[] *  Dof {d Volt} ,  {d Volt} ] ;  
                 In Bar ; Jacobian JVol ; Integration GradGrad ; }

When I run the Dirichlet at both ends I receive a Warning: Null right Hand Side in Linear Equation
which I assume is correct given the problem formulations.  I also get an indication that NaN is
being generated in the SOlver (GMRES).

Not surprisingly the Warning goes away when I use the Neumann BC.  THe Solver does not report any
NaN problems. <<simpleBar.geo>>  <<simpleBar.pre>>  <<solver.par>> 

Thanks,  the mail archives have been a big help

Laird Bolt

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060228/d11f2bd6/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: simpleBar.geo
Type: application/octet-stream
Size: 515 bytes
Desc: simpleBar.geo
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060228/d11f2bd6/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: simpleBar.pre
Type: application/octet-stream
Size: 10334 bytes
Desc: simpleBar.pre
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060228/d11f2bd6/attachment.pre>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: solver.par
Type: application/octet-stream
Size: 5031 bytes
Desc: solver.par
URL: <http://www.geuz.org/pipermail/getdp/attachments/20060228/d11f2bd6/attachment.par>

More information about the getdp mailing list