1 : C PRE5 SOURCE BECC 05/12/02 21:16:34 5258
2 : SUBROUTINE PRE5()
3 : C************************************************************************
4 : C
5 : C PROJET : CASTEM 2000
6 : C
7 : C NOM : PRE5
8 : C
9 : C DESCRIPTION : Voir PRET
10 : C
11 : C Resolution des Equations de la modele bifluide
12 : C
13 : C Calculation of the left and right states
14 : C at the interface for eact t step
15 : C
16 : C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
17 : C
18 : C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
19 : C Modified for two-fluid flow by
20 : C Jose R. Garcia Cascales
21 : C
22 : C************************************************************************
23 : C
24 : C
25 : C APPELES (Outils) : LIRENT, ERREUR
26 : C
27 : C APPELES (Calcul) :
28 : C
29 : C PRE1-------------- PRE12F (1st order)
30 : C |
31 : C |
32 : C --------------- PRE22F (2nd order)
33 : C
34 : C************************************************************************
35 : C
36 : C PHRASE D'APPEL (GIBIANE) :
37 : C
38 : C 1) two-fluid flow
39 : C
40 : C a) 1st order
41 : C
42 : C IALPHF, IUVF, IULF, IPF, ITVF, ITLF, IRVF, IRLF = 'PRET'
43 : C 'MULTIFLU' ORDESP ORDTPS $DOMA
44 : C VCH1 VCH2 VCH3 VCH4 VCH5 VCH6 VCH7 VCH8 ;
45 : C
46 : C CHAM1S CHAM2S CHAM3S CHAM4S CHAM5S CHAM6S CHAM7S CHAM8S = 'PRET'
47 : C 'MULTIFLU' ENTI1 ENTI2 TAB1
48 : C CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 CHPO7 CHPO8 ;
49 : C
50 : C
51 : C ENTREES :
52 : C
53 : C 'MULTIFLU' : key name (two phase flow)
54 : C
55 : C ENTI1 : ordre en espace (=1)
56 : C
57 : C ENTI2 : ordre en temps (=1 si ENTI1 = 1)
58 : C
59 : C TAB1 : TABLE de sous type "DOMAINE"
60 : C
61 : C CHPO1 : CHPOINT "CENTRE" void fraction
62 : C
63 : C CHPO2 : CHPOINT "CENTRE" vapour velocity
64 : C (NDIM components)
65 : C
66 : C CHPO3 : CHPOINT "CENTRE" liquid velocity
67 : C (NDIM components)
68 : C
69 : C CHPO4 : CHPOINT "CENTRE" pressure
70 : C
71 : C CHPO5 : CHPOINT "CENTRE" vapour temperature
72 : C
73 : C CHPO6 : CHPOINT "CENTRE" liquid temperature
74 : C
75 : C CHPO7 : CHPOINT "CENTRE" vapour density
76 : C
77 : C CHPO8 : CHPOINT "CENTRE" liquid density
78 : C
79 : C SORTIES :
80 : C
81 : C CHAM1S : MCHAML "FACEL" void fraction
82 : C (a gauche et a droite de chaque face).
83 : C Only one component ('SCAL')
84 : C
85 : C CHAM2S : MCHAML "FACEL" vapour velocity and the
86 : C director cosines (n,t) in the corresponding face;
87 : C in the 2D case 6 composantes:
88 : C 'UVN' = normal velocity
89 : C 'UVT' = tangent velocity
90 : C 'NX' = n.x
91 : C 'NY' = n.y
92 : C 'TX' = t.x
93 : C 'TY' = t.y
94 : C
95 : C CHAM3S : MCHAML "FACEL" liquid velocity
96 : C in the 2D case 2 composantes:
97 : C 'ULN' = normal velocity
98 : C 'ULT' = tangent velocity
99 : C
100 : C CHAM4S : MCHAML "FACEL" pressure
101 : C Only one component ('SCAL')
102 : C
103 : C CHAM5S : MCHAML "FACEL" vapour temperature
104 : C Only one component ('SCAL')
105 : C
106 : C CHAM6S : MCHAML "FACEL" liquid temperature
107 : C Only one component ('SCAL')
108 : C
109 : C CHAM7S : MCHAML "FACEL" vapour density
110 : C Only one component ('SCAL')
111 : C
112 : C CHAM8S : MCHAML "FACEL" liquid density
113 : C Only one component ('SCAL')
114 : C
115 : C************************************************************************
116 : C
117 : C HISTORIQUE (Anomalies et modifications éventuelles)
118 : C
119 : C HISTORIQUE : Créée le 11.6.98.(Modified February 26th, 2002)
120 : C
121 : C************************************************************************
122 : C
123 : C
124 : C**** Variables de COOPTIO
125 : C
126 : INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
127 : & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
128 : & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
129 : & ,IECHO, IIMPI, IOSPI
130 : & ,IDIM, IFICLE, IPREFI
131 : & ,MCOORD
132 : & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
133 : & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
134 : & ,NORINC,NORVAL,NORIND,NORVAD
135 : & ,NUCROU, IPSAUV
136 : C
137 : C**** Les variables
138 : C
139 : INTEGER ICOND, IRETOU
140 : & ,ORDESP, ORDTEM
141 : C
142 : C**** Les Includes
143 : C
144 : -INC CCOPTIO
145 : C
146 : C**** Les mots clés
147 : C
148 : C
149 : C**** Initialisation des parametres d'erreur
150 : C
151 : MOTERR(1:40) = ' '
152 : C
153 : C**** Lecture ORDESP
154 : C
155 : ICOND = 1
155 : CALL LIRENT(ORDESP,ICOND,IRETOU)
156 : IF(IERR .NE. 0) GOTO 9999
157 : IF((ORDESP .LT. 1) .OR. (ORDESP .GT. 2))THEN
158 : C
159 : C******* Message d'erreur standard
160 : C -301 0 %m1:40
161 : C
162 : MOTERR(1:40) = 'ENTI1 = ??? '
162 : CALL ERREUR(-301)
163 : C
164 : C******* Message d'erreur standard
165 : C Entier valant: %i1
166 : C -2 0
167 : C
168 : INTERR(1) = ORDESP
168 : CALL ERREUR(-2)
169 : C
170 : C******* Message d'erreur standard
171 : C 21 2
172 : C Données incompatibles
173 : C
173 : CALL ERREUR(21)
174 : GOTO 9999
175 : ENDIF
176 : C
177 : C**** Lecture ORDTEM
178 : C
179 : ICOND = 1
179 : CALL LIRENT(ORDTEM,ICOND,IRETOU)
180 : IF(IERR .NE. 0) GOTO 9999
181 : IF((ORDTEM .LT. 1) .OR. (ORDTEM .GT. ORDESP))THEN
182 : C
183 : C******* Message d'erreur standard
184 : C -301 0 %m1:40
185 : C
186 : MOTERR(1:40) = 'ENTI2 = ??? '
186 : CALL ERREUR(-301)
187 : C
188 : C******* Message d'erreur standard
189 : C Entier valant: %i1
190 : C -2 0
191 : C
192 : INTERR(1) = ORDTEM
192 : CALL ERREUR(-2)
193 : C
194 : C******* Message d'erreur standard
195 : C 21 2
196 : C Données incompatibles
197 : C
197 : CALL ERREUR(21)
198 : GOTO 9999
199 : ENDIF
200 : C
201 : IF((ORDESP .EQ. 1) .AND. (ORDTEM .EQ. 1))THEN
202 : C
203 : C******* ORDESP = 1, ORDTEM = 1
204 : C
204 : CALL PRE12F()
205 : ELSEIF((ORDESP .EQ. 2))THEN
206 : C
207 : C******* ORDESP = 2, ORDTEM = 1 ou 2
208 : C
208 : CALL PRE42F(ORDTEM)
209 : ELSE
210 : C
211 : C******* Message d'erreur standard
212 : C 5 3
213 : C Erreur anormale. Contactez votre support
214 : C
214 : CALL ERREUR(5)
215 : ENDIF
216 : C
217 : 9999 CONTINUE
218 : C
219 : RETURN
220 : END
221 :
222 :
223 :
224 :
225 :
© Cast3M 2003 - All rights reserved.
Disclaimer