* 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 * * * * 'ARGUMENT' mquaf*'MAILLAGE' ; 'ARGUMENT' un*'CHPOINT' ; 'ARGUMENT' clphi/'CHPOINT' ; * 'SI' ('NON' ('EXISTE' clphi)) ; 'FINSI' ; * 'ARGUMENT' rvm/'TABLE' ; 'SI' ('NON' ('EXISTE' rvm)) ; rvm = rv . 'METHINV' ; 'FINSI' ; * * 'SI' ('NEG' dim 2) ; 'ERREUR' 219 ; 'FINSI' ; * * 'SI' laxi ; cdpr = '*' ('COORDONNEE' 1 mquaf) ('*' PI 2.D0) ; * numop = 2 ; numvar = 1 ; numdat = 0 ; numcof = 0 ; numder = 2 ; * A . 'VAR' . 1 . 'DISC' = disc ; * numdat = 1 ; numcof = 1 ; B . 'VAR' . 1 . 'DISC' = disc ; B . 'DAT' . 1 . 'DISC' = disc ; B . 'DAT' . 1 . 'VALEUR' = cdpr ; B . 'COF' . 1 . 'COMPOR' = 'IDEN' ; * * numvar = 2 ; 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 . 'DISC' = disc ; C . 'DAT' . 1 . 'VALEUR' = cdpr ; C . 'COF' . 1 . 'COMPOR' = 'IDEN' ; * * * 'SINON' ; * * numop = 2 ; numvar = 1 ; numdat = 0 ; numcof = 0 ; numder = 2 ; * A . 'VAR' . 1 . 'DISC' = disc ; * numvar = 2 ; 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 ; * * 'FINSI' ; * 'RESPRO' psi ; * 'FINPROC' ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales