Search found 280 matches

by Luciano Garelli
Sat Oct 03, 2020 3:26 am
Forum: code_saturne usage
Topic: LES with conjugate heat transfer
Replies: 7
Views: 3316

Re: LES with conjugate heat transfer

Hello,

I did a test in a small case using CS 6.0.4 and the temperature field in celsius. I used LES dyn model restarting from K-omega SST and I don't get this behaviour in the temperature iterations. Maybe it was a bug that was fixed in the subsequent versions.


Regards,
Luciano
run_solver.log
(534.85 KiB) Downloaded 171 times
by Luciano Garelli
Wed Aug 19, 2020 4:25 am
Forum: code_saturne usage
Topic: Different convective schemes for different zones
Replies: 8
Views: 4926

Re: Different convective schemes for different zones

Hello Ok, but setting different blencv in your domain you are mixing first and second order schemes. May be the ischcv = 3 ( Blending between Second Order Linear Upwind and Centered scheme) is what you need, because your will be using two second order schemes and the unphysical oscillations can be r...
by Luciano Garelli
Tue Aug 18, 2020 4:11 pm
Forum: code_saturne usage
Topic: Different convective schemes for different zones
Replies: 8
Views: 4926

Re: Different convective schemes for different zones

Hello,

Did you try with another convective scheme like SOLU (ischcv). Also, you can set different limiters (isstpc)

Regards,
Luciano
by Luciano Garelli
Fri Jun 12, 2020 2:33 am
Forum: code_saturne usage
Topic: Store surface normal vector to array in parallel simulation
Replies: 4
Views: 3097

Re: Store surface normal vector to array in parallel simulation

Hello, I think that before call cs_parall_allgather_r() you have to do a parallel sum (parcpt) in order to sum "nlelt1" to get the total "n_g_elts". After you get "n_g_elts" you can call to cs_parall_allgather_r() Let me know if I'm right, because I never use it in CS, ...
by Luciano Garelli
Mon Jun 08, 2020 9:54 pm
Forum: code_saturne usage
Topic: Store surface normal vector to array in parallel simulation
Replies: 4
Views: 3097

Re: Store surface normal vector to array in parallel simulation

Hello, After partitioning the mesh each process will have a portion or not of the boundary WALL1, so each process has to create a local array and then you need to join all these array with a mpi gather. You can take a look to the file cs_user_extra_operations-parallel_operations.f90 in the example d...
by Luciano Garelli
Sun May 17, 2020 11:32 pm
Forum: code_saturne usage
Topic: Using fans definitions
Replies: 1
Views: 3652

Re: Using fans definitions

Hello,

You can check the Doxygen documentation

https://www.code-saturne.org/cms/sites/ ... an_8c.html

and also this old post

viewtopic.php?f=2&t=1831&start=10

Regards,
Luciano
by Luciano Garelli
Thu Mar 19, 2020 2:03 pm
Forum: code_saturne usage
Topic: ibeta
Replies: 2
Views: 2652

Re: ibeta

Hello, ibeta it is a (integer) pointer for thermal expansion coefficient. To access to the thermal expansion coeff you can use this example used inphyvar.f90, I think. double precision, dimension(:), pointer :: cpro_beta .... .... if (ibeta.ge.0) then iok1 = 0 call field_get_val_s(ibeta, cpro_beta) ...
by Luciano Garelli
Tue Mar 17, 2020 3:01 pm
Forum: code_saturne usage
Topic: thermal expansion coefficient
Replies: 4
Views: 3337

Re: thermal expansion coefficient

Hello, In the file /turb/divrit.f90 are the same lines that you mention, ! AFM model ! "-C_theta*k/eps*( xi* uT'.Grad u + eta*beta*g_i*T'^2)" if (iturt(iscal).eq.20) then if (itt.gt.0.and.ibeta.ge.0) then temp(ii) = temp(ii) - ctheta(iscal)*xtt* & etaafm*cpro_beta(iel)*grav(ii)*cvara_t...
by Luciano Garelli
Tue Mar 17, 2020 1:46 pm
Forum: code_saturne usage
Topic: Breaking the GUI
Replies: 3
Views: 2987

Re: Breaking the GUI

Hello,

I download your xml and can't reproduce you error. I make some changes and save it without any error or warning.

What version of CS are you using?

Regards,
Luciano
by Luciano Garelli
Tue Mar 17, 2020 1:33 pm
Forum: code_saturne usage
Topic: thermal expansion coefficient
Replies: 4
Views: 3337

Re: thermal expansion coefficient

Hello, First of all you have to activate the thermal model, then in the fluid property tab change density from constant to user law. In the mathematical expression editor you can write an equation for the density which include the thermal expansion coefficient, like the Boussinesq approximation for ...