Hi all,

I have troubles modeling a non-premixed laminar axissymetric co-flow methane/air combustion. During the simulation the temperature is getting too high, order 3300 K. Adiabatic methane/air temperature is 2279 K. Could anyone give me some advice?

Thanks in advance, Peter

P.s: the inlet profile and journals below arent really readable, but I can post them again. Think I reached maximum chars or something.

Im using the segregated solver, and also tried the coupled now. The first bit described the segregated solver set up.

- Species model uses the Stiff-chemistry, inlet Diffusion, Diffusion Energy Source, Full multicomponent diffusion, thermal diffusion

- I am importing a chemkin DRM19 mechanism. As a consequence my material properties are: in FLUID o Species Cp are piecewise-polynomial o Thermal conductivity & viscosity kinetic theory In Mixture: o Density is volume weighted mixing law o Thermal conductivity and viscosity is mass weighted mixing law. Also tried the ideal-gas-mixing-law here. Both gave the same issues.

- Boundary conditions at the inlet can be found below. The inlet of the methane in the co-flow is 4e-3 (2e-3 axissymetric) while the rest is the air co-flow to 5e-2

- Starting off with first order upwind and going to second order. Under relaxation of pressure and temperature to 0.9. While species to 0.8

- Starting up with just the flow. During these few iterations I use a pressure-outlet of -20 gage pressure to start up the flow. Changing this to outflow with the cold flow

- Then starting up with cold flow (so incl. species) with no volumetric reactions

- Patching a hot region of 2300 K to ignite the flow and turning on the volumetric reactions. Using the ISAT table I decrease error tolerance). While doing iteration I clean and create a new ISAT table now and then.

- At the end also doing direct integration.

- Mesh size is good. Looking at OH,CH3 species and its source give smooth results.

Coupled server;

The same approach as the segregated server. Difference is that during iteration I go from courant number of 0.5 to 5 and the high temperature patching is more slowly.

Under here one can find my boundary conditions and my journal files to illustrate what I am using.

