|
[Sponsors] |
Poisson equation: non-symmetric matrix with non-uniform spacing? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 25, 2012, 04:22 |
Poisson equation: non-symmetric matrix with non-uniform spacing?
|
#1 |
New Member
Daniel Weber
Join Date: Nov 2009
Posts: 4
Rep Power: 17 |
Hello,
I'm trying solve a poisson equation with finite differences. According to "Computational Methods for Fluid Dynamics" the finite difference approximation of the laplace operator in the 1D poisson equation is p_i+1 (x-i - x_i-1) + p_i-1(x_i+1 - x_i) - p_i (x_i+1 - x_i) / [ 1/2 * (x_i+1 - x_i-1) * (x_i+1 - x_i ) * ( x_i - x_i-1) ] If the grid spacing is uniform the matrix is symmetric, but not for non-uniform grid spacing. With non-uniformity the requirements is that the distance between each pair of odd and even unknowns (x_i - x_i+2, x_i+1 - x_i+3) must be equal. I could artificially move the unknows a bit, but generally that limits the transformation function of the grid. I would like to use a CG method, so is there anything I can do modify the scheme. (or maybe use finite volumes?) How is the pressure correction done in Finite Volume methods? Are the matrices as well non-symmetric? Any hint or a pointer to literature would be very helpful! Thank you! Daniel |
|
July 25, 2012, 04:54 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73 |
Quote:
Yes, on non-uniform grids the matrix is symmetric in its shape but not in the entry values |
||
July 25, 2012, 09:21 |
|
#3 | |
Member
G. S.
Join Date: Nov 2010
Posts: 54
Rep Power: 16 |
Quote:
Best regards... |
||
July 26, 2012, 04:03 |
|
#4 |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,291
Rep Power: 35 |
Poisson operator should give you symmetric matrix ie
Aij = Aji It is true in case of non uniform meshes too. |
|
July 26, 2012, 04:10 |
|
#5 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73 |
||
July 26, 2012, 06:52 |
|
#6 |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,291
Rep Power: 35 |
||
July 26, 2012, 12:16 |
|
#7 |
New Member
Daniel Weber
Join Date: Nov 2009
Posts: 4
Rep Power: 17 |
Hi,
thanks for the answers. Yes, I meant non-symmetric in the values. But I'm still curios if there are as well non-symmetric matrices(in the values) when I use a finite volume discretization for the pressure correction. Does anybody have i hint? best regards, Daniel |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Free jet simulation | msarkar | OpenFOAM Running, Solving & CFD | 39 | May 11, 2021 14:21 |
[Commercial meshers] fluentMeshToFoam multidomain mesh conversion problem | Attesz | OpenFOAM Meshing & Mesh Conversion | 12 | May 2, 2013 11:52 |
buoyantSimpleFoam divergence | TomS | OpenFOAM Running, Solving & CFD | 0 | November 24, 2011 15:16 |
Neumann Boundary Condition for Poisson Equation solution in Polar Coordinates | prapanj | Main CFD Forum | 2 | July 30, 2011 20:07 |
Problem with Poisson equation with Neumann BC | Jeff67 | Main CFD Forum | 0 | March 27, 2011 12:23 |