SteadyState and Transient Solvers
Could someone please explain the differences between steadystate solvers and Transient solvers also when is best to use each (specifically pisoFoam and SimpleFoam).
So far I understand that SteadyState is much faster and transient uses a timestep deltaT. But I am struggling to find any information past this. Thanks in advance. 
Greetings Thomas,
The terms steady state and transient have to do on if time is a parameter or not (very roughly). This has to do with the phenomena that you are studying.....if your physical flow does not change over time then you can see it as steady state. So in either case you must know before you do the simulation if your "real" flow is steady or transient and then procceed with the solvers you talked about. Poiseuille flow for example is treated as steady state. Regards, Nickolas 
Thanks for clearing this up for me. One thing remains though, how can turbulent flow be considered steady state?
Tom 
Well in the way you are posing this is very general. Maybe you can tell some information on what you have in mind to simulate. For example, if you process is a continuous process but the phenomena created are turbulent (very high Reynolds number) you can approximate it as steady state (since it is continuous) and insert a turbulent model such the Reynolds averaged stresses.
But all these are very general assumptions. With simpleFoam solver, which is steady state, you can switch between laminar and turbulent flow. Nickolas 
I am doing air flow over a small air vehicle. Reynolds 120,000 15deg angle of attack. Would I be right in saying that due to the air flow separation and turbulence this would be a noncontinuous simulation?

time evolution PisoFoam
Hello there,
I am struggling a little bit with the meaning of time in transient solvers like PisoFoam. I am new to openfoam and I am checking pisofoam for a transient newtonian flow. The problem is that whenever I decrease the time step on a pressuredriven flow the steady state is reached faster and the time scale changes dramatically. When I compare my results with analytic ones I find drastic differences, like if time would be scaled by a factor. Does someone know if the time that openFoam prints has a physical direct meaning?:confused::confused::confused: 
Greetings to all!
@Heliana: I got your private message and went looking for any posts you had made about this question. There are a few very important details that I need to ask, in order to assess what might be the problem:
Bruno 
Hello Bruno!!
I am so sorry for the late reply, i was on holidays :S My case is a pressure driven Newtonian flow. This is how I defined the mesh: Code:
/** C++ **\ These are the answers: 1. I studied the three tutorials :) 2. in the cavity case I compared also time scales between simpleFoam and pisoFoam and I found time differences regarding the time needed to get to steady state. 3. I checked both simpleFoam and pisoFoam for this case. Always simpleFoam would take importantly more time to get to steady state than pisoFoam (by a factor of approx. 100). 4. what do you mean by flow speeds? 5. these are my files more or less for pisoFoam controlDict Code:
FoamFile Code:
FoamFile Code:
/** C++ **\ Code:
/** C++ **\ Sorry uff it all sounds like too much. From simpleFoam to pisoFoam I changed the ddtSchemes from steadyState to Euler. in simpleFoam the div(phi,U) is bounded Gauss upwind and in pisoFoam it is just gauss upwind and fvSolution is the same. Actually the most different is fvSchemes for simpleFoam: Code:
/** C++ **\ sorry for the overload!! :s have a nice day!! heliana 
Hi Heliana,
Quote:
Code:
convertToMeters 0.01; With a "MaxCo" set to 0.7, the maximum velocity should be... I haven't done the math, but likely the maximum allowed velocity in the "U" field would be something like 0.1 to 0.01 m/s. Quote:
I don't have time to setup these cases myself, but I assume you're basing yourself on the tutorial "incompressible/pisoFoam/ras/cavity", correct? If you could share the cases, it would it easier for anyone to see and test this. Best regards, Bruno 
Hello Bruno,
Thank you for answering. Well I am using a better mesh, more uniform, and I am actually working more studying a pressure driven flow through a simple pipe. I am also very much aware of the convertToMetres command, though I tried to simplify more my mesh file. Now that the mesh is better I still find differences regarding the evolution of velocity when I decrease the timeStep / Courant number. Does that make sense? I mean for every timestep I use the velocity I am waiting for is eventually reached, but with low timesteps (1e6) the steady state is reached faster which tends to be similar to the analytic result. What do you think? Thank you for your patience :) 
Hi Heliana,
Quote:
Because if you are enforcing a fixed time step and are not careful with keeping the Courant Number below the 0.5 mark, then your solution is simply wrong. There is a specific line that tells you this, namely the one mentioning the word "continuity". If you look at that line in each time step, what values are you getting? Best regards, Bruno 
Hi Bruno!
Well, I have been using really low time steps, in the order of 10^6 amd for the time continuity errors I get something like this: time step continuity errors : sum local = 1.17661e22, global = 9.15827e25, cumulative = 3.50451e17 DICPCG: Solving for p, Initial residual = 1.33849e07, Final residual = 1.33849e07, No Iterations 0 time step continuity errors : sum local = 3.74667e22, global = 9.26575e25, cumulative = 3.50451e17 Which sounds pretty low. I am going to adjustTimeStep to see what happens :) I will let you know then ! Cheers, Heliana 
A more basic way of classifying the differences is this:
If you care about observing timeaccurate flow development, then you would use a transient solver (pisoFoam, pimpleFoam, etc) with a 2ndorder temporal scheme. If you only care about reaching a "steadystate" and do not wish to observe anything that goes on between the beginning and end states, then it is much better to choose one of the steadystate solvers (simpleFoam etc). Transient solvers generally go handinhand with DES and LES simulations, while steadystate solvers are better coupled with RANS simulations. 
Quote:

Quote:
One example where having a 2ndorder temporal scheme is important is when predicting the Strouhal number (St) of oscillation in the flow past a sphere using DES or LES. This is simply a measure of wake shedding frequencies. If you choose Euler, at the very least you'll certainly experience a lag/errors in the time history of St with respect to experiment. 
Quote:

Quote:
In cases where backward is unstable, you may desire to use CrankNicolson and blend it with Euler (using the coefficient 0 to 1) if that doesn't work either. 
Hi guys, thank you so much for answering, it has been really helpful... Well, I was checking my files and I found a mistake I have made (i feel really embarrassed). I was using relaxation factors so my transients were affected... when I removed them it worked just fine.
My goal is actually use the viscoelasticFluidFoam solver using a modified Maxwell model, and that is suposed to be a transient solver as far as I understood. Is any of you familiar with that solver?? cheers ! 
Is that a general principlethat relaxation factors should not be used on unsteady problems?

Quote:

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