Search found 239 matches

by Luciano Garelli
Wed Aug 19, 2020 4:25 am
Forum: code_saturne usage
Topic: Different convective schemes for different zones
Replies: 8
Views: 395

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: 395

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: 409

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, I have used something similar ...
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: 409

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: 1038

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: 750

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: 690

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_tt(iel) endif S...
by Luciano Garelli
Tue Mar 17, 2020 1:46 pm
Forum: code_saturne usage
Topic: Breaking the GUI
Replies: 3
Views: 611

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: 690

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 ...
by Luciano Garelli
Fri Dec 06, 2019 4:47 pm
Forum: code_saturne usage
Topic: How to get the coordinates of specific cells?
Replies: 4
Views: 557

Re: How to get the coordinates of specific cells?

Hello, In the Doxygen documentation that I posted previously is the data about your question. If you use ncel (cs_glob_mesh->n_cells) you will get the local number of cells in case of parallelism. If you use ncelet (cs_glob_mesh->n_cells_with_ghosts) you will get the local number of cells plus the g...