Compute Stresses at Boundary (condli.f90)

[vshinde]

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

[Yvan Fournier]

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