Hi Hrvoje,
Sorry to not rep
Hi Hrvoje,
Sorry to not replying you for the last time (about a_N coefficients. I reread your thesys, but it stil not ok for me), but in emmergency, I want to know where I'm wrong in the implementation of a Stokes solver: I begin with the simpleFOam.C file: so I just removed the div(phi,U) and changed some "minus" by "plus", according to he discretisation of NavierStokes, explained in your thesys. (in my simulation, I put turbulence OFF) **************************** tmp<fvvectormatrix> UEqn (turbulence>divR(U)); UEqn().relax(); solve(UEqn() == fvc::grad(p)); p.boundaryField().updateCoeffs(); volScalarField AU = UEqn().A(); U = UEqn().H()/AU; UEqn.clear(); phi = fvc::interpolate(U) & mesh.Sf(); adjustPhi(phi, U, p); for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) { fvScalarMatrix pEqn ( // I put a "minus" here, because now: Up=[H(U)+grad(p)]/ap, so the discretisation of div(U)=0 gives me the equation of pressure below: fvm::laplacian(1.0/AU, p) == fvc::div(phi)); pEqn.setReference(pRefCell, pRefValue); pEqn.solve(); if (nonOrth == nNonOrthCorr) { // "+=" instead of "=" because F = S.Uf = S.[H(U)/ap]f + S.[grad(p)/ap]f phi += pEqn.flux(); } } # include "continuityErrs.H" p.relax(); // "+=" instead of "=" because Up=[H(U)+grad(p)]/ap U += fvc::grad(p)/AU; U.correctBoundaryConditions(); } ******************** the end of the program is the same. It seems so easy, but I can't get by myself. I just need your patience to tell me where I'm wrong. Thanks in advance, Lam 
just take a little look please
just take a little look please http://www.cfdonline.com/OpenFOAM_D...part/happy.gif

? You have messed up the signs
? You have messed up the signs: just delete the convection terms and use a symmetric solver for the momentum equation.
I haven't got time to go through this in detail, but it should really be as simple as that. Hrv 
Hi Hrvoje,
Thanks to reply
Hi Hrvoje,
Thanks to reply me. Ok, I'm going to delete the convection term... even if I don't understand why the signs are wrong. Hope that you will explain me when you get more time. Cheers, Lam 
All times are GMT 4. The time now is 03:27. 