/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "rho.*"
    {
        solver          diagonal;
    }
    
    p_rgh
    {
        solver           GAMG;
        tolerance        1e-4; // changed from 1e-7
        relTol           0.1; // changed from0.01

        smoother         DIC;

    }
    
    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }
    
    UFinal // this section is newly added
   {
    solver          PBiCGStab; // newly added
    preconditioner  DILU;     
    tolerance       1e-04;  // Maximum allowed residual for velocity components . newly added  changed from 1e-05
    relTol          0.1;    // Relative tolerance for velocity field convergence . newly added
   }

    hFinal // this section is newly added
   {
    solver          PBiCGStab;
    preconditioner  DILU;
    tolerance       1e-04;  // Maximum allowed residual for enthalpy (or temperature) changed from 1e-05
    relTol          0.1;    // Relative tolerance for enthalpy (or temperature) convergence
   }
   
   epsilonFinal  // this section is newly added
   {
    solver          PBiCGStab;
    preconditioner  DILU;
    tolerance      1e-04;  // Maximum allowed residual for epsilon (dissipation rate) changed from 1e-05
    relTol         0.1;    // Relative tolerance for epsilon convergence
   }
   
   kFinal  // this section is newly added
   {
    solver          PBiCGStab;
    preconditioner  DILU;
    tolerance      1e-04;  // Maximum allowed residual for epsilon (dissipation rate)
    relTol         0.1;    // Relative tolerance for epsilon convergence
   }

    "(U|h|k|epsilon|omega)"
    {
        solver          PBiCGStab; // if converge is'nt happening try GAMG
        preconditioner  DILU;
        tolerance       1e-4; // changed from 1e-8
        relTol          0.1;  
        smoother        DIC; // newly added
    }
}

PIMPLE
{
    momentumPredictor  yes; // no;
    nNonOrthogonalCorrectors 1; // originally was 0
    nCorrectors       2;
    nInnerCorrectors 10;   // Number of inner correctors for each outer corrector . newly added
    tolerance 1e-04;       // Tolerance for convergence . newly added , changed on 10/04 from 1e-06
    pRefCell          0;
    pRefValue         0;
    correctPhi        yes;
}

relaxationFactors
{
    fields
    {
        p               0.5;     // newly added
        rho             0.7;      // originally was 1
        p_rgh           0.7;
    }
    equations
    {
        U               0.3;
        h               0.3;
        "(k|epsilon|omega)" 0.5;     // originally was 0.7
    }
}


// ************************************************************************* //
