[bbita]

Dear Foamers,

I added a source term to pressure equation at the interface in interFoam as follows:
fvScalarMatrix p_rghEqn
(
fvm::laplacian(rAUf, p_rgh) == fvc::div(phiHbyA)+pTerm
);
pTerm is attached. I observed the attached pressure distribution. Both pTerm and pressure are at the first timestep.
Pressure and velocity boundary conditions are attached.
I am wondering if this behaviour makes sense or I should observe the pressure change only around interface.

Many thanks,