# This file was generated by Barracuda Virtual Reactor. It will be re-written 
# the next time Barracuda Virtual Reactor is run. If you modify this file,    
# expect your changes to be lost.                                   
#
version 2300
# This file was generated by GUI version 25.0.0.i057.
projectUniqueIdStart

projectUniqueIdEnd
storeProjectFileAsUniqueId      1
projectCommentStart


projectCommentEnd
#
# =========== grid
grid
admesh 0
gridLastUsedTotalNumberCells 20000
rm_small_cells 3
rm_small_vol_I 0.040000
rm_small_vol_E 0.040000
max_aspect_ratio_I 30.000000
max_aspect_ratio_E 15.000000
open_closed_faces 1
rm_small_max_pass_combine 150
poly_col 0
draw_poly_faces 0
eps_scale -1.000000
stl 1
# x-Min and x-Max all stl files
-5.000000e+01  5.000000e+01
# y-Min and y-Max all stl files
-4.996954e+01  4.996954e+01
# z-Min and z-Max all stl files
 0.000000e+00  1.000000e+03
#
# =========== start stl file properties
BFB_reactor.stl
# header 
header ""Geom<stl unit=MM>""
# number facets 
180
# x-Min and x-Max 
-5.000000e+01  5.000000e+01
# y-Min and y-Max 
-4.996954e+01  4.996954e+01
# z-Min and z-Max 
 0.000000e+00  1.000000e+03
# =========== end stl file properties
x0 -5.000000e+01
y0 -4.996950e+01
z0  0.000000e+00
xdist 1
12  5.000000e+01   0.0000
ydist 1
12  4.996950e+01   0.0000
zdist 1
126  1.000000e+03   0.0000
units mm
endgrid
# =========== endgrid

# Variable units for solver output
InputOutputUnits
  time                    "s"               
  length                  "m"               
  mass                    "kg"              
  temperature             "K"               
  mole                    "mol"             
  pressure                "Pa"              
  pressureGradient        "Pa/m"            
  velocity                "m/s"             
  massFlux                "kg/(s*m^2)"      
  density                 "kg/m^3"          
  heatTransferCoef        "J/(s*m^2*K)"     
  heatTransferFlux        "J/(s*m^2)"       
  heatTransferRate        "J/s"             
  volAveReactionRate      "kmol/(s*m^3)"    
  discreteReactionRate    "kg/s"            
  particleSize            "micron-diameter" 
  massFlowRate            "kg/s"            
  volumeFlowRate          "m^3/s"           
  massConcentration       "kg/m^3"          
  moleConcentration       "mol/m^3"         
  viscosity               "Pa*s"            
  kinematicViscosity      "m^2/s"           
  thermalConductivity     "J/(s*m*K)"       
  specificHeat            "J/(kg*K)"        
  specificEnthalpy        "J/kg"            
  molecularWeight         "g/mol"           
  area                    "m^2"             
  volume                  "m^3"             
  energy                  "J"               
  wear                    "kg/(m^2*year)"   
  attrition               "kg^a*m^b/s^b"    
  cpuTime                 "hr"              
  simSpeed                "s/day"           
endInputOutputUnits

# Post-processor graphics mesh viewer to use for project
outputGmv 0
outputTecplot 1

# Thermal energy
isothermal      1
toggleThermal   on
radiationOn     0
isothermalT     8.731500e+02
presTermThermal 0
enableFluidTempSolver 0
outputMinMaxTemp 0

pHeatConduction
  enabled   Off
  iterMin   1
  iterMax   100
  errMax    1.000000e-12
end_pHeatConduction

toggleChem  on   0.000000e+00   0.000000e+00

# Put the number of fluids near start of input
numFluids  2

implicitParticleBuoyCalc on

# Put the use particles flag near start of input
useParticles 1

# Put the noParticleMomentum flag after useParticles
noParticleMomentum 0
allowPartialClouds on
virtualMassForce off
liftForce off
subgridWallModel off
wallFriction 0.000000e+00
coalescenceBreakupModel
  enabled off
  coalescenceFactor 1.000000e+00
  maximumBubbleSize 1.000000e-01
end_coalescenceBreakupModel
compressibleBubbles on

# Min and max temperature limit
minTempLimit  1.000000e+02
maxTempLimit  6.000000e+03

#Advanced GPU Options
decompSliceNormal automatic
tagParticlesByDecomposition OFF

# Put the name of grid file near start of input but after numFluids
gridFileName 00grid.grd

# minAreaBc used for minimum area percentage on flow/pres BC
minAreaBc    0.040000

# Time step control
time 1
#         dt       end time
#  ------------- -------------
    5.000000e-04  1.800000e+03
end_timeControl

# write the ray trace flag
rayTrace 1

# Iterations and convergance
maxErrPres     1.000000e-06
maxErrVel      1.000000e-07
maxErrVol      1.000000e-07
maxErrEnergy   1.000000e-06
maxErrFluidT   1.000000e-01
maxErrRadiate  1.000000e-09
maxIterPres    2000
maxIterPres0   2000
maxIterVel     50
maxIterVol     10
maxIterEnergy  100
maxIterFluidT  20
maxIterRadiate 1500
maxIterLoop    1

# Print, plot and restart frequency
printFreq         1.000000e-07
restartFreq       1.000000e+02
restartClockFreq  1.200000e+02

vizOutputInterval
  fileFlag     Off
  interval     6.000000e+01
  file         ""
end_vizOutputInterval

systemBalanceLogs
  logMass         Off
  massInterval    1.000000e-07
  logEnergy       Off
  energyInterval  1.000000e-07
end_systemBalanceLogs

# Time interval adjustment flags
adjDtHighVolFrac  1
adjDtHighCFL      1
adjDtHighVolErr   0
adjDtHighT        1
adjDtHighChemProd 1

# Min and max CFL
minCfl  8.000000e-01
maxCfl  1.500000e+00

# maxDifTemp - maximum cell difference in temperature
maxDifTemp       3.000000e+02
maxChangeChem 1.000000e-01

# Turn on/off pressure work & drag work in energy equation - 1=on, 0=off
presWork 0

