Hi all,

is it possible to block the velocity components u,v for a predefined cell id so that only w exists??

thx

Check the usersubroutine "sormom.f" for more info.

the usersubroutine sormom.f can't be used for velocity blocking!!!!

That depends on what you mean with blocking. I assumed you wanted 0-velocity for some cells.

Perhaps you could turn that cell in question to porous and define the coeficients such that there is no resistance in the w direction and infinite resistance in the other 2 directions?

Yes ,thank you for your help. I have adjusted the values so that nearly no resitance appears. I've tried the some with sormon.f in combination with sorent.f but this leads to strange results. So I think it is not possible to define a flow-equalizer(u,v blocked, w no restriction) with the subroutines.

I do not know what is your problem about but you can set arbitrary velocity at any cell using the sormom.f subroutine just write

if (ictid.eq.prostar_cell_type_number_where_the_subro utine_should_apply) then

s1u=velocity_you_like*great

s2u=great

and similarly for the rest

endif

instead of ictid you can use ip or something like that for cell number

using this procedure may appear unstable

regards

Why even bother with user subroutines? Can't you just use prostar to define which cells are porous and adjust the porous material properties accordingly?