TALK=F;RUN( 1, 1)

TEXT( Library case Y607: Combustion chamber with varying absorption

  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    PLANT information :
     * Data input groups used:  9
     * Ground groups planted :  1, 9-7, 19-6
     * Headings used  : PRPT??, SC06??
     * Functions used : None
     * Commands used  : IF.
  <<<<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

 LOAD(492)

NAMSAT=MOSG
STORE(ABSR)
    >>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     The above command provides the 3D-storage for aborption
     coefficient.
    <<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
RG(1)=0.2;RG(2)=0.3
    >>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     The use of RG's, as above, in PLANTed expression is especially
     convenient for parametric studies.
    <<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
   ABSR=RG(1)*FUEL+RG(2)*PROD
    >>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     The above statement makes the absorbtion coefficient to be the
     function of the fuel and combustion products mass fractions. It is
     done for the whole domain at the end of each iz-slab.
    <<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
PRNDTL(RADX)=-GRND
   LAMPR(RADX)=1/(:SCAT:+ABSR)
PRNDTL(RADZ)=-GRND
   LAMPR(RADZ)=1/(:SCAT:+ABSR)
    >>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     The above two statements followed by their pointers instruct PLANT
     to make the diffusivities of radiation fluxes in X- and
     Z-directions reciprocal to the sum of scattering and absorption
     coefficients.
    <<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
PRNDTL(RADY)=-GRND
   LAMPR(RADY)=1/(:SCAT:+ABSR)
  IF(CARTES)
   LAMPR(RADY)=RV2D**2/(1.+RV2D*(:SCAT:+ABSR))
  IF(.NOT.CARTES)
    >>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     PLANT makes the provision for the calculations of diffusivity of
     lateral radiation flux. The first statement of the above two is
     applied for cartesian coordinates, while the second will be active
     for cylindrical polar ones.
    <<<<<<<<<<<<<<<<<<<<  Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
STOP