courant number explosion on good mesh not shown in visual output

geth03 · October 15, 2021, 04:11

hi all,

i am running a simulation with reactingMultiphaseEulerFoam on OF version 7.

my Co keeps increasing from a certain time on, i could also say that my timestep is decreasing for a constant Co, i think the conclusion is the same. i did both, keeping Co constant and checking timestep and vice versa.

the strange part is that my mesh is really really good. also, when i try to visualize the high Co, the maximum value shown in paraview is really small, around 1, compared to what i see in the terminal output, which is 10, 20 or higher (depends on how far i let it go before stopping it, it increases with time).

i tried to modify the mesh, ran the case again and same thing happens.

does anybody have any idea why this happens?

regards.

geth03 · October 20, 2021, 10:54

the Co for multiphase cases is calculated for each phase and also for the total flux. the lib for Co takes only account the total flux, which is the reason that Co shown in the terminal does not match the ouput value from the lib.

October 15, 2021, 04:11	courant number explosion on good mesh not shown in visual output	#1
geth03 Senior Member Join Date: Dec 2019 Location: Cologne, Germany Posts: 369 Rep Power: 8	hi all, i am running a simulation with reactingMultiphaseEulerFoam on OF version 7. my Co keeps increasing from a certain time on, i could also say that my timestep is decreasing for a constant Co, i think the conclusion is the same. i did both, keeping Co constant and checking timestep and vice versa. the strange part is that my mesh is really really good. also, when i try to visualize the high Co, the maximum value shown in paraview is really small, around 1, compared to what i see in the terminal output, which is 10, 20 or higher (depends on how far i let it go before stopping it, it increases with time). i tried to modify the mesh, ran the case again and same thing happens. does anybody have any idea why this happens? regards.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[General] Extracting ParaView Data into Python Arrays	Jeffzda	ParaView	30	November 6, 2023 22:00
[Other] Can't Shake Erros: patch type 'patch' not constraint type 'empty'	BrendaEM	OpenFOAM Meshing & Mesh Conversion	12	April 3, 2022 19:32
[snappyHexMesh] No layers in a small gap	bobburnquist	OpenFOAM Meshing & Mesh Conversion	6	August 26, 2015 10:38
How to install CGNS under windows xp?	lzgwhy	Main CFD Forum	1	January 11, 2011 19:44
CGNS lib and Fortran compiler	manaliac	Main CFD Forum	2	November 29, 2010 07:25

October 20, 2021, 10:54		#2
geth03 Senior Member Join Date: Dec 2019 Location: Cologne, Germany Posts: 369 Rep Power: 8	the Co for multiphase cases is calculated for each phase and also for the total flux. the lib for Co takes only account the total flux, which is the reason that Co shown in the terminal does not match the ouput value from the lib.