1# 2# Run with ./extchem -options_file gasoline.inp 3# 4# The two input files come from OpenFoam tutorial for their 0-D chemistry solver chemFoam 5# 6# The file names are too generic for this example so make sure you don't have other files with these names already in the directory 7# 8# You may need to do 9# cp $PETSC_DIR/$PETSC_ARCH/share/periodictable.dat . 10# 11# 12-chem https://raw.githubusercontent.com/OpenFOAM/OpenFOAM-dev/master/tutorials/combustion/chemFoam/ic8h18/chemkin/chem.inp 13-thermo https://raw.githubusercontent.com/OpenFOAM/OpenFOAM-dev/master/tutorials/combustion/chemFoam/ic8h18/chemkin/therm.dat 14 15-initial_species IC8H18,O2,N2 16#initial concentration in moles 17-initial_mole 1,12.5,47.0 18 19#Temp in kelvin 20-Tini 800 21##pressure is in N/m2 22-pressure 5.06625e+06 23 24-ts_max_time 6.5e-3 25-ts_max_steps 100000000 26 27-monitor_temp 28 29# glee insists on very small timesteps so never makes any progress 30#-ts_adapt_type glee 31-ts_arkimex_fully_implicit 32-ts_max_snes_failures -1 33-ts_adapt_monitor 34-ts_adapt_dt_max 1e-3 35-ts_arkimex_type 4 36 37# the chemistry dictates reasonable bounds on the timestep 38-ts_adapt_dt_max 1e-4 39-ts_adapt_dt_min 1e-9 40 41-snes_max_it 10 42 43#-snes_stol 1.e-6 44 45# use direct solver (treats Jacobian as dense) 46-pc_type lu 47 48#-ts_monitor_lg_solution 49#-draw_pause -2 50-lg_use_markers false 51-ts_monitor_lg_solution_variables H2,O2,H2O,CH4,CO,CO2,C2H2,N2 52#-ts_monitor_envelope 53