CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

an error on PBiCG

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   January 9, 2013, 13:59
Default an error on PBiCG
  #1
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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 is offline   Reply With Quote

Old   January 9, 2013, 14:16
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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

Last edited by wyldckat; December 12, 2013 at 18:01. Reason: Fixed broken [CODE][/CODE]
immortality is offline   Reply With Quote

Old   December 12, 2013, 06:17
Default
  #3
Senior Member
 
ok
Join Date: Oct 2013
Posts: 257
Rep Power: 4
Maimouna is on a distinguished road
Dear Ehsan,

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

Regards
Maimouna is offline   Reply With Quote

Old   December 12, 2013, 09:21
Default
  #4
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   December 12, 2013, 09:35
Default
  #5
Senior Member
 
ok
Join Date: Oct 2013
Posts: 257
Rep Power: 4
Maimouna is on a distinguished road
Dear Ehsan,

many thanks for your reply.

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

Many thanks
Attached Files
File Type: gz TurekHron2.tar.gz (11.0 KB, 4 views)
Maimouna is offline   Reply With Quote

Old   December 28, 2013, 14:25
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,314
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
FYI: Maimouna's case is being addressed on this thread: Fluid structure interaction - Turek and Hron case
wyldckat is online now   Reply With Quote

Old   December 28, 2013, 18:43
Default
  #7
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Hi dear Maimouna,
sorry I haven't had time for looking into the case,have you found the source of the issue?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   December 30, 2013, 07:39
Default
  #8
Senior Member
 
ok
Join Date: Oct 2013
Posts: 257
Rep Power: 4
Maimouna is on a distinguished road
Dear Ehsan,

thank you for your reply. I still searching for the source of the issue.
Maimouna is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
smoothSolver diverges - solution in using PBiCG solver? makaveli_lcf OpenFOAM Running, Solving & CFD 3 September 11, 2013 12:44
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 9 January 3, 2012 09:06
Velocity blows up suddenly after 30,000+ iterations lordvon OpenFOAM 14 October 27, 2010 11:13
Problems with simulating TurbFOAM barath.ezhilan OpenFOAM 13 July 16, 2009 05:55


All times are GMT -4. The time now is 17:28.