- temperatures de paroi basses et hautes : 
t1=300
t2=900

- donnees initiales :
rho=1.3e-7
t=t1
u=0

- donnees physiques :
masse moleculaire : masse=0.663e-25
constante de Boltzmann : kboltz=1.3806503e-23
constante du gaz :   r_gaz=kboltz/masse
taux de relaxation:  tau=(2.117e-5*(t1/273.15))/(rho*r_gaz*t1)
(ou t1 et rho sont les valeurs donnees ci-dessus)

- temps de simulation : 
  -> temps necessaire pour atteindre un etat stationnaire
  -> 200 000 iterations avec dt=0.9*dt_maximum, ou dt_maximum est le pas de temps maximum assurant la positivite du schema
  -> necessite environ 16h de calcul en sequentiel sur un processeur a 2.4Ghz

- geometrie du domaine de calcul :
 
        
-----------------
|               |
|		|
|reservoir     	|    ______ 
|	    	|   / ____ \
|		|   ||    || tube de jonction
|		|   ||    ||
|		|   ||	  ||
|		|   ||	  ||
|------------  -|   ||    ||  	      	  
             ||     ||   -  --------------    	 
             ||     ||	 |                |   	 
             ||_____||	 |		  |      
             \_______/	 |		  |   	
                      	 |		  |   	
                         | reservoir	  | 	
		      	 |		  |   	
		      	 |		  |   	
		      	 |		  |   	
		      	 |----------------|   	
		      			      	



- geometrie du tube : il se compose de quatre parties separees ci-dessous 

    	    ______                     
    	   / ____ \                  	
    	   ||    ||    deux demi-cercles concentriques de centre (2,3) et de rayon 0.5 et 1.5
    	   ||    ||
    	   ||    ||
           
    	   
	   ||	 ||    quatre segments de longueur 3
    	   ||	 ||     
           ||	 ||
                 

    ||	   ||	       deux demi-cercles concentriques de centre (0,0) et de rayon 0.5 et 1.5
    ||     ||	 
    ||_____||    
    \_______/        
		      		            	


- temperature de paroi : 

 -> T=T1 sur les parois des reservoirs
 -> T(s) varie selon une loi affine le long de chacune des 4 parties
du tube, de T1 a T2 ou de T2 a T1, ou s est l'abscisse curviligne
associee a chaque paroi : 



        T1
-----------------
|               |
|		|
|	     	|    ______ 
|T1	    T1	|   / ____ \
|		|   ||    ||
|		| T1||    ||T2
|		|   ||	  ||
|		|   ||	  ||
|------------  -|   ||    || T1 	      	  
          T1 ||   T2||   -  --------------    	 
             ||     ||	 |                |   	 
             ||_____||	 |		  |      
             \_______/	 |		  |   	
                      	 |		  |   	
                      T1 |		  | T1	
		      	 |		  |   	
		      	 |		  |   	
		      	 |		  |   	
		      	 |----------------|   	
		      			      	
		      		T1            	
		      	
		      	
plus precisement, voici la repartition de la temperature pour chaque
partie du tube : 
    	    ______  
    	   / ____ \ 
    	   ||    || 
    	   ||    || 
    	 T1||    ||T2 
           	    
    	   	    
	 T1||	 ||T2 
    	   ||	 || 
         T2||	 ||T1 
                    
		    
  T1||	   ||T2	    
    ||     ||	    
    ||_____||       
    \_______/       
		    


- references des noeuds frontieres (pour le traitement des conditions
aux limites) : 

 -> 1 pour les noeuds des parois des reservoirs
 -> de 20 a 70 pour les parois du tube, comme indique ci-dessous :

              60
    	    ______  
    	   / ____ \ 
    	   ||    || 
    	   || 50 || 
    	   ||    || 
           	    
    	   	    
	   ||  70||70 
    	 40||40  || 
           ||	 || 
                    
		    
    ||	   ||	    
    ||  20 ||	    
    ||_____||       
    \_______/       

       30