# Material phases (and properties), must be before init. gasMixturePropAve: 0=mass-based, 1=mole-based
MaterialPhases 9
  nMaterials        9
  gasMixturePropAve 1
  equState          1
  compressible      1
  plotDeltaT        1.000000e+00
  plotTMin          0.000000e+00
  plotTMax          2.000000e+03
  evalT             3.000000e+02
  material
    chem_formula           "CH4_2"
    formula_description    "CH4 METHANE. ANHARMONIC"
    state                  G
    molecular_weight       1.604230e+01
    absorption_model OFF
    expression density
      polynomial  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression viscosity
      polynomial  1.596000e-06  3.439000e-08 -8.140000e-12  0.000000e+00  0.000000e+00 
      value_units "Pa*s"
      temperature_units "K"
      temp_limits 1 2.000000e+02 2.112000e+03
      value_limits 0 1.000000e-07 1.000000e-03
    end_expression
    delta_hf(298)          -4.650434e+06
    expression heat_capacity
      double_polynomial
        lowerCoefs  2.668848e+03 -7.081283e+00  2.547265e-02 -2.509909e-05  8.635257e-09 
        upperCoefs  8.569057e+02  5.196757e+00 -1.719040e-03  2.780657e-07 -1.631111e-11 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial -9.350000e-03  1.402800e-04  3.318000e-08  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 9.700000e+01 5.000000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
    absorption_coef         1.000000e-05
    scattering_coef         0.000000e+00
    refractive_index        1.000000e+00
    expression mass_diffusivity
      polynomial  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "m^2/s"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 1.000000e-06 1.000000e-04
    end_expression
    molarVolume 3.783200e-02 # [m3/kmol]
    henrysLaw 0.000000e+00
  end_material
  material
    chem_formula           "CO"
    formula_description    "CO CARBON MONOXIDE"
    state                  G
    molecular_weight       2.801010e+01
    absorption_model OFF
    expression density
      polynomial  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression viscosity
      polynomial  3.228000e-06  4.747000e-08 -9.648000e-12  0.000000e+00  0.000000e+00 
      value_units "Pa*s"
      temperature_units "K"
      temp_limits 1 6.800000e+01 2.460000e+03
      value_limits 0 1.000000e-07 1.000000e-03
    end_expression
    delta_hf(298)          -3.946267e+06
    expression heat_capacity
      double_polynomial
        lowerCoefs  1.062599e+03 -1.811859e-01  3.018454e-04  2.692483e-07 -2.684821e-10 
        upperCoefs  9.049553e+02  4.012659e-01 -1.442099e-04  2.340802e-08 -1.394642e-12 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial  5.066300e-04  9.123500e-05 -3.523800e-08  8.198200e-12  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
    absorption_coef         1.000000e-05
    scattering_coef         0.000000e+00
    refractive_index        1.000000e+00
    expression mass_diffusivity
      polynomial  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "m^2/s"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 1.000000e-06 1.000000e-04
    end_expression
    molarVolume 3.545600e-02 # [m3/kmol]
    henrysLaw 0.000000e+00
  end_material
  material
    chem_formula           "CO2"
    formula_description    "CO2 CARBON DIOXIDE"
    state                  G
    molecular_weight       4.400950e+01
    absorption_model OFF
    expression density
      polynomial  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression viscosity
      polynomial  2.545000e-06  4.549000e-08 -8.649000e-12  0.000000e+00  0.000000e+00 
      value_units "Pa*s"
      temperature_units "K"
      temp_limits 1 1.950000e+02 2.629784e+03
      value_limits 0 1.000000e-07 1.000000e-03
    end_expression
    delta_hf(298)          -8.941945e+06
    expression heat_capacity
      double_polynomial
        lowerCoefs  4.452849e+02  1.697418e+00 -1.345608e-03  4.642706e-07 -2.699610e-11 
        upperCoefs  8.760001e+02  5.179577e-01 -1.881601e-04  3.030268e-08 -1.731022e-12 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial -1.200000e-02  1.020800e-04 -2.240300e-08  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 1.950000e+02 2.278000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
    absorption_coef         1.000000e-05
    scattering_coef         0.000000e+00
    refractive_index        1.000000e+00
    expression mass_diffusivity
      polynomial  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "m^2/s"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 1.000000e-06 1.000000e-04
    end_expression
    molarVolume 3.727800e-02 # [m3/kmol]
    henrysLaw 0.000000e+00
  end_material
  material
    chem_formula           "CaCO3"
    formula_description    "CaCO3 Calcium Carbonate (S)"
    state                  S
    molecular_weight       1.000869e+02
    expression density
      polynomial  2.230000e+03  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    delta_hf(298)          -1.206515e+07
    expression heat_capacity
      double_polynomial
        lowerCoefs -9.367063e+01  5.750395e+00 -1.012179e-02  8.350020e-06 -2.613956e-09 
        upperCoefs  1.040782e+03  3.486184e-01 -1.560031e-04  4.185288e-08 -4.248083e-12 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 3.000000e+02 1.200000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial  9.017250e+01 -1.601400e-01  7.716100e-05  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 3.000000e+02 1.200000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
  end_material
  material
    chem_formula           "CaO"
    formula_description    "CaO Calcium Oxide (S)"
    state                  S
    molecular_weight       5.607740e+01
    expression density
      polynomial  1.560000e+03  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    delta_hf(298)          -1.491872e+07
    expression heat_capacity
      double_polynomial
        lowerCoefs -9.367063e+01  5.750395e+00 -1.012179e-02  8.350020e-06 -2.613956e-09 
        upperCoefs  1.040782e+03  3.486184e-01 -1.560031e-04  4.185288e-08 -4.248083e-12 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 3.000000e+02 1.200000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial  9.017250e+01 -1.601400e-01  7.716100e-05  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 3.000000e+02 1.200000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
  end_material
  material
    chem_formula           "H2"
    formula_description    "H2 HYDROGEN. REF ELEMENT"
    state                  G
    molecular_weight       2.015800e+00
    absorption_model OFF
    expression density
      polynomial  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression viscosity
      polynomial  2.775840e-06  2.119300e-08 -3.284000e-12  0.000000e+00  0.000000e+00 
      value_units "Pa*s"
      temperature_units "K"
      temp_limits 1 1.500000e+02 3.231000e+03
      value_limits 0 1.000000e-07 1.000000e-03
    end_expression
    delta_hf(298)           0.000000e+00
    expression heat_capacity
      double_polynomial
        lowerCoefs  9.669662e+03  3.291725e+01 -8.034153e-02  8.314244e-05 -3.042427e-08 
        upperCoefs  1.209719e+04  3.409510e+00 -6.038657e-04  6.356178e-08 -2.841112e-12 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial  3.951000e-02  4.591800e-04 -6.493300e-08  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 1.500000e+02 5.000000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
    absorption_coef         1.000000e-05
    scattering_coef         0.000000e+00
    refractive_index        1.000000e+00
    expression mass_diffusivity
      polynomial  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "m^2/s"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 1.000000e-06 1.000000e-04
    end_expression
    molarVolume 2.860400e-02 # [m3/kmol]
    henrysLaw 0.000000e+00
  end_material
  material
    chem_formula           "H2O"
    formula_description    "H2O STEAM"
    state                  G
    molecular_weight       1.801520e+01
    absorption_model OFF
    expression density
      polynomial  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression viscosity
      polynomial -3.682550e-06  4.291600e-08 -1.624000e-12  0.000000e+00  0.000000e+00 
      value_units "Pa*s"
      temperature_units "K"
      temp_limits 1 2.500000e+02 3.219000e+03
      value_limits 0 1.000000e-07 1.000000e-03
    end_expression
    delta_hf(298)          -1.342417e+07
    expression heat_capacity
      double_polynomial
        lowerCoefs  1.937891e+03 -9.399066e-01  3.009481e-03 -2.532967e-06  8.178565e-10 
        upperCoefs  1.235594e+03  1.372280e+00 -3.571351e-04  4.358604e-08 -1.970368e-12 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial  5.300000e-04  4.709300e-05  4.955100e-08  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 2.000000e+02 4.000000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
    absorption_coef         1.000000e-05
    scattering_coef         0.000000e+00
    refractive_index        1.000000e+00
    expression mass_diffusivity
      polynomial  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "m^2/s"
      temperature_units "K"
      temp_limits 1 2.000000e+02 6.000000e+03
      value_limits 0 1.000000e-06 1.000000e-04
    end_expression
    molarVolume 1.806900e-02 # [m3/kmol]
    henrysLaw 0.000000e+00
  end_material
  material
    chem_formula           "MgO"
    formula_description    "MgO MAGNESIUM OXIDE (S)"
    state                  S
    molecular_weight       4.030440e+01
    expression density
      polynomial  1.560000e+03  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    delta_hf(298)          -1.491872e+07
    expression heat_capacity
      double_polynomial
        lowerCoefs -9.367063e+01  5.750395e+00 -1.012179e-02  8.350020e-06 -2.613956e-09 
        upperCoefs  1.040782e+03  3.486184e-01 -1.560031e-04  4.185288e-08 -4.248083e-12 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 3.000000e+02 1.200000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial  9.017250e+01 -1.601400e-01  7.716100e-05  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 3.000000e+02 1.200000e+03
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
  end_material
  material
    chem_formula           "Ni_1"
    formula_description    "NI REF ELEMENT (S). (I)"
    state                  S
    molecular_weight       5.869340e+01
    expression density
      polynomial  2.200000e+03  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
      value_units "kg/m^3"
      temperature_units "K"
      temp_limits 0 0.000000e+00 6.000000e+03
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    delta_hf(298)           0.000000e+00
    expression heat_capacity
      double_polynomial
        lowerCoefs  5.554784e+02 -3.317657e+00  1.901617e-02 -3.909646e-05  2.812555e-08 
        upperCoefs  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00 
        breakTemperature 1.000000e+03
      end_double_polynomial
      value_units "J/kg/K"
      temperature_units "K"
      temp_limits 1 2.500000e+02 6.310000e+02
      value_limits 0 0.000000e+00 1.000000e+00
    end_expression
    expression thermal_conductivity
      polynomial  1.444370e+02 -2.279300e-01  1.643400e-04  0.000000e+00  0.000000e+00 
      value_units "W/m/K"
      temperature_units "K"
      temp_limits 1 2.500000e+02 6.310000e+02
      value_limits 0 1.000000e-04 1.000000e+00
    end_expression
  end_material
