CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Verification & Validation

Numerics in OpenFoam - PBiCG vs. PCG

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

Like Tree7Likes
  • 1 Post By meindert
  • 6 Post By thg

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 24, 2012, 10:47
Default Numerics in OpenFoam - PBiCG vs. PCG
  #1
New Member
 
Marwan Sabih
Join Date: Jun 2011
Posts: 8
Rep Power: 14
LeFrog is on a distinguished road
Dear 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
LeFrog is offline   Reply With Quote

Old   October 24, 2012, 15:27
Default
  #2
Member
 
Meindert de Groot
Join Date: Jun 2012
Location: Netherlands
Posts: 34
Rep Power: 13
meindert is on a distinguished road
The PCG method and PBicG method will give a different solution to the system of equations, say

\mathbf{A}\mathbf{x} = \mathbf{b}.

Given a certain tolerance \varepsilon, each method will give a solution \mathbf{x} that satisfies

||\mathbf{A}\mathbf{x}-\mathbf{b}|| < \varepsilon

for some norm ||\cdot||. So \mathbf{A}\mathbf{x} \approx \mathbf{b}.

If you would like to read a bit about these iterative methods, I would recommend reading 'Principles of computational fluid dynamics' by P. Wesseling. I should warn you, the book is written from a mathematical point of view. It might be quite challenging to read.

However, my guess is that your professor wants you to look at different numerical schemes, each giving a different system of equations.
Sourabhsup likes this.
meindert is offline   Reply With Quote

Old   October 25, 2012, 07:19
Default
  #3
thg
Member
 
Thorsten Grahs
Join Date: Oct 2009
Posts: 61
Rep Power: 16
thg is on a distinguished road
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 A

The 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 PBiCG
and 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.
thg is offline   Reply With Quote

Old   October 25, 2012, 11:00
Default
  #4
New Member
 
Marwan Sabih
Join Date: Jun 2011
Posts: 8
Rep Power: 14
LeFrog is on a distinguished road
Thank you very much meindert and thg,

this was helpfull.
LeFrog is offline   Reply With Quote

Reply

Tags
numerics pcg pbicg

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Floating point exception error Alan OpenFOAM Running, Solving & CFD 11 July 1, 2021 22:51
[Gmsh] 2D Mesh Generation Tutorial for GMSH aeroslacker OpenFOAM Meshing & Mesh Conversion 12 January 19, 2012 04:52
OpenFOAM 1.6.x, 1.7.0 and 1.7.x are not fully prepared to work with gcc-4.5.x wyldckat OpenFOAM Bugs 18 October 21, 2010 06:51
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 3 September 8, 2010 07:25
The OpenFOAM extensions project mbeaudoin OpenFOAM 16 October 9, 2007 10:33


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