CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Siemens (https://www.cfd-online.com/Forums/siemens/)
-   -   Call of Ratusr.f (https://www.cfd-online.com/Forums/siemens/54927-call-ratusr-f.html)

jojo October 12, 2005 13:36

Call of Ratusr.f
 
Hi everyone,

Although I have set C148 to -1, Ratusr.f is not called by star. What must I do to make it work, please?

Thanks a lot for your help.

JoJo

jojo October 12, 2005 13:37

Re: Call of Ratusr.f *NM*
 

elvis October 14, 2005 02:07

Re: Call of Ratusr.f
 
i don't know

Jamie October 14, 2005 08:48

Re: Call of Ratusr.f
 
Which version of the code are you using.

jojo October 14, 2005 09:31

Re: Call of Ratusr.f + sigsegv signal
 
Hello everybody.

Thanks to Jamie (--> I use 324) and Elvis for your message.

I have found the answer. It seems that once Ratusr.f has been written in the ufile directory, one has to enforce Prostar to save the chemical scheme again (by erasing and rewriting something in it). By doing this, Star will 'see' Ratusr.f at the next simulation and will link it.

I have another question: sometime, my simulation crashes with the message on the main window:

============================= [ 1 ] ***Received signal sigsegv - exiting PNP: Shutdown [2005-10-13-13:50:00] Execution aborted by request (SIGABRT) after 398 seconds (TOTAL ELAPSED TIME). =============================

I find nothing else in the .run and .info file. What does happen?


Jamie October 17, 2005 13:37

Re: Call of Ratusr.f + sigsegv signal
 
Not sure what the problem is. I'd double-check the user coding. I once had a typo in a variable name. The user coding would compile without an error but then the run would seg fault. Once I found the typo it ran fine.

jojo October 17, 2005 14:33

Re: Call of Ratusr.f + sigsegv signal
 
Thanks for your message.

The problem I have seems to be related to convergence with the VODE solver. I have solved it by imposing mass fraction exponents > 1 on every specie to relax the strength of the reaction (the pre-exponential factor must be kept to 1 as it enforces equilibrium).


All times are GMT -4. The time now is 11:03.