Search found 243 matches

by Luciano Garelli
Fri Oct 30, 2020 2:08 am
Forum: code_saturne usage
Topic: Internal coupling spurious velocity
Replies: 7
Views: 536

Re: Internal coupling spurious velocity

Hello Yvan, Thanks for the suggestion, I will try it. It is like to set a penalization for this term In the solid region. Another workaround that I found is to set the same density in both, fluid and solid, and in order to maintain the volumetric heat capacity in the solid part I have set the heat c...
by Luciano Garelli
Tue Oct 20, 2020 11:57 pm
Forum: code_saturne usage
Topic: Internal coupling spurious velocity
Replies: 7
Views: 536

Re: Internal coupling spurious velocity

Hello Yvan, Thanks for your quick reply. I have tried your suggestion but I get the same behavior with some iterations in the velocity field in the solid part. In the setup.log appears the text "Fluid-solid mode (disable dynamics in the solid part)" Regards, Luciano setup.log run_solver.log
by Luciano Garelli
Tue Oct 20, 2020 10:38 pm
Forum: code_saturne usage
Topic: Internal coupling spurious velocity
Replies: 7
Views: 536

Internal coupling spurious velocity

Hello, I'm running in CS 6.0.4 a simple example using the internal coupling feature, but when setting different densities (rho_air=1, rho_solid=100) spurious velocities appears at the solid region as seen in the picture. The densities and temperature are constant and the gravity is turned on. If gra...
by Luciano Garelli
Sat Oct 03, 2020 3:26 am
Forum: code_saturne usage
Topic: LES with conjugate heat transfer
Replies: 7
Views: 482

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 13 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: 922

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

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

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

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

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

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) ...