[lovecraft22]

Hi all;
I'm getting some strange results for a 2D circular cylinder run at 30 m/s, see this picture please:



Now, can anybody of you guys please check my case (attached) and tell me what I'm doing wrong?

Thank you!

[srakshit]

I'm getting some strange results for a 2D circular cylinder run at 30 m/s, see this picture please

In your setting,the velocity inlet is at 30 m/s but the outlet is at zeroGradient. This is creating the problem in your domain. keep the velocity constant at 30 m/s.

[MasterCooler]

Hey lovecraft22,

without having had the files opened and analysed (I am in a windows session) it reminds me of a problem I had some months ago. Then as a really new user I did the mistake to violate the rules to create blocks in blockMesh. My block was created instead of in a right-hand in a left-hand coordinate system. It was very strange that I could create the blockMesh without having problems. Also I could show it in paraView where it looked really nice and as assumed. As I didn't expect any problem I didn't execute checkMesh but started the calculation directly. And saw really strange results, inexplicable and senseless results like yours.

Therefore the question: have you performed a checkMesh run and says checkMesh that your mesh is okay?

Further I can't comprehend the hint of srakshit. I did several simulation with fixedValue for inlet velocity and zeroGradient for outlet velocity. As I know it should work if you use a zeroGradient condition for pressure at the inlet and fixedValue condition for pressure at the outlet.

Cheers
MasterCooler

[srakshit]

Yes, I agree with MasterCooler. Out of U and p, one value should be fixed for either inlet or outlet.

[lovecraft22]

Thanks to both of you.

I haven't performed a check mesh, I'll surely do that.
About my boundary conditions, how would you advise to set them as?

Thank you!

[Balakrshnan Ramakrishnan]

can you please explain why you have mentioned

frontAndBack as symmetryPlane ?

and upper and lowerWall as empty?

Are you trying to do a simulation in Z direction???

[lovecraft22]

I'm trying to simulate a flow in the x-y plane, so the planes parallel to this have to be empty…

[Balakrshnan Ramakrishnan]

are you sure??

can you try with frontAndBack with empty ,

upperWall and lowerWall with symmetry?

I think it will work..

[lovecraft22]

Actually, if you look at my blockMeshDict you'll see that frontAndBack are the planes parallel to the xz plane (red in the following picture) while upperWall and lowerWall are the planes parallel to the xy plane (yellow in the following picture).
So, I think this boundary conditions should be fine…

[srakshit]

If your domain is 1 cell thick in z direction, then keeping the frontAndBack patch as 'empty' will reduce the computational efforts (it wont solve U for that direction). Otherwise, i dont think it will matter. About the upper and lower wall, keeping them to 'symmetry' is a better choice. I have no idea what an empty patch will do at the wall, but again those area shouldn't effect the nature of flow you are working with.

Just use the suggestions mentioned earlier for your inlet and outlet boundary.

[lovecraft22]

Thank you Sukanata…
I haven't quite well understood what you mean with watch you suggested about the boundary conditions…

I'd really appreciate if you could give me some more detailed informations.

Thank you!

[srakshit]

ok...since you are doing flow past a cylinder with 30 m/s velocity. Your inlet boundary is consistent (since Velocity is initialized and pressure is kept floating).
Problem is with your outflow boundary. Both of your variables are floating (which may be the reason for waves bouncing inside the domain and velocity rising upto 900 m/s.) Initialize either pressure or velocity for the the outlet. Try using 30 m/s for velocity outlet if outlet pressure is unknown.

[lovecraft22]

Ok, this way thought I'll set the velocity at the outlet at 30 m/s but that's not the physical case… Am I correct?

[srakshit]

Are you validating some experiments?

If the pressure is unknown, you should define the velocity at the boundary. Keeping both of them floating is inconsistent. Since this is a steady state analysis of flow past a cylinder, you can assume a constant velocity.

[lovecraft22]

Ok, I'll set them as you suggested and we'll what comes out...

Just another quick question: which boundary conditions would you set if you were running a car in a free stream?

Thank you once again!

[joel.lehikoinen]

If you are running an incompressible simulation you can naturally fix pressure at the outlet.

[srakshit]

For an unsteady analysis, zero gradient for both pressure and velocity will be fine. But for steady state study, its better to have a steady state boundary for both inlet and outlet.

[lovecraft22]

Can't I set the pressure at the outlet to be the same as the freestream static pressure?

[srakshit]

Yes, you can try setting the outlet pressure at freeStream.

[lovecraft22]

Hi again;
I changed the boundary conditions but I'm still getting some really strange results, se the picture attached: the free stream should be coming from left to right at 30 m/s.

I've also attached my new boundary conditions.

Does anybody have a clue of what's going on here?


Thank you!