Steady vs. Unsteady solution
Dear foamers,
I do have a case and I would like to study both the steady and the unstedy solution. Something strange happens, since I perfectly get the steady solution, but the unsteady one suddenly diverges. I think it is a problem in the `controlDict' file, am I right? Do you have any suggestion for me? Do you need some files in order to help? Thanks a lot, Samuele |
It's a bit difficult to tell from the information provided what the cause could be. Therefore:
Which solvers are you using? How big is the courant number? Kind regards, L |
I am using bouyantPimpleFoam for the unsteady solution and bouyantSimpleFoam for the steady one.
These are my 2 controlDicts files: Code:
/*--------------------------------*- C++ -*----------------------------------*\ Code:
/*--------------------------------*- C++ -*----------------------------------*\ |
The controlDict seems fine. Maybe you can lower the maxCo to 0.1 or even lower and see what happens, just to exclude this as possible cause although I doubt this will make a difference.
How do you make the difference between the steady and the unsteady case exactly? Do you use a RANS turbulence model for the former and a LES model for the latter? If you use the same RANS model for both, the transient simulation will in the end result in the same solution as the steady simulation so you can always try the steady solution as initial condition for the transient simulation and see what happens. Can you should the output of the simulation where the solver starts to diverge? Regards, L |
Both for the steady and the unsteady simulation, I use the RANSModel (below my turbulenceProperties and RASProperties files):
Code:
/*--------------------------------*- C++ -*----------------------------------*\ Code:
/*--------------------------------*- C++ -*----------------------------------*\ Just one detail to make the simulation faster (so that I can give more tries): which are good `numbers' for the outer and non-orthogonal correctors? I guess that mine are too big (see below). Code:
/*--------------------------------*- C++ -*----------------------------------*\ Thanks a lot, Samuele |
15 is indeed really big. You should be able to see it from the solver output how many you need. 2 or 3 is often recommended but I don't have a lot of experience with this so you might want to check out other sources as well.
Since you're running a transient simulation, I would change the relaxation to 1.0 for all (although this will not help your problem, probably rather the opposite). Did you already try to run the case with the steady state solution as initial condition? Also post the solver output next time, this might tell us something as well... |
Here is the error message:
Code:
GAMG: Solving for p_rgh, Initial residual = 9.8126427e-07, Final residual = 9.8126427e-07, No Iterations 0 Quote:
Thanks a lot, Samuele |
All times are GMT -4. The time now is 07:02. |