Télécharger la plaquette Cast3M


fr en

Télécharger Dynasp.dgibi

Retour à la liste

Numérotation des lignes : 

  1. *****************************************************
  2. ************************************************************************
  3. ************************************************************************
  4. * fichier :  dynasp.dgibi                           *
  5. ** modifie le 15/06/2014 passage EQPR -> EQEX       *
  6. ** Il a été nécessaire de modifier aussi FROT       *
  7. ** soi yfrot.eso et yfrti.eso                       *
  8. *****************************************************
  9. * fichier :  dynasp.dgibi
  10. * Test DYNASP K-epsilon
  11.  
  12.   GRAPH    = 'N'                                                    ;
  13.   complet  = FAUX                                                   ;
  14.  
  15.   DISCR = 'LINE'  ;
  16.   KPRESS= 'CENTRE' ;
  17.   BETAP = 1.       ;
  18.  
  19.   titolo_g = chaine 'Dynasp-35 '                                    ;
  20.   titolo   = chaine titolo_g                                        ;
  21.   titr       titolo                                                 ;
  22.  
  23.   opti dime 2 elem qua4                                             ;
  24.  
  25. si ( complet )                                                      ;
  26.   n1   =          4 ; n2   =       4*n1 ; n3   =     -10*n1         ;
  27. sinon                                                               ;
  28.   n1   =          2 ; n2   =       4*n1 ; n3   =      -5*n1         ;
  29. finsi                                                               ;
  30.  
  31.   n4   =      -1*n1 ; n5   =      -1*n2                             ;
  32.   dmed = 2.000/n3   ; dmin = dmed/1.5   ; dmax = dmed*1.5           ;
  33.  
  34.   pa0 = 0.000 0.000 ; pb0 = 0.0435 0.000 ; pc0 = 0.200 0.000        ;
  35.   pa1 = 0.000 1.000 ; pb1 = 0.0435 1.000 ; pc1 = 0.200 1.000        ;
  36.   pa2 = 0.000 2.000 ; pb2 = 0.0435 2.000 ; pc2 = 0.200 2.000        ;
  37.   pa3 = 0.000 0.500 ; pb3 = 0.0435 0.500 ; pc3 = 0.200 0.500        ;
  38.   pa4 = 0.000 1.500 ; pb4 = 0.0435 1.500 ; pc4 = 0.200 1.500        ;
  39.  
  40.   axei  = pa0 droi n3 pa1 dini dmin dfin dmax                       ;
  41.   verti = pb0 droi n3 pb1 dini dmin dfin dmax                       ;
  42.   sidei = pc0 droi n3 pc1 dini dmin dfin dmax                       ;
  43.   axes  = pa1 droi n3 pa2 dini dmax dfin dmin                       ;
  44.   verts = pb1 droi n3 pb2 dini dmax dfin dmin                       ;
  45.   sides = pc1 droi n3 pc2 dini dmax dfin dmin                       ;
  46.  
  47.   weti  = axei  regl n1 verti                                       ;
  48.   dryi  = verti regl n2 sidei                                       ;
  49.   wets  = axes  regl n1 verts                                       ;
  50.   drys  = verts regl n2 sides                                       ;
  51.  
  52.   bot   = (cote 4 weti) et (cote 4 dryi)                            ;
  53.   hori  = (cote 4 wets) et (cote 4 drys)                            ;
  54.   in    = (cote 2 wets)                                             ;
  55.   dtop  = (cote 2 drys)                                             ;
  56.   top   = in            et  dtop                                    ;
  57.  
  58.   wet   = weti  et wets                                             ;
  59.   dry   = dryi  et drys                                             ;
  60.   axe   = axei  et axes                                             ;
  61.   side  = sidei et sides                                            ;
  62.   mt    = dry  et wet                                               ;
  63.   elim    mt  0.00001                                               ;
  64.   mt    = orie mt                                                   ;
  65.   cnt   = cont  mt                                                  ;
  66.   wall  = bot et side et dtop                                       ;
  67.   ppa3  = mt poin proc pa3                                          ;
  68.   ppb3  = mt poin proc pb3                                          ;
  69.   ppc3  = mt poin proc pc3                                          ;
  70.   ppa4  = mt poin proc pa4                                          ;
  71.   ppb4  = mt poin proc pb4                                          ;
  72.   ppc4  = mt poin proc pc4                                          ;
  73.   sez3  = ppa3 droi n1 ppb3 droi n2 ppc3                            ;
  74.   sez4  = ppa4 droi n1 ppb4 droi n2 ppc4                            ;
  75.   elim (mt et sez3 et sez4) 0.00001                                 ;
  76.   tass          mt                                                  ;
  77.  
  78.  
  79.   C        =        3.50e-2                                         ;
  80.   Vi       =        0.16                                            ;
  81.   Dp       =        4.70e-4                                         ;
  82.   S        =        0.400*0.087                                     ;
  83.   C_s      =        C/S                                             ;
  84.  
  85. * Options
  86.  
  88.  
  89. * Reading mesh
  90.  
  91.   yp            =    0.010                                          ;
  92.   rapp          =   (0.400 - yp)/0.400                              ;
  93.   depl mt affi rapp (0.0 0.0) (1.0 0.0)                             ;
  94. * trac mt                                                           ;
  95.  
  96. * Dominions
  97.  
  98.   Mmt    = 'CHAN' 'QUAF' mt                                         ;
  99.   Mbot   = 'CHAN' 'QUAF' bot                                        ;
  100.   Maxe   = 'CHAN' 'QUAF' axe                                        ;
  101.   Mside  = 'CHAN' 'QUAF' side                                       ;
  102.   Mtop   = 'CHAN' 'QUAF' top                                        ;
  103.   Min    = 'CHAN' 'QUAF' in                                         ;
  104.   Mcnt   = 'CHAN' 'QUAF' cnt                                        ;
  105.   Mwet   = 'CHAN' 'QUAF' wet                                        ;
  106.   Mdry   = 'CHAN' 'QUAF' dry                                        ;
  107.   Mwall  = 'CHAN' 'QUAF' wall                                       ;
  108.   Mdtop  = 'CHAN' 'QUAF' dtop                                       ;
  109.  
  110. 'ELIM' (Mmt et Mbot et Maxe et Mside et Mtop et Min et Mcnt et Mwet
  111.  et Mdry et Mwall et Mdtop) 1.e-5                                   ;
  112.  
  113.   $mt    = 'MODE' Mmt  'NAVIER_STOKES' DISCR                        ;
  114.   'DOMA' $mt 'IMPR'                                                 ;
  115.   $bot   = 'MODE' Mbot   'NAVIER_STOKES' DISCR                      ;
  116.   $axe   = 'MODE' Maxe   'NAVIER_STOKES' DISCR                      ;
  117.   $side  = 'MODE' Mside  'NAVIER_STOKES' DISCR                      ;
  118.   $top   = 'MODE' Mtop   'NAVIER_STOKES' DISCR                      ;
  119.   $in    = 'MODE' Min    'NAVIER_STOKES' DISCR                      ;
  120.   $cnt   = 'MODE' Mcnt   'NAVIER_STOKES' DISCR                      ;
  121.   $wet   = 'MODE' Mwet   'NAVIER_STOKES' DISCR                      ;
  122.   $dry   = 'MODE' Mdry   'NAVIER_STOKES' DISCR                      ;
  123.   $wall  = 'MODE' Mwall  'NAVIER_STOKES' DISCR                      ;
  124.   $dtop  = 'MODE' Mdtop  'NAVIER_STOKES' DISCR                      ;
  125.   volmt  = doma $mt         volume                                  ;
  126.  
  127.   mt    = 'DOMA' $mt    'MAILLAGE'                       ;
  128.   bot   = 'DOMA' $bot   'MAILLAGE'                       ;
  129.   axe   = 'DOMA' $axe   'MAILLAGE'                       ;
  130.   side  = 'DOMA' $side  'MAILLAGE'                       ;
  131.   top   = 'DOMA' $top   'MAILLAGE'                       ;
  132.   in    = 'DOMA' $in    'MAILLAGE'                       ;
  133.   cnt   = 'DOMA' $cnt   'MAILLAGE'                       ;
  134.   wet   = 'DOMA' $wet   'MAILLAGE'                       ;
  135.   dry   = 'DOMA' $dry   'MAILLAGE'                       ;
  136.   wall  = 'DOMA' $wall  'MAILLAGE'                       ;
  137.   dtop  = 'DOMA' $dtop  'MAILLAGE'                       ;
  138.  
  139. * Physical properties
  140.  
  141.   condf         =    2.62e-2                                        ;
  142.   nuf           =    1.50e-5                                        ;
  143.   rof           =    1.20e-0                                        ;
  144.   cpf           =    1.00e+3                                        ;
  145.   diff          =    condf/(rof*cpf)                                ;
  146.   nup           =    1.00e-20                                       ;
  147.   rop           =    2.50e+3                                        ;
  148.   cpp           =    8.40e+2                                        ;
  149.   difp          =    nup                                            ;
  150.   rrho          =    rof/rop                                        ;
  151.   qvol          =    0.00                                           ;
  152.   gx            =    0.00e+0                                        ;
  153.   gy            =   -9.81e+0                                        ;
  154.   gdr           =    0.0        (-9.81*(rop-rof)/rop)               ;
  155.   zero          =    0.0          0.0                               ;
  156.   beta          =    1.0          1.0                               ;
  157. * alpha         =   (3.14/6.0)*n_dot*(Dp**3.0)/Vi                   ;
  158.   alpha         =    C_s/rop/Vi                                     ;
  159.   n_dot         =   (6.0/3.14)*C_s/(rop*(Dp**3.0))                  ;
  160.   ustella       =    30.0*nuf/yp                                    ;
  161.  
  162.   L0            =     0.30                                          ;
  163.   K0            =     1.0e-4                                        ;
  164.   Rec           =    50.0                                           ;
  165.   Uref          =     4.                                            ;
  166.  
  167. * mess 'Cs' C_s 'Dp' Dp 'n_dot' n_dot ; opti donn 5                 ;
  168.  
  169.   tabbil           = table                                          ;
  170.   tabbil.'TF'      = mot 'TF'                                       ;
  171.   tabbil.'TP'      = mot 'TP'                                       ;
  172.   tabbil.'VP'      = mot 'VN'                                       ;
  173.   tabbil.'ALPHA'   = mot 'A'                                        ;
  174.   tabbil.'ROCPF'   = (rof*cpf)                                      ;
  175.   tabbil.'ROCPP'   = (rop*cpp)                                      ;
  176.   tabbil.'VOL'     = volmt                                          ;
  177.   tabbil.'CNT'     = $cnt                                           ;
  178.   tabbil.'TOP'     = $top                                           ;
  179.   tabbil.'BOT'     = $bot                                           ;
  180.  
  181. * Table for bif
  182.  
  183.   tabbif           = table                                          ;
  184.   tabbif.'RHOF'    = rof                                            ;
  185.   tabbif.'RHOP'    = rop                                            ;
  186.   tabbif.'DPART'   = Dp                                             ;
  187.   tabbif.'NUF'     = nuf                                            ;
  188.   tabbif.'ALPHA'   = mot 'A'                                        ;
  189.   tabbif.'UFLUID'  = mot 'UN'                                       ;
  190.   tabbif.'UPART'   = mot 'VN'                                       ;
  191.   tabbif.'KFLUID'  = mot 'KFL'                                      ;
  192.   tabbif.'KPART'   = mot 'KPR'                                      ;
  193.  
  194. * tabbif.'HPART'   = mot 'HP'                                       ;
  195. * tabbif.'HFLUID'  = mot 'HF'                                       ;
  196. * tabbif.'LAMBDAF' = condf                                          ;
  197. * tabbif.'ROCPP'   = (rop*cpp)                                      ;
  198. * tabbif.'ROCPF'   = (rof*cpf)                                      ;
  199. * tabbif.'TGASN'   = mot 'TF'                                       ;
  200. * tabbif.'TGASE'   = mot 'TFE'                                      ;
  201. * tabbif.'TPARTN'  = mot 'TP'                                       ;
  202. * tabbif.'TPARTE'  = mot 'TPE'                                      ;
  203.  
  204. * Points for time histories
  205.  
  206.   nc            =  10                                               ;
  207.   npti          =  nbno mt                                          ;
  208.   ps            =  npti  / nc                                       ;
  209.   pf            =  ps * nc                                          ;
  210.   lh       =  elem (doma $mt sommet) 'POI1' (lect  ps pas  ps  pf)  ;
  212.                          'UN' 1 lh 'UN' 2 lh
  213.                          'VN' 1 lh 'VN' 2 lh
  214.                          'KN'   lh 'EN'   lh
  215. *                        'TP'   lh 'TF'   lh
  216.                          'A'    lh
  217. ;
  218.  
  219. si ( complet )         ;
  220.    nbit = 10000        ;
  221. sinon                  ;
  222.    nbit =   800        ;
  223. *  nbit =   20         ;
  224. finsi                  ;
  225.  
  226. * Equations + boudary conditions
  227.  
  228. * Navier-Stokes + K-epsilon + particle mass conservation
  229.  
  230.   rv = EQEX  $mt   ITMA   nbit  ALFA  0.70
  231. *      ZONE  $mt   OPER NS   nuf           zero          INCO UN
  232.        OPTI 'NOCONS' 'SUPG' 'MMPG'
  233.        ZONE  $mt   OPER NSKE nuf  'NT'     zero          INCO UN KN EN
  234.        ZONE  $wall OPER FPU  nuf  'USTAR'  yp            INCO UN KN EN
  235.        ZONE  $mt   OPER NS   nup           gdr           INCO VN
  236.        OPTI 'CONS' 'SUPG'
  237.        ZONE  $mt   OPER TSCA nup  'VN' qvol 'TN'         INCO A
  238. ;
  239.  
  240. * Coupling terms on q.d.m equations
  241.  
  242.   rv = EQEX  rv
  243.        OPTI 'NOCONS' 'EFM1'
  244. *      OPTI 'IMPL'
  245.        ZONE  $mt   OPER FROT 'KFL' beta   'VN'           INCO UN
  246. *      OPTI 'IMPL'
  247.        ZONE  $mt   OPER FROT 'KPR' beta   'UN'           INCO VN
  248. ;
  249.   rv = EQEX  rv
  250.        ZONE  $mt   OPER FILTREKE Uref L0                 INCO KN EN
  251.        ZONE  $mt   OPER PPRE      'VN'     rrho          INCO VN
  252.        ZONE  $mt   OPER BIF  tabbif                      INCO VN
  253. ;
  254.  
  255. * Coupling terms on energy equations
  256.  
  257. * rv = EQEX  rv
  258. *      ZONE  $mt   OPER TSCA difp 'VN'     qvol          INCO TP
  259. *      ZONE  $mt   OPER TSCA diff 'UN'     qvol          INCO TF
  260. *      ZONE  $mt   OPER ECHI 'HP' 'TFE'                  INCO TP
  261. *      ZONE  $mt   OPER ECHI 'HF' 'TPE'                  INCO TF
  262. *      ZONE  $mt   OPER BILA tabbil
  263. ;
  264.   RV=EQEX RV OPTI EFM1 'CENTREE'
  265.   'ZONE' $MT      'OPER' 'DFDT' 1. 'UN' 'DELTAT' 'INCO' 'UN'
  266.   'ZONE' $MT      'OPER' 'DFDT' 1. 'VN' 'DELTAT' 'INCO' 'VN'
  267.   'ZONE' $MT      'OPER' 'DFDT' 1. 'KN' 'DELTAT' 'INCO' 'KN'
  268.   'ZONE' $MT      'OPER' 'DFDT' 1. 'EN' 'DELTAT' 'INCO' 'EN'
  269.   'ZONE' $MT      'OPER' 'DFDT' 1. 'A'  'DELTAT' 'INCO' 'A'
  270.   ;
  271.  
  272. * Boundary coditions
  273.  
  274.   rv = EQEX  rv
  275. *      CLIM  UN UIMP cnt  0.0     UN VIMP side 0.0
  276.        CLIM  UN VIMP top  0.0     UN VIMP bot  0.0
  277.        CLIM  UN UIMP axe  0.0     UN UIMP side 0.0
  278.        CLIM  VN UIMP axe  0.0     VN UIMP side 0.0
  279.        CLIM  VN VIMP in   (-1*Vi) VN VIMP dtop 0.0
  280.        CLIM  A  TIMP in   alpha   A  TIMP dtop 0.0
  281. *      CLIM  TP TIMP in  20.0
  282.  ;
  283.  
  284. * Pressure equation
  285.  
  286. RVP = EQEX 'OPTI' 'EF' KPRESS
  287.  'ZONE' $MT  OPER  KBBT  -1. betap  INCO 'UN' 'PRES'
  288. ;
  289.  
  290.     rvp.'METHINV'.TYPINV=1 ;
  291.     rvp.'METHINV'.IMPINV=0 ;
  292.     rvp.'METHINV'.NITMAX=300;
  293.     rvp.'METHINV'.PRECOND=3 ;                                       
  294.     rvp.'METHINV'.RESID  =1.e-8 ;
  295.     rvp.'METHINV' . 'FCPRECT'=100 ;
  296.     rvp.'METHINV' . 'FCPRECI'=100 ;
  297.  
  298.   RV.'PROJ'     =RVP ;
  299.  
  300. * Initial conditions
  301.  
  302.   alphaw          = kcht $wet  scal sommet (alpha  )                ;
  303.   tpwet           = kcht $wet  scal sommet (20.0   )                ;
  304.   rv.titolo       = titolo                                          ;
  305.   rv.INCO         = table INCO                                      ;
  306.   rv.INCO.'UN'    = kcht $mt   vect sommet (1.0e-5 1.0e-5)          ;
  307.   rv.INCO.'PRES'  = kcht $mt   scal kpress 0.                       ;
  308.   rv.INCO.'VN'    = kcht $mt   vect sommet (1.0e-5 (-1*Vi))         ;
  309.   rv.INCO.'KN'    = kcht $mt   scal sommet (1.0e-5       )          ;
  310.   rv.INCO.'EN'    = kcht $mt   scal sommet (1.0e-5       )          ;
  311.   rv.INCO.'NT'    = kcht $mt   scal centre                 nuf      ;
  312.   rv.INCO.'A'     = kcht $mt   scal sommet (0.0          ) alphaw   ;
  313.   rv.INCO.'TN'    = kcht $mt   scal sommet (0.0          )          ;
  314. * rv.INCO.'TP'    = kcht $mt   scal sommet (0.0          ) tpwet    ;
  315. * rv.INCO.'TF'    = kcht $mt   scal sommet (0.0          )          ;
  316. * rv.INCO.'TPE'   = kcht $mt   scal centre (0.0          )          ;
  317. * rv.INCO.'TFE'   = kcht $mt   scal centre (0.0          )          ;
  318.   rv.INCO.'KFL'   = kcht $mt   vect centre (0.0    0.0   )          ;
  319.   rv.INCO.'KPR'   = kcht $mt   vect centre (0.0    0.0   )          ;
  320. * rv.INCO.'HF'    = kcht $mt   scal centre (0.0          )          ;
  321. * rv.INCO.'HP'    = kcht $mt   scal centre (0.0          )          ;
  322.   rv.INCO.'USTAR' = kcht $wall scal centre                 ustella  ;
  323.   rv.'HIST'       = tran                                            ;
  324.   rv.'YP'         = yp                                              ;
  325.   rv.'NUF'        = nuf                                             ;
  326.  
  327. * Executing
  328.  
  329.   exec rv                                                           ;
  330.  
  331. * Saving results
  332.  
  333. * repix rv                                                          ;
  334. * sauv rv                                                           ;
  335.  
  336. * Procedura test
  337.  
  338.   debproc test                                                      ;
  339.   argu         rv*table complet*logique                             ;
  340.   pt1     = mt poin proc (0.0 1.260)                                ;
  341.   pt2     = mt poin proc (0.0 0.740)                                ;
  342.   u2      = exco uy (rv.inco.'UN') uy                               ;
  343.   v2      = exco uy (rv.inco.'VN') uy                               ;
  344.   upt1    = abs (extr (rv.inco.'UN') uy pt1)                        ;
  345.   upt2    = abs (extr (rv.inco.'UN') uy pt2)                        ;
  346.   vpt1    = abs (extr (rv.inco.'VN') uy pt1)                        ;
  347.   vpt2    = abs (extr (rv.inco.'VN') uy pt2)                        ;
  348.   apt1    = abs (extr (rv.inco.'A')  scal pt1)                      ;
  349.   apt2    = abs (extr (rv.inco.'A')  scal pt2)                      ;
  350.  
  351. si ( complet )                                                      ;
  352.   upt1_l  = 0.53580                                                 ;
  353.   upt2_l  = 0.59501                                                 ;
  354.   vpt1_l  = 3.0245                                                  ;
  355.   vpt2_l  = 3.5610                                                  ;
  356.   apt1_l  = 1.18535E-03                                             ;
  357.   apt2_l  = 1.00603E-03                                             ;
  358. sinon                                                               ;
  359.   upt1_l  = 0.43961                                                 ;
  360.   upt2_l  = 0.38635                                                 ;
  361.   vpt1_l  = 3.0859                                                  ;
  362.   vpt2_l  = 3.4686                                                  ;
  363.   apt1_l  = 6.78202E-04                                             ;
  364.   apt2_l  = 5.94036E-04                                             ;
  365. finsi                                                               ;
  366.  
  367.   upt1_e  = abs ((upt1 - upt1_l)/upt1_l)                            ;
  368.   upt2_e  = abs ((upt2 - upt2_l)/upt2_l)                            ;
  369.   vpt1_e  = abs ((vpt1 - vpt1_l)/vpt1_l)                            ;
  370.   vpt2_e  = abs ((vpt2 - vpt2_l)/vpt2_l)                            ;
  371.   apt1_e  = abs ((apt1 - apt1_l)/apt1_l)                            ;
  372.   apt2_e  = abs ((apt2 - apt2_l)/apt2_l)                            ;
  373.  
  374.   erlim = 5.00e-2                                                   ;
  375.   si (upt1_e > erlim)                                               ;
  376.      mess 'Error on Ugas : err_Upt1= ' upt1_e ' > ' erlim           ;
  377.      erreur 5                                                       ;
  378.   finsi                                                             ;
  379.  
  380.   si (upt2_e > erlim)                                               ;
  381.      mess 'Error on Ugas : err_Upt2= ' upt2_e ' > ' erlim           ;
  382.      erreur 5                                                       ;
  383.   finsi                                                             ;
  384.  
  385.   si (vpt1_e > erlim)                                               ;
  386.      mess 'Error on Vpart: err_Vpt1= ' vpt1_e ' > ' erlim           ;
  387.      erreur 5                                                       ;
  388.   finsi                                                             ;
  389.  
  390.   si (vpt2_e > erlim)                                               ;
  391.      mess 'Error on Vpart: err_Vpt2= ' vpt2_e ' > ' erlim           ;
  392.      erreur 5                                                       ;
  393.   finsi                                                             ;
  394.  
  395.   si (apt1_e > erlim)                                               ;
  396.      mess 'Error on alpha: err_Apt1= ' apt1_e ' > ' erlim           ;
  397.      erreur 5                                                       ;
  398.   finsi                                                             ;
  399.  
  400.   si (apt2_e > erlim)                                               ;
  401.      mess 'Error on alpha: err_Apt2= ' apt2_e ' > ' erlim            ;
  402.      erreur 5                                                       ;
  403.   finsi                                                             ;
  404.  
  405.   mess 'upt1'  upt1 'upt2' upt2 'vpt1' vpt1 'vpt2' vpt2             ;
  406.   mess 'apt1'  apt1 'apt2' apt2                                     ;
  407.   finproc                                                           ;
  408.  
  409. * Procedura visu
  410.  
  411.   debproc visu rv*table                                             ;
  412.   titr     rv.titolo                                                ;
  413.   tdom   = $mt                                                      ;
  414.   mesh   = doma $mt 'MAILLAGE'                                      ;
  415.   cmesh  = cont mesh                                                ;
  416.  
  417.   pt1    = poin mesh proc (0.0 0.740)                               ;
  418.   pt2    = poin mesh proc (0.4 0.740)                               ;
  419.   pt3    = poin mesh proc (0.0 1.260)                               ;
  420.   pt4    = poin mesh proc (0.4 1.260)                               ;
  421.  
  422.   pt12   = poin mesh droit pt2 pt1 0.0001                           ;
  423.   pt34   = poin mesh droit pt3 pt4 0.0001                           ;
  424.   meshl  = chan mesh lign                                           ;
  425.   line1  = elem meshl appuye pt12                                   ;
  426.   line2  = elem meshl appuye pt34                                   ;
  427.   line1  = inve line1                                               ;
  428.   line2  = inve line2                                               ;
  429.  
  430.   f1              = kcht $wet  scal sommet 1.0                      ;
  431.   filtro          = kcht $mt   scal sommet 0.0 f1                   ;
  432.  
  433.   un     = rv.inco.'UN'                                             ;
  434.   vn     = rv.inco.'VN'                                             ;
  435.   vn     = vn * filtro                                              ;
  436.   a      = rv.inco.'A'                                              ;
  437. * tp     = rv.inco.'TP'                                             ;
  438. * tf     = rv.inco.'TF'                                             ;
  439.   u1     = exco ux un             ux                                ;
  440.   u2     = exco uy un             uy                                ;
  441.   v1     = exco ux vn             ux                                ;
  442.   v2     = exco uy vn             uy                                ;
  443.  
  444.   unvt   = vect un 5.0e-2 ux uy jaun                                ;
  445.   vnvt   = vect vn 5.0e-3 ux uy roug                                ;
  446.  
  447.   titr (chaine titolo_g 'Uair')  ; trac unvt mesh cmesh             ;
  448.   titr (chaine titolo_g 'Upart') ; trac vnvt mesh cmesh             ;
  449.  
  450.   txa    = mot  'axe (m)'                                           ;
  451.   txo    = mot  'X   (m)'                                           ;
  452. * txl    = mot  'section 0.74 (m)'                                  ;
  453.   tyu    = mot  'velocity (m/s)'                                    ;
  454.   tyua   = mot  'air velocity (m/s)'                                ;
  455.   tyup   = mot  'u_I    (m/s)'                                      ;
  456.   tynt   = mot  'nu_t (m^2/s)'                                      ;
  457.   tya    = mot  'alpha'                                             ;
  458.   tyt    = mot  'temp. (C)'                                         ;
  459.   t1     = chaine titolo_g 'U;V /axe  '                             ;
  460.   t2     = chaine titolo_g 'Z=0.74 (m)'                             ;
  461.   t3     = chaine titolo_g 'Z=1.26 (m)'                             ;
  462.   t4     = chaine titolo_g 'Air veloc.'                             ;
  463.   t5     = chaine titolo_g 'Alpha/axe '                             ;
  464.   t6     = chaine titolo_g 'Alpha     '                             ;
  465.   t7     = chaine titolo_g 'Y+'                                     ;
  466.  
  467.   tab1   = table                                                    ;
  468.   tab1.1 = 'MARQ CROI REGU '                                        ;
  469.   tab1.2 = 'MARQ LOSA REGU '                                        ;
  470.   tab1.'TITRE'  = table                                             ;
  471.   tab1.'TITRE'. 1 = mot 'AIR'                                       ;
  472.   tab1.'TITRE'. 2 = mot 'PARTICLE'                                  ;
  473.  
  474.   tab2   = table                                                    ;
  475.   tab2.1 = 'MARQ LOSA REGU '                                        ;
  476.   tab2.2 = 'MARQ TRIA REGU '                                        ;
  477.   tab2.3 = 'MARQ CARR REGU '                                        ;
  478.   tab2.4 = 'MARQ ETOI REGU '                                        ;
  479.   tab2.'TITRE'   = table                                            ;
  480.   tab2.'TITRE'. 1 = mot 'H=1.50M'                                   ;
  481.   tab2.'TITRE'. 2 = mot 'H=1.00M'                                   ;
  482.   tab2.'TITRE'. 3 = mot 'H=0.50M'                                   ;
  483.   tab2.'TITRE'. 4 = mot 'H=0.74M'                                   ;
  484.  
  485.   tab3   = table                                                    ;
  486.   tab3.1 = 'MARQ LOSA REGU '                                        ;
  487.   tab3.2 = 'MARQ TRIA REGU '                                        ;
  488.   tab3.'TITRE'   = table                                            ;
  489.   tab3.'TITRE'. 1 = mot 'H=0.74M'                                   ;
  490.   tab3.'TITRE'. 2 = mot 'H=1.26M'                                   ;
  491.  
  492.   u2axe  = evol  chpo u2 uy    axe                                  ;
  493.   v2axe  = evol  chpo v2 uy    axe                                  ;
  494.   u2top  = evol  chpo u2 uy    top                                  ;
  495.   v2top  = evol  chpo v2 uy    top                                  ;
  496.   u2int  = evol  chpo u2 uy   (inve hori)                           ;
  497.   v2int  = evol  chpo v2 uy   (inve hori)                           ;
  498.   u2ln1  = evol  chpo u2 uy    line1                                ;
  499.   v2ln1  = evol  chpo v2 uy    line1                                ;
  500.   u2ln2  = evol  chpo u2 uy    line2                                ;
  501.   v2ln2  = evol  chpo v2 uy    line2                                ;
  502.  
  503. * tpaxe  = evol  chpo tp scal  axe                                  ;
  504. * tfaxe  = evol  chpo tf scal  axe                                  ;
  505. * tpln1  = evol  chpo tp scal  line1                                ;
  506. * tfln1  = evol  chpo tf scal  line1                                ;
  507.  
  508.   aaxe   = evol  chpo a  scal  axe                                  ;
  509.   atop   = evol  chpo a  scal (top      )                           ;
  510.   aln1   = evol  chpo a  scal  line1                                ;
  511.   aln2   = evol  chpo a  scal  line2                                ;
  512.   abot   = evol  chpo a  scal (inve bot )                           ;
  513.  
  514.   dess (u2axe et v2axe) mima titx txa tity tyu titr t1 gril lege tab1 ;
  515. * @excel1 u2axe 'pp'                              ;
  516.   dess (u2ln1 et v2ln1) mima titx txo tity tyu titr t2 gril lege tab1 ;
  517.   dess (u2ln2 et v2ln2) mima titx txo tity tyu titr t3 gril lege tab1 ;
  518.   dess (u2ln1 et u2ln2) mima titx txo tity tyua gril titr t4 lege tab3;
  519.   dess (aaxe          ) mima titx txa tity tya  gril titr t5        ;
  520.   dess (aln1  et aln2 ) mima titx txo tity tya  gril titr t6 lege tab3;
  521. * dess (tfaxe et tpaxe) mima titx txa tity tyt  gril lege tab1      ;
  522. * dess (tfln  et tpln ) mima titx txl tity tyt  gril lege tab1      ;
  523.  
  524.  
  525.   si (exist rv 'INCO' 'USTAR')                                      ;
  526.      yp    = rv.'YP'                                                ;
  527.      nuf   = rv.'NUF'                                               ;
  528.      ue    = rv.'INCO'.'USTAR'                                      ;
  529.      us    = elno $wall  ue                                         ;
  530.      Yplus = us * yp / nuf                                          ;
  531.      Ype   = evol chpo Yplus scal wall                              ;
  532.      dess Ype mima tity 'Y+' gril titr t7                           ;
  533.   finsi                                                             ;
  534.  
  535.   si (exist rv 'INCO' 'HP' )                                        ;
  536.      titr (chaine rv.titolo ' Hp'   )                               ;
  537.      trac (elno tdom (rv.'INCO'.'HP' )) mesh cmesh                  ;
  538.   finsi                                                             ;
  539.   si (exist rv 'INCO' 'HF' )                                        ;
  540.      titr (chaine rv.titolo ' Hf'   )                               ;
  541.      trac (elno tdom (rv.'INCO'.'HF' )) mesh cmesh                  ;
  542.   finsi                                                             ;
  543.   si (exist rv 'INCO' 'H'  )                                        ;
  544.      titr (chaine rv.titolo ' H '   )                               ;
  545.      trac (rv.'INCO'.'H')      mesh cmesh                           ;
  546.   finsi                                                             ;
  547.   si (exist rv 'INCO' 'TF')                                         ;
  548.      titr (chaine rv.titolo ' Tf'   )                               ;
  549.      trac (rv.inco.'TF')       mesh cmesh                           ;
  550.   finsi                                                             ;
  551.   si (exist rv 'INCO' 'TP')                                         ;
  552.      titr (chaine rv.titolo ' Tp'   )                               ;
  553.      trac (rv.inco.'TP')       mesh cmesh                           ;
  554.   finsi                                                             ;
  555.   si (exist rv 'INCO' 'A')                                          ;
  556.      titr (chaine rv.titolo ' alpha')                               ;
  557.      trac (rv.inco.'A' )       mesh cmesh                           ;
  558.   finsi                                                             ;
  559.   si (exist rv 'INCO' 'KN')                                         ;
  560.      kn = rv.inco.'KN'                                              ;
  561.      en = rv.inco.'EN'                                              ;
  562.      L0 = (kops (kops kn '**' 1.5) '*' (kops en '**' (-1.0)))       ;
  563.      titr (chaine rv.titolo ' kn')                                  ;
  564.      trac  kn                  mesh cmesh                           ;
  565.      titr (chaine rv.titolo ' en')                                  ;
  566.      trac  en                  mesh cmesh                           ;
  567.      mess 'min L0' (mini L0)  ; mess 'max L0' (maxi L0)             ;
  568.      titr (chaine rv.titolo ' L0')                                  ;
  569.      trac  L0                  mesh cmesh                           ;
  570. *    up   = kn * 2.0 / 3.0                                          ;
  571. *    up   = up ** 0.5                                               ;
  572. *    unor = kops un 'PSCA' un                                       ;
  573. **   up   = up / unor                                               ;
  574. *    upln = evol  chpo up scal  line1                               ;
  575. *    dess (upln          ) mima titx txl tity tyup gril             ;
  576.   finsi                                                             ;
  577.   si (exist rv 'INCO' 'NT')                                         ;
  578.      nt = elno tdom (rv.inco.'NT')                                  ;
  579.      titr (chaine rv.titolo ' nt')                                  ;
  580.      nte1   = evol  chpo nt   scal  line1                           ;
  581.      nte2   = evol  chpo nt   scal  line2                           ;
  582.      nte    = nte1 et nte2                                          ;
  583.      dess (nte         ) titx txo tity tynt gril lege tab3          ;
  584.      nr = nt / nuf                                                  ;
  585.      trac nt                   mesh cmesh                           ;
  586.      titr (chaine rv.titolo ' nr')                                  ;
  587.      trac nr                   mesh cmesh                           ;
  588.   finsi                                                             ;
  589.  
  590. * evoluzioni nel tempo
  591.  
  592.   si (exist rv 'HIST' '1UN')                                        ;
  593.      dess rv.hist.'1UN'                                             ;
  594.      dess rv.hist.'2UN'                                             ;
  595.   finsi                                                             ;
  596.   si (exist rv 'HIST' '1VN')                                        ;
  597.      dess rv.hist.'1VN'                                             ;
  598.      dess rv.hist.'2VN'                                             ;
  599.   finsi                                                             ;
  600.   si (exist rv 'HIST' 'KN')                                         ;
  601.      dess rv.hist.'KN'                                              ;
  602.      dess rv.hist.'EN'                                              ;
  603.   finsi                                                             ;
  604.   si (exist rv 'HIST' 'A' )                                         ;
  605.      dess rv.hist.'A'                                               ;
  606.   finsi                                                             ;
  607.   si (exist rv 'HIST' 'TF')                                         ;
  608.      dess rv.hist.'TF'                                              ;
  609.   finsi                                                             ;
  610.   si (exist rv 'HIST' 'TP')                                         ;
  611.      dess rv.hist.'TP'                                              ;
  612.   finsi                                                             ;
  613.  
  614.   finproc                                                           ;
  615.  
  616.  
  617.   test rv complet     ;
  618.  
  619. si ( 'EGA' graph 'O' );
  620.   visu rv             ;
  621. finsi                 ;
  622.   fin                                                               ;
  623.  
  624.  
  625.  
  626.  
  627.  
  628.  
  629.  
  630.  
  631.  
  632.  
  633.  
  634.  
  635.  
  636.  
  637.  
  638.  
  639.  
  640.  

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