#78 Immersed boundary and rhoPimpleFoam: immersedBoundary::gradientInternalCoeffs() not implemented

[Daniel Cardoso]

Hi!

I've been attempting to use the immersed boundary wall condition for a compressible case using rhoPimpleFoam. However, I get an error message:
"""
--> FOAM FATAL ERROR:
Not implemented

From function immersedBoundary::gradientInternalCoeffs()
in file lnInclude/fvPatchField.H at line 473.

FOAM aborting
"""

I believe it is caused by using the immersed boundary wall in the Temperature field. It seems to be some problem in class inheritance between IB and thermal patch field...
I also happened to find a solution that didn't work for me in a Chinese forum that mentioned that Mr. Jasak was already informed about the problem - I don't speak any Mandarin though :(. https://zhuanlan.zhihu.com/p/389383026

Is there any solution for this problem? What could I do to workaround it?
Any information would be really helpful. Thank you for your time!

Thanks,

Daniel Cardoso

[Hrvoje Jasak]

Hi, The code fails in the thermo package because h boundary condition on the immersed boundary patch is set to a default patch field type (immersed boundary) that should not be used. There are options on how to fix this, ie either specify the h field or change the hBoundaryTypes to support immersed boundary.
This is broken in the entire OpenFOAM, for all coupled implicit or special boundary types which are not spefically listed in functions such as hBoundaryTypes or heBoundaryTypes. A good virtual function needs to be added and the entire hBoundaryTypes system decentralised.

This is a 2-week job (for me) and I cannot do it without a support contract or similar.