Notices | Cast3M

$$$$ FCOURANT NOTICE  GOUNAND   18/10/25    21:15:01     9973           
                                             DATE     18/10/25

    Procedure FCOURANT
    ------------------

    CHPO1 = FCOURANT MAIL1 CHPO2 (RIG1 (CHPO3)) (TAB1)   ;


    Objet :
    _______

    Cette procedure calcule la fonction de courant scalaire, en 2D et 2D
    axisymetrique, correspondant au champ de vitesse a divergence nul
    donne.

    Commentaire :
    _____________

      MAIL1 : Domaine de definition de la vitesse (type MAILLAGE)

      CHPO2 : Champ de vitesse (type CHPOINT)

      RIG1  : Conditions sur la fonction de courant
      CHPO3   (Par defaut, on bloque le premier noeud de MAIL1)

      TAB1  : table optionnelle precisant les options pour le solveur
              de systeme lineaire KRES (cf. notice KRES)

      CHPO1 : Fonction de courant (type CHPOINT, nom de composante
              'PSI')

    Remarques :
    ___________

      1) Faute de modele, les conditions sur la fonction de courant
         (RIG1) doivent porter sur une inconnue de nom 'T'.

      2) Si TAB1 n'est pas presente, on utilise RESO pour resoudre le
         systeme lineaire.

    Notes :
    _______
    
    On utilise les definitions suivantes pour la fonction de courant psi
    en fonction du champ de vitesse u :

      1) En 2D plan : dpsi/dx =   u_y
                      dpsi/dy = - u_x

      2) En 2D axisymetrique : dpsi/dr = ( 2pi r) u_z
                               dpsi/dz = (-2pi r) u_r

    Ces definitions sont compatibles avec l'interpretation suivante :
    psi(B) - psi(A) = debit volumique (m^2.s-1 en 2D plan, m^3.s-1 en 2D 
    axi) traversant le segment [AB]. Le sens positif est celui du
    vecteur faisant un angle de +90 degres avec le vecteur AB.

    En 2D axi, psi est appelee fonction de courant de Stokes.
    Elle est egale a la composante hors plan du potentiel vecteur
    associe a u dans la decomposition de Helmholtz, divisee par la
    coordonnee radiale r.

    Si le champ de vitesse u n'est pas approximativement a divergence
    nulle, les isovaleurs de la fonction de courant psi ne seront pas
    necessairement tangentes aux vecteurs vitesse.