prescribing k and omega user profiles at inlet

[daniele]

Hello,

I simple question, can I use:

Code: Select all

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

[Yvan Fournier]

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

[daniele]

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