Télécharger clim3dj.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : clim3dj.dgibi
  2. *********************************************************************
  3. * VF, CLIM 3D *
  4. * Jacobians in a rotated mesh *
  5. * *
  6. * BECCANTINI A., DM2S/SFME/LTMF, JANVIER 2002 *
  7. *********************************************************************
  8.  
  9. TYEL = 'CU20' ;
  10.  
  11. 'OPTION' 'DIME' 3 'ELEM' TYEL 'ISOV' 'SULI' 'ECHO' 1 'TRAC' 'X' ;
  12.  
  13. GRAPH = VRAI ;
  14. GRAPH = FAUX ;
  15.  
  16. ******************
  17. **** MAILLAGE ****
  18. ******************
  19. *
  20. * P5 P4
  21. *
  22. * P2 P3
  23. *
  24.  
  25. DX = 0.5 ;
  26. P2 = 0.0 0.0 0.0 ;
  27. P3 = DX 0.0 0.0 ;
  28. P4 = DX DX 0.0 ;
  29. P5 = 0.0 DX 0.0 ;
  30. P2P3 = P2 'DROIT' 1 P3 ;
  31. P5P4 = P5 'DROIT' 1 P4 ;
  32. ENTREE = P2P3 'REGLER' 1 P5P4 ;
  33. DOMINT = ENTREE 'VOLU' 'TRAN' 4 (0.0 0.0 (4 '*' DX)) ;
  34.  
  35. ENTREE = ENTREE 'TOURNER' 15 (-1 0.0 0.0) (2.0 0.0 0.0) ;
  36. DOMINT = DOMINT 'TOURNER' 15 (-1 0.0 0.0) (2.0 0.0 0.0) ;
  37.  
  38. 'ELIMINATION' ENTREE DOMINT (DX '/' 100) ;
  39.  
  40. DOMINT = 'CHANGER' DOMINT 'QUAF' ;
  41. ENTREE = 'CHANGER' ENTREE 'QUAF' ;
  42.  
  43. 'ELIMINATION' DOMINT 0.1D-3 ENTREE ;
  44.  
  45. *
  46. **** Creation of DOMAINE tables via the MODEL object
  47. *
  48.  
  49. MDOMINT = 'MODELISER' DOMINT 'NAVIER_STOKES' 'LINE' ;
  50. MENTREE = 'MODELISER' ENTREE 'NAVIER_STOKES' 'LINE' ;
  51.  
  52. $DOMINT = 'DOMA' MDOMINT 'TABLE' ;
  53. $ENTREE = 'DOMA' MENTREE 'TABLE' ;
  54.  
  55. DOMINT = 'DOMA' MDOMINT 'MAILLAGE' ;
  56. ENTREE = 'DOMA' MENTREE 'MAILLAGE' ;
  57.  
  58. PCEN = ($DOMINT . 'CENTRE') 'POIN' 'PROC'
  59. ((DX '/' 2) (DX '/' 2) (DX '/' 2)) ;
  60.  
  61. 'SI' GRAPH ;
  62. 'TRACER' (DOMINT 'ET' PCEN) ;
  63. 'FINSI' ;
  64.  
  65.  
  66. **********************************************************************************
  67. **********************************************************************************
  68. *************************** Riemann invariants, inlet ***************************
  69. **********************************************************************************
  70. **********************************************************************************
  71.  
  72. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.15 ;
  73. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.011 ;
  74. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  75. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.19 'UX' 1.5
  76. 'UZ' 1.11 ;
  77.  
  78. GN = RN '*' VN ;
  79. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  80. ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  81.  
  82. CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 5 'RN' 0.91 'PN' 3.0
  83. 'UX' 0.35 'UY' 0.3 'UZ' 1.0 ;
  84.  
  85. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  86. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  87.  
  88. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  89. LISTC LISTP
  90. RN VN PN GAMN CHPLI 'INRI' ;
  91.  
  92. *
  93. **** Jacobian with respect to primitive variable
  94. *
  95. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  96. LISTC LISTP
  97. RN VN PN GAMN CHPLI 'INRI' ;
  98.  
  99. EPS=1.0D-6 ;
  100.  
  101. * RN
  102.  
  103. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  104. ELE1 = 'MANUEL' 'POI1' PCEN ;
  105. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  106.  
  107. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  108. LISTC LISTP
  109. (RN '+' DR) VN PN GAMN CHPLI 'INRI' ;
  110.  
  111. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  112. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  113. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  114.  
  115. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  116.  
  117. 'SI' (ERRO > 1.0D-3) ;
  118. 'ERREUR' 5 ;
  119. 'FINSI' ;
  120.  
  121.  
  122. * UX
  123.  
  124. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  125. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  126.  
  127. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  128. LISTC LISTP
  129. RN (VN '+' DUX) PN GAMN CHPLI 'INRI' ;
  130.  
  131. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  132. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  133. 'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  134.  
  135. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  136.  
  137. 'SI' (ERRO > 1.0D-3) ;
  138. 'ERREUR' 5 ;
  139. 'FINSI' ;
  140.  
  141. * UY
  142.  
  143. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  144. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  145.  
  146. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  147. LISTC LISTP
  148. RN (VN '+' DUY) PN GAMN CHPLI 'INRI' ;
  149.  
  150. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  151. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  152. 'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  153.  
  154. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  155.  
  156. 'SI' (ERRO > 1.0D-3) ;
  157. 'ERREUR' 5 ;
  158. 'FINSI' ;
  159.  
  160. * UZ
  161.  
  162. DUZSCAL = EPS '*' ('EXTRAIRE' VN 'UZ' PCEN) ;
  163. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  164.  
  165. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  166. LISTC LISTP
  167. RN (VN '+' DUZ) PN GAMN CHPLI 'INRI' ;
  168.  
  169. DRESDUZN = ( RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  170. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  171. 'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  172.  
  173. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  174.  
  175. 'SI' (ERRO > 1.0D-3) ;
  176. 'ERREUR' 5 ;
  177. 'FINSI' ;
  178.  
  179. * PN
  180.  
  181. DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  182. ELE1 = 'MANUEL' 'POI1' PCEN ;
  183. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  184.  
  185. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  186. LISTC LISTP
  187. RN VN (PN '+' DP) GAMN CHPLI 'INRI' ;
  188.  
  189. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  190. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  191. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  192.  
  193. ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  194.  
  195. 'SI' (ERRO > 1.0D-3) ;
  196. 'ERREUR' 5 ;
  197. 'FINSI' ;
  198.  
  199. *
  200. **** Jacobian with respect to conservative variables
  201. *
  202.  
  203. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  204. LISTC LISTP
  205. RN VN PN GAMN CHPLI 'INRI' ;
  206.  
  207. EPS=1.0D-6 ;
  208.  
  209. * RN
  210.  
  211. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  212. ELE1 = 'MANUEL' 'POI1' PCEN ;
  213. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  214.  
  215. RN1 = RN '+' DR ;
  216. GN1 = GN ;
  217. RETN1 = RETN ;
  218.  
  219. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  220.  
  221. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  222. LISTC LISTP
  223. RN1 VN1 PN1 GAMN CHPLI 'INRI' ;
  224.  
  225. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  226. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  227. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  228.  
  229. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  230.  
  231. 'SI' (ERRO > 1.0D-3) ;
  232. 'ERREUR' 5 ;
  233. 'FINSI' ;
  234.  
  235. * RUX
  236.  
  237. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  238. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  239.  
  240. RN1 = RN ;
  241. GN1 = GN '+' DUX ;
  242. RETN1 = RETN ;
  243.  
  244. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  245.  
  246. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  247. LISTC LISTP
  248. RN1 VN1 PN1 GAMN CHPLI 'INRI' ;
  249.  
  250. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  251. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  252. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  253.  
  254. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  255.  
  256. 'SI' (ERRO > 1.0D-3) ;
  257. 'ERREUR' 5 ;
  258. 'FINSI' ;
  259.  
  260. * RUY
  261.  
  262. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  263. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  264.  
  265. RN1 = RN ;
  266. GN1 = GN '+' DUY ;
  267. RETN1 = RETN ;
  268.  
  269. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  270.  
  271. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  272. LISTC LISTP
  273. RN1 VN1 PN1 GAMN CHPLI 'INRI' ;
  274.  
  275. DRESDUYN = (RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  276. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  277. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  278.  
  279. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  280.  
  281. 'SI' (ERRO > 1.0D-3) ;
  282. 'ERREUR' 5 ;
  283. 'FINSI' ;
  284.  
  285.  
  286. * RUZ
  287.  
  288. DUZSCAL = EPS '*' ('EXTRAIRE' GN 'UZ' PCEN) ;
  289. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  290.  
  291. RN1 = RN ;
  292. GN1 = GN '+' DUZ ;
  293. RETN1 = RETN ;
  294.  
  295. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  296.  
  297. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  298. LISTC LISTP
  299. RN1 VN1 PN1 GAMN CHPLI 'INRI' ;
  300.  
  301. DRESDUZN = (RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  302. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  303. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  304.  
  305. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  306.  
  307. 'SI' (ERRO > 1.0D-3) ;
  308. 'ERREUR' 5 ;
  309. 'FINSI' ;
  310.  
  311. * RETN
  312.  
  313. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  314. ELE1 = 'MANUEL' 'POI1' PCEN ;
  315. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  316.  
  317. RN1 = RN ;
  318. GN1 = GN ;
  319. RETN1 = RETN '+' DR ;
  320.  
  321. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  322.  
  323. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  324. LISTC LISTP
  325. RN1 VN1 PN1 GAMN CHPLI 'INRI' ;
  326.  
  327. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  328. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  329. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  330.  
  331. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  332.  
  333. 'SI' (ERRO > 1.0D-3) ;
  334. 'ERREUR' 5 ;
  335. 'FINSI' ;
  336.  
  337. **********************************************************************************
  338. **********************************************************************************
  339. *************************** Riemann invariants, outlet **************************
  340. **********************************************************************************
  341. **********************************************************************************
  342.  
  343. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.15 ;
  344. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.011 ;
  345. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  346. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 0.19 'UX' 0.15
  347. 'UZ' -0.11 ;
  348.  
  349. GN = RN '*' VN ;
  350. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  351. ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  352.  
  353. CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 5 'RN' 0.91 'PN' 1.10
  354. 'UX' 0.35 'UY' 0.3 'UZ' 0.011 ;
  355.  
  356. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  357. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  358.  
  359. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  360. LISTC LISTP
  361. RN VN PN GAMN CHPLI 'OUTRI' ;
  362.  
  363. *
  364. **** Jacobian with respect to primitive variable
  365. *
  366.  
  367. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  368. LISTC LISTP
  369. RN VN PN GAMN CHPLI 'OUTRI' ;
  370.  
  371. EPS=1.0D-6 ;
  372.  
  373. * RN
  374.  
  375. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  376. ELE1 = 'MANUEL' 'POI1' PCEN ;
  377. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  378.  
  379. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  380. LISTC LISTP
  381. (RN '+' DR) VN PN GAMN CHPLI 'OUTRI' ;
  382.  
  383. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  384. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  385. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  386.  
  387. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  388.  
  389. 'SI' (ERRO > 1.0D-3) ;
  390. 'ERREUR' 5 ;
  391. 'FINSI' ;
  392.  
  393. * UX
  394.  
  395. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  396. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  397.  
  398. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  399. LISTC LISTP
  400. RN (VN '+' DUX) PN GAMN CHPLI 'OUTRI' ;
  401.  
  402. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  403. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  404. 'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  405.  
  406. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  407.  
  408. 'SI' (ERRO > 1.0D-3) ;
  409. 'ERREUR' 5 ;
  410. 'FINSI' ;
  411.  
  412. * UY
  413.  
  414. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  415. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  416.  
  417. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  418. LISTC LISTP
  419. RN (VN '+' DUY) PN GAMN CHPLI 'OUTRI' ;
  420.  
  421. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  422. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  423. 'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  424.  
  425. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  426.  
  427. 'SI' (ERRO > 1.0D-3) ;
  428. 'ERREUR' 5 ;
  429. 'FINSI' ;
  430.  
  431. * UZ
  432.  
  433. DUZSCAL = EPS '*' ('EXTRAIRE' VN 'UZ' PCEN) ;
  434. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  435.  
  436. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  437. LISTC LISTP
  438. RN (VN '+' DUZ) PN GAMN CHPLI 'OUTRI' ;
  439.  
  440. DRESDUZN = ( RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  441. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  442. 'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  443.  
  444. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  445.  
  446. 'SI' (ERRO > 1.0D-3) ;
  447. 'ERREUR' 5 ;
  448. 'FINSI' ;
  449.  
  450. * PN
  451.  
  452. DPSCAL = EPS '' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  453. ELE1 = 'MANUEL' 'POI1' PCEN ;
  454. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  455.  
  456. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  457. LISTC LISTP
  458. RN VN (PN '+' DP) GAMN CHPLI 'OUTRI' ;
  459.  
  460. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  461. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  462. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  463.  
  464. ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  465.  
  466. 'SI' (ERRO > 1.0D-3) ;
  467. 'ERREUR' 5 ;
  468. 'FINSI' ;
  469.  
  470. *
  471. **** Jacobian with respect to conservative variables
  472. *
  473.  
  474. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  475. LISTC LISTP
  476. RN VN PN GAMN CHPLI 'OUTRI' ;
  477.  
  478. EPS=1.0D-6 ;
  479.  
  480. * RN
  481.  
  482. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  483. ELE1 = 'MANUEL' 'POI1' PCEN ;
  484. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  485.  
  486. RN1 = RN '+' DR ;
  487. GN1 = GN ;
  488. RETN1 = RETN ;
  489.  
  490. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  491.  
  492. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  493. LISTC LISTP
  494. RN1 VN1 PN1 GAMN CHPLI 'OUTRI' ;
  495.  
  496. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  497. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  498. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  499.  
  500. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  501.  
  502. 'SI' (ERRO > 1.0D-3) ;
  503. 'ERREUR' 5 ;
  504. 'FINSI' ;
  505.  
  506. * RUX
  507.  
  508. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  509. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  510.  
  511. RN1 = RN ;
  512. GN1 = GN '+' DUX ;
  513. RETN1 = RETN ;
  514.  
  515. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  516.  
  517. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  518. LISTC LISTP
  519. RN1 VN1 PN1 GAMN CHPLI 'OUTRI' ;
  520.  
  521. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  522. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  523. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  524.  
  525. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  526.  
  527. 'SI' (ERRO > 1.0D-3) ;
  528. 'ERREUR' 5 ;
  529. 'FINSI' ;
  530.  
  531. * RUY
  532.  
  533. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  534. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  535.  
  536. RN1 = RN ;
  537. GN1 = GN '+' DUY ;
  538. RETN1 = RETN ;
  539.  
  540. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  541.  
  542. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  543. LISTC LISTP
  544. RN1 VN1 PN1 GAMN CHPLI 'OUTRI' ;
  545.  
  546. DRESDUYN = (RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  547. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  548. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  549.  
  550. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  551.  
  552. 'SI' (ERRO > 1.0D-3) ;
  553. 'ERREUR' 5 ;
  554. 'FINSI' ;
  555.  
  556.  
  557. * RUZ
  558.  
  559. DUZSCAL = EPS '*' ('EXTRAIRE' GN 'UZ' PCEN) ;
  560. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  561.  
  562. RN1 = RN ;
  563. GN1 = GN '+' DUZ ;
  564. RETN1 = RETN ;
  565.  
  566. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  567.  
  568. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  569. LISTC LISTP
  570. RN1 VN1 PN1 GAMN CHPLI 'OUTRI' ;
  571.  
  572. DRESDUZN = (RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  573. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  574. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  575.  
  576. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  577.  
  578. 'SI' (ERRO > 1.0D-3) ;
  579. 'ERREUR' 5 ;
  580. 'FINSI' ;
  581.  
  582. * RETN
  583.  
  584. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  585. ELE1 = 'MANUEL' 'POI1' PCEN ;
  586. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  587.  
  588. RN1 = RN ;
  589. GN1 = GN ;
  590. RETN1 = RETN '+' DR ;
  591.  
  592. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  593.  
  594. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  595. LISTC LISTP
  596. RN1 VN1 PN1 GAMN CHPLI 'OUTRI' ;
  597.  
  598. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  599. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  600. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  601.  
  602. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  603.  
  604. 'SI' (ERRO > 1.0D-3) ;
  605. 'ERREUR' 5 ;
  606. 'FINSI' ;
  607.  
  608. **********************************************************************************
  609. **********************************************************************************
  610. *************************** Supersonic inlet ************************************
  611. **********************************************************************************
  612. **********************************************************************************
  613.  
  614. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  615. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  616. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  617. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 15. 'UX' 1. 'UZ' 1.7;
  618.  
  619. GN = RN '*' VN ;
  620. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  621. ('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
  622.  
  623. CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 5 'RN' 1.11 'PN' 2.13
  624. 'UX' -11. 'UY' 2. 'UZ' 1.8 ;
  625.  
  626. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  627. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  628.  
  629. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  630. LISTC LISTP
  631. RN VN PN GAMN CHPLI 'INSS' ;
  632.  
  633. ERRO = 'MAXIMUM' (RCHLIM '-' CHPLI) 'ABS' ;
  634.  
  635. 'SI' (ERRO > 1.0D-12) ;
  636. 'ERREUR' 5 ;
  637. 'FINSI' ;
  638.  
  639. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  640. LISTC LISTP
  641. RN VN PN GAMN CHPLI 'INSS' ;
  642.  
  643. AA = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 5 'VAR1' 1.11 'VAR2' 1.13
  644. 'VAR3' 1.15 'VAR4' 1.17 'VAR5' 1.37 ;
  645.  
  646. ERRO = 'MAXIMUM' ('KOPS' 'MULT' RJACO AA) 'ABS' ;
  647.  
  648. 'SI' (ERRO > 1.0D-12) ;
  649. 'ERREUR' 5 ;
  650. 'FINSI' ;
  651.  
  652. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  653. LISTC LISTP
  654. RN VN PN GAMN CHPLI 'INSS' ;
  655.  
  656. AA = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'RN' 1.11 'UX' 1.13
  657. 'UY' 1.15 'PN' 1.17 ;
  658.  
  659. ERRO = 'MAXIMUM' ('KOPS' 'MULT' RJACO AA) 'ABS' ;
  660.  
  661. 'SI' (ERRO > 1.0D-12) ;
  662. 'ERREUR' 5 ;
  663. 'FINSI' ;
  664.  
  665.  
  666. **********************************************************************************
  667. **********************************************************************************
  668. *************************** Supersonic outlet ***********************************
  669. **********************************************************************************
  670. **********************************************************************************
  671.  
  672. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  673. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  674. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  675. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' -15. 'UX' 1. 'UZ' -5.
  676. ;
  677.  
  678. GN = RN '*' VN ;
  679. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  680. ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  681.  
  682. CHPLI CACCA = 'KOPS' MATRIK ;
  683.  
  684. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  685. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  686.  
  687. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  688. LISTC LISTP
  689. RN VN PN GAMN CHPLI 'OUTSS' ;
  690.  
  691. *
  692. **** Jacobian with respect to primitive variable
  693. *
  694.  
  695. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  696. LISTC LISTP
  697. RN VN PN GAMN CHPLI 'OUTSS' ;
  698.  
  699. EPS=1.0D-6 ;
  700.  
  701. * RN
  702.  
  703. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  704. ELE1 = 'MANUEL' 'POI1' PCEN ;
  705. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  706.  
  707. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  708. LISTC LISTP
  709. (RN '+' DR) VN PN GAMN CHPLI 'OUTSS' ;
  710.  
  711. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  712. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  713. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  714.  
  715. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  716.  
  717.  
  718. 'SI' (ERRO > 1.0D-3) ;
  719. 'ERREUR' 5 ;
  720. 'FINSI' ;
  721.  
  722. * UX
  723.  
  724. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  725. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  726.  
  727. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  728. LISTC LISTP
  729. RN (VN '+' DUX) PN GAMN CHPLI 'OUTSS' ;
  730.  
  731. DRESDUXN = (RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  732. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  733. 'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  734.  
  735. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  736.  
  737. 'SI' (ERRO > 1.0D-3) ;
  738. 'ERREUR' 5 ;
  739. 'FINSI' ;
  740.  
  741. * UY
  742.  
  743. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  744. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  745.  
  746. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  747. LISTC LISTP
  748. RN (VN '+' DUY) PN GAMN CHPLI 'OUTSS' ;
  749.  
  750. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  751. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  752. 'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  753.  
  754. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  755.  
  756. 'SI' (ERRO > 1.0D-3) ;
  757. 'ERREUR' 5 ;
  758. 'FINSI' ;
  759.  
  760. * UZ
  761.  
  762. DUZSCAL = EPS '*' ('EXTRAIRE' VN 'UZ' PCEN) ;
  763. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  764.  
  765. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  766. LISTC LISTP
  767. RN (VN '+' DUZ) PN GAMN CHPLI 'OUTSS' ;
  768.  
  769. DRESDUZN = ( RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  770. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  771. 'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  772.  
  773. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  774.  
  775. 'SI' (ERRO > 1.0D-3) ;
  776. 'ERREUR' 5 ;
  777. 'FINSI' ;
  778.  
  779. * PN
  780.  
  781. DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  782. ELE1 = 'MANUEL' 'POI1' PCEN ;
  783. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  784.  
  785. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  786. LISTC LISTP
  787. RN VN (PN '+' DP) GAMN CHPLI 'OUTSS' ;
  788.  
  789. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  790. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  791. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  792.  
  793. ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  794.  
  795. 'SI' (ERRO > 1.0D-3) ;
  796. 'ERREUR' 5 ;
  797. 'FINSI' ;
  798.  
  799. *
  800. **** Jacobian with respect to conservative variables
  801. *
  802.  
  803. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  804. LISTC LISTP
  805. RN VN PN GAMN CHPLI 'OUTSS' ;
  806.  
  807. EPS=1.0D-6 ;
  808.  
  809. * RN
  810.  
  811. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  812. ELE1 = 'MANUEL' 'POI1' PCEN ;
  813. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  814.  
  815. RN1 = RN '+' DR ;
  816. GN1 = GN ;
  817. RETN1 = RETN ;
  818.  
  819. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  820.  
  821. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  822. LISTC LISTP
  823. RN1 VN1 PN1 GAMN CHPLI 'OUTSS' ;
  824.  
  825. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  826. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  827. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  828.  
  829. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  830.  
  831. 'SI' (ERRO > 1.0D-3) ;
  832. 'ERREUR' 5 ;
  833. 'FINSI' ;
  834.  
  835. * RUX
  836.  
  837. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  838. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  839.  
  840. RN1 = RN ;
  841. GN1 = GN '+' DUX ;
  842. RETN1 = RETN ;
  843.  
  844. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  845.  
  846. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  847. LISTC LISTP
  848. RN1 VN1 PN1 GAMN CHPLI 'OUTSS' ;
  849.  
  850. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  851. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  852. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  853.  
  854. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  855.  
  856. 'SI' (ERRO > 1.0D-3) ;
  857. 'ERREUR' 5 ;
  858. 'FINSI' ;
  859.  
  860. * RUY
  861.  
  862. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  863. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  864.  
  865. RN1 = RN ;
  866. GN1 = GN '+' DUY ;
  867. RETN1 = RETN ;
  868.  
  869. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  870.  
  871. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  872. LISTC LISTP
  873. RN1 VN1 PN1 GAMN CHPLI 'OUTSS' ;
  874.  
  875. DRESDUYN = (RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  876. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  877. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  878.  
  879. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  880.  
  881. 'SI' (ERRO > 1.0D-3) ;
  882. 'ERREUR' 5 ;
  883. 'FINSI' ;
  884.  
  885.  
  886. * RUZ
  887.  
  888. DUZSCAL = EPS '*' ('EXTRAIRE' GN 'UZ' PCEN) ;
  889. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  890.  
  891. RN1 = RN ;
  892. GN1 = GN '+' DUZ ;
  893. RETN1 = RETN ;
  894.  
  895. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  896.  
  897. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  898. LISTC LISTP
  899. RN1 VN1 PN1 GAMN CHPLI 'OUTSS' ;
  900.  
  901. DRESDUZN = (RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  902. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  903. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  904.  
  905. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  906.  
  907. 'SI' (ERRO > 1.0D-3) ;
  908. 'ERREUR' 5 ;
  909. 'FINSI' ;
  910.  
  911. * RETN
  912.  
  913. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  914. ELE1 = 'MANUEL' 'POI1' PCEN ;
  915. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  916.  
  917. RN1 = RN ;
  918. GN1 = GN ;
  919. RETN1 = RETN '+' DR ;
  920.  
  921. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  922.  
  923. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  924. LISTC LISTP
  925. RN1 VN1 PN1 GAMN CHPLI 'OUTSS' ;
  926.  
  927. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  928. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  929. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  930.  
  931. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  932.  
  933. 'SI' (ERRO > 1.0D-3) ;
  934. 'ERREUR' 5 ;
  935. 'FINSI' ;
  936.  
  937. **********************************************************************************
  938. **********************************************************************************
  939. *************************** Pressure at the outlet ******************************
  940. **********************************************************************************
  941. **********************************************************************************
  942. *
  943. * We solve a Riemann probleme between the state in the center and
  944. * the same state but the pressure equal to the pressure at the outlet
  945. * We use AUSMPLUS
  946. *
  947. *
  948. ***** Calcul du residu dans le cas particulier PF = PC
  949. *
  950.  
  951. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  952. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  953. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  954. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' -15. 'UX' 1.
  955. 'UZ' -3.1 ;
  956.  
  957. GN = RN '*' VN ;
  958. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  959. ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  960.  
  961. CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 1 'PN' 2.1 ;
  962.  
  963. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  964. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  965.  
  966. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  967. LISTC LISTP
  968. RN VN PN GAMN CHPLI 'OUTP' ;
  969.  
  970. *
  971. **** Jacobian with respect to primitive variable
  972. *
  973.  
  974. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  975. LISTC LISTP
  976. RN VN PN GAMN CHPLI 'OUTP' ;
  977.  
  978. EPS=1.0D-6 ;
  979.  
  980. * RN
  981.  
  982. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  983. ELE1 = 'MANUEL' 'POI1' PCEN ;
  984. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  985.  
  986. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  987. LISTC LISTP
  988. (RN '+' DR) VN PN GAMN CHPLI 'OUTP' ;
  989.  
  990. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  991. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  992. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  993.  
  994. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  995.  
  996. 'SI' (ERRO > 1.0D-3) ;
  997. 'ERREUR' 5 ;
  998. 'FINSI' ;
  999.  
  1000.  
  1001. * UX
  1002.  
  1003. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  1004. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1005.  
  1006. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1007. LISTC LISTP
  1008. RN (VN '+' DUX) PN GAMN CHPLI 'OUTP' ;
  1009.  
  1010. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1011. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1012. 'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1013.  
  1014. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1015.  
  1016. 'SI' (ERRO > 1.0D-3) ;
  1017. 'ERREUR' 5 ;
  1018. 'FINSI' ;
  1019.  
  1020. * UY
  1021.  
  1022. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  1023. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1024.  
  1025. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1026. LISTC LISTP
  1027. RN (VN '+' DUY) PN GAMN CHPLI 'OUTP' ;
  1028.  
  1029. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  1030. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1031. 'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1032.  
  1033. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1034.  
  1035. 'SI' (ERRO > 1.0D-3) ;
  1036. 'ERREUR' 5 ;
  1037. 'FINSI' ;
  1038.  
  1039. * UZ
  1040.  
  1041. DUZSCAL = EPS '*' ('EXTRAIRE' VN 'UZ' PCEN) ;
  1042. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1043.  
  1044. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1045. LISTC LISTP
  1046. RN (VN '+' DUZ) PN GAMN CHPLI 'OUTP' ;
  1047.  
  1048. DRESDUZN = ( RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  1049. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1050. 'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1051.  
  1052. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1053.  
  1054. 'SI' (ERRO > 1.0D-3) ;
  1055. 'ERREUR' 5 ;
  1056. 'FINSI' ;
  1057.  
  1058. * PN
  1059.  
  1060. DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1061. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1062. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1063.  
  1064. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1065. LISTC LISTP
  1066. RN VN (PN '+' DP) GAMN CHPLI 'OUTP' ;
  1067.  
  1068. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  1069. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1070. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1071.  
  1072. ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1073.  
  1074. 'SI' (ERRO > 1.0D-3) ;
  1075. 'ERREUR' 5 ;
  1076. 'FINSI' ;
  1077.  
  1078. *
  1079. **** Jacobian with respect to conservative variables
  1080. *
  1081.  
  1082. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  1083. LISTC LISTP
  1084. RN VN PN GAMN CHPLI 'OUTP' ;
  1085.  
  1086. EPS=1.0D-6 ;
  1087.  
  1088. * RN
  1089.  
  1090. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1091. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1092. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1093.  
  1094. RN1 = RN '+' DR ;
  1095. GN1 = GN ;
  1096. RETN1 = RETN ;
  1097.  
  1098. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1099.  
  1100. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1101. LISTC LISTP
  1102. RN1 VN1 PN1 GAMN CHPLI 'OUTP' ;
  1103.  
  1104. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1105. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1106. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1107.  
  1108. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1109.  
  1110. 'SI' (ERRO > 1.0D-3) ;
  1111. 'ERREUR' 5 ;
  1112. 'FINSI' ;
  1113.  
  1114. * RUX
  1115.  
  1116. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  1117. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1118.  
  1119. RN1 = RN ;
  1120. GN1 = GN '+' DUX ;
  1121. RETN1 = RETN ;
  1122.  
  1123. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1124.  
  1125. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1126. LISTC LISTP
  1127. RN1 VN1 PN1 GAMN CHPLI 'OUTP' ;
  1128.  
  1129. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1130. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1131. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1132.  
  1133. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1134.  
  1135. 'SI' (ERRO > 1.0D-3) ;
  1136. 'ERREUR' 5 ;
  1137. 'FINSI' ;
  1138.  
  1139. * RUY
  1140.  
  1141. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  1142. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1143.  
  1144. RN1 = RN ;
  1145. GN1 = GN '+' DUY ;
  1146. RETN1 = RETN ;
  1147.  
  1148. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1149.  
  1150. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1151. LISTC LISTP
  1152. RN1 VN1 PN1 GAMN CHPLI 'OUTP' ;
  1153.  
  1154. DRESDUYN = (RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  1155. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1156. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1157.  
  1158. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1159.  
  1160. 'SI' (ERRO > 1.0D-3) ;
  1161. 'ERREUR' 5 ;
  1162. 'FINSI' ;
  1163.  
  1164.  
  1165. * RUZ
  1166.  
  1167. DUZSCAL = EPS '*' ('EXTRAIRE' GN 'UZ' PCEN) ;
  1168. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1169.  
  1170. RN1 = RN ;
  1171. GN1 = GN '+' DUZ ;
  1172. RETN1 = RETN ;
  1173.  
  1174. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1175.  
  1176. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1177. LISTC LISTP
  1178. RN1 VN1 PN1 GAMN CHPLI 'OUTP' ;
  1179.  
  1180. DRESDUZN = (RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  1181. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1182. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  1183.  
  1184. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1185.  
  1186. 'SI' (ERRO > 1.0D-3) ;
  1187. 'ERREUR' 5 ;
  1188. 'FINSI' ;
  1189.  
  1190. * RETN
  1191.  
  1192. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1193. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1194. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1195.  
  1196. RN1 = RN ;
  1197. GN1 = GN ;
  1198. RETN1 = RETN '+' DR ;
  1199.  
  1200. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1201.  
  1202. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1203. LISTC LISTP
  1204. RN1 VN1 PN1 GAMN CHPLI 'OUTP' ;
  1205.  
  1206. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1207. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1208. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  1209.  
  1210. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1211.  
  1212. 'SI' (ERRO > 1.0D-3) ;
  1213. 'ERREUR' 5 ;
  1214. 'FINSI' ;
  1215.  
  1216. **********************************************************************************
  1217. **********************************************************************************
  1218. *************************** INSU ************************************************
  1219. **********************************************************************************
  1220. **********************************************************************************
  1221. *
  1222. * Total enthalpy, entropy, tangential velocity at the inlet
  1223. * Normal velocity inside
  1224. *
  1225. * We solve a Riemann probleme between the state at the interface thus
  1226. * reconstructed and the state at the center
  1227. *
  1228.  
  1229. *
  1230. **** Jacobian with respect to primitive variable
  1231. *
  1232.  
  1233. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1234. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1235. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  1236. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.13 'UX' 1.16
  1237. 'UZ' 1.11 ;
  1238. GN = RN * VN ;
  1239. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  1240. ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  1241.  
  1242. RHTN = RETN '+' PN ;
  1243. HTN = RHTN '/' RN ;
  1244.  
  1245. SN = PN '/' (RN '**' ('MAXI' GAMN)) ;
  1246.  
  1247. CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 2 'HT' (('MAXIMUM' HTN) *
  1248. 1.11) 'S' (('MAXIMUM' SN) '*' 1.13) ;
  1249.  
  1250.  
  1251. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  1252. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  1253.  
  1254. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1255. LISTC LISTP
  1256. RN VN PN GAMN CHPLI 'INSU' ;
  1257.  
  1258. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  1259. LISTC LISTP
  1260. RN VN PN GAMN CHPLI 'INSU' ;
  1261.  
  1262. EPS=1.0D-6 ;
  1263.  
  1264. * RN
  1265.  
  1266. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1267. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1268. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1269.  
  1270. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1271. LISTC LISTP
  1272. (RN '+' DR) VN PN GAMN CHPLI 'INSU' ;
  1273.  
  1274. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1275. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  1276. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1277.  
  1278. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1279.  
  1280. 'SI' (ERRO > 1.0D-3) ;
  1281. 'ERREUR' 5 ;
  1282. 'FINSI' ;
  1283.  
  1284.  
  1285. * UX
  1286.  
  1287. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  1288. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1289.  
  1290. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1291. LISTC LISTP
  1292. RN (VN '+' DUX) PN GAMN CHPLI 'INSU' ;
  1293.  
  1294. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1295. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1296. 'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1297.  
  1298. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1299.  
  1300. 'SI' (ERRO > 1.0D-3) ;
  1301. 'ERREUR' 5 ;
  1302. 'FINSI' ;
  1303.  
  1304. * UY
  1305.  
  1306. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  1307. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1308.  
  1309. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1310. LISTC LISTP
  1311. RN (VN '+' DUY) PN GAMN CHPLI 'INSU' ;
  1312.  
  1313. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  1314. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1315. 'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1316.  
  1317. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1318.  
  1319. 'SI' (ERRO > 1.0D-3) ;
  1320. 'ERREUR' 5 ;
  1321. 'FINSI' ;
  1322.  
  1323. * UZ
  1324.  
  1325. DUZSCAL = EPS '*' ('EXTRAIRE' VN 'UZ' PCEN) ;
  1326. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1327.  
  1328. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1329. LISTC LISTP
  1330. RN (VN '+' DUZ) PN GAMN CHPLI 'INSU' ;
  1331.  
  1332. DRESDUZN = ( RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  1333. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1334. 'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1335.  
  1336. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1337.  
  1338. 'SI' (ERRO > 1.0D-3) ;
  1339. 'ERREUR' 5 ;
  1340. 'FINSI' ;
  1341.  
  1342. * PN
  1343.  
  1344. DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1345. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1346. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1347.  
  1348. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1349. LISTC LISTP
  1350. RN VN (PN '+' DP) GAMN CHPLI 'INSU' ;
  1351.  
  1352. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  1353. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1354. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1355.  
  1356. ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1357.  
  1358. 'SI' (ERRO > 1.0D-3) ;
  1359. 'ERREUR' 5 ;
  1360. 'FINSI' ;
  1361.  
  1362. *
  1363. **** Jacobian with respect to conservative variables
  1364. *
  1365.  
  1366. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  1367. LISTC LISTP
  1368. RN VN PN GAMN CHPLI 'INSU' ;
  1369.  
  1370. EPS=1.0D-6 ;
  1371.  
  1372. * RN
  1373.  
  1374. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1375. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1376. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1377.  
  1378. RN1 = RN '+' DR ;
  1379. GN1 = GN ;
  1380. RETN1 = RETN ;
  1381.  
  1382. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1383.  
  1384. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1385. LISTC LISTP
  1386. RN1 VN1 PN1 GAMN CHPLI 'INSU' ;
  1387.  
  1388. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1389. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1390. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1391.  
  1392. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1393.  
  1394. 'SI' (ERRO > 1.0D-3) ;
  1395. 'ERREUR' 5 ;
  1396. 'FINSI' ;
  1397.  
  1398. * RUX
  1399.  
  1400. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  1401. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1402.  
  1403. RN1 = RN ;
  1404. GN1 = GN '+' DUX ;
  1405. RETN1 = RETN ;
  1406.  
  1407. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1408.  
  1409. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1410. LISTC LISTP
  1411. RN1 VN1 PN1 GAMN CHPLI 'INSU' ;
  1412.  
  1413. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1414. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1415. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1416.  
  1417. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1418.  
  1419. 'SI' (ERRO > 1.0D-3) ;
  1420. 'ERREUR' 5 ;
  1421. 'FINSI' ;
  1422.  
  1423. * RUY
  1424.  
  1425. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  1426. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1427.  
  1428. RN1 = RN ;
  1429. GN1 = GN '+' DUY ;
  1430. RETN1 = RETN ;
  1431.  
  1432. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1433.  
  1434. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1435. LISTC LISTP
  1436. RN1 VN1 PN1 GAMN CHPLI 'INSU' ;
  1437.  
  1438. DRESDUYN = (RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  1439. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1440. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1441.  
  1442. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1443.  
  1444. 'SI' (ERRO > 1.0D-3) ;
  1445. 'ERREUR' 5 ;
  1446. 'FINSI' ;
  1447.  
  1448.  
  1449. * RUZ
  1450.  
  1451. DUZSCAL = EPS '*' ('EXTRAIRE' GN 'UZ' PCEN) ;
  1452. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1453.  
  1454. RN1 = RN ;
  1455. GN1 = GN '+' DUZ ;
  1456. RETN1 = RETN ;
  1457.  
  1458. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1459.  
  1460. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1461. LISTC LISTP
  1462. RN1 VN1 PN1 GAMN CHPLI 'INSU' ;
  1463.  
  1464. DRESDUZN = (RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  1465. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1466. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  1467.  
  1468. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1469.  
  1470. 'SI' (ERRO > 1.0D-3) ;
  1471. 'ERREUR' 5 ;
  1472. 'FINSI' ;
  1473.  
  1474. * RETN
  1475.  
  1476. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1477. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1478. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1479.  
  1480. RN1 = RN ;
  1481. GN1 = GN ;
  1482. RETN1 = RETN '+' DR ;
  1483.  
  1484. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1485.  
  1486. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1487. LISTC LISTP
  1488. RN1 VN1 PN1 GAMN CHPLI 'INSU' ;
  1489.  
  1490. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1491. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1492. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  1493.  
  1494. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1495.  
  1496. 'SI' (ERRO > 1.0D-3) ;
  1497. 'ERREUR' 5 ;
  1498. 'FINSI' ;
  1499.  
  1500.  
  1501. **********************************************************************************
  1502. **********************************************************************************
  1503. *************************** INJE ************************************************
  1504. **********************************************************************************
  1505. **********************************************************************************
  1506. *
  1507. * We impose RT and the normal momontum at the inlet.
  1508.  
  1509. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1510. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1511. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  1512. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.11 'UX' -1.3 'UZ'
  1513. 1.5 ;
  1514.  
  1515. GN = RN '*' VN ;
  1516. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  1517. ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  1518.  
  1519. MOME = 0.111 ;
  1520. RT = 4.1 ;
  1521.  
  1522. CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 2 'RT' RT
  1523. 'MOME' MOME ;
  1524.  
  1525. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  1526. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  1527.  
  1528. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1529. LISTC LISTP
  1530. RN VN PN GAMN CHPLI 'INJE' ;
  1531.  
  1532. *
  1533. **** Jacobian with respect to primitive variable
  1534. *
  1535.  
  1536.  
  1537. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1538. LISTC LISTP
  1539. RN VN PN GAMN CHPLI 'INJE' ;
  1540.  
  1541. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  1542. LISTC LISTP
  1543. RN VN PN GAMN CHPLI 'INJE' ;
  1544.  
  1545. EPS=1.0D-6 ;
  1546.  
  1547. * RN
  1548.  
  1549. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1550. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1551. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1552.  
  1553. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1554. LISTC LISTP
  1555. (RN '+' DR) VN PN GAMN CHPLI 'INJE' ;
  1556.  
  1557. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1558. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  1559. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1560.  
  1561. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1562.  
  1563. 'SI' (ERRO > 1.0D-3) ;
  1564. 'ERREUR' 5 ;
  1565. 'FINSI' ;
  1566.  
  1567.  
  1568. * UX
  1569.  
  1570. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  1571. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1572.  
  1573. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1574. LISTC LISTP
  1575. RN (VN '+' DUX) PN GAMN CHPLI 'INJE' ;
  1576.  
  1577. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1578. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1579. 'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1580.  
  1581. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1582.  
  1583. 'SI' (ERRO > 1.0D-3) ;
  1584. 'ERREUR' 5 ;
  1585. 'FINSI' ;
  1586.  
  1587. * UY
  1588.  
  1589. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  1590. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1591.  
  1592. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1593. LISTC LISTP
  1594. RN (VN '+' DUY) PN GAMN CHPLI 'INJE' ;
  1595.  
  1596. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  1597. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1598. 'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1599.  
  1600. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1601.  
  1602. 'SI' (ERRO > 1.0D-3) ;
  1603. 'ERREUR' 5 ;
  1604. 'FINSI' ;
  1605.  
  1606. * UZ
  1607.  
  1608. DUZSCAL = EPS '*' ('EXTRAIRE' VN 'UZ' PCEN) ;
  1609. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1610.  
  1611. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1612. LISTC LISTP
  1613. RN (VN '+' DUZ) PN GAMN CHPLI 'INJE' ;
  1614.  
  1615. DRESDUZN = ( RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  1616. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1617. 'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1618.  
  1619. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1620.  
  1621. 'SI' (ERRO > 1.0D-3) ;
  1622. 'ERREUR' 5 ;
  1623. 'FINSI' ;
  1624.  
  1625. * PN
  1626.  
  1627. DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1628. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1629. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1630.  
  1631. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1632. LISTC LISTP
  1633. RN VN (PN '+' DP) GAMN CHPLI 'INJE' ;
  1634.  
  1635. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  1636. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1637. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1638.  
  1639. ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1640.  
  1641. 'SI' (ERRO > 1.0D-3) ;
  1642. 'ERREUR' 5 ;
  1643. 'FINSI' ;
  1644.  
  1645. *
  1646. **** Jacobian with respect to conservative variables
  1647. *
  1648.  
  1649. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  1650. LISTC LISTP
  1651. RN VN PN GAMN CHPLI 'INJE' ;
  1652.  
  1653. EPS=1.0D-6 ;
  1654.  
  1655. * RN
  1656.  
  1657. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1658. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1659. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1660.  
  1661. RN1 = RN '+' DR ;
  1662. GN1 = GN ;
  1663. RETN1 = RETN ;
  1664.  
  1665. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1666.  
  1667. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1668. LISTC LISTP
  1669. RN1 VN1 PN1 GAMN CHPLI 'INJE' ;
  1670.  
  1671. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1672. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1673. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1674.  
  1675. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1676.  
  1677. 'SI' (ERRO > 1.0D-3) ;
  1678. 'ERREUR' 5 ;
  1679. 'FINSI' ;
  1680.  
  1681. * RUX
  1682.  
  1683. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  1684. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1685.  
  1686. RN1 = RN ;
  1687. GN1 = GN '+' DUX ;
  1688. RETN1 = RETN ;
  1689.  
  1690. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1691.  
  1692. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1693. LISTC LISTP
  1694. RN1 VN1 PN1 GAMN CHPLI 'INJE' ;
  1695.  
  1696. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1697. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1698. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1699.  
  1700. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1701.  
  1702. 'SI' (ERRO > 1.0D-3) ;
  1703. 'ERREUR' 5 ;
  1704. 'FINSI' ;
  1705.  
  1706. * RUY
  1707.  
  1708. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  1709. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1710.  
  1711. RN1 = RN ;
  1712. GN1 = GN '+' DUY ;
  1713. RETN1 = RETN ;
  1714.  
  1715. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1716.  
  1717. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1718. LISTC LISTP
  1719. RN1 VN1 PN1 GAMN CHPLI 'INJE' ;
  1720.  
  1721. DRESDUYN = (RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  1722. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1723. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1724.  
  1725. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1726.  
  1727. 'SI' (ERRO > 1.0D-3) ;
  1728. 'ERREUR' 5 ;
  1729. 'FINSI' ;
  1730.  
  1731.  
  1732. * RUZ
  1733.  
  1734. DUZSCAL = EPS '*' ('EXTRAIRE' GN 'UZ' PCEN) ;
  1735. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1736.  
  1737. RN1 = RN ;
  1738. GN1 = GN '+' DUZ ;
  1739. RETN1 = RETN ;
  1740.  
  1741. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1742.  
  1743. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1744. LISTC LISTP
  1745. RN1 VN1 PN1 GAMN CHPLI 'INJE' ;
  1746.  
  1747. DRESDUZN = (RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  1748. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1749. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  1750.  
  1751. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1752.  
  1753. 'SI' (ERRO > 1.0D-3) ;
  1754. 'ERREUR' 5 ;
  1755. 'FINSI' ;
  1756.  
  1757. * RETN
  1758.  
  1759. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  1760. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1761. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1762.  
  1763. RN1 = RN ;
  1764. GN1 = GN ;
  1765. RETN1 = RETN '+' DR ;
  1766.  
  1767. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1768.  
  1769. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1770. LISTC LISTP
  1771. RN1 VN1 PN1 GAMN CHPLI 'INJE' ;
  1772.  
  1773. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1774. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1775. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  1776.  
  1777. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1778.  
  1779. 'SI' (ERRO > 1.0D-3) ;
  1780. 'ERREUR' 5 ;
  1781. 'FINSI' ;
  1782.  
  1783.  
  1784. **********************************************************************************
  1785. **********************************************************************************
  1786. *************************** INJELM **********************************************
  1787. **********************************************************************************
  1788. **********************************************************************************
  1789. *
  1790. * We impose RT and the normal momontum at the inlet.
  1791.  
  1792. RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
  1793. PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
  1794. GAMN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.37 ;
  1795. VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'UY' 1.11 'UX' -1.3 'UZ'
  1796. 1.5 ;
  1797.  
  1798. GN = RN '*' VN ;
  1799. RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
  1800. ('MOTS' 'UX' 'UY' 'UZ') ('MOTS' 'UX' 'UY' 'UZ'))) ;
  1801.  
  1802. MOME = 0.111 ;
  1803. RT = 4.1 ;
  1804.  
  1805. CHPLI = 'MANUEL' 'CHPO' ($ENTREE . 'CENTRE') 2 'RT' RT
  1806. 'MOME' MOME ;
  1807.  
  1808. LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
  1809. LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'UZ ' 'PN ' ;
  1810.  
  1811. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1812. LISTC LISTP
  1813. RN VN PN GAMN CHPLI 'INJELM' ;
  1814.  
  1815. *
  1816. **** Jacobian with respect to primitive variable
  1817. *
  1818.  
  1819.  
  1820. RCHLIM RCHRES = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1821. LISTC LISTP
  1822. RN VN PN GAMN CHPLI 'INJELM' ;
  1823.  
  1824. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOPRIM' MDOMINT MENTREE
  1825. LISTC LISTP
  1826. RN VN PN GAMN CHPLI 'INJELM' ;
  1827.  
  1828. EPS=1.0D-6 ;
  1829.  
  1830. * RN
  1831.  
  1832. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1833. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1834. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1835.  
  1836. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1837. LISTC LISTP
  1838. (RN '+' DR) VN PN GAMN CHPLI 'INJELM' ;
  1839.  
  1840. DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1841. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 1.0
  1842. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1843.  
  1844. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1845.  
  1846. 'SI' (ERRO > 1.0D-3) ;
  1847. 'ERREUR' 5 ;
  1848. 'FINSI' ;
  1849.  
  1850.  
  1851. * UX
  1852.  
  1853. DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
  1854. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1855.  
  1856. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1857. LISTC LISTP
  1858. RN (VN '+' DUX) PN GAMN CHPLI 'INJELM' ;
  1859.  
  1860. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1861. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1862. 'UX' 1.0 'UY' 0.0 'UZ' 0.0 'PN' 0.0) ;
  1863.  
  1864. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1865.  
  1866. 'SI' (ERRO > 1.0D-3) ;
  1867. 'ERREUR' 5 ;
  1868. 'FINSI' ;
  1869.  
  1870. * UY
  1871.  
  1872. DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
  1873. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1874.  
  1875. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1876. LISTC LISTP
  1877. RN (VN '+' DUY) PN GAMN CHPLI 'INJELM' ;
  1878.  
  1879. DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  1880. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1881. 'UX' 0.0 'UY' 1.0 'UZ' 0.0 'PN' 0.0) ;
  1882.  
  1883. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  1884.  
  1885. 'SI' (ERRO > 1.0D-3) ;
  1886. 'ERREUR' 5 ;
  1887. 'FINSI' ;
  1888.  
  1889. * UZ
  1890.  
  1891. DUZSCAL = EPS '*' ('EXTRAIRE' VN 'UZ' PCEN) ;
  1892. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  1893.  
  1894. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1895. LISTC LISTP
  1896. RN (VN '+' DUZ) PN GAMN CHPLI 'INJELM' ;
  1897.  
  1898. DRESDUZN = ( RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  1899. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1900. 'UX' 0.0 'UY' 0.0 'UZ' 1.0 'PN' 0.0) ;
  1901.  
  1902. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  1903.  
  1904. 'SI' (ERRO > 1.0D-3) ;
  1905. 'ERREUR' 5 ;
  1906. 'FINSI' ;
  1907.  
  1908. * PN
  1909.  
  1910. DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
  1911. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1912. DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
  1913.  
  1914. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1915. LISTC LISTP
  1916. RN VN (PN '+' DP) GAMN CHPLI 'INJELM' ;
  1917.  
  1918. DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
  1919. DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'RN' 0.0
  1920. 'UX' 0.0 'UY' 0.0 'UZ' 0.0 'PN' 1.0) ;
  1921.  
  1922. ERRO = 'MAXIMUM' (DRESDP '-' DRESDPN) 'ABS' ;
  1923.  
  1924. 'SI' (ERRO > 1.0D-3) ;
  1925. 'ERREUR' 5 ;
  1926. 'FINSI' ;
  1927.  
  1928. *
  1929. **** Jacobian with respect to conservative variables
  1930. *
  1931.  
  1932. RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'JACOCONS' MDOMINT MENTREE
  1933. LISTC LISTP
  1934. RN VN PN GAMN CHPLI 'INJELM' ;
  1935.  
  1936. EPS=1.0D-6 ;
  1937.  
  1938. * RN
  1939.  
  1940. DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
  1941. ELE1 = 'MANUEL' 'POI1' PCEN ;
  1942. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  1943.  
  1944. RN1 = RN '+' DR ;
  1945. GN1 = GN ;
  1946. RETN1 = RETN ;
  1947.  
  1948. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1949.  
  1950. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1951. LISTC LISTP
  1952. RN1 VN1 PN1 GAMN CHPLI 'INJELM' ;
  1953.  
  1954. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  1955. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 1.0
  1956. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1957.  
  1958. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  1959.  
  1960. 'SI' (ERRO > 1.0D-3) ;
  1961. 'ERREUR' 5 ;
  1962. 'FINSI' ;
  1963.  
  1964. * RUX
  1965.  
  1966. DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
  1967. DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
  1968.  
  1969. RN1 = RN ;
  1970. GN1 = GN '+' DUX ;
  1971. RETN1 = RETN ;
  1972.  
  1973. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1974.  
  1975. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  1976. LISTC LISTP
  1977. RN1 VN1 PN1 GAMN CHPLI 'INJELM' ;
  1978.  
  1979. DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
  1980. DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  1981. 'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0) ;
  1982.  
  1983. ERRO = 'MAXIMUM' (DRESDUX '-' DRESDUXN) 'ABS' ;
  1984.  
  1985. 'SI' (ERRO > 1.0D-3) ;
  1986. 'ERREUR' 5 ;
  1987. 'FINSI' ;
  1988.  
  1989. * RUY
  1990.  
  1991. DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
  1992. DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
  1993.  
  1994. RN1 = RN ;
  1995. GN1 = GN '+' DUY ;
  1996. RETN1 = RETN ;
  1997.  
  1998. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  1999.  
  2000. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  2001. LISTC LISTP
  2002. RN1 VN1 PN1 GAMN CHPLI 'INJELM' ;
  2003.  
  2004. DRESDUYN = (RCHRES1 '-' RCHRES) '/' DUYSCAL ;
  2005. DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  2006. 'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0) ;
  2007.  
  2008. ERRO = 'MAXIMUM' (DRESDUY '-' DRESDUYN) 'ABS' ;
  2009.  
  2010. 'SI' (ERRO > 1.0D-3) ;
  2011. 'ERREUR' 5 ;
  2012. 'FINSI' ;
  2013.  
  2014.  
  2015. * RUZ
  2016.  
  2017. DUZSCAL = EPS '*' ('EXTRAIRE' GN 'UZ' PCEN) ;
  2018. DUZ = ('MANUEL' 'CHPO' ELE1 1 'UZ' DUZSCAL) ;
  2019.  
  2020. RN1 = RN ;
  2021. GN1 = GN '+' DUZ ;
  2022. RETN1 = RETN ;
  2023.  
  2024. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  2025.  
  2026. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  2027. LISTC LISTP
  2028. RN1 VN1 PN1 GAMN CHPLI 'INJELM' ;
  2029.  
  2030. DRESDUZN = (RCHRES1 '-' RCHRES) '/' DUZSCAL ;
  2031. DRESDUZ = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  2032. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0) ;
  2033.  
  2034. ERRO = 'MAXIMUM' (DRESDUZ '-' DRESDUZN) 'ABS' ;
  2035.  
  2036. 'SI' (ERRO > 1.0D-3) ;
  2037. 'ERREUR' 5 ;
  2038. 'FINSI' ;
  2039.  
  2040. * RETN
  2041.  
  2042. DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
  2043. ELE1 = 'MANUEL' 'POI1' PCEN ;
  2044. DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
  2045.  
  2046. RN1 = RN ;
  2047. GN1 = GN ;
  2048. RETN1 = RETN '+' DR ;
  2049.  
  2050. VN1 PN1 = 'PRIM' 'PERFMONO' RN1 GN1 RETN1 GAMN ;
  2051.  
  2052. RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMONO' 'CLIM' 'RESI' MDOMINT MENTREE
  2053. LISTC LISTP
  2054. RN1 VN1 PN1 GAMN CHPLI 'INJELM' ;
  2055.  
  2056. DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
  2057. DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 5 'VAR1' 0.0
  2058. 'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0) ;
  2059.  
  2060. ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
  2061.  
  2062. 'SI' (ERRO > 1.0D-3) ;
  2063. 'ERREUR' 5 ;
  2064. 'FINSI' ;
  2065.  
  2066. 'FIN' ;
  2067.  
  2068.  
  2069.  
  2070.  
  2071.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales