Adding porous zones to compressible multiphase solvers
 

Dear all,

I am new to OpenFoam. I wish to add porous zones support to compressibleInterDyMFoam. By checking code of (interFoam, porousInterFoam, interDyMFoam, compressibleInterFoam & compressibleInterDyMFoam), it seems that I just have to include both the porousZones.H and createPorousZones.H headers to add porous zones support to the compressibleInterFoam or compressibleInterDyMFoam. Is this enough? Or should there be any other issues that I should attend to?

Best regards,
Hisham El Safti

I have such a question, too. I wish to add a porous zone support to MRFSimpleFoam. Is it possible?

Have you check your calculation?

I have not done it yet, but I plan to. When I do, I'll post the results. Of course it would be better if I would know if more work is to be anticipated than just including headers!

Do you know how to combine two solvers together?

Hi Hisham,

Did you succeeded in combining the two solvers? If yes, then could you share your brief method in doing so.

Thanks,

Manoj

Hello Manjo

I did not have time to do that just yet.

But I'm quite sure, now, that it is that straight forward (I think :confused:)

I may try to do it soon. I just got involved in other stuff

I know this is an old post but answering may help others.

You should also add the corresponding term to the momentum equation (UEqn). "pZones.ddt(rho, U)" on the left side.

This works on versions up to 2.1.
I think that 2.2 has a different way of treating porousZones.

hi marcelo, i tried with your suggestion but without succes unfortunately! I've tried to add in the Ueq both the term you suggested and, in a second time the term "pZones.addResistance(UEqn);" after Urelax (as is done in interporousfoam and in rhoPorousetc).
Have you other suggestions?actually the problem seems to fall in the create fields. I receive an error :
--> FOAM FATAL ERROR:
LHS and RHS of + have different dimensions
dimensions : [1 -3 0 0 0 0 0] + [0 0 0 0 0 0 0]

when the line:
volScalarField rho1("rho1", eos1->rho(p, T));

is reached.
I know that the post is very old but.... :D

Obviously....if someone else had good experiences with porous zones via fvOptions in more recent versions of compressibleInterFoam or twoPhaseEulerFoam is welcome!

thnx in advance to evrbdy ;)