|October 26, 2010, 05:09||
Strcture of pressure correction equation in PISO stable over time ?
Join Date: Feb 2010
Posts: 7Rep Power: 7
I have a problem in understanding why the things I'm doing work. I'm using the PISO algorithm to solve the Navier-Stokes equations. In there I have to solve the Poisson-like pressure equation which I'm solving with an algebraic multigrid solver.
Unfortunatly my PISO algorithm does not use a real Poisson pressure equation like -/LAPLACE/ p= RHS, but it includes the flux a caused by this correction into the operator.
This is, I have - /DIV/ a /DIV/ p = RHS.
For the standard PISO with a real Poisson pressure equation it is pretty clear that the structure of the matrix resulting from discretising this will never change over the time steps. But how is this in my case ?
I see that the matrix' structure doesn't really change so I can use the AMG coarsening once and once again without giving a new setup. But this is too heuristic in my eyes.
So my question is: Can I argue from any physical point of view that the structure of my matrix stays as it is ?
Does anyone have an idea or a hint where I could read about this ?
|October 26, 2010, 22:19||
Join Date: Mar 2009
Posts: 188Rep Power: 8
PISO algorithm was a method proposed for pressure-velocity coupling of incompressible flows. There are many ways in which the algorithm has been used in the literature. Check the paper based on which the code was written. It might be a good starting point to understand.
|Thread||Thread Starter||Forum||Replies||Last Post|
|calculation diverge after continue to run||zhajingjing||OpenFOAM||0||April 28, 2010 04:35|
|mass flow in is not equal to mass flow out||saii||CFX||2||September 18, 2009 08:07|
|Convergence moving mesh||lr103476||OpenFOAM Running, Solving & CFD||30||November 19, 2007 15:09|
|Transient simulation not converging||skabilan||OpenFOAM Running, Solving & CFD||12||September 17, 2007 17:48|
|Pressure Correction Equation||morteza||OpenFOAM Running, Solving & CFD||2||September 4, 2007 06:16|