|
[Sponsors] |
August 15, 2013, 11:51 |
The sum-up of forces on particles
|
#1 |
Member
X Meng
Join Date: Jun 2012
Location: Scotland
Posts: 89
Rep Power: 13 |
Hi, guys, I am using the solver 'icoUncoupledKinematicParcelFoam' to do some validations.
My question is, to a particle, or a parcel, the force on it should be from both environmental field (like drag force) and collision. Now, I can find the programme for drag force, gravity force, etc, I can also find the programme for collision spring force as well. But, still can't find where they are added together for the particle or parcel. Really confused... |
|
August 15, 2013, 14:26 |
|
#2 |
Senior Member
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 621
Rep Power: 0 |
Try looking at
Code:
$FOAM_SRC/lagrangian/intermediate/lnInclude/kinematicParcel.C Code:
const forceType& forces = td.cloud().forces(); // Momentum source due to particle forces const parcelType& p = static_cast<const parcelType&>(*this); const forceSuSp Fcp = forces.calcCoupled(p, dt, mass, Re, mu); const forceSuSp Fncp = forces.calcNonCoupled(p, dt, mass, Re, mu); const forceSuSp Feff = Fcp + Fncp; // New particle velocity //~~~~~~~~~~~~~~~~~~~~~~ Good luck. |
|
August 15, 2013, 16:04 |
|
#3 |
Member
X Meng
Join Date: Jun 2012
Location: Scotland
Posts: 89
Rep Power: 13 |
Thank you very much!
I found that and read it. By my understanding now, the velocity of parcel should be calculated twice, the one is in 'KinematicParcel.C' for particle forces, like drag force and gravity, while the other one is in 'CollidingParcel.C' for those forces due to collision. So, actually, the forces on the parcel are not collected altogether, but divided into two parts. Am I right? Thank you very much again! |
|
November 17, 2022, 12:48 |
|
#4 |
Member
Shah Akib Sarwar
Join Date: Mar 2021
Posts: 41
Rep Power: 5 |
Hello.. did you ever find an answer to this?
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
AMI interDyMFoam for mixer nu problem | danny123 | OpenFOAM Programming & Development | 8 | September 6, 2013 03:34 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 16:33 |
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 | bookie56 | OpenFOAM Installation | 8 | August 13, 2011 05:03 |
IcoFoam parallel woes | msrinath80 | OpenFOAM Running, Solving & CFD | 9 | July 22, 2007 03:58 |
Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 19:07 |