[angel_97]

Hello,

I am dealing with lagrangian particle and I have to insert a new force inside the model.
I went in src\lagrangian\parcel\submodels\Momentum\ParticleF orces and I found a lot of force but I do not understand the line of the code.

For example this drag force:

template<class CloudType>
Foam::WenYuDragForce<CloudType>::WenYuDragForce
(
CloudType& owner,
const fvMesh& mesh,
const dictionary& dict
)
:
DenseDragForce<CloudType>(owner, mesh, dict, typeName)
{}


template<class CloudType>
Foam::WenYuDragForce<CloudType>::WenYuDragForce
(
CloudType& owner,
const fvMesh& mesh,
const dictionary& dict,
const word& typeName
)
:
DenseDragForce<CloudType>(owner, mesh, dict, typeName)
{}


template<class CloudType>
Foam::WenYuDragForce<CloudType>::WenYuDragForce
(
const WenYuDragForce<CloudType>& df
)
:
DenseDragForce<CloudType>(df)
{}


// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //

template<class CloudType>
Foam::WenYuDragForce<CloudType>::~WenYuDragForce()
{}


// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //

template<class CloudType>
Foam::forceSuSp Foam::WenYuDragForce<CloudType>::calcCoupled
(
const typename CloudType:arcelType& p,
const typename CloudType:arcelType::trackingData& td,
const scalar dt,
const scalar mass,
const scalar Re,
const scalar muc
) const
{
const scalar alphac =
this->alphacInterp().interpolate
(
p.coordinates(),
p.currentTetIndices()
);
const scalar CdRe = SchillerNaumannDragForce<CloudType>::CdRe(alphac*R e);

return forceSuSp
(
Zero,
0.75*(mass/p.rho())*CdRe*muc*pow(alphac, -2.65)/(alphac*sqr(p.d()))
);
}

If someone can explain me line by line i such way I will be able to create a new force.

Thank you,

Best,

Angel