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

How do set a steady solution as an initial solution to an unsteady simulation?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 18, 2020, 06:39
Default How do set a steady solution as an initial solution to an unsteady simulation?
  #1
New Member
 
pro
Join Date: Mar 2020
Posts: 25
Rep Power: 6
pro_ is on a distinguished road
In the study of the aeroelastic Isogai Case A such as in this paper,
https://kiwi.oden.utexas.edu/papers/...es-Willcox.pdf

It states that, using SU2, they first initialized a steady state solution for an unsteady simulation. How do i go about in doing this (in general)?

Do i need to make 2 configuration files, one that is for the initial steady and the other for the unsteady?

I have tried this, and set the unsteady config file to restart from the steady config but it gives me an error, where it can't find the restart file. So i assume I can do this all in one config file?

I have checked the tutorial on Unsteady NACA 0012, where i assume they first initialized a steady state solution. But attached in the picture is where i dont understand how they did the following.


Attached is an example of my config file
Attached Images
File Type: jpg Capture.JPG (121.6 KB, 38 views)
Attached Files
File Type: txt 0_config.txt (10.2 KB, 11 views)
pro_ is offline   Reply With Quote

Old   April 18, 2020, 10:31
Default
  #2
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 13
pcg is on a distinguished road
I don't think you can do it in one config file.
I would just rename and duplicate the steady restart according to the time-domain naming convention, and then start the unsteady simulation.
Note that the first time steps of such an unsteady simulation are not useful, you still need to wait for transients to clear and for the simulation to start showing some kind of periodic behavior (starting from steady state hopefully makes this quicker...).

Regarding the passage you highlighted, I think it is not referring to this restart process, they are explaining the dual time-stepping strategy in simple terms (i.e. what goes on inside the solver, it is not something accessible from the user side).
pcg is offline   Reply With Quote

Old   April 18, 2020, 10:53
Default
  #3
New Member
 
pro
Join Date: Mar 2020
Posts: 25
Rep Power: 6
pro_ is on a distinguished road
Hi, I have renamed my steady state solutions (restart_flow) to solution_flow_00000. Running this would give me an error of 'negative iterator etc.' I mean i understand the DUAL_TIME stepping strategy would require 2 solution flows, I dont understand how i can create 2 steady state solutions (meaning i dont understand the logic behind this). So what i did was i copy pasted the solution_flow_00000.dat file and renamed it to solution_flow_00001.dat. I still get the error of a 'negative iterator'.

Ive attached my steady config file, hopefully you can help me out here.

Moreover, on the topic of the typical section model, I would like to ask. Is there anyway i can define the x position of the elastic axis in SU2? Based on the template only the Xcg is written.

Would the Elastic axis location be defined by the SU2 variable REF_ORIGIN_MOMENT_X? (for a typical section)

% Distance in semichords by which the center of gravity lies behind
% the elastic axis
CG_LOCATION= -0.2

If i can verify, does the CG_LOCATION defined above coincide with the attached picture?
Attached Images
File Type: jpg Capture.JPG (41.4 KB, 16 views)
Attached Files
File Type: txt 0_steady.txt (6.9 KB, 8 views)
pro_ is offline   Reply With Quote

Old   April 18, 2020, 11:23
Default
  #4
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 13
pcg is on a distinguished road
You need to set RESTART_ITER=2 (1 more than your last restart)
I cannot help with the aeroelastic features, I never used them.
pcg is offline   Reply With Quote

Old   April 20, 2020, 09:05
Default
  #5
New Member
 
pro
Join Date: Mar 2020
Posts: 25
Rep Power: 6
pro_ is on a distinguished road
Quote:
Originally Posted by pcg View Post
I don't think you can do it in one config file.
I would just rename and duplicate the steady restart according to the time-domain naming convention, and then start the unsteady simulation.
Note that the first time steps of such an unsteady simulation are not useful, you still need to wait for transients to clear and for the simulation to start showing some kind of periodic behavior (starting from steady state hopefully makes this quicker...).
.
Hi, you also mentioned removing the first few time steps
Referring to my attachment, would that mean removing the time history up to the blue line?
Attached Images
File Type: jpg Capture.JPG (53.6 KB, 24 views)
pro_ is offline   Reply With Quote

