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

steady for unsteady flow?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 19, 2013, 05:22
Default steady for unsteady flow?
  #1
Member
 
Shenren Xu
Join Date: Jan 2011
Location: London, U.K.
Posts: 67
Rep Power: 15
Shenren_CN is on a distinguished road
Hello,

I have a quite fundamental question about using steady flow solver on (physically) unsteady flow.

Let's take the example of laminar flow around cylinder, and I want to find the bifurcation point, i.e., the Reynold's number above which flow will be unsteady with periodic vortex shedding. I run the steady flow solver with gradually increasing Reynold's number until the steady flow solver does not fully converge any more, presumably it goes to limit cycle oscillation.

My question is, how much would the demarcation point I would find through this numerical example reflect the physical phenomenon. And how would using an explicit or implicit steady solver affect the demarcation point I find? And finally, my numerical experience tells me that demarcation point, i.e., the critical Reynold's number using implicit solver is higher than the explicit solver, why is there a difference in the demarcation point I find using explicit and implicit solvers?

Thanks in advance for any shared thoughts on this.
Cheers,
Shenren
Shenren_CN is offline   Reply With Quote

Old   July 19, 2013, 05:27
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,756
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Shenren_CN View Post
Hello,

I have a quite fundamental question about using steady flow solver on (physically) unsteady flow.

Let's take the example of laminar flow around cylinder, and I want to find the bifurcation point, i.e., the Reynold's number above which flow will be unsteady with periodic vortex shedding. I run the steady flow solver with gradually increasing Reynold's number until the steady flow solver does not fully converge any more, presumably it goes to limit cycle oscillation.

My question is, how much would the demarcation point I would find through this numerical example reflect the physical phenomenon. And how would using an explicit or implicit steady solver affect the demarcation point I find? And finally, my numerical experience tells me that demarcation point, i.e., the critical Reynold's number using implicit solver is higher than the explicit solver, why is there a difference in the demarcation point I find using explicit and implicit solvers?

Thanks in advance for any shared thoughts on this.
Cheers,
Shenren

I would consider a different stategy: use always time-marching scheme, changing the Re number you will find that the steady solution is no longer reached.

Implicit or explicit schemes of the same accuracy order must reach the same solution, provided that you are using the same time-step. Some difference can be due to the fact that implicit schemes require some particular BC at the new update and this task can introduce some approximation that does not appear for explicit scheme
FMDenaro is offline   Reply With Quote

Old   July 19, 2013, 05:34
Default
  #3
Member
 
Shenren Xu
Join Date: Jan 2011
Location: London, U.K.
Posts: 67
Rep Power: 15
Shenren_CN is on a distinguished road
Hi Filippo.

I think what I did is the same as what you proposed. I used a (pseudo) time-marching scheme and change the Re number until it does not converge.

My explicit solver uses Runge-Kutta+Multigrid+block-Jacobi preconditioner, so the time-accuracy is not preserved. By "time-marching scheme", are you actually suggesting to use a unsteady solver with physical time step?

As for implicit solver, for me, it's just a matrix preconditioner on the LHS, should it not affect the steady solution that's reached, if ever?

Please let me know if I didn't explain clearly.
Shenren



Quote:
Originally Posted by FMDenaro View Post
I would consider a different stategy: use always time-marching scheme, changing the Re number you will find that the steady solution is no longer reached.

Implicit or explicit schemes of the same accuracy order must reach the same solution, provided that you are using the same time-step. Some difference can be due to the fact that implicit schemes require some particular BC at the new update and this task can introduce some approximation that does not appear for explicit scheme
Shenren_CN is offline   Reply With Quote

Old   July 19, 2013, 05:53
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,756
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
yes, I mean a time-marching integration using a physical time-step.

Are you sure that in your implicit method the BCs for Un+1 is never required? However, implicit and explicit integration of the same accuracy, using the same time step, have the same magnitude order in the LTE, differences appearing in the sign of the coefficients...
FMDenaro is offline   Reply With Quote

Old   July 19, 2013, 06:06
Default
  #5
Member
 
Shenren Xu
Join Date: Jan 2011
Location: London, U.K.
Posts: 67
Rep Power: 15
Shenren_CN is on a distinguished road
Okay, then it's a different approach, and essentially you are saying I should use an unsteady solver. Maybe there is a fundamental problem with my approach of using only steady solver, I will think about it carefully.

I'm not using time-dependent boundary condition, so no, the BCs for Un+1 is never required. And by 'LTE', you mean leading truncation error?


Quote:
Originally Posted by FMDenaro View Post
yes, I mean a time-marching integration using a physical time-step.

Are you sure that in your implicit method the BCs for Un+1 is never required? However, implicit and explicit integration of the same accuracy, using the same time step, have the same magnitude order in the LTE, differences appearing in the sign of the coefficients...
Shenren_CN is offline   Reply With Quote

Old   July 19, 2013, 06:17
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,756
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
yes, LTE = local truncation error

Note that when the solution becomes unsteady at some Re number, you have to prescribe a suitable BC Un+1 at the outlet, otherwise you somehow will force the steady condition (due to elliptical part of the solution).
FMDenaro is offline   Reply With Quote

Old   July 19, 2013, 06:26
Default
  #7
Member
 
Shenren Xu
Join Date: Jan 2011
Location: London, U.K.
Posts: 67
Rep Power: 15
Shenren_CN is on a distinguished road
Thanks. I see. This approach probably will give me a better estimate
of the ctirical Re number. Then again, why using the explicit solver and
the Re number found by checking whether it converges or not does not
give me a satisfying answer? I'm more interested in this more 'theoretical'
problem than the more practical problem of predicting the critical Re #.

Any thought on that? :P


Quote:
Originally Posted by FMDenaro View Post
yes, LTE = local truncation error

Note that when the solution becomes unsteady at some Re number, you have to prescribe a suitable BC Un+1 at the outlet, otherwise you somehow will force the steady condition (due to elliptical part of the solution).
Shenren_CN is offline   Reply With Quote

Old   July 19, 2013, 06:32
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,756
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Shenren_CN View Post
Thanks. I see. This approach probably will give me a better estimate
of the ctirical Re number. Then again, why using the explicit solver and
the Re number found by checking whether it converges or not does not
give me a satisfying answer? I'm more interested in this more 'theoretical'
problem than the more practical problem of predicting the critical Re #.

Any thought on that? :P

Many issues can be involved...it is difficult for me thinking about without checking details ... also the chosen tolerance for the residuals can cause differences...
FMDenaro is offline   Reply With Quote

Reply

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
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
Steady state Ruben Main CFD Forum 43 May 7, 2011 04:32
Unsteady AND Steady mode for Fully Developped Flow Dominique FLUENT 5 April 16, 2004 18:58
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19
steady and unsteady flow shisheng Main CFD Forum 10 July 23, 1999 07:00


All times are GMT -4. The time now is 03:29.