CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Friction

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   May 28, 2012, 15:14
Default Friction
  #1
New Member
 
Join Date: Mar 2012
Posts: 25
Rep Power: 5
callahance is on a distinguished road
Hello all,

i am trying to induce a friction function on one patch of my mesh (exp. lowerWall). I am trying to add a friction term in the UEqn. I thought i can make a surfaceVectorField that acts only on the lower patch of my mesh but im not so good with C++. Anyone would have an idea about how to implement such field ?

Thanks in advance.
callahance is offline   Reply With Quote

Old   May 29, 2012, 10:44
Default
  #2
New Member
 
Join Date: Mar 2012
Posts: 25
Rep Power: 5
callahance is on a distinguished road
Nobody has a clue ?
callahance is offline   Reply With Quote

Old   June 2, 2012, 08:50
Default
  #3
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
Hi!

what type of function are you thinking of?

if you want to set a specific roughness then you have bcs for this. On this lowerwall if choose this rough bc and at the other you may set the standard bcs.

if you want to make a time dependent sort of bcs..then you can try the swak4Foam pack..

explain better what you want to do so we can try to help you
__________________
Best Regards
/calim

"Elune will grant us the strength"
calim_cfd is offline   Reply With Quote

Old   June 2, 2012, 10:48
Default
  #4
New Member
 
Join Date: Mar 2012
Posts: 25
Rep Power: 5
callahance is on a distinguished road
ok the problem is : a droplet sliding on an inclined plane. In nature, the droplet pinns till the inclination gets over a critical value (lets say 10 degrees). I tried in interFoam to implement different models of dynamic contact angles but no good results (the droplet starts sliding even at 1 degree inclination). I noticed that there is a friction force between the droplet and the surface which wont be considered by interFoam; F_c = f(theta_A , theta_R).

Now i want to add this force in the UEqn (Navier stokes equations). The problem is i dont know how to add a vector field that acts only on a certain patch (in this case the lowerWall).
Hope the problem is a bit clearer now
callahance is offline   Reply With Quote

Old   June 2, 2012, 12:44
Default
  #5
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
hey!


yeah it is.

if i were you i'd try using a roughwall bc to increase the reaction force of the surface.. what keeps the drop still is its surface tension and the normal component of the friction (am i missing sth?).. in this case the roughness of the wall should play an important role.. so i say double check this roughness issue..along with the value for the surface tension ur using... discretization associated with the surface tension value could trigger is phenomenom earlier i guess..large chunks of the water drop baddly refined....not sure thou..

i cant tell atm cuz im not on my desktop/of but have you gone though all the alphaContactAngle*.* files?

i guess you problem can be solved by implement at the worst case a new boundary condition instead of modifying a solver.. thw swak4foam could be an option.. idk.. im not an expert in multiphase flow

im kinda interested in your case now... if you could send ur mesh file with the 0 folder files it would be nice
__________________
Best Regards
/calim

"Elune will grant us the strength"
calim_cfd is offline   Reply With Quote

Old   June 2, 2012, 13:45
Default
  #6
New Member
 
Join Date: Mar 2012
Posts: 25
Rep Power: 5
callahance is on a distinguished road
Thanks for ur response ,

the only bcs i know for a wall are wall and slip. Do you know how i can increase the roughness ? by the way i dont know if the parameter "roughness" is implemented in the bc. I think modifying the solver by adding a force is easier than making a new bc. I went through the alpha..angle. but didnt find anything (i think im missing something in it; interesting is with a kistler dynamic contact angle model, the droplet slides a shorter path than with Seeberg dynamic contact angle model.) I attended a Symposium for Openfoam last weekend in Germany and such problems were discussed. The problem ist that on the surface of a droplet, there are two surface tensions of opposite directions. And now when u reach the wall, there is only one surface tension along the interface, while the opposite one wont be considered cz of the mesh (this is just a numerical mistake). This problem was solved by adding an artificial surface tension on the contact line (physically incorrect but the results are beautiful... ).
Anyway i dont have the case right now on this laptop. Ill upload it when i can. Ill keep u posted of any news or improvements
callahance is offline   Reply With Quote

Old   June 2, 2012, 15:18
Default
  #7
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
hihi!

slip and wall? check the folder with the bcs.. i dont recall the path right now but it's somewhere in /src/finitevolume/*fvpatchfields*/derived/ .. search for "uniformfixedvalue" bc inside src/finitevolume and you find the folder with all derived bcs.. you should find the ones for multiphase flow there as well.. you can set different alpha contact angle per patch or make it time dependent.. you should find the rough wall bcs there too.. search for rough inside the same folder.. when you use this bc you'll have to set an additional paramenter which is related to the average height/depth of the surface imperfections.. hence friction....you do know how to search for parameter and their definitions in openfoam right? cuz the user guide and the programmer guide tell you a lot but sometimes you need to dig deeper hehe

l8r and gl again
__________________
Best Regards
/calim

"Elune will grant us the strength"
calim_cfd is offline   Reply With Quote

Old   June 2, 2012, 17:39
Default
  #8
New Member
 
Join Date: Mar 2012
Posts: 25
Rep Power: 5
callahance is on a distinguished road
i couldnt find rough bc. in /src/finiteVolume/fields/fvPatchFields/derived/

im usin OF 2.0.

EDIT: ok i found the roughwallfunction in src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/ ... but does this work for interFoam ?
callahance is offline   Reply With Quote

Old   June 2, 2012, 20:08
Default
  #9
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
yeah.. all bcs work with all solvers in of.. as long as they are consistent physically.. the one you pointed is meant for turbulent flows.. and your case is laminar.. but i guess you can use one of those.. as i told you i dont have of atm... tomorrow night i might point you sth more specific and point you the bc i think better suit ur case.. if you odnt find it earlier
__________________
Best Regards
/calim

"Elune will grant us the strength"
calim_cfd is offline   Reply With Quote

Old   July 25, 2013, 06:46
Default
  #10
New Member
 
Join Date: Apr 2013
Posts: 24
Rep Power: 4
mebinitap is on a distinguished road
Hi,
I am trying to model a droplet on an oscillating surface. But the drop seems to slip out with time even for a fixed wall (not oscillating). I used zeroGradient BC for alpha1 and p_rgh and fixedValue (0 0 0) for U at the wall (fixed). The rest of domain are atmosphere.
As in previous posts, is it related to some kind of surface roughness, which is to be implemented to prevent the slipping in this case also. If so can anybody please tell me how to do that..

Regards
mebinitap is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Skin Friction and Shockwave evaluate problem mosessitio FLUENT 0 May 21, 2011 03:21
Skin friction in fluent is wrong syler3321 FLUENT 2 March 24, 2011 09:05
Openfoam friction coefficient differences. Peter88 OpenFOAM Post-Processing 0 February 8, 2011 11:07
Compute Skin friction coef by hand Francois FLUENT 1 February 10, 2006 07:21
caculation of friction factor Dilip Main CFD Forum 8 May 5, 2001 00:42


All times are GMT -4. The time now is 16:21.