diff --git a/applications/solvers/compressible/rhoPimpleFoam/UEqn.H b/applications/solvers/compressible/rhoPimpleFoam/UEqn.H index 8dfdcc5d9..3925650ec 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/UEqn.H +++ b/applications/solvers/compressible/rhoPimpleFoam/UEqn.H @@ -7,9 +7,16 @@ tmp UEqn + turbulence->divDevRhoReff(U) ); -if (oCorr == nOuterCorr-1) +if (oCorr == nOuterCorr - 1) { - UEqn().relax(1); + if (mesh.solutionDict().relax("UFinal")) + { + UEqn().relax(mesh.solutionDict().relaxationFactor("UFinal")); + } + else + { + UEqn().relax(1); + } } else { diff --git a/applications/solvers/compressible/rhoPorousMRFPimpleFoam/UEqn.H b/applications/solvers/compressible/rhoPorousMRFPimpleFoam/UEqn.H index 50291428a..ea80c2e7d 100644 --- a/applications/solvers/compressible/rhoPorousMRFPimpleFoam/UEqn.H +++ b/applications/solvers/compressible/rhoPorousMRFPimpleFoam/UEqn.H @@ -9,7 +9,14 @@ tmp UEqn if (oCorr == nOuterCorr-1) { - UEqn().relax(1); + if (mesh.solutionDict().relax("UFinal")) + { + UEqn().relax(mesh.solutionDict().relaxationFactor("UFinal")); + } + else + { + UEqn().relax(1); + } } else { diff --git a/applications/solvers/incompressible/pimpleDyMFoam/UEqn.H b/applications/solvers/incompressible/pimpleDyMFoam/UEqn.H index ae7c8edfe..aacdd5137 100644 --- a/applications/solvers/incompressible/pimpleDyMFoam/UEqn.H +++ b/applications/solvers/incompressible/pimpleDyMFoam/UEqn.H @@ -5,13 +5,20 @@ + turbulence->divDevReff(U) ); - if (ocorr != nOuterCorr - 1) + if (ocorr == nOuterCorr - 1) { - UEqn.relax(); + if (mesh.solutionDict().relax("UFinal")) + { + UEqn.relax(mesh.solutionDict().relaxationFactor("UFinal")); + } + else + { + UEqn.relax(1); + } } else { - UEqn.relax(mesh.solutionDict().relaxationFactor("UFinal")); + UEqn.relax(); } if (momentumPredictor) diff --git a/applications/solvers/incompressible/pimpleFoam/UEqn.H b/applications/solvers/incompressible/pimpleFoam/UEqn.H index 577745d7e..64ffbea02 100644 --- a/applications/solvers/incompressible/pimpleFoam/UEqn.H +++ b/applications/solvers/incompressible/pimpleFoam/UEqn.H @@ -7,9 +7,16 @@ tmp UEqn + turbulence->divDevReff(U) ); -if (oCorr == nOuterCorr-1) +if (oCorr == nOuterCorr - 1) { - UEqn().relax(1); + if (mesh.solutionDict().relax("UFinal")) + { + UEqn().relax(mesh.solutionDict().relaxationFactor("UFinal")); + } + else + { + UEqn().relax(1); + } } else { diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/UEqn.H b/applications/solvers/multiphase/twoLiquidMixingFoam/UEqn.H index 1bd6e2526..b6a7bb6ed 100644 --- a/applications/solvers/multiphase/twoLiquidMixingFoam/UEqn.H +++ b/applications/solvers/multiphase/twoLiquidMixingFoam/UEqn.H @@ -16,7 +16,14 @@ if (oCorr == nOuterCorr-1) { - UEqn.relax(1); + if (mesh.solutionDict().relax("UFinal")) + { + UEqn.relax(mesh.solutionDict().relaxationFactor("UFinal")); + } + else + { + UEqn.relax(1); + } } else {