
[Sponsors] 
December 22, 2012, 21:02 
Simplefoam, steady or unsteady

#1 
Member
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 6 
Hi FOAMERS
SIMPLEFOAM is steady, (due to the manual), but when I am running pitzDaily with it, why do I get a solution at times t = 0 to 1000 s? 

December 22, 2012, 23:29 

#2 
Senior Member

Hi,
The original simpleFoam solver is for the steady calculations. In steady calculations, deltaT has no effect on the results and is usually set to 1 so that the time means the number of iterations. Hope that helps, fumiya 

February 20, 2013, 05:59 

#3 
Member
Malik
Join Date: Dec 2012
Location: Austin, USA
Posts: 52
Rep Power: 5 
Hi,
If simpleFoam is for steady calculations, then why is it possible to change the ddtScheme into "Euler" rather than "steadyState" in system/fvSchemes ? is it ok to apply this change ? When I tested this on a cylinder at RE = 100 it gave me the same results for Cd and strouhal. By the way, the results are far from the experiments. With the solver piso and ico, the same mesh and boundary conditions, it is better. It makes me think that the simple Solver solver should never be used as a transient solver. Am i right ? 

February 20, 2013, 06:09 

#4 
Senior Member
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,211
Blog Entries: 1
Rep Power: 17 
simpleFoam is an steadystate solver and those changes has no effect on result
__________________
Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog (http://openfoam.blogfa.com/) Training Course on OpenFOAM at (http://www.isme.ir/) 

February 20, 2013, 06:15 

#5 
Member
Malik
Join Date: Dec 2012
Location: Austin, USA
Posts: 52
Rep Power: 5 
Thank you, it makes sense now.
It is quite confusing to let the user write Euler in the ddtSchemes. It would be better to fix this scheme as steadyState. 

February 20, 2013, 10:04 

#6  
Senior Member

Quote:
i don't know exactly how openFOAM deals with steady state but, one approach is to consider it as transient called (pseudo transient) it add transient term (d/dt) and start run from sensible values like starting from rest (it is easier to get these values for transient than steady state). it deals with transient term in different way from piso algorithm to accelerate solution. finally, you should ignore all values in developing time and only take into consider values after developing time when values become fixed. good luck 

March 14, 2014, 02:51 

#7 
Member
Fengjiao Bian
Join Date: Nov 2013
Location: beijing
Posts: 30
Rep Power: 5 
Hi ,Fumiya! I want to use LES turbulence model to simulate steady state of wind turbine with smagorinsky model. Could you tell me if I could choose the SimpleFoam solver?I see there are some articles with the method of LES in steady state but not in OF,So I am wondering how to achieve it. I know that LES method are usually used in trasient condition, I want to make sure does the LES inconsistent with simplefoam solver because when I run the MRFSimpleFoam ,it shows the error cannot find the file RASProperties. Could you help me ?


March 14, 2014, 03:07 

#8 
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 14 
It does not make any sense to do steady state LES simulations. I would be interested in any articles showing this. Of coarse, you can run a transient solver, on a stable, physically steady case, but this does not make it a steady state solver.


March 14, 2014, 03:31 

#9 
Member
Fengjiao Bian
Join Date: Nov 2013
Location: beijing
Posts: 30
Rep Power: 5 
《SteadyState LargeEddy Simulations to Study the Stratocumulus to Shallow Cumulus Cloud Transition 》 http://journals.ametsoc.org/doi/abs/...urnalCode=atsc


August 11, 2014, 07:37 

#10 
Senior Member
Join Date: Nov 2012
Location: Bavaria
Posts: 145
Rep Power: 6 
Hi jiaojiao,
did u find an answer on doing LES for a steady state flow? Ayla 

March 4, 2015, 05:43 
simpleFoam steady or unsteady

#11 
Member
Naresh Yathuru
Join Date: Feb 2015
Posts: 62
Rep Power: 3 
Hi Foamers,
I know this can be a silly question. I know that simple foam is a steady state solver and the delta T means the iteration number. When i run my simulation i get this Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 7.827025e13, No Iterations 20 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 5.0458368e13, No Iterations 22 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 3.9496433e13, No Iterations 25 GAMG: Solving for p, Initial residual = 1, Final residual = 3.6324613e11, No Iterations 1000 time step continuity errors : sum local = 6.3981153e14, global = 3.4832047e17, cumulative = 3.4832047e17 smoothSolver: Solving for omega, Initial residual = 0.17479046, Final residual = 5.3209563e13, No Iterations 20 smoothSolver: Solving for k, Initial residual = 1, Final residual = 7.9795701e13, No Iterations 22 ExecutionTime = 104.1 s ClockTime = 109 s So what does No iterations mean here? I think in trasient solvers there are inner iteration loops . This is confusing me could some one help me. I tried to go through the simpleFoam.c but i couldnt crack it. Regards, Naresh 

March 4, 2015, 06:00 

#12 
Senior Member
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,211
Blog Entries: 1
Rep Power: 17 
Those iteration numbers are for solution convergence in each time iteration
__________________
Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog (http://openfoam.blogfa.com/) Training Course on OpenFOAM at (http://www.isme.ir/) 

March 4, 2015, 07:22 

#13 
Member
Naresh Yathuru
Join Date: Feb 2015
Posts: 62
Rep Power: 3 
Thank you so much Nima for the quick reply.
Just to clarify , please correct me if i m wrong. 1. These iterations depends on the tolerance value i specify in the FVsolution. And the tolerance is also global so it stops the simulation when the required convergence is achevied. solvers { p { solver GAMG; tolerance 1e06; relTol 1e06; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; } "(UkepsilonomegaRnuTilda)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e06; relTol 1e06; } } now i have a doubt if "reltol" plays a role in the inner iterations. Thank you for your time. with regards, Naresh 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
conversion steady to unsteady  Clementhuon  CDadapco  7  January 26, 2012 05:33 
steady to unsteady  wlt_1985  FLUENT  6  December 4, 2010 17:17 
Steady needs unsteady.  nico  Main CFD Forum  0  September 21, 2007 04:50 
Turbulent: Steady or Unsteady: confusion  prem  FLUENT  0  March 30, 2006 10:40 
steady or unsteady? (in dpm)  winnie  FLUENT  1  April 28, 2003 11:30 