CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

porousSimpleFoam: A fundamental inconsistency to formula in Fluent

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 18, 2018, 08:13
Default porousSimpleFoam: A fundamental inconsistency to formula in Fluent
  #1
Member
 
Kai
Join Date: May 2010
Location: Shanghai
Posts: 61
Blog Entries: 1
Rep Power: 15
kaifu is on a distinguished road
Hi, Foamers

The formulation in solver porousSimpleFoam is based on the articles: Porous Media in OpenFOAM - Chalmers
http://www.tfd.chalmers.se/~hani/kur...ukurReport.pdf

in the article, the invicid Navier-Stokes Eqn is written down as
\frac{\partial}{\partial t}(\gamma \rho u_i)+u_j \frac{\partial}{\partial x_j}( \rho u_i)=-\frac{\partial p}{\partial x_i} +S_i \:\:\:[1]

I think the velocity here refers to the superficial velocity.

However, according to the formulation in Fluent and other ref:
[1]http://www.afs.enea.it/project/neptu...ug/node233.htm
[2]Novel porous media formulation for multiphase flow conservation equations, by Sha and Chao 2007
the Navier-Stokes Eqn should be
\frac{\partial}{\partial t}(\gamma \rho u_i)+\frac{\partial}{\partial x_j}( \gamma \rho u_i u_j )=-\gamma \frac{\partial p}{\partial x_i}+S_i
and here the velocity is physical velocity. The formula are derived according to time/volume averaging.

And if we substitute the above Eqn with u_{\rm superficial } = \gamma u_{\rm physical }, we have
\frac{\partial}{\partial t}(\gamma \rho u_i)+\frac{\partial}{\partial x_j}(  \rho u_i u_j )=-\gamma^2 \frac{\partial p}{\partial x_i}+ \gamma S_i \:\:\:[2]
how to explain the inconsistency of the coefficient in pressure term between Eqns [1] and [2]?
__________________
Kai

Last edited by kaifu; May 18, 2018 at 08:15. Reason: something wrong with solver's name
kaifu is offline   Reply With Quote

Old   May 18, 2018, 11:32
Default
  #2
Member
 
Kai
Join Date: May 2010
Location: Shanghai
Posts: 61
Blog Entries: 1
Rep Power: 15
kaifu is on a distinguished road
One of the possibilities may come from the temporal/spatial inertial term. Both of them may be negligible, compared to the pressure and viscous resistance in the pore. But still sounds something is wrong..
__________________
Kai
kaifu is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFX vs. FLUENT turbo CFX 4 April 13, 2021 09:08
heat transfer with RANS wall function, over a flat plate (validation with fluent) bruce OpenFOAM Running, Solving & CFD 6 January 20, 2017 07:22
batch mode fluent - beginner CompIng FLUENT 1 November 30, 2015 12:48
Fluent 12.0 is worst then Fluent 6.2 herntan FLUENT 5 December 14, 2009 03:57
Fluent formula for Nusselt number Alan FLUENT 1 May 29, 2001 13:50


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