Sample sludge code input file

*-------------------------------------------------------------
CONTROL
*-------------------------------------------------------------
  TITLE

CHXMIX: Illustrate heat transfer in composite heat sinks when
        the gas thermal conductivity in the gap between the
        heat sinks is derived from the gas mixture in a region.

  END TITLE
*
  TIMING
    TSTART      0.     !  START TIME, >0 FOR RESTART RUN
    TLAST   36000.     !  END TIME (Seconds) - 10 HOURS
    DTMIN       0.001  !  MIN TIMESTEP (Seconds)
    DTMAX              !  MAX TIMESTEP (Seconds)
        0.      3.0
    36000.      3.0
        1.0E6   3.0
    DTPRIN   3600.     !  PRINT INTERVAL (Seconds)
    PLTMIN     10.     !  MIN PLOT INTERVAL (Seconds)
    PLTMAX    100.     !  MAX INTERVAL WITHOUT PLOT (Seconds)
    DTRST    3600.     !  RESTART INTERVAL (Seconds)
    FTPCH       0.005  !  FRACTIONAL CHANGE IN T AND P
    FAECH       0.005  !  FRACTIONAL CHANGE IN AEROSOL MASS
    FPPLCH      0.03   !  FRACTIONAL CHANGE FOR PLOTTING
  END TIMING
*
  COMPOUNDS 2
    ARGON HELIUM
  END COMPOUNDS
*
  PRINT-T

    1,1  1,8  1,15       2,1  2,8  2,15

    3,1  3,8  3,15       4,1  4,8  4,15

    5,1  5,8  5,15       6,1  6,8  6,15

  END PRINT-T
*
  PLOT 1
    PRESSURE     4    1  2  3 10    ! Pressure
    GAS-T        4    1  2  3 10    ! Gas temperature
    HS-TI        6    1  2  3  4  5  6 ! Inner heat sink temperature
    HS-TA        6    1  2  3  4  5  6 ! Average heat sink temperature
    HS-TO        6    1  2  3  4  5  6 ! Outer heat sink temperature
    QSW          6    1  2  3  4  5  6 ! Heat input to sandwich heat sink
    QGAS-HSI     6    1  2  3  4  5  6 ! Heat input to heat sink from gas
    QGAS-HSO     6    1  2  3  4  5  6 ! Heat input to heat sink from gas
    QRAD-HSI     6    1  2  3  4  5  6 ! Rad heat input to heat sink
    QRAD-HSO     6    1  2  3  4  5  6 ! RAd heat input to heat sink
    GAS-X    HELIUM 4 1  2  3 10       ! Gas Mole Fraction
    GAS-X    ARGON  4 1  2  3 10       ! Gas Mole Fraction
    GAS-MASS HELIUM 4 1  2  3 10       ! Gas Mass (Species), kg
    GAS-MASS ARGON  4 1  2  3 10       ! Gas Mass (Species), kg
  END PLOT
*
  PLOT 2
    HS-T  1      5    1  2  3  4  5 ! Heat sink 1 detail
    HS-T  1      5    6  7  8  9 10 ! Heat sink 1 detail
    HS-T  1      5   11 12 13 14 15 ! Heat sink 1 detail
    HS-T  2      5    1  2  3  4  5 ! Heat sink 2 detail
    HS-T  2      5    6  7  8  9 10 ! Heat sink 2 detail
    HS-T  2      5   11 12 13 14 15 ! Heat sink 2 detail
    HS-T  3      5    1  2  3  4  5 ! Heat sink 3 detail
    HS-T  3      5    6  7  8  9 10 ! Heat sink 3 detail
    HS-T  3      5   11 12 13 14 15 ! Heat sink 3 detail
    HS-T  4      5    1  2  3  4  5 ! Heat sink 4 detail
    HS-T  4      5    6  7  8  9 10 ! Heat sink 4 detail
    HS-T  4      5   11 12 13 14 15 ! Heat sink 4 detail
    HS-T  5      5    1  2  3  4  5 ! Heat sink 5 detail
    HS-T  5      5    6  7  8  9 10 ! Heat sink 5 detail
    HS-T  5      5   11 12 13 14 15 ! Heat sink 5 detail
    HS-T  6      5    1  2  3  4  5 ! Heat sink 6 detail
    HS-T  6      5    6  7  8  9 10 ! Heat sink 6 detail
    HS-T  6      5   11 12 13 14 15 ! Heat sink 6 detail
  END PLOT
