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

 immortality January 9, 2013 13: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 14: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 06:17

Dear Ehsan,

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

Regards

 immortality December 12, 2013 09: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 09:35

1 Attachment(s)
Dear Ehsan,

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

Many thanks

 wyldckat December 28, 2013 14:25

 immortality December 28, 2013 18: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 07: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 15:24.