[dela]

Hello All,

Please I want to add the attached equations as a momentum source to my reactingTwoPhaseEulerFoam's Ueqn but I don't know how to write the code in OpenFOAM

Below is what I have done:

[volScalarField Kdg = 1.9*160*pow(as,2)*phase1.mu()/(36*0.95)+(1.37*0.16*as*alpha1*rho1*mag(U1)/(6*0.95)); //effective interfacial area per unit volume
volScalarField Kdl = 1.9*160*pow(as,2)*0.95*phase2.mu()/(36*pow(alpha2,2)+1)+(1.37*0.16*as*0.95*rho2*mag(U 2)/(6*alpha2+1));

volVectorField Udg = -f*mag(U1)*fvc::grad(alpha1)/(alpha1+1); //drift velocity, gas phase

volScalarField Kd(fluid.Kd());

volVectorField Udl = -f*mag(U2)*fvc::grad(alpha2)/(alpha2+1); // drift velocity, liquid phase

volVectorField F1 = alpha1*Kdg*Udg+0.95*Kd*(Udg-Udl); //gas dispersion force

volVectorField F2 = (alpha2+1)*Kdl*Udl;//-0.95*Kd*(Udg-Udl); //liquid dispersion force]