Télécharger enc2D-therco.dgibi
* fichier : enc2D-therco.dgibi ************************************************************************ ************************************************************************ * * Enceinte 2D Axisymetrique * Refroidissement d'un gaz chaud (AIR) * Interaction GAz/paroi/temperature exterieure * Teste thermique paroi et coefficient d'echange COMPLET= VRAI ; COMPLET= FAUX ; GRAPH = VRAI ; GRAPH = FAUX ; 'SI' COMPLET ; nbit=200 ; DT0 = 2. ; n1=3 ; n2=5; n3=10; ds1= 0.02 ; ds2= 0.2 ; 'SINON' ; nbit=20 ; DT0 = 2. ; n1=3 ; n2=5; n3=10; ds1= 0.02 ; ds2= 0.2 ; 'FINSI' ; p1=5.e-4 0 ; p2=0.5 0. ; p3=1 0 ; p4=1 2; p41=0.5 2. ; p5=5.e-4 2 ; mt= 'DALLER' bas paroid plaf axe ; Si GRAPH; trace (mt et vtp); Finsi ; rxt = 'TABLE' ; rxt.'VERSION'= 'V0' ; rxt.'vtf'=mt ; rxt.'axe' = axe ; rxt.'TF0' = 100. ; rxt.'PT0' = 1.0e5 ; rxt.'THERMP'=VRAI ; rxt.'THERCO'=VRAI ; rxt.'vtp'=vtp; rxt.'ROCP'=1.e4 ; rxt.'LAMBDA'=40. ; rxt.'Tp0' = 30. ; rxt.'ECHAN'=50. ; rxt.'ECHEXT' = VRAI ; rxt.'HEXT'=50. ; rxt.'TPEXT'=20. ; rxt.'parext'=parext ; rxt.'pi' = 0.5 0.5 ; rxt.'DT0' = DT0 ; rxt.'RENU' = 'RIEN' ; rxt.'DISCR'= 'QUAF'; rxt.'KPRE' = 'MSOMMET'; rxt.'MODTURB' = 'NUTURB' ; rxt.'NUT' = 0.005 ; rxt.'GRAPH' = GRAPH ; $vtp=rxt.'GEO'.'$vtp' ; $vtf=rxt.'GEO'.'$vtf' ; *rxt.'LAMBDA'=kcht $vtf scal sommet 40. ; EXECRXT nbit rxt ; un=rxt.'TIC'.'UN'; tf =rxt.'TIC'.'TF' ; $vtf=rxt.'GEO'.'$vtf' ; $vtp=rxt.'GEO'.'$vtp' ; rho=rxt.'TIC'.'RHO' ; un =rxt.'TIC'.'UN'; tp =rxt.'TIC'.'TP'; Si (NON COMPLET) ; ltfm=Prog 100.00 79.731 76.588 71.335 67.287 63.939 60.724 57.698 54.980 52.548 50.370 48.421 46.657 45.037 43.536 42.140 40.845 39.644 38.530 37.493 36.528 ; lPT =Prog 1.00000E+05 98275. 92267. 91693. 91538. 90412. 89390. 88653. 87949. 87272. 86679. 86162. 85693. 85262. 84864. 84492. 84145. 83823. 83523. 83246. 82987. ; Lmaxu=Prog 0.0000 0.0000 0.26029 0.44149 0.51786 0.48432 0.42517 0.35634 0.29129 0.23483 0.19800 0.17967 0.16889 0.16210 0.15651 0.15082 0.14500 0.13929 0.13407 0.12947 0.12546 ; tic=rxt.'TIC' ; Si (ERtf '>' 2.e-4) ; erreur 5 ; Finsi ; Si (ERPT '>' 2.e-4) ; erreur 5 ; Finsi ; Si (ERum '>' 2.e-3) ; erreur 5 ; Finsi ; list evuz ; lrr= prog 8.68197E-39 1.81018E-02 5.70173E-02 8.89906E-02 0.10936 0.12289 0.12546 0.12244 0.10940 9.93363E-02 8.06895E-02 7.14766E-02 5.53056E-02 5.39819E-02 4.18936E-02 4.67526E-02 3.46966E-02 3.43026E-02 2.18271E-02 3.34970E-03 -9.75217E-38 ; mess ' Ecart ' er ; Si ( er > 2.e-4) ; erreur 5 ; finsi ; Finsi ; Si GRAPH ; dess evauz 'TITRE' 'Vitesse sur l axe'; dess evarh 'TITRE' 'Densite melange sur l axe'; dess evatf 'TITRE' 'Temperature sur l axe'; trace ung vtf; trace tp vtp 'TITRE' ' Tp ' ; trace tf vtf 'TITRE' ' Tf ' ; trace rho vtf 'TITRE' ' Rho' ; Finsi ; FIN ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales