code_saturne User's Forum

Hello Yvan, Thank you very much for debugging. Unfortunately, this error occurred again. I use the new feature version 7.1. I found when I restart LES Smagorinsky from k-omega SST, it's OK. But when I restart LES Smagorinsky from Spalart-Allmaras, this error happened again. So there may be something...

Thank you Yvan! This is new knowledge for me.

Best regards,
Ruonan

Hello Yvan, Thank you very much for pointing out the mesh problem. At first I thought it was a warning, not an error, so I ignored it. Now I solved it. The original mesh was structured with multiple zones. I used Pointwise to join all the zones into one unstructured zone, then exported it. Then the ...

Hello Yvan, Thank you very much. Following your kind comments, I did some experiments: 1. First, I removed the source term subroutine, but ended up with the same error: divergence pressure. 2. Then, I decreased the timestep by a factor of 1000, from 1.9E-4s to 1.9E-7s, but again ended up with the sa...

Dear Code_Saturne Team, Could you please help me with this error? I am running a new case using the new feature version 7.1. But one error occurred at the first time step: Flexible Conjugate Gradient: error (divergence) solving for Pressure Conjugate Gradient [Pressure:coarse:10]: divergence after 1...

Hello Yvan, Thanks a lot for your suggestions. But sorry that I still can't get it work. I tried the first method you suggested. The compilation was ok, but there seems no change. The case still can't restart, with the same error (infinite density) as before. I then tried your second method, adding ...

Hello Yvan, Thank you very much, your help is much appreciated. Unfortunately, I tried both "read auxiliary restart file" and "standard I/O serial", but the error was still there. I still can't restart. But interestingly, when I remove the cs_user_physical_properties.f90, the cas...

Hello Yvan, Thank you very much. Yes, if I run it for only a few iterations and then restart, I won't have this issue either. Regarding this case, I ran 366390 iterations before I restart. I used "n_procs: 56, n_threads: 1". The run_solver.log for previous computation is too large to uploa...

Hello Yvan,

Attached is the second subroutine. Thank you!

Best regards,
Ruonan

Hello Yvan, Thank you very much for debugging! I attached the mesh, xml and subroutines. I use two subroutines, one for changing the density, the other one for adding centrifugal source term explicitly to the momentum equation. One strange thing is, when I run it for a shorter time and then restart,...