Can initial condition makes the solution diverge?

kim4545 · July 28, 2019, 16:02

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.

pete20r2 · July 30, 2019, 13: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, 04: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.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
decompose dependent solution	arionfard	OpenFOAM	3	December 10, 2018 10:36
Foam::error::printStack after large bounding epsilon events	StoneCFD727	OpenFOAM Running, Solving & CFD	1	January 1, 2018 18:04
Free surface issues with interDyMFoam for hydroturbine	oumnion	OpenFOAM Running, Solving & CFD	0	October 6, 2017 15:05
pressure in incompressible solvers e.g. simpleFoam	chrizzl	OpenFOAM Running, Solving & CFD	13	March 28, 2017 06:49
pimpleFoam: turbulence->correct(); is not executed when using residualControl	hfs	OpenFOAM Running, Solving & CFD	3	October 29, 2013 09:35

July 30, 2019, 13:46		#2
pete20r2 Senior Member Peter Baskovich Join Date: Jul 2014 Posts: 127 Rep Power: 12	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.

July 31, 2019, 04:56		#3
Carlo_P Senior Member Carlo_P Join Date: May 2019 Location: Italy Posts: 176 Rep Power: 8	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.