endMaterialPhases


#liquidCompressibilty beta p_ref
liquidCompressibilty 0.000000e+00 1.013250e+05


# Particles
Particles  2
  particleDescription
    particleSpcId   1
    massPerVol      0.000000e+00
    solidMaterial
      materialName  "Ni_1"
      massFraction  1.000000e+00
      ageFactor     1.000000e+00
    endSolidMaterial
    type            value
    fileName        ""
    radius1         1.50000e+02
    radius2         2.50000e+02
    units micron-diameter
    dragModel
      dragName              "WenYu-Ergun"
      multiplierFlag        CONSTANT
      multiplier            1.00000e+00
      multiplierFilePredef  ""
      multiplierFileUser    ""
      useLegacySlipVelInMult    off
      agglomeration
        onOff       off
        type        1
        radCut       1.80000e-05
        radBlend     5.00000e-06
        radBlendVar  2.00000e-01
        radFile      ""
        units micron-diameter
      endAgglomeration
    endDragModel

    # Species cloud resolution
    cloudResolution
      type              GLOBAL
      sliderLevel       3
      npPerCell         1.250000e+02
      npPerMass         5.000000e+03
    end_cloudResolution

    useSpeciesVoidCP off
    sphericity      1.00000e+00
    emissivity      1.00000e+00
    scatterFactor   0.00000e+00
    comment         "Nickel-Catalyst"
  end
  particleDescription
    particleSpcId   2
    massPerVol      0.000000e+00
    solidMaterial
      materialName  "CaO"
      massFraction  6.190000e-01
      ageFactor     1.000000e+00
    endSolidMaterial
    solidMaterial
      materialName  "MgO"
      massFraction  3.810000e-01
      ageFactor     1.000000e+00
    endSolidMaterial
    type            value
    fileName        ""
    radius1         1.25000e+02
    radius2         3.00000e+02
    units micron-diameter
    dragModel
      dragName              "WenYu-Ergun"
      multiplierFlag        CONSTANT
      multiplier            1.00000e+00
      multiplierFilePredef  ""
      multiplierFileUser    ""
      useLegacySlipVelInMult    off
      agglomeration
        onOff       off
        type        1
        radCut       1.80000e-05
        radBlend     5.00000e-06
        radBlendVar  2.00000e-01
        radFile      ""
        units micron-diameter
      endAgglomeration
    endDragModel

    # Species cloud resolution
    cloudResolution
      type              GLOBAL
      sliderLevel       3
      npPerCell         1.250000e+02
      npPerMass         5.000000e+03
    end_cloudResolution

    useSpeciesVoidCP off
    sphericity      1.00000e+00
    emissivity      1.00000e+00
    scatterFactor   0.00000e+00
    comment         "Dolomite"
  end
