# Another way to solve Stokes

 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 25, 2007, 09:57 Hi all, I want to use an ot #1 Member   Hoang-Lam Join Date: Mar 2009 Location: Lausanne, Switzerland Posts: 60 Rep Power: 8 Hi all, I want to use an other way to solve a steady Stokes equation for incompressible flow: div(U)=0 laplacian(nu,U)=Grad(p) (2) by div(U)=epsi*p, where epsi is very small i put this expression in (2), and I have: Laplacian(nu,U)=grad(div(U)/epsi) (2') But my problem is I have not managed to discretize "grad(div(U)/epsi)" I've tried : tmp UEqn (fvm::laplacian(nuEff(), U)); solve(UEqn()==fvc::grad(div(U))); But it tells me someting like in compilation: " error: no matching function for call to 'div(Foam::volVectorField&) " I think it a intressing way to test, even if it could be expensive for the memory. Hope taht anybody will give me some advice to make a good discretisation for (2') Thanks in advance, Lam

 June 25, 2007, 12:15 nobody has an idea? #2 Member   Hoang-Lam Join Date: Mar 2009 Location: Lausanne, Switzerland Posts: 60 Rep Power: 8 nobody has an idea?

 June 26, 2007, 05:54 I have tried this: tmp #3 Member   Hoang-Lam Join Date: Mar 2009 Location: Lausanne, Switzerland Posts: 60 Rep Power: 8 I have tried this: tmp UEqn (turbulence->divR(U)); solve(UEqn()==-fvc::grad((fvc::div(U))/epsi)); where epsi is a dimensioned scalar [0 -2 1 0 0 0 0] But my solution doesn't obey the continuity! Can somebody tell me where I am wrong?

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post ztdep Main CFD Forum 0 February 15, 2006 23:54 guru FLUENT 0 August 20, 2004 08:20 Rui Main CFD Forum 1 March 17, 2003 14:14 Jose Choy Main CFD Forum 2 May 18, 2000 05:45 Yangang Bao Main CFD Forum 0 September 29, 1999 10:26

All times are GMT -4. The time now is 10:10.