CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

time step of twophaseeulerfoam

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

Like Tree1Likes
  • 1 Post By rdbisme

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 2, 2020, 12:02
Question time step of twophaseeulerfoam
  #1
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Hi guys,


I am using twophaseeulerfoam to simulate the solid-liquid pipe flow. The total number of the cells are 180000 (unstructure mesh) but I need to set a very very small time step to keep simulation normally runnning. Do you have any suggestion that I can define a larger time step? I saw some people used this solver to simulate the same phonomenon. They used more cells than mine. I cannot imagine how much time they used to finish one simulation if the situation like mine.

Thanks,
qi.yang@polimi.it is offline   Reply With Quote

Old   July 2, 2020, 14:40
Default
  #2
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by qi.yang@polimi.it View Post
Hi guys,


I am using twophaseeulerfoam to simulate the solid-liquid pipe flow. The total number of the cells are 180000 (unstructure mesh) but I need to set a very very small time step to keep simulation normally runnning. Do you have any suggestion that I can define a larger time step? I saw some people used this solver to simulate the same phonomenon. They used more cells than mine. I cannot imagine how much time they used to finishone simulation if the situation like mine.

Thanks,
Hello qi.yang, I think you might get more helpful comments if you give more details about your configuration. You can check my old post for an example of which details might be interesting.
rdbisme is offline   Reply With Quote

Old   July 2, 2020, 20:47
Default
  #3
Member
 
Ardalan
Join Date: Jul 2012
Location: Atlanta, USA
Posts: 77
Rep Power: 14
Ardali is on a distinguished road
Quote:
Originally Posted by qi.yang@polimi.it View Post
Hi guys,


I am using twophaseeulerfoam to simulate the solid-liquid pipe flow. The total number of the cells are 180000 (unstructure mesh) but I need to set a very very small time step to keep simulation normally runnning. Do you have any suggestion that I can define a larger time step? I saw some people used this solver to simulate the same phonomenon. They used more cells than mine. I cannot imagine how much time they used to finish one simulation if the situation like mine.

Thanks,



If you do not need energy equation discard it in the solver.

If you need it specify correct boundary conditions in T files.

diameterModel in phaseProperties matters as well.
If comprisibility matters to you a lot, e.g.sonic jet or so on, It will be reflected in p_rgh, therefore it wont be a sudden crash. Anyhow, small courant number in the order of 0.1 may help too.
Good luck
Ardalan
Ardali is offline   Reply With Quote

Old   July 3, 2020, 03:15
Default
  #4
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Thanks a lot. As you said, I already neglected energy euqation. Even though I set maxCo 0.1, the maximum time step is 2e-5 or 2e-6. The totall cells are 260000. How about your time step used in the simulation?

Quote:
Originally Posted by Ardali View Post
If you do not need energy equation discard it in the solver.

If you need it specify correct boundary conditions in T files.

diameterModel in phaseProperties matters as well.
If comprisibility matters to you a lot, e.g.sonic jet or so on, It will be reflected in p_rgh, therefore it wont be a sudden crash. Anyhow, small courant number in the order of 0.1 may help too.
Good luck
Ardalan
qi.yang@polimi.it is offline   Reply With Quote

Old   July 3, 2020, 03:24
Default
  #5
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Thank you. I uploaded my test files and mesh here. Could you please help me to have a look? By the way, the parameters in slurry pipe I used are the following:

pipe diameter 0.1 m (only half)
particle diameter 1e-4 m
Uin=2.5 m/s
alphas=0.2


Quote:
Originally Posted by rdbisme View Post
Hello qi.yang, I think you might get more helpful comments if you give more details about your configuration. You can check my old post for an example of which details might be interesting.
Attached Images
File Type: png WeChat Image_20200703091844.png (65.3 KB, 11 views)
qi.yang@polimi.it is offline   Reply With Quote

Old   July 3, 2020, 03:29
Default
  #6
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
I uploaded the files here without mesh, due to the large size.
Attached Files
File Type: zip slurrypipe.zip (18.8 KB, 1 views)
qi.yang@polimi.it is offline   Reply With Quote

Old   July 3, 2020, 03:38
Default
  #7
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
I saw your previous post and I checked your test case uploaded there. The timestep in your case is 2e-4 s. So you need several days or more to finish one simulation?

