CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

min time step in transient CFD

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 29, 2017, 05:28
Default
  #21
Senior Member
 
A. Min
Join Date: Mar 2015
Posts: 305
Rep Power: 12
alimea is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
For an explicit time-marching scheme, the numerical stability region is determined by the curve (in the 1D case) cfl=f(Re_h), being Re_h the cell Reynolds number. Only for Re_h >>1 you can recover the constraint due to only the cfl value.
In your case, at Re=100 I suppose you are working at Re_h=O(1) so that the max cfl value for the stability is much lower than it would be for the inviscid case.
yes, my Re_h is in O(1). that's almost 0.85.
please introduce me a reference to see this curve (cfl=f(Re_h)) and how to obtain that.
Also, what's the meaning of "recovering the constraint due to only the cfl value" ?
alimea is offline   Reply With Quote

Old   December 29, 2017, 05:50
Default
  #22
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,793
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
That is a classical issue in CFD and numerical analysis and requires the use of the von Neumann analysis. Many texbooks give details about that.
An example for a 2D case is shown in https://www.researchgate.net/publica...y-driven_flows
FMDenaro is offline   Reply With Quote

Old   December 29, 2017, 05:55
Default
  #23
Senior Member
 
A. Min
Join Date: Mar 2015
Posts: 305
Rep Power: 12
alimea is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
That is a classical issue in CFD and numerical analysis and requires the use of the von Neumann analysis. Many texbooks give details about that.
An example for a 2D case is shown in https://www.researchgate.net/publica...y-driven_flows
Thank you for participating in this discussion and sharing your valuable experiences.
FMDenaro likes this.
alimea is offline   Reply With Quote

Old   December 29, 2017, 22:40
Default
  #24
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,278
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by alimea View Post
Thanks for your nice answer. because it will force me to learn some new concept. could you please explain about these concepts, or give me link or reference to read about them:
1- dissipation (physical or numerical?). especially what is Rhie Chow dissipation?
2- what are implicit and explicit under relaxations? what is the difference between them?

why did you think that in my problem, viscosity becomes dominant term?
I am really short on time to explain these two things. Also you would have to do some search on it. For now just take it as a hint and try to investigate.

About the Rhie and Chow term, which is difficult to find from books as to what my comment was talking about here is very short hint.

Rhie and chow term is inversely propertional to diagonal of momentum matrix. Diagonal of momentum matrix is inversely propertional to time step size. For euler time stepping it would be ( density * volume / delta_T ).

So when delta goes to 0, it shall go to infinity and 1/Ap goes to 0.

So basically as delta_T reduces the Rhie and Chow term becomes weaker and weaker. After certain value, it could be too small at some parts of simulation that solution can diverge.
alimea likes this.
arjun is offline   Reply With Quote

Old   January 9, 2018, 00:36
Default
  #25
Senior Member
 
A. Min
Join Date: Mar 2015
Posts: 305
Rep Power: 12
alimea is on a distinguished road
Hi all
I want to solve the flow around the cylinder in viscoelastic fluid. At first I did it for Newtonian fluid in Re= 10 , 40 , 100 and for viscoelastic fluid in Re= 10 , 100 and obtained exact results. But when I change the Re to 40 for viscoelastic fluid, it became diverged!
I changed many parameters of solution:
urf = from 0.1 to 0.9 --> diverged!
corunt Nu. = from 0.1 to 0.9 --> diverged!
changing the kind and size of mesh --> diverged!
changing the kind and size of mesh --> diverged!
increasing the number of solving pressure correction Eqn from 2 to 20 --> diverged!
decreasing the min residuals to 1e-8 --> diverged!
changing discretization schemes of div terms: Gauss linear (central), Gauss upwind(1st), Gauss upwind(2nd), limitedLinear , QUICK --> diverged!

But I have done some ways to solve that:
1- I saw in this page that sb proposed to the other that decrease the residuals to 1e-19 !! I did it and my solution became converged! I can't analyze that! Is it possible?
2- In the other way I increased the number of solving pressure correction Eqn to 20, decreasing the min residuals to 1e-8
3- setting the time step to 0.001 instead of setting Cr=0.3

Now I don't know that my results are reliable or not!
Could you please tell me what happend that these solutions are appropriate for solve it?

Thanks
alimea is offline   Reply With Quote

Reply

Tags
time step size, transient

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
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field lakeat OpenFOAM Community Contributions 58 December 23, 2021 02:36
p_rgh initial residual no change with different settings manuc OpenFOAM Running, Solving & CFD 3 June 26, 2018 15:53
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 05:28
simpleFoam error - "Floating point exception" mbcx4jc2 OpenFOAM Running, Solving & CFD 12 August 4, 2015 02:20
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58


All times are GMT -4. The time now is 17:24.