Télécharger pret_ther4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : pret_ther4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **********************************************************************
  5. **** APPROCHE VF "Cell-Centred Formulation" pour la solution des ****
  6. **** Equations d'Euler pour un gaz parfait. ****
  7. **** ****
  8. **** OPERATEUR PRET ****
  9. **** Operateur qui 'recontruit les variables primitives aux faces ****
  10. **** Cas gaz "thermally perfect" multi-especes ****
  11. **** Interieur et Mur ****
  12. **** Splitting des scalaires passifs ****
  13. **** ****
  14. **** A. BECCANTINI DRN/DMT/SEMT/TTMF FEVRIER 2000 ****
  15. **********************************************************************
  16.  
  17. 'OPTION' 'DIME' 2 ;
  18. 'OPTION' 'ELEM' QUA4 ;
  19. 'OPTION' 'ECHO' 0 ;
  20. 'OPTION' 'TRAC' 'X' ;
  21.  
  22. *
  23. *** GRAPH
  24. *
  25.  
  26. * GRAPH = VRAI ;
  27. GRAPH = FAUX ;
  28.  
  29.  
  30. ***************************
  31. *** PROCEDURE RUO1 ****
  32. ***************************
  33. *
  34. *** Precedure pour la rotation d'un tenseur de premier
  35. * ordre
  36. *
  37.  
  38. 'DEBPROC' RUO1 ;
  39. 'ARGUMENT' ALPHA*'FLOTTANT' UN*'FLOTTANT' UT*'FLOTTANT';
  40.  
  41. SINA = 'SIN' ALPHA ;
  42. COSA = 'COS' ALPHA ;
  43.  
  44. UX = (UN * COSA ) '-' (UT * SINA) ;
  45. UY = (UN * SINA ) '+' (UT * COSA) ;
  46.  
  47.  
  48. 'FINPROC' UX UY;
  49.  
  50. ***************************
  51. *** PROCEDURE RUO2 ****
  52. ***************************
  53. *
  54. *** Precedure pour la rotation d'un tenseur de deuxieme
  55. * ordre
  56. *
  57.  
  58. 'DEBPROC' RUO2 ;
  59. 'ARGUMENT' ALPHA*'FLOTTANT' UNN*'FLOTTANT' UNT*'FLOTTANT'
  60. UTN*'FLOTTANT' UTT*'FLOTTANT';
  61. *
  62. **** (n,t) -> (x,y)
  63. *
  64. * n = CA x '+' SA y ;
  65. * t = -SA x '+' CA y ;
  66. *
  67. * UNT = DUN/DT
  68. *
  69.  
  70. SA = 'SIN' ALPHA ;
  71. CA = 'COS' ALPHA ;
  72. CA2 = CA * CA ;
  73. CASA = CA * SA ;
  74. SA2 = SA * SA;
  75.  
  76. UXX = (CA2 * UNN) '-' (CASA * (UNT '+' UTN)) '+' (SA2 *UTT) ;
  77. UYX = (CASA * (UNN '-' UTT)) '-' (SA2 * UNT) '+' (CA2 * UTN ) ;
  78. UXY = (CASA * (UNN '-' UTT)) '+' (CA2 * UNT) '-' (SA2 * UTN ) ;
  79. UYY = (SA2 * UNN) '+' (CASA * (UNT '+' UTN)) '+' (CA2 *UTT) ;
  80.  
  81. 'FINPROC' UXX UXY UYX UYY;
  82.  
  83. *******************************************************
  84. **** Cas multi-especes : la table proprieté de gaz ****
  85. *******************************************************
  86.  
  87. PGAZ = 'TABLE' ;
  88. PGAZ2= 'TABLE' ;
  89.  
  90. *
  91. **** Especes qui sont dans les equations d'Euler
  92. *
  93.  
  94. PGAZ . 'ESPEULE' = 'MOTS' 'H2 ' 'N2 ' ;
  95. PGAZ2. 'ESPEULE' = 'MOTS' 'H2 ' 'N2 ' ;
  96.  
  97. *
  98. **** Espece qui n'y est pas
  99. *
  100.  
  101.  
  102. PGAZ . 'ESPNEULE' = 'O2 ';
  103. PGAZ2. 'ESPNEULE' = 'O2 ';
  104.  
  105. *
  106.  
  107. PGAZ . 'H2 ' = 'TABLE' ;
  108. PGAZ . 'N2 ' = 'TABLE' ;
  109. PGAZ . 'O2 ' = 'TABLE' ;
  110. PGAZ2. 'H2 ' = 'TABLE' ;
  111. PGAZ2. 'N2 ' = 'TABLE' ;
  112. PGAZ2. 'O2 ' = 'TABLE' ;
  113.  
  114. *
  115. **** R (J/Kg/K)
  116. *
  117.  
  118. PGAZ . 'H2 ' . 'R' = 4130.0 ;
  119. PGAZ . 'N2 ' . 'R' = 296.8 ;
  120. PGAZ . 'O2 ' . 'R' = 259.8 ;
  121. PGAZ2. 'H2 ' . 'R' = 4130.0 ;
  122. PGAZ2. 'N2 ' . 'R' = 296.8 ;
  123. PGAZ2. 'O2 ' . 'R' = 259.8 ;
  124.  
  125.  
  126. *
  127. **** Regressions polynomials
  128. *
  129.  
  130. PGAZ . 'H2 ' . 'A' = 'PROG' 9834.91866 0.0 0.0 0.0 0.0 ;
  131. PGAZ . 'O2 ' . 'A' = 'PROG' 575.012333 0.0 0.0 0.0 0. ;
  132. PGAZ . 'N2 ' . 'A' = 'PROG' 652.940766 0.0 0.0 0.0 0.0 ;
  133. PGAZ2. 'H2 ' . 'A' = 'PROG' 9834.91866 0.0 0.0 0.0 0.0 ;
  134. PGAZ2. 'O2 ' . 'A' = 'PROG' 575.012333 0.0 0.0 0.0 0. ;
  135. PGAZ2. 'N2 ' . 'A' = 'PROG' 652.940766 0.0 0.0 0.0 0.0 ;
  136.  
  137. *
  138. **** "Enthalpies" (ou energies) de formations a OK (J/Kg)
  139. * Note: ce sont des entites numeriques
  140. * h_i = h_i(T0) '-' \int_0^{T0} cp_i(x) dx =
  141. * h_i(T0) '-' (\int_0^{T0} cv_i(x) dx '+' R_i * T0)
  142. *
  143.  
  144. PGAZ . 'H2 ' . 'H0K' = -4.195D6 ;
  145. PGAZ . 'O2 ' . 'H0K' = -2.634D5 ;
  146. PGAZ . 'N2 ' . 'H0K' = -2.953D5 ;
  147. PGAZ2. 'H2 ' . 'H0K' = -4.195D6 ;
  148. PGAZ2. 'O2 ' . 'H0K' = -2.634D5 ;
  149. PGAZ2. 'N2 ' . 'H0K' = -2.953D5 ;
  150.  
  151. PGAZ2 . 'SCALPASS' = 'MOTS' 'TITI' 'TUTU' 'TATA' ;
  152.  
  153. ***************************
  154. ***** DOMAINE SPATIAL ****
  155. ***************************
  156.  
  157. *
  158. **** Deux carre
  159. *
  160.  
  161. A1 = 0.0D0 0.0D0 ;
  162. A2 = 1.0D0 0.0D0 ;
  163. A3 = 2.0D0 0.0D0 ;
  164. A4 = 2.0D0 1.0D0 ;
  165. A5 = 1.0D0 1.0D0 ;
  166. A6 = 0.0D0 1.0D0 ;
  167.  
  168. L12 = A1 'DROIT' 1 A2 ;
  169. L23 = A2 'DROIT' 1 A3 ;
  170. L34 = A3 'DROIT' 1 A4 ;
  171. L45 = A4 'DROIT' 1 A5 ;
  172. L56 = A5 'DROIT' 1 A6 ;
  173. L61 = A6 'DROIT' 1 A1 ;
  174. L25 = A2 'DROIT' 1 A5 ;
  175.  
  176.  
  177. DOM10 = 'DALL' L12 L25 L56 L61
  178. 'PLANE';
  179. DOM20 = 'DALL' L23 L34 L45 ('INVERSE' L25)
  180. 'PLANE';
  181.  
  182. *
  183. **** P10: Point face entre le deux carre, ou on fait les controlles
  184. **** P20: Point face sur le "mur"
  185. **
  186.  
  187. P10 = 1.0 0.5 ;
  188. P20 = 2.0 0.5 ;
  189.  
  190. DOM1 = DOM10 ;
  191. DOM2 = DOM20 ;
  192.  
  193.  
  194. DOMTOT = DOM1 ET DOM2;
  195. 'ELIMINATION' (DOM1 ET DOM2) 1D-6;
  196. DOMTOT = DOM1 ET DOM2;
  197.  
  198.  
  199. $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  200.  
  201. $DOM1 = 'MODELISER' DOM1 'EULER';
  202. $DOM2 = 'MODELISER' DOM2 'EULER';
  203.  
  204. TDOMTOT = 'DOMA' $DOMTOT 'VF';
  205.  
  206. TDOM1 = 'DOMA' $DOM1 'VF';
  207. TDOM2 = 'DOMA' $DOM2 'VF';
  208.  
  209. MDOM1 = TDOM1 . 'QUAF' ;
  210. MDOM2 = TDOM2 . 'QUAF' ;
  211. MDOMTOT = TDOMTOT . 'QUAF' ;
  212.  
  213. 'ELIMINATION' (MDOMTOT ET MDOM1) 0.0001 ;
  214. 'ELIMINATION' (MDOMTOT ET MDOM2) 0.0001 ;
  215.  
  216. 'SI' GRAPH;
  217. 'TRACER' (('DOMA' $DOMTOT 'MAILLAGE') 'ET' ('DOMA' $DOMTOT 'FACEL')
  218. 'ET' P10 'ET' P20) 'TITRE' 'Domaine et FACEL';
  219. 'FINSI' ;
  220.  
  221.  
  222. *
  223. **** TEST1: si les limiteurs sont nuls ou les gradients son nuls
  224. * le deuxieme ordre en espace 'et premier ordre
  225. * en temps degenere en premier ordre en espace
  226. *
  227.  
  228.  
  229. RN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 1.0D0;
  230.  
  231. PN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 2.0D0;
  232.  
  233. VIT = 'KCHT' $DOMTOT 'VECT' 'CENTRE' (4.0 5.0);
  234.  
  235. GAMMA = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 1.5D0;
  236.  
  237.  
  238. *
  239. **** Limiteurs nuls, gradients non-nuls
  240. *
  241.  
  242.  
  243.  
  244. GRADR = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  245. (1.0D0 1.0D0);
  246.  
  247.  
  248.  
  249. GRADP = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  250. (1.0D0 1.0D0);
  251.  
  252.  
  253.  
  254. GRADVX = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  255. (1.0D0 1.0D0);
  256.  
  257.  
  258. GRADVY = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  259. (1.0D0 1.0D0);
  260.  
  261. GRADV = GRADVX 'ET' GRADVY ;
  262.  
  263.  
  264.  
  265. ALR = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1' 0.0D0 ;
  266.  
  267.  
  268. ALP = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1' 0.0D0 ;
  269.  
  270.  
  271. ALV = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  272. (0.0 0.0D0) ;
  273.  
  274. * YN
  275.  
  276.  
  277. YG = 'PROG' 0.2 0.1 ;
  278. YD = 'PROG' 0.0 0.13 ;
  279.  
  280. YN1 = 'KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' ('MOTS' 'H2 ' 'N2 ' )
  281. (('EXTRAIRE' 1 YG) ('EXTRAIRE' 2 YG)) ;
  282. YN2 = 'KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' ('MOTS' 'H2 ' 'N2 ')
  283. (('EXTRAIRE' 1 YD) ('EXTRAIRE' 2 YD)) ;
  284. YN = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' ('MOTS' 'H2 ' 'N2 ')
  285. (YN1 'ET' YN2) ;
  286.  
  287. GRADYH2 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  288. (1.0D0 1.0D0);
  289.  
  290. GRADYN2 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  291. (1.0D0 1.0D0);
  292. GRADYN = GRADYH2 'ET' GRADYN2 ;
  293.  
  294. ALYN = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1' 'P2' (0.0 0.0) ;
  295.  
  296. *
  297.  
  298. scag = 11.342 ;
  299. scad = 34.2;
  300.  
  301. SCAC11 = 'KCHT' $DOM1 'SCAL' 'CENTRE' 'COMP' 'TITI' scag ;
  302. SCAC12 = 'KCHT' $DOM2 'SCAL' 'CENTRE' 'COMP' 'TITI' scad ;
  303. SCAC1 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'TITI'
  304. (SCAC11 'ET' SCAC12) ;
  305.  
  306. SCAC21 = 'KCHT' $DOM1 'SCAL' 'CENTRE' 'COMP' 'TUTU' scag ;
  307. SCAC22 = 'KCHT' $DOM2 'SCAL' 'CENTRE' 'COMP' 'TUTU' scad ;
  308. SCAC2 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'TUTU'
  309. (SCAC21 'ET' SCAC22) ;
  310.  
  311. SCAC31 = 'KCHT' $DOM1 'SCAL' 'CENTRE' 'COMP' 'TATA' scag ;
  312. SCAC32 = 'KCHT' $DOM2 'SCAL' 'CENTRE' 'COMP' 'TATA' scad ;
  313. SCAC3 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'TATA'
  314. (SCAC31 'ET' SCAC32) ;
  315.  
  316.  
  317. SCAC = SCAC1 'ET' SCAC2 'ET' SCAC3 ;
  318.  
  319.  
  320. GRADS1 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  321. (1.0D0 1.0D0);
  322. GRADS2 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  323. (1.0D0 1.0D0);
  324. GRADS3 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P3DX' 'P3DY'
  325. (1.0D0 1.0D0);
  326.  
  327. GRADS = GRADS1 'ET' GRADS2 'ET' GRADS3 ;
  328.  
  329. ALS1 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1' 0.0 ;
  330. ALS2 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P2' 0.0 ;
  331. ALS3 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P3' 0.0 ;
  332.  
  333. ALS = ALS1 'ET' ALS2 'ET' ALS3 ;
  334.  
  335. *
  336.  
  337. ORDESP = 1;
  338. ORDTEM = 1;
  339.  
  340.  
  341. ROF1 VITF1 PF1 YF1 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  342. $DOMTOT PGAZ RN VIT PN YN ;
  343.  
  344. ROF11 VITF11 PF11 YF11 SCAF11 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  345. $DOMTOT PGAZ2 RN VIT PN YN SCAC ;
  346.  
  347.  
  348. ORDESP = 2;
  349. ORDTEM = 1;
  350.  
  351.  
  352. ROF VITF PF YF = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  353. $DOMTOT PGAZ RN GRADR ALR
  354. VIT GRADV ALV
  355. PN GRADP ALP
  356. YN GRADYN ALYN;
  357.  
  358. ROF22 VITF22 PF22 YF22 SCAF22 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  359. $DOMTOT PGAZ2 RN GRADR ALR
  360. VIT GRADV ALV
  361. PN GRADP ALP
  362. YN GRADYN ALYN
  363. SCAC GRADS ALS;
  364.  
  365.  
  366. ERRO = 'MAXIMUM' ('PROG'
  367. ('MAXIMUM' (ROF '-' ROF1) 'ABS') ('MAXIMUM' (ROF '-' ROF11) 'ABS')
  368. ('MAXIMUM' (ROF '-' ROF22) 'ABS')
  369. ('MAXIMUM' (PF '-' PF1) 'ABS') ('MAXIMUM' (PF '-' PF11) 'ABS')
  370. ('MAXIMUM' (PF '-' PF22) 'ABS')
  371. ('MAXIMUM' (VITF '-' VITF1) 'ABS') ('MAXIMUM'(VITF '-' VITF11)'ABS')
  372. ('MAXIMUM' (VITF '-' VITF22) 'ABS') ) ;
  373.  
  374. ERRO = 'MAXIMUM' ('PROG' ERRO
  375. ('MAXIMUM' (YF '-' YF1) 'ABS') ('MAXIMUM' (YF '-' YF11) 'ABS')
  376. ('MAXIMUM' (YF '-' YF22) 'ABS') ('MAXIMUM' (SCAF11 '-' SCAF22) 'ABS'));
  377.  
  378.  
  379. 'SI' (ERRO > 1.0D-12);
  380. 'ERREUR' 5;
  381. 'FINSI' ;
  382.  
  383.  
  384. *
  385. **** Limiteurs non-nuls, gradients nuls
  386. *
  387.  
  388.  
  389. GRADR = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  390. (0.0D0 0.0D0);
  391.  
  392. GRADP = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  393. (0.0D0 0.0D0);
  394.  
  395. GRADVX = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  396. (0.0D0 0.0D0);
  397. GRADVY = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  398. (0.0D0 0.0D0);
  399. GRADV = GRADVX 'ET' GRADVY ;
  400.  
  401. ALR = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1' 0.4D0 ;
  402.  
  403.  
  404. ALP = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1' 0.8D0 ;
  405.  
  406.  
  407. ALV = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  408. (0.9 0.9D0) ;
  409.  
  410. GRADYN = ('KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  411. (0.0D0 0.0D0)) 'ET'
  412. ('KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  413. (0.0D0 0.0D0));
  414.  
  415. ALYN = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1' 'P2 ' (1.0 1.0) ;
  416.  
  417.  
  418.  
  419. GRADS1 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  420. (0.0D0 0.0D0);
  421. GRADS2 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  422. (0.0D0 0.0D0);
  423. GRADS3 = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P3DX' 'P3DY'
  424. (0.0D0 0.0D0);
  425.  
  426. GRADS = GRADS1 'ET' GRADS2 'ET' GRADS3 ;
  427.  
  428. ALS1 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1' 1.0 ;
  429. ALS2 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P2' 1.0 ;
  430. ALS3 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P3' 1.0 ;
  431.  
  432. ALS = ALS1 'ET' ALS2 'ET' ALS3 ;
  433.  
  434. *
  435.  
  436.  
  437. ORDESP = 1;
  438. ORDTEM = 1;
  439.  
  440.  
  441. ROF1 VITF1 PF1 YF1 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  442. $DOMTOT PGAZ RN VIT PN YN ;
  443.  
  444. ROF11 VITF11 PF11 YF11 SCAF11 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  445. $DOMTOT PGAZ2 RN VIT PN YN SCAC ;
  446.  
  447.  
  448. ORDESP = 2;
  449. ORDTEM = 1;
  450.  
  451.  
  452. ROF VITF PF YF = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  453. $DOMTOT PGAZ RN GRADR ALR
  454. VIT GRADV ALV
  455. PN GRADP ALP
  456. YN GRADYN ALYN;
  457.  
  458. ROF22 VITF22 PF22 YF22 SCAF22 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  459. $DOMTOT PGAZ2 RN GRADR ALR
  460. VIT GRADV ALV
  461. PN GRADP ALP
  462. YN GRADYN ALYN
  463. SCAC GRADS ALS;
  464.  
  465.  
  466. ERRO = 'MAXIMUM' ('PROG'
  467. ('MAXIMUM' (ROF '-' ROF1) 'ABS') ('MAXIMUM' (ROF '-' ROF11) 'ABS')
  468. ('MAXIMUM' (ROF '-' ROF22) 'ABS')
  469. ('MAXIMUM' (PF '-' PF1) 'ABS') ('MAXIMUM' (PF '-' PF11) 'ABS')
  470. ('MAXIMUM' (PF '-' PF22) 'ABS')
  471. ('MAXIMUM' (VITF '-' VITF1) 'ABS') ('MAXIMUM'(VITF '-' VITF11)'ABS')
  472. ('MAXIMUM' (VITF '-' VITF22) 'ABS') ) ;
  473.  
  474. ERRO = 'MAXIMUM' ('PROG' ERRO
  475. ('MAXIMUM' (YF '-' YF1) 'ABS') ('MAXIMUM' (YF '-' YF11) 'ABS')
  476. ('MAXIMUM' (YF '-' YF22) 'ABS') ('MAXIMUM' (SCAF11 '-' SCAF22) 'ABS'));
  477.  
  478.  
  479. 'SI' (ERRO > 1.0D-12);
  480. 'ERREUR' 5;
  481. 'FINSI' ;
  482.  
  483.  
  484. *
  485. **** TEST2 : Deuxieme ordre en espace, deuxieme ordre en
  486. * temps: invariance par rotation
  487. *
  488. * Il faut remarquer que:
  489. * densite, pression sont des tenseurs d'ordre 0,
  490. * donc leurs gradients sont des tenseurs d'ordre 1
  491. * la vitesse est un tenseur d'ordre 1,
  492. * donc sont gradient est un tenseur d'ordre 2
  493. *
  494. **** On considere un repere solidal avec la surface, (n,t),
  495. * et un repere (x,y)
  496. * On consider une probleme 1D dans le repaire (n,t)
  497. *
  498. * Dans le cas de transort de scalaires passifs
  499. * (PGAZ2 . 'SCALPASS' = 'MOTS' 'TITI' 'TUTU' 'TATA')
  500. * on impose TITI = N2, TOTO = H2, TATA = N2
  501. *
  502. *** Etats gauche et droite
  503. *
  504.  
  505.  
  506.  
  507. gamg = 1.4D0;
  508. rog = 1.00 ;
  509. grong = -1.0D0 ;
  510. grotg = 0.0D0;
  511. alrog = 0.1 ;
  512. pg = 2.0D0 ;
  513. gpng = -3.0D0 ;
  514. gptg = 0.0D0;
  515. alpg = 0.2 ;
  516. ung = 3.0D0 ;
  517. gunng = -5.0D0 ;
  518. guntg = 0.0D0 ;
  519. utg = 0.0D0 ;
  520. gutng = 0.0D0 ;
  521. guttg = 0.0D0;
  522. alug = 0.3 ;
  523. yh2g = 0.1 ;
  524. gyh2ng = -0.0111 ;
  525. gyh2tg = 0.0 ;
  526. alyh2g = 0.123 ;
  527. yn2g = 0.34 ;
  528. gyn2ng = -0.0311 ;
  529. gyn2tg = 0.0 ;
  530. alyn2g = 0.323 ;
  531.  
  532. c2g = gamg '*' pg '/' rog ;
  533.  
  534.  
  535. gamd = 1.8 ;
  536. rod = 1.50 ;
  537. grond = -1.5D0 ;
  538. grotd = 0.0D0 ;
  539. alrod = 0.25 ;
  540. pd = 2.5D0 ;
  541. gpnd = -3.5D0 ;
  542. gptd = 0.0D0 ;
  543. alpd = 0.3 ;
  544. und = 3.5D0 ;
  545. gunnd = -5.5D0 ;
  546. guntd = 0.0D0 ;
  547. utd = 0.0D0 ;
  548. gutnd = 0.0D0 ;
  549. guttd = 0.0D0;
  550. alud = 0.15 ;
  551. yh2d = 0.1131 ;
  552. gyh2nd = -0.03111 ;
  553. gyh2td = 0.0 ;
  554. alyh2d = 0.4123 ;
  555. yn2d = 0.134 ;
  556. gyn2nd = -0.01311 ;
  557. gyn2td = 0.0 ;
  558. alyn2d = 0.1323 ;
  559.  
  560. c2d = gamd '*' pd '/' rod ;
  561.  
  562. deltat = 1.0 ;
  563.  
  564. * Encrement predictive
  565.  
  566. drog = -1.0 '*' ((ung '*' grong '*' alrog) '+' (rog '*' gunng '*' alug))
  567. ;
  568. drog = drog '*' deltat ;
  569.  
  570. drod = -1.0 '*' ((und '*' grond '*' alrod) '+' (rod '*' gunnd '*' alud))
  571. ;
  572. drod = drod '*' deltat ;
  573.  
  574. dung = -1.0 '*' ((ung '*' gunng '*' alug) '+' (alpg '*' gpng '/' rog)) ;
  575. dung = dung '*' deltat ;
  576.  
  577. dund = -1.0 '*' ((und '*' gunnd '*' alud) '+' (alpd '*' gpnd '/' rod)) ;
  578. dund = dund '*' deltat ;
  579.  
  580. dpg = -1.0 '*' ((rog '*' c2g '*' gunng '*' alug) '+'
  581. (ung '*' gpng '*' alpg)) ;
  582. dpg = dpg '*' deltat ;
  583.  
  584. dpd = -1.0 '*' ((rod '*' c2d '*' gunnd '*' alud) '+'
  585. (und '*' gpnd '*' alpd)) ;
  586. dpd = dpd '*' deltat ;
  587.  
  588. dyh2g = -1.0 '*' ((ung '*' gyh2ng '*' alyh2g)) ;
  589. dyh2g = dyh2g '*' deltat ;
  590.  
  591. dyh2d = -1.0 '*' ((und '*' gyh2nd '*' alyh2d)) ;
  592. dyh2d = dyh2d '*' deltat ;
  593.  
  594. dyn2g = -1.0 '*' ((ung '*' gyn2ng '*' alyn2g)) ;
  595. dyn2g = dyn2g '*' deltat ;
  596.  
  597. dyn2d = -1.0 '*' ((und '*' gyn2nd '*' alyn2d)) ;
  598. dyn2d = dyn2d '*' deltat ;
  599.  
  600. * Encrement corrective; dxg = -dxd = 0.5 (point P1)
  601.  
  602. dsrog = 0.5 '*' grong '*' alrog ;
  603. dsrod = -0.5 '*' grond '*' alrod ;
  604.  
  605. dspg = 0.5 '*' gpng '*' alpg ;
  606. dspd = -0.5 '*' gpnd '*' alpd ;
  607.  
  608. dsung = 0.5 '*' gunng '*' alug ;
  609. dsund = -0.5 '*' gunnd '*' alud ;
  610.  
  611. dsyh2g = 0.5 '*' gyh2ng '*' alyh2g ;
  612. dsyh2d = -0.5 '*' gyh2nd '*' alyh2d ;
  613.  
  614. dsyn2g = 0.5 '*' gyn2ng '*' alyn2g ;
  615. dsyn2d = -0.5 '*' gyn2nd '*' alyn2d ;
  616.  
  617. * Encrement corrective; dxd = 0.5 (point P2)
  618.  
  619.  
  620. dsrod2 = -1.0 '*' dsrod ;
  621. dspd2 = -1.0 '*' dspd ;
  622. dsund2 = -1.0 '*' dsund ;
  623. dsyh2d2 = -1.0 '*' dsyh2d ;
  624. dsyn2d2 = -1.0 '*' dsyn2d ;
  625.  
  626. * Les etats (point P1)
  627.  
  628. etrog = rog '+' drog '+' dsrog ;
  629. etrod = rod '+' drod '+' dsrod ;
  630.  
  631. etung = ung '+' dung '+' dsung ;
  632. etund = und '+' dund '+' dsund ;
  633.  
  634. etpg = pg '+' dpg '+' dspg ;
  635. etpd = pd '+' dpd '+' dspd ;
  636.  
  637. etyh2g = yh2g '+' dyh2g '+' dsyh2g ;
  638. etyh2d = yh2d '+' dyh2d '+' dsyh2d ;
  639.  
  640. etyn2g = yn2g '+' dyn2g '+' dsyn2g ;
  641. etyn2d = yn2d '+' dyn2d '+' dsyn2d ;
  642.  
  643. * Les etats (point P2)
  644.  
  645. etrod2 = rod '+' drod '+' dsrod2 ;
  646. etund2 = und '+' dund '+' dsund2 ;
  647. etpd2 = pd '+' dpd '+' dspd2 ;
  648. etyh2d2 = yh2d '+' dyh2d '+' dsyh2d2 ;
  649. etyn2d2 = yn2d '+' dyn2d '+' dsyn2d2 ;
  650.  
  651.  
  652. *
  653. *** Etats a gauche et a droite du point face P1
  654. *
  655. * ro, un, ut, p, yh2, yn2
  656.  
  657. ETATG = 'PROG' etrog etung 0.0 etpg etyh2g etyn2g ;
  658. ETATD = 'PROG' etrod etund 0.0 etpd etyh2d etyn2d ;
  659.  
  660. *
  661. *** Etats a gauche et a droite du point face P1
  662. *
  663. * ro, un, ut, p, yh2, yn2
  664.  
  665. ETATG2 = 'PROG' etrod2 etund2 0.0 etpd2 etyh2d2 etyn2d2 ;
  666. ETATD2 = 'PROG' etrod2 (-1.0 '*'etund2) 0.0 etpd2 etyh2d2 etyn2d2 ;
  667.  
  668. *
  669. *** Rotation
  670. *
  671.  
  672. DANGLE = 85;
  673. ANGLE = -1.0 '*' DANGLE ;
  674. ORIG = 0.0D0 0.0D0;
  675.  
  676.  
  677. 'REPETER' BL1 5;
  678.  
  679. ANGLE = ANGLE '+' DANGLE ;
  680.  
  681. uxg uyg = RUO1 ANGLE UNG UTG ;
  682. uxd uyd = RUO1 ANGLE UND UTD ;
  683.  
  684. groxg groyg = RUO1 ANGLE grong grotg;
  685. groxd groyd = RUO1 ANGLE grond grotd;
  686.  
  687. gpxg gpyg = RUO1 ANGLE gpng gptg;
  688. gpxd gpyd = RUO1 ANGLE gpnd gptd;
  689.  
  690. guxxg guxyg guyxg guyyg = RUO2 ANGLE gunng guntg gutng guttg;
  691. guxxd guxyd guyxd guyyd = RUO2 ANGLE gunnd guntd gutnd guttd;
  692.  
  693. gyh2xg gyh2yg = RUO1 ANGLE gyh2ng gyh2tg ;
  694. gyh2xd gyh2yd = RUO1 ANGLE gyh2nd gyh2td ;
  695.  
  696. gyn2xg gyn2yg = RUO1 ANGLE gyn2ng gyn2tg ;
  697. gyn2xd gyn2yd = RUO1 ANGLE gyn2nd gyn2td ;
  698.  
  699. 'MESSAGE' ;
  700. 'MESSAGE' (CHAIN 'Angle de rotation= ' ANGLE);
  701. 'MESSAGE' ;
  702.  
  703. DOM1 = DOM10 'TOURNER' ANGLE ORIG;
  704. DOM2 = DOM20 'TOURNER' ANGLE ORIG;
  705. P1 = P10 'TOURNER' ANGLE ORIG;
  706. P2 = P20 'TOURNER' ANGLE ORIG;
  707.  
  708. DOMTOT = DOM1 ET DOM2;
  709. 'ELIMINATION' (DOM1 ET DOM2) 1D-6;
  710. DOMTOT = DOM1 ET DOM2;
  711.  
  712.  
  713. $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  714.  
  715. $DOM1 = 'MODELISER' DOM1 'EULER';
  716. $DOM2 = 'MODELISER' DOM2 'EULER';
  717.  
  718. TDOMTOT = 'DOMA' $DOMTOT 'VF';
  719.  
  720. TDOM1 = 'DOMA' $DOM1 'VF';
  721. TDOM2 = 'DOMA' $DOM2 'VF';
  722.  
  723. MDOM1 = TDOM1 . 'QUAF' ;
  724. MDOM2 = TDOM2 . 'QUAF' ;
  725. MDOMTOT = TDOMTOT . 'QUAF' ;
  726.  
  727. 'ELIMINATION' (MDOMTOT ET MDOM1) 0.0001 ;
  728. 'ELIMINATION' (MDOMTOT ET MDOM2) 0.0001 ;
  729.  
  730. 'SI' GRAPH;
  731. 'TRACER' (('DOMA' $DOMTOT 'MAILLAGE') 'ET' ('DOMA' $DOMTOT 'FACEL')
  732. 'ET' P1 'ET' P2) 'TITRE' 'Domaine et FACEL';
  733. 'FINSI' ;
  734.  
  735. *
  736. **** Redefinition de P1 dans $DOMTOT . 'FACE'
  737. *
  738.  
  739. P1 = ('DOMA' $DOMTOT 'FACE') 'POIN' 'PROC' P1;
  740. P2 = ('DOMA' $DOMTOT 'FACE') 'POIN' 'PROC' P2;
  741.  
  742.  
  743. ***********************
  744. **** Les CHPOINTs ****
  745. ***********************
  746.  
  747. GAM1 = 'KCHT' $DOM1 'SCAL' 'CENTRE' gamg ;
  748. GAM2 = 'KCHT' $DOM2 'SCAL' 'CENTRE' gamd ;
  749. GAMMA = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' (GAM1 'ET' GAM2) ;
  750.  
  751. RN1 = 'KCHT' $DOM1 'SCAL' 'CENTRE' rog;
  752. RN2 = 'KCHT' $DOM2 'SCAL' 'CENTRE' rod;
  753. RN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' (RN1 'ET' RN2);
  754.  
  755.  
  756. VIT1 = 'KCHT' $DOM1 'VECT' 'CENTRE' (uxg uyg);
  757. VIT2 = 'KCHT' $DOM2 'VECT' 'CENTRE' (uxd uyd);
  758. VIT = 'KCHT' $DOMTOT 'VECT' 'CENTRE' (VIT1 'ET' VIT2);
  759.  
  760. PN1 = 'KCHT' $DOM1 'SCAL' 'CENTRE' pg;
  761. PN2 = 'KCHT' $DOM2 'SCAL' 'CENTRE' pd;
  762. PN = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' (PN1 'ET' PN2);
  763.  
  764. YN1 = 'KCHT' $DOM1 'VECT' 'CENTRE'
  765. 'COMP' 'H2 ' 'N2 ' (yh2g yn2g) ;
  766. YN2 = 'KCHT' $DOM2 'VECT' 'CENTRE'
  767. 'COMP' 'H2 ' 'N2 ' (yh2d yn2d) ;
  768. YN = YN1 'ET' YN2 ;
  769.  
  770. *
  771. **** On impose les gradients et le limiteurs
  772. *
  773.  
  774. ALR1 = 'KCHT' $DOM1 'SCAL' 'CENTRE' 'COMP' 'P1' alrog ;
  775. ALR2 = 'KCHT' $DOM2 'SCAL' 'CENTRE' 'COMP' 'P1' alrod ;
  776. ALR = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1'
  777. (ALR1 'ET' ALR2) ;
  778.  
  779.  
  780. ALP1 = 'KCHT' $DOM1 'SCAL' 'CENTRE' 'COMP' 'P1' alpg ;
  781. ALP2 = 'KCHT' $DOM2 'SCAL' 'CENTRE' 'COMP' 'P1' alpd ;
  782. ALP = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'P1'
  783. (ALP1 'ET' ALP2) ;
  784.  
  785.  
  786. ALV1 = 'KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  787. (alug alug) ;
  788. ALV2 = 'KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  789. (alud alud) ;
  790. ALV = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  791. (ALV1 'ET' ALV2) ;
  792.  
  793.  
  794. ALY1 = 'KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  795. (alyh2g alyn2g) ;
  796. ALY2 = 'KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  797. (alyh2d alyn2d) ;
  798. ALYN = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1' 'P2'
  799. (ALY1 'ET' ALY2) ;
  800.  
  801.  
  802. GRADR1 = 'KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  803. (groxg groyg) ;
  804.  
  805. GRADR2 = 'KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  806. (groxd groyd) ;
  807.  
  808. GRADR = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  809. (GRADR1 'ET' GRADR2);
  810.  
  811.  
  812.  
  813. GRADP1 = 'KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  814. (gpxg gpyg) ;
  815.  
  816. GRADP2 = 'KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  817. (gpxd gpyd) ;
  818.  
  819. GRADP = 'KCHT' $DOMTOT 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  820. (GRADP1 'ET' GRADP2);
  821.  
  822.  
  823.  
  824. GRADVX1 = 'KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  825. (guxxg guxyg);
  826.  
  827. GRADVX2 = 'KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  828. (guxxd guxyd);
  829.  
  830. GRADVY1 = 'KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  831. (guyxg guyyg);
  832.  
  833. GRADVY2 = 'KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  834. (guyxd guyyd);
  835.  
  836. GRADVX = GRADVX1 'ET' GRADVX2 ;
  837. GRADVY = GRADVY1 'ET' GRADVY2 ;
  838. GRADV = GRADVX 'ET' GRADVY ;
  839.  
  840. GRADYN = ('KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  841. (gyh2xg gyh2yg )) 'ET'
  842. ('KCHT' $DOM1 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  843. (gyn2xg gyn2yg )) 'ET'
  844. ('KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P1DX' 'P1DY'
  845. (gyh2xd gyh2yd )) 'ET'
  846. ('KCHT' $DOM2 'VECT' 'CENTRE' 'COMP' 'P2DX' 'P2DY'
  847. (gyn2xd gyn2yd ));
  848.  
  849.  
  850. *N.B: PGAZ2 . 'SCALPASS' = 'MOTS' 'TITI' 'TUTU' 'TATA'
  851. * PGAZ2 . 'ESPEULE' = 'MOTS' 'H2' 'N2'
  852.  
  853.  
  854. SCAC = 'EXCO' ('MOTS' 'N2' 'H2' 'N2') YN
  855. ('MOTS' 'TITI' 'TUTU' 'TATA') 'NATURE' 'DISCRET' ;
  856.  
  857. GRADS = 'EXCO' ('MOTS' 'P2DX' 'P2DY' 'P1DX' 'P1DY' 'P2DX' 'P2DY')
  858. GRADYN
  859. ('MOTS' 'P1DX' 'P1DY' 'P2DX' 'P2DY' 'P3DX' 'P3DY' )
  860. 'NATURE' 'DISCRET' ;
  861.  
  862. ALS = 'EXCO' ('MOTS' 'P2' 'P1' 'P2') ALYN ('MOTS' 'P1' 'P2' 'P3')
  863. 'NATURE' 'DISCRET' ;
  864.  
  865.  
  866.  
  867.  
  868.  
  869. ORDESP = 2;
  870. ORDTEM = 1;
  871.  
  872. ROF VITF PF YF = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  873. $DOMTOT PGAZ RN GRADR ALR
  874. VIT GRADV ALV
  875. PN GRADP ALP
  876. YN GRADYN ALYN ;
  877.  
  878.  
  879. ROF1 VITF1 PF1 YF1 SCAF1 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  880. $DOMTOT PGAZ2 RN GRADR ALR
  881. VIT GRADV ALV
  882. PN GRADP ALP
  883. YN GRADYN ALYN
  884. SCAC GRADS ALS ;
  885.  
  886. YF11 = 'EXCO' 'H2' YF1 'SCAL' ;
  887. YF12 = 'EXCO' 'N2' YF1 'SCAL' ;
  888.  
  889. SCAF11 = 'EXCO' 'TITI' SCAF1 'SCAL' ;
  890. SCAF12 = 'EXCO' 'TUTU' SCAF1 'SCAL' ;
  891. SCAF13 = 'EXCO' 'TATA' SCAF1 'SCAL' ;
  892.  
  893. ERRO = 'MAXIMUM' ('PROG'
  894. ('MAXIMUM' (ROF '-' ROF) 'ABS')
  895. ('MAXIMUM' (PF '-' PF1) 'ABS')
  896. ('MAXIMUM' (VITF '-' VITF1) 'ABS')
  897. ('MAXIMUM' (YF '-' YF1) 'ABS') );
  898.  
  899. ERRO = 'MAXIMUM' ('PROG' ERRO
  900. ('MAXIMUM' (YF11 '-' SCAF12) 'ABS')
  901. ('MAXIMUM' (YF12 '-' SCAF11) 'ABS')
  902. ('MAXIMUM' (YF12 '-' SCAF13) 'ABS')) ;
  903.  
  904.  
  905. 'SI' (ERRO > 1.0D-12);
  906. 'ERREUR' 5;
  907. 'FINSI' ;
  908.  
  909.  
  910. *
  911.  
  912. ORDESP = 2;
  913. ORDTEM = 2;
  914.  
  915. ROF VITF PF YF = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  916. $DOMTOT PGAZ RN GRADR ALR
  917. VIT GRADV ALV
  918. PN GRADP ALP
  919. YN GRADYN ALYN
  920. GAMMA deltat ;
  921.  
  922.  
  923. ROF1 VITF1 PF1 YF1 SCAF1 = 'PRET' 'PERFTEMP' ORDESP ORDTEM
  924. $DOMTOT PGAZ2 RN GRADR ALR
  925. VIT GRADV ALV
  926. PN GRADP ALP
  927. YN GRADYN ALYN
  928. SCAC GRADS ALS
  929. GAMMA deltat ;
  930.  
  931. YF11 = 'EXCO' 'H2' YF1 'SCAL' ;
  932. YF12 = 'EXCO' 'N2' YF1 'SCAL' ;
  933.  
  934. SCAF11 = 'EXCO' 'TITI' SCAF1 'SCAL' ;
  935. SCAF12 = 'EXCO' 'TUTU' SCAF1 'SCAL' ;
  936. SCAF13 = 'EXCO' 'TATA' SCAF1 'SCAL' ;
  937.  
  938. ERRO = 'MAXIMUM' ('PROG'
  939. ('MAXIMUM' (ROF '-' ROF) 'ABS')
  940. ('MAXIMUM' (PF '-' PF1) 'ABS')
  941. ('MAXIMUM' (VITF '-' VITF1) 'ABS')
  942. ('MAXIMUM' (YF '-' YF1) 'ABS') );
  943.  
  944. ERRO = 'MAXIMUM' ('PROG' ERRO
  945. ('MAXIMUM' (YF11 '-' SCAF12) 'ABS')
  946. ('MAXIMUM' (YF12 '-' SCAF11) 'ABS')
  947. ('MAXIMUM' (YF12 '-' SCAF13) 'ABS')) ;
  948.  
  949.  
  950. 'SI' (ERRO > 1.0D-12);
  951. 'ERREUR' 5;
  952. 'FINSI' ;
  953.  
  954.  
  955.  
  956. *********************************************************
  957. *** Control des etats sur la surface qui contient P1 ****
  958. *********************************************************
  959.  
  960.  
  961. GEOP1 = ('DOMA' $DOMTOT 'FACEL') 'ELEM' 'APPUYE' 'LARGEMENT' P1 ;
  962. GEOP2 = ('DOMA' $DOMTOT 'FACE') 'ELEM' 'APPUYE' 'LARGEMENT' P1 ;
  963. GEOP12 = ('DOMA' $DOMTOT 'FACEL') 'ELEM' 'APPUYE' 'LARGEMENT' P2 ;
  964.  
  965.  
  966.  
  967. ROGEOP1 = 'REDU' ROF GEOP1;
  968. VGEOP1 = 'REDU' VITF GEOP1;
  969. PGEOP1 = 'REDU' PF GEOP1;
  970. REFGEOP1 = 'REDU' VITF GEOP2;
  971. YNGEOP1 = 'REDU' YF GEOP1 ;
  972.  
  973. ROGEOP2 = 'REDU' ROF GEOP12 ;
  974. VGEOP2 = 'REDU' VITF GEOP12 ;
  975. PGEOP2 = 'REDU' PF GEOP12 ;
  976. YNGEOP2 = 'REDU' YF GEOP12 ;
  977.  
  978. rofag = 'EXTRAIRE' ROGEOP1 'SCAL' 1 1 1;
  979. rofad = 'EXTRAIRE' ROGEOP1 'SCAL' 1 1 3;
  980. unfag = 'EXTRAIRE' VGEOP1 'UN ' 1 1 1;
  981. unfad = 'EXTRAIRE' VGEOP1 'UN ' 1 1 3;
  982. utfag = 'EXTRAIRE' VGEOP1 'UT ' 1 1 1;
  983. utfad = 'EXTRAIRE' VGEOP1 'UT ' 1 1 3;
  984. pfag = 'EXTRAIRE' PGEOP1 'SCAL' 1 1 1;
  985. pfad = 'EXTRAIRE' PGEOP1 'SCAL' 1 1 3;
  986. yh2fag = 'EXTRAIRE' YNGEOP1 'H2 ' 1 1 1;
  987. yh2fad= 'EXTRAIRE' YNGEOP1 'H2 ' 1 1 3;
  988. yn2fag = 'EXTRAIRE' YNGEOP1 'N2 ' 1 1 1;
  989. yn2fad= 'EXTRAIRE' YNGEOP1 'N2 ' 1 1 3;
  990.  
  991. rofag2 = 'EXTRAIRE' ROGEOP2 'SCAL' 1 1 1;
  992. rofad2 = 'EXTRAIRE' ROGEOP2 'SCAL' 1 1 3;
  993. unfag2 = 'EXTRAIRE' VGEOP2 'UN ' 1 1 1;
  994. unfad2 = 'EXTRAIRE' VGEOP2 'UN ' 1 1 3;
  995. utfag2 = 'EXTRAIRE' VGEOP2 'UT ' 1 1 1;
  996. utfad2 = 'EXTRAIRE' VGEOP2 'UT ' 1 1 3;
  997. pfag2 = 'EXTRAIRE' PGEOP2 'SCAL' 1 1 1;
  998. pfad2 = 'EXTRAIRE' PGEOP2 'SCAL' 1 1 3;
  999. yh2fag2 = 'EXTRAIRE' YNGEOP2 'H2 ' 1 1 1;
  1000. yh2fad2 = 'EXTRAIRE' YNGEOP2 'H2 ' 1 1 3;
  1001. yn2fag2 = 'EXTRAIRE' YNGEOP2 'N2 ' 1 1 1;
  1002. yn2fad2 = 'EXTRAIRE' YNGEOP2 'N2 ' 1 1 3;
  1003.  
  1004. *
  1005. **** Orientation de la normal n de castem par raport a la
  1006. * notre; t est par consequence (sur P1)
  1007. *
  1008.  
  1009. NCOS = 'EXTRAIRE' REFGEOP1 'NX' 1 1 1;
  1010. NSIN = 'EXTRAIRE' REFGEOP1 'NY' 1 1 1;
  1011.  
  1012. 'SI' (('ABS' NCOS) > ('ABS' NSIN));
  1013. ORIENT = ('COS' ANGLE) '/' NCOS;
  1014. 'SINON';
  1015. ORIENT = ('SIN' ANGLE) '/' NSIN;
  1016. 'FINSI' ;
  1017.  
  1018. ORIENT = 'SIGNE' 'FLOTTANT' ORIENT;
  1019.  
  1020.  
  1021. *
  1022. **** ORIENT = -1 -> Mon etat gauche est son etat droite
  1023. *
  1024.  
  1025. 'SI' (ORIENT > 0);
  1026.  
  1027. ERRLIG = 'PROG' rofag (unfag '*' ORIENT)
  1028. (utfag '*' ORIENT) pfag yh2fag yn2fag ;
  1029.  
  1030. ERRLID = 'PROG' rofad (unfad '*' ORIENT)
  1031. (utfad '*' ORIENT) pfad yh2fad yn2fad ;
  1032.  
  1033. 'SINON' ;
  1034.  
  1035. ERRLID = 'PROG' rofag (unfag '*' ORIENT)
  1036. (utfag '*' ORIENT) pfag yh2fag yn2fag;
  1037.  
  1038. ERRLIG = 'PROG' rofad (unfad '*' ORIENT)
  1039. (utfad '*' ORIENT) pfad yh2fad yn2fad ;
  1040.  
  1041. 'FINSI' ;
  1042.  
  1043. ERRO = 'MAXIMUM' ('PROG'
  1044. ('MAXIMUM' (ETATG '-' ERRLIG) 'ABS')
  1045. ('MAXIMUM' (ETATD '-' ERRLID) 'ABS')
  1046. );
  1047.  
  1048. 'SI' (ERRO > 1.0D-14)
  1049. 'MESSAGE' 'Ordre en espace = 2';
  1050. 'MESSAGE' 'Ordre en temps = 2';
  1051. 'ERREUR' 5 ;
  1052. 'FINSI' ;
  1053.  
  1054. *
  1055. **** Sir le mur (P2)
  1056. *
  1057.  
  1058. ERRLIG = 'PROG' rofag2 unfag2 utfag2 pfag2 yh2fag2 yn2fag2 ;
  1059.  
  1060. ERRLID = 'PROG' rofad2 unfad2 utfad2 pfad2 yh2fad2 yn2fad2 ;
  1061.  
  1062. ERRO = 'MAXIMUM' ('PROG'
  1063. ('MAXIMUM' (ETATG2 '-' ERRLIG) 'ABS')
  1064. ('MAXIMUM' (ETATD2 '-' ERRLID) 'ABS')
  1065. );
  1066.  
  1067.  
  1068. 'SI' (ERRO > 1.0D-14)
  1069. 'Message' 'Mur' ;
  1070. 'MESSAGE' 'Ordre en espace = 2';
  1071. 'MESSAGE' 'Ordre en temps = 2';
  1072. 'ERREUR' 5 ;
  1073. 'FINSI' ;
  1074.  
  1075. 'FIN' BL1;
  1076.  
  1077.  
  1078.  
  1079. 'FIN' ;
  1080.  
  1081.  
  1082.  
  1083.  
  1084.  
  1085.  
  1086.  
  1087.  
  1088.  

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