[pappo1890]

Hi Foamers.

The momentum source terms due to the dispersed phase can be identified as those due to the mass transfer and the forces acting on the particle. In literature, this term is expressed as in the figure in the attachment (Experiments and Numerical Simulations of Diluted Spray Combustion, B. Merci, D. Roekaerts, A. Saidiki., p.17)

I would like to have feedback about this expression and the one is used in OpenFOAM.

In OpenFOAM the momentum source term is written in the UEqn.H as parcels.SU(U).
SU(volVectorField& u) const is defined in line 563 of ~/OpenFOAM2.4.0/src/lagrangianba/intermediate/clouds/Templates/KinematicCloud/KinematicCloudI.H:

tmp<fvVectorMatrix> tfvm(new fvVectorMatrix(U, dimForce));
fvVectorMatrix& fvm = tfvm();

fvm.source() = -UTrans()/(this->db().time().deltaT());

return tfvm;

UTrans() is updated twice when solving the lagrangian libraries. First in the KinematicCloud, for the momentum spray source term due to the force acting on the particle and later in parcels/Templates/ReactingParcel/ReactingParcel.C (line 689):

td.cloud().UTrans()[cellI] += dm*U0;

U0 is (line 579):

const vector& U0 = this->U_;

and at line 731 the this-> U_ pointer is calculated as the particle velocity.

The expression does not completely look like that found in literature, where the difference between the carrier and the particle velocity is considered. Could, please someone double-check if my reading of the following lines of code is correct?

Thanks in advance