|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Join Date: Aug 2011
Posts: 30
Rep Power: 15 ![]() |
Hi everyone!
I am simulating flow over cylinder at Re=100 with my own code. Setting time step to dt=1s, I monitor initial values of drag coef. and it solves it without problem. But when I set it to say dt=0.001 then initially drag coef. rises to very high values (hundreds) and then it becomes negative even. It is not important whether it will be better as time passes since I am also trying to understand why decreasing time step causes this kind of problem. Other specifications:
Any idea? |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 414
Rep Power: 17 ![]() |
You are starting the flow from a standstill, so that takes a big pressure drop to accomplish. That is one reason for the high drag. Also, you are only doing 10 inner iterations, which probably means it is not converging within each time step. In that case, the initial results don't mean much anyways.
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Join Date: Aug 2011
Posts: 30
Rep Power: 15 ![]() |
Dear Chris
OK but this should be true for both time steps. That is why I decreased time step however it became worse unexpectedly. Also I tried to make inner iterations 100 but that only caused to reach high values faster. Last edited by orxan.shibli; June 3, 2012 at 17:45. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 414
Rep Power: 17 ![]() |
Well, your small timestep is 1000 times smaller than your large time step, so I would guess that the large time step just "steps over" a lot of the initial transient stuff.
Also, imagine in reality having completely still fluid and then instantly accelerating it to whatever value you have specified at the inlet. Not going to happen. It will take some (maybe small) amount of time to accelerate. In your simulation that time to accelerate is your first time step. Thus, I believe a smaller timestep will imply a larger initial pressure drop and also drag forces. In any case, these initial transients will die out eventually. |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Member
Join Date: Aug 2011
Posts: 30
Rep Power: 15 ![]() |
I found the problem. It was due to time-step dependency of original Rhie and Chow MIM. Hence, it is better to use modified version:
"Discussion on Momentum Interpolation Method for Collocated Grids of Incompressible Flow", Yu et al. |
|
![]() |
![]() |
![]() |
![]() |
#6 |
Senior Member
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 414
Rep Power: 17 ![]() |
Very interesting. Thanks for sharing.
|
|
![]() |
![]() |
![]() |
Tags |
cross, cylinder, drag, flow |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Extrusion with OpenFoam problem No. Iterations 0 | Lord Kelvin | OpenFOAM Running, Solving & CFD | 8 | March 28, 2016 11:08 |
Interfoam blows on parallel run | danvica | OpenFOAM Running, Solving & CFD | 16 | December 22, 2012 02:09 |
Negative value of k causing simulation to stop | velan | OpenFOAM Running, Solving & CFD | 1 | October 17, 2008 05:36 |
Unknown error | sivakumar | OpenFOAM Pre-Processing | 9 | September 9, 2008 12:53 |
Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 18:07 |