
[Sponsors] 
March 9, 2018, 05:48 
Steady state solution > transient

#1 
Member
Nat K
Join Date: Oct 2017
Posts: 68
Rep Power: 7 
Hi,
I am looking to use a laminar steady state solution as an initial condition for a laminar transient case. I know it is possible to do so on CFX, but unsure how to do it on OpenFOAM. I have done the steady state case with simpleFoam, and I guess I need to change the solver to pisoFoam for transient, and the ddt scheme to a transient scheme. I am a bit unsure about the startFrom and startTime, because the steady state case is in iterations, but a transient case is in time steps. I have read through some posts but I can't find one where someone has managed to do this (please let me know if there is a thread about it) I think I can use mapFields but unsure how to use if I have the same mesh. Most probably more straighforward than I am thinking. Also is it possible to do this with a decomposed case? 

March 9, 2018, 12:00 

#2 
Member
Arvind Jay
Join Date: Sep 2012
Posts: 96
Rep Power: 13 
just change the startFrom flag in the controlDict also with the solver.
Code:
startFrom latestTime; Jay 

March 12, 2018, 09:25 

#3  
Senior Member

Quote:
OP; You need to first change the delta T to a smaller value representing courant number <1 if pisoFoam is your choice of solver. YOu need to change the fvSchems time marching scheme from steady state to Euler, backward or Crank nickolson (or etc). If it is the same mesh you don't need to map the results as well just take the saved data for the latest time step in simpleFoam. Inside the data files there will be a file called uniform > time. Edit the index of time to one, time delta value to the desired delta T value. From there I think your simulation can start. 

March 12, 2018, 09:50 

#4  
Member
Nat K
Join Date: Oct 2017
Posts: 68
Rep Power: 7 
Quote:
Thanks for your answer. I've change fvSchemes to make it transient, and fvSolution for PISO loop. I have set my required timestep 1e5 and changed the solver to pisoFoam. So in the last time directory e.g. 100. I change index  1, deltaT  1e5 and deltaT0  1e5? And in the controlDict I just change startTime to latestTime (which should be 100 directory)? 

March 12, 2018, 11:49 

#5  
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,246
Rep Power: 63 
Quote:
Only thing you might want to do is change 100 to 1 or 0. But this is not absolutely necessary because OF will just think that the current time is 100s. When you do startFrom latestTime; OF runs a utility that scans your case directory and looks for the biggest time. You can also do: Code:
startFrom startTime; startTime 100; This is same for transient as it is for steady state except that the meaning of time is iteration for steady state and time for transient. 

Tags 
laminar, map fields, steady state 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
2way FSi, Initialize with steady solution, Fluent, Transient Sturcural, System Coupli  mmkkeshavarzi  FLUENT  0  June 22, 2016 08:26 
Steady state solution in chtMultiRegionSimpleFoam depends on kappa (t. conductivity)?  karoltomek  OpenFOAM Running, Solving & CFD  3  May 8, 2014 01:54 
is it possible to predict how long it takes to reach steady state solution in unstead  Alimohamadi_nasr  CFX  4  November 11, 2013 06:11 
time to reach steady state with transient simulati  prabhu  CFX  2  April 21, 2008 10:18 
About the difference between steady and unsteady problems  Lisa  Main CFD Forum  11  July 5, 2000 14:37 