* CALCDISP PROCEDUR GBM 04/12/22 21:15:01 4973 *--------------------------------------------------------------------- *---------PROCEDURE DE CALCUL DE LA DISPERSIVITE---------------------- *--------------------------------------------------------------------- * |-----------------------------------------------------------------| * | Généralités : CALCDISP calcule le tenseur de dispersion | * | du problème de transport convection-diffusion. | * |-----------------------------------------------------------------| * | | * |-----------------------------------------------------------------| * | ENTREES | * |-----------------------------------------------------------------| * | | * | ALPHAL : coefficient longitudinal de dispersivité | * | | * | ALPHAT : coefficient transverse de dispersivité | * | | * | QELEM : vitesse au centre de chaque élément | * |-----------------------------------------------------------------| * | SORTIES | * |-----------------------------------------------------------------| * | | * | | * | dif_disp : C'est le tenseur de dispersion | * | | * | | ********************************************************************** * Seulement si présence de convection et présence de * dispersivité et disp_l et disp_t 'SI' (DISPERSI ) ; * vpetit = 1.D-20 sur le meme nombre de zone que DISPL * on calcul la dispersivité * dimension 2 VC = (V1 * V1) '+' (V2 * V2) ; VC = (VC**0.5) ; VC = VC + vpetit ; VCM1 = VC '**' -1 ; D11P = ((DISPL - DISPT) * V1*V1*VCM1) + (DISPT * VC) ; D22P = ((DISPL - DISPT) * V2*V2*VCM1) + (DISPT * VC) ; D21P = (DISPL - DISPT) * V1*V2*VCM1 ; dif_disp = D11P '+' D22P '+' D21P ; 'FINSI' ; * dimension 3 VC = (V1 * V1) '+' (V2 * V2) '+' (V3 * V3) ; VC = VC**0.5 ; VC = VC + vpetit ; VCM1 = VC '**' -1 ; D11P = ((DISPL - DISPT) * V1*V1*VCM1) + (DISPT * VC) ; D22P = ((DISPL - DISPT) * V2*V2*VCM1) + (DISPT * VC) ; D33P = ((DISPL - DISPT) * V3*V3*VCM1) + (DISPT * VC) ; D21P = (DISPL - DISPT) * V1*V2*VCM1 ; D31P = (DISPL - DISPT) * V1*V3*VCM1 ; D32P = (DISPL - DISPT) * V2*V3*VCM1 ; dif_disp = D11P '+' D22P '+' D33P '+' D21P '+' D31P '+' D32P ; 'FINSI' ; 'FINSI' ; 'FINPROC' dif_disp ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales