velocity blocking
 

Hi all,

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

thx

Re: velocity blocking
 

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

Re: velocity blocking
 

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

Re: velocity blocking
 

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

Re: velocity blocking
 

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?

Re: velocity blocking
 

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.

Re: velocity blocking
 

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

Re: velocity blocking
 

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?