*
* Run diagnostics
*
  PLOT 3
    DELT
    ITDX
    ITERJN
    CPU
  END PLOT
*
  REPORT 1
    QSW      INDEX=1   TIME=1800
    QSW      INDEX=3   TIME=1800
    QSW      INDEX=5   TIME=1800

    QSW      INDEX=1   TIME=7200
    QSW      INDEX=3   TIME=7200
    QSW      INDEX=5   TIME=7200

    QSW      INDEX=1   END
    QSW      INDEX=3   END
    QSW      INDEX=5   END

    GAS-X    HELIUM INDEX=3 TIME=1800
    GAS-X    ARGON  INDEX=3 TIME=1800

    GAS-X    HELIUM INDEX=3 TIME=7200
    GAS-X    ARGON  INDEX=3 TIME=7200

    GAS-X    HELIUM INDEX=3 END
    GAS-X    ARGON  INDEX=3 END
  END REPORT
*
  ACTIVE MODELS
    IJUNC   1       ! Junction flow model
    ICCFLW  0       ! Counter-current flow model
    IHSINK  1       ! Heat sinks
    ICNDS   0       ! Condensation
    IHXPOL  0       ! Heat and Mass transfer on liquid pool
    IASED   0       ! Aerosol Sedimentation
    IALEAK  0       ! Aerosol Leakage
    IFOG    0       ! Fog formation
    ISRC    1       ! User-defined sources
    IMCO    0       ! MCO models
    ISENS   0       ! Sensitivity runs
  END ACTIVE MODELS
*
  OPTIONS
    CELSIUS
    IPLTYP  2       ! 1=wrap around, 2=no wrap (spreadsheet format)
  END OPTIONS
*
*-----------------------------------------------------------------------
* SOURCE GROUP:
* ENTER: TIME, TEMP, FLOWRATES, POWER
*
  SOURCES
    SOURCE 1
      REGION  3  GASES 1
        PLUME  RADIUS 1.0 ELEVATION 10.0 DOWN
        ARGON
*     Time      Temp    kg/s    Power
          0.0   50.0    0.0       0.E0
       3600.0   50.0    0.08      0.E0
      32400.0   50.0    0.30      0.E0
      END REGION
  END SOURCES
*
END CONTROL
*
*-------------------------------------------------------------------
VOLUMES
*-------------------------------------------------------------------
*
  REGIONS         1           10
    LABEL     HE_TANK        AMBIENT
    VOLUME     1000.0          1.0E6
    ZTOP         10.0        100.0
    ELEVATION     0.0         -1.0
    SED_AREA     50.0       1000.0
    TEMP_GAS     50.0         20.0
    PRESSURE      1.01325E5    1.01325E5
  END REGIONS
*
  REGIONS           2            3
    LABEL     AR_TANK   PURGE_TANK
    SAME_AS         1            1
  END REGIONS
*
* Gas composition of each region; specify mole fraction of each
* gas. No more than five columns at a time.
*
  GASES          1       2       3      10
    HELIUM       1.0     0.0     1.0     1.0
    ARGON        0.0     1.0     0.0     0.0
  END GASES
