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

IcoFoam Is it a NavierStokesEquation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 10, 2008, 06:24
Default Hi, I try to analyse the ic
  #1
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
Hi,

I try to analyse the icoFoam Solver.
The Solver solves the Equation:
ddt(U)+div(phi,U)-Laplacian(nu,U)=-grad(p)

this themes seems to be a Navier-Stokes-Equation but I miss an U at div(phi,U).
The Navier-Stokes Equation I know sould be:
ddt(U)+U*div(phi,U)-Laplacian(nu,U)=-grad(p)

Does the U appear later in the discretisation or is it a variation of the Navier-Stokes-Equation?

Thx,
Jörn
joern is offline   Reply With Quote

Old   September 10, 2008, 08:33
Default No: the N-S equation that _you
  #2
Senior Member
 
Gavin Tabor
Join Date: Mar 2009
Posts: 181
Rep Power: 8
grtabor is on a distinguished road
No: the N-S equation that _you_ know is

ddt(U)+U*div(U)-Laplacian(nu,U)=-grad(p)

Since div(u)=0 we can rearrange this as

ddt(U)+div(U*U)-Laplacian(nu,U)=-grad(p)

The finite volume discretisation however linearises this by representing one of the U terms in div(U*U) as the flux phi, hence

ddt(U)+div(phi,U)-Laplacian(nu,U)=-grad(p)

We then need to iterate around a bit to take account of the fact that when U changes from this solution, phi must be updated...

Clear now?

Gavin
grtabor is offline   Reply With Quote

Old   September 10, 2008, 08:52
Default Thanks for the fast help. Thi
  #3
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
Thanks for the fast help.
This helps a lot.

I didn't see the step:
ddt(U)+U*div(U)-Laplacian(nu,U)=-grad(p)
to
ddt(U)+div(U*U)-Laplacian(nu,U)=-grad(p)

Joern
joern is offline   Reply With Quote

Old   September 10, 2008, 09:24
Default Here it is, with subscript x,y
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,592
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Here it is, with subscript x,y,z meaning differentiation and u,v,w being velocity components and U the vector:

In the x-direction div(U * U) reads:

(u * u)_x + (u * v)_y + (u * w)_z =
u_x u + u u_x + u_x u + u v_y + u_z w + u w_z =
u_x u + u_x u + u_z w + u*(u_x + v_y + w_z) =
u_x u + u_x u + u_z w

where the term in the bracket are identical to zero for incompressible flows, as it is your continuity equation.

The same can be done for the other directions and you will have understood the step above.

Best regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   September 10, 2008, 10:29
Default Thx for your Help. I think
  #5
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
Thx for your Help.

I think I now understand the solver algorithm. Now I have to understand the discretisation ;)

Joern
joern is offline   Reply With Quote

Old   November 26, 2008, 18:26
Default Hallo, i have a new problem
  #6
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
Hallo,

i have a new problem with the same topic

I try to understand the sover sonicFoam.
This solver solves the Equations:
ddt(rho, U)+div(phi, U)-laplacian(mu, U)=-grad(p)
and
ddt(rho, e)+div(phi, e)-laplacian(mu, e)=-div(phi/rho)+mu*magSqr(symm(grad(U)))

For me it looks like some kind of impulse and energy equation from the compressible Navier Stokes equations.

Does someone know which equations sonicFoam solves or where i can find some literature about the algorithm an his mathematic roots?


sorry for posting nearly the same question again, but i need the infos for my further work and i don't see it alone.

thx,
Joern
joern is offline   Reply With Quote

Old   November 27, 2008, 18:43
Default sonicfoam is a solver for comp
  #7
New Member
 
Mark Michael
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 5
Rep Power: 8
mark_michael is on a distinguished road
sonicfoam is a solver for compressible ideal gas flow. That why the solver use the continuity-,momentum-, energy-, and ideal gas equation in its algorithm.
look in the programmer's guide ! there is an nice tutorial named: "Supersonic flow over a forward-facing step" on page 58 !
http://foam.sourceforge.net/doc/Guid...mmersGuide.pdf
mark_michael is offline   Reply With Quote

Old   June 15, 2009, 13:44
Default
  #8
Member
 
Joern Bader
Join Date: Mar 2009
Posts: 33
Rep Power: 8
joern is on a distinguished road
I reactivate this old thread, cause i have a question about the sonicFoam solver.

sonicFoam is a solver with PISO loop for the compressible Navier-Stokes equations. I think i understand most of the solver but one think is not clear:

if i set mu=0 in sonicFoam it solves the Euler Equations for a compressible, inviscid, perfect gas. (thats what i want to simulate)
But there is a little difference. sonicFoam just looks at the inner energy e=T*Cv (just thermal energy). The Euler Equations (and also the Navier-Stokes) work with the energy e=T*Cv+1/2*mag(v)^2 (thermal energy + kinetic energy).

why can the kinetic energy be ignored? And why is the the solver still consisten to the Euler-/Navier-Stokes Equations?


thx for help,
Joern
joern 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
IcoFoam aap OpenFOAM Running, Solving & CFD 15 May 28, 2012 08:30
Density in icoFoam Densidad en icoFoam manuel OpenFOAM Running, Solving & CFD 8 September 22, 2010 04:10
About phi in icoFoam kar OpenFOAM Running, Solving & CFD 3 February 20, 2008 06:20
Possible bug in icoFoam msrinath80 OpenFOAM Bugs 6 November 19, 2007 18:35
IcoFoam on AIX 53 ds2taieb OpenFOAM Installation 1 March 24, 2006 04:22


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