Entropy generation

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Yvan Fournier
Posts: 4208
Joined: Mon Feb 20, 2012 3:25 pm

Re: Entropy generation

Post by Yvan Fournier »

Hello,

There seems to be an extra "iccocg" and a missing "ilved" argument to grdvec in your example.
Check (grep) for grdvec in other files for the correct syntax.

In addition, you setup seems to be based on 3.0-rc1, rather than a 3.0.0 or more current version. I recommend updating, as release candidate versions are not supported. We recommend using the latest release of a series, to avoid known and solved bugs, but if an issue seems to have no relation to some fix mentioned in the NEWS file, we usually try to help with slightly older releases. On the other hand, we do not bother with -rc or -beta versions once ".0" versions are released.

Best regards,

Yvan
Hakim
Posts: 37
Joined: Tue Sep 18, 2012 9:05 am

Re: Entropy generation

Post by Hakim »

hello,

Can someone explain me why I get an evolution of the variable on surface only?

best regards,

hakim.
Attachments
processus.tar.gz
(525.64 KiB) Downloaded 444 times
Yvan Fournier
Posts: 4208
Joined: Mon Feb 20, 2012 3:25 pm

Re: Entropy generation

Post by Yvan Fournier »

Hello,

Are you sure this is the case you ran ? One of your user subroutines does not compile :
<....>/CASE1/SRC/cs_user_postprocess_var.f90:284.2:

visct = max (visct,rom*sqdu)
1
Error: Incompatible ranks 0 and 1 in assignment at (1)
(you need sqdu(iel) here).

I'll let you check this first.

Regards,

Yvan
Hakim
Posts: 37
Joined: Tue Sep 18, 2012 9:05 am

Re: Entropy generation

Post by Hakim »

Hello,

even with the correction I still get the same problem, namely the evolution at the surface only,

best regards,

Hakim.
Hakim
Posts: 37
Joined: Tue Sep 18, 2012 9:05 am

Re: Entropy generation

Post by Hakim »

Hello Ivan,

I'd like to extract the maximum values of ​​SFBR (ILOC) calculated in cs_postprocess_var to a datafile.dat and that for each time step.
I tried with the following program but the problem is that I still do not get the same values ​​if I redo the simulation, of course after convergence. but it works for the temperature, I always get the same temperature values.

impout = impusr(1)
open(impout,file='profile.dat')
write(impout,99) max(rtp(iel,ivar),1.d-10),max(sfbr(nfbrps),1.d-10), ntcabs
99 format(6g17.9)

Best regards.

Hakim.
Hakim
Posts: 37
Joined: Tue Sep 18, 2012 9:05 am

Re: Entropy generation

Post by Hakim »

the same thing hapens with the following:
impout = impusr(1)
open(impout,file='profile.dat')
write(impout,99) max(rtp(iel,ivar),1.d-10),max(sfbr(iloc),1.d-10), ntcabs
99 format(6g17.9)
Yvan Fournier
Posts: 4208
Joined: Mon Feb 20, 2012 3:25 pm

Re: Entropy generation

Post by Yvan Fournier »

Hello,

Are you running in parallel or serial ? Your program includes no precautions for parallellism.

In addition, you're not looping on "iel", so it is either not initialized, or uses the last loop variable (probably past-the end).

Regards,

Yvan
Post Reply