|
[Sponsors] |
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: 17 |
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<fvvectormatrix> 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: 17 |
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: 17 |
I have tried this:
tmp<fvvectormatrix> 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? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
what is the "stokes problems " | ztdep | Main CFD Forum | 0 | February 15, 2006 22:54 |
stokes number | guru | FLUENT | 0 | August 20, 2004 08:20 |
Stokes algorithm? | Rui | Main CFD Forum | 1 | March 17, 2003 13:14 |
help: I am trying to solve Navier Stokes compressible and viscid flow | Jose Choy | Main CFD Forum | 2 | May 18, 2000 05:45 |
any better way to solve 3D stokes flow? | Yangang Bao | Main CFD Forum | 0 | September 29, 1999 10:26 |