Quote:
Originally Posted by rdbisme View Post
Hello qi.yang, I think you might get more helpful comments if you give more details about your configuration. You can check my old post for an example of which details might be interesting.
qi.yang@polimi.it is offline   Reply With Quote

Old   July 6, 2020, 10:30
Default
  #8
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Nobody is curious about this issue?
qi.yang@polimi.it is offline   Reply With Quote

Old   July 6, 2020, 11:38
Default
  #9
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
Quote:
Originally Posted by qi.yang@polimi.it View Post
I saw your previous post and I checked your test case uploaded there. The timestep in your case is 2e-4 s. So you need several days or more to finish one simulation?
Hello,


actually.. no. Running the test case in parallel the test would finish in a matter of hours... Why do you say days? Are you running in parallel?
rdbisme is offline   Reply With Quote

Old   July 7, 2020, 08:27
Default
  #10
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Thanks. I did not try parallel but in the end I will do it. I am curious that there is the way to increase the time step.


Quote:
Originally Posted by rdbisme View Post
Hello,


actually.. no. Running the test case in parallel the test would finish in a matter of hours... Why do you say days? Are you running in parallel?
qi.yang@polimi.it is offline   Reply With Quote

Old   July 7, 2020, 08:32
Default
  #11
Senior Member
 
rdbisme's Avatar
 
Ruben Di Battista
Join Date: May 2013
Location: Paris
Posts: 137
Rep Power: 13
rdbisme is on a distinguished road
With a very basic explanation: The time step of a simulation can't be arbitrarily changed. Based on the time and space schemes you are using, there's a so-called CFL condition that restricts the maximum time step you can achieve. This CFL condition takes also in consideration the cell size: that means if you have a very fine mesh, you're generally gonna have a small time step to keep the numerical scheme stable.

What you can do to increase the time step is to compromise on some aspects:
- Use a coarse mesh
- Use low order numerical schemes that have larger stability regions
- Use implicit time schemes
- Increase the Courant number (that will reduce your stability)
granzer likes this.
rdbisme is offline   Reply With Quote

Old   July 7, 2020, 08:42
Default
  #12
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Hi guys, I met a new problem that the converency issue between the unstrcture mesh created by workbench and structure mesh generated by blockMesh. Why the iteration of P_rg is 1000, which means unconverged, right? Do you habe any suggestion? Thank u.
Attached Images
File Type: png structuremesh.png (26.2 KB, 4 views)
File Type: png unstructuremesh.png (31.5 KB, 3 views)
File Type: png structuremesh1.png (51.8 KB, 8 views)
File Type: png unstructuremesh1.png (70.0 KB, 10 views)
qi.yang@polimi.it is offline   Reply With Quote

Old   July 7, 2020, 08:52
Default
  #13
Senior Member
 
kimy
Join Date: Mar 2019
Location: https://t.me/pump_upp
Posts: 164
Rep Power: 7
qi.yang@polimi.it is on a distinguished road
Send a message via ICQ to qi.yang@polimi.it Send a message via AIM to qi.yang@polimi.it Send a message via Yahoo to qi.yang@polimi.it
Thanks a lot. I already considered what you mentioned below however, the time step is the samll one. Now I regard that issue as normal one. Really thank u. Gentilissimo!

Quote:
Originally Posted by rdbisme View Post
With a very basic explanation: The time step of a simulation can't be arbitrarily changed. Based on the time and space schemes you are using, there's a so-called CFL condition that restricts the maximum time step you can achieve. This CFL condition takes also in consideration the cell size: that means if you have a very fine mesh, you're generally gonna have a small time step to keep the numerical scheme stable.

What you can do to increase the time step is to compromise on some aspects:
- Use a coarse mesh
- Use low order numerical schemes that have larger stability regions
- Use implicit time schemes
- Increase the Courant number (that will reduce your stability)
qi.yang@polimi.it is offline   Reply With Quote

Reply

Tags
timestep, twophaseeulerfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
courant number increases to rather large values 6863523 OpenFOAM Running, Solving & CFD 22 July 5, 2023 23:48
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 13:40
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM Running, Solving & CFD 8 March 28, 2016 11:08
Floating point exception error lpz_michele OpenFOAM Running, Solving & CFD 53 October 19, 2015 02:50
Star cd es-ice solver error ernarasimman STAR-CD 2 September 12, 2014 00:01


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