diff --git a/tutorials/multiphase/interIsoFoam/damBreak/system/fvSchemes b/tutorials/multiphase/interIsoFoam/damBreak/system/fvSchemes index 0974ab0dca..3471c6fb90 100644 --- a/tutorials/multiphase/interIsoFoam/damBreak/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/damBreak/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss limitedLinearV 1; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss interfaceCompression; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } diff --git a/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/constant/dynamicMeshDict b/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/constant/dynamicMeshDict index e9f00c12b9..ee75b90d65 100644 --- a/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/constant/dynamicMeshDict +++ b/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/constant/dynamicMeshDict @@ -47,7 +47,7 @@ correctFluxes (phi none) (nHatf none) (rhoPhi none) - (alphaPhi none) + (alphaPhi_ none) (ghf none) (phi0 none) (dVf_ none) diff --git a/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSchemes b/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSchemes index 880a89205c..cf772a5117 100644 --- a/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss upwind; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; diff --git a/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSolution b/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSolution index 0f414e2cea..25fc8f1703 100644 --- a/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/damBreakWithObstacle/system/fvSolution @@ -29,9 +29,8 @@ solvers writeIsoFaces false; reconstructionScheme plicRDF; // isoAlpha - nAlphaCorr 1; nAlphaSubCycles 1; - cAlpha 1; + cAlpha 1; // Read by interIsoFoam but not used } p_rgh diff --git a/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSchemes b/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSchemes index f3fb2d3ee6..2d9aaf9341 100644 --- a/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rho*phi,U) Gauss limitedLinearV 1; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss interfaceCompression; div((muEff*dev(T(grad(U))))) Gauss linear; } diff --git a/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSolution b/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSolution index 8f14cbb4d8..e9a5e58049 100644 --- a/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/discInConstantFlow/system/fvSolution @@ -19,12 +19,12 @@ solvers { "alpha.*" { - isoFaceTol 1e-8; - surfCellTol 1e-8; - nAlphaBounds 3; + isoFaceTol 1e-10; + surfCellTol 1e-10; + nAlphaBounds 5; snapTol 0; clip false; - cAlpha 1; + cAlpha 1; //Read by interIsoFoam but not used nAlphaSubCycles 1; reconstructionScheme plicRDF; // isoAlpha } diff --git a/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSchemes b/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSchemes index 27539b3d29..85c1ddb24f 100644 --- a/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss linearUpwind grad(U); - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } diff --git a/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSolution b/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSolution index 5bae5e24da..559f8d146d 100644 --- a/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/discInConstantFlowCyclicBCs/system/fvSolution @@ -27,17 +27,8 @@ solvers clip true; reconstructionScheme isoAlpha; // isoAlpha - nAlphaCorr 2; nAlphaSubCycles 1; - cAlpha 1; - - MULESCorr yes; - nLimiterIter 3; - - solver smoothSolver; - smoother symGaussSeidel; - tolerance 1e-8; - relTol 0; + cAlpha 1; // Read but not used by interIsoFoam } "pcorr.*" diff --git a/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSchemes b/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSchemes index d9dd7e1883..982056b31f 100644 --- a/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss limitedLinearV 1; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; } laplacianSchemes diff --git a/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSolution b/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSolution index 05fe08f1c6..629a2b6615 100644 --- a/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/discInReversedVortexFlow/system/fvSolution @@ -26,7 +26,7 @@ solvers clip false; nAlphaSubCycles 1; - cAlpha 1; + cAlpha 1; // Read by interIsoFoam but not used reconstructionScheme plicRDF; // isoAlpha } } diff --git a/tutorials/multiphase/interIsoFoam/iobasin/system/fvSchemes b/tutorials/multiphase/interIsoFoam/iobasin/system/fvSchemes index 675c85755a..6fcc0b01b5 100644 --- a/tutorials/multiphase/interIsoFoam/iobasin/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/iobasin/system/fvSchemes @@ -29,8 +29,6 @@ divSchemes { default none; div(rhoPhi,U) Gauss limitedLinearV 1; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; "div\(phi,(k|omega)\)" Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; diff --git a/tutorials/multiphase/interIsoFoam/iobasin/system/fvSolution b/tutorials/multiphase/interIsoFoam/iobasin/system/fvSolution index ab46cad2e4..c3246331ec 100644 --- a/tutorials/multiphase/interIsoFoam/iobasin/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/iobasin/system/fvSolution @@ -20,7 +20,7 @@ solvers "alpha.water.*" { nAlphaSubCycles 1; - cAlpha 1; + cAlpha 1; // Read by interIsoFoam but not used reconstructionScheme isoAlpha; // plicRDF isoFaceTol 1e-8; surfCellTol 1e-6; diff --git a/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSchemes b/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSchemes index e18a11069b..908966f393 100644 --- a/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss limitedLinearV 1; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } diff --git a/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSolution b/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSolution index f2db19bd4c..f9c1e7a451 100644 --- a/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/notchedDiscInSolidBodyRotation/system/fvSolution @@ -25,7 +25,7 @@ solvers nAlphaBounds 3; clip false; nAlphaSubCycles 1; - cAlpha 1; + cAlpha 1; // Read by interIsoFoam but not used reconstructionScheme isoAlpha; // plicRDF } } diff --git a/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSchemes b/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSchemes index 1ea5308e02..88493c5c6a 100644 --- a/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSchemes @@ -27,8 +27,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss vanLeerV; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss vanLeer; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } diff --git a/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSolution b/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSolution index dea662b317..85c405cb1f 100644 --- a/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/sloshingTank2D/system/fvSolution @@ -18,9 +18,8 @@ solvers { alpha.water { - nAlphaCorr 1; nAlphaSubCycles 3; - cAlpha 1.5; + cAlpha 1; // Read by interIsoFoam but not used isoFaceTol 1e-10; surfCellTol 1e-6; diff --git a/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSchemes b/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSchemes index 74eda64d19..67ebdf926f 100644 --- a/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss limitedLinearV 1; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; } laplacianSchemes diff --git a/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSolution b/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSolution index bb303b1d33..6216c70b2c 100644 --- a/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/sphereInReversedVortexFlow/system/fvSolution @@ -25,7 +25,7 @@ solvers nAlphaBounds 8; clip false; nAlphaSubCycles 1; - cAlpha 1; + cAlpha 1; // Read but not used by interIsoFoam reconstructionScheme isoAlpha; } } diff --git a/tutorials/multiphase/interIsoFoam/standingWave/system/fvSolution b/tutorials/multiphase/interIsoFoam/standingWave/system/fvSolution index f0db5bfd89..141b14462e 100644 --- a/tutorials/multiphase/interIsoFoam/standingWave/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/standingWave/system/fvSolution @@ -26,7 +26,7 @@ solvers clip true; nAlphaSubCycles 1; - cAlpha 1; + cAlpha 1; // Read by interIsoFoam but not used reconstructionScheme isoAlpha; // plicRDF } diff --git a/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSchemes b/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSchemes index b0d1975449..8d08d76b21 100644 --- a/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSchemes @@ -27,8 +27,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss linearUpwind grad(U); - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } diff --git a/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSolution b/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSolution index 7c55e4775d..f15e1510da 100644 --- a/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSolution +++ b/tutorials/multiphase/interIsoFoam/waveExampleStreamFunction/system/fvSolution @@ -24,9 +24,8 @@ solvers snapTol 1e-12; clip true; - nAlphaCorr 1; nAlphaSubCycles 1; - cAlpha 1; + cAlpha 1; // Read by interIsoFoam but not used reconstructionScheme plicRDF; // plicRDF } diff --git a/tutorials/multiphase/interIsoFoam/weirOverflow/system/fvSchemes b/tutorials/multiphase/interIsoFoam/weirOverflow/system/fvSchemes index 1d3baa2721..1d715d5a48 100644 --- a/tutorials/multiphase/interIsoFoam/weirOverflow/system/fvSchemes +++ b/tutorials/multiphase/interIsoFoam/weirOverflow/system/fvSchemes @@ -28,8 +28,6 @@ gradSchemes divSchemes { div(rhoPhi,U) Gauss linear; - div(phi,alpha) Gauss vanLeer; - div(phirb,alpha) Gauss linear; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div(phi,R) Gauss upwind;