diff --git a/tutorials/multiphase/icoReactingMultiPhaseInterFoam/evaporationMultiComponent/system/fvSolution b/tutorials/multiphase/icoReactingMultiPhaseInterFoam/evaporationMultiComponent/system/fvSolution index fa57dcc179ccddf5b6b38dd1596c9cf2656006ff..b8d923c24dba7a09d05b22775371195a36543153 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 7cba53c1131b561311f80dbef6e1ee43d6a43141..06be31bba1068cb9e801e7ad238d53ef0fab500b 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 e07ef333cee50b7ea225fb20c31cbec87f6b6713..1a6c64832489229a49d12af469e382ade2fc2bad 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 820815d6468c27f024ca69bcc2ee71e4b62588fb..eaf78c1a7e3dfa54093cf8dd4ca6121057b436c5 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 0fd5ca804a7d7b6fd2112609f7cd0e64b18ffdab..f142c63606e6a4b5baa8f29c066a0bb30977e0e7 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 c7a9c0ee26244af6cc8fba4e48661528b43efb3c..863ae999ec7d4be42a297fed8d4c93957f1febd5 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;