c@a
c@versb
C-----------------------------------------------------------------------
C
CVERS Code_Saturne version 1.3
C ------------------------
C
C This file is part of the Code_Saturne Kernel, element of the
C Code_Saturne CFD tool.
C
C Copyright (C) 1998-2007 EDF S.A., France
C
C contact: saturne-support@edf.fr
C
C The Code_Saturne Kernel is free software; you can redistribute it
C and/or modify it under the terms of the GNU General Public License
C as published by the Free Software Foundation; either version 2 of
C the License, or (at your option) any later version.
C
C The Code_Saturne Kernel is distributed in the hope that it will be
C useful, but WITHOUT ANY WARRANTY; without even the implied warranty
C of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
C GNU General Public License for more details.
C
C You should have received a copy of the GNU General Public License
C along with the Code_Saturne Kernel; if not, write to the
C Free Software Foundation, Inc.,
C 51 Franklin St, Fifth Floor,
C Boston, MA 02110-1301 USA
C
C-----------------------------------------------------------------------
c@verse
C cpincl.h
C
C***********************************************************************
C
C INCLUDE POUR LA PHYSIQUE PARTICULIERE RELATIF
C A LA COMBUSTION DU CHARBON PULVERISE
C
C Necessite ppppar.h et ppthch
C-----------------------------------------------------------------------
C
C--> POINTEURS VARIABLES COMBUSTION CHARBON PULVERISE cpincl, ppincl
C
C -------voir ppppar
CC NCHARM --> Nombre maximal de charbons
CC NCPCMX --> Nombre maximal de classes par charbon
CC NCLCPM --> Nombre total de classes
CC
C INTEGER NCHARM , NCPCMX , NCLCPM
C PARAMETER (NCHARM=3, NCPCMX=10, NCLCPM=NCHARM*NCPCMX)
C
CC ------
C EPSICP : Precision pour les tests
C
DOUBLE PRECISION EPSICP
PARAMETER ( EPSICP = 1.D-8 )
C
C
C--> DONNEES RELATIVES AU CHARBON
C
C NCHARB --> Nombre de charbons
INTEGER NCHARB
C
C ---- PAR CHARBON (grandeurs fournies)
C
C - Distribution granulometrique
C NCLPCH(CH) --> Nombre de classes par charbon
C
INTEGER NCLPCH(NCHARM)
C
C - Proprietes sur charbon sec
C CCH(CH) --> Composition elementaire en C, H, O sur sec (%)
C HCH(CH) du charbon
C OCH(CH)
C ALPHA(CH) --> Composition du charbon reactif
C BETA(CH) sous la forme CH(ALPHA)O(BETA)
C ALPHA(CH) = HCH(CH)/CCH(CH)
C BETA(CH) = OCH(CH)/CCH(CH)
C ALPHAM --> ALPHA moyen
C BETAM --> BETA moyen
C PCICH(CH) --> PCI (J/kg) charbon
C RHO0CH(CH) --> Masse volumique initiale (kg/m3)
C - Proprietes sur charbon sec du coke
C CCK(CH) --> Composition elementaire en C, H, O sur sec (%)
C HCK(CH) du coke
C OCK(CH)
C GAMMA(CH) --> Composition du coke
C DELTA(CH) sous la forme CH(GAMMA)O(DELTA)
C GAMMA(CH) = HCK(CH)/CCK(CH)
C DELTA(CH) = OCK(CH)/CCK(CH)
C PCICK(CH) --> PCI (J/kg) coke
C RHOCK(CH) --> Masse volumique coke
C - Proprietes sur charbon sec des cendres
C XASHCH(CH) --> Taux de cendre (kg/kg)
C CPASHC(CH) --> CP des cendres (J/kg/K)
C H0ASHC(CH) --> Enthalpie de formation des cendres (J/kg)
C H02CH --> H0 du Charbon
C CPCH --> CP du Charbon
C
DOUBLE PRECISION CCH (NCHARM), HCH (NCHARM), OCH (NCHARM),
& ALPHA (NCHARM), BETA (NCHARM), ALPHAM, BETAM ,
& PCICH (NCHARM), RHO0CH(NCHARM),
& CCK (NCHARM), HCK (NCHARM), OCK (NCHARM),
& GAMMA (NCHARM), DELTA (NCHARM),
& RHOCK (NCHARM), PCICK (NCHARM),
& XASHCH(NCHARM), CPASHC(NCHARM),
& H0ASHC(NCHARM),
& H02CH (NCHARM), CP2CH (NCHARM)
C
C - Parametres cinetiques pour la devolatilisation
C (Modele de Kobayashi)
C IY1CH(CH) --> Indicateur : 0 si MVl = {CH4;CO}
C 1 si MVl = {CHz;CO}
C Y1CH(CH) --> Coefficient stoechiometrique (adim)
C calcule si IY1CH = 0 ; donne si IY1CH = 1
C A1CH(CH) --> Facteur pre-exponetielle (1/s)
C E1CH(CH) --> Energie d'activation (J/mol)
C IY2CH(CH) --> Indicateur : 0 si MVL = {C2H4;CO}
C 1 si MVL = {CxHy;CO}
C Y2CH(CH) --> Coefficient stoechiometrique (adim)
C calcule si IY2CH = 0 ; donne si IY2CH = 1
C A2CH(CH) --> Constante preexponetielle (1/s)
C E2CH(CH) --> Energie d'activation (J/mol)
C - Parametres cinetiques pour la combustion heterogene du coke
C (Modele a sphere retrecissante)
C AHETCH(CH) --> Constante pre-exponentielle (kg/m2/s/atm)
C EHETCH(CH) --> Energie d'activation (kcal/mol)
C IOCHET(CH) --> Ordre de la reaction 0.5 si = 0 1 si = 1
C
INTEGER IY1CH (NCHARM), IY2CH (NCHARM)
INTEGER IOCHET (NCHARM)
DOUBLE PRECISION Y1CH (NCHARM), A1CH (NCHARM), E1CH (NCHARM),
& Y2CH (NCHARM), A2CH (NCHARM), E2CH (NCHARM),
& AHETCH(NCHARM), EHETCH(NCHARM)
C
C - Enthalpie du charbon reactif, coke et cendres
C ICH(CH) --> Pointeur dans le tableau EHSOLI pour
C le Charbon Reactif
C ICK(CH) --> Pointeur dans le tableau EHSOLI pour
C le Coke
C IASH(CH) --> Pointeur dans le tableau EHSOLI pour
C les cendres
C NSOLID --> Nb constituants solides (Ch.Reactif, Coke, Ash)
C NSOLIM --> Nb maximal de constituants solides
C EHSOLI(S,IT) --> Enthalpie massique (J/kg) du constituant solide
C no S a la temperature T(IT)
C WMOLS(S) --> Masse molaire du constituant solide
C
INTEGER NSOLIM
PARAMETER( NSOLIM = 3*NCHARM )
C
INTEGER NSOLID, ICH(NCHARM), ICK(NCHARM), IASH(NCHARM)
DOUBLE PRECISION EHSOLI(NSOLIM,NPOT), WMOLS(NSOLIM)
C
C ---- PAR CLASSES (grandeurs deduites)
C
C NCLACP --> Nb de classes
C
INTEGER NCLACP
C
C - Proprietes
C ICHCOR(CL) --> = ICH si la classe consideree appartient
C au charbon ICH (1, 2, ...)
C DIAM20(CL) --> Diametre initial (m)
C DIA2MN(CL) --> Diametre minimum (m)
C RHO20(CL) --> Masse volumique initiale (kg/m3)
C RHO2MN(CL) --> Masse volumique minimale (kg/m3)
C XMP0(CL) --> Masse initiale de la particule (m)
C XMASH(CL) --> Masse de cendres de la particule (m)
C
INTEGER ICHCOR(NCLCPM)
DOUBLE PRECISION DIAM20(NCLCPM), DIA2MN(NCLCPM),
& RHO20 (NCLCPM), RHO2MN(NCLCPM),
& XMP0 (NCLCPM), XMASH (NCLCPM)
C
C
C--> DONNEES RELATIVES A LA COMBUSTION DES ESPECES GAZEUSES
C
C ICHX1C(CH) --> Pointeur CHx1 pour EHGAZE et WMOLE
C ICHX2C(CH) --> Pointeur CHx2 pour EHGAZE et WMOLE
C ICHX1 --> Pointeur CHx1m pour EHGAZE et WMOLE
C ICHX2 --> Pointeur CHx2m pour EHGAZE et WMOLE
C ICO --> Pointeur CO pour EHGAZE et WMOLE
C IO2 --> Pointeur O2 pour EHGAZE et WMOLE
C ICO2 --> Pointeur CO2 pour EHGAZE et WMOLE
C IH2O --> Pointeur H2O pour EHGAZE et WMOLE
C IN2 --> Pointeur N2 pour EHGAZE et WMOLE
C XSI --> XSI = 3,76 pour de l'air
C F3MAX --> Maximum pour le traceur F3
C CHX1(CH) --> Composition de l'hydrocarbure relatif
C au MVl : CH(X1)
C CHX2(CH) --> Composition de l'hydrocarbure relatif
C au MVL : CH(X2)
C A1(CH), --> Coefficients stoechiometriques molaires pour
C B1(CH) la reaction de devolatilisation a basses T
C A2(CH), --> Coefficients stoechiometriques molaires pour
C B2(CH) la reaction de devolatilisation a basses T
C
INTEGER ICHX1C(NCHARM), ICHX2C(NCHARM),
& ICHX1, ICHX2, ICO, IO2, ICO2, IH2O, IN2
DOUBLE PRECISION XSI, F3MAX,
& CHX1(NCHARM), CHX2(NCHARM),
& A1(NCHARM), B1(NCHARM),
& A2(NCHARM), B2(NCHARM)
C
C--> DONNEES COMPLEMENTAIRES RELATIVES AU CALCUL DE RHO
C SUR LES FACETTES DE BORD
C
C IENTAT(IENT) --> Indicateur air par type de facette d'entree
C IENTCP(IENT) --> Indicateur CP par type de facette d'entree
C TIMPAT(IENT) --> Temperature en K pour l'air relative
C a l'entree IENT
C X20(IENT, --> Fraction massique dans le melange de charbon
C ICLA ) de la classe ICLA relative a l'entree IENT
C
INTEGER IENTAT(NOZPPM), IENTCP(NOZPPM)
DOUBLE PRECISION TIMPAT(NOZPPM), X20(NOZPPM,NCLCPM)
C
C--> POINTEURS DANS LE TABLEAU TBMCR
C
INTEGER IF1MC(NCHARM) , IF2MC(NCHARM)
INTEGER IX1MC ,IX2MC, ICHX1F1, ICHX2F2
INTEGER ICOF1, ICOF2
C
C--> DEFINITION DES COMMONS
C
COMMON / ICPCOM / NCHARB, NCLPCH, NCLACP, IY1CH , IY2CH ,
& ICHCOR,
& NSOLID, ICH , ICK , IASH ,
& ICHX1C, ICHX2C,
& ICHX1 , ICHX2 ,
& ICO , IO2 , ICO2 , IH2O , IN2 ,
& IENTAT, IENTCP,
& IF1MC , IF2MC ,
& IX1MC , IX2MC , ICHX1F1, ICHX2F2,
& ICOF1 , ICOF2 , IOCHET
C
COMMON / RCPCOM / CCH , HCH , OCH ,
& ALPHA , BETA , ALPHAM, BETAM ,
& PCICH , RHO0CH,
& CCK , HCK , OCK ,
& GAMMA , DELTA ,
& RHOCK , PCICK ,
& XASHCH, CPASHC, H0ASHC,
& H02CH , CP2CH ,
& Y1CH , A1CH , E1CH , Y2CH , A2CH , E2CH ,
& AHETCH, EHETCH,
& EHSOLI, WMOLS ,
& DIAM20, DIA2MN, RHO20 , RHO2MN,
& XMP0 , XMASH ,
& XSI , F3MAX ,
& CHX1 , CHX2 ,
& A1 , B1 , A2 , B2 ,
& TIMPAT, X20
C
C--> GRANDEURS FOURNIES PAR L'UTILISATEUR EN CONDITIONS AUX LIMITES
C PERMETTANT DE CALCULER AUTOMATIQUEMENT LA VITESSE, LA TURBULENCE,
C L'ENTHALPIE D'ENTREE.
C
C POUR LES ENTREES UNIQUEMENT , IENT ETANT LE NUMERO DE ZONE FRONTIERE
C
C QIMPAT(IENT) --> Debit Air en kg/s
C TIMPAT(IENT) --> Temperature Air en K
C QIMPCP(IENT,ICHA) --> Debit Charbon ICHA en kg/s
C TIMPCP(IENT,ICHA) --> Temperature Charbon ICHA en K
C DISTCH(IENT,ICHA,ICLA) --> Distribution en %masse de la classe ICLA
C pour le charbon ICHA
C
DOUBLE PRECISION QIMPAT(NOZPPM)
DOUBLE PRECISION QIMPCP(NOZPPM,NCHARM), TIMPCP(NOZPPM,NCHARM)
DOUBLE PRECISION DISTCH(NOZPPM,NCHARM,NCPCMX)
c DOUBLE PRECISION COEFE(NGAZEM), XSOLID(NSOLIM)
c DOUBLE PRECISION F1MC(NCHARM) , F2MC(NCHARM)
C
COMMON / RCPCLI / QIMPAT ,
& QIMPCP , TIMPCP ,
& DISTCH
c & COEFE , XSOLID ,
c & F1MC , F2MC
C
C FIN
c@z
syntax highlighted by Code2HTML, v. 0.9.1