Pressure Correction
 

Hello,

I am modelling a cylinder oscillating in a uniform current due to VIV. The simulation will run up to about 25 seconds flow time, but then it fails due to a pressure correction.
I have had a read through all the other threads available about trying to solve this problem. I am running the simulation with PISO pressure-velocity coupling with skewness-Neighbour Coupling disabled. My URF values are 0.3,1,1 and 0.7 - Pressure, Density, Body Forces and Momentum respectively. Under the Discretization tab the pressure is set to Standard and Momentum as Second Order Upwind.
I have tried to change all these values but with still no success.

However.

When I change my time step from 0.001 s to 0.003 s the analysis will keep iterating. If I let it run at this time step I run into a negative cell volume in my mesh, but If I quickly change the time step back to 0.001 s after 1 or 2 iterations the analysis will run for about 10 more time steps then fail due to a negative cell volume.

Can anyone please advise me on how to resolve this problem?

Many thanks.

P.

Could it be anything to do with the number of iterations per time step? Currently I run the simulation at 1 iteration per time step, and plan to crank it up to 20-30 when the cylinder settles into a regular motion.

Cheers

Currently I am running Fluent in 2d, would running it in 2ddp be better?

Regards.

The double precision solver does help in getting improved results however I don't think its the solution to your problem. Check your mesh. Also try using the green gauss node based scheme for gradient. Maybe this is your problem.

Neville,

Thanks for your advise. I am running the simulations as we speak, they take out 24 hours to run so I will post how I get on.

Thanks!

P

Neville,

This may seem trivial. But when you say check the mesh, What am I checking for?

Peter

Sorry to flag his up again people.

Can anyone explain to me what I should be checking for in my mesh to try stop getting these pressure convergence errors?

Thanks for your help

Peter

Hello Peter,
Sorry for the late reply but By mesh check I mean you should see whether your mesh is uniform and that cells are not clustered in particular region unless necessary refinement is required in that region.

Also the min cell volume for the mesh should be positive. I had an experience the mesh being discontinuous in my work. this is due to the fact that the mashing algorithm failed in that region but I was lucky to notice it visually.

Neville,

Thanks very much for that info there.

One thing I have found with my problem seems pretty weird. If I reduce the mass of my cylinder, the simulations run for a lot longer. I ran a simulation a few months ago with an extremely light cylinder that was critically damped. And it seemed to settle into a steady state.

Any idea why this could be?

Thanks again,

Peter

Sorry Peter, but I can't help you on this as I have not worked on vibrations.