[qwerty13]

Hi all,


I am simulating the transport of a passive scalar with OpenFOAM 2.2.x. using a modified solver with adaptive mesh refinement just to check how well AMR works. There are two things I wanted to ask in case someone can help:

1. If the boundary condition for the scalar at the inlet is set to fixedValue everything works well. However, if it set to totalFlowRateAdvectiveDiffusive then the simulation always crashes right after the lowerRefineLevel is reached in the dynamicRefineFvMesh file. Essentially the passive scalar equation cannot be solved (No Iterations 1001) and the scalar values just explode. Does somebody perhaps have any clue regarding this matter?
2. During the simulation I track the mass flow rate coming inside the domain from the inlet patch with the following piece of code:

massFlowRate
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl runTime; //outputTime;
writeInterval 0.005;
log false;
valueOutput false;
source patch;
sourceName inlet;
operation sum;
fields
(
phi
);
}



Regardless what the boundary condition is for the scalar at the inlet

(fixedValue or totalFlowRateAdvectiveDiffusive) when the adaptive mesh
refinement occurs then the calculated mass flow rate at the inlet is 4

times less than it should be which obviously has to do with the mesh

refinement. Any ideas on how to circumvent this problem?


Thanks in advance for any hint or help.