[Gmsh] Problem with free charge and floating potential

[krzysztof.pelczar at gmail.com]

  Dear All,

my gratitude for the software is beyond imagination. Nevertheless let
me first thank you for your wonderful job.
  As in the subject, I have a problem in formulation of electrostatic
problem with floating potential. I tried to follow the examples from
GetDP manual but with no visible result. The problem is as follows.
Using formulation without taking into account floating potentials and
free charges the solution I get for scalar potential and E-field seems
(still I can not judge if it correct as it is unexpected to work at
all...) to work (resolution EleSta_v in the project files attached).
For resolution with definition required by floating potentials, based
on the manual, namely EleSta_vf, I get results file full of zeroes,
and postoperation fails with memory segmentation. I believe there is
my mistake in formulation hidden, but I can not figure out what
exactly is wrong and how to get charge distribution properly.
  I am kindly waiting for your reply.

PS: the mesh was done using twice refinement by splitting with default
settings. Windows, Gmsh version 2.5.0, GetDP 2.1.0.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: float.zip
Type: application/zip
Size: 47898 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20110329/8791fc1f/attachment.zip>