code_saturne User's Forum

Dear Code Saturne team,

I am using Code Saturne 3.0.
I use ALE (fluid structure interaction) module where forces on the structure faces are given by forbr(ndim, ifac). There seems to be a typo error in calculating this force in step 1 over 5 from 'condli.f90' file. Below is the piece of code from the same file. Should not the term in bold be 'velipb(ifac,3)' ? Please confirm.

------------------------------------------------------------------------------------------------
condli.f90
------------------------------------------------------------------------------------------------
! The implicit term is added after having updated the velocity
forbr(1,ifac) = ( cofafu(1,ifac) &
+ cofbfu(1,1,ifac) * velipb(ifac,1) &
+ cofbfu(1,2,ifac) * velipb(ifac,2) &
+ cofbfu(1,3,ifac) * velipb(ifac,3) )*srfbnf
forbr(2,ifac) = ( cofafu(2,ifac) &
+ cofbfu(2,1,ifac) * velipb(ifac,1) &
+ cofbfu(2,2,ifac) * velipb(ifac,2) &
+ cofbfu(2,3,ifac) * velipb(ifac,2))*srfbnf
forbr(3,ifac) = ( coefau(3,ifac) &
+ cofbfu(3,1,ifac) * velipb(ifac,1) &
+ cofbfu(3,2,ifac) * velipb(ifac,2) &
+ cofbfu(3,3,ifac) * velipb(ifac,3) )*srfbnf
------------------------------------------------------------------------------------------------

Thanks,
Vilas

Hello,

Yes, you are right. Thanks for reporting this. The fix has just been commited to version control, so it will be in the next patch releases (3.1.1 and 3.0.2).

Regards,

Yvan