Segregated journal: file read-case /home/s021018/Diffusioncase_full/100/CH4_grav_2orde_100_lenonunity_mesh007035.cas read-data /home/s021018/Diffusioncase_full/100/CH4_grav_2orde_100_lenonunity_mesh007035.dat /file/read-profile/ /home/s021018/Diffusioncase_full/100/inlet.prof /file/import/chemkin-mechanism "chemkin-import" "/home/s021018/Diffusioncase_full/drm19.mech" "/home/s021018/Diffusioncase_full/60/thermo211.dat" no yes "/home/s021018/Diffusioncase_full/transport.dat" /define/materials/change-create chemkin-import chemkin-import no no no yes volume-weighted-mixing-law no yes ideal-gas-mixing-law yes ideal-gas-mixing-law no no no no /define/operating-conditions/gravity yes -9.81 0 /solve/set/limits 10132 5e+10 280 3800 /solve/monitors/residual/print yes /define/models/species/volumetric-reactions no /solve/set/equations/flow yes /solve/set/equations/temperature no /solve/set/equations/species-0 no /solve/set/equations/species-1 no /solve/set/equations/species-2 no /solve/set/equations/species-3 no /solve/set/equations/species-4 no /solve/set/equations/species-5 no /solve/set/equations/species-6 no /solve/set/equations/species-7 no /solve/set/equations/species-8 no /solve/set/equations/species-9 no /solve/set/equations/species-10 no /solve/set/equations/species-11 no /solve/set/equations/species-12 no /solve/set/equations/species-13 no /solve/set/equations/species-14 no /solve/set/equations/species-15 no /solve/set/equations/species-16 no /solve/set/equations/species-17 no /solve/set/equations/species-18 no /solve/set/equations/species-19 no /solve/iterate/ 400 /solve/set/equations/temperature yes /solve/set/equations/species-0 yes /solve/set/equations/species-1 yes /solve/set/equations/species-2 yes /solve/set/equations/species-3 yes /solve/set/equations/species-4 yes /solve/set/equations/species-5 yes /solve/set/equations/species-6 yes /solve/set/equations/species-7 yes /solve/set/equations/species-8 yes /solve/set/equations/species-9 yes /solve/set/equations/species-10 yes /solve/set/equations/species-11 yes /solve/set/equations/species-12 yes /solve/set/equations/species-13 yes /solve/set/equations/species-14 yes /solve/set/equations/species-15 yes /solve/set/equations/species-16 yes /solve/set/equations/species-17 yes /solve/set/equations/species-18 yes /solve/set/equations/species-19 yes /solve/iterate/ 600 /define/boundary-conditions/zone-type outflow outflow /solve/iterate/ 600 /define/boundary-conditions/zone-type radius outflow /solve/iterate/ 600 /file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_initialfl ow.cas /define/models/species/volumetric-reactions yes /adapt/mark-inout-rectangle yes no 0 0.01 0 0.005 /solve/patch () hexahedron-r0 () temperature 2300 /define/models/species/clear-isat-table /define/models/species/integration-parameters yes 1e-08 1e-09 0.001 100 0.001 1 0 /solve/iterate/ 3000 /define/models/species/clear-isat-table /define/models/species/integration-parameters yes 1e-08 1e-09 0.0001 100 0.001 1 0 /solve/iterate/ 3000 /solve/set/discretization-scheme/mom 1 /solve/set/discretization-scheme/pressure 12 /solve/set/discretization-scheme/species-0 1 /solve/set/discretization-scheme/species-1 1 /solve/set/discretization-scheme/species-2 1 /solve/set/discretization-scheme/species-3 1 /solve/set/discretization-scheme/species-4 1 /solve/set/discretization-scheme/species-5 1 /solve/set/discretization-scheme/species-6 1 /solve/set/discretization-scheme/species-7 1 /solve/set/discretization-scheme/species-8 1 /solve/set/discretization-scheme/species-9 1 /solve/set/discretization-scheme/species-10 1 /solve/set/discretization-scheme/species-11 1 /solve/set/discretization-scheme/species-12 1 /solve/set/discretization-scheme/species-13 1 /solve/set/discretization-scheme/species-14 1 /solve/set/discretization-scheme/species-15 1 /solve/set/discretization-scheme/species-16 1 /solve/set/discretization-scheme/species-17 1 /solve/set/discretization-scheme/species-18 1 /solve/set/discretization-scheme/species-19 1 /solve/set/discretization-scheme/temperature 1 /solve/set/p-v-coupling 21 /solve/iterate/ 500 /define/models/species/clear-isat-table /define/models/species/integration-parameters yes 1e-08 1e-09 0.0001 100 0.001 1 0 /solve/set/limits 10132 5e+10 280 3500 /solve/iterate/ 2000 /file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_ISAT.cas /define/models/species/clear-isat-table /define/models/species/integration-parameters no 1e-08 1e-09 /solve/iterate/ 500 /file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_DIRECT.ca s q q exit yes