endParticles

# GUI run settings (solver does not read these, only for GUI to pass command line)
runSettingsBegin
  useOpenMP               no
  openMPThreadCnt         19
  useGPU                  yes
  runSettingsGpuBegin
    gpuSelectMode         "prompt"
    gpuSelectDevices      "all"
    gpuAutoSelectMode     "availableMemory"
    gpuAutoSelectCount    "1"
    gpuFallBack           ""
  runSettingsGpuEnd
  quitOnPrompt            no
  quitOnLostLicense       no
  quitOnFE_INVALID        no
  tempLimitMode           ""
  chemistryBalanceMode    ""
  minDiskPercentMode      ""
  minDiskSpacePercentage  1
  advancedOptions         ""
runSettingsEnd

# Advanced GPU Settings
particleLoadBalanceBegin
  enabled ON
  decompositionType autoDecomp
  manualSettings OFF
  threshold 0.900000
  checkFreq 100
  dFactor 0.020000
  minCount 2000000
particleLoadBalanceEnd

guiFluidLibPath "C:\Program Files\CPFD\Barracuda\25.0.0.beta\props\cpfd_prop.prp"

inviscid 0
inviscidWall 0

# Particle close pack vol frac.
voidCP 0.600000

# Local close pack vol frac. Vol frac are fraction of voidCP 
#        on/off ResetVolFrac MinVolFrac  MaxVolFrac      C1        MinImpact
localVoidCP  0   0.8000000   1.0000000   1.1000000  1.000000e+03  -1.000000e+00 

# Gravity (body force)
gravityBegin
  type   CONSTANT   # CONSTANT, ROTATING, or SFF

# If type = CONSTANT or ROTATING, these values are used
  xGrav  0.000000e+00
  yGrav  0.000000e+00
  zGrav -9.810000e+00

# If type = ROTATING, the following information is used
  rotationAxis X
  rotationVel  0.000000e+00

# If type = SFF, the following file is used
  gravityFile ""
gravityEnd

# Particle-to-wall bounce
wallReflNRetain    0.300000
wallReflTRetain    0.990000
wallReflSpecular   5
wallImpactModels 1
  wallImpactModel 0
    type             legacy   # legacy, accretion, oka
    enabled          Off   # on, off
    startTime        0.000000e+00
    location
        x1                [
        x2                ]
        y1                [
        y2                ]
        z1                [
        z2                ]
        direction         normal
        normalLimit       normal
        vectorDir         0.00000e+00  0.00000e+00  1.00000e+00
        forceAbsoluteDir  0
        varAngleVec       1.50000e+01
        useExpression     0
    end_location
    wallMaterial     user_defined
    particleMaterial user_defined
    modelParameters
      impactAngles 11
        0.0               0.000000e+00
        0.1               7.100000e-01
        0.2               1.000000e+00
        0.3               9.700000e-01
        0.4               9.000000e-01
        0.5               8.000000e-01
        0.6               7.200000e-01
        0.7               6.100000e-01
        0.8               5.200000e-01
        0.9               3.800000e-01
        1.0               3.000000e-01
      endImpactAngles
      coefK             1.000000e+00
      powMass           1.500000e+00
      powVel            3.500000e+00
      minImpact  0.000000e+00
    endModelParameters
    comment          ""
  endWallImpactModel
endWallImpactModels

# Particle-to-wall impact
#          attritionFlag  powVel     powMass     minTrauma    startTime 
wallTrauma     0          2          1                 0      0.000000e+00 

# Particle-to-particle restitution
p_pRestitution   4.000000e-01
modelRedirection 0

# Small control volume control
smallCv          0.000000e+00

# particle volume interpolation at cut cells on a PresBC
usePolatePBcFlag on

# Boundary Conditions Properties
typeParticleFeed                       1
particleOverfeed                       1.00000e+00
particleFeedDeficit                    1.00000e-01
particleFeedDeficitFreq               -1.00000e+01
particleExitSmoothing                  0


# Pressure BC
presBC 1

presBCdescription  0
  fileFlag           1
  dir                2
  x1                 -5.000000e-02
  x2                 5.000000e-02
  y1                 -4.996950e-02
  y2                 4.996950e-02
  z1                 9.905660e-01
  z2                 1.000000e+00
  useExpression     0
  areaFrac           1.00000e+00
  kFact              0.00000e+00
  limit              2
  pBcBehavior        EXIT_UNRESTRICT
  calcHydrostatic    Off
  pres               0.00000e+00
  temp               3.00000e+02
  file               "reactor_outlet.sff"
  maxExitRad         1.00000e+10
  minExitRad         0.00000e+00
  units micron-diameter
  comment            "TOP PRESSURE BC"
  fluxPlaneName      "FLUX_REACTOR_OUT"
  showRad             1
  nFluxRad            5
  showRadFile         ""
  gasSpcFluxPlane     MASS_FRACTION
  rawParticleData     on
  rawTracerData       off
  typeGasFrac  MASS_FRACTION
  fluidSpecies
    gasSpc    "H2O" 0.750000
    gasSpc    "CH4_2" 0.250000
  endFluidSpecies
  tracerFeed
    enabled           OFF
    species           2000
    tracerResolution  1.000000e+02
    useBcConnector    OFF
    useSffFile        OFF
    sffFile           ""
  endTracerFeed
  usePipeFeedback   Off
