Télécharger la plaquette Cast3M


fr en

Télécharger cl_E_7.dgibi

Retour à la liste

Numérotation des lignes : 

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

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