
[Sponsors] 
October 11, 2008, 17:55 
some question on matrix in FVM

#1 
Guest
Posts: n/a

I am new to CFD and ask some simple questions. when NS equations is solved over unstructured grid based on SIMPLE coupling method, matrices for the discretized momentum and pressure correction equations are symmetrical or nonsymmetrical? When we may have symmetrical or nonsymmetrical matrices?


October 14, 2008, 23:58 
Re: some question on matrix in FVM

#2 
Guest
Posts: n/a

I would appreciate your help.


October 15, 2008, 02:14 
Re: some question on matrix in FVM

#3 
Guest
Posts: n/a

pressure can be symm.
but mom. can be unsym. due to usung upwinding (if you use central schemes, not feasible for convective terms, you mom. can be symm. too) 

October 15, 2008, 03:00 
Re: some question on matrix in FVM

#4 
Guest
Posts: n/a

Than you very much. Here is my understanding of your replies, symm has nothing to with mesh topology, right? I mean for hybrid grids containing mixing elements such as hex, tetra, and prism, and for grids containing hex only, pressure equations are both symm. This sounds a little strange. For hybrid grids, different volumes may have different unknown variable, e.g. 5 for a tetra and y for a hexa. it seems symm is not about how many unknown variables a discretised equation for a control volume has. How should I to understand symm?


October 15, 2008, 03:18 
Re: some question on matrix in FVM

#5 
Guest
Posts: n/a

> symm has nothing to with mesh topology, right?
yes, first (you use simple method )i assume you mean incompressible flow; pressure eq. is resulted fro discritization of Poisson equation like this Del^P = rhs in fact a steady diffusion equation for which the best schemes are central schems, so each quantity has same share from each cell and so system of eq. will be symm. but we should be careful for boundary cells, but it is possible to apply bc so system will be symmetric > How should I to understand symm? try to derive your system on a small uniform cartesin grid, then you will see how momentum eq. will be obviously unsymmetric and how it would be possible to have symm. pressure eq. assume neighbour cells A and B, and a shared face between them, assume flow direction is from cell A to B, so using first order upwind method (consider just convectie terms), momentum equation of cell B has a contribution from cell A, but cell A does not have any contribution from cell B which leads to unsymmetric system of equation 

October 15, 2008, 04:40 
Re: some question on matrix in FVM

#6 
Guest
Posts: n/a

thx


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Unanswered question  niklas  OpenFOAM  2  July 31, 2013 16:03 
PISO Question  titio  OpenFOAM Running, Solving & CFD  1  July 6, 2009 05:33 
Accessing the motionU matrix in OF ver 14  philippose  OpenFOAM Running, Solving & CFD  2  April 16, 2007 18:16 
Elemtary matrix to CSR global matrix  xueying  Main CFD Forum  2  September 24, 2002 09:44 
question  K.L.Huang  CDadapco  1  March 29, 2000 04:57 