diff --git a/tutorials/combustion/chemFoam/README b/tutorials/combustion/chemFoam/README index 50e94283bcefb7958faaaf29660e35b9882b669b..2cd8c89273293ba04ca4114a55f86b1ac3c62feb 100644 --- a/tutorials/combustion/chemFoam/README +++ b/tutorials/combustion/chemFoam/README @@ -1,13 +1,11 @@ Stiff chemistry solver validation test cases gri : GRI-Mech 3.0. CH4 combustion, 53 species, 325 reactions - h2 : H2 combustion, 10 species, 27 reactions - nc7h16 : n-Heptane combustion, 544 species, 2446 reactions - ic8h18 : iso-Octane combustion, 874 species, 3796 reactions Results interpreted in 'validation' sub-folder, where OpenFOAM results -are compared against those predicted by Senkin (CHEMKIN II) +are compared against those predicted by CHEMKIN II. +Overall the best performing ODE solver is seulex followed closely by rodas32. diff --git a/tutorials/combustion/chemFoam/gri/constant/chemistryProperties b/tutorials/combustion/chemFoam/gri/constant/chemistryProperties index 314f83cc2da0f7328af993b2d37d22543bcf6ac4..bf097e0ba39679d2ec60ded2d59bb8ffe51fa227 100644 --- a/tutorials/combustion/chemFoam/gri/constant/chemistryProperties +++ b/tutorials/combustion/chemFoam/gri/constant/chemistryProperties @@ -33,7 +33,7 @@ EulerImplicitCoeffs odeCoeffs { - solver rodas32; + solver seulex; absTol 1e-12; relTol 1e-1; } diff --git a/tutorials/combustion/chemFoam/h2/constant/chemistryProperties b/tutorials/combustion/chemFoam/h2/constant/chemistryProperties index 314f83cc2da0f7328af993b2d37d22543bcf6ac4..b00cac494d4cf528afe6947d96595c7a8058ce49 100644 --- a/tutorials/combustion/chemFoam/h2/constant/chemistryProperties +++ b/tutorials/combustion/chemFoam/h2/constant/chemistryProperties @@ -33,7 +33,7 @@ EulerImplicitCoeffs odeCoeffs { - solver rodas32; + solver seulex; //rodas32; absTol 1e-12; relTol 1e-1; } diff --git a/tutorials/combustion/chemFoam/ic8h18/constant/chemistryProperties b/tutorials/combustion/chemFoam/ic8h18/constant/chemistryProperties index 314f83cc2da0f7328af993b2d37d22543bcf6ac4..bf097e0ba39679d2ec60ded2d59bb8ffe51fa227 100644 --- a/tutorials/combustion/chemFoam/ic8h18/constant/chemistryProperties +++ b/tutorials/combustion/chemFoam/ic8h18/constant/chemistryProperties @@ -33,7 +33,7 @@ EulerImplicitCoeffs odeCoeffs { - solver rodas32; + solver seulex; absTol 1e-12; relTol 1e-1; } diff --git a/tutorials/combustion/chemFoam/nc7h16/constant/chemistryProperties b/tutorials/combustion/chemFoam/nc7h16/constant/chemistryProperties index aee69a3239ea0ba38d36b2315cbb9331a78c73a5..a2a586199750a45e70f0b577a05c075a57b213e0 100644 --- a/tutorials/combustion/chemFoam/nc7h16/constant/chemistryProperties +++ b/tutorials/combustion/chemFoam/nc7h16/constant/chemistryProperties @@ -33,7 +33,7 @@ EulerImplicitCoeffs odeCoeffs { - solver rodas32; + solver seulex; absTol 1e-14; relTol 1e-1; }