Hi Nicolas,

On 04/27/2016 04:49 PM, Nicolas DELAN wrote:

Hi everyone,

I am a new user of sfepy, i need it to compute the diffusions equations on a turbine blade. But i am a little lost with the boundary conditions. I have tried to impose some dirichlet boundary conditions on some points of the edges of the airfoil but it keeps crashing with exceptions. For example in my .mesh file i have this at the line where my point is defined

0.1 0.5 25 (x) (y)

So I have written this in the region definition

regions = { 'Omega' : 'all', 'point25' : ('vertices of group 25', 'edge'), }

So the group 25 vertices are forming an edge, right? Then this should work. If this is a single vertex, use 'vertex' kind instead of 'edge', as the 'edge' region would be empty.

I have tried before with :

refions = { 'Omega' : 'all', 'extrados' : ('vertices in (y>0.1) ', 'facet'), 'intrados' : ('vertices in (y<0)', 'facet'), }

and it was giving resuls (it was only for test purpose)

I don't understand very well how to define the boundary conditions (and to be honest the hole configuration file despite the documentation)

For example, to impose Dirichlet BC u = 1 in the point25 region to a scalar variable u, do:

ebcs = { 'a name' : ('point25', {'u.0' : 1.0}), }

For more examples, check [1], if you have not done that already.

Does that help?

r.

[1] http://sfepy.org/doc-devel/users_guide.html#essential-boundary-conditions-an...