[leroyv]

Dear OpenFOAM users,
I recently wrote a simple custom boundary condition which computes the fluid velocity at the wall using a formula of the type:

v = f(rho) * normal

While my code yields the expected results using the rhoPimpleFoam solver, everything is messed up when using it with chtMultiRegionFoam, even though I use the exact same parameters for the BC. The near-wall velocity has an inverted orientation, and the upstream flow has a stagnation point, even if the wall velocity is oriented downstream. This is obviously a physical and numerical nonsense ...

Is there something in particular that should be taken care of when using chtMultiRegionFoam with custom BCs?

Parameters with rhoPimpleFoam:

squareWall
    {
        type            pyrolysisVelocity;
        refValue        0.001;
        massPR          0.001;
        rhoS            -1;
        rhoName         rho;
    }

Result: https://www.dropbox.com/s/n7cfws80w3...Foam_outer.png

With chtMultiRegionFoam:

fluid_to_solid
            {
                type pyrolysisVelocity;
                refValue 0.001;
                massPR 0.001;
                rhoS -1;
                rhoName rho;
            }

Result: https://www.dropbox.com/s/8t2re4piih...Foam_outer.png