end


# Inlet flow BC
flowBC 1

flowBCdescription  0
  fileFlag          1
  fileFlagParticle  0
  typeFlowBc        VELOCITY
  typeDirFlow       zDir
  normalLimit       normal
  vectorDir         0.00000e+00  0.00000e+00  1.00000e+00
  forceAbsoluteDir  0
  varAngleVec       1.50000e+01
  useExpression     0
  x1                -5.000000e-02
  x2                5.000000e-02
  y1                -4.996950e-02
  y2                4.996950e-02
  z1                0.000000e+00
  z2                9.433960e-03
  areaFrac          1.00000e+00
  pBcBehavior       BLOCKED
  flowBc            0.00000e+00
  temp              3.00000e+02
  pres              0.00000e+00
  file              "reactor_inlet_fluidization.sff"
  fileParticle      "reactor_inlet_fluidization.sff"
  maxExitRad        1.00000e+10
  minExitRad        0.00000e+00
  units micron-diameter
  comment           ""
  fluxPlaneName     "FLUX_REACTOR_INLET"
  showRad             0
  nFluxRad            100
  showRadFile         ""
  gasSpcFluxPlane     MASS_FRACTION
  rawParticleData     off
  rawTracerData       off
  typeGasFrac  MOLE_FRACTION
  fluidSpecies
    gasSpc    "H2O" 0.750000
    gasSpc    "CH4_2" 0.250000
  endFluidSpecies
  particleExit
    onOff                   OFF
    typeRateControl         RATE_METERED
    particleFlowRate        0
    particleVolumeFraction  0
  endParticleExit
  particleFeed
    particleResolutionFlag  off

    # Particle feed cloud resolution
    cloudResolution
      type              NP_PER_CELL
      sliderLevel       3
      npPerCell         1.250000e+02
      npPerMass         5.000000e+03
    end_cloudResolution

    particleResolutionFile  ""
    typeMixtureFrac         massFrac
    particleSpecies
    endParticleSpecies
    flux                    0.00000e+00
    flowRate                0.00000e+00
    pVolFrac                3.00000e-01
    slip                    5.00000e-01
  endParticleFeed
  particleFeedControl
    mode             MODE_NONE
    control          CTRL_NONE
    i                -1
    j                -1
    k                -1
    fluxPlaneName    ""
    fluxPlaneVarType CTRL_FPVAR_NONE
    test1            TST_NOTEST
      testVal1         0.00000e+00
      action1          ACT_NOACTION
      fluxPlaneName1   ""
      actionVal11      1.00000e+00
      actionVal12      1.00000e+00
      dtimeRange1      0.00000e+00
    test2            TST_NOTEST
      testVal2         0.00000e+00
      action2          ACT_NOACTION
      fluxPlaneName2   ""
      actionVal21      1.00000e+00
      actionVal22      1.00000e+00
      dtimeRange2      0.00000e+00
  endFeedControl
  tracerFeed
    enabled           OFF
    species           2000
    tracerResolution  1.000000e+02
    useBcConnector    OFF
    useSffFile        OFF
    sffFile           ""
  endTracerFeed
end


# Normal stress
normStres   0            1             3          1e-08 

# Initial conditions ... must be entered after fluid prop.
initCond 1

ICdescription
  enabled     ON
  type        SPECIFIED
  x1          -5.000000e-02
  x2          5.000000e-02
  y1          -4.996950e-02
  y2          4.996950e-02
  z1          0.000000e+00
  z2          1.000000e+00
  u            0.00000e+00
  v            0.00000e+00
  w            0.00000e+00
  pres         1.01325e+05
  temp         3.00000e+02
  comment      ""
  typeGasFrac  MASS_FRACTION
  fluidSpecies
    gasSpc    "H2O" 0.750000
    gasSpc    "CH4_2" 0.250000
  endFluidSpecies
  fluidICfile     ""
  neighborCount   8
  distanceFactor  1.000000e+00
  calcHydrostatic Off
end


# Initial particle velocity.
InitVelP 0

# Global cloud resolution
cloudResolution
  type              SLIDER
  sliderLevel       3
  npPerCell         1.250000e+02
  npPerMass         5.000000e+03
end_cloudResolution


# Initial particle distribution
Distribution 2

distribution 0
  enabled                 ON
  species                 1
  x1                      -5.000000e-02
  x2                      5.000000e-02
  y1                      -4.996950e-02
  y2                      4.996950e-02
  z1                      0.000000e+00
  z2                      3.000000e-01
  volFrac                 1.617000e-01
  mass                    2.200000e+00
  temperature             3.000000e+02
  noMomentum              off
  type                    REGION
  typeVolume              VOLUME_MASS

  # Particle cloud resolution
  cloudResolution
    type              LEGACY_SLIDER
    sliderLevel       3
    npPerCell         1.250000e+02
    npPerMass         5.000000e+03
    npPerRegion       1.000000e+04
    legacySlider      3.000000e+00
    legacyNpPerCell   1.250000e+02
  end_cloudResolution

  randomize               1
  dVolFracRand            1.000000e+00
  comment                 "Initial Mass of Nickel-Catalyst in the reactor"
  particleResolutionFlag  off
  particleResolutionFile  ""
  particleICfile          ""
  resetResidence          OFF
  splitByGridRatio        OFF
endDistribution

distribution 1
  enabled                 ON
  species                 2
  x1                      -5.000000e-02
  x2                      5.000000e-02
  y1                      -4.996950e-02
  y2                      4.996950e-02
  z1                      0.000000e+00
  z2                      3.000000e-01
  volFrac                 4.382000e-01
  mass                    9.000000e-01
  temperature             3.000000e+02
  noMomentum              off
  type                    REGION
  typeVolume              VOLUME_MASS

  # Particle cloud resolution
  cloudResolution
    type              LEGACY_SLIDER
    sliderLevel       3
    npPerCell         1.250000e+02
    npPerMass         5.000000e+03
    npPerRegion       1.000000e+04
    legacySlider      3.000000e+00
    legacyNpPerCell   1.250000e+02
  end_cloudResolution

  randomize               1
  dVolFracRand            1.000000e+00
  comment                 "Initial Mass of Dolomite in the reactor"
  particleResolutionFlag  off
  particleResolutionFile  ""
  particleICfile          ""
  resetResidence          OFF
  splitByGridRatio        OFF
