|
[Sponsors] |
![]() |
![]() |
#1 |
Guest
Posts: n/a
|
Hi everybody,
I want to set a fixed velocity in some fluid cells in the middle of my mesh. I have tried to do this in "the subroutine "sormom.f" without success. Does anyone know how to do this? Thanks, Lennart |
|
![]() |
![]() |
![]() |
#2 |
Guest
Posts: n/a
|
Lennart,
Hi and greetings again! sormom is the exact right way to do this: C------------------------------------------------------------------------- C This subroutine enables the user to specify the momentum source C term (per unit volume) in linearised form: C C Source in x direction = S1U-S2U*U, (N/m3) C Source in y direction = S1V-S2V*V, (N/m3) C Source in z direction = S1W-S2W*W, (N/m3) C C in an arbitray manner. C C ** Parameters to be returned to STAR: S1U,S2U,S1V,S2V,S1W,S2W, C POROS C C------------------------------------------------------------------------- What may not be obvious is what to set s1u to etc., this is best shown in sorent.f: C----------------------------------------------------------------------- C C This subroutine enables the user to specify a source term (per unit C volume) for enthalpy in linearized form: C C Source = S1P-S2P*T, (W/m3) C C in an arbitrary manner. C C If temperature is to be fixed to a given value T, then the C following may be used: C C S1P=GREAT*T C S2P=GREAT, C C where T can be a constant or an arbitrary function of the C parameters in the parameter list. C C ** Parameters to be returned to STAR: S1P,S2P C C------------------------------------------------------------------------- so combining the two means that if you set SU1 to great*desired velocity and SU2 to great it fixes the velocity. You need to set all three velocity components. Generally, this is done in and IF block based upon cell table (ICTID). This is the info that would be supplied very quickly by an email to support@cd.co.uk regards Philip |
|
![]() |
![]() |
![]() |
#3 |
Guest
Posts: n/a
|
Thanks alot Philip. It works perfectly.
Regards, Lennart |
|
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to modify interFoam to keep a constant velocity field | thibault_pringuey | OpenFOAM | 0 | January 9, 2009 05:59 |
velocity-decay constant | hungryfish | CFX | 1 | December 18, 2008 16:27 |
wall boundary condition with constant velocity | Rogerio Fernandes Brito | CFX | 0 | August 17, 2008 10:01 |
transient temperature field with constant velocity | Törnquist | CFX | 0 | September 16, 2003 04:22 |
Specifying constant-velocity fans | Craig B | Phoenics | 3 | October 9, 2001 05:04 |