code_saturne User's Forum

Hello,

I'm trying to use the radiative heat transfer module in version 3.2.4 together with some faces with itypfb(ifac) = iindef.

Checking raycll.f90, there is no standard treatment of such faces. Is there anything somewhere else? Is this a bug? Or am I supposed to edit some user subroutine with such treatment? I'm already using usray2.f90 and I define a group also for those faces, but the code crashes anyway.


Thanks

Hello,

As you suggested there is no treatment of such faces in the boundary conditions treatment for radiative heat transfer. I don't exactly know why, but I suspect it is linked to the fact that there is no icodcl, rcodcl arrays for the radiative intensity.

For a standard variable, when you set an undefined type, the code then directly transposes the information given in icodcl (Dirichlet, Neumann to make it short), rcodcl (Dirichlet value, Neumann flux value).

Why do you need to have boundary faces with undefined type ?
Would it be possible to prescribe one of the existing type and then adapt the boundary conditions values in usray2 ?

Regards,
Erwan.