1 : C HEA010 SOURCE SERRE 12/07/10 21:15:18 7431 2 : SUBROUTINE HEA010(T,TC,HEAS,HEAF,DTCSDT,MP,TRAV) 3 : C--------------------------------------------------------------------- 4 : C Calcul du coefficient 'ki' de RK4 (voir HEA011) 5 : C--------------------------------------------------------------------- 6 : C 7 : C--------------------------- 8 : C Parametres Entree/Sortie : 9 : C--------------------------- 10 : C 11 : C E/ T : flottant : Temps courant au pas du RK4 considéré (s) 12 : C E/ TC : flottant : Température des plaques au pas du RK4 13 : C considéré (K) 14 : C /S HEAS : flottant : Puissance générée dans le solide (J/s) 15 : C /S HEAF : flottant : Puissance transmise au fluide (J/s) 16 : C /S DTCSDT : flottant : dTc/dt (K/s) 17 : C /S MP : flottant : Débit de mélange traversant le HEATER (kg/s) 18 : C E/ TRAV : entier : Pointeur du segment de travail 19 : C 20 : C------------------------------ 21 : C Variables de TRAV utilisées : 22 : C------------------------------ 23 : C 24 : C E/ PRESSION : Pression à l'entrée du HEATER (Pa) 25 : C E/ TEMPENT : Température à l'entrée du HEATER (K) 26 : C E/ XiMOY : Fraction molaire à l'entrée du HEATER 27 : C E/ EPS_MH2 : Seuil en deça duquel on suppose que la consommation 28 : C d'H2 max est mise à 0. (kg/s) 29 : C E/ EPS_DT : Ecart entre la température des plaques et celle du 30 : C gaz (Tc-T) en deça duquel on suppose qu'il n'y a pas 31 : C d'échange convectif (K) 32 : C E/ U : Vitesse minimale dans le HEATER (= 0.01 m/s) 33 : C E/ S : Surface des plaques (m2) 34 : C E/ SP : Section de passage entrée/sortie fluide du HEATER (m2) 35 : C E/ MC : Masse des plaques (kg) 36 : C E/ CPC : Chaleur spécifique des plaques (J/kg/K) 37 : C 38 : C------------------- 39 : C Autres variables : 40 : C------------------- 41 : C 42 : C MH2 : Consommation fictive d'H2 par les plaques associée à la 43 : C puissance thermique délivrée (kg/s) 44 : C 45 : C--------------------------------------------------------------------- 46 : C 47 : C Langage : ESOPE + FORTRAN 77 48 : C 49 : C Mise en oeuvre : H. Paillère (1997, TTMF) 50 : C 51 : C--------------------------------------------------------------------- 52 : IMPLICIT INTEGER(I-N) 53 : IMPLICIT REAL*8 (A-H,O-Z) 54 : REAL*8 XE(7) 55 : REAL*8 MH2,MP 56 : segment trav 57 : integer iKALP,iMODEL 58 : real*8 e,L,Lch,Dh,S,sp,Ck 59 : real*8 mc,Cpc 60 : real*8 g,R,deltah 61 : real*8 M(nbesp),cstmod(ncst) 62 : real*8 Cpi(nbesp),al 63 : real*8 eps_mh2,eps_dt,eps_con,u 64 : real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT 65 : real*8 XHEMOY,XCO2MOY,XCOMOY 66 : endsegment 67 : C 68 : C Densité du mélange 69 : PE = PRESSION 70 : TE = TEMPENT 71 : XE(1) = XN2MOY 72 : XE(2) = XO2MOY 73 : XE(3) = XH2MOY 74 : XE(4) = XH2OMOY 75 : XE(5) = XHEMOY 76 : XE(6) = XCO2MOY 77 : XE(7) = XCOMOY 78 : ROE = FLA001(PE,TE,XE,TRAV) 79 : C 80 : C Débit minimum dans le PAR thermique (valeur par défaut) 81 : MP = U*ROE*SP 82 : C 83 : C dmh2/dt équivalent à la puissance thermique du HEATER 84 : MH2 = CSTMOD(1) / deltah 85 : C 86 : C 1. Cas sans comsommation d'hydrogène (puissance thermique résiduelle) 87 : IF (MH2 .LT. EPS_MH2) THEN 88 : IF ((TC-TE) .LT. EPS_DT) THEN 89 : C 1.1 et sans echange convectif 90 : HEAS = 0.D0 91 : HEAF = 0.D0 92 : ELSE 93 : C 1.2 et avec echange convectif (bilan d'énergie des plaques) 94 : HEAS = 0.D0 95 : HEAF = -S*FLA009(PE,TE,TC,TE,XE,TRAV)*(TC-TE) 96 : ENDIF 97 : ELSE 98 : C 2. Cas avec comsommation d'hydrogène 99 : IF ((TC-TE) .LT. EPS_DT) THEN 100 : C 2.1 et sans echange convectif 101 : HEAS = DELTAH * MH2 102 : HEAF = 0.D0 103 : ELSE 104 : C 2.2 et avec echange convectif (bilan d'énergie des plaques) 105 : MP = FLA015(T,TC,MH2,TRAV) 106 : HEAS = DELTAH * MH2 107 : HEAF = -S*FLA009(PE,TE,TC,TE,XE,TRAV)*(TC-TE) 108 : ENDIF 109 : ENDIF 110 : HEAS = HEAS * cstmod(2) 111 : DTCSDT = (HEAS + HEAF) / (MC*CPC) 112 : C 113 : RETURN 114 : END 115 : 116 : 117 : 118 : 119 : 120 : 121 : 122 : 123 : 124 :
© Cast3M 2003 - All rights reserved.
Disclaimer