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

solidParticleFoam reflection?

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

Like Tree3Likes
  • 1 Post By andrewryan
  • 2 Post By Sylvain

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 5, 2011, 08:59
Default solidParticleFoam reflection?
  #1
Member
 
Andrew Ryan
Join Date: Mar 2009
Posts: 47
Rep Power: 17
andrewryan is on a distinguished road
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.
randolph likes this.
andrewryan is offline   Reply With Quote

Old   April 18, 2011, 05:02
Default
  #2
Member
 
Sylvain Aguinaga
Join Date: Feb 2010
Posts: 41
Rep Power: 16
Sylvain is on a distinguished road
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
randolph and rezaeimahdi like this.
Sylvain is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
shock tube reflection boundary conditions euler 1d micheletuttafesta Main CFD Forum 1 May 8, 2010 12:01
boundary reflection shuo OpenFOAM 1 March 9, 2010 16:10
boundary reflection shuo Main CFD Forum 0 March 9, 2010 02:18
solidParticleFoam particle positions caw OpenFOAM Running, Solving & CFD 1 October 6, 2009 04:50
Reflection coefficient of particle in DPM model S.J.R FLUENT 1 June 9, 2007 20:03


All times are GMT -4. The time now is 08:03.