Coupled journal: file read-case /home/s021018/Diffusioncase_full/100/mesh00502100400.cas read-data /home/s021018/Diffusioncase_full/100/mesh00502100400.dat /file/read-profile/ /home/s021018/Diffusioncase_full/100/inlet.prof /file/import/chemkin-mechanism "chemkin-import" "/home/s021018/Diffusioncase_full/drm19.mech" "/home/s021018/Diffusioncase_full/60/thermo211.dat" no yes "/home/s021018/Diffusioncase_full/transport.dat" /define/materials/change-create chemkin-import chemkin-import no no no yes volume-weighted-mixing-law no yes mass-weighted-mixing-law yes mass-weighted-mixing-law no no no no /define/operating-conditions/gravity yes -9.81 0 /solve/set/limits 10132 5e+10 285 3000 0.2 /solve/monitors/residual/print yes /define/models/species/volumetric-reactions no /solve/iterate/ 40 /define/boundary-conditions/zone-type exit outflow /solve/iterate/ 100 /file/write-case-data/100/mesh00502100400_begin.cas /define/models/species/volumetric-reactions yes /solve/set/stiff-chemistry yes 0.2 0.25 0.9 /adapt/mark-inout-rectangle yes no 0 0.01 0 0.005 /solve/patch () hexahedron-r0 () temperature 1200 /solve/set/courant-number 0.5 /solve/iterate/ 50 /adapt/mark-inout-rectangle yes no 0 0.01 0 0.005 /solve/patch () hexahedron-r0 () temperature 1600 /solve/set/courant-number 1 /solve/iterate/ 50 /adapt/mark-inout-rectangle yes no 0 0.01 0 0.005 /solve/patch () hexahedron-r0 () temperature 1800 /solve/set/courant-number 3 /solve/iterate/ 50 /adapt/mark-inout-rectangle yes no 0 0.01 0 0.005 /solve/patch () hexahedron-r0 () temperature 2000 /solve/set/courant-number 5 /solve/iterate/ 50 /adapt/mark-inout-rectangle yes no 0 0.01 0 0.005 /solve/patch () hexahedron-r0 () temperature 2279 /solve/iterate/ 6000 /file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_coupled_w all.cas /define/boundary-conditions/zone-type radius outflow /solve/iterate/ 6000 /file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_coupled_a lloutflow.cas q q exit yes

Im trying to simulate results from already done measurement. The inlet of the methane is 4e-3 wide while my co-flow and domain are 5e-2 wide. Boundary inlet conditions can be found below.

((inletprofile radial 29 ) (r 0.000000e+00 5.000000e-04 1.000000e-03 1.500000e-03 2.000000e-03 2.500000e-03 3.000000e-03 3.500000e-03 3.600000e-03 3.850000e-03 4.000000e-03 4.500000e-03 5.000000e-03 5.500000e-03 6.000000e-03 6.500000e-03 7.000000e-03 7.500002e-03 7.999999e-03 8.500003e-03 8.856000e-03 8.928000e-03 9.100000e-03 9.450000e-03 9.499997e-03 1.000000e-02 1.050000e-02 1.515000e-02 5.000000e-02 ) (vax 1.326291e-01 1.326291e-01 1.326291e-01 1.326291e-01 1.326291e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 ) (CH4 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 ) (O2 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 ) (N2 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 ) (CO2 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 ) (H2O 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 ) (T 2.940000e+02 2.940000e+02 2.940000e+02 2.940000e+02 2.940000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 ) )

The inlet profile in a better template

((inletprofile radial 29 )

(r

0.000000e+00 5.000000e-04 1.000000e-03 1.500000e-03 2.000000e-03 2.500000e-03 3.000000e-03 3.500000e-03 3.600000e-03 3.850000e-03 4.000000e-03 4.500000e-03 5.000000e-03 5.500000e-03 6.000000e-03 6.500000e-03 7.000000e-03 7.500002e-03 7.999999e-03 8.500003e-03 8.856000e-03 8.928000e-03 9.100000e-03 9.450000e-03 9.499997e-03 1.000000e-02 1.050000e-02 1.515000e-02 5.000000e-02 )

(vax

1.326291e-01 1.326291e-01 1.326291e-01 1.326291e-01 1.326291e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 1.537728e-01 )

(CH4

1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 )

(O2

0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 2.311300e-01 )

(N2

0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 7.625000e-01 )

(CO2

0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 5.500000e-04 )

(H2O

0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 5.810000e-03 )

(T

2.940000e+02 2.940000e+02 2.940000e+02 2.940000e+02 2.940000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 2.910000e+02 ) )

The journal file of the segregated solver in a better template

file

read-case

