Trouble modelling fluid flow (air) between two roofs of a biogas plant
Posted Sep 19, 2019, 4:41 AM EDT Computational Fluid Dynamics (CFD), Heat Transfer & Phase Change Version 5.4 3 Replies
Please login with a confirmed email address before reporting spam
I am currently modelling heat loss in a biogas plant (manure digester) placed on a farm. The digester is cyllindrical in shape, with a double domed roof. Between the two roofs, air is pumped through to assure an outflow of potentially dangerous gases leaking through the lower roof, where the inflowing air has the outside temperature, much lower than the temperatures inside the digester (10 degC for the inflow, 42 degC for the manure). I included a snapshot of the geometry under this post, as well as the model's MPH file. I use a combination of the CFD module (turbulent flow) and the possibilities for Heat Flow (Heat transfer in solids and fluids) in the basic licence (I do not have the Heat Flux module to my disposal). The following are the problems I encounter:
- I am having trouble modelling (or possibly, understanding the result of) the air flow between the roofs. As can be seen in the snapshots included (one in the plane of the inlet and outlet and one perpendicular to that), the airflow directly at the inlet and outlet is correct, but somehow the flow between the two roofs is in the direction of the inlet (negative x-components of velocity field), rather than the outlet, which seems odd.
- Also, it seems like the airflow from the inlet somehow penetrates the inner roof, which is a closed surface of a PVC foil, which by no means is permeable to such an extent.
- Thirdly, it seems like the velocity of the airflow between the roofs has the highest magnitude directly under the upper roof, where due to the no-slip condition imposed, I would expect the highest magnitude in the middle of the distance between the two roofs.
- On top of that, the result does not display a correct circulation of the biogas below the inner roof, or at least one would expect. The result shows a cold 'bubble' forming behind the lower roof near the inlet. This can be seen in the Temperature snapshot.
- Lastly, when integrating the Velocity magnitude over the surface of the inlet and doing the same for the outlet, shows a different value for the two. This should not be correct, because there is no other way for the inflowing air to exit the roof than the outlet on the other side.
Can anyone elaborate on what I might be doing incorrectly, or if the model is correct, where my thinking and expectations are off?
Any help would be greatly appreciated. Thank you in advance!