1 : C TURB5 SOURCE CHAT 05/01/13 03:56:59 5004
2 : C***********************************************************************
3 : C* *
4 : C* PROJET : Opérateur LIMI *
5 : C* NOM : TURB5 *
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 Cousteix. *
9 : C* LANGAGE : Esope *
10 : C* AUTEUR : Guillaume VENCO - DRN/DMT/SEMT/LTMF *
11 : C* *
12 : C***********************************************************************
13 : C* *
14 : C* APPELES : lEBAl2.eso, sous-programme *
15 : C* RACFGA.eso, sous-programme *
16 : C* JACOB.eso, sous-programme *
17 : C* equa1 , fonction *
18 : C* equa2 , fonction *
19 : C* equa3 , fonction *
20 : C* *
21 : C***********************************************************************
22 : C* *
23 : C* ENTREES : D2N : -Epaisseur de quantité de mouvement, *
24 : C* correspond à l'inconnue de l'équation intégrale*
25 : C* de quantité de mouvement *
26 : C* Y1N : -Différence entre delta et delta1, *
27 : C* correspond à l'inconnue de l'équation intégrale*
28 : C* d'entrainement *
29 : C* UEN : -Champ de vitesse à l'extérieur de la couche *
30 : C* limite *
31 : C* N : -Nombre de noeuds du maillage *
32 : C* *
33 : C* SORTIES : HH : -Facteur H*, correspond au rapport Y1N/D2N *
34 : C* HN : -Facteur de forme *
35 : C* CFN : -Coefficient de frottement *
36 : C* CEN : -Coefficient d'entraînement *
37 : C* D1N : -Epaisseur de déplacement *
38 : C* DeN : -Epaisseur de couche limite *
39 : C* *
40 : C***********************************************************************
41 : C* *
42 : C* VERSION : 20/04/2000 *
43 : C* CREATION : *
44 : C* *
45 : C***********************************************************************
46 :
47 : SUBROUTINE TURB5(D2N,Y1N,UEN,HH,HN,CFN,CEN,D1N,DeN,N)
48 :
49 : C* *********************************
50 : C* *** Déclaration des variables ***
51 : C* *********************************
52 :
53 : IMPLICIT INTEGER(I-N)
54 : IMPLICIT REAL*8 (A-H,O-Z)
55 : REAL*8 D2N(N),Y1N(N),UEN(N)
56 : REAL*8 HH(N),HN(N),CFN(N),CEN(N),D1N(N),DeN(N)
57 : REAL*8 F,P,NU,GAN,GN,A0,B0,be
58 : REAL*8 WX(3,3),QX(3,3)
59 : REAL*8 RD1,X1,X2,X3,XP1,XP2,XP3,f1P,f2P,f3P,r1,r2,r3,R
60 : PARAMETER (eps=1.d-5)
61 :
62 :
63 : C* *********************************
64 : C* ******** Sous-programme *********
65 : C* *********************************
66 :
67 : C---- Viscosité du fluide
68 : NU = 1e-6
69 :
70 : DO 10 I=1,N
71 : c write(6,*) 'Ue(',I,')=',UEN(I)
72 : C---- Facteur H*
73 : HH(I)=Y1N(I)*(D2N(I)**(-1))
74 : c write(6,*) 'HH(',I,')=',HH(I)
75 :
76 : C---- Calcul des coefficients H,gamma et G par itérations
77 : C---- Méthode de Newton
78 :
79 : C---- Valeurs initiales calculées à partir de la
80 : C---- formule approchée
81 : A0=1
82 : B0=2
83 : be=0.613
83 : CALL LEBAL2(HH(I),be,HN(I),A0,B0)
84 : D1N(I)=D2N(I)*HN(I)
85 : A0=1.D-3
86 : B0= 10
86 : CALL RACFGA(UEN(I),D1N(I),HN(I),GAN,A0,B0)
87 : GN=(HN(I)-1)/(HN(I)*GAN)
88 : RD1=UEN(I)*D1N(I)/NU
89 : c write(6,*) 'D1N(',I,')=',D1N(I)
90 : c write(6,*) 'H0(',I,')=',HN(I)
91 : c write(6,*) 'GN0(',I,')=',GN
92 : c write(6,*) 'Gamma0(',I,')=',GAN
93 :
94 : C---- Calcul des coefficients réels par itérations
95 : X1=HN(I)
96 : X2=GN
97 : X3=GAN
98 : K=0
98 : 40 CALL JACOB(X1,X2,X3,WX,QX)
99 : K=K+1
100 : c write(6,*) 'it=',K
101 : f1P= equa1(X1,X2,HH(I))
102 : f2P= equa2(X2,X3,RD1)
103 : f3P= equa3(X1,X2,X3)
104 :
105 : XP1=X1-(QX(1,1)*f1P+QX(1,2)*f2P+QX(1,3)*f3P)
106 : XP2=X2-(QX(2,1)*f1P+QX(2,2)*f2P+QX(2,3)*f3P)
107 : XP3=X3-(QX(3,1)*f1P+QX(3,2)*f2P+QX(3,3)*f3P)
108 : c write(6,*) 'X1(P+1)=',XP1
109 : c write(6,*) 'X2(P+1)=',XP2
110 : c write(6,*) 'X3(P+1)=',XP3
111 : r1=abs(XP1-X1)
112 : r2=abs(XP2-X2)
113 : r3=abs(XP3-X3)
114 :
115 : R=max(r1,r2,r3)
116 : IF (R.GE.eps) GOTO 20
117 : GOTO 30
118 :
119 : 20 X1=XP1
120 : X2=XP2
121 : X3=XP3
122 : GOTO 40
123 :
124 : 30 HN(I)=XP1
125 : GN=XP2
126 : GAN=XP3
127 : c write(6,*) 'H',HN(I)
128 : c write(6,*) 'gamma=',GAN
129 : c write(6,*) 'G=',GN
130 :
131 : C---- Valeurs des formules de similitude F(G) et P(G)
132 : P=0.074*GN-1.0957*(GN**(-1))
133 : F=0.613*GN-((3.6+76.86*((GN**(-1))-0.154)**2)*(GN**(-1)))
134 :
135 : C---- Coefficient de frottement
136 : CFN(I)=2*(GAN**2)
137 : c write(6,*) 'Cf=',CFN(I)
138 :
139 : C---- Coefficient d'entrainement Ce
140 : CEN(I) = P*GAN
141 : c write(6,*) 'CE=',CEN(I)
142 :
143 : C---- Epaisseur de couche limite
144 : DeN(I)=D1N(I)*((F*GAN)**(-1))
145 : c write(6,*) 'delta=',DeN(I)
146 :
147 : 10 CONTINUE
148 :
149 : RETURN
150 : END
151 : C* *********************************
152 : C* **** Fin Sous-Programme *******
153 : C* *********************************
154 :
155 :
156 :
157 :
158 :
159 :
160 :
© Cast3M 2003 - All rights reserved.
Disclaimer