Télécharger clmult2D.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : clmult2D.dgibi
  2. *********************************************************************
  3. * VF, Condition Limites en 2D pour les Equations d'Euler *
  4. * multiespeces *
  5. * Kudriakov S., DM2S/SFME/LTMF, FEVRIER 2003 *
  6. * *
  7. * Beccantini A., DM2S/SFME/LTMF, FEVRIER 2006: *
  8. * BC of RESErvoir type *
  9. *********************************************************************
  10.  
  11. TYEL = 'QUA8' ;
  12. * TYEL = 'TRI6' ;
  13.  
  14. 'OPTION' 'DIME' 2 'ELEM' TYEL 'ISOV' 'SULI'
  15. 'ECHO' 1 'TRAC' 'X' ;
  16.  
  17.  
  18. GRAPH = VRAI ;
  19. GRAPH = FAUX ;
  20.  
  21. ******************
  22. **** MAILLAGE ****
  23. ******************
  24. *
  25. * P2 P3
  26. *
  27.  
  28. DX = 0.5 ;
  29. P2 = 0.0 0.0 ;
  30. P3 = DX DX ;
  31. P2P3 = P2 'DROIT' 1 P3 ;
  32.  
  33. LIGHAU = (P2P3) 'PLUS' (4 '*' ((-1*DX) DX)) ;
  34.  
  35. DOMINT = P2P3 'REGLER' 4 LIGHAU ;
  36.  
  37. DOMINT = 'CHANGER' DOMINT 'QUAF' ;
  38.  
  39. *
  40. **** Creation of DOMAINE tables via the MODEL object
  41. *
  42.  
  43. MDOMINT = 'MODELISER' DOMINT 'NAVIER_STOKES' 'LINE' ;
  44. MP2P3 = 'MODELISER' P2P3 'NAVIER_STOKES' 'LINE' ;
  45.  
  46. $DOMINT = 'DOMA' MDOMINT 'TABLE' ;
  47. $P2P3 = 'DOMA' MP2P3 'TABLE' ;
  48.  
  49. DOMINT = 'DOMA' MDOMINT 'MAILLAGE' ;
  50. P2P3 = 'DOMA' MP2P3 'MAILLAGE' ;
  51.  
  52. PCEN = ($DOMINT . 'CENTRE') 'POIN' 'PROC' ((DX '/' 2) (DX /2)) ;
  53.  
  54. MDICEN = $DOMINT . 'CENTRE' ;
  55. MFCEN = $P2P3 . 'CENTRE' ;
  56.  
  57. 'SI' GRAPH ;
  58. 'TRACER' (DOMINT 'ET' PCEN) ;
  59. 'FINSI' ;
  60.  
  61. **********************************************************************************
  62. **********************************************************************************
  63. *************************** Riemann invariants, inlet ***************************
  64. **********************************************************************************
  65. **********************************************************************************
  66. NESP = 4 ;
  67.  
  68. *-----------------------------------------------------
  69. *** GAS: H_2, O_2, H_2O, N_2
  70. *
  71. * CP, CV en J/Kg/K @ T = 3000
  72. *-----------------------------------------------------
  73.  
  74. PGAS = 'TABLE' ;
  75.  
  76. PGAS . 'CP' = 'TABLE' ;
  77. PGAS . 'CP' . 'H2 ' = .18729066D+05 ;
  78. PGAS . 'CP' . 'O2 ' = .11886820D+04 ;
  79. PGAS . 'CP' . 'H2O ' = .31209047D+04 ;
  80. PGAS . 'CP' . 'N2 ' = .12993995D+04 ;
  81.  
  82.  
  83.  
  84. PGAS . 'CV' = 'TABLE' ;
  85. PGAS . 'CV' . 'H2 ' = .14571861D+05 ;
  86. PGAS . 'CV' . 'O2 ' = .92885670D+03 ;
  87. PGAS . 'CV' . 'H2O ' = .26589930D+04 ;
  88. PGAS . 'CV' . 'N2 ' = .10024563D+04;
  89.  
  90. *----------------------------------------------------
  91. **** Especes qui sont dans les equations d'Euler
  92. *----------------------------------------------------
  93.  
  94. PGAS . 'ESPEULE' = 'MOTS' 'H2 ' 'O2 ' 'H2O ' ;
  95.  
  96. *----------------------------------------------------
  97. **** Espece qui n'y est pas
  98. *----------------------------------------------------
  99.  
  100.  
  101. PGAS . 'ESPNEULE' = 'N2 ' ;
  102.  
  103. *-----------------
  104. *** CP ****
  105. *-----------------
  106. CP = ('MANUEL' 'CHPO' MDICEN 1 'H2'
  107. (PGAS . 'CP' . 'H2 ') 'NATU' 'DISCRET') 'ET'
  108. ('MANUEL' 'CHPO' MDICEN 1 'O2'
  109. (PGAS . 'CP' . 'O2 ') 'NATU' 'DISCRET') 'ET'
  110. ('MANUEL' 'CHPO' MDICEN 1 'H2O'
  111. (PGAS . 'CP' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  112. ('MANUEL' 'CHPO' MDICEN 1 'N2'
  113. (PGAS . 'CP' . 'N2 ') 'NATU' 'DISCRET') ;
  114. ***************************
  115. CPCEL = ('MANUEL' 'CHPO' MFCEN 1 'H2'
  116. (PGAS . 'CP' . 'H2 ') 'NATU' 'DISCRET') 'ET'
  117. ('MANUEL' 'CHPO' MFCEN 1 'O2'
  118. (PGAS . 'CP' . 'O2 ') 'NATU' 'DISCRET') 'ET'
  119. ('MANUEL' 'CHPO' MFCEN 1 'H2O'
  120. (PGAS . 'CP' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  121. ('MANUEL' 'CHPO' MFCEN 1 'N2'
  122. (PGAS . 'CP' . 'N2 ') 'NATU' 'DISCRET') ;
  123. *-----------------
  124. *** CV ****
  125. *-----------------
  126. CV = ('MANUEL' 'CHPO' MDICEN 1 'H2'
  127. (PGAS . 'CV' . 'H2 ') 'NATU' 'DISCRET') 'ET'
  128. ('MANUEL' 'CHPO' MDICEN 1 'O2'
  129. (PGAS . 'CV' . 'O2 ') 'NATU' 'DISCRET') 'ET'
  130. ('MANUEL' 'CHPO' MDICEN 1 'H2O'
  131. (PGAS . 'CV' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  132. ('MANUEL' 'CHPO' MDICEN 1 'N2'
  133. (PGAS . 'CV' . 'N2 ') 'NATU' 'DISCRET') ;
  134. *****************************
  135. CVCEL = ('MANUEL' 'CHPO' MFCEN 1 'H2'
  136. (PGAS . 'CV' . 'H2 ') 'NATU' 'DISCRET') 'ET'
  137. ('MANUEL' 'CHPO' MFCEN 1 'O2'
  138. (PGAS . 'CV' . 'O2 ') 'NATU' 'DISCRET') 'ET'
  139. ('MANUEL' 'CHPO' MFCEN 1 'H2O'
  140. (PGAS . 'CV' . 'H2O ') 'NATU' 'DISCRET') 'ET'
  141. ('MANUEL' 'CHPO' MFCEN 1 'N2'
  142. (PGAS . 'CV' . 'N2 ') 'NATU' 'DISCRET') ;
  143. ***********************************************************
  144. ***********************************************************
  145.  
  146. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.110 ;
  147. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.41 ;
  148. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 0.1 'UX' 1.65 ;
  149. YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
  150. 'H2O' 0.1 ;
  151. YN1 = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
  152. 'H2O' 0.1 'NATU' 'DISCRET' ;
  153. **********************************************************
  154. UNCH = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.0 ;
  155. YNESP = UNCH '-' ('PSCAL' UNCH YN ('MOTS' 'SCAL' 'SCAL' 'SCAL')
  156. (PGAS . 'ESPEULE')) ;
  157. YNESP = 'NOMC' 'N2 ' YNESP 'NATU' 'DISCRET' ;
  158. YNTOT = YN1 'ET' YNESP ;
  159. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  160. CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  161. CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  162. GAMN = CPTOT '/' CVTOT ;
  163. **********************************************************
  164.  
  165. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' 2.11 'PN' 111.4
  166. 'UX' 2.0 'UY' 1.0 'H2' 0.1 'O2' 0.7 'H2O' 0.05 ;
  167.  
  168.  
  169. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
  170. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
  171.  
  172. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  173. PGAS LISTC LISTP
  174. RN VN PN YN CHPLI 'INRI' ;
  175. ************************************************************
  176. **** On calcule les variables conservatives **************
  177. ************************************************************
  178. RAZM = NESP '+' 3 ;
  179.  
  180. GN = RN '*' VN ('MOTS' 'SCAL' 'SCAL')
  181. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY') ;
  182. LMOT1 = 'MOTS' 'UX' 'UY' ;
  183. VVV = 'PSCAL' VN VN LMOT1 LMOT1 ;
  184. ECIN = 0.5D0 '*' RN '*' VVV ;
  185.  
  186. RETN = 'NOMC' 'SCAL' ((PN '/' (GAMN '-' 1.0)) '+' ECIN)
  187. 'NATU' 'DISCRET' ;
  188.  
  189. LMOT1 = (PGAS . 'ESPEULE') ;
  190. LMOT2 = 'MOTS' 'SCAL' 'SCAL' 'SCAL' ;
  191. LMOT3 = (PGAS . 'ESPEULE') ;
  192.  
  193. RYN = YN '*' RN LMOT1 LMOT2 LMOT3 ;
  194.  
  195. VIT PRES TTEMP FRAM GAMM = 'PRIM' 'PERFMULT' PGAS RN GN RETN RYN ;
  196.  
  197.  
  198. ERR1 = 'MAXIMUM' (PRES '-' PN) 'ABS' ;
  199. ERR2 = 'MAXIMUM' (GN '-' (RN '*' VIT ('MOTS' 'SCAL' 'SCAL' )
  200. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) 'ABS' ;
  201. XXTOL = 1.D-12;
  202. 'SI' (ERR1 > XXTOL) ;
  203. CHAI1 = CHAI 'ERREUR NUMERO 1 : ' ERR1 'POUR ' XXTOL;
  204. 'MESS' CHAI1;
  205. 'ERREUR' 5 ;
  206. 'FINSI' ;
  207.  
  208. 'SI' (ERR2 > XXTOL) ;
  209. CHAI2 = CHAI 'ERREUR NUMERO 2 : ' ERR2 'POUR ' XXTOL;
  210. 'MESS' CHAI2;
  211. 'ERREUR' 5 ;
  212. 'FINSI' ;
  213.  
  214. *****************************************
  215. **** On calcule RCHLIM en gibiane
  216. *****************************************
  217.  
  218. UNC = (('EXCO' VN 'UY' 'SCAL') '-' ('EXCO' VN 'UX' 'SCAL')) '/' (2
  219. '**' 0.5) ;
  220.  
  221. GAMMA = 'MAXIMUM' GAMN ;
  222. GM1 = '-' GAMMA 1. ;
  223. USGM1 = 1. '/' GM1 ;
  224. DSGM1 = 2. '*' USGM1 ;
  225.  
  226. ASONC = GAMMA * (PN '/' RN) ;
  227. ASONC = ASONC '**' 0.5 ;
  228.  
  229. G1 = UNC '-' (DSGM1 * ASONC) ;
  230. G1 = 'EXTRAIRE' G1 PCEN 'SCAL' ;
  231. *******************************************************
  232. RNCEL= 'EXCO' CHPLI 'RN' ;
  233. PNCEL = 'EXCO' CHPLI 'PN' ;
  234. UNCELC = (('EXCO' CHPLI 'UY' 'SCAL') '-' ('EXCO' CHPLI 'UX' 'SCAL'))
  235. '/' (2 '**' 0.5) ;
  236. UTCELC = (('EXCO' CHPLI 'UY' 'SCAL') '+' ('EXCO' CHPLI 'UX' 'SCAL'))
  237. '/' (-1*(2 '**' 0.5)) ;
  238. YH2CEL = 'EXCO' CHPLI 'H2' 'NATU' 'DISCRET' ;
  239. YO2CEL = 'EXCO' CHPLI 'O2' 'NATU' 'DISCRET' ;
  240. YHOCEL = 'EXCO' CHPLI 'H2O' 'NATU' 'DISCRET' ;
  241. YNCEL = YH2CEL 'ET' YO2CEL ;
  242. YNCEL = YNCEL 'ET' YHOCEL ;
  243. **********************************************************
  244. YNTCEL = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 4 'H2' 0.1 'O2' 0.7
  245. 'H2O' 0.05 'N2' 0.15 'NATU' 'DISCRET' ;
  246. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  247. CVTOT = 'PSCAL' YNTCEL CVCEL LMOT1 LMOT1 ;
  248. CPTOT = 'PSCAL' YNTCEL CPCEL LMOT1 LMOT1 ;
  249. GAMCEL = CPTOT '/' CVTOT ;
  250. ***********************************************************
  251. GAMF = 'MAXIMUM' GAMCEL ;
  252. GM1F = '-' GAMF 1. ;
  253. USGM1F = 1. '/' GM1F ;
  254. DSGM1F = 2. '*' USGM1F ;
  255.  
  256. ASONF = GAMF * (PNCEL '/' RNCEL) ;
  257. ASONF = ASONF '**' 0.5 ;
  258.  
  259. G3 = UNCELC '+' (DSGM1F * ASONF) ;
  260. G3 = 'EXTRAIRE' G3 ('POIN' 1 ($P2P3 . 'CENTRE')) 'SCAL' ;
  261. S = PNCEL '/' (RNCEL** GAMF) ;
  262. S = 'EXTRAIRE' S ('POIN' 1 ($P2P3 . 'CENTRE')) 'SCAL' ;
  263.  
  264. GAMID = (GAMMA '+' GAMF) '/' 2.0D0 ;
  265. GM1M = '-' GAMID 1. ;
  266. USGM1M = 1. '/' GM1M ;
  267. DSGM1M = 2. '*' USGM1M ;
  268.  
  269. ASON2=(0.5D0*(G3-G1)) ;
  270. ASON2=ASON2/DSGM1F ;
  271. ASON2=ASON2*ASON2 ;
  272. RHO=ASON2/(GAMF*S) ;
  273. RHO=RHO**USGM1F ;
  274. P=(RHO*ASON2)/GAMF ;
  275. UNCEL = 0.5 '*' (G3 '+' G1) ;
  276. UTCEL = 'EXTRAIRE' UTCELC ('POIN' 1 ($P2P3 . 'CENTRE')) 'SCAL' ;
  277.  
  278. UX = (UNCEL '+' UTCEL) '/' (-1*(2 '**' 0.5)) ;
  279. UY = (UNCEL '-' UTCEL) '/' (2 '**' 0.5) ;
  280.  
  281.  
  282. CHLIM2 = 'MANU' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' RHO 'UX'
  283. UX 'UY' UY 'PN' P 'H2' 0.1 'O2' 0.7 'H2O' 0.05 ;
  284.  
  285. ERRO = 'MAXIMUM' (CHLIM2 '-' RCHLIM) ;
  286.  
  287. XXTOL = 1.D-12;
  288. 'SI' (ERRO > XXTOL) ;
  289. CHAI1 = CHAI 'ERREUR NUMERO 3 : ' ERRO 'POUR ' XXTOL;
  290. 'MESS' CHAI1;
  291. 'ERREUR' 5 ;
  292. 'FINSI' ;
  293.  
  294. 'FINSI' ;
  295.  
  296. ******************************************************************
  297. **** Jacobian with respect to primitive variable
  298. ******************************************************************
  299.  
  300. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
  301. PGAS LISTC LISTP
  302. RN VN PN YN CHPLI 'INRI' ;
  303. EPS=1.0D-6 ;
  304. ERROM = 1.0D-3 ;
  305. * RN ****************** RN ******* RN *********************
  306.  
  307. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  308. ELE1 = 'MANUEL' 'POI1' PCEN ;
  309. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  310.  
  311. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  312. PGAS LISTC LISTP
  313. (RN '+' DR) VN PN YN CHPLI 'INRI' ;
  314.  
  315. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  316. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
  317. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  318. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  319.  
  320. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  321.  
  322. XXTOL = ERROM;
  323. ERR1 = ERRO;
  324. 'SI' (ERR1 > XXTOL) ;
  325. CHAI1 = CHAI 'ERREUR NUMERO 4 : ' ERR1 'POUR ' XXTOL;
  326. 'MESS' CHAI1;
  327. 'ERREUR' 5 ;
  328. 'FINSI' ;
  329.  
  330.  
  331. * UX ***** UX ****** UX ****** UX ************************
  332.  
  333. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  334. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  335.  
  336. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  337. PGAS LISTC LISTP
  338. RN (VN '+' DUX) PN YN CHPLI 'INRI' ;
  339.  
  340. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  341. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  342. 'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0) ;
  343. MMM = 'MAXIMUM' DRESDUX 'ABS' ;
  344.  
  345. ERRO = 'MAXIMUM' ((DRESDUX '-' DRESDUXN) '/' MMM) 'ABS' ;
  346.  
  347. XXTOL = ERROM;
  348. ERR1 = ERRO;
  349. 'SI' (ERR1 > XXTOL) ;
  350. CHAI1 = CHAI 'ERREUR NUMERO 5 : ' ERR1 'POUR ' XXTOL;
  351. 'MESS' CHAI1;
  352. 'ERREUR' 5 ;
  353. 'FINSI' ;
  354.  
  355.  
  356. *** UY *************** UY ************* UY *******************
  357.  
  358. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  359. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  360.  
  361. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  362. PGAS LISTC LISTP
  363. RN (VN '+' DUY) PN YN CHPLI 'INRI' ;
  364.  
  365. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  366. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  367. 'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  368.  
  369. MMM = 'MAXIMUM' DRESDUY 'ABS' ;
  370.  
  371. ERRO = 'MAXIMUM' ((DRESDUY '-' DRESDUYN) '/' MMM) 'ABS' ;
  372.  
  373. XXTOL = ERROM;
  374. ERR1 = ERRO;
  375. 'SI' (ERR1 > XXTOL) ;
  376. CHAI1 = CHAI 'ERREUR NUMERO 6 : ' ERR1 'POUR ' XXTOL;
  377. 'MESS' CHAI1;
  378. 'ERREUR' 5 ;
  379. 'FINSI' ;
  380.  
  381.  
  382. **** PN ********** PN ************ PN ***********************
  383.  
  384. DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  385. ELE1 = 'MANUEL' 'POI1' PCEN ;
  386. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  387.  
  388. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  389. PGAS LISTC LISTP
  390. RN VN (PN '+' DP) YN CHPLI 'INRI' ;
  391.  
  392. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  393. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  394. 'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  395.  
  396. MMM = 'MAXIMUM' DRESDP 'ABS' ;
  397.  
  398. ERRO = 'MAXIMUM' ((DRESDP '-' DRESDPN)'/' MMM) 'ABS' ;
  399.  
  400.  
  401. XXTOL = ERROM;
  402. ERR1 = ERRO;
  403. 'SI' (ERR1 > XXTOL) ;
  404. CHAI1 = CHAI 'ERREUR NUMERO 7 : ' ERR1 'POUR ' XXTOL;
  405. 'MESS' CHAI1;
  406. 'ERREUR' 5 ;
  407. 'FINSI' ;
  408.  
  409.  
  410. **** YH2 ********** YH2 ************ YH2 ***********************
  411.  
  412. DYH2SCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
  413. ELE1 = 'MANUEL' 'POI1' PCEN ;
  414. DYH2 = ('MANUEL' 'CHPO' ELE1 1 'H2' DYH2SCAL ) ;
  415.  
  416. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  417. PGAS LISTC LISTP
  418. RN VN PN (YN1 '+' DYH2) CHPLI 'INRI' ;
  419.  
  420. DRESDYN = ( RCHRES1 '-' RCHRES) '/' DYH2SCAL ;
  421. DRESDY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  422. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
  423.  
  424. MMM = 'MAXIMUM' DRESDY 'ABS' ;
  425.  
  426. ERRO = 'MAXIMUM' ((DRESDY '-' DRESDYN)'/' MMM) 'ABS' ;
  427.  
  428. XXTOL = ERROM;
  429. ERR1 = ERRO;
  430. 'SI' (ERR1 > XXTOL) ;
  431. CHAI1 = CHAI 'ERREUR NUMERO 8 : ' ERR1 'POUR ' XXTOL;
  432. 'MESS' CHAI1;
  433. 'ERREUR' 5 ;
  434. 'FINSI' ;
  435.  
  436.  
  437. ***** YO2 ****** YO2 *********** YO2 ******************
  438.  
  439. DYO2SCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
  440.  
  441. ELE1 = 'MANUEL' 'POI1' PCEN ;
  442. DYO2 = ('MANUEL' 'CHPO' ELE1 1 'O2' DYO2SCAL ) ;
  443.  
  444. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  445. PGAS LISTC LISTP
  446. RN VN PN (YN1 '+' DYO2) CHPLI 'INRI' ;
  447.  
  448. DRESDYN = ( RCHRES1 '-' RCHRES) '/' DYO2SCAL ;
  449. DRESDY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  450. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
  451.  
  452. MMM = 'MAXIMUM' DRESDY 'ABS' ;
  453.  
  454. ERRO = 'MAXIMUM' ((DRESDY '-' DRESDYN)'/' MMM) 'ABS' ;
  455.  
  456. XXTOL = ERROM;
  457. ERR1 = ERRO;
  458. 'SI' (ERR1 > XXTOL) ;
  459. CHAI1 = CHAI 'ERREUR NUMERO 9 : ' ERR1 'POUR ' XXTOL;
  460. 'MESS' CHAI1;
  461. 'ERREUR' 5 ;
  462. 'FINSI' ;
  463.  
  464.  
  465. ***** YHO ****** YHO *********** YHO ******************
  466.  
  467. DYHOSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
  468.  
  469. ELE1 = 'MANUEL' 'POI1' PCEN ;
  470. DYHO = ('MANUEL' 'CHPO' ELE1 1 'H2O' DYHOSCAL ) ;
  471.  
  472. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  473. PGAS LISTC LISTP
  474. RN VN PN (YN1 '+' DYHO) CHPLI 'INRI' ;
  475.  
  476. DRESDYN = ( RCHRES1 '-' RCHRES) '/' DYHOSCAL ;
  477. DRESDY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  478. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
  479.  
  480. MMM = 'MAXIMUM' DRESDY 'ABS' ;
  481.  
  482. ERRO = 'MAXIMUM' ((DRESDY '-' DRESDYN)'/' MMM) 'ABS' ;
  483.  
  484. XXTOL = ERROM;
  485. ERR1 = ERRO;
  486. 'SI' (ERR1 > XXTOL) ;
  487. CHAI1 = CHAI 'ERREUR NUMERO 10: ' ERR1 'POUR ' XXTOL;
  488. 'MESS' CHAI1;
  489. 'ERREUR' 5 ;
  490. 'FINSI' ;
  491.  
  492.  
  493. **************************************************************
  494. **** Jacobian with respect to conservative variables
  495. **************************************************************
  496.  
  497. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
  498. PGAS LISTC LISTP
  499. RN VN PN YN CHPLI 'INRI' ;
  500.  
  501. ******** RN ************* RN ********** RN **************
  502.  
  503. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  504. ELE1 = 'MANUEL' 'POI1' PCEN ;
  505. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  506.  
  507. RN1 = RN '+' DR ;
  508. GN1 = 'COPIER' GN ;
  509. RETN1 = 'COPIER' RETN ;
  510. RYN1 = 'COPIER' RYN ;
  511.  
  512. VN1 PN1 TTEMP YN1 GAMM =
  513. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  514.  
  515. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  516. PGAS LISTC LISTP
  517. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  518.  
  519. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  520. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
  521. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
  522. 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  523.  
  524. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  525. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  526.  
  527. XXTOL = ERROM;
  528. ERR1 = ERRO;
  529. 'SI' (ERR1 > XXTOL) ;
  530. CHAI1 = CHAI 'ERREUR NUMERO 11: ' ERR1 'POUR ' XXTOL;
  531. 'MESS' CHAI1;
  532. 'ERREUR' 5 ;
  533. 'FINSI' ;
  534.  
  535. ***** RUX ******* RUX ******** RUX *******************
  536.  
  537. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  538. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  539.  
  540. RN1 = 'COPIER' RN ;
  541. GN1 = GN '+' DUX ;
  542. RETN1 = 'COPIER' RETN ;
  543. RYN1 = 'COPIER' RYN ;
  544.  
  545. VN1 PN1 TTEMP YN1 GAMM =
  546. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  547.  
  548. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  549. PGAS LISTC LISTP
  550. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  551.  
  552. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  553. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  554. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0
  555. 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  556.  
  557. MMM = 'MAXIMUM' DRESDUX 'ABS' ;
  558. ERRO = 'MAXIMUM' ((DRESDUX '-' DRESDUXN) '/' MMM) 'ABS' ;
  559.  
  560. XXTOL = ERROM;
  561. ERR1 = ERRO;
  562. 'SI' (ERR1 > XXTOL) ;
  563. CHAI1 = CHAI 'ERREUR NUMERO 12: ' ERR1 'POUR ' XXTOL;
  564. 'MESS' CHAI1;
  565. 'ERREUR' 5 ;
  566. 'FINSI' ;
  567.  
  568. ***** RUY ******* RUY ******** RUY *******************
  569.  
  570. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  571. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  572.  
  573. RN1 = 'COPIER' RN ;
  574. GN1 = GN '+' DUY ;
  575. RETN1 = 'COPIER' RETN ;
  576. RYN1 = 'COPIER' RYN ;
  577.  
  578. VN1 PN1 TTEMP YN1 GAMM =
  579. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  580.  
  581. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  582. PGAS LISTC LISTP
  583. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  584.  
  585. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  586. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  587. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0
  588. 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  589.  
  590. MMM = 'MAXIMUM' DRESDUY 'ABS' ;
  591. ERRO = 'MAXIMUM' ((DRESDUY '-' DRESDUYN) '/' MMM) 'ABS' ;
  592.  
  593. XXTOL = ERROM;
  594. ERR1 = ERRO;
  595. 'SI' (ERR1 > XXTOL) ;
  596. CHAI1 = CHAI 'ERREUR NUMERO 13: ' ERR1 'POUR ' XXTOL;
  597. 'MESS' CHAI1;
  598. 'ERREUR' 5 ;
  599. 'FINSI' ;
  600.  
  601. ******** RETN ************ RETN ************************
  602.  
  603. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  604. ELE1 = 'MANUEL' 'POI1' PCEN ;
  605. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  606.  
  607. RN1 = 'COPIER' RN ;
  608. GN1 = 'COPIER' GN ;
  609. RETN1 = RETN '+' DR ;
  610. RYN1 = 'COPIER' RYN ;
  611.  
  612. VN1 PN1 TTEMP YN1 GAMM =
  613. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  614.  
  615. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  616. PGAS LISTC LISTP
  617. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  618.  
  619. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  620. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  621. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0
  622. 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  623.  
  624. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  625. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  626.  
  627. XXTOL = ERROM;
  628. ERR1 = ERRO;
  629. 'SI' (ERR1 > XXTOL) ;
  630. CHAI1 = CHAI 'ERREUR NUMERO 14: ' ERR1 'POUR ' XXTOL;
  631. 'MESS' CHAI1;
  632. 'ERREUR' 5 ;
  633. 'FINSI' ;
  634.  
  635. ******** RYN ************ RYN ************************
  636.  
  637. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
  638. ELE1 = 'MANUEL' 'POI1' PCEN ;
  639. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  640.  
  641. RN1 = 'COPIER' RN ;
  642. GN1 = 'COPIER' GN ;
  643. RETN1 = RETN ;
  644. RYN1 = RYN '+' DR ;
  645.  
  646. VN1 PN1 TTEMP YN1 GAMM =
  647. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  648.  
  649. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  650. PGAS LISTC LISTP
  651. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  652.  
  653. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  654. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  655. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
  656. 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  657.  
  658. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  659. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  660.  
  661. XXTOL = ERROM;
  662. ERR1 = ERRO;
  663. 'SI' (ERR1 > XXTOL) ;
  664. CHAI1 = CHAI 'ERREUR NUMERO 15: ' ERR1 'POUR ' XXTOL;
  665. 'MESS' CHAI1;
  666. 'ERREUR' 5 ;
  667. 'FINSI' ;
  668.  
  669. ******** RYN ************ RYN ************************
  670.  
  671. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
  672. ELE1 = 'MANUEL' 'POI1' PCEN ;
  673. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  674.  
  675. RN1 = 'COPIER' RN ;
  676. GN1 = 'COPIER' GN ;
  677. RETN1 = RETN ;
  678. RYN1 = RYN '+' DR ;
  679.  
  680. VN1 PN1 TTEMP YN1 GAMM =
  681. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  682.  
  683. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  684. PGAS LISTC LISTP
  685. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  686.  
  687. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  688. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  689. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
  690. 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  691.  
  692. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  693. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  694.  
  695. XXTOL = ERROM;
  696. ERR1 = ERRO;
  697. 'SI' (ERR1 > XXTOL) ;
  698. CHAI1 = CHAI 'ERREUR NUMERO 16: ' ERR1 'POUR ' XXTOL;
  699. 'MESS' CHAI1;
  700. 'ERREUR' 5 ;
  701. 'FINSI' ;
  702.  
  703. ******** RYNH2 ************ RYNH2 ************************
  704.  
  705. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
  706. ELE1 = 'MANUEL' 'POI1' PCEN ;
  707. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  708.  
  709. RN1 = 'COPIER' RN ;
  710. GN1 = 'COPIER' GN ;
  711. RETN1 = RETN ;
  712. RYN1 = RYN '+' DR ;
  713.  
  714. VN1 PN1 TTEMP YN1 GAMM =
  715. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  716.  
  717. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  718. PGAS LISTC LISTP
  719. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  720.  
  721. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  722. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  723. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
  724. 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  725.  
  726. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  727. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  728.  
  729. XXTOL = ERROM;
  730. ERR1 = ERRO;
  731. 'SI' (ERR1 > XXTOL) ;
  732. CHAI1 = CHAI 'ERREUR NUMERO 17: ' ERR1 'POUR ' XXTOL;
  733. 'MESS' CHAI1;
  734. 'ERREUR' 5 ;
  735. 'FINSI' ;
  736.  
  737. ******** RYNO2 ************ RYNO2 ************************
  738.  
  739. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
  740. ELE1 = 'MANUEL' 'POI1' PCEN ;
  741. DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
  742.  
  743. RN1 = 'COPIER' RN ;
  744. GN1 = 'COPIER' GN ;
  745. RETN1 = RETN ;
  746. RYN1 = RYN '+' DR ;
  747.  
  748. VN1 PN1 TTEMP YN1 GAMM =
  749. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  750.  
  751. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  752. PGAS LISTC LISTP
  753. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  754.  
  755. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  756. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  757. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
  758. 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
  759.  
  760. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  761. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  762.  
  763. XXTOL = ERROM;
  764. ERR1 = ERRO;
  765. 'SI' (ERR1 > XXTOL) ;
  766. CHAI1 = CHAI 'ERREUR NUMERO 18: ' ERR1 'POUR ' XXTOL;
  767. 'MESS' CHAI1;
  768. 'ERREUR' 5 ;
  769. 'FINSI' ;
  770.  
  771.  
  772. ******** RYNH20 ************ RYNH20 ************************
  773.  
  774. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
  775. ELE1 = 'MANUEL' 'POI1' PCEN ;
  776. DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
  777.  
  778. RN1 = 'COPIER' RN ;
  779. GN1 = 'COPIER' GN ;
  780. RETN1 = RETN ;
  781. RYN1 = RYN '+' DR ;
  782.  
  783. VN1 PN1 TTEMP YN1 GAMM =
  784. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  785.  
  786. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  787. PGAS LISTC LISTP
  788. RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
  789.  
  790. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  791. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  792. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
  793. 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
  794.  
  795. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  796. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  797.  
  798. XXTOL = ERROM;
  799. ERR1 = ERRO;
  800. 'SI' (ERR1 > XXTOL) ;
  801. CHAI1 = CHAI 'ERREUR NUMERO 19: ' ERR1 'POUR ' XXTOL;
  802. 'MESS' CHAI1;
  803. 'ERREUR' 5 ;
  804. 'FINSI' ;
  805.  
  806. ******************************************************************************
  807. ******************************************************************************
  808. *************************** INJE ********************************************
  809. ******************************************************************************
  810. ******************************************************************************
  811. *
  812. * We impose RT and the normal momentum at the inlet.
  813. * (les fractions massiques sont impsees aussi)
  814. ***** Calcul du residu dans le cas particulier UF = UC
  815. *****************************************************************************
  816. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  817. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  818. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 15. 'UX' -15. ;
  819. YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
  820. 'H2O' 0.1 ;
  821.  
  822. GN = RN '*' VN ;
  823. MOME = (15 '*' (2. '**' 0.5)) * 1.1 ;
  824. RT = 2.1 '/' 1.1 ;
  825.  
  826. UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
  827. UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
  828. PC = ('EXTRAIRE' PN 'SCAL' PCEN) ;
  829. RC = ('EXTRAIRE' RN 'SCAL' PCEN) ;
  830. YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
  831. YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
  832. YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
  833.  
  834. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'RT' RT
  835. 'MOME' MOME 'H2' 0.3 'O2' 0.5 'H2O' 0.1 ;
  836.  
  837.  
  838. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
  839. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
  840.  
  841. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  842. PGAS LISTC LISTP
  843. RN VN PN YN CHPLI 'INJE' ;
  844. *---------------------------------------
  845. ** Testing RCHLIM !!!!!
  846. *---------------------------------------
  847. RCHLIM1 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' RC 'PN' PC 'UX' UXC
  848. 'UY' UYC 'H2' 0.3 'O2' 0.5 'H2O' 0.1 ;
  849.  
  850. ERRO = 'MAXIMUM' (RCHLIM1 '-' RCHLIM) 'ABS' ;
  851.  
  852. XXTOL = 1.D-12;
  853. ERR1 = ERRO;
  854. 'SI' (ERR1 > XXTOL) ;
  855. CHAI1 = CHAI 'ERREUR NUMERO 20: ' ERR1 'POUR ' XXTOL;
  856. 'MESS' CHAI1;
  857. 'ERREUR' 5 ;
  858. 'FINSI' ;
  859.  
  860. *-------------------------------------
  861. ** Testing RCHRES !!!!!
  862. *-------------------------------------
  863. 'SI' ('EGA' TYEL 'QUA8') ;
  864. SURSVOL = 1. '/' ((2 '**' 0.5) * 0.5) ;
  865. 'SINON' ;
  866. SURSVOL = 2. '/' ((2 '**' 0.5) * 0.5) ;
  867. 'FINSI' ;
  868. *-------------------------------------------------------------
  869. YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
  870. YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
  871. 'H2O' YH2O 'N2' YN2 ;
  872. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  873. CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  874. CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  875. GAMN = CPTOT '/' CVTOT ;
  876. GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
  877. *--------------------------------------------------------------
  878. FMASS = MOME ;
  879. FUX = (MOME * UXC) '-' (PC '/' (2. '**' 0.5)) ;
  880. FUY = (MOME * UYC) '+' (PC '/' (2. '**' 0.5)) ;
  881. HIN = (GAMC/(GAMC-1)) '*' (PC '/' RC) ;
  882. ECIN = 0.5 '*' ((UXC '*' UXC) '+' (UYC '*' UYC)) ;
  883. FEN = MOME '*' (HIN '+' ECIN)
  884. ;
  885. RCHRES1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' PCEN) 7
  886. 'VAR1' (SURSVOL '*' FMASS) 'VAR2' (FUX '*' SURSVOL) 'VAR3' (FUY '*'
  887. SURSVOL) 'VAR4' (FEN '*' SURSVOL) 'VAR5' (FMASS '*' YH2 '*' SURSVOL)
  888. 'VAR6' (FMASS '*' YO2'*'SURSVOL) 'VAR7' (FMASS '*' YH2O'*'SURSVOL);
  889.  
  890. ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
  891.  
  892. XXTOL = 1.D-10;
  893. ERR1 = ERRO;
  894. 'SI' (ERR1 > XXTOL) ;
  895. CHAI1 = CHAI 'ERREUR NUMERO 21: ' ERR1 'POUR ' XXTOL;
  896. 'MESS' CHAI1;
  897. 'ERREUR' 5 ;
  898. 'FINSI' ;
  899.  
  900. ****************************************************************
  901. ************* JACOBIANS ***************************************
  902. ****************************************************************
  903.  
  904. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  905. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  906. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 1.11 'UX' -1.3 ;
  907. YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
  908. 'H2O' 0.1 ;
  909. YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
  910. YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
  911. YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
  912. YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
  913. YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
  914. 'H2O' YH2O 'N2' YN2 ;
  915. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  916. CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  917. CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  918. GAMN = CPTOT '/' CVTOT ;
  919. GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
  920.  
  921. LMOT1 = (PGAS . 'ESPEULE') ;
  922. LMOT2 = 'MOTS' 'SCAL' 'SCAL' 'SCAL' ;
  923. LMOT3 = (PGAS . 'ESPEULE') ;
  924.  
  925. RYN = YN '*' RN LMOT1 LMOT2 LMOT3 ;
  926.  
  927. ************************************************
  928.  
  929. GN = RN '*' VN ;
  930. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  931. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
  932. MOME = 0.111 ;
  933. RT = 4.1 ;
  934.  
  935. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'RT' RT
  936. 'MOME' MOME 'H2' 0.1 'O2' 0.5 'H2O' 0.3 ;
  937.  
  938.  
  939. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
  940. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
  941.  
  942. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  943. PGAS LISTC LISTP
  944. RN VN PN YN CHPLI 'INJE' ;
  945.  
  946. ***************************************************************
  947. ***** Jacobian with respect to the primitive variables ********
  948. ***************************************************************
  949.  
  950. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
  951. PGAS LISTC LISTP
  952. RN VN PN YN CHPLI 'INJE' ;
  953.  
  954. * EPS=1.0D-6 ;
  955. * ERROM = 1.0D-5 ;
  956.  
  957. ************ RN *************** RN *****************************
  958.  
  959. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  960. ELE1 = 'MANUEL' 'POI1' PCEN ;
  961. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  962.  
  963. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  964. PGAS LISTC LISTP
  965. (RN '+' DR) VN PN YN CHPLI 'INJE' ;
  966.  
  967. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  968. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
  969. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  970.  
  971. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  972.  
  973. XXTOL = ERROM;
  974. ERR1 = ERRO;
  975. 'SI' (ERR1 > XXTOL) ;
  976. CHAI1 = CHAI 'ERREUR NUMERO 22: ' ERR1 'POUR ' XXTOL;
  977. 'MESS' CHAI1;
  978. 'ERREUR' 5 ;
  979. 'FINSI' ;
  980.  
  981.  
  982. ************ UX *************** UX *****************************
  983.  
  984. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  985. ELE1 = 'MANUEL' 'POI1' PCEN ;
  986. DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
  987.  
  988. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  989. PGAS LISTC LISTP
  990. RN (VN '+' DR) PN YN CHPLI 'INJE' ;
  991.  
  992. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  993. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  994. 'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  995.  
  996. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  997.  
  998. XXTOL = ERROM;
  999. ERR1 = ERRO;
  1000. 'SI' (ERR1 > XXTOL) ;
  1001. CHAI1 = CHAI 'ERREUR NUMERO 23: ' ERR1 'POUR ' XXTOL;
  1002. 'MESS' CHAI1;
  1003. 'ERREUR' 5 ;
  1004. 'FINSI' ;
  1005.  
  1006. ************ UY *************** UY *****************************
  1007.  
  1008. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  1009. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1010. DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
  1011.  
  1012. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1013. PGAS LISTC LISTP
  1014. RN (VN '+' DR) PN YN CHPLI 'INJE' ;
  1015.  
  1016. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1017. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1018. 'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  1019.  
  1020. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1021.  
  1022. XXTOL = ERROM;
  1023. ERR1 = ERRO;
  1024. 'SI' (ERR1 > XXTOL) ;
  1025. CHAI1 = CHAI 'ERREUR NUMERO 24: ' ERR1 'POUR ' XXTOL;
  1026. 'MESS' CHAI1;
  1027. 'ERREUR' 5 ;
  1028. 'FINSI' ;
  1029.  
  1030. ************ PN *************** PN *****************************
  1031.  
  1032. DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1033. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1034. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1035.  
  1036. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1037. PGAS LISTC LISTP
  1038. RN VN (PN '+' DR) YN CHPLI 'INJE' ;
  1039.  
  1040. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1041. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1042. 'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  1043.  
  1044. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1045.  
  1046. XXTOL = ERROM;
  1047. ERR1 = ERRO;
  1048. 'SI' (ERR1 > XXTOL) ;
  1049. CHAI1 = CHAI 'ERREUR NUMERO 25: ' ERR1 'POUR ' XXTOL;
  1050. 'MESS' CHAI1;
  1051. 'ERREUR' 5 ;
  1052. 'FINSI' ;
  1053.  
  1054. ************ H2 *************** H2 *****************************
  1055.  
  1056. DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
  1057. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1058. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  1059.  
  1060. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1061. PGAS LISTC LISTP
  1062. RN VN PN (YN '+' DR) CHPLI 'INJE' ;
  1063.  
  1064. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1065. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1066. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
  1067.  
  1068. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1069.  
  1070. XXTOL = ERROM;
  1071. ERR1 = ERRO;
  1072. 'SI' (ERR1 > XXTOL) ;
  1073. CHAI1 = CHAI 'ERREUR NUMERO 26: ' ERR1 'POUR ' XXTOL;
  1074. 'MESS' CHAI1;
  1075. 'ERREUR' 5 ;
  1076. 'FINSI' ;
  1077.  
  1078. ************ O2 *************** O2 *****************************
  1079.  
  1080. DRSCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
  1081. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1082. DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
  1083.  
  1084. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1085. PGAS LISTC LISTP
  1086. RN VN PN (YN '+' DR) CHPLI 'INJE' ;
  1087.  
  1088. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1089. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1090. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
  1091.  
  1092. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1093.  
  1094. XXTOL = ERROM;
  1095. ERR1 = ERRO;
  1096. 'SI' (ERR1 > XXTOL) ;
  1097. CHAI1 = CHAI 'ERREUR NUMERO 27: ' ERR1 'POUR ' XXTOL;
  1098. 'MESS' CHAI1;
  1099. 'ERREUR' 5 ;
  1100. 'FINSI' ;
  1101.  
  1102. ************ H2O *************** H2O *****************************
  1103.  
  1104. DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
  1105. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1106. DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
  1107.  
  1108. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1109. PGAS LISTC LISTP
  1110. RN VN PN (YN '+' DR) CHPLI 'INJE' ;
  1111.  
  1112. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1113. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1114. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
  1115.  
  1116. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1117.  
  1118. XXTOL = ERROM;
  1119. ERR1 = ERRO;
  1120. 'SI' (ERR1 > XXTOL) ;
  1121. CHAI1 = CHAI 'ERREUR NUMERO 28: ' ERR1 'POUR ' XXTOL;
  1122. 'MESS' CHAI1;
  1123. 'ERREUR' 5 ;
  1124. 'FINSI' ;
  1125.  
  1126. ************************************************************
  1127. **** Jacobian with respect to conservative variables ******
  1128. ************************************************************
  1129.  
  1130. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1131. PGAS LISTC LISTP
  1132. RN VN PN YN CHPLI 'INJE' ;
  1133.  
  1134. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
  1135. PGAS LISTC LISTP
  1136. RN VN PN YN CHPLI 'INJE' ;
  1137.  
  1138. ***************** RN ************************************
  1139.  
  1140. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1141. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1142. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1143.  
  1144. RN1 = RN '+' DR ;
  1145. GN1 = 'COPIER' GN ;
  1146. RETN1 = 'COPIER' RETN ;
  1147. RYN1 = 'COPIER' RYN ;
  1148.  
  1149. VN1 PN1 TTEMP YN1 GAMM =
  1150. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1151.  
  1152. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1153. PGAS LISTC LISTP
  1154. RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
  1155.  
  1156. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1157. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
  1158. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1159.  
  1160. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1161.  
  1162. XXTOL = ERROM;
  1163. ERR1 = ERRO;
  1164. 'SI' (ERR1 > XXTOL) ;
  1165. CHAI1 = CHAI 'ERREUR NUMERO 29: ' ERR1 'POUR ' XXTOL;
  1166. 'MESS' CHAI1;
  1167. 'ERREUR' 5 ;
  1168. 'FINSI' ;
  1169.  
  1170. ***************** RUX ************************************
  1171.  
  1172. DRSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  1173. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1174. DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
  1175.  
  1176. RN1 = 'COPIER' RN ;
  1177. GN1 = GN '+' DR ;
  1178. RETN1 = 'COPIER' RETN ;
  1179. RYN1 = 'COPIER' RYN ;
  1180.  
  1181. VN1 PN1 TTEMP YN1 GAMM =
  1182. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1183.  
  1184. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1185. PGAS LISTC LISTP
  1186. RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
  1187.  
  1188. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1189. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1190. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1191.  
  1192. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1193.  
  1194. XXTOL = ERROM;
  1195. ERR1 = ERRO;
  1196. 'SI' (ERR1 > XXTOL) ;
  1197. CHAI1 = CHAI 'ERREUR NUMERO 30: ' ERR1 'POUR ' XXTOL;
  1198. 'MESS' CHAI1;
  1199. 'ERREUR' 5 ;
  1200. 'FINSI' ;
  1201.  
  1202. ***************** RUY ************************************
  1203.  
  1204. DRSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  1205. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1206. DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
  1207.  
  1208. RN1 = 'COPIER' RN ;
  1209. GN1 = GN '+' DR ;
  1210. RETN1 = 'COPIER' RETN ;
  1211. RYN1 = 'COPIER' RYN ;
  1212.  
  1213. VN1 PN1 TTEMP YN1 GAMM =
  1214. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1215.  
  1216. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1217. PGAS LISTC LISTP
  1218. RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
  1219.  
  1220. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1221. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1222. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1223.  
  1224. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1225.  
  1226. XXTOL = ERROM;
  1227. ERR1 = ERRO;
  1228. 'SI' (ERR1 > XXTOL) ;
  1229. CHAI1 = CHAI 'ERREUR NUMERO 31: ' ERR1 'POUR ' XXTOL;
  1230. 'MESS' CHAI1;
  1231. 'ERREUR' 5 ;
  1232. 'FINSI' ;
  1233.  
  1234.  
  1235. ***************** RETN ************************************
  1236.  
  1237. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1238. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1239. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1240.  
  1241. RN1 = 'COPIER' RN ;
  1242. GN1 = 'COPIER' GN ;
  1243. RETN1 = RETN '+' DR ;
  1244. RYN1 = 'COPIER' RYN ;
  1245.  
  1246. VN1 PN1 TTEMP YN1 GAMM =
  1247. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1248.  
  1249. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1250. PGAS LISTC LISTP
  1251. RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
  1252.  
  1253. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1254. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1255. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1256.  
  1257. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1258.  
  1259. XXTOL = ERROM;
  1260. ERR1 = ERRO;
  1261. 'SI' (ERR1 > XXTOL) ;
  1262. CHAI1 = CHAI 'ERREUR NUMERO 32: ' ERR1 'POUR ' XXTOL;
  1263. 'MESS' CHAI1;
  1264. 'ERREUR' 5 ;
  1265. 'FINSI' ;
  1266.  
  1267.  
  1268. ***************** RYH2 ************************************
  1269.  
  1270. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
  1271. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1272. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  1273.  
  1274. RN1 = 'COPIER' RN ;
  1275. GN1 = 'COPIER' GN ;
  1276. RETN1 = 'COPIER' RETN ;
  1277. RYN1 = RYN '+' DR ;
  1278.  
  1279. VN1 PN1 TTEMP YN1 GAMM =
  1280. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1281.  
  1282. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1283. PGAS LISTC LISTP
  1284. RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
  1285.  
  1286. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1287. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1288. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1289.  
  1290. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1291.  
  1292. XXTOL = ERROM;
  1293. ERR1 = ERRO;
  1294. 'SI' (ERR1 > XXTOL) ;
  1295. CHAI1 = CHAI 'ERREUR NUMERO 33: ' ERR1 'POUR ' XXTOL;
  1296. 'MESS' CHAI1;
  1297. 'ERREUR' 5 ;
  1298. 'FINSI' ;
  1299.  
  1300.  
  1301. ***************** RYO2 ************************************
  1302.  
  1303. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
  1304. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1305. DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
  1306.  
  1307. RN1 = 'COPIER' RN ;
  1308. GN1 = 'COPIER' GN ;
  1309. RETN1 = 'COPIER' RETN ;
  1310. RYN1 = RYN '+' DR ;
  1311.  
  1312. VN1 PN1 TTEMP YN1 GAMM =
  1313. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1314.  
  1315. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1316. PGAS LISTC LISTP
  1317. RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
  1318.  
  1319. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1320. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1321. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
  1322.  
  1323. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1324.  
  1325. XXTOL = ERROM;
  1326. ERR1 = ERRO;
  1327. 'SI' (ERR1 > XXTOL) ;
  1328. CHAI1 = CHAI 'ERREUR NUMERO 34: ' ERR1 'POUR ' XXTOL;
  1329. 'MESS' CHAI1;
  1330. 'ERREUR' 5 ;
  1331. 'FINSI' ;
  1332.  
  1333.  
  1334. ***************** RYH2O ************************************
  1335.  
  1336. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
  1337. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1338. DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
  1339.  
  1340. RN1 = 'COPIER' RN ;
  1341. GN1 = 'COPIER' GN ;
  1342. RETN1 = 'COPIER' RETN ;
  1343. RYN1 = RYN '+' DR ;
  1344.  
  1345. VN1 PN1 TTEMP YN1 GAMM =
  1346. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1347.  
  1348. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1349. PGAS LISTC LISTP
  1350. RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
  1351.  
  1352. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1353. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1354. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
  1355.  
  1356. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1357.  
  1358. XXTOL = ERROM;
  1359. ERR1 = ERRO;
  1360. 'SI' (ERR1 > XXTOL) ;
  1361. CHAI1 = CHAI 'ERREUR NUMERO 35: ' ERR1 'POUR ' XXTOL;
  1362. 'MESS' CHAI1;
  1363. 'ERREUR' 5 ;
  1364. 'FINSI' ;
  1365.  
  1366.  
  1367. ************************************************************************
  1368. ************************************************************************
  1369. *************************** Pressure at the outlet ********************
  1370. ************************************************************************
  1371. ************************************************************************
  1372. *
  1373. * We solve a Riemann probleme between the state in the center and
  1374. * the same state but the pressure equal to the pressure at the outlet
  1375. * We use AUSMPLUS
  1376. *
  1377. *
  1378. ***** Calcul du residu dans le cas particulier PF = PC
  1379.  
  1380. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1381. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1382. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' -15. 'UX' 1. ;
  1383. YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
  1384. 'H2O' 0.1 ;
  1385.  
  1386. GN = RN '*' VN ;
  1387. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 1 'PN' 2.1 ;
  1388. *----------------------------------------------------------
  1389. YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
  1390. YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
  1391. YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
  1392. YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
  1393. YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
  1394. 'H2O' YH2O 'N2' YN2 ;
  1395. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  1396. CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  1397. CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  1398. GAMN = CPTOT '/' CVTOT ;
  1399. GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
  1400. *----------------------------------------------------------
  1401. LMOT1 = (PGAS . 'ESPEULE') ;
  1402. LMOT2 = 'MOTS' 'SCAL' 'SCAL' 'SCAL' ;
  1403. LMOT3 = (PGAS . 'ESPEULE') ;
  1404.  
  1405. RYN = YN '*' RN LMOT1 LMOT2 LMOT3 ;
  1406.  
  1407. ************************************************
  1408.  
  1409. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  1410. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
  1411.  
  1412. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
  1413. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
  1414.  
  1415. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1416. PGAS LISTC LISTP
  1417. RN VN PN YN CHPLI 'OUTP' ;
  1418.  
  1419. *-----------------------------------------------
  1420. 'SI' ('EGA' TYEL 'QUA8') ;
  1421. SURSVOL = -1. '/' ((2 '**' 0.5) * 0.5) ;
  1422. 'SINON' ;
  1423. SURSVOL = -2. '/' ((2 '**' 0.5) * 0.5) ;
  1424. 'FINSI' ;
  1425. UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
  1426. UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
  1427. UNORM = (UXC '-' UYC) '/' (2. '**' 0.5) ;
  1428. UT = (UXC '+' UYC) '/' (2. '**' 0.5) ;
  1429. PC = 'EXTRAIRE' PN PCEN 'SCAL' ;
  1430. RC = 'EXTRAIRE' RN PCEN 'SCAL' ;
  1431.  
  1432. FMASS = RC '*' UNORM ;
  1433. FUN = (RC '*' UNORM '*' UNORM) '+' PC ;
  1434. FUT = (RC '*' UNORM '*' UT) ;
  1435. FUX = (FUN '+' FUT) '/' (2. '**' 0.5) ;
  1436. FUY = (FUT '-' FUN) '/' (2. '**' 0.5) ;
  1437. FEN = (UNORM '*' ((GAMC/(GAMC-1)) '*' PC)) '+'
  1438. (0.5 '*' UNORM '*' RC '*' ((UNORM '*' UNORM) '+' (UT '*' UT)))
  1439. ;
  1440. FYH2 = YH2 '*' RC '*' UNORM ;
  1441. FYO2 = YO2 '*' RC '*' UNORM ;
  1442. FYH2O = YH2O '*' RC '*' UNORM ;
  1443.  
  1444. RCHRES1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' PCEN) 7
  1445. 'VAR1' (SURSVOL '*' FMASS) 'VAR2' (FUX '*' SURSVOL) 'VAR3' (FUY '*'
  1446. SURSVOL) 'VAR4' (FEN '*' SURSVOL) 'VAR5' (FYH2 '*' SURSVOL)
  1447. 'VAR6' (FYO2 '*' SURSVOL) 'VAR7' (FYH2O '*' SURSVOL) ;
  1448.  
  1449. ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
  1450.  
  1451. XXTOL = 1.D-12;
  1452. ERR1 = ERRO;
  1453. 'SI' (ERR1 > XXTOL) ;
  1454. CHAI1 = CHAI 'ERREUR NUMERO 36: ' ERR1 'POUR ' XXTOL;
  1455. 'MESS' CHAI1;
  1456. 'ERREUR' 5 ;
  1457. 'FINSI' ;
  1458.  
  1459.  
  1460. *------------------------------------------------------------
  1461. ***** Jacobian with respect to the primitive variables ******
  1462. *------------------------------------------------------------
  1463. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 1 'PN' 4.1 ;
  1464. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1465. PGAS LISTC LISTP
  1466. RN VN PN YN CHPLI 'OUTP' ;
  1467.  
  1468. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
  1469. PGAS LISTC LISTP
  1470. RN VN PN YN CHPLI 'OUTP' ;
  1471.  
  1472.  
  1473. ************* RN *********** RN ***********************
  1474.  
  1475. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1476. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1477. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1478.  
  1479. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1480. PGAS LISTC LISTP
  1481. (RN '+' DR) VN PN YN CHPLI 'OUTP' ;
  1482.  
  1483. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1484. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
  1485. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  1486.  
  1487. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1488.  
  1489. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1490.  
  1491. XXTOL = ERROM;
  1492. ERR1 = ERRO;
  1493. 'SI' (ERR1 > XXTOL) ;
  1494. CHAI1 = CHAI 'ERREUR NUMERO 37: ' ERR1 'POUR ' XXTOL;
  1495. 'MESS' CHAI1;
  1496. 'ERREUR' 5 ;
  1497. 'FINSI' ;
  1498.  
  1499.  
  1500. ************* UX *********** UX ***********************
  1501.  
  1502. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  1503. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1504. DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
  1505.  
  1506. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1507. PGAS LISTC LISTP
  1508. RN (VN '+' DR) PN YN CHPLI 'OUTP' ;
  1509.  
  1510. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1511. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1512. 'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  1513.  
  1514. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1515.  
  1516. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1517.  
  1518. XXTOL = ERROM;
  1519. ERR1 = ERRO;
  1520. 'SI' (ERR1 > XXTOL) ;
  1521. CHAI1 = CHAI 'ERREUR NUMERO 38: ' ERR1 'POUR ' XXTOL;
  1522. 'MESS' CHAI1;
  1523. 'ERREUR' 5 ;
  1524. 'FINSI' ;
  1525.  
  1526.  
  1527. ************* UY *********** UY ***********************
  1528.  
  1529. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  1530. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1531. DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
  1532.  
  1533. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1534. PGAS LISTC LISTP
  1535. RN (VN '+' DR) PN YN CHPLI 'OUTP' ;
  1536.  
  1537. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1538. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1539. 'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  1540.  
  1541. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1542.  
  1543. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1544.  
  1545. XXTOL = ERROM;
  1546. ERR1 = ERRO;
  1547. 'SI' (ERR1 > XXTOL) ;
  1548. CHAI1 = CHAI 'ERREUR NUMERO 39: ' ERR1 'POUR ' XXTOL;
  1549. 'MESS' CHAI1;
  1550. 'ERREUR' 5 ;
  1551. 'FINSI' ;
  1552.  
  1553.  
  1554. ************* PN *********** PN ***********************
  1555.  
  1556. DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1557. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1558. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1559.  
  1560. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1561. PGAS LISTC LISTP
  1562. RN VN (PN '+' DR) YN CHPLI 'OUTP' ;
  1563.  
  1564. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1565. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1566. 'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  1567.  
  1568. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1569.  
  1570. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1571.  
  1572. XXTOL = ERROM;
  1573. ERR1 = ERRO;
  1574. 'SI' (ERR1 > XXTOL) ;
  1575. CHAI1 = CHAI 'ERREUR NUMERO 40: ' ERR1 'POUR ' XXTOL;
  1576. 'MESS' CHAI1;
  1577. 'ERREUR' 5 ;
  1578. 'FINSI' ;
  1579.  
  1580.  
  1581. ************* H2 *********** H2 ***********************
  1582.  
  1583. DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
  1584. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1585. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  1586.  
  1587. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1588. PGAS LISTC LISTP
  1589. RN VN PN (YN '+' DR) CHPLI 'OUTP' ;
  1590.  
  1591. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1592. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1593. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
  1594.  
  1595. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1596.  
  1597. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1598.  
  1599. XXTOL = ERROM;
  1600. ERR1 = ERRO;
  1601. 'SI' (ERR1 > XXTOL) ;
  1602. CHAI1 = CHAI 'ERREUR NUMERO 41: ' ERR1 'POUR ' XXTOL;
  1603. 'MESS' CHAI1;
  1604. 'ERREUR' 5 ;
  1605. 'FINSI' ;
  1606.  
  1607.  
  1608. ************* O2 *********** O2 ***********************
  1609.  
  1610. DRSCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
  1611. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1612. DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
  1613.  
  1614. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1615. PGAS LISTC LISTP
  1616. RN VN PN (YN '+' DR) CHPLI 'OUTP' ;
  1617.  
  1618. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1619. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1620. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
  1621.  
  1622. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1623.  
  1624. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1625.  
  1626. XXTOL = ERROM;
  1627. ERR1 = ERRO;
  1628. 'SI' (ERR1 > XXTOL) ;
  1629. CHAI1 = CHAI 'ERREUR NUMERO 42: ' ERR1 'POUR ' XXTOL;
  1630. 'MESS' CHAI1;
  1631. 'ERREUR' 5 ;
  1632. 'FINSI' ;
  1633.  
  1634.  
  1635.  
  1636. ************* H2O *********** H2O ***********************
  1637.  
  1638. DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
  1639. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1640. DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
  1641.  
  1642. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1643. PGAS LISTC LISTP
  1644. RN VN PN (YN '+' DR) CHPLI 'OUTP' ;
  1645.  
  1646. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1647. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  1648. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
  1649.  
  1650. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1651.  
  1652. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1653.  
  1654. XXTOL = ERROM;
  1655. ERR1 = ERRO;
  1656. 'SI' (ERR1 > XXTOL) ;
  1657. CHAI1 = CHAI 'ERREUR NUMERO 43: ' ERR1 'POUR ' XXTOL;
  1658. 'MESS' CHAI1;
  1659. 'ERREUR' 5 ;
  1660. 'FINSI' ;
  1661.  
  1662.  
  1663. *----------------------------------------------------------
  1664. **** Jacobian with respect to conservative variables
  1665. *---------------------------------------------------------
  1666.  
  1667. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
  1668. PGAS LISTC LISTP
  1669. RN VN PN YN CHPLI 'OUTP' ;
  1670.  
  1671. *************** RN **************** RN *******************
  1672.  
  1673. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1674. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1675. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1676.  
  1677. RN1 = RN '+' DR ;
  1678. GN1 = 'COPIER' GN ;
  1679. RETN1 = 'COPIER' RETN ;
  1680. RYN1 = 'COPIER' RYN ;
  1681.  
  1682. VN1 PN1 TTEMP YN1 GAMM =
  1683. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1684.  
  1685. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1686. PGAS LISTC LISTP
  1687. RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
  1688.  
  1689. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1690. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
  1691. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1692.  
  1693. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1694.  
  1695. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
  1696.  
  1697. XXTOL = ERROM;
  1698. ERR1 = ERRO;
  1699. 'SI' (ERR1 > XXTOL) ;
  1700. CHAI1 = CHAI 'ERREUR NUMERO 44: ' ERR1 'POUR ' XXTOL;
  1701. 'MESS' CHAI1;
  1702. 'ERREUR' 5 ;
  1703. 'FINSI' ;
  1704.  
  1705.  
  1706. ***************** RUX ************************************
  1707.  
  1708. DRSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  1709. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1710. DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
  1711.  
  1712. RN1 = 'COPIER' RN ;
  1713. GN1 = GN '+' DR ;
  1714. RETN1 = 'COPIER' RETN ;
  1715. RYN1 = 'COPIER' RYN ;
  1716.  
  1717. VN1 PN1 TTEMP YN1 GAMM =
  1718. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1719.  
  1720. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1721. PGAS LISTC LISTP
  1722. RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
  1723.  
  1724. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1725. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1726. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1727.  
  1728. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1729.  
  1730. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1731.  
  1732. XXTOL = ERROM;
  1733. ERR1 = ERRO;
  1734. 'SI' (ERR1 > XXTOL) ;
  1735. CHAI1 = CHAI 'ERREUR NUMERO 45: ' ERR1 'POUR ' XXTOL;
  1736. 'MESS' CHAI1;
  1737. 'ERREUR' 5 ;
  1738. 'FINSI' ;
  1739.  
  1740.  
  1741. ***************** RUY ************************************
  1742.  
  1743. DRSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  1744. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1745. DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
  1746.  
  1747. RN1 = 'COPIER' RN ;
  1748. GN1 = GN '+' DR ;
  1749. RETN1 = 'COPIER' RETN ;
  1750. RYN1 = 'COPIER' RYN ;
  1751.  
  1752. VN1 PN1 TTEMP YN1 GAMM =
  1753. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1754.  
  1755. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1756. PGAS LISTC LISTP
  1757. RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
  1758.  
  1759. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1760. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1761. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1762.  
  1763. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1764.  
  1765. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1766.  
  1767. XXTOL = ERROM;
  1768. ERR1 = ERRO;
  1769. 'SI' (ERR1 > XXTOL) ;
  1770. CHAI1 = CHAI 'ERREUR NUMERO 46: ' ERR1 'POUR ' XXTOL;
  1771. 'MESS' CHAI1;
  1772. 'ERREUR' 5 ;
  1773. 'FINSI' ;
  1774.  
  1775.  
  1776. ***************** RETN ************************************
  1777.  
  1778. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1779. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1780. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1781.  
  1782. RN1 = 'COPIER' RN ;
  1783. GN1 = 'COPIER' GN ;
  1784. RETN1 = RETN '+' DR ;
  1785. RYN1 = 'COPIER' RYN ;
  1786.  
  1787. VN1 PN1 TTEMP YN1 GAMM =
  1788. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1789.  
  1790. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1791. PGAS LISTC LISTP
  1792. RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
  1793.  
  1794. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1795. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1796. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1797.  
  1798. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1799.  
  1800. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1801.  
  1802. XXTOL = ERROM;
  1803. ERR1 = ERRO;
  1804. 'SI' (ERR1 > XXTOL) ;
  1805. CHAI1 = CHAI 'ERREUR NUMERO 47: ' ERR1 'POUR ' XXTOL;
  1806. 'MESS' CHAI1;
  1807. 'ERREUR' 5 ;
  1808. 'FINSI' ;
  1809.  
  1810.  
  1811. ***************** RYH2 ************************************
  1812.  
  1813. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
  1814. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1815. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  1816.  
  1817. RN1 = 'COPIER' RN ;
  1818. GN1 = 'COPIER' GN ;
  1819. RETN1 = 'COPIER' RETN ;
  1820. RYN1 = RYN '+' DR ;
  1821.  
  1822. VN1 PN1 TTEMP YN1 GAMM =
  1823. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1824.  
  1825. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1826. PGAS LISTC LISTP
  1827. RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
  1828.  
  1829. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1830. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1831. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  1832.  
  1833. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1834.  
  1835. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1836.  
  1837. XXTOL = ERROM;
  1838. ERR1 = ERRO;
  1839. 'SI' (ERR1 > XXTOL) ;
  1840. CHAI1 = CHAI 'ERREUR NUMERO 48: ' ERR1 'POUR ' XXTOL;
  1841. 'MESS' CHAI1;
  1842. 'ERREUR' 5 ;
  1843. 'FINSI' ;
  1844.  
  1845.  
  1846. ***************** RYO2 ************************************
  1847.  
  1848. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
  1849. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1850. DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
  1851.  
  1852. RN1 = 'COPIER' RN ;
  1853. GN1 = 'COPIER' GN ;
  1854. RETN1 = 'COPIER' RETN ;
  1855. RYN1 = RYN '+' DR ;
  1856.  
  1857. VN1 PN1 TTEMP YN1 GAMM =
  1858. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1859.  
  1860. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1861. PGAS LISTC LISTP
  1862. RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
  1863.  
  1864. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1865. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1866. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
  1867.  
  1868. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1869.  
  1870. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1871.  
  1872. XXTOL = ERROM;
  1873. ERR1 = ERRO;
  1874. 'SI' (ERR1 > XXTOL) ;
  1875. CHAI1 = CHAI 'ERREUR NUMERO 49: ' ERR1 'POUR ' XXTOL;
  1876. 'MESS' CHAI1;
  1877. 'ERREUR' 5 ;
  1878. 'FINSI' ;
  1879.  
  1880.  
  1881. ***************** RYH2O ************************************
  1882.  
  1883. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
  1884. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1885. DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
  1886.  
  1887. RN1 = 'COPIER' RN ;
  1888. GN1 = 'COPIER' GN ;
  1889. RETN1 = 'COPIER' RETN ;
  1890. RYN1 = RYN '+' DR ;
  1891.  
  1892. VN1 PN1 TTEMP YN1 GAMM =
  1893. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  1894.  
  1895. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1896. PGAS LISTC LISTP
  1897. RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
  1898.  
  1899. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1900. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  1901. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
  1902.  
  1903. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  1904.  
  1905. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  1906.  
  1907. XXTOL = ERROM;
  1908. ERR1 = ERRO;
  1909. 'SI' (ERR1 > XXTOL) ;
  1910. CHAI1 = CHAI 'ERREUR NUMERO 50: ' ERR1 'POUR ' XXTOL;
  1911. 'MESS' CHAI1;
  1912. 'ERREUR' 5 ;
  1913. 'FINSI' ;
  1914.  
  1915.  
  1916. *******************************************************************
  1917. *******************************************************************
  1918. *************************** INSU *********************************
  1919. *******************************************************************
  1920. *******************************************************************
  1921. * Total enthalpy, entropy, tangential velocity at the inlet
  1922. * Normal velocity inside
  1923. *
  1924. * We solve a Riemann probleme between the state at the interface thus
  1925. * reconstructed and the state at the center
  1926. *
  1927. ***** Calcul du residu dans le cas particulier UF = UC
  1928. *
  1929.  
  1930. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1931. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1932. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 15. 'UX' -15. ;
  1933. YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
  1934. 'H2O' 0.1 ;
  1935.  
  1936. GN = RN '*' VN ;
  1937. *----------------------------------------------------------
  1938. YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
  1939. YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
  1940. YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
  1941. YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
  1942. YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
  1943. 'H2O' YH2O 'N2' YN2 ;
  1944. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  1945. CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  1946. CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  1947. GAMN = CPTOT '/' CVTOT ;
  1948. GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
  1949. *----------------------------------------------------------
  1950. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  1951. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
  1952. RHTN = RETN '+' PN ;
  1953. HTN = RHTN '/' RN ;
  1954. UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
  1955. UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
  1956. UTC = -1 '*' (UXC '+' UYC) '/' (2. '**' 0.5) ;
  1957.  
  1958. SN = PN '/' (RN '**' ('MAXI' GAMN)) ;
  1959.  
  1960. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'HT' ('MAXIMUM' HTN)
  1961. 'S' ('MAXIMUM' SN) 'H2' 0.3 'O2' 0.5 'H2O' 0.1;
  1962.  
  1963. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
  1964. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
  1965.  
  1966. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  1967. PGAS LISTC LISTP
  1968. RN VN PN YN CHPLI 'INSU' ;
  1969.  
  1970. *-------------------------------------------------------------
  1971.  
  1972. 'SI' ('EGA' TYEL 'QUA8') ;
  1973. SURSVOL = -1. '/' ((2 '**' 0.5) * 0.5) ;
  1974. 'SINON' ;
  1975. SURSVOL = -2. '/' ((2 '**' 0.5) * 0.5) ;
  1976. 'FINSI' ;
  1977.  
  1978. UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
  1979. UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
  1980. UNORM = (UXC '-' UYC) '/' (2. '**' 0.5) ;
  1981. UT = (UXC '+' UYC) '/' (2. '**' 0.5) ;
  1982. PC = 'EXTRAIRE' PN PCEN 'SCAL' ;
  1983. RC = 'EXTRAIRE' RN PCEN 'SCAL' ;
  1984. GAMC = 'EXTRAIRE' GAMN PCEN 'SCAL' ;
  1985.  
  1986. RCHLIM2 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' RC 'PN' PC 'UX'
  1987. UXC 'UY' UYC 'H2' YH2 'O2' YO2 'H2O' YH2O ;
  1988.  
  1989. ERRO = 'MAXIMUM' (RCHLIM '-' RCHLIM2) 'ABS' ;
  1990.  
  1991. XXTOL = 1.D-6;
  1992. ERR1 = ERRO;
  1993. 'SI' (ERR1 > XXTOL) ;
  1994. CHAI1 = CHAI 'ERREUR NUMERO 51: ' ERR1 'POUR ' XXTOL;
  1995. 'MESS' CHAI1;
  1996. 'ERREUR' 5 ;
  1997. 'FINSI' ;
  1998.  
  1999.  
  2000. FMASS = RC '*' UNORM ;
  2001. FUN = (RC '*' UNORM '*' UNORM) '+' PC ;
  2002. FUT = (RC '*' UNORM '*' UT) ;
  2003. FUX = (FUN '+' FUT) '/' (2. '**' 0.5) ;
  2004. FUY = (FUT '-' FUN) '/' (2. '**' 0.5) ;
  2005. FEN = (UNORM '*' ((GAMC/(GAMC-1)) '*' PC)) '+'
  2006. (0.5 '*' UNORM '*' RC '*' ((UNORM '*' UNORM) '+'
  2007. (UT '*' UT))) ;
  2008. FH2 = RC '*' UNORM '*' YH2 ;
  2009. FO2 = RC '*' UNORM '*' YO2 ;
  2010. FHO = RC '*' UNORM '*' YH2O ;
  2011. RCHRES1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' PCEN) 7
  2012. 'VAR1' (SURSVOL '*' FMASS) 'VAR2' (FUX '*' SURSVOL)
  2013. 'VAR3' (FUY '*' SURSVOL) 'VAR4' (FEN '*' SURSVOL)
  2014. 'VAR5' (FH2 '*' SURSVOL) 'VAR6' (FO2 '*' SURSVOL)
  2015. 'VAR7' (FHO '*' SURSVOL) ;
  2016.  
  2017. ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
  2018.  
  2019.  
  2020. XXTOL = 1.D-6;
  2021. ERR1 = ERRO;
  2022. 'SI' (ERR1 > XXTOL) ;
  2023. CHAI1 = CHAI 'ERREUR NUMERO 52: ' ERR1 'POUR ' XXTOL;
  2024. 'MESS' CHAI1;
  2025. 'ERREUR' 5 ;
  2026. 'FINSI' ;
  2027.  
  2028. *------------------------------------------------------------
  2029. ***** Jacobian with respect to the primitive variables ******
  2030. *------------------------------------------------------------
  2031. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'HT' 250.0
  2032. 'S' 2.0 'H2' 0.1 'O2' 0.7 'H2O' 0.2;
  2033.  
  2034. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2035. PGAS LISTC LISTP
  2036. RN VN PN YN CHPLI 'INSU' ;
  2037.  
  2038. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
  2039. PGAS LISTC LISTP
  2040. RN VN PN YN CHPLI 'INSU' ;
  2041.  
  2042. ************* RN *********** RN ***********************
  2043.  
  2044. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  2045. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2046. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2047.  
  2048. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2049. PGAS LISTC LISTP
  2050. (RN '+' DR) VN PN YN CHPLI 'INSU' ;
  2051.  
  2052. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2053. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
  2054. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2055.  
  2056. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2057.  
  2058. 'SI' (MMM < ERROM) ;
  2059. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2060. 'SINON' ;
  2061. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2062. 'FINSI' ;
  2063.  
  2064.  
  2065. XXTOL = ERROM;
  2066. ERR1 = ERRO;
  2067. 'SI' (ERR1 > XXTOL) ;
  2068. CHAI1 = CHAI 'ERREUR NUMERO 53: ' ERR1 'POUR ' XXTOL;
  2069. 'MESS' CHAI1;
  2070. 'ERREUR' 5 ;
  2071. 'FINSI' ;
  2072.  
  2073.  
  2074. ************* UX *********** UX ***********************
  2075.  
  2076. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  2077. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2078. DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
  2079.  
  2080. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2081. PGAS LISTC LISTP
  2082. RN (VN '+' DR) PN YN CHPLI 'INSU' ;
  2083.  
  2084. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2085. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2086. 'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2087.  
  2088. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2089.  
  2090.  
  2091. 'SI' (MMM < ERROM) ;
  2092. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2093. 'SINON' ;
  2094. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2095. 'FINSI' ;
  2096.  
  2097. XXTOL = ERROM;
  2098. ERR1 = ERRO;
  2099. 'SI' (ERR1 > XXTOL) ;
  2100. CHAI1 = CHAI 'ERREUR NUMERO 54: ' ERR1 'POUR ' XXTOL;
  2101. 'MESS' CHAI1;
  2102. 'ERREUR' 5 ;
  2103. 'FINSI' ;
  2104.  
  2105.  
  2106. ************* UY *********** UY ***********************
  2107.  
  2108. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  2109. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2110. DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
  2111.  
  2112. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2113. PGAS LISTC LISTP
  2114. RN (VN '+' DR) PN YN CHPLI 'INSU' ;
  2115.  
  2116. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2117. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2118. 'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2119.  
  2120. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2121.  
  2122. 'SI' (MMM < ERROM) ;
  2123. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2124. 'SINON' ;
  2125. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2126. 'FINSI' ;
  2127.  
  2128. XXTOL = ERROM;
  2129. ERR1 = ERRO;
  2130. 'SI' (ERR1 > XXTOL) ;
  2131. CHAI1 = CHAI 'ERREUR NUMERO 55: ' ERR1 'POUR ' XXTOL;
  2132. 'MESS' CHAI1;
  2133. 'ERREUR' 5 ;
  2134. 'FINSI' ;
  2135.  
  2136.  
  2137. ************* PN *********** PN ***********************
  2138.  
  2139. DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  2140. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2141. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2142.  
  2143. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2144. PGAS LISTC LISTP
  2145. RN VN (PN '+' DR) YN CHPLI 'INSU' ;
  2146.  
  2147. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2148. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2149. 'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2150.  
  2151. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2152.  
  2153. 'SI' (MMM < ERROM) ;
  2154. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2155. 'SINON' ;
  2156. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2157. 'FINSI' ;
  2158.  
  2159. XXTOL = ERROM;
  2160. ERR1 = ERRO;
  2161. 'SI' (ERR1 > XXTOL) ;
  2162. CHAI1 = CHAI 'ERREUR NUMERO 56: ' ERR1 'POUR ' XXTOL;
  2163. 'MESS' CHAI1;
  2164. 'ERREUR' 5 ;
  2165. 'FINSI' ;
  2166.  
  2167.  
  2168. ************* H2 *********** H2 ***********************
  2169.  
  2170. DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
  2171. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2172. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  2173.  
  2174. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2175. PGAS LISTC LISTP
  2176. RN VN PN (YN '+' DR) CHPLI 'INSU' ;
  2177.  
  2178. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2179. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2180. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
  2181.  
  2182. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2183.  
  2184. 'SI' (MMM < ERROM) ;
  2185. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2186. 'SINON' ;
  2187. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2188. 'FINSI' ;
  2189.  
  2190. XXTOL = ERROM;
  2191. ERR1 = ERRO;
  2192. 'SI' (ERR1 > XXTOL) ;
  2193. CHAI1 = CHAI 'ERREUR NUMERO 57: ' ERR1 'POUR ' XXTOL;
  2194. 'MESS' CHAI1;
  2195. 'ERREUR' 5 ;
  2196. 'FINSI' ;
  2197.  
  2198. ************* O2 *********** O2 ***********************
  2199.  
  2200. DRSCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
  2201. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2202. DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
  2203.  
  2204. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2205. PGAS LISTC LISTP
  2206. RN VN PN (YN '+' DR) CHPLI 'INSU' ;
  2207.  
  2208. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2209. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2210. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
  2211.  
  2212. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2213.  
  2214. 'SI' (MMM < ERROM) ;
  2215. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2216. 'SINON' ;
  2217. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2218. 'FINSI' ;
  2219.  
  2220. XXTOL = ERROM;
  2221. ERR1 = ERRO;
  2222. 'SI' (ERR1 > XXTOL) ;
  2223. CHAI1 = CHAI 'ERREUR NUMERO 58: ' ERR1 'POUR ' XXTOL;
  2224. 'MESS' CHAI1;
  2225. 'ERREUR' 5 ;
  2226. 'FINSI' ;
  2227.  
  2228.  
  2229.  
  2230. ************* H2O *********** H2O ***********************
  2231.  
  2232. DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
  2233. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2234. DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
  2235.  
  2236. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2237. PGAS LISTC LISTP
  2238. RN VN PN (YN '+' DR) CHPLI 'INSU' ;
  2239.  
  2240. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2241. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2242. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
  2243.  
  2244. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2245.  
  2246. 'SI' (MMM < ERROM) ;
  2247. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2248. 'SINON' ;
  2249. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2250. 'FINSI' ;
  2251.  
  2252.  
  2253. XXTOL = ERROM;
  2254. ERR1 = ERRO;
  2255. 'SI' (ERR1 > XXTOL) ;
  2256. CHAI1 = CHAI 'ERREUR NUMERO 59: ' ERR1 'POUR ' XXTOL;
  2257. 'MESS' CHAI1;
  2258. 'ERREUR' 5 ;
  2259. 'FINSI' ;
  2260.  
  2261.  
  2262. *----------------------------------------------------------
  2263. **** Jacobian with respect to conservative variables
  2264. *---------------------------------------------------------
  2265.  
  2266. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
  2267. PGAS LISTC LISTP
  2268. RN VN PN YN CHPLI 'INSU' ;
  2269.  
  2270.  
  2271. *************** RN **************** RN *******************
  2272.  
  2273. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  2274. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2275. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2276.  
  2277. RN1 = RN '+' DR ;
  2278. GN1 = 'COPIER' GN ;
  2279. RETN1 = 'COPIER' RETN ;
  2280. RYN1 = 'COPIER' RYN ;
  2281.  
  2282. VN1 PN1 TTEMP YN1 GAMM =
  2283. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  2284.  
  2285. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2286. PGAS LISTC LISTP
  2287. RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
  2288.  
  2289. NNN = 'MAXIMUM' RCHRES 'ABS' ;
  2290.  
  2291. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2292. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
  2293. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  2294.  
  2295. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2296.  
  2297. 'SI' (MMM < ERROM) ;
  2298. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2299. 'SINON' ;
  2300. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2301. 'FINSI' ;
  2302.  
  2303. * 'OPTION' DONN 5 ;
  2304.  
  2305. XXTOL = ERROM;
  2306. ERR1 = ERRO;
  2307. 'SI' (ERR1 > XXTOL) ;
  2308. CHAI1 = CHAI 'ERREUR NUMERO 60: ' ERR1 'POUR ' XXTOL;
  2309. 'MESS' CHAI1;
  2310. 'ERREUR' 5 ;
  2311. 'FINSI' ;
  2312.  
  2313.  
  2314. ***************** RUX ************************************
  2315.  
  2316. DRSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  2317. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2318. DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
  2319.  
  2320. RN1 = 'COPIER' RN ;
  2321. GN1 = GN '+' DR ;
  2322. RETN1 = 'COPIER' RETN ;
  2323. RYN1 = 'COPIER' RYN ;
  2324.  
  2325. VN1 PN1 TTEMP YN1 GAMM =
  2326. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  2327.  
  2328. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2329. PGAS LISTC LISTP
  2330. RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
  2331.  
  2332. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2333. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  2334. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  2335.  
  2336. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2337.  
  2338. 'SI' (MMM < ERROM) ;
  2339. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2340. 'SINON' ;
  2341. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2342. 'FINSI' ;
  2343.  
  2344. * 'OPTION' DONN 5 ;
  2345.  
  2346. XXTOL = ERROM;
  2347. ERR1 = ERRO;
  2348. 'SI' (ERR1 > XXTOL) ;
  2349. CHAI1 = CHAI 'ERREUR NUMERO 61: ' ERR1 'POUR ' XXTOL;
  2350. 'MESS' CHAI1;
  2351. 'ERREUR' 5 ;
  2352. 'FINSI' ;
  2353.  
  2354.  
  2355. ***************** RUY ************************************
  2356.  
  2357. DRSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  2358. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2359. DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
  2360.  
  2361. RN1 = 'COPIER' RN ;
  2362. GN1 = GN '+' DR ;
  2363. RETN1 = 'COPIER' RETN ;
  2364. RYN1 = 'COPIER' RYN ;
  2365.  
  2366. VN1 PN1 TTEMP YN1 GAMM =
  2367. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  2368.  
  2369. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2370. PGAS LISTC LISTP
  2371. RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
  2372.  
  2373. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2374. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  2375. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  2376.  
  2377. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2378.  
  2379. 'SI' (MMM < ERROM) ;
  2380. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2381. 'SINON' ;
  2382. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2383. 'FINSI' ;
  2384.  
  2385. XXTOL = ERROM;
  2386. ERR1 = ERRO;
  2387. 'SI' (ERR1 > XXTOL) ;
  2388. CHAI1 = CHAI 'ERREUR NUMERO 62: ' ERR1 'POUR ' XXTOL;
  2389. 'MESS' CHAI1;
  2390. 'ERREUR' 5 ;
  2391. 'FINSI' ;
  2392.  
  2393.  
  2394. ***************** RETN ************************************
  2395.  
  2396. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  2397. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2398. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2399.  
  2400. RN1 = 'COPIER' RN ;
  2401. GN1 = 'COPIER' GN ;
  2402. RETN1 = RETN '+' DR ;
  2403. RYN1 = 'COPIER' RYN ;
  2404.  
  2405. VN1 PN1 TTEMP YN1 GAMM =
  2406. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  2407.  
  2408. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2409. PGAS LISTC LISTP
  2410. RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
  2411.  
  2412. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2413. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  2414. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  2415.  
  2416. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2417.  
  2418. 'SI' (MMM < ERROM) ;
  2419. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2420. 'SINON' ;
  2421. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2422. 'FINSI' ;
  2423.  
  2424. XXTOL = ERROM;
  2425. ERR1 = ERRO;
  2426. 'SI' (ERR1 > XXTOL) ;
  2427. CHAI1 = CHAI 'ERREUR NUMERO 63: ' ERR1 'POUR ' XXTOL;
  2428. 'MESS' CHAI1;
  2429. 'ERREUR' 5 ;
  2430. 'FINSI' ;
  2431.  
  2432. ***************** RYH2 ************************************
  2433.  
  2434. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
  2435. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2436. DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
  2437.  
  2438. RN1 = 'COPIER' RN ;
  2439. GN1 = 'COPIER' GN ;
  2440. RETN1 = 'COPIER' RETN ;
  2441. RYN1 = RYN '+' DR ;
  2442.  
  2443. VN1 PN1 TTEMP YN1 GAMM =
  2444. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  2445.  
  2446. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2447. PGAS LISTC LISTP
  2448. RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
  2449.  
  2450. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2451. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  2452. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
  2453.  
  2454. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2455.  
  2456. 'SI' (MMM < ERROM) ;
  2457. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2458. 'SINON' ;
  2459. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2460. 'FINSI' ;
  2461.  
  2462. XXTOL = ERROM;
  2463. ERR1 = ERRO;
  2464. 'SI' (ERR1 > XXTOL) ;
  2465. CHAI1 = CHAI 'ERREUR NUMERO 64: ' ERR1 'POUR ' XXTOL;
  2466. 'MESS' CHAI1;
  2467. 'ERREUR' 5 ;
  2468. 'FINSI' ;
  2469.  
  2470.  
  2471. ***************** RYO2 ************************************
  2472.  
  2473. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
  2474. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2475. DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
  2476.  
  2477. RN1 = 'COPIER' RN ;
  2478. GN1 = 'COPIER' GN ;
  2479. RETN1 = 'COPIER' RETN ;
  2480. RYN1 = RYN '+' DR ;
  2481.  
  2482. VN1 PN1 TTEMP YN1 GAMM =
  2483. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  2484.  
  2485. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2486. PGAS LISTC LISTP
  2487. RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
  2488.  
  2489. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2490. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  2491. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
  2492.  
  2493. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2494.  
  2495. 'SI' (MMM < ERROM) ;
  2496. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2497. 'SINON' ;
  2498. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2499. 'FINSI' ;
  2500.  
  2501. XXTOL = ERROM;
  2502. ERR1 = ERRO;
  2503. 'SI' (ERR1 > XXTOL) ;
  2504. CHAI1 = CHAI 'ERREUR NUMERO 65: ' ERR1 'POUR ' XXTOL;
  2505. 'MESS' CHAI1;
  2506. 'ERREUR' 5 ;
  2507. 'FINSI' ;
  2508.  
  2509.  
  2510. ***************** RYH2O ************************************
  2511.  
  2512. DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
  2513. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2514. DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
  2515.  
  2516. RN1 = 'COPIER' RN ;
  2517. GN1 = 'COPIER' GN ;
  2518. RETN1 = 'COPIER' RETN ;
  2519. RYN1 = RYN '+' DR ;
  2520.  
  2521. VN1 PN1 TTEMP YN1 GAMM =
  2522. 'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
  2523.  
  2524. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2525. PGAS LISTC LISTP
  2526. RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
  2527.  
  2528. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2529. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
  2530. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
  2531.  
  2532. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2533.  
  2534. 'SI' (MMM < ERROM) ;
  2535. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2536. 'SINON' ;
  2537. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2538. 'FINSI' ;
  2539.  
  2540. XXTOL = ERROM;
  2541. ERR1 = ERRO;
  2542. 'SI' (ERR1 > XXTOL) ;
  2543. CHAI1 = CHAI 'ERREUR NUMERO 66: ' ERR1 'POUR ' XXTOL;
  2544. 'MESS' CHAI1;
  2545. 'ERREUR' 5 ;
  2546. 'FINSI' ;
  2547.  
  2548.  
  2549. *************************************************************
  2550. *************************************************************
  2551. *------- 'INOU' -------*
  2552. *************************************************************
  2553. *************************************************************
  2554. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  2555. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  2556. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 15. 'UX' -15. ;
  2557. VN1 = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' -15. 'UX' 25. ;
  2558. VN2 = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' -15. 'UX' 25. ;
  2559. YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
  2560. 'H2O' 0.1 ;
  2561.  
  2562. GN = RN '*' VN ;
  2563. GN2 = RN '*' VN1 ;
  2564. *---------------------------------------------------------
  2565. YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
  2566. YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
  2567. YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
  2568. YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
  2569. YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
  2570. 'H2O' YH2O 'N2' YN2 ;
  2571. LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
  2572. CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
  2573. CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
  2574. GAMN = CPTOT '/' CVTOT ;
  2575. GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
  2576. *----------------------------------------------------------
  2577. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  2578. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
  2579. RETN2 =(('INVERSE' (GAMN '-' 1)) '*' PN) '+'
  2580. (0.5 '*' ('PSCAL' VN1 GN2
  2581. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
  2582.  
  2583. RHTN = RETN '+' PN ;
  2584. HTN = RHTN '/' RN ;
  2585. UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
  2586. UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
  2587. UTC = -1 '*' (UXC '+' UYC) '/' (2. '**' 0.5) ;
  2588.  
  2589. SN = PN '/' (RN '**' ('MAXI' GAMN)) ;
  2590.  
  2591. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 6 'HT' ('MAXIMUM' HTN)
  2592. 'S' ('MAXIMUM' SN) 'H2' 0.3 'O2' 0.5 'H2O' 0.1 'PN' 50.0 ;
  2593.  
  2594. CHPLI1 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'HT' ('MAXIMUM' HTN)
  2595. 'S' ('MAXIMUM' SN) 'H2' 0.3 'O2' 0.5 'H2O' 0.1 ;
  2596.  
  2597. CHPLI2 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 1 'PN' 50.0 ;
  2598.  
  2599. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
  2600. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
  2601.  
  2602. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2603. PGAS LISTC LISTP
  2604. RN VN PN YN CHPLI 'INOU' ;
  2605.  
  2606. RCHLIM1 RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2607. PGAS LISTC LISTP
  2608. RN VN PN YN CHPLI1 'INSU' ;
  2609.  
  2610. ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
  2611. ERRO1 = 'MAXIMUM' (RCHLIM1 '-' RCHLIM) 'ABS' ;
  2612.  
  2613. * 'OPTION' DONN 5 ;
  2614.  
  2615. XXTOL = 1.D-9;
  2616. ERR1 = ERRO;
  2617. 'SI' (ERR1 > XXTOL) ;
  2618. CHAI1 = CHAI 'ERREUR NUMERO 67: ' ERR1 'POUR ' XXTOL;
  2619. 'MESS' CHAI1;
  2620. 'ERREUR' 5 ;
  2621. 'FINSI' ;
  2622.  
  2623. XXTOL = 1.D-12;
  2624. ERR1 = ERRO1;
  2625. 'SI' (ERR1 > XXTOL) ;
  2626. CHAI1 = CHAI 'ERREUR NUMERO 68: ' ERR1 'POUR ' XXTOL;
  2627. 'MESS' CHAI1;
  2628. 'ERREUR' 5 ;
  2629. 'FINSI' ;
  2630. *--------------------------------------------------------------
  2631. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2632. PGAS LISTC LISTP
  2633. RN VN1 PN YN CHPLI 'INOU' ;
  2634.  
  2635.  
  2636. RCHLIM2 RCHRES2 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2637. PGAS LISTC LISTP
  2638. RN VN1 PN YN CHPLI2 'OUTP' ;
  2639.  
  2640. ERRO = 'MAXIMUM' (RCHRES2 '-' RCHRES) 'ABS' ;
  2641. ERRO1 = 'MAXIMUM' (RCHLIM2 '-' RCHLIM) 'ABS' ;
  2642.  
  2643. XXTOL = 1.D-12;
  2644. ERR1 = ERRO;
  2645. 'SI' (ERR1 > XXTOL) ;
  2646. CHAI1 = CHAI 'ERREUR NUMERO 69: ' ERR1 'POUR ' XXTOL;
  2647. 'MESS' CHAI1;
  2648. 'ERREUR' 5 ;
  2649. 'FINSI' ;
  2650.  
  2651. XXTOL = 1.D-12;
  2652. ERR1 = ERRO1;
  2653. 'SI' (ERR1 > XXTOL) ;
  2654. CHAI1 = CHAI 'ERREUR NUMERO 70: ' ERR1 'POUR ' XXTOL;
  2655. 'MESS' CHAI1;
  2656. 'ERREUR' 5 ;
  2657. 'FINSI' ;
  2658.  
  2659. *------------------------------------------------------------
  2660. ***** Jacobian with respect to the primitive variables ******
  2661. *------------------------------------------------------------
  2662. CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 6 'HT' 250.0
  2663. 'S' 2.0 'H2' 0.1 'O2' 0.7 'H2O' 0.2 'PN' 50.0 ;
  2664.  
  2665. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2666. PGAS LISTC LISTP
  2667. RN VN PN YN CHPLI 'INOU' ;
  2668.  
  2669. RCHLIM2 RCHRES2 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2670. PGAS LISTC LISTP
  2671. RN VN1 PN YN CHPLI 'INOU' ;
  2672.  
  2673. RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
  2674. PGAS LISTC LISTP
  2675. RN VN PN YN CHPLI 'INOU' ;
  2676.  
  2677.  
  2678. RJACO2 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
  2679. PGAS LISTC LISTP
  2680. RN VN1 PN YN CHPLI 'INOU' ;
  2681. ************* RN *********** RN ***********************
  2682.  
  2683. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  2684. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2685. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2686.  
  2687. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2688. PGAS LISTC LISTP
  2689. (RN '+' DR) VN PN YN CHPLI 'INOU' ;
  2690.  
  2691. RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2692. PGAS LISTC LISTP
  2693. (RN '+' DR) VN1 PN YN CHPLI 'INOU' ;
  2694.  
  2695. DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL ;
  2696. DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
  2697. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2698.  
  2699. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2700. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
  2701. 'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2702.  
  2703. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2704. 'SI' (MMM < ERROM) ;
  2705. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2706. 'SINON' ;
  2707. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2708. 'FINSI' ;
  2709. *---------------------------
  2710. MMM = 'MAXIMUM' DRESDR2 'ABS' ;
  2711. 'SI' (MMM < ERROM) ;
  2712. ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
  2713. 'SINON' ;
  2714. ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
  2715. 'FINSI' ;
  2716.  
  2717. XXTOL = ERROM;
  2718. ERR1 = ERRO;
  2719. 'SI' (ERR1 > XXTOL) ;
  2720. CHAI1 = CHAI 'ERREUR NUMERO 71: ' ERR1 'POUR ' XXTOL;
  2721. 'MESS' CHAI1;
  2722. 'ERREUR' 5 ;
  2723. 'FINSI' ;
  2724.  
  2725. XXTOL = ERROM;
  2726. ERR1 = ERRO2;
  2727. 'SI' (ERR1 > XXTOL) ;
  2728. CHAI1 = CHAI 'ERREUR NUMERO 72: ' ERR1 'POUR ' XXTOL;
  2729. 'MESS' CHAI1;
  2730. 'ERREUR' 5 ;
  2731. 'FINSI' ;
  2732.  
  2733.  
  2734. ************* UX *********** UX ***********************
  2735.  
  2736. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  2737. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2738. DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
  2739. *-------------------------------------------
  2740. DRSCAL2 = EPS '*' ('EXTRAIRE' VN2 'UX' PCEN) ;
  2741. DR2 = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL2) ;
  2742. *-------------------------------------------
  2743.  
  2744. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2745. PGAS LISTC LISTP
  2746. RN (VN '+' DR) PN YN CHPLI 'INOU' ;
  2747.  
  2748. RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2749. PGAS LISTC LISTP
  2750. RN (VN1 '+' DR2) PN YN CHPLI 'INOU' ;
  2751.  
  2752. DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL2 ;
  2753. DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2754. 'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2755.  
  2756. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2757. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2758. 'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2759.  
  2760. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2761. 'SI' (MMM < ERROM) ;
  2762. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2763. 'SINON' ;
  2764. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2765. 'FINSI' ;
  2766. *---------------------------
  2767. MMM = 'MAXIMUM' DRESDR2 'ABS' ;
  2768. 'SI' (MMM < ERROM) ;
  2769. ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
  2770. 'SINON' ;
  2771. ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
  2772. 'FINSI' ;
  2773.  
  2774. XXTOL = ERROM;
  2775. ERR1 = ERRO;
  2776. 'SI' (ERR1 > XXTOL) ;
  2777. CHAI1 = CHAI 'ERREUR NUMERO 73: ' ERR1 'POUR ' XXTOL;
  2778. 'MESS' CHAI1;
  2779. 'ERREUR' 5 ;
  2780. 'FINSI' ;
  2781.  
  2782. XXTOL = ERROM;
  2783. ERR1 = ERRO2;
  2784. 'SI' (ERR1 > XXTOL) ;
  2785. CHAI1 = CHAI 'ERREUR NUMERO 74: ' ERR1 'POUR ' XXTOL;
  2786. 'MESS' CHAI1;
  2787. 'ERREUR' 5 ;
  2788. 'FINSI' ;
  2789.  
  2790. ************* UY *********** UY ***********************
  2791.  
  2792. DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  2793. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2794. DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
  2795. *-------------------------------------------
  2796. DRSCAL2 = EPS '*' ('EXTRAIRE' VN2 'UY' PCEN) ;
  2797. DR2 = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL2) ;
  2798. *-------------------------------------------
  2799.  
  2800. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2801. PGAS LISTC LISTP
  2802. RN (VN '+' DR) PN YN CHPLI 'INOU' ;
  2803.  
  2804. RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2805. PGAS LISTC LISTP
  2806. RN (VN1 '+' DR2) PN YN CHPLI 'INOU' ;
  2807.  
  2808. DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL2 ;
  2809. DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2810. 'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2811.  
  2812. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2813. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2814. 'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2815.  
  2816. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2817. 'SI' (MMM < ERROM) ;
  2818. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2819. 'SINON' ;
  2820. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2821. 'FINSI' ;
  2822. *---------------------------
  2823. MMM = 'MAXIMUM' DRESDR2 'ABS' ;
  2824. 'SI' (MMM < ERROM) ;
  2825. ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
  2826. 'SINON' ;
  2827. ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
  2828. 'FINSI' ;
  2829.  
  2830. XXTOL = ERROM;
  2831. ERR1 = ERRO;
  2832. 'SI' (ERR1 > XXTOL) ;
  2833. CHAI1 = CHAI 'ERREUR NUMERO 75: ' ERR1 'POUR ' XXTOL;
  2834. 'MESS' CHAI1;
  2835. 'ERREUR' 5 ;
  2836. 'FINSI' ;
  2837.  
  2838. XXTOL = ERROM;
  2839. ERR1 = ERRO2;
  2840. 'SI' (ERR1 > XXTOL) ;
  2841. CHAI1 = CHAI 'ERREUR NUMERO 76: ' ERR1 'POUR ' XXTOL;
  2842. 'MESS' CHAI1;
  2843. 'ERREUR' 5 ;
  2844. 'FINSI' ;
  2845.  
  2846.  
  2847. ************* PN *********** PN ***********************
  2848.  
  2849. DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  2850. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2851. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2852.  
  2853. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2854. PGAS LISTC LISTP
  2855. RN VN (PN '+' DR) YN CHPLI 'INOU' ;
  2856.  
  2857. RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
  2858. PGAS LISTC LISTP
  2859. RN VN1 (PN '+' DR) YN CHPLI 'INOU' ;
  2860.  
  2861. DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL ;
  2862. DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2863. 'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2864.  
  2865. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2866. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
  2867. 'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
  2868.  
  2869. MMM = 'MAXIMUM' DRESDR 'ABS' ;
  2870. 'SI' (MMM < ERROM) ;
  2871. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2872. 'SINON' ;
  2873. ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
  2874. 'FINSI' ;
  2875. *---------------------------
  2876. MMM = 'MAXIMUM' DRESDR2 'ABS' ;
  2877. 'SI' (MMM < ERROM) ;
  2878. ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
  2879. 'SINON' ;
  2880. ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
  2881. 'FINSI' ;
  2882.  
  2883. XXTOL = ERROM;
  2884. ERR1 = ERRO;
  2885. 'SI' (ERR1 > XXTOL) ;
  2886. CHAI1 = CHAI 'ERREUR NUMERO 77: ' ERR1 'POUR ' XXTOL;
  2887. 'MESS' CHAI1;
  2888. 'ERREUR' 5 ;
  2889. 'FINSI' ;