Accueil | Cast3M

   1 : $$$$ KMAB     NOTICE  CHAT      11/09/12    21:16:44     7124           
   2 :                                              DATE     11/09/12
   3 :                                                                                 
   4 :   Operateur KMAB                           Voir aussi : NAVI KBBT
   5 :     --------------                                        DUDW  EQEX  
   6 :                                                            
   7 :     SYNTAXE ( EQEX ) :   Cf operateur EQEX                 
   8 :     __________________
   9 : 
  10 : 
  11 :       'OPER'  'KMAB' coef <beta>     'INCO' 'UN' 'PRES'
  12 : 
  13 : 
  14 :                                                                                 
  15 :                                                                                 
  16 :     OBJET :
  17 :     -----
  18 : 
  19 :  L'operateur KMAB discretise le terme Div U du systeme d'equations de
  20 : Stokes (ou Navier-Stokes) par une methode d'elements finis.
  21 : 
  22 :     Commentaires
  23 :     ____________
  24 : 
  25 :      coef    coefficent multiplicateur
  26 :              FLOTTANT
  27 :              ou CHPOINT SCAL SOMMET  (porosite volumique)
  28 :              ou CHPOINT VECT SOMMET  (porosite directionnelle)
  29 :              ou MOT
  30 : 
  31 :      beta    parametre de stabilisation pour les elements lineaires
  32 :              FLOTTANT ou MOT
  33 : 
  34 :      UN      Champ de vitesse
  35 :              CHPOINT VECT SOMMET ou MOT
  36 : 
  37 :      PRES    Champ de pression
  38 :              CHPOINT SCAL CENTRE   ou MOT
  39 :              CHPOINT SCAL CENTREP1 ou MOT
  40 :              CHPOINT SCAL CENTREP0 ou MOT
  41 :              le type doit etre precise dans les options mot cle INCOD
  42 : 
  43 :  Un coefficient de type MOT indique que l'operateur va chercher le
  44 :  champ dans la table INCO a l'indice MOT.
  45 : 
  46 : 
  47 :     Complement d'information :
  48 :     __________________________
  49 : 
  50 : 
  51 :   1/ La formulation variationnelle est :
  52 : 
  53 :    /               
  54 :    | q*Div U dv =  
  55 :    |v            
  56 : 
  57 :   2/ Si toutes les vitesses normales sont imposees (nulles ou non) sur
  58 :      les frontieres, il faut IMPERATIVEMENT imposer la pression en un
  59 :      point. C'est le cas pour tout ecoulement d'un fluide incompressible
  60 :      en cavite fermee.
  61 : 
  62 :   3/ Si le coefficient est de type CHPOINT SCAL SOMMET
  63 :     (Porosite volumique H) On calcule :
  64 : 
  65 :    /                 /                 /
  66 :    | q*Div(H U) dv = | q H Div U dv  + | q U Grad H  dv
  67 :    |v                |v                |v
  68 : 
  69 : 
  70 :   4/ Si le coefficient est de type CHPOINT VECT SOMMET
  71 :     (Porosite surfacique ou directionnelle Hi ) On calcule :
  72 : 
  73 :    /
  74 :    | q*Div((HU)i) dv
  75 :    |v
  76 : 
  77 : 
  78 :     Options : (EQEX)
  79 :     _________
  80 : 
  81 :     OPTI INCOD CENTRE
  82 :                CENTREP1
  83 :                CENTREP0
  84 :