endDistribution


# Data plane variables. Data written each plot freq time
# gridResolution 0 or 1. 0 = specified resolution, 1 = use grid
dataPlanes 0

   timeInterval 1.000000e-01
   nonCellValue NAN
   useOldFileNames OFF

end_dataPlanes


# Minimum bc area reduction for aFrac (volFrac*aFrac*area >= areaFracMin*aFrac*area
areaFracMin 1.000000
#

# Cell variables
draw_f  prho ijk volume fCompMassConc fCompMassFrac fCompMoleConc fCompMoleFrac frho fTemp velf void pres

# Particle mesh variables
draw_p  pid density mass size speed temp velocity void npCloud species

# Save average data
SaveAveData start 0.000000e+00 
SaveAveDataByVolFrac  on  # default is on, other option is off

# Raw Data Output
rawDataOutput
  timeInterval 1.000000e+20
  particleDataEnabled off
  particleData  pid
  cellDataEnabled off
  cellData  xyz cellId cellVol
endRawDataOutput

# Population Data
population
  enabled               ON
  outputFrequency       1.000000e+02
  radiusDivisions       35
  typeRadiusDivision    TYPE_MINMAX_LEGACY
  radiusMin             0.000000e+00
  radiusMax             0.000000e+00
  units micron-diameter
endPopulation


# Write flux plane
fluxPlane 0

# Uniform cloud volume
uniformCloudVol 1

# Collision model
collisionRestitCoef       9.800000e-01
collisionRelaxEquil       off
collisionRelaxIsotropy    off
collisionLiqMassTransfer  on
liquidCollisionParams
    kus 1.000000e-02
    kul 1.000000e-01
    R 1.000000e+00
endLiquidCollisionParams
#Collision Blend Acceleration model
collisionBlendAccel off  1.000000e-01  6.000000e-01

# Turbulence model
turbulence
  turbModel           2
  SmagorinskyConst    1.000000e-02
  keProductionConst   9.000000e-02
  keDissipationConst  0.000000e+00
  thermalModel        0
  particleInTurbEddy  off
endTurbulence

# Blend 1st and 2nd order advection
advection
  type       QSOU
  alphaQSOU  0.300000
  betaQSOU   1.000000
  blendDSOU  0.000000
  userWarned 
end

# Blend explicit/implicit advection, diffusion and lumped resistance. 50=Crank-Nicolson
implicitWt 100

# Blend explicit/implicit Up in continuum drag
# (Note:  Not currently used)

# k-factor adjustment 
kFactAdjust 0.000000

# Close pack limit coefficient for adjusting initial volume fractions
volFracCoefToClosePack 1.020000

# Set initial seed with or without random variation 
randLoc 0

# damp particle normal stress induced velocity 
wacker 1

# ------------ Energy equation parameters ------------
thermal

  # cap total particle projected area in cells for radiation P1 model
  capParticleAreaRadiateP1    on

  interpParticleAreaRadiateP1    on


  # Thermal wall BCs
  thermalWallBc  0

  endThermalWallBc


  # Heat transfer coefficient constants

  # Fluid-to-wall heat transfer coefficient: c0 c1 c2 n1 n2
  htc_f2w_c0 4.600000e-01
  htc_f2w_c1 3.660000e+00
  htc_f2w_c2 0.000000e+00
  htc_f2w_n1 5.000000e-01
  htc_f2w_n2 3.300000e-01

  # Dense phase heat transfer coefficient: c0 n1
  htc_f2w_dense_c0  5.250000e-01
  htc_f2w_dense_n1  7.500000e-01

  # Fluid-to-particle heat transfer coefficient: c0 c1 c2 n1
  htc_f2p_c0 3.700000e-01
  htc_f2p_c1 1.000000e-01
  htc_f2p_c2 0.000000e+00
  htc_f2p_n1 6.000000e-01

endthermal
# ---------- End energy equation parameters ----------

# wall heat output format (controls all output): RATE (J/s) or FLUX (J/sm^2)
wallHeatOutputFormat  RATE
multBaseVel   1.000000e+00


