CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   an error on PBiCG (http://www.cfd-online.com/Forums/openfoam-solving/111550-error-pbicg.html)

immortality January 9, 2013 12:59

an error on PBiCG
 
its strange.why rhoPimpleFoam shows an error on its default solver as I copied?U and others are not symmetric but why openFoam wants a symmetric matrix solver?please submit your idea.thanks.
Code:

Starting time loop

Mean and max Courant Numbers = 0.1422220012 0.5709693515
deltaT = 9.807881957e-09
Time = 9.80788e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0


--> FOAM FATAL IO ERROR:
Unknown symmetric matrix solver PBiCG

Valid symmetric matrix solvers are :

4
(
GAMG
ICCG
PCG
smoothSolver
)


file: /home/thesis/Desktop/WR_Pimple/system/fvSolution::solvers::(U|h|R|k|epsilon|omega) from line 43 to line 46.

    From function lduMatrix::solver::New
    in file matrices/lduMatrix/lduMatrix/lduMatrixSolver.C at line 77.

Code:

solvers
{
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance      1e-06;
        relTol          0;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    "rho.*"
    {
        $p;
        tolerance      1e-05;
        relTol          0;
    }

    "(U|h|R|k|epsilon|omega)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance      1e-05;
        relTol          0;
    }

    "(U|h|R|k|epsilon|omega)Final"
    {
        $U;
        relTol          0;
    }
}

PIMPLE
{
    momentumPredictor yes;
    nOuterCorrectors 1;
    nCorrectors    2;
    nNonOrthogonalCorrectors 0;
    rhoMin          rhoMin [ 1 -3 0 0 0 ] 0.5;
    rhoMax          rhoMax [ 1 -3 0 0 0 ] 2.0;
}


immortality January 9, 2013 13:16

very strange behaviour of openFoam.when I changed it the error became reciprocal:
Code:

Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1        0.85034;
    alphaK2        1;
    alphaOmega1    0.5;
    alphaOmega2    0.85616;
    Prt            1;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1          0.075;
    beta2          0.0828;
    betaStar        0.09;
    a1              0.31;
    c1              10;
}

Reading thermophysicalProperties


Starting time loop

Mean and max Courant Numbers = 0.1422220012 0.5709693515
deltaT = 9.807881957e-09
Time = 9.80788e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for Ux, Initial residual = 0.001234410983, Final residual = 1.598743971e-09, No Iterations 1
DICPCG:  Solving for Uy, Initial residual = 0.9999999996, Final residual = 4.396979854e-07, No Iterations 1
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
swak4Foam: Allocating new repository for sampledGlobalVariables
DICPCG:  Solving for e, Initial residual = 0.001547344897, Final residual = 2.015913876e-09, No Iterations 1


--> FOAM FATAL IO ERROR:
Unknown asymmetric matrix solver PCG

Valid asymmetric matrix solvers are :

4
(
BICCG
GAMG
PBiCG
smoothSolver
)


file: /home/thesis/Desktop/WR_Pimple/system/fvSolution::solvers::(U|e|h|R|k|epsilon|omega) from line 43 to line 46.

    From function lduMatrix::solver::New
    in file matrices/lduMatrix/lduMatrix/lduMatrixSolver.C at line 106.

FOAM exiting

thesis@thesis-X58A-UD7:~/Desktop/WR_Pimple$ rhoCentralFoam | tee log
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.1.0-0bc225064152
Exec  : rhoCentralFoam
Date  : Feb 09 2012
Time  : 21:41:49
Host  : "thesis-X58A-UD7"
PID    : 6985
Case  : /home/thesis/Desktop/WR_Pimple
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field U

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1        0.85034;
    alphaK2        1;
    alphaOmega1    0.5;
    alphaOmega2    0.85616;
    Prt            1;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1          0.075;
    beta2          0.0828;
    betaStar        0.09;
    a1              0.31;
    c1              10;
}

Reading thermophysicalProperties


Starting time loop

Mean and max Courant Numbers = 0.1422220012 0.5709693515
deltaT = 9.807881957e-09
Time = 9.80788e-09

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0


--> FOAM FATAL IO ERROR:
Unknown symmetric matrix solver PBiCG

Valid symmetric matrix solvers are :

4
(
GAMG
ICCG
PCG
smoothSolver
)


file: /home/thesis/Desktop/WR_Pimple/system/fvSolution::solvers::(U|e|h|R|k|epsilon|omega) from line 43 to line 46.

    From function lduMatrix::solver::New
    in file matrices/lduMatrix/lduMatrix/lduMatrixSolver.C at line 77.

FOAM exiting


Maimouna December 12, 2013 05:17

Dear Ehsan,

I got the same error that you posted above. Do you know how to solve it?

Regards

immortality December 12, 2013 08:21

Hi Maimouna,
how are you?
I remember the issue but not sure what was the origin,I think I had copied fvSolution of another solver return of rhoPimpleFoam and it didn't match to rhoPimpleFoam,
if you can send your fvSolution dictionary maybe I remember what was the cause exactly.

Maimouna December 12, 2013 08:35

1 Attachment(s)
Dear Ehsan,

many thanks for your reply.

I'm good. I attached the case that I'm working.

Many thanks

wyldckat December 28, 2013 13:25

FYI: Maimouna's case is being addressed on this thread: http://www.cfd-online.com/Forums/ope...hron-case.html

immortality December 28, 2013 17:43

Hi dear Maimouna,
sorry I haven't had time for looking into the case,have you found the source of the issue?

Maimouna December 30, 2013 06:39

Dear Ehsan,

thank you for your reply. I still searching for the source of the issue.


All times are GMT -4. The time now is 14:10.