1 : C TURB6 SOURCE CHAT 05/01/13 03:57:03 5004
2 : C***********************************************************************
3 : C* *
4 : C* PROJET : Opérateur LIMI *
5 : C* NOM : TURB6 *
6 : C* DESCRIPTION : Ce sous-programme calcule les relations de fermeture *
7 : C* des méthodes intégrales pour les couches limites *
8 : C* turbulentes. Utilisation des formules de Head et *
9 : C* Ludwieg&Tillmann. *
10 : C* LANGAGE : Esope *
11 : C* AUTEUR : Guillaume VENCO - DRN/DMT/SEMT/LTMF *
12 : C* *
13 : C***********************************************************************
14 : C* *
15 : C* APPELES : Aucun *
16 : C* *
17 : C***********************************************************************
18 : C* *
19 : C* ENTREES : D2N : -Epaisseur de quantité de mouvement, *
20 : C* correspond à l'inconnue de l'équation intégrale*
21 : C* de quantité de mouvement *
22 : C* Y1N : -Différence entre delta et delta1, *
23 : C* correspond à l'inconnue de l'équation intégrale*
24 : C* d'entrainement *
25 : C* UEN : -Champ de vitesse à l'extérieur de la couche *
26 : C* limite *
27 : C* N : -Nombre de noeuds du maillage *
28 : C* *
29 : C* SORTIES : HH : -Facteur H*, correspond au rapport Y1N/D2N *
30 : C* HN : -Facteur de forme *
31 : C* CFN : -Coefficient de frottement *
32 : C* CEN : -Coefficient d'entraînement *
33 : C* D1N : -Epaisseur de déplacement *
34 : C* *
35 : C***********************************************************************
36 : C* *
37 : C* VERSION : 20/04/2000 *
38 : C* CREATION : *
39 : C* *
40 : C***********************************************************************
41 :
42 : SUBROUTINE TURB6(D2N,Y1N,UEN,HH,HN,CFN,CEN,D1N,N)
43 :
44 : C* *********************************
45 : C* *** Déclaration des variables ***
46 : C* *********************************
47 :
48 : IMPLICIT INTEGER(I-N)
49 : IMPLICIT REAL*8 (A-H,O-Z)
50 : REAL*8 D2N(N),Y1N(N),UEN(N)
51 : REAL*8 HH(N),HN(N),CFN(N),CEN(N),D1N(N),NU
52 :
53 : C* *********************************
54 : C* ******** Sous-Programme *********
55 : C* *********************************
56 :
57 : C---- Viscosité du fluide
58 : NU = 1e-6
59 :
60 :
61 : DO 10 I=1,N
62 : c write(6,*) 'Ue(',I,')=',UEN(I)
63 : C---- facteur H*
64 : HH(I)=Y1N(I)*(D2N(I)**(-1))
65 : c write(6,*) 'HH(',I,')=',HH(I)
66 : C---- Coefficient d'entrainement
67 : CEN(I)=0.0306*((HH(I)-3.)**(-0.653))
68 : c write(6,*) 'CEN(',I,')=',CEN(I)
69 : C---- Détermination du facteur de forme H
70 : HN(I)=0.7+(((1./1.535)*(HH(I)-3.3))**(-1./2.715))
71 : c write(6,*) 'H(',I,')=',HN(I)
72 : C---- Coefficient de frottement
73 : CFN(I)=0.246*(10**(-0.678*HN(I)))
74 : & *((UEN(I)*D2N(I)/NU)**(-0.268))
75 : c write(6,*) 'CFN(',I,')=',CFN(I)
76 : C---- Calcul de D1N
77 : D1N(I) = HN(I)*D2N(I)
78 : 10 CONTINUE
79 :
80 :
81 : RETURN
82 : END
83 :
84 : C* *********************************
85 : C* **** Fin Sous-Programme *******
86 : C* *********************************
87 :
88 :
89 :
90 :
91 :
92 :
93 :
94 :
95 :
© Cast3M 2003 - All rights reserved.
Disclaimer