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

Steady State Discretisation and Time Step

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

Like Tree1Likes
  • 1 Post By LuckyTran

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 25, 2021, 09:12
Default Steady State Discretisation and Time Step
  #1
Senior Member
 
Brett
Join Date: May 2013
Posts: 212
Rep Power: 14
Bdew8556 is on a distinguished road
Hi everyone,

So I know when solving a steady state problem we use a kind of pseudo time step to march towards convergence. What I want to know is how do we specifically discretise our equations and where is the timestep implemented in this process?
If all the differentials with respect to time can be cancelled due to the problem being steady state where does the time step come in??

B
Bdew8556 is offline   Reply With Quote

Old   August 25, 2021, 15:36
Default
  #2
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,739
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
There are unsteady/transient solvers that use real time-stepping. If you march these in time, eventually you may arrive at a steady state solution (if it even exists). These always have temporal discretization because they are unsteady solvers and terms cancel out in the solution whenever they feel like canceling out.

There are steady solvers that solve problems that don't depend on time. Steady solvers don't need pseudo time-stepping and most of them don't have any temporal discretization because the terms are cancelled out. But sometimes a steady solver will use a pseudo time-step. Are you asking how these work? The pseudo time-step is a form of implicit underrelaxation and isn't actually a form of a temporal discretization.
LuckyTran is offline   Reply With Quote

Old   August 25, 2021, 15:43
Default
  #3
Senior Member
 
Brett
Join Date: May 2013
Posts: 212
Rep Power: 14
Bdew8556 is on a distinguished road
Hey LuckyTran.

If you could explain how they work that would be cool!

So theres no temporal discretisation in steady problems?
Bdew8556 is offline   Reply With Quote

Old   August 25, 2021, 16:41
Default
  #4
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,739
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
Alright so for steady problems you have a discretized linear system that needs to be solved. You have an initial guess serving as an old solution. If you solve the linear system you get an updated solution.


With implicit underrelaxation, you can inject (by adding) a certain amount of the old solution into the linear system. When you solve this modified underrelaxed system instead, the updated solution is biased and tends more towards the old solution. This is equivalent to limiting the solution change by a local time-step. You can try and figure out how the weights are related to the local time-step (hint: it's related to the flow Courant number). Or.... you can just put the time-step into the weight factor and that's pseudo transient.
LuckyTran is offline   Reply With Quote

Old   August 25, 2021, 18:40
Default
  #5
Senior Member
 
Brett
Join Date: May 2013
Posts: 212
Rep Power: 14
Bdew8556 is on a distinguished road
care to indulge me?
Bdew8556 is offline   Reply With Quote

Old   August 25, 2021, 18:46
Default
  #6
Senior Member
 
Brett
Join Date: May 2013
Posts: 212
Rep Power: 14
Bdew8556 is on a distinguished road
For example when I set a physical timescale value in CFX, what exactly am I setting and where does it present itself in the equations (Simplified, say the Euler equations)
Bdew8556 is offline   Reply With Quote

Old   August 26, 2021, 00:05
Default
  #7
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,739
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
The first rule of talking about pseudo transient is that we don't talk about what the real transient PDE is. We start with the presumption and restrict ourselves to only steady pde's. The pseudo transient term looks like so, and you can do this for any steady pde. You just artificially add a transient term into the equation. And any time you inject stuff like this into the linear system, it is a form of implicit under/over-relaxation.



Formally there is a temporal discretization going on, but then you will ask me what is the difference between the transient equation and this steady one with an artificial timestep. Well in this example there isn't a lot. But you should forget that there is a transient form of the N-S equations as long as you are talking about pseudo transient unless you want to confuse yourself. Fluent/CFX uses a 1st order relaxation that looks very much like 1st order backward Euler being applied to a transient PDE. And that's true. But in general, you can artificially add higher order temporal derivatives out of nowhere to further stabilize your solver or accelerate convergence.


The reason the pseudo time-step is useful, is because you can choose the pseudo time-step size based on physical phenomenon that you know exists in your problem to stabilize particular physics (convection, diffusion, buoyancy, etc.) rather than using under-relaxation factors that range somewhere from 0 to 1.
ari003 likes this.
LuckyTran is offline   Reply With Quote

Reply

Tags
discretisation, numerical methods, steady, timestep

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
CrankNicolson time discretisation GerhardHolzinger OpenFOAM Programming & Development 3 November 15, 2019 08:27
higher order time and space discretisation for elliptic eq. with a transport eq. emjay Main CFD Forum 1 August 20, 2019 16:50
Time Step Problem in NS equations code eulerian Main CFD Forum 0 June 18, 2013 20:12
Large time Step Continuity Errors ozzythewise OpenFOAM Running, Solving & CFD 1 September 29, 2010 20:55
Time discretisation scheme for steady state case Andrey CFX 7 September 16, 2008 03:40


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