Hello,
I am, too, facing an issue regarding an imposed concentration and solid discharges at boundaries.
I would like to impose a concentration at one of the liquid boundary. First I tried to use the key word "CONCENTRATIONS PAR CLASSES AUX FRONTIERES" and the graphic printout variable CS1. No concentration can be seen in the domain.
Then I tried to implement a fortran subroutine 'conlit.f',changing features in the right loop =>> CBOR%ADR(1)%P%R(K)=0.001 for instance
But it does not work either.
Regarding the imposed solid discharge at the boundary, I used the same subroutine and added QBOR%ADR(1)%P%R(K)=0.001 where K is the boundary points
It works only for the piece of mesh at the very near of the boundary line. Solid discharge does not spread into the domain.
I managed to modify the fourth column in the BOUNDARY CONDITION FILE, which is linked to Q2BOR:
4 5 5 1E-4 0 0 0 4 0 0 0
LIHBOR LIQBOR Q2BOR LIEBOR
Results are unfortunately the same.
Attached are the sisyphe and telemac steering file. I use the V6P1 version. There is also the fortran file.
I would be very glad if you could help me.
Best regards
Nicolas LORRAIN
ENSEEIHT TOULOUSE