|
[Sponsors] |
April 5, 2011, 08:59 |
solidParticleFoam reflection?
|
#1 |
Member
Andrew Ryan
Join Date: Mar 2009
Posts: 47
Rep Power: 17 |
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; Code:
vector n = wpp.faceAreas()[wpp.whichFace(face())]; n /= mag(n); U_ -= 2 * (U_ & n) * n; |
|
April 18, 2011, 05:02 |
|
#2 |
Member
Sylvain Aguinaga
Join Date: Feb 2010
Posts: 41
Rep Power: 16 |
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 | Search this Thread |
Display Modes | |
|
|
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 |