*
END VOLUMES
*
*-------------------------------------------------------------------
HEAT_SINKS
*-------------------------------------------------------------------
* Syntax:
* IGEOM   1 for plane, 0 or 2 for cylinder
* IMATHS  material type
* XRI     Inner Radius (m)
* XRO     Outer Radius (m) for cylindrical,
*         thickness (m) for planar
* AHS     One-sided average heat sink area (m^2)
* XZHS    Axial length for conduction ( m)
* TIINIT  Initial inside surface temperature (C)
* TOINIT  Initial outside surface temperature (C)
* IMSLAB  Number of slabs, 3 is minimum. Use 1 for lumped heat
*         sink
* IREGI   Region index for inner surface or 0 (insulated)
*         or -1 for constant temperature
* TIHS    Region surface temperature when IREGI = -1  (C)
* IREGO   Region index for outer surface or 0 (insulated)
*         or -1 for constant temperature
* TOHS    Region surface temperature when IREGO = -1  (C)
* XLHS    Characteristic length for natural convection (m)
* If different from material type:
* RHO     Density (kg/m^3),
* KHS     Thermal Conductivity (W/m/K)
* CPHS    Specific Heat (J/kg/K)
* QV      Volumetric Heat Generation (W/m^3)
* EHSI    Emissivity of inner surface
*         User has an option to input a temperature dependent
*         emissivity by inputting a negative integer for the
*         emissivity and supplying the corresponding temperature
*         versus emissivity look-up table. See TTABLE & ETABLE
*         keywords.
* EHSO    Emissivity of outer surface
******************************************************************
*
  SINKS        1        2
    LABEL      HOT_1    COLD_1
    IGEOM      1        1      ! Planar geometry
    IMATHS     1        1      ! Copper
    IORIHS     0        0      ! Vertical orientation
    IMSLAB    15       15      ! 14 nodes plus 2 half-nodes
    IREGO     -1        0      ! Outer boundary condition: 50C vs Reg. 1
    IREGI      0       -1      ! Inner boundary condition: Reg 1. vs 20C
    TOHS      50.0      0.0    ! Fixed outer surface temperature
    TIHS       0.0     20.0    ! Fixed inner surface temperature
    XRI        0.0      0.0    ! Inner radius
    XRO        0.01     0.01   ! Thickness
    AHS        1.0      1.0    ! Mean cross-sectional area
    TOINIT    50.0     20.0    ! Outside surface initial temperature
    TIINIT    50.0     20.0    ! Inside surface initial temperature
    XZHS       1.0      1.0    ! Axial conduction length
    ZTHS       2.0      2.0    ! Upper coordinate
    ZBHS       1.0      1.0    ! Lower coordinate
    VHS        0.01     0.01   ! Volume
  END
*
  SINKS        3        4
    LABEL      HOT_2    COLD_2
    SAME_AS    1        2
    IREGO     -1        0      ! Outer boundary: 50C vs Reg. 2
    IREGI      0       -1      ! Inner boundary: Reg 2. vs 20C
  END
*
  SINKS        5        6
    LABEL      HOT_3    COLD_3
    SAME_AS    1        2
    IREGO     -1        0      ! Outer boundary: 50C vs Reg. 3
    IREGI      0       -1      ! Inner boundary: Reg 3. vs 20C
  END
*
  MATERIAL LIBRARY
*   Name            IMATHS   RHO   KHS      CP    QV  EHSI    EHSO
    COPPER          1      8954.  398.0   384.    0.   0.7     0.7
    STAINLESS-STEEL 2      8000.   16.0   500.    0.   0.3     0.3
  END MATERIAL LIBRARY
*
* Composite heat sinks:
*
  COMPOSITE NET=1  1             2
    IFLIP          0             0
    FMULT          1.0           1.0
    XGAP                0.01
    KGAP                HELIUM
  END COMPOSITE
*
  COMPOSITE NET=2  3             4
    IFLIP          0             0
    FMULT          1.0           1.0
    XGAP                0.01
    KGAP                ARGON
  END COMPOSITE
*
  COMPOSITE NET=3  5             6
    IFLIP          0             0
    FMULT          1.0           1.0
    XGAP                0.01
    KGAP                REGION-3
  END COMPOSITE
*
END HEAT_SINKS
*
*-------------------------------------------------------------------
JUNCTIONS
*-------------------------------------------------------------------
  PATHS      1
    LABEL    VENT
    IJTYP    1
    ICCJN    1
    IHORIZ   0
    AJN      0.01
    CJN      1.0
    IR1      3
    IR2     10
    Z1JN    10.0
    Z2JN    11.0
  END PATHS
END JUNCTIONS
*-------------------------------------------------------------------
* __END__