[futerxiaosage]

Dear members of the forum,

I have modified UEqn in UEqn.H as follows:

volScalarField alpha11 = mesh.lookupObject<volScalarField>("alpha.water");
volScalarField nu_sum=nu1*alpha11+nu2*(1-alpha11);
fvVectorMatrix UEqn
(
fvm::ddt(U) + fvm::div(phi, U)
+ MRF.DDt(U)
- fvm::laplacian(nu_sum,U)
==
fvOptions(U)
);

And I also have modified the pEqn.H as follows:

surfaceScalarField phig
(
(
// mixture.surfaceTensionForce()
- ghf*fvc::snGrad(rho)
)*fvc::interpolate(1/UEqn.A())*mesh.magSf()/fvc::interpolate(rho)
);

Ignoring the influence of surface tension, I have modified the fvm::laplacian in fvm.upper() by dividing the term rho[faci]. However, I have noticed significant differences in the pressure field results compared to the original interFoam solver after the first time step. I would like to understand where the main issue lies and how to make the necessary modifications. Thank you！！！

[futerxiaosage]

If anyone is interested, I would be glad to engage in a discussion and collaborate on further exploring this issue. It would be beneficial to exchange ideas and insights to find potential solutions and improve the modifications made to UEqn and pEqn.