Hi Francesca,

On 12/10/2015 11:41 AM, Francesca Camorrino wrote:

Hello! I'm trying to solve Poisson equation with Sfepy in a very simple domain, but I found some problem:

1. I need triangular mesh on a rectangular/square domain

For this, you can try gmsh. Why do you need triangles, by the way?

2. I have null von neumann conditions ( and I omitted them in boundary condition) like in the discussion https://groups.google.com/forum/#!searchin/sfepy-devel/null$20von$20neumann/...

OK

3)I know that this problem have many convergence problem, in one of this case is the simulation successful? Or I find some error?

You fix the potential in some, so it should work fine.

But your definition for region 'p' selects nothing for the mesh you generate - try

 'p' : ('vertices in (x <= 0.99951 )', 'vertex'),

instead.

Always verify your regions using --save-regions-as-groups option:

./simple.py prob_desc_file.py --save-regions-as-groups

• then check the 0_regions.vtk file using postproc.py.

I don't know because the simulation is successful but in the postprocessing the result isn't correct I attach the problem description file too. If you want I would like some advice about my problem

When I run your code, it does not do anything, because there is no load - maybe some right-hand side term is missing? Or provide a Dirichlet BC also in the region 'n' (again, fix the region definition).

r. PS: I missed your message, so my reply is late.