/home/s021018/Diffusioncase_full/100/CH4_grav_2orde_100_lenonunity_mesh007035.cas

read-data

/home/s021018/Diffusioncase_full/100/CH4_grav_2orde_100_lenonunity_mesh007035.dat

/file/read-profile/ /home/s021018/Diffusioncase_full/100/inlet.prof

/file/import/chemkin-mechanism "chemkin-import" "/home/s021018/Diffusioncase_full/drm19.mech" "/home/s021018/Diffusioncase_full/60/thermo211.dat"

no

yes "/home/s021018/Diffusioncase_full/transport.dat"

/define/materials/change-create chemkin-import

chemkin-import

no

no

no

yes volume-weighted-mixing-law

no

yes ideal-gas-mixing-law

yes ideal-gas-mixing-law

no

no

no

no

/define/operating-conditions/gravity yes -9.81 0

/solve/set/limits 10132 5e+10 280 3800

/solve/monitors/residual/print yes

/define/models/species/volumetric-reactions no

/solve/set/equations/flow yes

/solve/set/equations/temperature no

/solve/set/equations/species-0 no

/solve/set/equations/species-1 no

/solve/set/equations/species-2 no

/solve/set/equations/species-3 no

/solve/set/equations/species-4 no

/solve/set/equations/species-5 no

/solve/set/equations/species-6 no

/solve/set/equations/species-7 no

/solve/set/equations/species-8 no

/solve/set/equations/species-9 no

/solve/set/equations/species-10 no

/solve/set/equations/species-11 no

/solve/set/equations/species-12 no

/solve/set/equations/species-13 no

/solve/set/equations/species-14 no

/solve/set/equations/species-15 no

/solve/set/equations/species-16 no

/solve/set/equations/species-17 no

/solve/set/equations/species-18 no

/solve/set/equations/species-19 no

/solve/iterate/ 400

/solve/set/equations/temperature yes

/solve/set/equations/species-0 yes

/solve/set/equations/species-1 yes

/solve/set/equations/species-2 yes

/solve/set/equations/species-3 yes

/solve/set/equations/species-4 yes

/solve/set/equations/species-5 yes

/solve/set/equations/species-6 yes

/solve/set/equations/species-7 yes

/solve/set/equations/species-8 yes

/solve/set/equations/species-9 yes

/solve/set/equations/species-10 yes

/solve/set/equations/species-11 yes

/solve/set/equations/species-12 yes

/solve/set/equations/species-13 yes

/solve/set/equations/species-14 yes

/solve/set/equations/species-15 yes

/solve/set/equations/species-16 yes

/solve/set/equations/species-17 yes

/solve/set/equations/species-18 yes

/solve/set/equations/species-19 yes

/solve/iterate/ 600

/define/boundary-conditions/zone-type outflow outflow

/solve/iterate/ 600

/define/boundary-conditions/zone-type radius outflow

/solve/iterate/ 600

/file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_initialfl ow.cas

/define/models/species/volumetric-reactions yes

/adapt/mark-inout-rectangle yes no 0 0.01 0 0.005

/solve/patch () hexahedron-r0 () temperature 2300

/define/models/species/clear-isat-table

/define/models/species/integration-parameters yes 1e-08 1e-09 0.001 100 0.001 1 0

/solve/iterate/ 2000

/define/models/species/clear-isat-table

/define/models/species/integration-parameters yes 1e-08 1e-09 0.0001 100 0.001 1 0

/solve/iterate/ 2000

/solve/set/discretization-scheme/mom 1

/solve/set/discretization-scheme/pressure 12

/solve/set/discretization-scheme/species-0 1

/solve/set/discretization-scheme/species-1 1

/solve/set/discretization-scheme/species-2 1

/solve/set/discretization-scheme/species-3 1

/solve/set/discretization-scheme/species-4 1

/solve/set/discretization-scheme/species-5 1

/solve/set/discretization-scheme/species-6 1

