[Tempest]

Dear Foamers,

I have been trying to simulate flow across a porous square heated cylinder, subjected to buoyancy. I am using BuoyantBoussinesqPimpleFoam with explicitPorositySource via fvOptions. I have attached the case file for reference. (BlockMesh + TopoSet + SetFields before trying simulation)

Now, I am getting negative lift variation, as per the literature, for solid cylinder (see solid.png in attachment). When I introduce porous cylinder instead of solid, the drag variation remains same, but the lift coefficient values (see porous.png in attachment) jump near zero, which is way off the literture value.

The case that I have attached is at Ri = 0.5 and Re = 40 and the coefficient of lift (average) should be near -1.

I have also attached the lift coefficient variation.

What can be the reason for this variation in Cl values and how can I solve it? I have tried almost all variations in putting the Boussinesq approximation into the code, but nothing seems to work out.

Please help me out. Thanks in advance.

P.S. The axis for Cl curve is on RHS.

[damu1414]

Dear Anirudh,

I have issues while simulating flow over a heated cylinder. Since you have successfully simulated the solid cylinder case, I would be grateful to have your suggestions.

buoyantBoussinesqPimpleFoam as such is not giving me negative lift coefficients. Also, the wake is not getting captured properly. Did you have any such issues? Please let me know.

Thank you