code_saturne User's Forum

Posted: **Tue Oct 17, 2023 1:13 pm**

Hello,

I simple question, can I use:

call getfbr('INLET', nlelt, lstelt)

do ilelt = 1, nlelt

  ifac = lstelt(ilelt)

  rcodcl(ifac,iu,1) = 0.5d0
  rcodcl(ifac,iv,1) = 0.d0
  rcodcl(ifac,iw,1) = 0.d0
  
  rcodcl(ifac,ik,1) = 0.5d0
  rcodcl(ifac,iomg,1) = 10000.d0

enddo

(values prescribed just as an example) inside cs_user_boundary.f90 to impose the desired profiles of k and omega at inlet?
When I do that, the code seems to ignore the subroutine and takes the values specified in the GUI (I've set "Calculation by formula" for the Turbulence).

Thank you very much in advance.
Kind regards,
Daniele

Posted: **Wed Oct 18, 2023 4:19 pm**

Hello,

Which version of the code are you using, and is this with a specific physical model ?
Is the velocity itself taken into account ?

Regards,

Yvan

Posted: **Wed Oct 25, 2023 9:52 am**

Hello Yvan,

After double checking, the routine actually works correctly, sorry for the useless loss of time.

Thank you for your help.
Best regards,
Daniele

code_saturne User's Forum

prescribing k and omega user profiles at inlet

prescribing k and omega user profiles at inlet

Re: prescribing k and omega user profiles at inlet

Re: prescribing k and omega user profiles at inlet