CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Another way to solve Stokes

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 25, 2007, 09:57
Default Hi all, I want to use an ot
  #1
lam
Member
 
Hoang-Lam
Join Date: Mar 2009
Location: Lausanne, Switzerland
Posts: 60
Rep Power: 8
lam is on a distinguished road
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
lam is offline   Reply With Quote

Old   June 25, 2007, 12:15
Default nobody has an idea?
  #2
lam
Member
 
Hoang-Lam
Join Date: Mar 2009
Location: Lausanne, Switzerland
Posts: 60
Rep Power: 8
lam is on a distinguished road
nobody has an idea?
lam is offline   Reply With Quote

Old   June 26, 2007, 05:54
Default I have tried this: tmp
  #3
lam
Member
 
Hoang-Lam
Join Date: Mar 2009
Location: Lausanne, Switzerland
Posts: 60
Rep Power: 8
lam is on a distinguished road
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?
lam is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
what is the "stokes problems " ztdep Main CFD Forum 0 February 15, 2006 23:54
stokes number guru FLUENT 0 August 20, 2004 08:20
Stokes algorithm? Rui Main CFD Forum 1 March 17, 2003 14: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


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