2. 1 -  Définir les CL de flux

Définir les CL de flux dans le programme principal :   

On procède comme pour les conditions de Dirichlet en décrivant les conditions de flux par un tableau clne, de dimension (nclne,2) :

  • la 1ère colonne de clne contient les numéros des noeuds d’extrémités où le flux est imposé,

  • la 2ème colonne de clne contient les valeurs imposées,

  • le nombre de ligne est le nombre des extrémités où le flux est imposé

Ici :
une seule extrémité est concernée : $nclne=1$ et $clne(1,2)=q_ o$ est la valeur imposée à l’extrémité de numéro $clne(1,1)=nx$

sous scilab  : qo =1 ;  clne =  [  nx , qo ] ;

remarque : ne pas oublier d’adapter les CL de Dirichlet.