|October 27, 2008, 11:04||
Dear all, I'm new to OF and
Join Date: Mar 2009
Location: Stockholm, Stockholm, Sweden
Posts: 73Rep Power: 9
I'm new to OF and trying to implement a Poisson solver to calculate a pressure field from a velocity field. For the simplest case, inviscid and incompressible that would be, laplace(p) = -div(div(UU)).
After reading U and setting pP (pressure to be derived) initially to p, using the cavity case in IcoFoam as a test case, the Poisson equation is implemented as
fvm::laplacian(pP) == -fvc::div(fvc::div(phi, U))
and solved with the GAMG solver. The problem is that I'm getting strange results. Is there anything wrong with the implementation? What is the difference with the following implementation, (which does not want to work)
solve(pPEqn == fvm::laplacian(pP));
Is there a direct way of setting the boundary conditions for pP to fixed values, something similar to boundaryField().types()?
|April 17, 2009, 16:59||
did you success to implement your Poisson solver?
I'm managing to implement exactly the same solver, to get a pressure field from an incompressible velocity field. Do you obtain good results?
|Thread||Thread Starter||Forum||Replies||Last Post|
|looking for a 3D Poisson Solver||phsieh2005||Main CFD Forum||2||March 1, 2007 20:33|
|Different ans of poisson eqn with different solver||Quarkz||Main CFD Forum||0||October 6, 2005 11:49|
|Poisson solver||Steve||Main CFD Forum||6||July 22, 2004 21:16|
|solver for Poisson eqn||billboard||Main CFD Forum||3||May 22, 2002 00:07|
|Solver for poisson equ.||Capi||Main CFD Forum||1||May 3, 2001 11:20|