solidParticleFoam reflection?

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 April 5, 2011, 07:59 solidParticleFoam reflection? #1 Member   Andrew Ryan Join Date: Mar 2009 Posts: 47 Rep Power: 8 Dear Foamers Looking at solidParticleFoam I didn't understand what the code in hitWallPatch (in solidParticle.C) is supposed to do, it looks like it just inverts the direction of the particle?? Code: ``` vector nw = wpp.faceAreas()[wpp.whichFace(face())]; nw /= mag(nw); scalar Un = U_ & nw; vector Ut = U_ - Un*nw; if (Un > 0) { U_ -= (1.0 + td.spc().e())*Un*nw; } U_ -= td.spc().mu()*Ut;``` Anyways I just replaced it a formula for specular reflection. Code: ``` vector n = wpp.faceAreas()[wpp.whichFace(face())]; n /= mag(n); U_ -= 2 * (U_ & n) * n;``` The result now looks much better in paraview.

 April 18, 2011, 04:02 #2 Member   Sylvain Aguinaga Join Date: Feb 2010 Posts: 33 Rep Power: 7 Hi Andrew, I think that the way the reflection is implemented in OpenFOAM takes into account: -the decrease of the normal velocity due to inelastic rebound, using the energy loss coefficient "e" -the decrease of the tangential velocity due to the friction of the particle with the wall, using the coefficient mu I think that if you take mu=0 (no friction) and e=1(perfect elastic rebound without energy loss) your formulation is the same as the one implemented in OpenFOAM. Looks like anyway… If you specify e=1 and mu=0, the results are very different between your implementation and the one used in OpenFOAM? Hope it will be helpfull… Sylvain

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post micheletuttafesta Main CFD Forum 1 May 8, 2010 11:01 shuo OpenFOAM 1 March 9, 2010 16:10 shuo Main CFD Forum 0 March 9, 2010 02:18 caw OpenFOAM Running, Solving & CFD 1 October 6, 2009 03:50 S.J.R FLUENT 1 June 9, 2007 19:03

All times are GMT -4. The time now is 05:22.

 Contact Us - CFD Online - Top