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

The sum-up of forces on particles

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 15, 2013, 10:51
Default The sum-up of forces on particles
  #1
Member
 
Xiangyin Meng
Join Date: Jun 2012
Location: London
Posts: 69
Rep Power: 5
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, 13:26
Default
  #2
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 546
Rep Power: 18
chegdan 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.
__________________
Dan

Find me on twitter @dancombest and LinkedIn
chegdan is offline   Reply With Quote

Old   August 15, 2013, 15:04
Default
  #3
Member
 
Xiangyin Meng
Join Date: Jun 2012
Location: London
Posts: 69
Rep Power: 5
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

Reply

Thread Tools
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 On
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 02:34
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02: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 12:03.