code_saturne User's Forum

hello;
it seems that the mesh is good and the rotor is so far from the outlet, i think that the problem occurs at the interface boundary, calculations don't start if the face joining is not added in the GUI.
regards.

hello, thank you, now it work. but during rotation i have some problems: incoming flow at outlet and the mesh is destructed at the boundary rotor/stator (please see attached file). for that calculations stop by casting this message in the listing file: "Jacobi: error (divergence) solving for ve...

hello,
please find attached the "cs_user_turbomachinery.c" subroutine.
i update to the 3.2.1 version, and now i have a new message error:
Signal SIGSEGV (accès à une zone mémoire interdite) intercepté !

best regards.

hello, thank you for the clarification; now i have an error message: **************************************************************************** cs_solver: cs_turbomachinery.c :619 : _select_rotor_cells: l'assertion « tbm->rotor_cells_c != ((void *)0) » a échoué. cs_solver: cs_turbomachinery.c :619...

hello, i have a problem with cell_criteria definition in cs_user_turbomachinery.c: i didn't understand what means the first raw "cylinder[0.0, 0.0, 0.0," for the second & the third i think that corresponds to rotation axis and velocity. note: my rotor is a cylindrical domain with a rad...

thank's it's always the same thing even without the if (. true.) return I send you my subroutine and file compil.log when installing my cs3, it given an error in cgns instal and it proposed to installed - without-cgns!!! is that it is causing problem ??? thank you very very much compile.log cs_user_...

hello thank you for your response for term source I used a constant source term S=AU+B A=0 B=cst crvimp(iel) = 0.d0 crvexp(1,iel)=cst if (ivar.eq.iu .or. ivar.eq.iv .or. ivar.eq.iw) then qdm = 100.d0 call getcel('rot', nlelt, lstelt) do iel = 1, nlet crvimp(1,1,iel) = 0.0d0 enddo do iel = 1, nlet cr...

Hello I ask you for help because I have problems with the new version of saturne CS3 in the applications of ustsns in cs-usersouceterm and the user-boundary. I have adapted the subroutine to this new version, they work in first sight! but after (1000 iteration)get back with my results but without th...

hello,
you can try to give an hight stiffens to the two direction that you want to fix, and eliminate forces applied in these directions. (ex: for a needed motion in y direction: fluid_fx=fluid_fz=0 and kx,kz>>ky)
best regards.

Hello ;
thank you luca ,
do you mean modifying the routine called (usstr2) ?
And what about defining a rigid buddy around the structure ? I usually use a hight mesh viscosity but for large displacement I had negative volume.
good luck.