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

The sum-up of forces on particles

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 15, 2013, 11:51
Default The sum-up of forces on particles
  #1
Member
 
X Meng
Join Date: Jun 2012
Location: Scotland
Posts: 89
Rep Power: 13
mxylondon is on a distinguished road
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...
mxylondon is offline   Reply With Quote

Old   August 15, 2013, 14:26
Default
  #2
Senior Member
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 621
Rep Power: 0
chegdan will become famous soon enoughchegdan will become famous soon enough
Try looking at

Code:
$FOAM_SRC/lagrangian/intermediate/lnInclude/kinematicParcel.C
and look at the method calcVelocity. In there, you should see where the forces are collected to calculate the new parcel velocity in the lines

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.
chegdan is offline   Reply With Quote

Old   August 15, 2013, 16:04
Default
  #3
Member
 
X Meng
Join Date: Jun 2012
Location: Scotland
Posts: 89
Rep Power: 13
mxylondon is on a distinguished road
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!
mxylondon is offline   Reply With Quote

Old   November 17, 2022, 12:48
Default
  #4
Member
 
Shah Akib Sarwar
Join Date: Mar 2021
Posts: 41
Rep Power: 5
Shah Akib Sarwar is on a distinguished road
Hello.. did you ever find an answer to this?
Shah Akib Sarwar 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
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


All times are GMT -4. The time now is 09:29.