On 12/09/2015 03:24 PM, Jacky Ko wrote:

The example http://sfepy.org/doc-devel/examples/navier_stokes/navier_stokes.html clearly illustrated a NS problem with inlet velocity normal to y axis.

Is it possible to set the inlet with arbitrary plane that is not normal to any axis? e.g. inlet normal velocity = (u1, u2, u3)

I see that in this example the inlet and outlet boundaries are defined by the cinc functions:

cinc_name = 'cinc_' + op.splitext(op.basename(filename_mesh))[0]cinc = getattr(utils, cinc_name) functions = { 'cinc0' : (lambda coors, domain=None: cinc(coors, 0),), 'cinc1' : (lambda coors, domain=None: cinc(coors, 1),),}

how would this utility work with my arbitrary inlet/outlet boundaries?

Concerning the question about the "region by function" definition, the docs are here [1]. You need to select vertices of the mesh that belong to the region, based on their coordinates. So either you can express a condition for inlet/outlet point coordinates mathematically, or define the regions using some other way. For example, you can assign an id number to each node and select a region based on that id ('vertices of group <integer>').

Best regards, r.

[1] http://sfepy.org/doc-devel/users_guide.html#functions