- **OpenFOAM Verification & Validation**
(*http://www.cfd-online.com/Forums/openfoam-verification-validation/*)

- - **Numerics in OpenFoam - PBiCG vs. PCG**
(*http://www.cfd-online.com/Forums/openfoam-verification-validation/108464-numerics-openfoam-pbicg-vs-pcg.html*)

Numerics in OpenFoam - PBiCG vs. PCGDear all,
to give an accurate validation for some channelFlow-case, my Prof wants me to estimate how good the numeric methods are, that i'm using in OpenFoam. But honestly not having much Background in numerical-Math but a bit in physics, it's hard for me to grasp even the concept of the piso solver. Using Gauss linear-Approximation everywhere for spacial discretisation shemes and using Chrank Nicholson shemes for time discretisation, I think this Methods are used to solve the euqations of the piso-loop - what i don't know where does PCG and PBiCG come into play - as i found there are in use for every Field U p... and if i use Chrank Nicholoson 0.5 which is implicit - what kind of sheme is used therefore - more or less every field should be involved in this linear system? Is it that for every gradient divergence and so on... a PCG or PBicG method will be used? Sorry if this question are not precise, but this might be cause i don't understand everything completely, if you could give recommand some liturature it would be nice. Kind regards Le Frog |

Hi Frog,
Meindert already mentioned, the difference between a PCG or a PBiCG,... you can find in any textbook about linear iterative solvers like Wessling, Saad, Trefethen, Stran,.... I think the book from Saad you find free on the www. But I think your problems are a bit different. Sine OF is a segregated solver, for every field (variable) there has to be solved a linear system (i.e. for U, p,...) you have AU=b, Bp=c, A;B matrices.Your choice of the numerical schemes (entries for the operator discretization in fvSchemes, i.e. (linear, upwind, limitedLinear,...) for the operators is respobsible of the coefficients of this matrices. So the choice of operators regarding f.i. U are responsible for the properties/coefficients of the matrix AThe choice of your solution/solving algorithm (i.e. entries in fvSolution) is responsible which iterative solver is used to solve your linear system. So choosing f.i. PCG for U and PBiCG for p, AU=b is solved iteratively bc PCG, Bp=c with PBiCGand so forth. Being also a mathematician, as a book enlighting CFD and numerics, I would also recommend Ferziger/Peric: Computational methods for fluid dynamics Th. |

Thank you very much meindert and thg,
this was helpfull. |

All times are GMT -4. The time now is 12:15. |