CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Can initial condition makes the solution diverge? (https://www.cfd-online.com/Forums/openfoam-solving/219471-can-initial-condition-makes-solution-diverge.html)

kim4545 July 28, 2019 15:02

Can initial condition makes the solution diverge?
 
1 Attachment(s)
Hello Foamer! I'm beginner of CFD.

I'm trying to simulate Supersonic nozzle with various Mach number.

The inlet pressure is 10bar, and outlet pressure is vacuum.

Inside of the nozzle is initially vacuum too.

Imagine the air is dumped into the vacuum tank from the pressure tank through the supersonic nozzle.

*The problem is that if i initialize the values to vacuum condition(P=0 or 1), solution diverges.

When i set the initial condition same as the inlet values(P=10bar), solution works.

But i want the processes of reaching steadystate from vacuum condtition, i want to set the initial values as vacuum.

I used Openfoam and Fluent, both of them diverged in vacuum initial condition. (SonicFoam, rhocentralFoam in openfoam)

If you have any ideas, please let me know.

Thanks for your attention. Iím looking forward to your reply.


https://www.cfd-online.com/Forums/at...1&d=1564339334

pete20r2 July 30, 2019 12:46

Yes, initial conditions can and often do cause divergence (usually in combination with other factors).
I would not set outlet to 0 as this is non-physical (bad start if you want real physics out).
Which pressure are you setting?
Can you share your solver output, it usually tells you what went wrong.

Carlo_P July 31, 2019 03:56

Hey Kim,
maybe the problem con be the relaxation factors.
In a case so difficult, you can have maybe some numerical issues.
At which iteration would diverge?
Can you try with very small relaxation factor?


Are you initializing all the domain with p=0?
This can be a problem, since not all the domain in the result is p=0.
If you want, you can try to initialize in different way, using TopoSetDict .
Cheers,
Carlo



p.s. can you post your BC files? I'm really interested in it.


All times are GMT -4. The time now is 16:49.