[sturgeon]

I was hoping someone could clear something up for me. I tried searching the forums but not finding out much.

As I understood it, fvc::flux converts a volume field to a surface scalar field. In the case of things like fvc::flux(HbyA), I understand its use in the algorithm. I also thought fvc:reconstruct essentially reverses the process.

However, I'm looking for errors in a solver I'm making and checking out the numbers produced and notice something strange. I have phiHbyA = fvc::flux(HbyA), and then I added some lines to output both HbyA and fvc::reconstruct(phiHbyA), and the values do not match up.

I'd understand small differences due to numerical error, but the difference between the two is fairly significant.

Am I misunderstanding the function of flux and reconstruct? Or are these errors indicative of some error in my code?