[sepp.zell]

Originally Posted by krndv

sorry , i am new to openfoam. Please tell me where to copy this (means which file and where in file)

thank you

Copy the source files of porousSimpleFoam to your case directory. The files can be found by using

cd $FOAM_SOLVERS/incompressible/simpleFoam

Copy the whole directory of porousSimpleFoam (as mentioned for icoFoam in the tutorial).
As mentioned in the bug-report you have to replace

U -= rUA*fvc::grad(p)

by

U = fvc::reconstruct(phi)

in the file pEqn.H. Also changes in

 .../Make/files

have to be done as mentioned in the icoFoam tutorial. Do not forget to give your solver a new name. Finally run

wmake

to compile the new solver.

[krndv]

Originally Posted by sepp.zell

Copy the source files of porousSimpleFoam to your case directory. The files can be found by using

cd $FOAM_SOLVERS/incompressible/simpleFoam

Copy the whole directory of porousSimpleFoam (as mentioned for icoFoam in the tutorial).
As mentioned in the bug-report you have to replace

U -= rUA*fvc::grad(p)

by

U = fvc::reconstruct(phi)

in the file pEqn.H. Also changes in

 .../Make/files

have to be done as mentioned in the icoFoam tutorial. Do not forget to give your solver a new name. Finally run

wmake

to compile the new solver.

Thanks for your Reply.
I am unable to find this equation , Can you please tell in which file it belongs to?
peq.h or ueq.h ?

[sepp.zell]

Originally Posted by sepp.zell

As mentioned in the bug-report you have to replace

U -= rUA*fvc::grad(p)

by

U = fvc::reconstruct(phi)

in the file pEqn.H.

As I wrote you have to look into the file pEqn.H.

Since I do not know, which version of openFoam you are using, I can not give you all the details. In the newest version it should correspond to

if (pressureImplicitPorosity)
{
    U = HbyA - (trTU()&fvc::grad(p));
}
else
{
    U = HbyA - (trAU()*fvc::grad(p));
}

I guess here you have to change the lines.