Error run LES when restart from RANS

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Ruonan
Posts: 136
Joined: Mon Dec 14, 2020 11:38 am

Error run LES when restart from RANS

Post by Ruonan »

Dear developers,

Please could you help me with this problem:

I'm trying to run an LES solution, restart from a converged RANS solution, to make the LES solution converge quickly. But one error occurred: "Invalid (not-a-number) values detected for a field".
error.png
I checked the log file, guess maybe it's because of an error when restarting LES from the previous RANS result. I have no problem when restarting an LES solution from a previous LES result, or restarting RANS from RANS.

Another thing need to mention is that, the previous RANS mesh is coarse, and the current LES mesh is refined. But the CFD model is the same. So I ticked the "Different mesh" box, and choose the .csm file from the previous RANS checkpoint folder. Also, in "Advanced options", I didn't tick the "Read the auxiliary restart file" box.
Start Restart.PNG
Could you please suggest to me, if the setting is correct or not? And how to deal with this error? Thank you so much for your help!

Best regards,
Ruonan
Ruonan
Posts: 136
Joined: Mon Dec 14, 2020 11:38 am

Re: Error run LES when restart from RANS

Post by Ruonan »

Please find the xml and log files here. Thank you!

Best regards,
Ruonan
Attachments
setup.xml
(16.75 KiB) Downloaded 110 times
run_solver.log
(30.68 KiB) Downloaded 111 times
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: Error run LES when restart from RANS

Post by Yvan Fournier »

Hello,

I received another bug report and test case for a similar issue, probably the same one as here. I located the error (exponential in RANS to LES term, with negative R_00 value in some cells; in src/turb/cs_les_synthetic_eddy_method.c, line 1670 for current master branch), and an checking with colleagues which formula should be used there instead of the current one. So I should have a fix very soon.

Best regards,

Yvan
Ruonan
Posts: 136
Joined: Mon Dec 14, 2020 11:38 am

Re: Error run LES when restart from RANS

Post by Ruonan »

Hello Yvan,

Thank you very much! I will follow your updates on GitHub.

Best regards,
Ruonan
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: Error run LES when restart from RANS

Post by Yvan Fournier »

Hello,

This has been fixed on the master and v7.0 branches a few days ago.

Best regards,

Yvan
Gen
Posts: 17
Joined: Sat Sep 05, 2020 7:39 am

Re: Error run LES when restart from RANS

Post by Gen »

Hello Yvan,

I got the fixed version from Git but I still received the same error as Ruonan mentioned before, when I tried to restart LES from RANS.
src/base/cs_log_iteration.c:842: Fatal error.

Invalid (not-a-number) values detected for a field.
Did I make a mistake? Looking forward to your reply, thank you so much!

Best regards,
Gen
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: Error run LES when restart from RANS

Post by Yvan Fournier »

Hello,

There might be another bug. Do you have a small test case on which I could try to reproduce this ?
If your cas is large, the other option would be to run it on your side with a debug build, which would detect the cause of the "NaN" earlier and help locate the issue.

Best regards,

Yvan
Gen
Posts: 17
Joined: Sat Sep 05, 2020 7:39 am

Re: Error run LES when restart from RANS

Post by Gen »

Hello Yvan,

Please kindly find the setup file and the mesh of my test case in the attachment. It's a simple channel flow, I first tested it with k-omega model than I tried to run Smagorinsky model based on the previous checkpoint result.
In the meanwhile I will try to fix it by my side with the debug mode.
Thank you so much for your help!Bonne journée!

Best regards,
Gen
Attachments
Mesh_590_2.zip
(2.48 MiB) Downloaded 105 times
setup.xml
(12.16 KiB) Downloaded 101 times
Gen
Posts: 17
Joined: Sat Sep 05, 2020 7:39 am

Re: Error run LES when restart from RANS

Post by Gen »

As for debug build, shall I perform a separate installation? Could I just modify setup file like this when I install saturne.py?
# Install Code_Saturne with debugging symbols
#--------------------------------------------------------
debug yes
And if I use other code_saturne versions, I will always get the following error:
READING THE MAIN RESTART FILE

Start reading

Warning: the number of variables or properties has been changed
relative to the restart file.

current: 4 model and 0 user variables
17 model and 0 user properties
previous: 5 model and 0 user variables
12 model and 0 user properties

The computation continues, with a partial restart.
Reading dimensions complete
Reading the previous time step number (restarting computation) NTPABS = 80000
Reading the previous time step number (restarting computation) TTPABS = 0.8000E+01
Reading options complete
restart/main.csc: section "wall_yplus::vals::0" not present.

Warning: the turbulence model has been changed
relative to the restart file.

current model: 40
previous model: 60

The computation continues, with a partial an/or adapted restart.
restart/main.csc: section "epsilon::vals::0" not present.

SIGFPE signal (floating point exception) intercepted!
I will be really thankful if you could give me some help. Merci beaucoup!

Best regards,
Gen
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: Error run LES when restart from RANS

Post by Yvan Fournier »

Hello,

Yes, to install a debug version, you can change the setting in the installer.

Debug builds are usually about 3x slower than regular builds, so I personally prefer to do a second install an another (adjacent) path, so as to always be able to access both versions.

In any case, the SIGFPE you observe and the "Not a Number" mentioned in a previous post are probably related.
Running your case under a debugger, it seems "epsilon" is zero, and e divide by that.

So I assume the omega->epsilon conversion is missing here, and the restart would probably work from a k-epsilon computation. I'll check if I can fix this myself for k-omega or see with colleagues more exert in this part of the code.

Best regards,

Yvan

Edit: I fixed to bug on the master and v7.0 branches. The GitHub mirror has been synchronized.
Post Reply