Procédures | Cast3M

* FCOURANT  PROCEDUR  GOUNAND   07/12/21    21:15:01     6010           
************************************************************************
* NOM         : FCOURANT
* DESCRIPTION : Calcul la fonction de courant en 2D et 2D Axi
*               par une méthode d'éléments finis moindres carrés
*
* LANGAGE     : GIBIANE-CAST3M
* AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
*               mél : gounand@semt2.smts.cea.fr
**********************************************************************
* VERSION    : v1, 22/11/2007, version initiale
* HISTORIQUE : v1, 22/11/2007, création
* HISTORIQUE :
* HISTORIQUE :
************************************************************************
*
* Procedure fonction de courant 2D et 2D axi
* - modns : modele Navier-Stokes
* - disc  : discretisation des champs au sommet (QUAF ou LINE) 
* - un    : champ de vitesses
* - clphi : conditions de Dirichlet sur phi
*           (nom de composante 'SCAL')
*
* en 2D, on résout : dpsi/dx = u_y
*                    dpsi/dy = -u_x
* en 2D axi , on résout : dpsi/dr = 2pi r u_z
*                         dpsi/dz = - 2 pi r u_r
*
*
*
'DEBPROC' FCOURANT ;
*
'ARGUMENT' mquaf*'MAILLAGE' ;
'ARGUMENT' disc*'MOT' ;
'ARGUMENT' un*'CHPOINT' ;
'ARGUMENT' clphi/'CHPOINT' ;
*
'SI' ('NON' ('EXISTE' clphi)) ;
   pref  = 'POIN' ('CHANGER' mquaf 'POI1') 1 ;
   clphi = 'MANUEL' 'CHPO' pref 1 'SCAL' 0.D0 ;
'FINSI' ;
*
'ARGUMENT' rvm/'TABLE' ;
'SI' ('NON' ('EXISTE' rvm)) ;
   rv = 'EQEX' ;
   rvm = rv . 'METHINV' ;
'FINSI' ;   
*
dim = 'VALEUR' 'DIME' ;
*
'SI' ('NEG' dim 2) ;
   'ERREUR' 219 ;
'FINSI' ;
*
laxi = 'EGA' ('VALEUR' 'MODE') 'AXIS' ;
ncx = 'MOTS' 'UX' ; ncy = 'MOTS' 'UY' ;
*
'SI' laxi ;
   cdpr  = '*' ('COORDONNEE' 1 mquaf) ('*' PI 2.D0) ;
   ur = 'EXCO' ncx un ;
   uz = 'EXCO' ncy un ;
*
   numop  = 2 ; numvar = 1 ; numdat = 0 ;
   numcof = 0 ; numder = 2 ;
*
   A = ININLIN numop numvar numdat numcof numder ;
   A . 'VAR' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
   A . 'VAR' . 1 . 'DISC'   = disc ;
   A . 1 . 1 . 1 = 'LECT' ;
   A . 2 . 1 . 2 = 'LECT' ;
*
   numdat = 1 ; numcof = 1 ;
   B = ININLIN numop numvar numdat numcof numder ;
   B . 'VAR' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
   B . 'VAR' . 1 . 'DISC'   = disc ;
   B . 'DAT' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
   B . 'DAT' . 1 . 'DISC'   = disc ;
   B . 'DAT' . 1 . 'VALEUR' = cdpr ;
   B . 'COF' . 1 . 'COMPOR' = 'IDEN' ;
   B . 'COF' . 1 . 'LDAT'   = 'LECT' 1 ;
*   
   B . 1 . 1 . 1 = 'LECT' 1 ;
   B . 2 . 1 . 2 = 'LECT' 1 ;
*
   numvar = 2 ;
   C = ININLIN numop numvar numdat numcof numder ;
   C . 'VAR' . 1 . 'NOMDDL' = ncx ;
   C . 'VAR' . 1 . 'DISC'   = disc ;
   C . 'VAR' . 1 . 'VALEUR' = '*' ur -1. ;
   C . 'VAR' . 2 . 'NOMDDL' = ncy ;
   C . 'VAR' . 2 . 'DISC'   = disc ;
   C . 'VAR' . 2 . 'VALEUR' = uz ;
   C . 'DAT' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
   C . 'DAT' . 1 . 'DISC'   = disc ;
   C . 'DAT' . 1 . 'VALEUR' = cdpr ;
   C . 'COF' . 1 . 'COMPOR' = 'IDEN' ;
   C . 'COF' . 1 . 'LDAT'   = 'LECT' 1 ;
*
   C . 1 . 2 . 0 = 'LECT' 1 ;
   C . 2 . 1 . 0 = 'LECT' 1 ;
*
   mat = NLIN mquaf disc A B 'GAU7' ;
   mat = 'KOPS' 'RIMA' mat ;
   smb = NLIN mquaf disc C B 'GAU7' ;
*
'SINON'  ;
*
   ux = 'EXCO' ncx un ;
   uy = 'EXCO' ncy un ;
*
   numop  = 2 ; numvar = 1 ; numdat = 0 ;
   numcof = 0 ; numder = 2 ;
*
   A = ININLIN numop numvar numdat numcof numder ;
   A . 'VAR' . 1 . 'NOMDDL' = 'MOTS' 'SCAL' ;
   A . 'VAR' . 1 . 'DISC'   = disc ;
   A . 1 . 1 . 1 = 'LECT' ;
   A . 2 . 1 . 2 = 'LECT' ;
*
   numvar = 2 ;
   B = ININLIN numop numvar numdat numcof numder ;
   B . 'VAR' . 1 . 'NOMDDL' = ncx ;
   B . 'VAR' . 1 . 'DISC'   = disc ;
   B . 'VAR' . 1 . 'VALEUR' = '*' ux -1. ;
   B . 'VAR' . 2 . 'NOMDDL' = ncy ;
   B . 'VAR' . 2 . 'DISC'   = disc ;
   B . 'VAR' . 2 . 'VALEUR' = uy ;
*
   B . 1 . 2 . 0 = 'LECT' ;
   B . 2 . 1 . 0 = 'LECT' ;
*
   mat = NLIN mquaf disc A A 'GAU7' ;
   mat = 'KOPS' 'RIMA' mat ;
   smb = NLIN mquaf disc B A 'GAU7' ;
'FINSI' ;
*
psi = 'KRES' mat 'TYPI' rvm 'CLIM' clphi 'SMBR' smb ;
'RESPRO' psi ;
*
'FINPROC' ;