Télécharger la plaquette Cast3M


fr en

Télécharger cl_B_2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  cl_B_2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *C***********************************************************************
  5. *C*                                                                     *
  6. *C* PROJET      : Opérateur CLMI                                        *
  7. *C* NOM         : profil_B_2.dgibi                                      *
  8. *C* DESCRIPTION : Jeu de données pour le calcul d'un écoulement en      *
  9. *C*               fluide visqueux entre 2 plaques                       *
  10. *C* LANGAGE     : Gibiane                                               *
  11. *C* AUTEUR      : Guillaume VENCO - DRN/DMT/SEMT/LTMF                   *
  12. *C*                                                                     *
  13. *C***********************************************************************
  14. *C*                                                                     *
  15. *C* APPELES     : Opérateur CLMI                                        *
  16. *C*                                                                     *
  17. *C***********************************************************************
  18. *C*                                                                     *
  19. *C*                                                                     *
  20. *C***********************************************************************
  21. *C*                                                                     *
  22. *C* VERSION     : 15/05/2000                                            *
  23. *C* CREATION    : 18/02/2000                                            *
  24. *C*                                                                     *
  25. *C***********************************************************************
  26.  
  27. 'OPTION' dime 2 elem seg2 ;
  28. 'OPTION'  trace x ;
  29. GRAPH=FAUX;
  30. COMPLET=FAUX;
  31.  
  32. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  33. *>>>>>>>>>>>>>Maillage et domaine>>>>>>>>>>>>>>>>>>>>>>>>>>>
  34. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  35. P1 = 0.1 0.;
  36. p2 =  10. 0. ;
  37. nb1 = 10;
  38. p1p2 = 'DROIT' p1 p2 'DINI' 0.00001 'DFIN' 1. ;
  39. *p1p2 = 'DROIT' p1 p2 20;
  40. mt = p1p2 ;
  41. * 'TRACER'  mt ;
  42.  
  43. *Définition du modèle
  44. Mmt = 'CHANGER'  mt 'QUAF' ;
  45. Mail = 'CHANGER' Mmt 'POI1';
  46. Pt0 = Mail 'ELEM' 1;
  47. $mt = 'MODELISER'  Mmt 'NAVIER_STOKES' 'LINE' ;
  48.  
  49. *Données
  50. Ue0 = 0.001;
  51. NU = 1.d-6 ;
  52.  
  53. *valeurs initiales: valeurs données par Falkner-Skan
  54. D20 =0.66411*((0.1*(NU/Ue0))**(0.5));
  55. D30 = 1.57257*D20;
  56.  
  57.  
  58. *Pas de temps et nombre de pas
  59. DT = 0.1;
  60. NBIT = 30;
  61. Si Complet ;
  62. NBIT = 300;
  63. Finsi ;
  64. *Champ de vitesse constant
  65. X = 'COORDONNEE' 1 mt;
  66. UE = 'KCHT' $mt scal 'SOMMET' Ue0;
  67. DUE = 'KCHT' $mt scal 'SOMMET' 0.;
  68.  
  69. *Coefficients de l'équation
  70. A = 'KCHT' $mt scal sommet 1.;
  71. B = 'KCHT' $mt scal sommet 1.;
  72. SM  = 'KCHT'  $mt scal sommet 1. ;
  73.  
  74. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  75. *>>>>>>>>>>>>Création de la table RV et de la table RK>>>>>>
  76. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  77. *               dY         Cf               >>>>>>>>>>>>>>>
  78. *CLMI résout A ---- + B = ---- = SM         >>>>>>>>>>>>>>>>
  79. *               dX          2               >>>>>>>>>>>>>>>>
  80. * Les coef A, B et SM sont calculés dans    >>>>>>>>>>>>>>>>
  81. *        CLMI suivant les cas               >>>>>>>>>>>>>>>>
  82. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  83.  
  84. rv = 'EQEX'  'OMEGA' 1.  'NITER' 1  'ITMA' 0   'FIDT' 1
  85.      'OPTI' 'EF' 'IMPL' 'CENTREE' 'BDF2' 
  86.      'ZONE' $mt
  87.      'OPER' 'DFDT' 1. 'D2NM' 'D2N2' DT  'INCO' 'D2' ;
  88. *     'OPER' 'DFDT' 1. 'D2NM' DT  'INCO' 'D2' ;
  89.  
  90. rv = 'EQEX' rv
  91.      'OPTI' 'EF' 'IMPL' 'SUPG'
  92.      'ZONE'  $mt
  93.      'OPER'  'CLMI' 3 1 UE DUE 'D2NM' 'D3NM' 'INCO' 'D2';
  94.  
  95. rv = 'EQEX' rv
  96.      'OPTI' 'EF' 'IMPL' 'CENTREE' 'BDF2' 
  97.      'ZONE' $mt
  98. *    'OPER' 'DFDT' 1. 'D3NM' DT  'INCO' 'D3' ;
  99.      'OPER'  'DFDT' 1. 'D3NM' 'D3N2' DT  'INCO' 'D3';
  100.  
  101. rv = 'EQEX' rv
  102.      'OPTI' 'EF' 'IMPL' 'SUPG'
  103.      'ZONE'  $mt
  104.      'OPER'  'CLMI' 3 2 UE DUE 'D3NM' 'D2NM' 'INCO' 'D3';
  105.  
  106. rv = 'EQEX' rv
  107.       'CLIM' 'D2' 'TIMP' pt0 D20
  108.       'CLIM' 'D3' 'TIMP' pt0 D30;
  109.  
  110.  
  111. rv.inco = 'TABLE'  inco ;
  112. rv.inco.'D2'=  'KCHT'  $mt scal sommet D20;
  113. rv.inco.'D2NM'=  'KCHT'  $mt scal sommet D20;
  114. rv.inco.'D2N2'=  'KCHT'  $mt scal sommet D20;
  115. rv.inco.'D3'=  'KCHT'  $mt scal sommet D30;
  116. rv.inco.'D3NM'=  'KCHT'  $mt scal sommet D30;
  117. rv.inco.'D3NM'=  'KCHT'  $mt scal sommet D30;
  118. rv.inco.'D3N2'=  'KCHT'  $mt scal sommet D30;
  119.  
  120. ltps = 'PROG' 0. ;
  121. lteta = 'PROG' 0. ;
  122. TPS = 0. ;
  123. nupt = 0 ;
  124. mts1 = 'DOMA'  $mt 'SOMMET' ;
  125. pt1 = 'POIN'  mts1 2 ;
  126. 'LISTE'  pt1 ;
  127.  
  128. *'OPTION'  donn 5 ;
  129.  
  130. 'REPETER'  BCLT nbit ;
  131. nupt=nupt+1;
  132. tps=tps+dt ;
  133. mess ' Pas de temps DT=' dt ' Temps = ' tps ;
  134. ltps=ltps et (prog tps) ;
  135. EXEC  rv ;
  136. teta='EXTRAIRE'  (rv.inco.'D2') 'SCAL' pt1 ;
  137. mess ' Teta=' teta ;
  138. lteta= lteta et ('PROG' teta);
  139.  
  140.  
  141. YNM = rv.'INCO'.'D2NM';
  142. ZNM = rv.'INCO'.'D3NM';
  143. rv.'INCO'.'D2N2' = rv.'INCO'.'D2NM';
  144. rv.'INCO'.'D2NM' = rv.'INCO'.'D2';
  145. rv.'INCO'.'D3N2' = rv.'INCO'.'D3NM';
  146. rv.'INCO'.'D3NM' = rv.'INCO'.'D3';
  147.  
  148. 'FIN'  bclt ;
  149.  
  150. li4 = 'PROG' 'EXPO' 'A' -1. 'B' 0. ltps;
  151. ev2 = 'EVOL' 'ROUGE' manu ltps li4 ;
  152. evt= evol manu ltps lteta;
  153.  
  154.  
  155.  
  156. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  157. *>>>>>>>>>>>>>>>>>>>>>>> résultats>>>>>>>>>>>>>>>>>>>>>>>>>>
  158. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  159. *quantité de mouvement
  160. d21 = (rv.inco.'D2');
  161. ev = 'EVOL' 'ROUGE' 'CHPO' d21 mt;
  162. RD2 = (UE*d21)*(1./NU);
  163. ev = 'EVOL' 'ROUGE' 'CHPO' RD2 mt;
  164. *épaisseur d'énergie cinétique
  165. d31= (rv.inco.'D3');
  166. ev10 = 'EVOL' 'ROUGE' 'CHPO' d31 mt;
  167. RD3 = UE*d31/NU;
  168. ev10 = 'EVOL' 'ROUGE' 'CHPO' RD3 mt;
  169. *Epaisseur de déplacement
  170. d11 = (rv.inco.'D1');
  171. RD1 = UE*d11/NU;
  172. evd1 = 'EVOL' 'ROUGE' 'CHPO' RD1 mt;
  173. *évolution de la vitesse extérieure et du gradient
  174. vit = 'EVOL' 'CHPO' UE mt;
  175. gvit = 'EVOL' 'ROUG' 'CHPO' DUE mt;
  176. *évolution du frottement à la paroi
  177. Cf=2*(0.22052D-6)*((YNM*UE)**(-1));
  178. f = 'EVOL' 'ROUGE' 'CHPO' Cf mt;
  179. *évolution du facteur de forme H final
  180. *hev = 'EVOL' 'ROUGE''CHPO' H mt;
  181. *évolution du facteur H32
  182. H32=d31*(d21**(-1));
  183. ev32 = 'EVOL' 'ROUGE' 'CHPO' H32 mt;
  184.  
  185. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  186. *solution théorique de FALKNER SKAN>>>>>>>>>>>>>>>>>>>>>>>>
  187. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  188. *caractéristique des solutions de FALKNER-SKAN
  189. cteta=0.66411;
  190. cd1 = 1.72079;
  191. ccf=0.33206;
  192. Hfs=2.5911;
  193. *coefficient de frottement
  194. RX = X*UE*(1./NU);
  195. Cfb = 2*ccf*(RX**(-0.5));
  196. ev3 = 'EVOL' 'BLEU' 'CHPO' Cfb mt;
  197. *quantité de mvt
  198. * Rteta = 0.664(Rx)^0.5
  199. RD2b = (cteta*(RX**0.5)); 
  200. ev4 =  'EVOL' 'BLEU' 'CHPO' RD2b mt;
  201. *facteur de forme
  202. Hb = 'KCHT' $mt scal sommet Hfs;
  203. ev5 = 'EVOL' 'BLEU' 'CHPO' Hb mt;
  204. *coef H32 théorique
  205. H32b='KCHT' $mt scal sommet 1.57257;
  206. ev11= 'EVOL' 'BLEU' 'CHPO' H32b mt;
  207. *épaisseur de'énergie cinétique
  208. RD3t=H32b*RD2b;
  209. ev12='EVOL' 'BLEU' 'CHPO' RD3t mt;
  210. *épaisseur de déplacement théorique
  211. RD1b = (cd1*(RX**0.5));
  212. evdb =  'EVOL' 'BLEU' 'CHPO' RD1b mt;
  213. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  214. *Courbes>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  215. *>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  216. *'OPTION' trace psc;
  217.  
  218. Si GRAPH ;
  219. 'DESSIN'  (evt) 'TITR' 'Convergence de la valeur au dernier noeud'
  220.            'TITY' 'Valeurs de Y' 'TITX' 'Temps';
  221.  
  222. 'DESSIN' (vit 'ET' gvit)
  223.      'TITRE' 'Vitesse et gradient de vitesse extérieure';
  224.  
  225. TAB1 = 'TABLE' ;
  226. TAB1.'TITRE' = 'TABLE' ;
  227. TAB1.'TITRE' . 1 = MOT ' Simulation';
  228. TAB1.'TITRE' . 2 = MOT ' Falkner-Skan';
  229.  
  230. 'DESSIN' (f 'ET' ev3) 'TITRE' 'Coefficient de frottement à la paroi'
  231.            'TITX' 'Abscisses'
  232.            'TITY' 'Cf'
  233.            'GRIL'
  234.            'LEGE' TAB1;     
  235.  
  236. 'DESSIN' (ev 'ET' ev4) 'TITRE' 'Epaisseur de quantité de mouvement'
  237.        'TITY' 'Rteta'
  238.        'TITX' 'X'
  239.        'GRIL'
  240.        'LEGE' TAB1;
  241.  
  242. 'DESSIN' (ev10 'ET' ev12) 'TITRE' 'Epaisseur d énergie cinétique'
  243.        'TITY' 'Rdelta3'
  244.        'TITX' 'X'
  245.        'GRIL'
  246.        'LEGE' TAB1;
  247.  
  248. 'DESSIN' (evd1 'ET' evdb) 'TITRE' 'Epaisseur de déplacement'
  249.        'TITY' 'Rdelta1'
  250.        'TITX' 'X'
  251.        'GRIL'
  252.        'LEGE' TAB1;         
  253.  
  254. *'DESSIN' (hev 'ET' ev5) 'TITRE' 'Evolution du favteur de forme'
  255. *        'TITX' 'Abscisses'
  256. *        'TITY' 'Cf'
  257. *        'LEGE' TAB1;
  258.  
  259. 'DESSIN' (ev32 'ET' ev11) 'TITRE' 'Evolution du facteur H32'
  260.         'TITX' 'Abscisses'
  261.         'TITY' 'H32'
  262.         'LEGE' TAB1;
  263.  
  264. Finsi ;
  265. *calcul de l'erreur au sens de la norme L2
  266. diag = 'DOMA' $mt 'XXDIAGSI';
  267. r = (RD2-RD2B)**2;
  268. erl2 = 'SOMT'(diag*r);
  269. 'MESSAGE' 'Erreur au sens de la norme L2=';
  270. 'LISTE'  erl2;
  271. **************************
  272. *autres verif ************
  273. **************************
  274. w=RD2*(Rx**(-0.5));
  275. wev = 'EVOL' 'ROUGE' 'CHPO' mt w;
  276. wtheo = 'KCHT'$mt 'SCAL' 'SOMMET' cteta;
  277. wtev = 'EVOL' 'BLEU' 'CHPO' mt wtheo;
  278.  
  279. Si GRAPH ;
  280. 'DESSIN' (wev 'ET' wtev) 'TITRE' 'comparaison theo/simul coef 1'
  281.         'TITX' 'Abscisses'
  282.         'TITY' 'Rteta/(Rx**0.5)'
  283.         'LEGE' TAB1;
  284. Finsi ;
  285.  
  286. *q=(CfKP/2.)*(Rx**(0.5));
  287. *qev = 'EVOL' 'ROUGE' 'CHPO' mt q;
  288. *qtheo = 'KCHT'$mt 'SCAL' 'SOMMET' ccf;
  289. *qtev = 'EVOL' 'BLEU' 'CHPO' mt qtheo;
  290. *'DESSIN' (qev 'ET' qtev) 'TITRE' 'comparaison theo/simul coef 2'
  291. *        'TITX' 'Abscisses'
  292. *        'TITY' '0.5*Cf*(Rx**0.5)'
  293. *        'LEGE' TAB1;         
  294.  'FIN';
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  
  303.  
  304.  
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  

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