[hiep.nguyentrong]

I look in to steadyCompressibleFoam in foam-extend-4.1 and run the simple tutorial test like transonicbump.

The base code used U instead of field HbyA run well.
When i create new volVectorField HbyA with same definition like below and replace U with HbyA, the code still run but divergence with wrong result.

What is the trick here? Can anyone explain it?

        volVectorField HbyA = rUA*UEqn().H();
        U =  rUA*UEqn().H();

        // Calculate phi for boundary conditions
        phi = rhof*(fvc::interpolate(HbyA) & mesh.Sf());