[nuttita]

Hi, I need to couple InterFoam with porousInterFoam for flows through deforming porous medium. InterFoam needs pressure field from porousInterFoam and in turn it will give alpha and U to porousInterFoam. The steps in the coupled model will be:

At time i, p —> solve InterFoam —> alpha, U
Time i+dt/2, alpha, U —> solve porousInterFoam —> p
Time i+dt, p —> solve InterFoam —> alpha, U
...

Does anyone know how to couple solvers in this way? Is there any utilities or libraries for this coupling? Please help. Thanks!

[Tobi]

Hi,
is there a reason why you cannot use porousInterFoam directly? You mentioned some deformation - I am not familiar with OpenFOAM porousInterFoam - however, if there is the mesh motion missing, I suggest to add this to the solver. In principle, coupling two solvers as you mentioned is not a good idea; however, there are different opinions.

[ebringley]

It sound to me like you're tyring to do some strang operator splitting type routine, but you didn't include your time steps for each step, so it's impossible to tell exactly what you're doing... To the best of my knowledge, there's no out-of-the box way to "just couple two solvers".

If porousInterFoam doesn't do what you want, you'll have to create your own copy of the source code and implement your own operator splitting, likely starting from porousInterFoam.