# Initial drag for cylinder in cross-flow

 Register Blogs Members List Search Today's Posts Mark Forums Read

 May 28, 2012, 15:26 Initial drag for cylinder in cross-flow #1 New Member   Orxan Shibliyev Join Date: Aug 2011 Posts: 18 Rep Power: 8 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: Unsteady Fully Implicit, 1st order Euler temporal discre. 2nd order upwind spatial discre. Max. Inner Iterations = 10 2D, laminar Initially velocity and pressure field is zero Any idea?

 May 28, 2012, 19:46 #2 Senior Member   Chris DeGroot Join Date: Nov 2011 Location: Canada Posts: 388 Rep Power: 9 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.

 May 29, 2012, 05:44 #3 New Member   Orxan Shibliyev Join Date: Aug 2011 Posts: 18 Rep Power: 8 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.

 May 29, 2012, 10:06 #4 Senior Member   Chris DeGroot Join Date: Nov 2011 Location: Canada Posts: 388 Rep Power: 9 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.

 June 3, 2012, 17:43 #5 New Member   Orxan Shibliyev Join Date: Aug 2011 Posts: 18 Rep Power: 8 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.

June 4, 2012, 01:10
#6
Senior Member

Chris DeGroot
Join Date: Nov 2011
Posts: 388
Rep Power: 9
Quote:
 Originally Posted by orxan.shibli 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.
Very interesting. Thanks for sharing.

 Tags cross, cylinder, drag, flow

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Lord Kelvin OpenFOAM Running, Solving & CFD 8 March 28, 2016 11:08 danvica OpenFOAM Running, Solving & CFD 16 December 22, 2012 03:09 velan OpenFOAM Running, Solving & CFD 1 October 17, 2008 05:36 sivakumar OpenFOAM Pre-Processing 9 September 9, 2008 12:53 liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 21:28.