Old   April 20, 2020, 12:14
Default
  #6
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 13
pcg is on a distinguished road
A lot more than that, the characteristic time (period) seems to be around 0.08s, I would disregard around 10 periods, after which maybe you start seeing a clearly repeating cycle (same peak to peak amplitude, same average, etc.).
That is assuming what you are studying is periodic. If you want to study a short transient that is not periodic, e.g. a quick movement of a control surface, you still need an "unsteady initialization", the gist of it is, until the initial particles (from a steady state of freestream initialization) leave the domain the results are not good.
pcg is offline   Reply With Quote

Old   April 26, 2020, 09:15
Default
  #7
New Member
 
pro
Join Date: Mar 2020
Posts: 25
Rep Power: 6
pro_ is on a distinguished road
Quote:
Originally Posted by pcg View Post
A lot more than that, the characteristic time (period) seems to be around 0.08s, I would disregard around 10 periods, after which maybe you start seeing a clearly repeating cycle (same peak to peak amplitude, same average, etc.).
That is assuming what you are studying is periodic. If you want to study a short transient that is not periodic, e.g. a quick movement of a control surface, you still need an "unsteady initialization", the gist of it is, until the initial particles (from a steady state of freestream initialization) leave the domain the results are not good.
It seems that when i try restarting from my steady solutions, the unsteady solution quickly diverges. Any clue as to why this occurs? and how shall i get around this
pro_ is offline   Reply With Quote

Old   April 26, 2020, 19:23
Default
  #8
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 13
pcg is on a distinguished road
You have a moving grid problem, you are probably missing the grid velocities when you restart (similar to this problem Restart causes FGMRES orthogonalization failed).
Try using RESTART_STEADY_STATE= YES and be sure to use v7.0.3.
pcg is offline   Reply With Quote

Old   April 27, 2020, 12:54
Default
  #9
New Member
 
pro
Join Date: Mar 2020
Posts: 25
Rep Power: 6
pro_ is on a distinguished road
Quote:
Originally Posted by pcg View Post
You have a moving grid problem, you are probably missing the grid velocities when you restart (similar to this problem Restart causes FGMRES orthogonalization failed).
Try using RESTART_STEADY_STATE= YES and be sure to use v7.0.3.
I have tried that but obtained a ' linear solver diverged'. I think it has something to do with my initial steady state solution. So to clear things this is what i have done,

Ran the SAME steady state case TWICE and renamed each restart_flow file to solution_flow_00000 and solution_flow_00001, twice because my unsteady simulation is in DUAL_TIME.

I wonder if this is actually a correct procedure?
pro_ is offline   Reply With Quote

Old   April 28, 2020, 16:24
Default
  #10
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 13
pcg is on a distinguished road
I played around with your case and I think there is a bug somewhere preventing aeroelastic restarts from steady state. Likely because the grid velocities are not in steady restart files.
I'll investigate this a bit but for now the conclusion is there is a bug.

On another note your case will run much faster if you use ILU as the preconditioner for deformation (for the reasons explained here https://su2code.github.io/docs_v7/Li...econditioners/)
pcg is offline   Reply With Quote

Old   April 28, 2020, 18:05
Default
  #11
New Member
 
pro
Join Date: Mar 2020
Posts: 25
Rep Power: 6
pro_ is on a distinguished road
Quote:
Originally Posted by pcg View Post
I played around with your case and I think there is a bug somewhere preventing aeroelastic restarts from steady state. Likely because the grid velocities are not in steady restart files.
I'll investigate this a bit but for now the conclusion is there is a bug.

On another note your case will run much faster if you use ILU as the preconditioner for deformation (for the reasons explained here https://su2code.github.io/docs_v7/Li...econditioners/)
Thanks for taking a look into the case. Will take your advice
pro_ 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
[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend" bigphil OpenFOAM CC Toolkits for Fluid-Structure Interaction 686 December 22, 2022 10:10
pressure in incompressible solvers e.g. simpleFoam chrizzl OpenFOAM Running, Solving & CFD 13 March 28, 2017 06:49
Cannot run the code properly: very large time step continuity error crst15 OpenFOAM Running, Solving & CFD 9 December 14, 2014 19:17
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
pimpleFoam: turbulence->correct(); is not executed when using residualControl hfs OpenFOAM Running, Solving & CFD 3 October 29, 2013 09:35


All times are GMT -4. The time now is 08:15.