[renyun0511]

Hello,
I am studying pressure correction equation in SIMPLE algorithm by OpenFOAM-1.5-dev, and now I encounter a tough problem. So, I need your help.

The SIMPLE algorithm I studied is in the solver of MRFSimpleFoam and I want to add a new pressure prediction between the mass conserving pressure correction equation and the application of under-relaxation. Actually, the new pressure is the average pressure evaluated by multidimensional averaging. Thus
……………………………………………………………….….①
Where, is a previous iteration of pressure field. is computed via pressure equation in SIMPLE of MRFSimpleFOAM. is evaluated by multidimensional averaging. is the distance between the owner point and neighbor cell centers. is the neighbor nodes

Then, the new pressure field distribution () is used for correcting pressure and velocity by applying the under-relaxation, and this step is the same as it in SIMPLE .
My problem is how to calculate and represents them in OpenFOAM codes.

Thanks!