[Lena_K]

Hello!

Sorry for disturbing, i am not familiar with OpenFOAM, but currently i need to simulate the density distribution inside the gas target (i'm using sonicFoam because of high speed and kOmegaSST model because of high pressure gradient). I've wrote geometry in blockMeshDict, and it seems to be ok, and there's no error when i'm running the simulation. But the problem is that final density distribution (after there's almost no changes in distributions at different timesteps) does not looks reasonable.
I've tried different pressure for the inlets (i've changed time step respectively), and simulation "stuck" at different stage. I would be very grateful for any help.

Regards,
Lena

[Lena_K]

Looks like problem is lying under proper mesh. I'm getting very high courant number for big timesteps, therefore i reduce the timestep up to 10^-11 s and i got reasonable max courant number, but very small mean courant number. Could that cause a problem?

[Lena_K]

I was running a 3D simulation with timestep of 10^-9 s at fine mesh, and it does not reach steady state, therefore i build 2D model and test it with the same timestep. Simulation has been running for two weeks with more progress, than 3D simulation, but the fluctuations in the resulting fields are very big, i guess due to high pressure (for low pressures and therefore bigger timesteps it works good). Therefore i have a problem: to avoid the fluctuation i need do increase time step, this lead to higher Courant number -> simulation crash. What can i do to solve this problem?