* fichier : clmult2D.dgibi
*********************************************************************
* VF, Condition Limites en 2D pour les Equations d'Euler *
* multiespeces *
* Kudriakov S., DM2S/SFME/LTMF, FEVRIER 2003 *
* *
* Beccantini A., DM2S/SFME/LTMF, FEVRIER 2006: *
* BC of RESErvoir type *
*********************************************************************
TYEL = 'QUA8' ;
* TYEL = 'TRI6' ;
'OPTION' 'DIME' 2 'ELEM' TYEL 'ISOV' 'SULI'
'ECHO' 1 'TRAC' 'X' ;
GRAPH = VRAI ;
GRAPH = FAUX ;
******************
**** MAILLAGE ****
******************
*
* P2 P3
*
DX = 0.5 ;
P2 = 0.0 0.0 ;
P3 = DX DX ;
P2P3 = P2 'DROIT' 1 P3 ;
LIGHAU = (P2P3) 'PLUS' (4 '*' ((-1*DX) DX)) ;
DOMINT = P2P3 'REGLER' 4 LIGHAU ;
DOMINT = 'CHANGER' DOMINT 'QUAF' ;
*
**** Creation of DOMAINE tables via the MODEL object
*
MDOMINT = 'MODELISER' DOMINT 'NAVIER_STOKES' 'LINE' ;
MP2P3 = 'MODELISER' P2P3 'NAVIER_STOKES' 'LINE' ;
$DOMINT = 'DOMA' MDOMINT 'TABLE' ;
$P2P3 = 'DOMA' MP2P3 'TABLE' ;
DOMINT = 'DOMA' MDOMINT 'MAILLAGE' ;
P2P3 = 'DOMA' MP2P3 'MAILLAGE' ;
PCEN = ($DOMINT . 'CENTRE') 'POIN' 'PROC' ((DX '/' 2) (DX /2)) ;
MDICEN = $DOMINT . 'CENTRE' ;
MFCEN = $P2P3 . 'CENTRE' ;
'SI' GRAPH ;
'TRACER' (DOMINT 'ET' PCEN) ;
'FINSI' ;
**********************************************************************************
**********************************************************************************
*************************** Riemann invariants, inlet ***************************
**********************************************************************************
**********************************************************************************
NESP = 4 ;
*-----------------------------------------------------
*** GAS: H_2, O_2, H_2O, N_2
*
* CP, CV en J/Kg/K @ T = 3000
*-----------------------------------------------------
PGAS = 'TABLE' ;
PGAS . 'CP' = 'TABLE' ;
PGAS . 'CP' . 'H2 ' = .18729066D+05 ;
PGAS . 'CP' . 'O2 ' = .11886820D+04 ;
PGAS . 'CP' . 'H2O ' = .31209047D+04 ;
PGAS . 'CP' . 'N2 ' = .12993995D+04 ;
PGAS . 'CV' = 'TABLE' ;
PGAS . 'CV' . 'H2 ' = .14571861D+05 ;
PGAS . 'CV' . 'O2 ' = .92885670D+03 ;
PGAS . 'CV' . 'H2O ' = .26589930D+04 ;
PGAS . 'CV' . 'N2 ' = .10024563D+04;
*----------------------------------------------------
**** Especes qui sont dans les equations d'Euler
*----------------------------------------------------
PGAS . 'ESPEULE' = 'MOTS' 'H2 ' 'O2 ' 'H2O ' ;
*----------------------------------------------------
**** Espece qui n'y est pas
*----------------------------------------------------
PGAS . 'ESPNEULE' = 'N2 ' ;
*-----------------
*** CP ****
*-----------------
CP = ('MANUEL' 'CHPO' MDICEN 1 'H2'
(PGAS . 'CP' . 'H2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MDICEN 1 'O2'
(PGAS . 'CP' . 'O2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MDICEN 1 'H2O'
(PGAS . 'CP' . 'H2O ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MDICEN 1 'N2'
(PGAS . 'CP' . 'N2 ') 'NATU' 'DISCRET') ;
***************************
CPCEL = ('MANUEL' 'CHPO' MFCEN 1 'H2'
(PGAS . 'CP' . 'H2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MFCEN 1 'O2'
(PGAS . 'CP' . 'O2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MFCEN 1 'H2O'
(PGAS . 'CP' . 'H2O ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MFCEN 1 'N2'
(PGAS . 'CP' . 'N2 ') 'NATU' 'DISCRET') ;
*-----------------
*** CV ****
*-----------------
CV = ('MANUEL' 'CHPO' MDICEN 1 'H2'
(PGAS . 'CV' . 'H2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MDICEN 1 'O2'
(PGAS . 'CV' . 'O2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MDICEN 1 'H2O'
(PGAS . 'CV' . 'H2O ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MDICEN 1 'N2'
(PGAS . 'CV' . 'N2 ') 'NATU' 'DISCRET') ;
*****************************
CVCEL = ('MANUEL' 'CHPO' MFCEN 1 'H2'
(PGAS . 'CV' . 'H2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MFCEN 1 'O2'
(PGAS . 'CV' . 'O2 ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MFCEN 1 'H2O'
(PGAS . 'CV' . 'H2O ') 'NATU' 'DISCRET') 'ET'
('MANUEL' 'CHPO' MFCEN 1 'N2'
(PGAS . 'CV' . 'N2 ') 'NATU' 'DISCRET') ;
***********************************************************
***********************************************************
RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.110 ;
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.41 ;
VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 0.1 'UX' 1.65 ;
YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
'H2O' 0.1 ;
YN1 = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
'H2O' 0.1 'NATU' 'DISCRET' ;
**********************************************************
UNCH = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.0 ;
YNESP = UNCH '-' ('PSCAL' UNCH YN ('MOTS' 'SCAL' 'SCAL' 'SCAL')
(PGAS . 'ESPEULE')) ;
YNESP = 'NOMC' 'N2 ' YNESP 'NATU' 'DISCRET' ;
YNTOT = YN1 'ET' YNESP ;
LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
GAMN = CPTOT '/' CVTOT ;
**********************************************************
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' 2.11 'PN' 111.4
'UX' 2.0 'UY' 1.0 'H2' 0.1 'O2' 0.7 'H2O' 0.05 ;
LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INRI' ;
************************************************************
**** On calcule les variables conservatives **************
************************************************************
RAZM = NESP '+' 3 ;
GN = RN '*' VN ('MOTS' 'SCAL' 'SCAL')
('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY') ;
LMOT1 = 'MOTS' 'UX' 'UY' ;
VVV = 'PSCAL' VN VN LMOT1 LMOT1 ;
ECIN = 0.5D0 '*' RN '*' VVV ;
RETN = 'NOMC' 'SCAL' ((PN '/' (GAMN '-' 1.0)) '+' ECIN)
'NATU' 'DISCRET' ;
LMOT1 = (PGAS . 'ESPEULE') ;
LMOT2 = 'MOTS' 'SCAL' 'SCAL' 'SCAL' ;
LMOT3 = (PGAS . 'ESPEULE') ;
RYN = YN '*' RN LMOT1 LMOT2 LMOT3 ;
VIT PRES TTEMP FRAM GAMM = 'PRIM' 'PERFMULT' PGAS RN GN RETN RYN ;
ERR1 = 'MAXIMUM' (PRES '-' PN) 'ABS' ;
ERR2 = 'MAXIMUM' (GN '-' (RN '*' VIT ('MOTS' 'SCAL' 'SCAL' )
('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) 'ABS' ;
'SI' (ERR1 > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI' ;
'SI' (ERR2 > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI' ;
*****************************************
**** On calcule RCHLIM en gibiane
*****************************************
UNC = (('EXCO' VN 'UY' 'SCAL') '-' ('EXCO' VN 'UX' 'SCAL')) '/' (2
'**' 0.5) ;
GAMMA = 'MAXIMUM' GAMN ;
GM1 = '-' GAMMA 1. ;
USGM1 = 1. '/' GM1 ;
DSGM1 = 2. '*' USGM1 ;
ASONC = GAMMA * (PN '/' RN) ;
ASONC = ASONC '**' 0.5 ;
G1 = UNC '-' (DSGM1 * ASONC) ;
G1 = 'EXTRAIRE' G1 PCEN 'SCAL' ;
*******************************************************
RNCEL= 'EXCO' CHPLI 'RN' ;
PNCEL = 'EXCO' CHPLI 'PN' ;
UNCELC = (('EXCO' CHPLI 'UY' 'SCAL') '-' ('EXCO' CHPLI 'UX' 'SCAL'))
'/' (2 '**' 0.5) ;
UTCELC = (('EXCO' CHPLI 'UY' 'SCAL') '+' ('EXCO' CHPLI 'UX' 'SCAL'))
'/' (-1*(2 '**' 0.5)) ;
YH2CEL = 'EXCO' CHPLI 'H2' 'NATU' 'DISCRET' ;
YO2CEL = 'EXCO' CHPLI 'O2' 'NATU' 'DISCRET' ;
YHOCEL = 'EXCO' CHPLI 'H2O' 'NATU' 'DISCRET' ;
YNCEL = YH2CEL 'ET' YO2CEL ;
YNCEL = YNCEL 'ET' YHOCEL ;
**********************************************************
YNTCEL = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 4 'H2' 0.1 'O2' 0.7
'H2O' 0.05 'N2' 0.15 'NATU' 'DISCRET' ;
LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
CVTOT = 'PSCAL' YNTCEL CVCEL LMOT1 LMOT1 ;
CPTOT = 'PSCAL' YNTCEL CPCEL LMOT1 LMOT1 ;
GAMCEL = CPTOT '/' CVTOT ;
***********************************************************
GAMF = 'MAXIMUM' GAMCEL ;
GM1F = '-' GAMF 1. ;
USGM1F = 1. '/' GM1F ;
DSGM1F = 2. '*' USGM1F ;
ASONF = GAMF * (PNCEL '/' RNCEL) ;
ASONF = ASONF '**' 0.5 ;
G3 = UNCELC '+' (DSGM1F * ASONF) ;
G3 = 'EXTRAIRE' G3 ('POIN' 1 ($P2P3 . 'CENTRE')) 'SCAL' ;
S = PNCEL '/' (RNCEL** GAMF) ;
S = 'EXTRAIRE' S ('POIN' 1 ($P2P3 . 'CENTRE')) 'SCAL' ;
GAMID = (GAMMA '+' GAMF) '/' 2.0D0 ;
GM1M = '-' GAMID 1. ;
USGM1M = 1. '/' GM1M ;
DSGM1M = 2. '*' USGM1M ;
ASON2=(0.5D0*(G3-G1)) ;
ASON2=ASON2/DSGM1F ;
ASON2=ASON2*ASON2 ;
RHO=ASON2/(GAMF*S) ;
RHO=RHO**USGM1F ;
P=(RHO*ASON2)/GAMF ;
UNCEL = 0.5 '*' (G3 '+' G1) ;
UTCEL = 'EXTRAIRE' UTCELC ('POIN' 1 ($P2P3 . 'CENTRE')) 'SCAL' ;
UX = (UNCEL '+' UTCEL) '/' (-1*(2 '**' 0.5)) ;
UY = (UNCEL '-' UTCEL) '/' (2 '**' 0.5) ;
CHLIM2 = 'MANU' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' RHO 'UX'
UX 'UY' UY 'PN' P 'H2' 0.1 'O2' 0.7 'H2O' 0.05 ;
ERRO = 'MAXIMUM' (CHLIM2 '-' RCHLIM) ;
'SI' (ERRO > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI' ;
'FINSI' ;
******************************************************************
**** Jacobian with respect to primitive variable
******************************************************************
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INRI' ;
EPS=1.0D-6 ;
ERROM = 1.0D-3 ;
* RN ****************** RN ******* RN *********************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
(RN '+' DR) VN PN YN CHPLI 'INRI' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
* UX ***** UX ****** UX ****** UX ************************
DUXSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DUX) PN YN CHPLI 'INRI' ;
DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0) ;
MMM = 'MAXIMUM' DRESDUX 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDUX '-' DRESDUXN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*** UY *************** UY ************* UY *******************
DUYSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DUY) PN YN CHPLI 'INRI' ;
DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDUY 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDUY '-' DRESDUYN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
**** PN ********** PN ************ PN ***********************
DPSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DP = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DPSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN (PN '+' DP) YN CHPLI 'INRI' ;
DRESDPN = ( RCHRES1 '-' RCHRES) '/' DPSCAL ;
DRESDP = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDP 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDP '-' DRESDPN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
**** YH2 ********** YH2 ************ YH2 ***********************
DYH2SCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DYH2 = ('MANUEL' 'CHPO' ELE1 1 'H2' DYH2SCAL ) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN1 '+' DYH2) CHPLI 'INRI' ;
DRESDYN = ( RCHRES1 '-' RCHRES) '/' DYH2SCAL ;
DRESDY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDY 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDY '-' DRESDYN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***** YO2 ****** YO2 *********** YO2 ******************
DYO2SCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DYO2 = ('MANUEL' 'CHPO' ELE1 1 'O2' DYO2SCAL ) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN1 '+' DYO2) CHPLI 'INRI' ;
DRESDYN = ( RCHRES1 '-' RCHRES) '/' DYO2SCAL ;
DRESDY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDY 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDY '-' DRESDYN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***** YHO ****** YHO *********** YHO ******************
DYHOSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DYHO = ('MANUEL' 'CHPO' ELE1 1 'H2O' DYHOSCAL ) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN1 '+' DYHO) CHPLI 'INRI' ;
DRESDYN = ( RCHRES1 '-' RCHRES) '/' DYHOSCAL ;
DRESDY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
MMM = 'MAXIMUM' DRESDY 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDY '-' DRESDYN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
**************************************************************
**** Jacobian with respect to conservative variables
**************************************************************
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INRI' ;
******** RN ************* RN ********** RN **************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = RN '+' DR ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***** RUX ******* RUX ******** RUX *******************
DUXSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
DUX = ('MANUEL' 'CHPO' ELE1 1 'UX' DUXSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DUX ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDUXN = ( RCHRES1 '-' RCHRES) '/' DUXSCAL ;
DRESDUX = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0
'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDUX 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDUX '-' DRESDUXN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***** RUY ******* RUY ******** RUY *******************
DUYSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
DUY = ('MANUEL' 'CHPO' ELE1 1 'UY' DUYSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DUY ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDUYN = ( RCHRES1 '-' RCHRES) '/' DUYSCAL ;
DRESDUY = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0
'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDUY 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDUY '-' DRESDUYN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
******** RETN ************ RETN ************************
DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN '+' DR ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0
'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
******** RYN ************ RYN ************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
******** RYN ************ RYN ************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
******** RYNH2 ************ RYNH2 ************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
******** RYNO2 ************ RYNO2 ************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
******** RYNH20 ************ RYNH20 ************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INRI' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0
'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
******************************************************************************
******************************************************************************
*************************** INJE ********************************************
******************************************************************************
******************************************************************************
*
* We impose RT and the normal momentum at the inlet.
* (les fractions massiques sont impsees aussi)
***** Calcul du residu dans le cas particulier UF = UC
*****************************************************************************
RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 15. 'UX' -15. ;
YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
'H2O' 0.1 ;
GN = RN '*' VN ;
MOME = (15 '*' (2. '**' 0.5)) * 1.1 ;
RT = 2.1 '/' 1.1 ;
UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
PC = ('EXTRAIRE' PN 'SCAL' PCEN) ;
RC = ('EXTRAIRE' RN 'SCAL' PCEN) ;
YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'RT' RT
'MOME' MOME 'H2' 0.3 'O2' 0.5 'H2O' 0.1 ;
LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INJE' ;
*---------------------------------------
** Testing RCHLIM !!!!!
*---------------------------------------
RCHLIM1 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' RC 'PN' PC 'UX' UXC
'UY' UYC 'H2' 0.3 'O2' 0.5 'H2O' 0.1 ;
ERRO = 'MAXIMUM' (RCHLIM1 '-' RCHLIM) 'ABS' ;
'SI' (ERRO > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI' ;
*-------------------------------------
** Testing RCHRES !!!!!
*-------------------------------------
'SI' ('EGA' TYEL 'QUA8') ;
SURSVOL = 1. '/' ((2 '**' 0.5) * 0.5) ;
'SINON' ;
SURSVOL = 2. '/' ((2 '**' 0.5) * 0.5) ;
'FINSI' ;
*-------------------------------------------------------------
YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
'H2O' YH2O 'N2' YN2 ;
LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
GAMN = CPTOT '/' CVTOT ;
GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
*--------------------------------------------------------------
FMASS = MOME ;
FUX = (MOME * UXC) '-' (PC '/' (2. '**' 0.5)) ;
FUY = (MOME * UYC) '+' (PC '/' (2. '**' 0.5)) ;
HIN = (GAMC/(GAMC-1)) '*' (PC '/' RC) ;
ECIN = 0.5 '*' ((UXC '*' UXC) '+' (UYC '*' UYC)) ;
FEN = MOME '*' (HIN '+' ECIN)
;
RCHRES1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' PCEN) 7
'VAR1' (SURSVOL '*' FMASS) 'VAR2' (FUX '*' SURSVOL) 'VAR3' (FUY '*'
SURSVOL) 'VAR4' (FEN '*' SURSVOL) 'VAR5' (FMASS '*' YH2 '*' SURSVOL)
'VAR6' (FMASS '*' YO2'*'SURSVOL) 'VAR7' (FMASS '*' YH2O'*'SURSVOL);
ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
'SI' (ERRO > 1.0D-10) ;
'ERREUR' 5 ;
'FINSI' ;
****************************************************************
************* JACOBIANS ***************************************
****************************************************************
RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 1.11 'UX' -1.3 ;
YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
'H2O' 0.1 ;
YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
'H2O' YH2O 'N2' YN2 ;
LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
GAMN = CPTOT '/' CVTOT ;
GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
LMOT1 = (PGAS . 'ESPEULE') ;
LMOT2 = 'MOTS' 'SCAL' 'SCAL' 'SCAL' ;
LMOT3 = (PGAS . 'ESPEULE') ;
RYN = YN '*' RN LMOT1 LMOT2 LMOT3 ;
************************************************
GN = RN '*' VN ;
RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
MOME = 0.111 ;
RT = 4.1 ;
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'RT' RT
'MOME' MOME 'H2' 0.1 'O2' 0.5 'H2O' 0.3 ;
LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INJE' ;
***************************************************************
***** Jacobian with respect to the primitive variables ********
***************************************************************
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INJE' ;
* EPS=1.0D-6 ;
* ERROM = 1.0D-5 ;
************ RN *************** RN *****************************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
(RN '+' DR) VN PN YN CHPLI 'INJE' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************ UX *************** UX *****************************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'INJE' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************ UY *************** UY *****************************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'INJE' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************ PN *************** PN *****************************
DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN (PN '+' DR) YN CHPLI 'INJE' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************ H2 *************** H2 *****************************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INJE' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************ O2 *************** O2 *****************************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INJE' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************ H2O *************** H2O *****************************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INJE' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************************************************************
**** Jacobian with respect to conservative variables ******
************************************************************
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INJE' ;
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INJE' ;
***************** RN ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = RN '+' DR ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RUX ************************************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RUY ************************************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RETN ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN '+' DR ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYH2 ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYO2 ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYH2O ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INJE' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************************************************************************
************************************************************************
*************************** Pressure at the outlet ********************
************************************************************************
************************************************************************
*
* We solve a Riemann probleme between the state in the center and
* the same state but the pressure equal to the pressure at the outlet
* We use AUSMPLUS
*
*
***** Calcul du residu dans le cas particulier PF = PC
RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' -15. 'UX' 1. ;
YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
'H2O' 0.1 ;
GN = RN '*' VN ;
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 1 'PN' 2.1 ;
*----------------------------------------------------------
YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
'H2O' YH2O 'N2' YN2 ;
LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
GAMN = CPTOT '/' CVTOT ;
GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
*----------------------------------------------------------
LMOT1 = (PGAS . 'ESPEULE') ;
LMOT2 = 'MOTS' 'SCAL' 'SCAL' 'SCAL' ;
LMOT3 = (PGAS . 'ESPEULE') ;
RYN = YN '*' RN LMOT1 LMOT2 LMOT3 ;
************************************************
RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'OUTP' ;
*-----------------------------------------------
'SI' ('EGA' TYEL 'QUA8') ;
SURSVOL = -1. '/' ((2 '**' 0.5) * 0.5) ;
'SINON' ;
SURSVOL = -2. '/' ((2 '**' 0.5) * 0.5) ;
'FINSI' ;
UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
UNORM = (UXC '-' UYC) '/' (2. '**' 0.5) ;
UT = (UXC '+' UYC) '/' (2. '**' 0.5) ;
PC = 'EXTRAIRE' PN PCEN 'SCAL' ;
RC = 'EXTRAIRE' RN PCEN 'SCAL' ;
FMASS = RC '*' UNORM ;
FUN = (RC '*' UNORM '*' UNORM) '+' PC ;
FUT = (RC '*' UNORM '*' UT) ;
FUX = (FUN '+' FUT) '/' (2. '**' 0.5) ;
FUY = (FUT '-' FUN) '/' (2. '**' 0.5) ;
FEN = (UNORM '*' ((GAMC/(GAMC-1)) '*' PC)) '+'
(0.5 '*' UNORM '*' RC '*' ((UNORM '*' UNORM) '+' (UT '*' UT)))
;
FYH2 = YH2 '*' RC '*' UNORM ;
FYO2 = YO2 '*' RC '*' UNORM ;
FYH2O = YH2O '*' RC '*' UNORM ;
RCHRES1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' PCEN) 7
'VAR1' (SURSVOL '*' FMASS) 'VAR2' (FUX '*' SURSVOL) 'VAR3' (FUY '*'
SURSVOL) 'VAR4' (FEN '*' SURSVOL) 'VAR5' (FYH2 '*' SURSVOL)
'VAR6' (FYO2 '*' SURSVOL) 'VAR7' (FYH2O '*' SURSVOL) ;
ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
'SI' (ERRO > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI';
*------------------------------------------------------------
***** Jacobian with respect to the primitive variables ******
*------------------------------------------------------------
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 1 'PN' 4.1 ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'OUTP' ;
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'OUTP' ;
************* RN *********** RN ***********************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
(RN '+' DR) VN PN YN CHPLI 'OUTP' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* UX *********** UX ***********************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'OUTP' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* UY *********** UY ***********************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'OUTP' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* PN *********** PN ***********************
DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN (PN '+' DR) YN CHPLI 'OUTP' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* H2 *********** H2 ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'OUTP' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* O2 *********** O2 ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'OUTP' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* H2O *********** H2O ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'OUTP' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------------------------------------
**** Jacobian with respect to conservative variables
*---------------------------------------------------------
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'OUTP' ;
*************** RN **************** RN *******************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = RN '+' DR ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN)'/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RUX ************************************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RUY ************************************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RETN ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN '+' DR ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYH2 ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYO2 ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYH2O ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'OUTP' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*******************************************************************
*******************************************************************
*************************** INSU *********************************
*******************************************************************
*******************************************************************
* Total enthalpy, entropy, tangential velocity at the inlet
* Normal velocity inside
*
* We solve a Riemann probleme between the state at the interface thus
* reconstructed and the state at the center
*
***** Calcul du residu dans le cas particulier UF = UC
*
RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 15. 'UX' -15. ;
YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
'H2O' 0.1 ;
GN = RN '*' VN ;
*----------------------------------------------------------
YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
'H2O' YH2O 'N2' YN2 ;
LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
GAMN = CPTOT '/' CVTOT ;
GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
*----------------------------------------------------------
RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
RHTN = RETN '+' PN ;
HTN = RHTN '/' RN ;
UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
UTC = -1 '*' (UXC '+' UYC) '/' (2. '**' 0.5) ;
SN = PN '/' (RN '**' ('MAXI' GAMN)) ;
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'HT' ('MAXIMUM' HTN)
'S' ('MAXIMUM' SN) 'H2' 0.3 'O2' 0.5 'H2O' 0.1;
LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INSU' ;
*-------------------------------------------------------------
'SI' ('EGA' TYEL 'QUA8') ;
SURSVOL = -1. '/' ((2 '**' 0.5) * 0.5) ;
'SINON' ;
SURSVOL = -2. '/' ((2 '**' 0.5) * 0.5) ;
'FINSI' ;
UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
UNORM = (UXC '-' UYC) '/' (2. '**' 0.5) ;
UT = (UXC '+' UYC) '/' (2. '**' 0.5) ;
PC = 'EXTRAIRE' PN PCEN 'SCAL' ;
RC = 'EXTRAIRE' RN PCEN 'SCAL' ;
GAMC = 'EXTRAIRE' GAMN PCEN 'SCAL' ;
RCHLIM2 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 7 'RN' RC 'PN' PC 'UX'
UXC 'UY' UYC 'H2' YH2 'O2' YO2 'H2O' YH2O ;
ERRO = 'MAXIMUM' (RCHLIM '-' RCHLIM2) 'ABS' ;
'SI' (ERRO > 1.0D-6) ;
'ERREUR' 5 ;
'FINSI' ;
FMASS = RC '*' UNORM ;
FUN = (RC '*' UNORM '*' UNORM) '+' PC ;
FUT = (RC '*' UNORM '*' UT) ;
FUX = (FUN '+' FUT) '/' (2. '**' 0.5) ;
FUY = (FUT '-' FUN) '/' (2. '**' 0.5) ;
FEN = (UNORM '*' ((GAMC/(GAMC-1)) '*' PC)) '+'
(0.5 '*' UNORM '*' RC '*' ((UNORM '*' UNORM) '+'
(UT '*' UT))) ;
FH2 = RC '*' UNORM '*' YH2 ;
FO2 = RC '*' UNORM '*' YO2 ;
FHO = RC '*' UNORM '*' YH2O ;
RCHRES1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' PCEN) 7
'VAR1' (SURSVOL '*' FMASS) 'VAR2' (FUX '*' SURSVOL)
'VAR3' (FUY '*' SURSVOL) 'VAR4' (FEN '*' SURSVOL)
'VAR5' (FH2 '*' SURSVOL) 'VAR6' (FO2 '*' SURSVOL)
'VAR7' (FHO '*' SURSVOL) ;
ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
'SI' (ERRO > 1.0D-6) ;
'ERREUR' 5 ;
'FINSI' ;
*------------------------------------------------------------
***** Jacobian with respect to the primitive variables ******
*------------------------------------------------------------
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'HT' 250.0
'S' 2.0 'H2' 0.1 'O2' 0.7 'H2O' 0.2;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INSU' ;
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INSU' ;
************* RN *********** RN ***********************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
(RN '+' DR) VN PN YN CHPLI 'INSU' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* UX *********** UX ***********************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'INSU' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* UY *********** UY ***********************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'INSU' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* PN *********** PN ***********************
DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN (PN '+' DR) YN CHPLI 'INSU' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* H2 *********** H2 ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INSU' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* O2 *********** O2 ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INSU' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* H2O *********** H2O ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INSU' ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------------------------------------
**** Jacobian with respect to conservative variables
*---------------------------------------------------------
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INSU' ;
*************** RN **************** RN *******************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = RN '+' DR ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
NNN = 'MAXIMUM' RCHRES 'ABS' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
* 'OPTION' DONN 5 ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RUX ************************************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
* 'OPTION' DONN 5 ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RUY ************************************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RETN ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = RETN '+' DR ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYH2 ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYO2 ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
***************** RYH2O ************************************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INSU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*************************************************************
*************************************************************
*------- 'INOU' -------*
*************************************************************
*************************************************************
RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 1.1 ;
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' 2.1 ;
VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' 15. 'UX' -15. ;
VN1 = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' -15. 'UX' 25. ;
VN2 = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UY' -15. 'UX' 25. ;
YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 3 'H2' 0.3 'O2' 0.5
'H2O' 0.1 ;
GN = RN '*' VN ;
GN2 = RN '*' VN1 ;
*---------------------------------------------------------
YH2 = ('EXTRAIRE' YN 'H2' PCEN) ;
YO2 = ('EXTRAIRE' YN 'O2' PCEN) ;
YH2O = ('EXTRAIRE' YN 'H2O' PCEN) ;
YN2 = 1.0D0 '-' (YH2 '+' YO2 '+' YH2O) ;
YNTOT = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 4 'H2' YH2 'O2' YO2
'H2O' YH2O 'N2' YN2 ;
LMOT1 = (PGAS . 'ESPEULE') 'ET' ('MOTS' (PGAS . 'ESPNEULE')) ;
CVTOT = 'PSCAL' YNTOT CV LMOT1 LMOT1 ;
CPTOT = 'PSCAL' YNTOT CP LMOT1 LMOT1 ;
GAMN = CPTOT '/' CVTOT ;
GAMC = ('EXTRAIRE' GAMN 'SCAL' PCEN) ;
*----------------------------------------------------------
RETN = (('INVERSE' (GAMN '-' 1)) '*' PN) '+' (0.5 '*' ('PSCAL' VN GN
('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
RETN2 =(('INVERSE' (GAMN '-' 1)) '*' PN) '+'
(0.5 '*' ('PSCAL' VN1 GN2
('MOTS' 'UX' 'UY') ('MOTS' 'UX' 'UY'))) ;
RHTN = RETN '+' PN ;
HTN = RHTN '/' RN ;
UXC = ('EXTRAIRE' VN 'UX' PCEN) ;
UYC = ('EXTRAIRE' VN 'UY' PCEN) ;
UTC = -1 '*' (UXC '+' UYC) '/' (2. '**' 0.5) ;
SN = PN '/' (RN '**' ('MAXI' GAMN)) ;
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 6 'HT' ('MAXIMUM' HTN)
'S' ('MAXIMUM' SN) 'H2' 0.3 'O2' 0.5 'H2O' 0.1 'PN' 50.0 ;
CHPLI1 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 5 'HT' ('MAXIMUM' HTN)
'S' ('MAXIMUM' SN) 'H2' 0.3 'O2' 0.5 'H2O' 0.1 ;
CHPLI2 = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 1 'PN' 50.0 ;
LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' 'VAR6' 'VAR7' ;
LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'H2 ' 'O2 ' 'H2O ' ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INOU' ;
RCHLIM1 RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI1 'INSU' ;
ERRO = 'MAXIMUM' (RCHRES1 '-' RCHRES) 'ABS' ;
ERRO1 = 'MAXIMUM' (RCHLIM1 '-' RCHLIM) 'ABS' ;
* 'OPTION' DONN 5 ;
'SI' (ERRO > 1.0D-9) ;
'ERREUR' 5 ;
'FINSI';
'SI' (ERRO1 > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI';
*--------------------------------------------------------------
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN YN CHPLI 'INOU' ;
RCHLIM2 RCHRES2 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN YN CHPLI2 'OUTP' ;
ERRO = 'MAXIMUM' (RCHRES2 '-' RCHRES) 'ABS' ;
ERRO1 = 'MAXIMUM' (RCHLIM2 '-' RCHLIM) 'ABS' ;
'SI' (ERRO > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI';
'SI' (ERRO1 > 1.0D-12) ;
'ERREUR' 5 ;
'FINSI';
*------------------------------------------------------------
***** Jacobian with respect to the primitive variables ******
*------------------------------------------------------------
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 6 'HT' 250.0
'S' 2.0 'H2' 0.1 'O2' 0.7 'H2O' 0.2 'PN' 50.0 ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INOU' ;
RCHLIM2 RCHRES2 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN YN CHPLI 'INOU' ;
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INOU' ;
RJACO2 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOPRIM' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN YN CHPLI 'INOU' ;
************* RN *********** RN ***********************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
(RN '+' DR) VN PN YN CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
(RN '+' DR) VN1 PN YN CHPLI 'INOU' ;
DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 1.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* UX *********** UX ***********************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
*-------------------------------------------
DRSCAL2 = EPS '*' ('EXTRAIRE' VN2 'UX' PCEN) ;
DR2 = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL2) ;
*-------------------------------------------
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN1 '+' DR2) PN YN CHPLI 'INOU' ;
DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL2 ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 1.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* UY *********** UY ***********************
DRSCAL = EPS '*' ('EXTRAIRE' VN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
*-------------------------------------------
DRSCAL2 = EPS '*' ('EXTRAIRE' VN2 'UY' PCEN) ;
DR2 = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL2) ;
*-------------------------------------------
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN '+' DR) PN YN CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN (VN1 '+' DR2) PN YN CHPLI 'INOU' ;
DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL2 ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 1.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* PN *********** PN ***********************
DRSCAL = EPS '*' ('EXTRAIRE' PN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN (PN '+' DR) YN CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 (PN '+' DR) YN CHPLI 'INOU' ;
DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 1.0 'H2' 0.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* H2 *********** H2 ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN (YN '+' DR) CHPLI 'INOU' ;
DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 1.0 'O2' 0.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* O2 *********** O2 ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN (YN '+' DR) CHPLI 'INOU' ;
DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 1.0 'H2O' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
************* H20 *********** H20 ***********************
DRSCAL = EPS '*' ('EXTRAIRE' YN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN (YN '+' DR) CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN (YN '+' DR) CHPLI 'INOU' ;
DRESDRN2 = ( RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
DRESDRN = ( RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'RN' 0.0
'UX' 0.0 'UY' 0.0 'PN' 0.0 'H2' 0.0 'O2' 0.0 'H2O' 1.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------------------------------------
**** Jacobian with respect to conservative variables
*---------------------------------------------------------
RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'INOU' ;
RJACO2 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'JACOCONS' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN1 PN YN CHPLI 'INOU' ;
*************** RN **************** RN *******************
DRSCAL = EPS '*' ('EXTRAIRE' RN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
RN1 = RN '+' DR ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
VN2 PN2 TTEMP YN2 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN2 RETN2 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN2 PN2 YN2 CHPLI 'INOU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
DRESDRN2 = (RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 1.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*************** RUX **************** RUX *******************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL) ;
*---------------------------------------------
DRSCAL2 = EPS '*' ('EXTRAIRE' GN2 'UX' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR2 = ('MANUEL' 'CHPO' ELE1 1 'UX' DRSCAL2) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
GN3 = GN2 '+' DR2 ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
VN2 PN2 TTEMP YN2 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN3 RETN2 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN2 PN2 YN2 CHPLI 'INOU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
DRESDRN2 = (RCHRES3 '-' RCHRES2) '/' DRSCAL2 ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 1.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*************** RUY **************** RUY *******************
DRSCAL = EPS '*' ('EXTRAIRE' GN 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL) ;
*---------------------------------------------
DRSCAL2 = EPS '*' ('EXTRAIRE' GN2 'UY' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR2 = ('MANUEL' 'CHPO' ELE1 1 'UY' DRSCAL2) ;
RN1 = 'COPIER' RN ;
GN1 = GN '+' DR ;
GN3 = GN2 '+' DR2 ;
RETN1 = 'COPIER' RETN ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
VN2 PN2 TTEMP YN2 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN3 RETN2 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN2 PN2 YN2 CHPLI 'INOU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
DRESDRN2 = (RCHRES3 '-' RCHRES2) '/' DRSCAL2 ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 1.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*************** RETN **************** RETN *******************
DRSCAL = EPS '*' ('EXTRAIRE' RETN 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL) ;
*---------------------------------------------
DRSCAL2 = EPS '*' ('EXTRAIRE' RETN2 'SCAL' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR2 = ('MANUEL' 'CHPO' ELE1 1 'SCAL' DRSCAL2) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = (RETN '+' DR) ;
RETN3 = (RETN2 '+' DR2) ;
RYN1 = 'COPIER' RYN ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
VN2 PN2 TTEMP YN2 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN2 RETN3 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN2 PN2 YN2 CHPLI 'INOU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
DRESDRN2 = (RCHRES3 '-' RCHRES2) '/' DRSCAL2 ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 1.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
*----------------------------
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*************** RH2 **************** RH2 *******************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
VN2 PN2 TTEMP YN2 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN2 RETN2 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN2 PN2 YN2 CHPLI 'INOU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
DRESDRN2 = (RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 1.0 'VAR6' 0.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
*----------------------------
* 'OPTION' DONN 5 ;
* 'SI' (ERRO > ERROM) ;
* 'ERREUR' 5 ;
* 'FINSI' ;
*----------------------------
* 'SI' (ERRO2 > ERROM) ;
* 'ERREUR' 5 ;
* 'FINSI' ;
*************** RO2 **************** RO2 *******************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'O2' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'O2' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
VN2 PN2 TTEMP YN2 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN2 RETN2 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN2 PN2 YN2 CHPLI 'INOU' ;
NNN = 'MAXIMUM' RCHRES 'ABS' ;
VERH = (RCHRES1 '-' RCHRES) '/' NNN ;
DRESDRN = VERH '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
DRESDRN2 = (RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 1.0 'VAR7' 0.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
* 'OPTION' DONN 5 ;
*----------------------------
'SI' (ERRO > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*************** RH2O **************** RH2O *******************
DRSCAL = EPS '*' ('EXTRAIRE' RYN 'H2O' PCEN) ;
ELE1 = 'MANUEL' 'POI1' PCEN ;
DR = ('MANUEL' 'CHPO' ELE1 1 'H2O' DRSCAL) ;
RN1 = 'COPIER' RN ;
GN1 = 'COPIER' GN ;
RETN1 = 'COPIER' RETN ;
RYN1 = RYN '+' DR ;
VN1 PN1 TTEMP YN1 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN1 RETN1 RYN1 ;
VN2 PN2 TTEMP YN2 GAMM =
'PRIM' 'PERFMULT' PGAS RN1 GN2 RETN2 RYN1 ;
RCHLIM RCHRES1 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN1 PN1 YN1 CHPLI 'INOU' ;
RCHLIM RCHRES3 = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN1 VN2 PN2 YN2 CHPLI 'INOU' ;
DRESDRN = (RCHRES1 '-' RCHRES) '/' DRSCAL ;
DRESDR = 'KOPS' RJACO '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
DRESDRN2 = (RCHRES3 '-' RCHRES2) '/' DRSCAL ;
DRESDR2 = 'KOPS' RJACO2 '*' ('MANUEL' 'CHPO' ELE1 7 'VAR1' 0.0
'VAR2' 0.0 'VAR3' 0.0 'VAR4' 0.0 'VAR5' 0.0 'VAR6' 0.0 'VAR7' 1.0 ) ;
MMM = 'MAXIMUM' DRESDR 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO = 'MAXIMUM' (DRESDR '-' DRESDRN) 'ABS' ;
'SINON' ;
ERRO = 'MAXIMUM' ((DRESDR '-' DRESDRN) '/' MMM) 'ABS' ;
'FINSI' ;
*---------------------------
MMM = 'MAXIMUM' DRESDR2 'ABS' ;
'SI' (MMM < ERROM) ;
ERRO2 = 'MAXIMUM' (DRESDR2 '-' DRESDRN2) 'ABS' ;
'SINON' ;
ERRO2 = 'MAXIMUM' ((DRESDR2 '-' DRESDRN2) '/' MMM) 'ABS' ;
'FINSI' ;
*----------------------------
* 'OPTION' DONN 5 ;
* 'SI' (ERRO > ERROM) ;
* 'ERREUR' 5 ;
* 'FINSI' ;
*----------------------------
'SI' (ERRO2 > ERROM) ;
'ERREUR' 5 ;
'FINSI' ;
*
* BC of 'RESE'rvoir type
*
SURSVOL = ('MAXIMUM' ('DOMA' MP2P3 'VOLUME')) '/'
('MAXIMUM' ('DOMA' MDOMINT . 'VOLUME')) ;
*
* We give the pressure of the reservoir, the density, the homogeneous
* composition. We suppose that the tangential velocity is zero.
* We compute the inlet boundary conditions by supposing that the inlet
* is at the throat, namely the inlet pressure is equal to the maximum
* between the critical and the inlet pressure.
* The inlet density and velocity are found using the Riemann invariant.
*
* Two species
* HELIUM
* AIR
*
PGAS = 'TABLE' ;
PGAS . 'CP' = 'TABLE' ;
PGAS . 'CP' . 'AIR ' = 1.00D+03 ;
PGAS . 'CP' . 'HE ' = 5.18D+03 ;
PGAS . 'CV' = 'TABLE' ;
PGAS . 'CV' . 'AIR ' = .717D+03 ;
PGAS . 'CV' . 'HE ' = 3.10D+03 ;
*----------------------------------------------------
**** Especes qui sont dans les equations d'Euler
*----------------------------------------------------
PGAS . 'ESPEULE' = 'MOTS' 'AIR ' ;
*----------------------------------------------------
**** Espece qui n'y est pas
*----------------------------------------------------
PGAS . 'ESPNEULE' = 'HE ' ;
*
* Reservoir conditions
*
P0 = 200.0D5 ;
RHO0 = 32.1D0 ;
Y0 = 1.0D-5 ;
gamf = ((PGAS . 'CP' . 'HE ') '*' (1. '-' Y0)) '+'
((PGAS . 'CP' . 'AIR ') '*' (Y0)) ;
gamf = gamf '/'
(((PGAS . 'CV' . 'HE ') '*' (1. '-' Y0)) '+'
((PGAS . 'CV' . 'AIR ') '*' (Y0))) ;
PCR = (2. '/' (gamf '+' 1.)) ;
PCR = PCR '**' (gamf '/' (gamf '-' 1.)) ;
PCR = PCR '*' P0 ;
RCR = (PCR '/' P0) '**' (1. '/' gamf) ;
RCR = RCR '*' RHO0 ;
UCR = (2 '*' gamf '/' (gamf '-' 1.0D0)) '*' ((P0 '/' RHO0) '-'
(PCR '/' RCR)) ;
UCR = UCR '**' 0.5 ;
'SI' ('NEG' (ucr * ucr) (gamf * pcr '/' rcr) (1.0D-12 *ucr)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
* Case 1: internal pressure lower than critical pressure
*
P1 = 1.0D5 ;
RHO1 = 0.954D0 ;
Y1 = 0.99999 ;
*
YN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'AIR' Y1 ;
RN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' RHO1 ;
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' P1 ;
VN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 2 'UX' 0.0 'UY' 0.0 ;
CHPLI = 'MANUEL' 'CHPO' ($P2P3 . 'CENTRE') 3 'RN' RHO0 'PN' P0
'AIR ' y0 ;
LISTC = 'MOTS' 'VAR1' 'VAR2' 'VAR3' 'VAR4' 'VAR5' ;
LISTP = 'MOTS' 'RN ' 'UX ' 'UY ' 'PN ' 'AIR ' ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'RESE' ;
aa = 'MAXIMUM' ('EXCO' RCHLIM 'RN') ;
'SI' ('NEG' (rcr) (aa) (1.0D-12 * aa)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
bb = 'MAXIMUM' ('EXCO' RCHLIM 'PN') ;
'SI' ('NEG' (pcr) (bb) (1.0D-12 * bb)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
cc = 'MAXIMUM' ('EXCO' RCHLIM 'UX') ;
dd = 'MAXIMUM' ('EXCO' RCHLIM 'UY') ;
ee = (cc * cc) '+' (dd * dd) ;
ee = ee '**' 0.5 ;
'SI' ('NEG' (ee) (ucr) (1.0D-12 * ee)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
aa = 'MAXIMUM' ('EXCO' RCHRES 'VAR1') ;
bb = SURSVOL '*' rcr * ucr ;
'SI' ('NEG' (bb) (aa) (1.0D-12 * aa)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
htot = (gamf '/' (gamf '-' 1.0)) '*' p0 '/' rho0 ;
aa = htot * rcr * ucr * sursvol ;
bb = 'MAXIMUM' ('EXCO' RCHRES 'VAR4') ;
'SI' ('NEG' (bb) (aa) (1.0D-12 * aa)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
* Case 2: internal pressure larger than critical pressure
*
P1 = 1.0D7 ;
*
R1 = (P1 '/' P0) '**' (1. '/' gamf) ;
R1 = R1 '*' RHO0 ;
U1 = (2 '*' gamf '/' (gamf '-' 1.0D0)) '*' ((P0 '/' RHO0) '-'
(P1 '/' R1)) ;
U1 = U1 '**' 0.5 ;
*
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' P1 ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'RESE' ;
aa = 'MAXIMUM' ('EXCO' RCHLIM 'RN') ;
'SI' ('NEG' (r1) (aa) (1.0D-12 * aa)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
bb = 'MAXIMUM' ('EXCO' RCHLIM 'PN') ;
'SI' ('NEG' (p1) (bb) (1.0D-12 * bb)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
cc = 'MAXIMUM' ('EXCO' RCHLIM 'UX') ;
dd = 'MAXIMUM' ('EXCO' RCHLIM 'UY') ;
ee = (cc * cc) '+' (dd * dd) ;
ee = ee '**' 0.5 ;
'SI' ('NEG' (ee) (u1) (1.0D-12 * ee)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
aa = 'MAXIMUM' ('EXCO' RCHRES 'VAR1') ;
bb = SURSVOL '*' r1 * u1 ;
'SI' ('NEG' (bb) (aa) (1.0D-12 * aa)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
aa = htot * r1 * u1 * sursvol ;
bb = 'MAXIMUM' ('EXCO' RCHRES 'VAR4') ;
'SI' ('NEG' (bb) (aa) (1.0D-12 * aa)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
* Case 3: internal pressure larger than reservoir pressure
*
P1 = 210.0D5 ;
*
*
PN = 'MANUEL' 'CHPO' ($DOMINT . 'CENTRE') 1 'SCAL' P1 ;
RCHLIM RCHRES = 'KONV' 'VF' 'PERFMULT' 'CLIM' 'RESI' MDOMINT MP2P3
PGAS LISTC LISTP
RN VN PN YN CHPLI 'RESE' ;
aa = 'MAXIMUM' ('EXCO' RCHLIM 'RN') ;
'SI' ('NEG' (rho0) (aa) (1.0D-12 * aa)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
bb = 'MAXIMUM' ('EXCO' RCHLIM 'PN') ;
'SI' ('NEG' (p0) (bb) (1.0D-12 * bb)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
cc = 'MAXIMUM' ('EXCO' RCHLIM 'UX') ;
dd = 'MAXIMUM' ('EXCO' RCHLIM 'UY') ;
ee = (cc * cc) '+' (dd * dd) ;
ee = ee '**' 0.5 ;
'SI' ('NEG' (ee) (0.0) (1.0D-12 * ucr)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
aa = 'MAXIMUM' ('EXCO' RCHRES 'VAR1') ;
bb = 0.0 ;
'SI' ('NEG' (bb) (aa) (1.0D-12 * rcr * ucr)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
aa = 0.0 ;
bb = 'MAXIMUM' ('EXCO' RCHRES 'VAR4') ;
'SI' ('NEG' (bb) (aa) (1.0D-12 * rcr * (ucr ** 3))) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
*
cc = 'MAXIMUM' ('EXCO' RCHRES 'VAR2') ;
dd = 'MAXIMUM' ('EXCO' RCHRES 'VAR3') ;
ee = (cc * cc) '+' (dd * dd) ;
ee = ee '**' 0.5 ;
ff = SURSVOL * p1 ;
'SI' ('NEG' (ee) (ff) (1.0D-12 * ee)) ;
'MESSAGE' 'RESE BC' ;
'ERREUR' 5 ;
'FINSI' ;
'FIN' ;
© Cast3M 2003 - All rights reserved.
Disclaimer