i would like to use a Coulomb-Type BC for the simulation of a collapse of a granular column with interFoam.
Therefore, i want to mimic a kind of stick-slip BC by utilizing the partSlip BC and introducing a calculated valueFraction.
The flow at the boundary should feel an additional friction term F_c=-u_t/|u_t|*tan(delta)*N, where delta is the bottom friction angle, N the normal force acting on the Boundary and u_t is the velocity at the patch.
The change in velocity due to this force could be callculated, by assuming free slip, and correcting the velocity by:
u_t=u_fs - deltaT*u_fs/|u_fs| *tan(delta)*pressure*area
so i would like to calculte the valueFraction by:
1) Is it legal to handle the bottom friction in this way? Or is there a better way to introduce an bottom frition law?
2) I'm new to OpenFoam and, therefore, have some problems with the fundamental basics of implementing the BC.
a) Can I use
b) And to calculate the valueFraction could I use:
Thanks a lot!
The code is now partly working, i have now used:
But I will try to check if the BC is somehow meaningfull :rolleyes:
I am also interested in this coulomb-type slip bc. I guess you started from the partialSlipfvPatchField.C/.H and modified the "evaluate" member function? are there any caveats with this code, i.e. is it "meaningful", as you decribe in your last post?
thanks for your answer
|All times are GMT -4. The time now is 10:07.|