Search found 245 matches

by Luciano Garelli
Thu Jan 21, 2021 3:36 am
Forum: code_saturne usage
Topic: Hagen–Poiseuille flow comparison between CFD and Theory
Replies: 3
Views: 20

Re: Hagen–Poiseuille flow comparison between CFD and Theory

Hello If you set an average/constant inlet velocity of 0.01m/s you will never have a maximum speed lower than this value (0.0074778) due to mass conservation, near the wall the velocity is lower than 0.01, so in the center has to be greater than 0,01. The theory for 2D flow between parallel flat pla...
by Luciano Garelli
Wed Dec 23, 2020 3:36 am
Forum: code_saturne usage
Topic: Pressure-driven flow for pipe
Replies: 6
Views: 218

Re: Pressure-driven flow for pipe

Hello, The line "call getfbr('2 and x < 0.01', nlelt, lstelt)" is used for the selection of boundary faces satisfying a specific criteria. The criteria is that belongs to a group of face number/color "2" and the coordinates of the face centers x-coordinate "x<0.01". In ...
by Luciano Garelli
Fri Oct 30, 2020 2:08 am
Forum: code_saturne usage
Topic: Internal coupling spurious velocity
Replies: 7
Views: 730

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

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

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

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 22 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: 1050

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

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

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

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