From 8cd5301c634e98b5b60bb893e5b463256bdffb9b Mon Sep 17 00:00:00 2001 From: danielque Date: Wed, 17 Jan 2024 11:42:07 +0100 Subject: [PATCH] add pMinFactor to fvSolution of compressible cases for OF>4 we are using pressure control instead of density control in cfdemSolverRhoPimple and cfdemSolverRhoPimpleChem which require pMin/pMax or pMinFactor/pMaxFactor instead of rhoMin/rhoMax; also corrected rhoMin in single particle chemistry tutorials for expected rho --- .../HeatTransferRanzMarshall/CFD/system/fvSolution | 8 ++++---- .../PolydisperseFluidizedBed/R1_FB/CFD/system/fvSolution | 6 ++++-- .../PolydisperseFluidizedBed/R2_FB/CFD/system/fvSolution | 6 ++++-- .../PolydisperseFluidizedBed/R3_FB/CFD/system/fvSolution | 6 ++++-- .../30CO6H2_Hematite/CFD/system/fvSolution | 6 ++++-- .../33CO3H2_Hematite/CFD/system/fvSolution | 8 +++++--- .../GOD0k1_Hematite/CFD/system/fvSolution | 8 +++++--- .../SingleParticleCases/ISO4695/CFD/system/fvSolution | 8 +++++--- 8 files changed, 35 insertions(+), 21 deletions(-) diff --git a/tutorials/cfdemSolverRhoPimple/HeatTransferRanzMarshall/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimple/HeatTransferRanzMarshall/CFD/system/fvSolution index 3bc33d5f..e9dc23a8 100644 --- a/tutorials/cfdemSolverRhoPimple/HeatTransferRanzMarshall/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimple/HeatTransferRanzMarshall/CFD/system/fvSolution @@ -69,10 +69,10 @@ PIMPLE nCorrectors 1; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.5; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; - pMinFactor 0.5; - pMaxFactor 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.5; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors diff --git a/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R1_FB/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R1_FB/CFD/system/fvSolution index 121600f9..47107880 100755 --- a/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R1_FB/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R1_FB/CFD/system/fvSolution @@ -95,8 +95,10 @@ PIMPLE nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors diff --git a/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R2_FB/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R2_FB/CFD/system/fvSolution index a6754a58..70b7dbbe 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R2_FB/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R2_FB/CFD/system/fvSolution @@ -94,8 +94,10 @@ PIMPLE nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors diff --git a/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R3_FB/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R3_FB/CFD/system/fvSolution index 121600f9..47107880 100755 --- a/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R3_FB/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimpleChem/PolydisperseFluidizedBed/R3_FB/CFD/system/fvSolution @@ -95,8 +95,10 @@ PIMPLE nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/CFD/system/fvSolution index 47b89b11..c9298e70 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/CFD/system/fvSolution @@ -109,8 +109,10 @@ PIMPLE nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.1; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/CFD/system/fvSolution index 2ebb617d..d5c88c23 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/CFD/system/fvSolution @@ -22,7 +22,7 @@ solvers preconditioner DIC; tolerance 1e-6; relTol 0.01; - /*solver GAMG; + /*solver GAMG; tolerance 1e-06; relTol 0.01; smoother GaussSeidel; @@ -109,8 +109,10 @@ PIMPLE nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.1; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/CFD/system/fvSolution index 2b3f9225..92a33cf5 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/CFD/system/fvSolution @@ -41,7 +41,7 @@ solvers tolerance 1e-06; relTol 0; } - + "(rho|G)" { solver PCG; @@ -109,8 +109,10 @@ PIMPLE nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.1; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/CFD/system/fvSolution index 2b3f9225..92a33cf5 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/CFD/system/fvSolution +++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/CFD/system/fvSolution @@ -41,7 +41,7 @@ solvers tolerance 1e-06; relTol 0; } - + "(rho|G)" { solver PCG; @@ -109,8 +109,10 @@ PIMPLE nOuterCorrectors 5; nCorrectors 2; nNonOrthogonalCorrectors 0; - rhoMin rhoMin [ 1 -3 0 0 0 ] 0.4; - rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; + rhoMin rhoMin [ 1 -3 0 0 0 ] 0.1; // OF4 + rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0; // OF4 + pMinFactor 0.5; // OF5, OF6 + pMaxFactor 2.0; // OF5, OF6 } relaxationFactors