Télécharger cl_D_4.dgibi

Retour à la liste

Numérotation des lignes :

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

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