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...