# Chemistry module section
chemistry

  # Volume-Average reaction type, stored user setting, 0=Stoich, 1=Species
  reactionTypeVolAve   0

  # kg/m3, g-mol/m3, or kg-mol/m3
  concentrationType    kg/m3

  # volumeFraction or surfaceArea
  solidConcentrationType  volumeFraction

  # reaction boundary layer, on/off
  reactionBoundaryLayer  off

  # site activation is a multiplier times solid surface area
  siteActivationCoef     1

  heatDistribution       AUTO
  wtSolidGasHeatReact    1.000000e+00
  heatPartitionVolAve    MANUAL
  wtSolidHeatRxVolAve    0.000000e+00
  relativeToleranceOde   1.000000e-09
  absoluteToleranceOde   1.000000e-09
  maxNumStepsOde         200000
  maxNumIterOde          4
  # Minimum value used in the ODE solver
  minValueConcenOde      1.000000e-09
  resetMinValueThreshold 0
  negConcRecErr          0
  maxOrderBdf            5
  maxNonlinIters         3
  enableVolumeRatioLimit 0
  volumeRatioLimit       1.000000e+00

  chemicalRateCoef 11
    chemRate
      id           0
      typeRate     0
      wT           5.000000e-01
      c0           5.117000e-08
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG   -1.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E            2.887900e+04    1.000000e+00    0.000000e+00
      E0           4.456700e+01
      solidSpc     "Ni_1"
        type       mass
        pow        1.000000e+00    1.000000e+00    0.000000e+00
      endSolidSpc
      unitsSolidMassPerVol  "kg/m^3"
      comment               "k1 (kin const) - (kmol*bar^0.5)/(m^3*sec)"
    end

    chemRate
      id           1
      typeRate     0
      wT           5.000000e-01
      c0           6.091700e-09
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG   -1.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E            2.933610e+04    1.000000e+00    0.000000e+00
      E0           4.527000e+01
      solidSpc     "Ni_1"
        type       mass
        pow        1.000000e+00    1.000000e+00    0.000000e+00
      endSolidSpc
      unitsSolidMassPerVol  "kg/m^3"
      comment               "k2 (kin const) - (kmol*bar^0.5)/(m^3*sec)"
    end

    chemRate
      id           2
      typeRate     0
      wT           5.000000e-01
      c0           2.094440e-03
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG   -1.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E            8.074330e+03    1.000000e+00    0.000000e+00
      E0           1.246040e+01
      solidSpc     "Ni_1"
        type       mass
        pow        1.000000e+00    1.000000e+00    0.000000e+00
      endSolidSpc
      unitsSolidMassPerVol  "kg/m^3"
      comment               "k3 (kin const) - (kmol)/(m^3*sec*bar)"
    end

    chemRate
      id           3
      typeRate     0
      wT           5.000000e-01
      c0           4.707000e+12
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG    0.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E            2.694251e+04    1.000000e+00    0.000000e+00
      E0           0.000000e+00
      comment               "K_I (eq const) - bar^2"
    end

    chemRate
      id           4
      typeRate     0
      wT           5.000000e-01
      c0           5.375394e+10
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG    0.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E            2.245610e+04    1.000000e+00    0.000000e+00
      E0           0.000000e+00
      comment               "K_II (eq const) - bar^2"
    end

    chemRate
      id           5
      typeRate     0
      wT           5.000000e-01
      c0           1.142000e-02
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG    0.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E           -4.486410e+03    1.000000e+00    0.000000e+00
      E0           0.000000e+00
      comment               "K_III (eq const) - unitless"
    end

    chemRate
      id           6
      typeRate     0
      wT           5.000000e-01
      c0           1.790000e-01
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG    0.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E           -4.604280e+03    1.000000e+00    0.000000e+00
      E0          -5.594500e+00
      comment               "K_CH4 (ad const) - bar^-1"
    end

    chemRate
      id           7
      typeRate     0
      wT           5.000000e-01
      c0           4.152000e-01
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG    0.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E            1.066634e+04    1.000000e+00    0.000000e+00
      E0           1.296030e+01
      comment               "K_H2O (ad const) - unitless"
    end

    chemRate
      id           8
      typeRate     0
      wT           5.000000e-01
      c0           4.091000e+01
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG    0.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E           -8.497710e+03    1.000000e+00    0.000000e+00
      E0          -1.311370e+01
      comment               "K_CO (ad const) - bar^-1"
    end

    chemRate
      id           9
      typeRate     0
      wT           5.000000e-01
      c0           2.960000e-02
      nTemp        0.000000e+00    1.000000e+00    0.000000e+00
      nPres        0.000000e+00    1.000000e+00    0.000000e+00
      nRhoF        0.000000e+00    1.000000e+00    0.000000e+00
      nVolFracG    0.000000e+00    1.000000e+00    0.000000e+00
      nNpPerVol    0.000000e+00    1.000000e+00    0.000000e+00
      nReynolds    0.000000e+00    1.000000e+00    0.000000e+00
      E           -9.971130e+03    1.000000e+00    0.000000e+00
      E0          -1.538700e+01
      comment               "K_H2 (ad const) - bar^-1"
    end

    rateCoefExpression
      id           10
      typeRate     1
      expression
        // User-defined rate coefficient expression generated from Arrhenius rate coefficient
        MW_CaO = 56.0774; //kg/kmol
        MW_MgO= 40.3044; //kg/kmol
        S_0 = 31*10^3; //initial specific area of CaO in m^2/kg_sor
        p_CO2_eq_lowT = 4.1918*10^12*exp(-20474/tempF); //Pa (for T<=1,173.15K)
        p_CO2_eq_highT = 1.216*10^12*exp(-19130/tempF); //Pa (for T>1,173.15K)
        k_carb_lowP = 1.04*10^-10*exp(-2453.69/tempF);//(kmol/(m^-2*pa^-1*s^-1)) f(or P <=10,000pa)
        k_carb_highP = 1.04*10^-6*exp(-2453.69/tempF);//(kmol/(m^-2*s^-1)) f(or P >10,000pa)
        n_low = 1;
        n_high = 0;
        P_CO2_eq = if(tempF<=1173.15,p_CO2_eq_lowT,p_CO2_eq_highT);
        k_carb = if((pressure*nf_CO2-P_CO2_eq)<=10000,k_carb_lowP,k_carb_highP); 
        n = if((pressure*nf_CO2-P_CO2_eq)<=10000,n_low,n_high); 
        r_carb = -k_carb*(pressure*nf_CO2-P_CO2_eq)^n*S_0*mass_CaO; //kmol*kg_sor^-1*sec^-1*kg_sor
        return r_carb;
      endExpression
      comment               "Discrete Carbonation reaction"
    end

  endChemicalCoef


  chemicalReactionSolid 1

    chemEquation 1
      unitsRate               "kmol/s"
      unitsConcentration      "kg/m^3"
      shrinkingCore
        state                 off
        diffusionCoef         1.000000e-04
        diffusionUnits        "m^2/s"
      endShrinkingCore
      solidSpc "CaO"
      comment                 ""
      rate  1
        superGroup 1
          operator     "+"
          power        1.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId     (1.00000e+00 k10) 
            materialSpecies
            endMaterialSpecies
          endChemGroup
        endSuperGroup
      endRate
      product
      gasSpc "CO2"  (1.00000e+00)
      solidSpc "CaCO3"  (-1.00000e+00)
      endProduct
      enforceMaxRate          off
      maxRate                 1.000000e+30
      enforceMinRate          off
      minRate                 -1.000000e+30
    endChemEquation

  endChemicalReaction

  stoichiometric

    stoichiometricEquations
      00:  CH4_2 + H2O => CO + 3.000000 H2
      01:  CO + 3.000000 H2 => CH4_2 + H2O
      02:  CH4_2 + 2.000000 H2O => CO2 + 4.000000 H2
      03:  CO2 + 4.000000 H2 => CH4_2 + 2.000000 H2O
      04:  CO + H2O => CO2 + H2
      05:  CO2 + H2 => CO + H2O
    endStoichiometricEquations

    stoichiometricRates

    rateEqu  00
      unitsRate               "kmol/s*m^3"
      unitsConcentration      "bar"
      comment                 "R1 - forward reaction"
      rate  3
        superGroup 1
          operator     "+"
          power        1.00000e+00
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k0) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
              gasSpc "H2O"   1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 1
            const      1.00000e+00
            coefId    
            materialSpecies
              gasSpc "H2"   2.50000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 5
          operator     "*"
          power       -2.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId    
            materialSpecies
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k8) 
            materialSpecies
              gasSpc "CO"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k9) 
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k6) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k7) 
            materialSpecies
              gasSpc "H2O"   1.00000e+00
              gasSpc "H2"  -1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
      endRate
    endRateEqu

    rateEqu  01
      unitsRate               "kmol/s*m^3"
      unitsConcentration      "bar"
      comment                 "R1 - backward reaction"
      rate  4
        superGroup 1
          operator     "+"
          power        1.00000e+00
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k0) 
            materialSpecies
              gasSpc "H2"   3.00000e+00
              gasSpc "CO"   1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 1
            const      1.00000e+00
            coefId    
            materialSpecies
              gasSpc "H2"   2.50000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 5
          operator     "*"
          power       -2.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId    
            materialSpecies
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k8) 
            materialSpecies
              gasSpc "CO"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k9) 
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k6) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k7) 
            materialSpecies
              gasSpc "H2O"   1.00000e+00
              gasSpc "H2"  -1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId     (1.00000e+00 k3) 
            materialSpecies
            endMaterialSpecies
          endChemGroup
        endSuperGroup
      endRate
    endRateEqu

    rateEqu  02
      unitsRate               "kmol/s*m^3"
      unitsConcentration      "bar"
      comment                 "R2 - forward reaction"
      rate  3
        superGroup 1
          operator     "+"
          power        1.00000e+00
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k1) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
              gasSpc "H2O"   2.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 1
            const      1.00000e+00
            coefId    
            materialSpecies
              gasSpc "H2"   3.50000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 5
          operator     "*"
          power       -2.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId    
            materialSpecies
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k8) 
            materialSpecies
              gasSpc "CO"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k9) 
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k6) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k7) 
            materialSpecies
              gasSpc "H2O"   1.00000e+00
              gasSpc "H2"  -1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
      endRate
    endRateEqu

    rateEqu  03
      unitsRate               "kmol/s*m^3"
      unitsConcentration      "bar"
      comment                 "R2 - backward reaction"
      rate  4
        superGroup 1
          operator     "+"
          power        1.00000e+00
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k1) 
            materialSpecies
              gasSpc "H2"   4.00000e+00
              gasSpc "CO2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 1
            const      1.00000e+00
            coefId    
            materialSpecies
              gasSpc "H2"   3.50000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 5
          operator     "*"
          power       -2.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId    
            materialSpecies
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k8) 
            materialSpecies
              gasSpc "CO"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k9) 
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k6) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k7) 
            materialSpecies
              gasSpc "H2O"   1.00000e+00
              gasSpc "H2"  -1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId     (1.00000e+00 k4) 
            materialSpecies
            endMaterialSpecies
          endChemGroup
        endSuperGroup
      endRate
    endRateEqu

    rateEqu  04
      unitsRate               "kmol/s*m^3"
      unitsConcentration      "bar"
      comment                 "R3 - forward reaction"
      rate  3
        superGroup 1
          operator     "+"
          power        1.00000e+00
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k2) 
            materialSpecies
              gasSpc "CO"   1.00000e+00
              gasSpc "H2O"   1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 1
            const      1.00000e+00
            coefId    
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 5
          operator     "*"
          power       -2.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId    
            materialSpecies
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k8) 
            materialSpecies
              gasSpc "CO"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k9) 
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k6) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k7) 
            materialSpecies
              gasSpc "H2O"   1.00000e+00
              gasSpc "H2"  -1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
      endRate
    endRateEqu

    rateEqu  05
      unitsRate               "kmol/s*m^3"
      unitsConcentration      "bar"
      comment                 "R3 - backward reaction"
      rate  4
        superGroup 1
          operator     "+"
          power        1.00000e+00
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k2) 
            materialSpecies
              gasSpc "H2"   1.00000e+00
              gasSpc "CO2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 1
            const      1.00000e+00
            coefId    
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 5
          operator     "*"
          power       -2.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId    
            materialSpecies
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k8) 
            materialSpecies
              gasSpc "CO"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k9) 
            materialSpecies
              gasSpc "H2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 1
            const      1.00000e+00
            coefId     (1.00000e+00 k6) 
            materialSpecies
              gasSpc "CH4_2"   1.00000e+00
            endMaterialSpecies
          endChemGroup
          chemGroup 2
            const      1.00000e+00
            coefId     (1.00000e+00 k7) 
            materialSpecies
              gasSpc "H2O"   1.00000e+00
              gasSpc "H2"  -1.00000e+00
            endMaterialSpecies
          endChemGroup
        endSuperGroup
        superGroup 1
          operator     "*"
          power       -1.00000e+00
          chemGroup 0
            const      1.00000e+00
            coefId     (1.00000e+00 k5) 
            materialSpecies
            endMaterialSpecies
          endChemGroup
        endSuperGroup
      endRate
    endRateEqu

    endStoichiometricRates

  endStoichiometric

endChemistry

# transient plot variables 
transientPlot

wallBlockHeatTransientOut  1

# nXyz write flush      Output file    
# ---- ----- ----- --------------------
  2   1   0    trans.data  

# 0=loc 1=node 

# type                    loc        x             y             z        comment
# ----------------------- --- ------------- ------------- ------------- ------------
 pres                      0  0.000000e+00  0.000000e+00  5.000000e-03  "bed bottom (z = 0.005m)"
 pres                      0  0.000000e+00  0.000000e+00  8.000000e-01  "above expanded bed (z = 0.8m)"

end_transientPlot
# End transient plot variables 


# --- development --- features in solver not yet in gui
developmentBegin
developmentEnd
# --------------- End of development ---------------

# End
