CFD Online Discussion Forums - The sum-up of forces on particles

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)

- - The sum-up of forces on particles (https://www.cfd-online.com/Forums/openfoam/122259-sum-up-forces-particles.html)

mxylondon

August 15, 2013 10:51

The sum-up of forces on particles

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...

chegdan

August 15, 2013 13:26

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.

mxylondon

August 15, 2013 15:04

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!

Shah Akib Sarwar

November 17, 2022 11:48

Hello.. did you ever find an answer to this?

All times are GMT -4. The time now is 13:28.