CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > Siemens > STAR-CCM+

Doubts about time step and iteration number

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 26, 2018, 07:58
Default Doubts about time step and iteration number
  #1
New Member
 
Carlos Cordeiro
Join Date: Mar 2018
Location: Denmark
Posts: 10
Rep Power: 8
ccordeiro is on a distinguished road
Good afternoon to all.
First let me be really clear: I'm really a beginner in CFD simulations and in fact never had any lectures on it. I'm writing my Master thesis and I need to perform some CFD simulations regarding the effect of wind in a suspended bridge. The main goal is to look into turbulence and wake eddies.

Details about the model:
- Main domain size = 400x140m
- Bridge deck domain size = 50x15m
- Bridge deck dimensions = 31x4m

- Base size for mesh in main domain = 3m
- Base size for mesh in bridge domain = 1m
- Volumetric control created for overset mech

- Solvers: implicit unsteady, turbulent, constant density, RANS, SST k-omega,

For this initial solution, body motion is defined as stationary. However further I'll need to simulate it with rotational motion.

I choose a time step of 0,1s with 15 iterations. After 24h running the plots looks like the ones attached, but the simulation is not completed yet.

My question is: with these conditions, there are some methods (rule of thumbs) to estimate if my values for time step and iterations make some sense? I'm afraid these long simulations can take too long and after the result will be trash.

Any feedback will be much appreciated.
Attached Images
File Type: jpg Mesh.jpg (186.3 KB, 54 views)
File Type: jpg Residuals.jpg (80.7 KB, 59 views)
File Type: jpg Solver iteration CPU time monitor plot.jpg (71.0 KB, 49 views)
ccordeiro is offline   Reply With Quote

Old   March 26, 2018, 08:42
Smile
  #2
Senior Member
 
ashokac7's Avatar
 
Ashok Chaudhari
Join Date: Aug 2016
Location: Pune, India
Posts: 260
Rep Power: 10
ashokac7 is on a distinguished road
Send a message via Skype™ to ashokac7
Quote:
Originally Posted by ccordeiro View Post
Good afternoon to all.
First let me be really clear: I'm really a beginner in CFD simulations and in fact never had any lectures on it. I'm writing my Master thesis and I need to perform some CFD simulations regarding the effect of wind in a suspended bridge. The main goal is to look into turbulence and wake eddies.

Details about the model:
- Main domain size = 400x140m
- Bridge deck domain size = 50x15m
- Bridge deck dimensions = 31x4m

- Base size for mesh in main domain = 3m
- Base size for mesh in bridge domain = 1m
- Volumetric control created for overset mech

- Solvers: implicit unsteady, turbulent, constant density, RANS, SST k-omega,

For this initial solution, body motion is defined as stationary. However further I'll need to simulate it with rotational motion.

I choose a time step of 0,1s with 15 iterations. After 24h running the plots looks like the ones attached, but the simulation is not completed yet.

My question is: with these conditions, there are some methods (rule of thumbs) to estimate if my values for time step and iterations make some sense? I'm afraid these long simulations can take too long and after the result will be trash.

Any feedback will be much appreciated.
Don't know why you are using polyhedral mesh. For Overset mesh behavior, CCM recommends trimmer, isn't it? Please do make sure. The mesh count will be lot less with trimmer. Again for complex flow behavior, polyhedral performs better than trimmer.

The domain size you are dealing with quite large. So Plot Courant number and try to keep it near 1. Time step size though depend on type of problem, I would recommend using a script for time step such as for some initial hundred odd steps it is small then gradually increase the time step. Like $Time<0.1?0.001: ($Time<0.02?0.01:0.1).

The residual for continuity or momentum should fall by order of 3 or 4 (10^-3 or 10^-4) during each time step. If it does you can increase your time step. Otherwise reduce time step or increase inner iterations. The residual should tend to converge at end of time step which looks better in snapshots you have provided. Hope this helps!!
ashokac7 is offline   Reply With Quote

Old   March 26, 2018, 08:58
Default
  #3
New Member
 
Carlos Cordeiro
Join Date: Mar 2018
Location: Denmark
Posts: 10
Rep Power: 8
ccordeiro is on a distinguished road
Thanks for the feedback ashokac7.

Most of the parameters I took from another similar thesis which my supervisor gave me with the remarks the guys made a really good job. I guess in this case we can assume the flow behaviour is complex. Maybe in the screenshots is not visible but I've modelled also the railings since they're relatively close to the deck and the flow has a lot of disturbance there.

Can you formulate a bit more in detail how I can add a script for time step? Or point me out in the right direction.

Seems to me the iterations are quite ok. Especially after 5-6 iterations the values seem to converge. Do you agree or should I try to increase it?

Again, thanks a lot for the valuable inputs.
ccordeiro is offline   Reply With Quote

Old   March 26, 2018, 09:59
Smile
  #4
Senior Member
 
ashokac7's Avatar
 
Ashok Chaudhari
Join Date: Aug 2016
Location: Pune, India
Posts: 260
Rep Power: 10
ashokac7 is on a distinguished road
Send a message via Skype™ to ashokac7
Quote:
Originally Posted by ccordeiro View Post
Thanks for the feedback ashokac7.

Most of the parameters I took from another similar thesis which my supervisor gave me with the remarks the guys made a really good job. I guess in this case we can assume the flow behavior is complex. Maybe in the screenshots is not visible but I've modeled also the railings since they're relatively close to the deck and the flow has a lot of disturbance there.

Can you formulate a bit more in detail how I can add a script for time step? Or point me out in the right direction.

Seems to me the iterations are quite ok. Especially after 5-6 iterations the values seem to converge. Do you agree or should I try to increase it?

Again, thanks a lot for the valuable inputs.
For the script of time step

$Time<0.1 ? 0.001 : 0.01

This means for time less than 0.1 sec, time step is 0.001 sec (100 time steps with 0.001 s time step) otherwise 0.01 sec (After that it is 0.01 sec). You can add bracket as I have mentioned in previous comment for more gradual increase.

Residuals looks okay to me. Some more details about the problem is needed for further comment though.
ashokac7 is offline   Reply With Quote

Reply


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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
AMI speed performance danny123 OpenFOAM 21 October 24, 2020 04:13
Inconsistencies in reading .dat file during run time in new injection model Scram_1 OpenFOAM 0 March 23, 2018 22:29
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 05:28
Help for the small implementation in turbulence model shipman OpenFOAM Programming & Development 25 March 19, 2014 10:08
AMI interDyMFoam for mixer nu problem danny123 OpenFOAM Programming & Development 8 September 6, 2013 02:34


All times are GMT -4. The time now is 07:55.