Télécharger pente3D.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : pente3D.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***********************************************************
  5. **** Finite Volume, "Cell-Centred Formulation". ****
  6. **** PENT, operator to compute gradients and limiters ****
  7. **** ****
  8. **** 3D VALIDATION ****
  9. **** ****
  10. **** T. KLOCZKO, LTMF JULY 2005 ****
  11. **** ****
  12. ***********************************************************
  13. *
  14. *** Choix du type d'élément
  15. *
  16. TYEL = TET4 ;
  17. * TYEL = PYR5 ;
  18. * TYEL = PRI6 ;
  19. * TYEL = CUB8 ;
  20. *
  21. *** Options gibiane
  22. *
  23. 'OPTION' 'DIME' 3 'ELEM' TYEL 'TRAC' 'X'
  24. 'ISOV' 'SULI' 'ECHO' 1 ;
  25. *
  26. *** Option de tracé
  27. *
  28. * GRAPH = VRAI ;
  29. GRAPH = FAUX ;
  30. *
  31. ***************************************************
  32. *** CANAL SECTION CARRÉE AVEC BOSSE SINUSOÏDALE ***
  33. ***************************************************
  34. *
  35. RAF = 1 ;
  36. VRAF = 9. '*' RAF ;
  37. Lref = 1. ;
  38. Size = Lref '/' VRAF ;
  39.  
  40. Lx = 4. ;
  41. Ly = 0.5 ;
  42. Lz = 1. ;
  43. DX = Size ;
  44. DY = Size ;
  45. DZ = Size ;
  46. NX = 'ENTIER' (Lx '/' Size '+' 0.1) ;
  47. NY = 'ENTIER' (Ly '/' Size '+' 0.1) ;
  48. NZ = 'ENTIER' (Lz '/' Size '+' 0.1) ;
  49. *
  50. *** Sommets du domaine de calcul
  51. *
  52. P0 = 0.0 0.0 0.0 ;
  53. P1 = (Lx '/' 4.) 0.0 0.0 ;
  54. P2 = (3. '*' Lx '/' 4.) 0.0 0.0 ;
  55. P3 = Lx 0.0 0.0 ;
  56. P4 = Lx 0.0 Lz ;
  57. P5 = 0.0 0.0 Lz ;
  58.  
  59. P6 = 0.0 Ly 0.0 ;
  60. P7 = (Lx '/' 4.) Ly 0.0 ;
  61. P8 = (3. '*' Lx '/' 4.) Ly 0.0 ;
  62. P9 = Lx Ly 0.0 ;
  63. P10 = Lx Ly Lz ;
  64. P11 = 0.0 Ly Lz ;
  65.  
  66. *
  67. *** Calcul des arêtes
  68. *
  69. P0P1 = 'DROIT' P0 P1 DINI DX DFIN DX ;
  70. P1P2 = 'DROIT' P1 P2 DINI DX DFIN DX ;
  71. P2P3 = 'DROIT' P2 P3 DINI DX DFIN DX ;
  72. P0P3 = P0P1 'ET' P1P2 'ET' P2P3 ;
  73. P4P5 = 'DROIT' P4 P5 DINI DX DFIN DX ;
  74.  
  75. P6P7 = 'DROIT' P6 P7 DINI DX DFIN DX ;
  76. P7P8 = 'DROIT' P7 P8 DINI DX DFIN DX ;
  77. P8P9 = 'DROIT' P8 P9 DINI DX DFIN DX ;
  78. P6P9 = P6P7 'ET' P7P8 'ET' P8P9 ;
  79. P10P11 = 'DROIT' P10 P11 DINI DX DFIN DX ;
  80.  
  81. P0P6 = 'DROIT' P0 P6 DINI DY DFIN DY ;
  82. P5P11 = 'DROIT' P5 P11 DINI DY DFIN DY ;
  83. P3P9 = 'DROIT' P3 P9 DINI DY DFIN DY ;
  84. P4P10 = 'DROIT' P4 P10 DINI DY DFIN DY ;
  85.  
  86. P0P5 = 'DROIT' P0 P5 DINI DZ DFIN DZ ;
  87. P3P4 = 'DROIT' P3 P4 DINI DZ DFIN DZ ;
  88. P9P10 = 'DROIT' P9 P10 DINI DZ DFIN DZ ;
  89. P6P11 = 'DROIT' P6 P11 DINI DZ DFIN DZ ;
  90. *
  91. *** Calcul des surfaces
  92. *
  93. * Inlet
  94. SURFI = 'SURFACE' (('INVERSE' P0P6) 'ET' P0P5 'ET' P5P11
  95. 'ET' ('INVERSE' P6P11)) 'PLAN' ;
  96. SURFI = 'COULEUR' SURFI 'VERT' ;
  97.  
  98. * Outlet
  99. SURFO = 'SURFACE' (('INVERSE' P3P9) 'ET' P3P4 'ET' P4P10
  100. 'ET' ('INVERSE' P9P10)) 'PLAN' ;
  101. SURFO = 'COULEUR' SURFO 'ROUG' ;
  102.  
  103. * Up
  104. SURFU = 'SURFACE' (('INVERSE' P4P5) 'ET' P4P10 'ET' P10P11
  105. 'ET' ('INVERSE' P5P11)) 'PLAN' ;
  106. SURFU = 'COULEUR' SURFU 'TURQ' ;
  107.  
  108. * Down
  109. *
  110. * On trace en premier lieu une surface plane que l'on déforme
  111. * par un champ de déplacement ad hoc.
  112. SURFD = 'SURFACE' (P0P3 'ET' P3P9 'ET' ('INVERSE'P6P9)
  113. 'ET' ('INVERSE' P0P6)) 'PLAN' ;
  114. SURFD = 'COULEUR' SURFD 'BLEU' ;
  115. *
  116. *** Création du déplacement pour générer la bosse
  117. *
  118. NBPTD = 'NBNO' SURFD ;
  119. SURD1 = 'CHANGER' SURFD 'POI1' ;
  120. CHPMOUV CACA = 'KOPS' MATRIK ;
  121. 'REPETER' DBCL NBPTD ;
  122. PSD = 'POINT' SURD1 &DBCL ;
  123. XSD = ('COORDONNEE' 1 PSD) ;
  124. 'SI' (XSD < (Lx '/' 4.)) ;
  125. CHPI = 'MANU' 'CHPO' PSD 3
  126. 'UX' 0.D0 'UY' 0.D0 'UZ' 0.D0 ;
  127. 'SINON' ;
  128. 'SI' (XSD > (3. '*' Lx '/' 4.)) ;
  129. CHPI = 'MANU' 'CHPO' PSD 3
  130. 'UX' 0.D0 'UY' 0.D0 'UZ' 0.D0 ;
  131. 'SINON' ;
  132. ZSD = (0.1D0 '*' Lx '/' 4.) '*' (1.D0 +
  133. ('COS' (180.D0 '*' ((4. '*' XSD '/' Lx) '-' 2.)))) ;
  134. CHPI = 'MANU' 'CHPO' PSD 3
  135. 'UX' 0.D0 'UY' 0.D0 'UZ' ZSD ;
  136. 'FINSI' ;
  137. 'FINSI' ;
  138. CHPMOUV = CHPMOUV '+' CHPI ;
  139. 'FIN' DBCL ;
  140. *
  141. * Modification de la forme de la surface du bas
  142. FSURD = 'FORM' CHPMOUV ;
  143.  
  144. * Front
  145. SURFF = 'SURFACE' (P0P3 'ET' P3P4 'ET' P4P5
  146. 'ET' ('INVERSE' P0P5)) 'PLAN' ;
  147. * Back
  148. SURFB = 'SURFACE' (P6P9 'ET' P9P10 'ET' P10P11
  149. 'ET' ('INVERSE' P6P11)) 'PLAN' ;
  150. * Total Face
  151. SURTOT = (SURFF 'ET' SURFB 'ET' SURFI
  152. 'ET' SURFO 'ET' SURFU 'ET' SURFD) ;
  153. 'ELIMINATION' SURTOT (Size '/' 1.D3) ;
  154.  
  155. 'SI' GRAPH ;
  156. 'TRACER' SURTOT 'TITRE'
  157. ('CHAINE' 'Enveloppe, nbel=' ('NBEL' SURTOT)) ;
  158. 'FINSI' ;
  159. *
  160. *** Volume
  161. *
  162. DOMTOT = 'VOLUME' SURTOT ;
  163. *
  164. *** Surface latérale
  165. *
  166. SURTOT = 'ENVELOPPE' DOMTOT ;
  167. si graph;'TRACER' surtot ;finsi;
  168.  
  169. MAICON = DOMTOT 'ELEM' 'APPUYE' 'LARGEMENT' SURTOT ;
  170. DOMINT = DOMTOT 'DIFF' MAICON 'COULEUR' 'ROUGE' ;
  171. *
  172. **** Creation of MODELS
  173. *
  174.  
  175. MDOMTOT = 'MODELISER' DOMTOT 'EULER' ;
  176. MDOMINT = 'MODELISER' DOMINT 'EULER' ;
  177. MSURTOT = 'MODELISER' SURTOT 'EULER' ;
  178.  
  179. TDOMTOT = 'DOMA' MDOMTOT 'VF' ;
  180. TDOMINT = 'DOMA' MDOMINT 'VF' ;
  181. TSURTOT = 'DOMA' MSURTOT 'VF' ;
  182.  
  183. QDOMTOT = TDOMTOT . 'QUAF' ;
  184. QDOMINT = TDOMINT . 'QUAF' ;
  185. QSURTOT = TSURTOT . 'QUAF' ;
  186.  
  187. 'ELIMINATION' (QDOMTOT 'ET' QDOMINT 'ET' QSURTOT)
  188. (1.d-3 '/' (RAF '*' 1.)) ;
  189.  
  190. 'SI' GRAPH ;
  191. 'TRACER' DOMTOT 'TITRE'
  192. ('CHAINE' 'Domaine, nbel=' ('NBEL' DOMTOT)) ;
  193. 'TRACER' MAICON 'TITRE'
  194. ('CHAINE' 'Mailles extérieures, nbel=' ('NBEL' MAICON)) ;
  195. 'TRACER' DOMINT 'TITRE'
  196. ('CHAINE' 'Mailles intérieures, nbel=' ('NBEL' DOMINT)) ;
  197. 'FINSI' ;
  198. *
  199. **** The interface points (for limit condition)
  200. *
  201. SURC = ('DOMA' MSURTOT 'CENTRE') 'COULEUR' 'BLEU';
  202. 'REPETER' BLLIM (('NBEL' SURC)) ;
  203. POIN0 = 'POIN' SURC &BLLIM ;
  204. XFAC YFAC ZFAC = 'COORDONNEE' POIN0 ;
  205. PFAC = ('DOMA' MDOMTOT 'FACE') 'POIN' 'PROC' (XFAC YFAC ZFAC);
  206. 'SI' (&BLLIM 'EGA' 1) ;
  207. GEOLIM = 'MANUEL' 'POI1' PFAC 'BLEU' ;
  208. 'SINON' ;
  209. GEOLIM = GEOLIM 'ET' ('MANUEL' 'POI1' PFAC 'BLEU') ;
  210. 'FINSI' ;
  211. 'FIN' BLLIM ;
  212. *
  213. 'SI' GRAPH ;
  214. 'TRACER' (DOMTOT 'ET' GEOLIM)
  215. 'TITRE' 'DOMTOT et points faces aux bords';
  216. 'FINSI' ;
  217. *
  218. ***********************************************
  219. ********* CHP: scalar linear field *********
  220. ********* EULESCAL *********
  221. ***********************************************
  222. *
  223. * (A11 '*' XX) '+' (A12 '*' YY) '+' (A13 '*' ZZ) '+' A0;
  224. *
  225.  
  226. A11 = 12.01517 ;
  227. A12 = 13.1421 ;
  228. A13 = 14.0412 ;
  229. A0 = -3.21 ;
  230.  
  231. B11 = 22.01517 ;
  232. B12 = 43.1421 ;
  233. B13 = 35.5463 ;
  234. B0 = -5.21 ;
  235.  
  236. C11 = 24.01517 ;
  237. C12 = 46.1421 ;
  238. C13 = 39.5463 ;
  239. C0 = -6.21 ;
  240.  
  241. XX YY ZZ = 'COORDONNEE' ('DOMA' MDOMTOT 'CENTRE') ;
  242.  
  243. CHP1 = (A11 '*' XX) '+' (A12 '*' YY) '+' (A13 '*' ZZ) '+' A0 ;
  244.  
  245. CHP2 = (B11 '*' XX) '+' (B12 '*' YY) '+' (B13 '*' ZZ) '+' B0 ;
  246.  
  247. CHP3 = (C11 '*' XX) '+' (C12 '*' YY) '+' (C13 '*' ZZ) '+' C0 ;
  248.  
  249. CHP1 = 'NOMC' 'C1' CHP1 'NATU' 'DISCRET' ;
  250. CHP2 = 'NOMC' 'C2' CHP2 'NATU' 'DISCRET' ;
  251. CHP3 = 'NOMC' 'C3' CHP3 'NATU' 'DISCRET' ;
  252.  
  253. CHP = CHP1 'ET' CHP2 'ET' CHP3 ;
  254.  
  255. MOTC = ('MOTS' 'C1' 'C2' 'C3') ;
  256.  
  257. GRCHP LIMCH GRG = 'PENT' MDOMTOT 'CENTRE'
  258. 'EULESCAL' 'NOLIMITE' MOTC CHP ;
  259.  
  260. LIMCH1 = 'MANU' 'CHPO' ('DOMA' MDOMTOT 'CENTRE')
  261. 3 'P1' 1.0 'P2' 1.0 'P3' 1.0 ;
  262.  
  263. ERRO = 'MAXIMUM' (LIMCH1 '-' LIMCH) 'ABS' ;
  264.  
  265. 'SI' (ERRO > 1.0D-8) ;
  266. 'MESSAGE' ;
  267. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  268. 'MESSAGE' ;
  269. 'ERREUR' 5;
  270. 'FINSI' ;
  271.  
  272. *
  273. * TEST1 : we check that we have the same values if we use
  274. * GRG to compute gradients
  275. *
  276.  
  277. MOTC = ('MOTS' 'C1' 'C2' 'C3') ;
  278. GRCHP1 LIMCH = 'PENT' MDOMTOT 'CENTRE'
  279. 'EULESCAL' 'LIMITEUR' MOTC
  280. CHP 'GRADGEO' GRG ;
  281. 'LISTE' MOTC ;
  282.  
  283. ERRO = 'MAXIMUM' (GRCHP '-' GRCHP1) 'ABS' ;
  284.  
  285. 'SI' (ERRO > 1.0D-8) ;
  286. 'MESSAGE' ;
  287. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  288. 'MESSAGE' ;
  289. 'ERREUR' 5;
  290. 'FINSI' ;
  291.  
  292.  
  293. *
  294. * TEST2 : we check that the computation is exact on the internal domain
  295. *
  296.  
  297. GRCHPI = 'REDU' GRCHP ('DOMA' MDOMINT 'CENTRE') ;
  298.  
  299.  
  300. GR1DX = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P1DX' A11
  301. 'NATU' 'DISCRET' ;
  302.  
  303. GR1DY = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P1DY' A12
  304. 'NATU' 'DISCRET';
  305.  
  306. GR1DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P1DZ' A13
  307. 'NATU' 'DISCRET';
  308.  
  309. GR2DX = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P2DX' B11
  310. 'NATU' 'DISCRET';
  311.  
  312. GR2DY = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P2DY' B12
  313. 'NATU' 'DISCRET';
  314.  
  315. GR2DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P2DZ' B13
  316. 'NATU' 'DISCRET';
  317.  
  318. GR3DX = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P3DX' C11
  319. 'NATU' 'DISCRET';
  320.  
  321. GR3DY = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P3DY' C12
  322. 'NATU' 'DISCRET';
  323.  
  324. GR3DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P3DZ' C13
  325. 'NATU' 'DISCRET';
  326.  
  327. GREXA = GR1DX 'ET' GR1DY 'ET' GR1DZ 'ET'
  328. GR2DX 'ET' GR2DY 'ET' GR2DZ 'ET'
  329. GR3DX 'ET' GR3DY 'ET' GR3DZ ;
  330.  
  331. ERRCH = GREXA '-' GRCHPI ;
  332.  
  333. LMOT = 'MOTS' 'P1DX' 'P1DY' 'P1DZ'
  334. 'P2DX' 'P2DY' 'P2DZ'
  335. 'P3DX' 'P3DY' 'P3DZ' ;
  336. MERRCH = ('PSCAL' ERRCH ERRCH LMOT LMOT) '**' 0.5 ;
  337.  
  338.  
  339. *
  340. *** Error graphics (we change the error CHAMPOINT into a CHAMELEM
  341. *
  342.  
  343. 'SI' GRAPH ;
  344. MOD1 = 'MODELISER' ('DOMA' MDOMINT 'MAILLAGE') 'THERMIQUE' ;
  345. CHM_MER = 'KCHA' MDOMINT 'CHAM' MERRCH ;
  346. 'TRACER' CHM_MER MOD1 'TITR'
  347. ('CHAINE' 'EULESCAL: erreur') ;
  348. 'FINSI' ;
  349.  
  350. ERRO = 'MAXIMUM' ERRCH 'ABS' ;
  351.  
  352. 'SI' (ERRO > 1.0D-8) ;
  353. 'MESSAGE' ;
  354. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  355. 'MESSAGE' ;
  356. 'ERREUR' 5;
  357. 'FINSI' ;
  358.  
  359.  
  360. *
  361. * TEST3 : we impose the linear field values on the border.
  362. * In that case error must be "zero" everywhere.
  363. *
  364.  
  365. XLIM YLIM ZLIM = 'COORDONNEE' GEOLIM ;
  366.  
  367. CHPLI1 = (A11 '*' XLIM) '+' (A12 '*' YLIM) '+' (A13 '*' ZLIM) '+' A0;
  368.  
  369. CHPLI2 = (B11 '*' XLIM) '+' (B12 '*' YLIM) '+' (B13 '*' ZLIM) '+' B0;
  370.  
  371. CHPLI3 = (C11 '*' XLIM) '+' (C12 '*' YLIM) '+' (C13 '*' ZLIM) '+' C0;
  372.  
  373. CHPLI1 = 'NOMC' 'C1' CHPLI1 'NATU' 'DISCRET' ;
  374. CHPLI2 = 'NOMC' 'C2' CHPLI2 'NATU' 'DISCRET' ;
  375. CHPLI3 = 'NOMC' 'C3' CHPLI3 'NATU' 'DISCRET' ;
  376.  
  377. CHPLI = CHPLI1 'ET' CHPLI2 'ET' CHPLI3 ;
  378.  
  379. MOTC = ('MOTS' 'C1' 'C2' 'C3') ;
  380.  
  381. GRCHPL LIMCH GRGL = 'PENT' MDOMTOT 'CENTRE'
  382. 'EULESCAL' 'NOLIMITE' MOTC CHP 'CLIM' CHPLI ;
  383.  
  384. MOTC = ('MOTS' 'C1' 'C2' 'C3') ;
  385. GRCHP1 LIMCH = 'PENT' MDOMTOT 'CENTRE'
  386. 'EULESCAL' 'NOLIMITE' MOTC
  387. CHP 'CLIM' CHPLI 'GRADGEO' GRGL ;
  388.  
  389. ERRO = 'MAXIMUM' (GRCHPL '-' GRCHP1) 'ABS' ;
  390.  
  391. 'SI' (ERRO > 1.0D-8) ;
  392. 'MESSAGE' ;
  393. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  394. 'MESSAGE' ;
  395. 'ERREUR' 5;
  396. 'FINSI' ;
  397.  
  398. GR1DX = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P1DX' A11
  399. 'NATU' 'DISCRET' ;
  400.  
  401. GR1DY = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P1DY' A12
  402. 'NATU' 'DISCRET';
  403.  
  404. GR1DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P1DZ' A13
  405. 'NATU' 'DISCRET';
  406.  
  407. GR2DX = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P2DX' B11
  408. 'NATU' 'DISCRET';
  409.  
  410. GR2DY = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P2DY' B12
  411. 'NATU' 'DISCRET';
  412.  
  413. GR2DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P2DZ' B13
  414. 'NATU' 'DISCRET';
  415.  
  416. GR3DX = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P3DX' C11
  417. 'NATU' 'DISCRET';
  418.  
  419. GR3DY = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P3DY' C12
  420. 'NATU' 'DISCRET';
  421.  
  422. GR3DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P3DZ' C13
  423. 'NATU' 'DISCRET';
  424.  
  425. GREXA = GR1DX 'ET' GR1DY 'ET' GR1DZ 'ET'
  426. GR2DX 'ET' GR2DY 'ET' GR2DZ 'ET'
  427. GR3DX 'ET' GR3DY 'ET' GR3DZ ;
  428.  
  429. ERRCH = GREXA '-' GRCHPL ;
  430.  
  431. LMOT = 'MOTS' 'P1DX' 'P1DY' 'P1DZ'
  432. 'P2DX' 'P2DY' 'P2DZ'
  433. 'P3DX' 'P3DY' 'P3DZ' ;
  434.  
  435. MERRCH = ('PSCAL' ERRCH ERRCH LMOT LMOT) '**' 0.5 ;
  436.  
  437. 'SI' GRAPH ;
  438. MOD1 = 'MODELISER' ('DOMA' MDOMTOT 'MAILLAGE') 'THERMIQUE' ;
  439. CHM_MER = 'KCHA' MDOMTOT 'CHAM' MERRCH ;
  440. 'TRACER' CHM_MER MOD1 'TITR'
  441. ('CHAINE' 'EULESCAL avec C.L. : Erreur') ;
  442. 'FINSI' ;
  443.  
  444. ERRO = 'MAXIMUM' ERRCH 'ABS' ;
  445.  
  446. 'SI' (ERRO > 1.0D-8) ;
  447. 'MESSAGE' ;
  448. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  449. 'MESSAGE' ;
  450. 'ERREUR' 5;
  451. 'FINSI' ;
  452.  
  453.  
  454. *
  455. **** TEST4: boudary condition for a scalar field (EULESCAL)
  456. *
  457.  
  458. *** Symétrie de référence et objet modèle
  459. *
  460. VOLI = DOMTOT 'ELEM' 'APPUYE' 'LARGEMENT' SURFI ;
  461. $VOLI = 'MODELISER' VOLI 'EULER' ;
  462. TVOLI = 'DOMA' $VOLI 'VF' ;
  463. QVOLI = TVOLI . 'QUAF' ;
  464. *
  465. *** Champoints quelconques
  466. *
  467. XX YY ZZ = 'COORDONNEE' ('DOMA' $VOLI 'CENTRE') ;
  468.  
  469. CHP1 = ('NOMC' 'C1' (yy * ('EXP'(XX) '+' ZZ)) 'NATU' 'DISCRET') ;
  470. CHP2 = ('NOMC' 'C2' (yy '+' (3 '*' ZZ)) 'NATU' 'DISCRET') ;
  471. CHP3 = ('NOMC' 'C3' (XX * ('EXP'(ZZ) '+' ZZ)) 'NATU' 'DISCRET') ;
  472. *
  473. *** Plan de symétrie
  474. *
  475. PP0 = P0 'PLUS' ((-1.d-11) 0.0 0.0 ) ;
  476. PP5 = P5 'PLUS' ((-1.d-11) 0.0 0.0 ) ;
  477. PP6 = P6 'PLUS' ((-1.d-11) 0.0 0.0 ) ;
  478. *
  479. *** Champoints symétriques
  480. *
  481. CNEW1 CHP11 = ('EXTRAIRE' CHP1 'MAILLAGE') CHP1 'SYMETRIE'
  482. 'PLAN' PP0 PP5 PP6 ;
  483. CNEW2 CHP21 = ('EXTRAIRE' CHP2 'MAILLAGE') CHP2 'SYMETRIE'
  484. 'PLAN' PP0 PP5 PP6 ;
  485. CNEW3 CHP31 = ('EXTRAIRE' CHP3 'MAILLAGE') CHP3 'SYMETRIE'
  486. 'PLAN' PP0 PP5 PP6 ;
  487. 'ELIMINATION' (CNEW1 'ET' CNEW2 'ET' CNEW3) 1.0D-9 ;
  488. *
  489. *** Visualisation des maillages
  490. *
  491. VSYM = (VOLI 'SYMETRIE' 'PLAN' PP0 PP5 PP6) 'COULEUR' 'ROUG' ;
  492. VTOT = (VOLI 'ET' VSYM) ;
  493. 'ELIMINATION' VTOT 1.d-5 ;
  494. *
  495. 'SI' GRAPH ;
  496. 'TRACER' (VSYM 'ET' VTOT) 'TITRE' 'Traitement de C.L.' ;
  497. 'FINSI' ;
  498. *
  499. *** Objet modèle
  500. *
  501. $VTOT = 'MODELISER' VTOT 'EULER' ;
  502. TVTOT = 'DOMA' $VTOT 'VF' ;
  503. QVTOT = TVTOT . 'QUAF' ;
  504. 'ELIMINATION' QVOLI 1.D-5 QVTOT ;
  505. 'ELIMINATION' QVTOT 1.D-5 CNEW1 ;
  506.  
  507. *
  508. *** Calcul du gradient
  509. *
  510.  
  511. CHPPET = CHP1 'ET' CHP2 'ET' CHP3 ;
  512. CHPSYM = CHP11 'ET' CHP21 'ET' CHP31 ;
  513. CHPTOT = CHPSYM 'ET' CHPPET ;
  514.  
  515.  
  516. GRCHP LIMCHP GRG = 'PENT' $VTOT 'CENTRE' 'EULESCAL' 'NOLIMITE'
  517. ('MOTS' 'C1' 'C2' 'C3') CHPTOT ;
  518.  
  519. GRCHP1 LIMCHP GRG1 = 'PENT' $VOLI 'CENTRE' 'EULESCAL' 'NOLIMITE'
  520. ('MOTS' 'C1' 'C2' 'C3')
  521. ('REDU' CHPTOT ('DOMA' $VOLI 'CENTRE')) ;
  522.  
  523.  
  524. ERRCH = GRCHP1 '-' ('REDU' GRCHP ('DOMA' $VOLI 'CENTRE')) ;
  525.  
  526. LMOT = 'MOTS' 'P1DX' 'P1DY' 'P1DZ'
  527. 'P2DX' 'P2DY' 'P2DZ'
  528. 'P3DX' 'P3DY' 'P3DZ' ;
  529.  
  530. MERRCH = ('PSCAL' ERRCH ERRCH LMOT LMOT) '**' 0.5 ;
  531.  
  532. 'SI' GRAPH ;
  533. MOD1 = 'MODELISER' ('DOMA' $VOLI 'MAILLAGE') 'THERMIQUE' ;
  534. CHM_MER = 'KCHA' $VOLI 'CHAM' MERRCH ;
  535. 'TRACER' CHM_MER MOD1
  536. 'TITR' ('CHAINE' 'EULESCAL sur les bords: erreur') ;
  537. 'FINSI' ;
  538.  
  539. ERRO = 'MAXIMUM' MERRCH 'ABS' ;
  540. 'LISTE' ERRO ;
  541.  
  542.  
  543. 'SI' (ERRO > 1.0D-8) ;
  544. 'MESSAGE' ;
  545. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  546. 'MESSAGE' ;
  547. 'ERREUR' 5 ;
  548. 'FINSI' ;
  549. *
  550. ***************************************************
  551. ********* CHP: vectorial linear field ************
  552. ********* EULEVECT ****************
  553. ***************************************************
  554. *
  555. * (A11 '*' XX) '+' (A12 '*' YY) '+' (A13 '*' ZZ) '+' A0;
  556. *
  557.  
  558. A11 = 11.3 ;
  559. A12 = 0.23 ;
  560. A13 = 2.46 ;
  561. A0 = 1.0 ;
  562.  
  563. B11 = 77.0 ;
  564. B12 = 0.6 ;
  565. B13 = 9.2 ;
  566. B0 = 9.11 ;
  567.  
  568. C11 = 1.89 ;
  569. C12 = 44.7 ;
  570. C13 = 0.87 ;
  571. C0 = 4.0 ;
  572.  
  573. XX YY ZZ = 'COORDONNEE' ('DOMA' MDOMTOT 'CENTRE') ;
  574.  
  575. CHP1 = (A11 '*' XX) '+' (A12 '*' YY) '+' (A13 '*' ZZ) '+' A0 ;
  576.  
  577. CHP2 = (B11 '*' XX) '+' (B12 '*' YY) '+' (B13 '*' ZZ) '+' B0 ;
  578.  
  579. CHP3 = (C11 '*' XX) '+' (C12 '*' YY) '+' (C13 '*' ZZ) '+' C0 ;
  580.  
  581. CHP1 = 'NOMC' 'UX' CHP1 'NATU' 'DISCRET' ;
  582. CHP2 = 'NOMC' 'UY' CHP2 'NATU' 'DISCRET' ;
  583. CHP3 = 'NOMC' 'UZ' CHP3 'NATU' 'DISCRET' ;
  584.  
  585. CHP = CHP1 'ET' CHP2 'ET' CHP3 ;
  586.  
  587. MOTC = ('MOTS' 'UX' 'UY' 'UZ') ;
  588. GRCHP LIMCHP GRG = 'PENT' MDOMTOT 'CENTRE'
  589. 'EULEVECT' 'NOLIMITE' MOTC CHP ;
  590. 'LISTE' MOTC ;
  591.  
  592. LIMCH1 = 'MANU' 'CHPO' ('DOMA' MDOMTOT 'CENTRE')
  593. 3 'P1' 1.0 'P2' 1.0 'P3' 1.0 ;
  594.  
  595. ERRO = 'MAXIMUM' (LIMCH1 '-' LIMCHP) 'ABS' ;
  596.  
  597. 'SI' (ERRO > 1.0D-8) ;
  598. 'MESSAGE' ;
  599. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  600. 'MESSAGE' ;
  601. 'ERREUR' 5;
  602. 'FINSI' ;
  603.  
  604. *
  605. * TEST1 : we check that we have the same values if we use
  606. * GRG to compute gradients
  607. *
  608.  
  609. GRCHP1 LIMCHP = 'PENT' MDOMTOT 'CENTRE'
  610. 'EULEVECT' 'NOLIMITE' MOTC
  611. CHP 'GRADGEO' GRG ;
  612.  
  613. ERRO = 'MAXIMUM' (GRCHP '-' GRCHP1) 'ABS' ;
  614.  
  615. 'SI' (ERRO > 5.0D-8) ;
  616. 'MESSAGE' ;
  617. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  618. 'MESSAGE' ;
  619. 'ERREUR' 5;
  620. 'FINSI' ;
  621.  
  622. *
  623. * TEST2 : we check that the computation is exact on the internal domain
  624. *
  625.  
  626. GRCHPI = 'REDU' GRCHP ('DOMA' MDOMINT 'CENTRE') ;
  627.  
  628.  
  629. GR1DX = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P1DX' A11
  630. 'NATU' 'DISCRET' ;
  631.  
  632. GR1DY = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P1DY' A12
  633. 'NATU' 'DISCRET';
  634.  
  635. GR1DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P1DZ' A13
  636. 'NATU' 'DISCRET';
  637.  
  638. GR2DX = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P2DX' B11
  639. 'NATU' 'DISCRET';
  640.  
  641. GR2DY = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P2DY' B12
  642. 'NATU' 'DISCRET';
  643.  
  644. GR2DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P2DZ' B13
  645. 'NATU' 'DISCRET';
  646.  
  647. GR3DX = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P3DX' C11
  648. 'NATU' 'DISCRET';
  649.  
  650. GR3DY = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P3DY' C12
  651. 'NATU' 'DISCRET';
  652.  
  653. GR3DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMINT 'CENTRE') 1 'P3DZ' C13
  654. 'NATU' 'DISCRET';
  655.  
  656. GREXA = GR1DX 'ET' GRqDY 'sp`n style="color: #0000BA;">ET' GR1DZ 'ET'
  657. GR2DX 'ET' GR2DY 'ET' GR2DZ 'ET'
  658. GR3DX 'ET' GR3DY 'ET' GR3DZ ;
  659.  
  660. ERRCH = GREXA '-' GRCHPI ;
  661.  
  662. LMOT = 'MOTS' 'P1DX' 'P1DY' 'P1DZ'
  663. 'P2DX' 'P2DY' 'P2DZ'
  664. 'P3DX' 'P3DY' 'P3DZ' ;
  665.  
  666. MERRCH = ('PSCAL' ERRCH ERRCH LMOT LMOT) '**' 0.5 ;
  667.  
  668. *
  669. *** CREATION DE 'MODEL' POUR GRAPHIQUER LE CHAMELEM ***
  670. *
  671.  
  672. 'SI' GRAPH ;
  673. MOD1 = 'MODELISER' ('DOMA' MDOMINT 'MAILLAGE') 'THERMIQUE' ;
  674. CHM_MER = 'KCHA' MDOMINT 'CHAM' MERRCH ;
  675. 'TRACER' CHM_MER MOD1 'TITR'
  676. ('CHAINE' 'EULEVECT: Erreur') ;
  677. 'FINSI' ;
  678.  
  679. ERRO = 'MAXIMUM' ERRCH 'ABS' ;
  680.  
  681. 'SI' (ERRO > 5.0D-8) ;
  682. 'MESSAGE' ;
  683. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  684. 'MESSAGE' ;
  685. 'ERREUR' 5;
  686. 'FINSI' ;
  687.  
  688.  
  689. *
  690. * TEST3 : we impose the linear field values on the border.
  691. * In that case error must be "zero" everywhere.
  692. *
  693.  
  694. XLIM YLIM ZLIM = 'COORDONNEE' GEOLIM ;
  695.  
  696. CHPLI1 = (A11 '*' XLIM) '+' (A12 '*' YLIM) '+' (A13 '*' ZLIM) '+' A0 ;
  697.  
  698. CHPLI2 = (B11 '*' XLIM) '+' (B12 '*' YLIM) '+' (B13 '*' ZLIM) '+' B0 ;
  699.  
  700. CHPLI3 = (C11 '*' XLIM) '+' (C12 '*' YLIM) '+' (C13 '*' ZLIM) '+' C0 ;
  701.  
  702. CHPLI1 = 'NOMC' 'UX' CHPLI1 'NATU' 'DISCRET' ;
  703. CHPLI2 = 'NOMC' 'UY' CHPLI2 'NATU' 'DISCRET' ;
  704. CHPLI3 = 'NOMC' 'UZ' CHPLI3 'NATU' 'DISCRET' ;
  705.  
  706. CHPLI = CHPLI1 'ET' CHPLI2 'ET' CHPLI3 ;
  707.  
  708.  
  709. GRCHPL LIMCHP GRGL = 'PENT' MDOMTOT 'CENTRE'
  710. 'EULEVECT' 'NOLIMITE' ('MOTS' 'UX' 'UY' 'UZ')
  711. CHP 'CLIM' CHPLI ;
  712.  
  713. GRCHP1 LIMCHP = 'PENT' MDOMTOT 'CENTRE'
  714. 'EULEVECT' 'NOLIMITE' ('MOTS' 'UX' 'UY' 'UZ')
  715. CHP 'CLIM' CHPLI 'GRADGEO' GRGL ;
  716.  
  717. ERRO = 'MAXIMUM' (GRCHPL '-' GRCHP1) 'ABS' ;
  718.  
  719. 'SI' (ERRO > 1.0D-8) ;
  720. 'MESSAGE' ;
  721. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  722. 'MESSAGE' ;
  723. 'ERREUR' 5;
  724. 'FINSI' ;
  725.  
  726. GR1DX = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P1DX' A11
  727. 'NATU' 'DISCRET' ;
  728.  
  729. GR1DY = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P1DY' A12
  730. 'NATU' 'DISCRET';
  731.  
  732. GR1DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P1DZ' A13
  733. 'NATU' 'DISCRET';
  734.  
  735. GR2DX = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P2DX' B11
  736. 'NATU' 'DISCRET';
  737.  
  738. GR2DY = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P2DY' B12
  739. 'NATU' 'DISCRET';
  740.  
  741. GR2DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P2DZ' B13
  742. 'NATU' 'DISCRET';
  743.  
  744. GR3DX = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P3DX' C11
  745. 'NATU' 'DISCRET';
  746.  
  747. GR3DY = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P3DY' C12
  748. 'NATU' 'DISCRET';
  749.  
  750. GR3DZ = 'MANUEL' 'CHPO' ('DOMA' MDOMTOT 'CENTRE') 1 'P3DZ' C13
  751. 'NATU' 'DISCRET';
  752.  
  753.  
  754. GREXA = GR1DX 'ET' GR1DY 'ET' GR1DZ 'ET'
  755. GR2DX 'ET' GR2DY 'ET' GR2DZ 'ET'
  756. GR3DX 'ET' GR3DY 'ET' GR3DZ ;
  757.  
  758. ERRCH = GREXA '-' GRCHPL ;
  759.  
  760. LMOT = 'MOTS' 'P1DX' 'P1DY' 'P1DZ'
  761. 'P2DX' 'P2DY' 'P2DZ'
  762. 'P3DX' 'P3DY' 'P3DZ' ;
  763. MERRCH = ('PSCAL' ERRCH ERRCH LMOT LMOT) '**' 0.5 ;
  764.  
  765. 'SI' GRAPH ;
  766. MOD1 = 'MODELISER' ('DOMA' MDOMTOT 'MAILLAGE') 'THERMIQUE' ;
  767. CHM_MER = 'KCHA' MDOMTOT 'CHAM' MERRCH ;
  768. 'TRACER' CHM_MER MOD1 'TITR'
  769. ('CHAINE' 'EULEVECT avec C.L.: erreur') ;
  770. 'FINSI' ;
  771.  
  772. ERRO = 'MAXIMUM' ERRCH 'ABS' ;
  773.  
  774. 'SI' (ERRO > 1.0D-8) ;
  775. 'MESSAGE' ;
  776. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  777. 'MESSAGE' ;
  778. 'ERREUR' 5;
  779. 'FINSI' ;
  780.  
  781. *
  782. **** TEST4: boudary condition for a vectorial field (EULEVECT)
  783. *
  784.  
  785. *** Géométrie de base et objet modèle
  786. *
  787. VOLI = DOMTOT 'ELEM' 'APPUYE' 'LARGEMENT' SURFI ;
  788. $VOLI = 'MODELISER' VOLI 'EULER' ;
  789. TVOLI = 'DOMA' $VOLI 'VF' ;
  790. QVOLI = TVOLI . 'QUAF' ;
  791. *
  792. *** Construction d'un champoint
  793. *
  794. XX YY ZZ = 'COORDONNEE' ('DOMA' $VOLI 'CENTRE') ;
  795.  
  796. CHP1 = ('NOMC' 'UX' (yy * ('EXP'(XX) '+' ZZ)) 'NATU' 'DISCRET') ;
  797. CHP2 = ('NOMC' 'UY' (yy '+' (3 '*' ZZ)) 'NATU' 'DISCRET') ;
  798. CHP3 = ('NOMC' 'UZ' (XX * ('EXP'(ZZ) '+' ZZ)) 'NATU' 'DISCRET') ;
  799. *
  800. *** Plan de symétrie
  801. *
  802. PP0 = P0 'PLUS' ((-1.d-11) 0.0 0.0 ) ;
  803. PP5 = P5 'PLUS' ((-1.d-11) 0.0 0.0 ) ;
  804. PP6 = P6 'PLUS' ((-1.d-11) 0.0 0.0 ) ;
  805. *
  806. *** Construction du champoint symétrique par rapport au plan
  807. *
  808. CNEW1 CHP11 = ('EXTRAIRE' CHP1 'MAILLAGE') CHP1 'SYMETRIE'
  809. 'PLAN' PP0 PP5 PP6 ;
  810. CHP11 = CHP11 '*' (-1.D0) ;
  811.  
  812. CNEW2 CHP21 = ('EXTRAIRE' CHP2 'MAILLAGE') CHP2 'SYMETRIE'
  813. 'PLAN' PP0 PP5 PP6 ;
  814.  
  815. CNEW3 CHP31 = ('EXTRAIRE' CHP3 'MAILLAGE') CHP3 'SYMETRIE'
  816. 'PLAN' PP0 PP5 PP6 ;
  817. 'ELIMINATION' (CNEW1 'ET' CNEW2 'ET' CNEW3) 1.0D-9 ;
  818. *
  819. *** Visualisation des maillages
  820. *
  821. VSYM = (VOLI 'SYMETRIE' 'PLAN' PP0 PP5 PP6) 'COULEUR' 'ROUG' ;
  822. VTOT = (VOLI 'ET' VSYM) ;
  823. 'ELIMINATION' VTOT 1.d-5 ;
  824. *
  825. 'SI' GRAPH ;
  826. 'TRACER' (VSYM 'ET' VTOT) 'TITRE' 'Traitement de C.L.' ;
  827. 'FINSI' ;
  828. *
  829. *** Objet modèle
  830. *
  831. $VTOT = 'MODELISER' VTOT 'EULER' ;
  832. TVTOT = 'DOMA' $VTOT 'VF' ;
  833. QVTOT = TVTOT . 'QUAF' ;
  834. 'ELIMINATION' QVOLI 1.D-5 QVTOT ;
  835. 'ELIMINATION' QVTOT 1.D-5 CNEW1 ;
  836. *
  837. *** Calcul des gradients
  838. *
  839.  
  840. CHPPET = CHP1 'ET' CHP2 'ET' CHP3 ;
  841. CHPSYM = CHP11 'ET' CHP21 'ET' CHP31 ;
  842. CHPTOT = CHPSYM 'ET' CHPPET ;
  843.  
  844.  
  845. GRCHP LIMCHP GRG = 'PENT' $VTOT 'CENTRE' 'EULEVECT' 'NOLIMITE'
  846. ('MOTS' 'UX' 'UY' 'UZ') CHPTOT ;
  847.  
  848. GRCHP1 LIMCHP GRG1 = 'PENT' $VOLI 'CENTRE' 'EULEVECT' 'NOLIMITE'
  849. ('MOTS' 'UX' 'UY' 'UZ')
  850. ('REDU' CHPTOT ('DOMA' $VOLI 'CENTRE')) ;
  851.  
  852.  
  853. ERRCH = GRCHP1 '-' ('REDU' GRCHP ('DOMA' $VOLI 'CENTRE'));
  854.  
  855. LMOT = 'MOTS' 'P1DX' 'P1DY' 'P1DZ'
  856. 'P2DX' 'P2DY' 'P2DZ'
  857. 'P3DX' 'P3DY' 'P3DZ' ;
  858.  
  859. MERRCH = ('PSCAL' ERRCH ERRCH LMOT LMOT) '**' 0.5 ;
  860.  
  861. 'SI' GRAPH ;
  862. MOD1 = 'MODELISER' ('DOMA' $VOLI 'MAILLAGE') 'THERMIQUE' ;
  863. CHM_MER = 'KCHA' $VOLI 'CHAM' MERRCH ;
  864. 'TRACER' CHM_MER MOD1
  865. 'TITR' ('CHAINE' 'EULEVECT sur les bords: erreur') ;
  866. 'FINSI' ;
  867.  
  868. ERRO = 'MAXIMUM' MERRCH 'ABS' ;
  869. 'LISTE' ERRO ;
  870.  
  871.  
  872. 'SI' (ERRO > 1.0D-8) ;
  873. 'MESSAGE' ;
  874. 'MESSAGE' ('CHAINE' 'Probleme dans PENT');
  875. 'MESSAGE' ;
  876. 'ERREUR' 5 ;
  877. 'FINSI' ;
  878.  
  879. *
  880. *** Fin des tests
  881. *
  882. 'MESSAGE' ;
  883. 'MESSAGE' ('CHAINE' 'OK DOKY') ;
  884. 'MESSAGE' ;
  885.  
  886. 'FIN' ;
  887.  
  888.  
  889.  
  890.  
  891.  
  892.  
  893.  
  894.  
  895.  
  896.  
  897.  
  898.  
  899.  
  900.  
  901.  
  902.  
  903.  
  904.  
  905.  

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