From 5c4d714b4531c86019f5f778e78c25f59c5eabf1 Mon Sep 17 00:00:00 2001 From: sergio <s.ferraris@opencfd.co.uk> Date: Thu, 18 Mar 2021 11:16:27 -0700 Subject: [PATCH] ENH: Update icoReactingMultiPhaseInterFoam tutorials --- .../system/fvSolution | 8 ++++++ .../system/controlDict | 2 +- .../system/fvSchemes | 2 +- .../system/fvSolution | 28 +++++++++++-------- .../poolEvaporation/system/fvSolution | 5 ++++ .../solidMelting2D/system/fvSolution | 2 +- 6 files changed, 33 insertions(+), 14 deletions(-) diff --git a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/evaporationMultiComponent/system/fvSolution b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/evaporationMultiComponent/system/fvSolution index fa57dcc179c..b8d923c24db 100644 --- a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/evaporationMultiComponent/system/fvSolution +++ b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/evaporationMultiComponent/system/fvSolution @@ -39,6 +39,14 @@ solvers cYi 1; } + "pcorr.*" + { + solver PCG; + preconditioner DIC; + tolerance 1e-5; + relTol 0; + } + "(U.*|epsilon.*|k.*)" { solver smoothSolver; diff --git a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/controlDict b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/controlDict index 7cba53c1131..06be31bba10 100644 --- a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/controlDict +++ b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/controlDict @@ -29,7 +29,7 @@ deltaT 1e-5; writeControl adjustable; -writeInterval 0.02; +writeInterval 0.01; purgeWrite 0; diff --git a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSchemes b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSchemes index e07ef333cee..1a6c6483248 100644 --- a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSchemes +++ b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSchemes @@ -29,7 +29,7 @@ divSchemes { default none; - div(rhoPhi,U) Gauss linearUpwind gradU; + div(rhoPhi,U) Gauss limitedLinearV 1; "div\(phi,alpha.*\)" Gauss vanLeer; "div\(phir,alpha.*\)" Gauss linear; diff --git a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSolution b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSolution index 820815d6468..eaf78c1a7e3 100644 --- a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSolution +++ b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/inertMultiphaseMultiComponent/system/fvSolution @@ -39,7 +39,7 @@ solvers ); nAlphaCorr 2; - nAlphaSubCycles 2; + nAlphaSubCycles 3; // Compressiion factor for species in each alpha phase // NOTE: It should be similar to cAlpha @@ -74,19 +74,16 @@ solvers p_rgh { - solver PCG; - preconditioner DIC; - tolerance 1e-9; - relTol 0.01; - - nPreSweeps 0; - nPostSweeps 2; + solver GAMG; + smoother DICGaussSeidel; + tolerance 1e-11; + relTol 0.01; - cacheAgglomeration true; + cacheAgglomeration true; nCellsInCoarsestLevel 10; - agglomerator faceAreaPair; - mergeLevels 1; + agglomerator faceAreaPair; + mergeLevels 1; }; p_rghFinal @@ -94,10 +91,19 @@ solvers $p_rgh; relTol 0; }; + + "pcorr.*" + { + solver PCG; + preconditioner DIC; + tolerance 1e-5; + relTol 0; + } } PIMPLE { + correctPhi false; momentumPredictor false; nOuterCorrectors 1; nCorrectors 4; diff --git a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/poolEvaporation/system/fvSolution b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/poolEvaporation/system/fvSolution index 0fd5ca804a7..f142c63606e 100644 --- a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/poolEvaporation/system/fvSolution +++ b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/poolEvaporation/system/fvSolution @@ -49,6 +49,11 @@ solvers relTol 0; } + "pcorr.*" + { + $p_rgh; + } + "U.*" { solver smoothSolver; diff --git a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/solidMelting2D/system/fvSolution b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/solidMelting2D/system/fvSolution index c7a9c0ee262..863ae999ec7 100644 --- a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/solidMelting2D/system/fvSolution +++ b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/solidMelting2D/system/fvSolution @@ -36,7 +36,7 @@ solvers relTol 0; } - pcorr + "pcorr.*" { solver PCG; preconditioner DIC; -- GitLab