From bc2d925d61f85d4d373e526ec93ccf9d614f009d Mon Sep 17 00:00:00 2001 From: Henry <Henry> Date: Thu, 7 Nov 2013 18:38:58 +0000 Subject: [PATCH] chemFoam tutorials: set ODE solver to seulex which performs best overall --- tutorials/combustion/chemFoam/README | 6 ++---- .../combustion/chemFoam/gri/constant/chemistryProperties | 2 +- .../combustion/chemFoam/h2/constant/chemistryProperties | 2 +- .../combustion/chemFoam/ic8h18/constant/chemistryProperties | 2 +- .../combustion/chemFoam/nc7h16/constant/chemistryProperties | 2 +- 5 files changed, 6 insertions(+), 8 deletions(-) diff --git a/tutorials/combustion/chemFoam/README b/tutorials/combustion/chemFoam/README index 50e94283bce..2cd8c892732 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 314f83cc2da..bf097e0ba39 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 314f83cc2da..b00cac494d4 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 314f83cc2da..bf097e0ba39 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 aee69a3239e..a2a58619975 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; } -- GitLab