[tbhanderi]

Hi,

I am new to OpenFOAM. I am trying to solve 2D steady state heat conduction problem as per case setup in URL link.

I have postprocessed the results in Paraview, got the correct temperature profile "but heat flux is not conserved between both the boundaries."

Anyone has an idea what could be the issue?

However in Ansys heat flux is conserved with the same mesh and boundary conditions.


Note :
Solver-laplacianFoam.
To calculate heat flux at the boundaries in Paraview , I am using extract block and integrate variable options.



casefile download link : https://drive.google.com/open?id=1Bf...BDSFMrtYEyPR07

[geth03]

hi,
i did face a similar problem, and i believe that the calculations with openfoam are correct but paraview postprocessing misses the point somewhat.

i have a pipe with two fluids, i compare the integrated velocities, which are the volume fluxes, and strangely i don't get the same values at the inlet and outlet of the pipe with paraview for those two fluids. when i postprocess w/ openfoam functions with correct weighting i get the correct results and mass is conserved.

i think paraview's integrated variables doesn't use correct weighting and creates wrong answers.

[tbhanderi]

Hi,

Thank for the reply.

I will try with OpenFOAM postprocessing functions. I also believe that the issue is with paraview integrating functions because i got all the other results correct.

As far as your case is concern (Pipe flow with two fluids), Paraview calculates correct fluxes as I have checked flux conservation in some pipe flow cases. Could you provide your case just to check flux conservation in paraview ?

[geth03]

in my case i'm using a multiphase solver, so there is the void fraction involved.
when i integrate the velocity, normally i should get the flux, but the void fraction weighting is not considered, so just by comparing those integrated velocities i don't get the correct volume flux. paraview doesn't take this into account.