12.07.2015 Views

Tamtam Proceedings - lamsin

Tamtam Proceedings - lamsin

Tamtam Proceedings - lamsin

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

368 Younes et al.La matrice B étant inversible, l’équation (4) s’écrit3A 1A Aiijjj1Q a B P TP(5)Avec (voir [9])rjkrjk rjkrki rjkrij 1 1 111 1B rjkrki rkirki rij rki 1 1 1A 3 rjk rij rij rki rij r ij 1 1 1 (6) étant le coefficient défini par [9],3ij ik jk Bij 48 A312(7)j1avec rijle vecteur partant du nœud i au nœud j , ijson module ( ij rij).L’équation (1) est discrétisée avec un schéma implicite en temps :3sA1 PA,nA,n A A P Qj A fA Qs(8)tj1En utilisant (5), on obtient3AA,n1 aA,n1 A,n QsP iTPi P (9) i1 3 3sA1 1 1 , i B ij , i3 , a .tavec L’équation (5) devient j1 j13 3AA,n 1 ai A,n 1 1 A,n 1 i A,n iQ sQi ajTPj Bij TPj P j1 j1 ALe système final est obtenu en écrivant la continuité des pressions TPi TPi (10)des flux entre deux éléments adjacents A et B :Q Q 0(11)A,n1 B,n1iiPour une arête i , le terme diagonal de la matrice lui correspondant est :2A BA1aimii mii miiavec miiaBii (12) BetTAMTAM –Tunis– 2005

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!