Télécharger cl_D_7.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : cl_D_7.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *C***********************************************************************
  5. *C* *
  6. *C* PROJET : Opérateur CLMI *
  7. *C* NOM : cl_D_6 .dgibi *
  8. *C* DESCRIPTION : Jeu de données pour le calcul de la couche limite *
  9. *C* turbulente en utilisant les relations de fermeture de *
  10. *C* Cousteix - Cas de la plaque plane avec gradient de *
  11. *C* pression *
  12. *C* LANGAGE : Gibiane *
  13. *C* AUTEUR : Guillaume VENCO - DRN/DMT/SEMT/LTMF *
  14. *C* *
  15. *C***********************************************************************
  16. *C* *
  17. *C* APPELES : Opérateur CLMI *
  18. *C* *
  19. *C***********************************************************************
  20. *C* *
  21. *C* *
  22. *C***********************************************************************
  23. *C* *
  24. *C* VERSION : 05/05/2000 *
  25. *C* CREATION : 05/04/2000 *
  26. *C* *
  27. *C***********************************************************************
  28.  
  29. 'OPTION' dime 2 elem seg2 ;
  30. 'OPTION' trace x ;
  31. GRAPH=FAUX;
  32. COMPLET=FAUX;
  33.  
  34. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  35. *>>>>>>>>>>>>>Maillage et domaine>>>>>>>>>>>>>>>>>>>>>>>>>>>
  36. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  37. P1 = 0.15 0.;
  38. p2 = 0.8 0. ;
  39. p3 = 10. 0;
  40. nb1 = 10;
  41. *p1p2 = 'DROIT' p1 p2 'DINI' 0.00001 'DFIN' 0.01;
  42. p1p2 = 'DROIT' p1 p2 500;
  43. *p2p3 = 'DROIT' p2 p3 'DINI' 0.01 'DFIN' 1;
  44. *mt = p1p2 'ET' p2p3;
  45. mt = p1p2;
  46. *'ELIMINATION' 1E-3 mt;
  47. * 'TRACER' mt ;
  48.  
  49. *Définition du modèle
  50. Mmt = 'CHANGER' mt 'QUAF' ;
  51. Mail = 'CHANGER' Mmt 'POI1';
  52. Pt0 = Mail 'ELEM' 1;
  53. $mt = 'MODELISER' Mmt 'NAVIER_STOKES' 'LINE' ;
  54.  
  55. *Données
  56. Ue0 = 1.;
  57. NU = 1.d-6 ;
  58. *Valeurs initiales données par la loi de puissance
  59. D20=4.335d-4;
  60. Q10=7.85*D20;
  61. *Pas de temps et nombre de pas
  62. DT = 1.d-1;
  63.  
  64. NBIT = 10;
  65. Si Complet;
  66. NBIT = 100;
  67. Finsi;
  68.  
  69. DT2= 1.d-1;
  70. NBI2=1;
  71. *Champ de vitesse à l'extérieur
  72.  
  73. * vitesse exponentielle
  74. m1= 'DOMA' $mt 'SOMMET';
  75. X = 'COORDONNEE' 1 m1;
  76. *X2 = (-5.)*(X-0.15);
  77. *X1 = 'EXP' (X2);
  78. UEX = 2-(X-0.15);
  79. UE = 'KCHT' $mt scal 'SOMMET' UEX;
  80. *DUEX = (-5.)*X1;
  81. DUEX = -1.;
  82. DUE = 'KCHT' $mt scal 'SOMMET' DUEX;
  83.  
  84. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  85. *>>>>>>>>>>>>Création de la table RV et de la table RK>>>>>>
  86. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  87.  
  88. *Dans CLMI le 1er coef correspond au type de fermeture
  89. *le 2ème coef correspond au n°de l'équation à traiter
  90. *1=QDM 2=Energie cinétique
  91.  
  92. rv = 'EQEX' 'OMEGA' 1. 'NITER' 1 'ITMA' 0 'FIDT' 1
  93. 'OPTI' 'EF' 'IMPL' 'CENTREE' 'BDF2'
  94. 'ZONE' $mt
  95. 'OPER' 'DFDT' 1. 'D2NM' 'D2N2' 'DT' 'INCO' 'D2' ;
  96. * 'OPER' 'DFDT' 1. 'D2NM' DT 'INCO' 'D2' ;
  97.  
  98. rv = 'EQEX' rv
  99. 'OPTI' 'EF' 'IMPL' 'SUPG'
  100. 'ZONE' $mt
  101. 'OPER' 'CLMI' 4 1 UE DUE 'D2NM' 'Q1NM'
  102. 'INCO' 'D2' ;
  103.  
  104. rv = 'EQEX' rv
  105. 'OPTI' 'EF' 'IMPL' 'CENTREE' 'BDF2'
  106. 'ZONE' $mt
  107. 'OPER' 'DFDT' 1. 'Q1NM' 'Q1N2' 'DT' 'INCO' 'Q1' ;
  108. * 'OPER' 'DFDT' 1. 'D1NM' DT 'INCO' 'D1' ;
  109.  
  110. rv = 'EQEX' rv
  111. 'OPTI' 'EF' 'IMPL' 'SUPG'
  112. 'ZONE' $mt
  113. 'OPER' 'CLMI' 4 3 UE DUE 'Q1NM' 'D2NM'
  114. 'INCO' 'Q1' ;
  115.  
  116. rv = 'EQEX' rv
  117. 'CLIM' 'D2' 'TIMP' pt0 D20
  118. 'CLIM' 'Q1' 'TIMP' pt0 Q10;
  119.  
  120.  
  121. rv.inco = 'TABLE' inco ;
  122. rv.inco.'D2'= 'KCHT' $mt scal sommet D20;
  123. rv.inco.'D2NM'= 'KCHT' $mt scal sommet D20;
  124. rv.inco.'D2N2'= 'KCHT' $mt scal sommet D20;
  125. rv.inco.'Q1'= 'KCHT' $mt scal sommet Q10;
  126. rv.inco.'Q1NM'= 'KCHT' $mt scal sommet Q10;
  127. rv.inco.'Q1N2'= 'KCHT' $mt scal sommet Q10;
  128. rv.inco.'DT' = dt;
  129. ltps = 'PROG' 0. ;
  130. lteta = 'PROG' 0. ;
  131. TPS = 0. ;
  132. nupt = 0 ;
  133. mts1 = 'DOMA' $mt 'SOMMET' ;
  134. pt1 = 'POIN' mts1 2 ;
  135. 'LISTE' pt1 ;
  136.  
  137. *'OPTION' donn 5 ;
  138.  
  139. 'REPETER' BCLT nbit ;
  140. nupt=nupt+1;
  141. tps=tps+dt ;
  142. mess ' Pas de temps DT=' dt ' Temps = ' tps ;
  143. ltps=ltps et (prog tps) ;
  144. EXEC rv ;
  145. teta='EXTRAIRE' (rv.inco.'D2') 'SCAL' pt1 ;
  146. mess ' Teta=' teta ;
  147. lteta= lteta et ('PROG' teta);
  148.  
  149.  
  150. YNM = rv.'INCO'.'D2NM';
  151. ZNM = rv.'INCO'.'Q1NM';
  152.  
  153. rv.'INCO'.'D2N2' = rv.'INCO'.'D2NM';
  154. rv.'INCO'.'D2NM' = rv.'INCO'.'D2';
  155. rv.'INCO'.'Q1N2' = rv.'INCO'.'Q1NM';
  156. rv.'INCO'.'Q1NM' = rv.'INCO'.'Q1';
  157.  
  158. 'FIN' bclt ;
  159.  
  160.  
  161. *******************************************************
  162.  
  163. 'REPETER' BCLT2 nbi2 ;
  164. rv.inco.'DT'=1.d-2;
  165. DT=1.d-2;
  166. nupt=nupt+1;
  167. tps=tps+dt2 ;
  168. mess ' Pas de temps DT=' dt ' Temps = ' tps ;
  169. ltps=ltps et (prog tps) ;
  170. EXEC rv ;
  171. teta='EXTRAIRE' (rv.inco.'D2') 'SCAL' pt1 ;
  172. mess ' Teta=' teta ;
  173. lteta= lteta et ('PROG' teta);
  174.  
  175.  
  176. YNM = rv.'INCO'.'D2NM';
  177. ZNM = rv.'INCO'.'Q1NM';
  178.  
  179. rv.'INCO'.'D2N2' = rv.'INCO'.'D2NM';
  180. rv.'INCO'.'D2NM' = rv.'INCO'.'D2';
  181. rv.'INCO'.'Q1N2' = rv.'INCO'.'Q1NM';
  182. rv.'INCO'.'Q1NM' = rv.'INCO'.'Q1';
  183.  
  184. 'FIN' bclt2 ;
  185.  
  186. li4 = 'PROG' 'EXPO' 'A' -1. 'B' 0. ltps;
  187. evt2 = 'EVOL' 'ROUGE' manu ltps li4 ;
  188. evt= evol manu ltps lteta;
  189.  
  190.  
  191.  
  192. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  193. *>>>>>>>>>>>>>>>>>>>>>>> résultats>>>>>>>>>>>>>>>>>>>>>>>>>>
  194. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  195. *évolution de la vitesse extérieure et du gradient
  196. vit = 'EVOL' 'CHPO' UE mt;
  197. gvit = 'EVOL' 'ROUG' 'CHPO' DUE mt;
  198. *quantité de mouvement
  199. D2S = (rv.inco.'D2');
  200. RD2S = UE*D2S/NU;
  201. ev1s = 'EVOL' 'ROUGE' 'CHPO' RD2S mt;
  202. *différence entre (d-d1)
  203. Q1S = (rv.inco.'Q1');
  204. RQ1S = UE*Q1S/NU;
  205. ev2s = 'EVOL' 'ROUGE' 'CHPO' RQ1S mt;
  206. *Epaisseur de déplacement
  207. D1S = (rv.'INCO'.'D1');
  208. RD1S = UE*D1S/NU;
  209. ev3s = 'EVOL' 'ROUGE' 'CHPO' RD1S mt;
  210. *évolution du frottement à la paroi
  211. CfS = rv.'INCO'.'CF';
  212. ev4s = 'EVOL' 'ROUGE' 'CHPO' CfS mt;
  213. *évolution du facteur de forme H final
  214. HS = rv.'INCO'.'H';
  215. ev5s = 'EVOL' 'ROUGE''CHPO' HS mt;
  216. *évolution du facteur H*
  217. HHS = Q1S*(D2S**(-1));
  218. ev6s = 'EVOL' 'ROUGE' 'CHPO' HHS mt;
  219.  
  220. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  221. *solution théorique cf:SCHLICHTING >>>>>>>>>>>>>>>>>>>>>>>>
  222. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  223. *RX = X*UE*(1./NU);
  224. RX = 2.*X*(1./NU);
  225. *Quantité de mvt
  226. *----Rteta donné par Nikuradze
  227. RD2N = (0.01738/1.304)*(RX**(0.861));
  228. ev1n = 'EVOL' 'BLEU' 'CHPO' RD2N mt;
  229. *----Rteta par la loi de puissance
  230. RD2P = (1./1.e-6)*(0.036*X*(RX**(-1./5)))*UE;
  231. ev1p = 'EVOL' 'VERT' 'CHPO' RD2P mt;
  232. *différence (d-d1) théorique
  233. Q1P=(0.37-(1.2857*0.036))*(X*(RX**(-1./5)));
  234. RQ1P=UE*Q1P/NU;
  235. ev2p='EVOL' 'VERT' 'CHPO' RQ1P mt;
  236. *Epaisseur de déplacement
  237. *----Rdelta1 par la Nikuradze
  238. RD1N = 1.304*RD2N;
  239. ev3n = 'EVOL' 'BLEU' 'CHPO' RD1N mt;
  240. *----Rdelta1 par la loi de puissance
  241. RD1P = (9./7.)*RD2P;
  242. ev3p = 'EVOL' 'VERT' 'CHPO' RD1P mt;
  243. *Coefficient de frottement
  244. *----Formule par Nikuradze
  245. CfN = 0.02296*(RX**(-0.139));
  246. ev4n = 'EVOL' 'BLEU' 'CHPO' CfN mt;
  247. *----Formule par loi de puissance
  248. CfP = 0.0576*(RX**(-1./5.));
  249. ev4p = 'EVOL' 'VERT' 'CHPO' CfP mt;
  250. *facteur de forme
  251. *----H par Nikuradze
  252. HN = 'KCHT' $mt scal sommet 1.304;
  253. ev5n = 'EVOL' 'BLEU' 'CHPO' HN mt;
  254. *----H par loi de puissance
  255. HP = 'KCHT' $mt scal sommet 1.2857;
  256. ev5p = 'EVOL' 'VERT' 'CHPO' HP mt;
  257. *coef H* théorique
  258. HHP='KCHT' $mt scal sommet 7.89;
  259. ev6p= 'EVOL' 'BLEU' 'CHPO' HHP mt;
  260.  
  261. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  262. *Courbes>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  263. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  264. Si GRAPH ;
  265.  
  266. 'DESSIN' (evt) 'TITR' 'Convergence de la valeur au dernier noeud'
  267. 'TITY' 'Valeurs de Y' 'TITX' 'Temps';
  268.  
  269. TAB0 = 'TABLE' ;
  270. TAB0.'TITRE' = 'TABLE' ;
  271. TAB0.'TITRE' . 1 = MOT ' Vitesse';
  272. TAB0.'TITRE' . 2 = MOT ' Gradient';
  273.  
  274. 'DESSIN' (vit 'ET' gvit)
  275. 'TITRE' 'Vitesse et gradient de vitesse extérieure';
  276.  
  277. TAB1 = 'TABLE' ;
  278. TAB1.'TITRE' = 'TABLE' ;
  279. TAB1.'TITRE' . 1 = MOT ' Simulation';
  280. TAB1.'TITRE' . 2 = MOT ' Nikuradze';
  281. TAB1.'TITRE' . 3 = MOT ' Loi de puissance';
  282.  
  283. TAB2 = 'TABLE' ;
  284. TAB2.'TITRE' = 'TABLE' ;
  285. TAB2.'TITRE' . 1 = MOT ' Simulation';
  286. TAB2.'TITRE' . 2 = MOT ' Loi de puissance';
  287.  
  288. 'DESSIN' (ev1s 'ET' ev1n 'ET' ev1p)
  289. 'TITRE' 'Epaisseur de quantité de mouvement'
  290. 'TITY' 'Rteta'
  291. 'TITX' 'Abscisses'
  292. 'GRIL'
  293. 'LEGE' TAB1;
  294.  
  295. 'DESSIN' (ev2s 'ET' ev2p) 'TITRE' 'différence entre delta et delta1'
  296. 'TITY' 'Rdelta-Rdelta1'
  297. 'TITX' 'Abscisses'
  298. 'GRIL'
  299. 'LEGE' TAB2;
  300.  
  301. 'DESSIN' (ev3s 'ET' ev3n 'ET' ev3p)
  302. 'TITRE' 'Epaisseur de déplacement'
  303. 'TITY' 'Rdelta1'
  304. 'TITX' 'Abscisses'
  305. 'GRIL'
  306. 'LEGE' TAB1;
  307.  
  308. 'DESSIN' (ev4s 'ET' ev4n 'ET' ev4p)
  309. 'TITRE' 'Coefficient de frottement à la paroi'
  310. 'TITX' 'Abscisses'
  311. 'TITY' 'Cf'
  312. 'GRIL'
  313. 'LEGE' TAB1;
  314.  
  315. 'DESSIN' (ev5s 'ET' ev5n 'ET' ev5p)
  316. 'TITRE' 'Evolution du facteur de forme'
  317. 'TITX' 'Abscisses'
  318. 'TITY' 'H=delta1/theta'
  319. 'LEGE' TAB1;
  320.  
  321. 'DESSIN' (ev6s 'ET' ev6p) 'TITRE' 'Evolution du facteur H*'
  322. 'TITX' 'Abscisses'
  323. 'TITY' 'H'
  324. 'LEGE' TAB1;
  325. Finsi;
  326.  
  327. *calcul de l'erreur au sens de la norme L2
  328. diag = 'DOMA' $mt 'XXDIAGSI';
  329. r = (RD2S-RD2N)**2;
  330. r2= (RD2S-RD2P)**2;
  331. EL2N = 'SOMT'(diag*r);
  332. EL2P = 'SOMT'(diag*r2);
  333. 'MESSAGE' 'Erreur sur teta au sens de la norme L2=';
  334. 'LISTE' EL2N;
  335. 'LISTE' EL2P;
  336. r3=(RQ1S-RQ1P)**2;
  337. EL3P= 'SOMT'(diag*r3);
  338. 'MESSAGE' 'Erreur sur (delta-delta1) au sens de la norme L2=';
  339. 'LISTE' EL3P;
  340. *autre dessin
  341. a = 'EXTRAIRE' ev1s 'SCAL';
  342. log = 'LOG' a;
  343. l10= 'LOG' 10;
  344. loga = (1./l10)*log;
  345. b = 'EXTRAIRE' ev5s 'SCAL';
  346. evH = 'EVOL' 'MANU' 'Rtheta' a 'H' b;
  347. SI Graph;
  348. 'DESSIN' evH
  349. 'LOGX'
  350. 'TITRE' 'facteur de forme en couche limite turbulente'
  351. 'TITX' 'Log(Rtheta)'
  352. 'TITY' 'H'
  353. 'GRIL';
  354. Finsi;
  355. 'FIN';
  356.  
  357.  
  358.  
  359.  
  360.  
  361.  
  362.  
  363.  
  364.  
  365.  
  366.  
  367.  
  368.  
  369.  
  370.  
  371.  
  372.  

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