From a67dfdb23a88226d5a4a8b9417b33d3d03096bba Mon Sep 17 00:00:00 2001 From: Henry Date: Mon, 3 Feb 2014 18:04:18 +0000 Subject: [PATCH] interFoam/ras/damBreak: updated for predictor-corrector MULES and changes to phi in pEqn --- .../interFoam/ras/damBreak/system/controlDict | 4 +- .../interFoam/ras/damBreak/system/fvSchemes | 7 +-- .../interFoam/ras/damBreak/system/fvSolution | 50 +++++++++++-------- .../ras/damBreakPorousBaffle/0/p_rgh | 8 +-- .../damBreakPorousBaffle/system/controlDict | 2 +- .../system/createBafflesDict | 4 +- .../ras/damBreakPorousBaffle/system/fvSchemes | 7 +-- .../damBreakPorousBaffle/system/fvSolution | 50 ++++++++++--------- 8 files changed, 68 insertions(+), 64 deletions(-) diff --git a/tutorials/multiphase/interFoam/ras/damBreak/system/controlDict b/tutorials/multiphase/interFoam/ras/damBreak/system/controlDict index 8d4bf8a10b..b507f78b95 100644 --- a/tutorials/multiphase/interFoam/ras/damBreak/system/controlDict +++ b/tutorials/multiphase/interFoam/ras/damBreak/system/controlDict @@ -47,8 +47,8 @@ runTimeModifiable yes; adjustTimeStep on; -maxCo 0.2; -maxAlphaCo 0.2; +maxCo 1; +maxAlphaCo 1; maxDeltaT 1; diff --git a/tutorials/multiphase/interFoam/ras/damBreak/system/fvSchemes b/tutorials/multiphase/interFoam/ras/damBreak/system/fvSchemes index 86fad259ec..3bc5347598 100644 --- a/tutorials/multiphase/interFoam/ras/damBreak/system/fvSchemes +++ b/tutorials/multiphase/interFoam/ras/damBreak/system/fvSchemes @@ -27,14 +27,11 @@ gradSchemes divSchemes { - div(rhoPhi,U) Gauss linear; + div(rhoPhi,U) Gauss linearUpwind grad(U); 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; - div(R) Gauss linear; - div(phi,nuTilda) Gauss upwind; div((muEff*dev(T(grad(U))))) Gauss linear; } @@ -58,7 +55,7 @@ fluxRequired default no; p_rgh; pcorr; - alpha; + alpha.water; } diff --git a/tutorials/multiphase/interFoam/ras/damBreak/system/fvSolution b/tutorials/multiphase/interFoam/ras/damBreak/system/fvSolution index f6725a947b..972bb3ec25 100644 --- a/tutorials/multiphase/interFoam/ras/damBreak/system/fvSolution +++ b/tutorials/multiphase/interFoam/ras/damBreak/system/fvSolution @@ -17,18 +17,27 @@ FoamFile solvers { - alpha.water + "alpha.water.*" { - nAlphaCorr 1; - nAlphaSubCycles 4; - cAlpha 2; + nAlphaCorr 2; + nAlphaSubCycles 1; + alphaOuterCorrectors yes; + cAlpha 1; + + MULESCorr yes; + nLimiterIter 3; + + solver PBiCG; + preconditioner DILU; + tolerance 1e-8; + relTol 0; } pcorr { solver PCG; preconditioner DIC; - tolerance 1e-10; + tolerance 1e-5; relTol 0; } @@ -42,35 +51,34 @@ solvers p_rghFinal { - solver PCG; - preconditioner DIC; - tolerance 1e-07; + $p_rgh; relTol 0; } - "(U|k|epsilon)" + "(U|k|epsilon).*" { - solver PBiCG; - preconditioner DILU; + solver smoothSolver; + smoother symGaussSeidel; tolerance 1e-06; relTol 0; - } - - "(U|k|epsilon)Final" - { - solver PBiCG; - preconditioner DILU; - tolerance 1e-08; - relTol 0; + minIter 1; } } PIMPLE { - momentumPredictor no; - nCorrectors 3; + momentumPredictor no; + nOuterCorrectors 1; + nCorrectors 3; nNonOrthogonalCorrectors 0; } +relaxationFactors +{ + equations + { + ".*" 1; + } +} // ************************************************************************* // diff --git a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/0/p_rgh b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/0/p_rgh index 0b3ae667f0..f24e6e1ede 100644 --- a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/0/p_rgh +++ b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/0/p_rgh @@ -58,18 +58,18 @@ boundaryField patchType cyclic; jump uniform 0; value uniform 0; - D 700; + D 1000; I 500; - length 1.05; + length 0.15; } porous_half1 { type porousBafflePressure; patchType cyclic; value uniform 0; - D 700; + D 1000; I 500; - length 1.05; + length 0.15; } } diff --git a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/controlDict b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/controlDict index bf6c452b46..53bb12ff93 100644 --- a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/controlDict +++ b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/controlDict @@ -29,7 +29,7 @@ deltaT 0.001; writeControl adjustableRunTime; -writeInterval 0.1; +writeInterval 0.05; purgeWrite 0; diff --git a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/createBafflesDict b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/createBafflesDict index 46e0c66bf8..789f61dbc4 100644 --- a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/createBafflesDict +++ b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/createBafflesDict @@ -49,9 +49,9 @@ baffles { type porousBafflePressure; patchType cyclic; - D 700; + D 1000; I 500; - length 1.05; + length 0.15; jump uniform 0; value uniform 0; } diff --git a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSchemes b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSchemes index 86fad259ec..3bc5347598 100644 --- a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSchemes +++ b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSchemes @@ -27,14 +27,11 @@ gradSchemes divSchemes { - div(rhoPhi,U) Gauss linear; + div(rhoPhi,U) Gauss linearUpwind grad(U); 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; - div(R) Gauss linear; - div(phi,nuTilda) Gauss upwind; div((muEff*dev(T(grad(U))))) Gauss linear; } @@ -58,7 +55,7 @@ fluxRequired default no; p_rgh; pcorr; - alpha; + alpha.water; } diff --git a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSolution b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSolution index a7103f6d62..972bb3ec25 100644 --- a/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSolution +++ b/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/fvSolution @@ -17,18 +17,27 @@ FoamFile solvers { - alpha.water + "alpha.water.*" { - nAlphaCorr 1; - nAlphaSubCycles 4; - cAlpha 2; + nAlphaCorr 2; + nAlphaSubCycles 1; + alphaOuterCorrectors yes; + cAlpha 1; + + MULESCorr yes; + nLimiterIter 3; + + solver PBiCG; + preconditioner DILU; + tolerance 1e-8; + relTol 0; } pcorr { solver PCG; preconditioner DIC; - tolerance 1e-10; + tolerance 1e-5; relTol 0; } @@ -42,41 +51,34 @@ solvers p_rghFinal { - solver PCG; - preconditioner DIC; - tolerance 1e-08; + $p_rgh; relTol 0; } - "(U|k|epsilon)" + "(U|k|epsilon).*" { - solver PBiCG; - preconditioner DILU; + solver smoothSolver; + smoother symGaussSeidel; tolerance 1e-06; relTol 0; - } - - "(U|k|epsilon)Final" - { - solver PBiCG; - preconditioner DILU; - tolerance 1e-07; - relTol 0; + minIter 1; } } PIMPLE { - momentumPredictor yes; - nCorrectors 3; + momentumPredictor no; + nOuterCorrectors 1; + nCorrectors 3; nNonOrthogonalCorrectors 0; } relaxationFactors { - "U.*" 1; - "k.*" 1; - "epsilon.*" 1; + equations + { + ".*" 1; + } } // ************************************************************************* //