CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > SU2

Unsteady 2nd order Dual Time-Step Diverges When Restarted From Steady Solution

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

Like Tree2Likes
  • 1 Post By pcg
  • 1 Post By pcg

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 26, 2021, 11:27
Default Unsteady 2nd order Dual Time-Step Diverges When Restarted From Steady Solution
  #1
Member
 
Sangeet
Join Date: Jun 2017
Location: India
Posts: 43
Rep Power: 8
sangeet is on a distinguished road
Hello,


I am currently running trying to run a unsteady case in SU2 v7.1.1 and I want to restart from a converged steady solution that I already have. I used the 2nd order dual time-step method. Since I need two solutions to restart, I just made two copies of the steady solution, but SU2 diverges before any iterations start. Below is the unsteady paramaters part of the config file.


This google drive link has my case files.
https://drive.google.com/drive/folde...J8SK9aW2xwMNeq



Code:
% ------------------------- UNSTEADY SIMULATION -------------------------------%
%
TIME_DOMAIN = YES
%
% Numerical Method for Unsteady simulation(NO, TIME_STEPPING, DUAL_TIME_STEPPING-1ST_ORDER, DUAL_TIME_STEPPING-2ND_ORDER, TIME_SPECTRAL)
TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER
%
% Time Step for dual time stepping simulations (s) -- Only used when UNST_CFL_NUMBER = 0.0
%TIME_STEP= 5e-4
%
% Maximum Number of physical time steps.
TIME_ITER= 100000

% Total Physical Time for dual time stepping simulations (s)
MAX_TIME= 0.1

% Unsteady Courant-Friedrichs-Lewy number of the finest grid
UNST_CFL_NUMBER= 1
%
% Number of internal iterations (dual time method)
INNER_ITER= 100
%
% Restart after the transient phase has passed
RESTART_SOL = YES
%
% Specify unsteady restart iter
RESTART_ITER = 2

Last edited by sangeet; April 26, 2021 at 12:14. Reason: Added SU2 version
sangeet is offline   Reply With Quote

Old   April 26, 2021, 18:37
Default
  #2
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 465
Rep Power: 13
pcg is on a distinguished road
This UNST_CFL_NUMBER= 1 gives you a variable time step which is not correct for the second order dual time approach, because the time derivative is based on finite differences.
Set UNST_CFL_NUMBER to 0 and use TIME_STEP instead.
pcg is offline   Reply With Quote

Old   April 26, 2021, 18:41
Default
  #3
Member
 
Sangeet
Join Date: Jun 2017
Location: India
Posts: 43
Rep Power: 8
sangeet is on a distinguished road
Hello,
Oh, so specifying the CFL number makes it like local time stepping? But then what is its role for an unsteady problem?


EDIT : and then why does this config work when I am not restarting?

Last edited by sangeet; April 26, 2021 at 18:51. Reason: Added question
sangeet is offline   Reply With Quote

Old   April 26, 2021, 19:32
Default
  #4
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 465
Rep Power: 13
pcg is on a distinguished road
Hi,
You have a point... I looked at the code again and it looks like a bug.
For TIME_STEPPING the unsteady CFL is used at each time iteration to determine the time step.

For DUAL_TIME_x it looks like the intention is to determine an automatic time step on the first iteration, and then fix it (because it cannot be variable for these methods) however the logic for this seems to be broken when the simulation is restarted.

I'll see about fixing this, but for now set it to zero and use TIME_STEP.

In dual time methods, the pseudo time step determined for inner iterations is based on the normal CFL, this UNST_CFL_NUMBER option is only for the "real" time step.
That is, dual time is just steady state with a source term.
sangeet likes this.
pcg is offline   Reply With Quote

Old   April 27, 2021, 03:10
Default
  #5
Member
 
Sangeet
Join Date: Jun 2017
Location: India
Posts: 43
Rep Power: 8
sangeet is on a distinguished road
Hello,
Thank you for the clarification. I had a question. When you say that the time step is determined by the cfl number, do you mean that su2 chooses the minimum timestep available as the cell sizes are different and the timesteps calculated from them would all be different?
sangeet is offline   Reply With Quote

Old   April 27, 2021, 13:04
Default
  #6
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 465
Rep Power: 13
pcg is on a distinguished road
Yes correct, minimum timestep over the entire problem.
sangeet likes this.
pcg is offline   Reply With Quote

Old   May 18, 2021, 05:56
Default
  #7
Member
 
Ravi
Join Date: May 2017
Posts: 31
Rep Power: 8
raviramesh10 is on a distinguished road
Quote:
Originally Posted by pcg View Post
Yes correct, minimum timestep over the entire problem.
Hi pcg, as I am also working on unsteady flow problems, I had a question regarding this too. But wouldn't this time step depend on the local grid size as well through the CFL number? If yes, does that mean that it takes the minimum grid size into account for calculating the time step?
raviramesh10 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
Time Step Continuity Errors simpleFoam Dorian1504 OpenFOAM Running, Solving & CFD 1 October 9, 2022 10:23
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field lakeat OpenFOAM Community Contributions 58 December 23, 2021 03:36
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 06:28
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 14:40
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 08:47


All times are GMT -4. The time now is 22:54.