Anybody has advice? Thanks a lot.

Thanks a lot, Dongyue. Very helpful. I found in SN model CdRe was multiplied by another Re. Thus, the formula is consistent.

Hi, everyone. I am confused that the expression of drag force used in Openfoam, as shown in the following codes. Why there is nu and divided by square of particle diameter. I saw the formula is without nu and only divided by diameter in some literature. Do you have any tip for me? Thanks a lot. By the way, After I deleted pair_.continuous().nu() and only divided by diameter. The results I obtained are totally the same.

Foam::tmp<Foam::volScalarField> Foam::dragModel::Ki() const
{
    return
        0.75
       *CdRe()
       *swarmCorrection_->Cs()
       *pair_.continuous().rho()
       *pair_.continuous().nu()
       /sqr(pair_.dispersed().d());
}

I added a new scalar in twophasesystem/phaseModel.C and I compiled successfully. The problem is that there is an error when I run the case. My objective is to use this constant in phasePair.C to calculate Re(). In oder to change the constant myBeta directly from dictionary (phaseProperties) without compiling every time, I think it is necessary to add a new scalar like alphaMax in phaseModel.C. Do you have any tip?

Thanks a lot. The problem was solved.

@东岳 Hi, dongyue. I tried many ways to reach my goal, such as create a newtwophasesystem.C (I found this file in the reactingtwophaseeulerfoam, so it should be possible to use a new one) and compile. But I still cannot link the new file to the solver if I put it in the user's library. Do you have any suggestion? Many Thanks.

Anyone has suggestion about how to link the revised twophasesystem.C to the solver as a user's library?

@xpqiu Thanks for your reply. In fact, There is no error for the compilation and simulation. I don't know how to link the new modified twophasesystem.C to the solver.

@东岳 Hi Dongyue, it works. Now I am trying to add it into user's library. I changed the Make/files as following and I include "libmycompressibleTwoPhaseSystem" in system/controDict, but the solver neglect it. I really don't know how to make the solver to call the modified "twophasesystem.C", or I need to change the name of "twophasesystem.C". Do you have any suggestion? Thanks.

LIB = $(FOAM_USER_LIBBIN)/libmycompressibleTwoPhaseSystem

I see. Great!
Grazie mille.