/solve/set/discretization-scheme/species-7 1

/solve/set/discretization-scheme/species-8 1

/solve/set/discretization-scheme/species-9 1

/solve/set/discretization-scheme/species-10 1

/solve/set/discretization-scheme/species-11 1

/solve/set/discretization-scheme/species-12 1

/solve/set/discretization-scheme/species-13 1

/solve/set/discretization-scheme/species-14 1

/solve/set/discretization-scheme/species-15 1

/solve/set/discretization-scheme/species-16 1

/solve/set/discretization-scheme/species-17 1

/solve/set/discretization-scheme/species-18 1

/solve/set/discretization-scheme/species-19 1

/solve/set/discretization-scheme/temperature 1

/solve/set/p-v-coupling 21

/solve/iterate/ 500

/define/models/species/clear-isat-table

/define/models/species/integration-parameters yes 1e-08 1e-09 0.0001 100 0.001 1 0

/solve/set/limits 10132 5e+10 280 3500

/solve/iterate/ 2000

/file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_ISAT.cas

/define/models/species/clear-isat-table

/define/models/species/integration-parameters no 1e-08 1e-09

/solve/iterate/ 500

/file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_DIRECT.ca s

q

q

exit

yes

The journal file for the coupled solver in a better template

file

read-case

/home/s021018/Diffusioncase_full/100/mesh00502100400.cas

read-data

/home/s021018/Diffusioncase_full/100/mesh00502100400.dat

/file/read-profile/ /home/s021018/Diffusioncase_full/100/inlet.prof

/file/import/chemkin-mechanism "chemkin-import" "/home/s021018/Diffusioncase_full/drm19.mech" "/home/s021018/Diffusioncase_full/60/thermo211.dat"

no

yes "/home/s021018/Diffusioncase_full/transport.dat"

/define/materials/change-create chemkin-import

chemkin-import

no

no

no

yes volume-weighted-mixing-law

no

yes mass-weighted-mixing-law

yes mass-weighted-mixing-law

no

no

no

no

/define/operating-conditions/gravity yes -9.81 0

/solve/set/limits 10132 5e+10 285 3000 0.2

/solve/monitors/residual/print yes

/define/models/species/volumetric-reactions no

/solve/iterate/ 40

/define/boundary-conditions/zone-type exit outflow

/solve/iterate/ 100

/file/write-case-data/100/mesh00502100400_begin.cas

/define/models/species/volumetric-reactions yes

/solve/set/stiff-chemistry yes 0.2 0.25 0.9

/adapt/mark-inout-rectangle yes no 0 0.01 0 0.005

/solve/patch () hexahedron-r0 () temperature 1200

/solve/set/courant-number 0.5

/solve/iterate/ 50

/adapt/mark-inout-rectangle yes no 0 0.01 0 0.005

/solve/patch () hexahedron-r0 () temperature 1600

/solve/set/courant-number 1

/solve/iterate/ 50

/adapt/mark-inout-rectangle yes no 0 0.01 0 0.005

/solve/patch () hexahedron-r0 () temperature 1800

/solve/set/courant-number 3

/solve/iterate/ 50

/adapt/mark-inout-rectangle yes no 0 0.01 0 0.005

/solve/patch () hexahedron-r0 () temperature 2000

/solve/set/courant-number 5

/solve/iterate/ 50

/adapt/mark-inout-rectangle yes no 0 0.01 0 0.005

/solve/patch () hexahedron-r0 () temperature 2279

/solve/iterate/ 6000

/file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_coupled_w all.cas

/define/boundary-conditions/zone-type radius outflow

/solve/iterate/ 6000

/file/write-case-data/100/CH4_grav_2orde_100_lenonunity_mesh007035_coupled_a lloutflow.cas

q

q

exit

yes

Does someone experience the same problem?

Hi,

I have the same problem. Have you found a solution?

Greetings, Karel