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

Mistype in Darcy-Forchheimer porosity model description

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 10, 2009, 04:27
Default Mistype in Darcy-Forchheimer porosity model description
  #1
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Hi all!

I noticed a mistype in description of Darcy-Forchheimer model at
http://foam.sourceforge.net/doc/Doxy...ce.html#l00036

it is

S = - (mu*d*U + rho*|U|*f/2) * U

but should be

S = - (mu*d + rho*|U|*f/2) * U

(according to the literature and source code implementation).
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 5, 2011, 13:09
Default
  #2
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
acordingo to the header porousZone.H it has been fixed, if there was this problem in first place:
*************************************8
Implemented porosity models:

powerLaw (\e C0 and \e C1 parameters)
\f[
S = - \rho C_0 |U|^{(C_1 - 1)/2} U
\f]

Darcy-Forchheimer (@e d and \e f parameters)
\f[
S = - (\mu \, d + \frac{\rho |U|}{2} \, f) U
\f]
**************************************************
calim_cfd is offline   Reply With Quote

Old   August 5, 2011, 16:26
Default
  #3
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Yeah, I posted this 2 years before)))
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 6, 2011, 12:16
Default
  #4
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
hi maka!,glad to see you're still around! i'm new to OP and i'm trying to use the porosity solver and i bumped into ur post :P


since u've dealt with this topic b4 i'm taking this chance and i'll ask you.. do you know how to apply a permeability using in a cylindrical CS? i got a filter in a shape of a tube and i need to set its permeability in the radial direction.. any clues??

i found a few topics in the forums but the ones that seem to fit my need seemed to indicate some major changes in the codes and that's a lengthy path... do you know any shortcuts??

thanks a lot!!
calim_cfd is offline   Reply With Quote

Old   August 8, 2011, 02:25
Default
  #5
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Hi!

Don't know what do you really need, but you can use e.g. a source term in your momentum equation and adjust it accordin to the local radial direction (just calculate some vector field field at the beginning). But that is apart of porous class. You can also hack it to meet your needs.
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   August 8, 2011, 07:56
Default
  #6
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
Hello again Maka!

I'll try to be clearer. I want to model a porous zone. To achieve so, i need to set resistences in space, at all directions. Thing is, how do i "pick" a radial direction? i need a cylindrical coordinate system so that i can, lets say, set a vector similar to (100 1e7 1e7) representing a small resistence in direction e1=er and infinite in the others, i. e., e2(etheta)=e3(ez). I guess, It all boils down to: can i, and if so, how do i build a Cylindrical CS in OP.... so that i can set sth like:

coordinateSystem cylindrical_123;
// porosity 0.781;
Darcy
{
d d [0 -2 0 0 0] (x y z);
f f [0 -1 0 0 0] ( a b c);
}

http://imageshack.us/photo/my-images...radeteladn.png


thanks a lot :P


***********************

i found this post with the link and i guess i may have found what i needed.. it has a whole lof of other things..it's worth checking it out

Porous Case

Last edited by calim_cfd; August 8, 2011 at 10:23. Reason: photo
calim_cfd is offline   Reply With Quote

Old   December 11, 2014, 17:07
Default
  #7
Member
 
Hao Chen
Join Date: Aug 2014
Posts: 62
Rep Power: 3
hchen is on a distinguished road
Hi Mauricio:

I know it was a post long time ago...
But I am wondering if you manage to introduce a cylindrical coordinate into porous media? Since I am doing that right now!

Thanks!

Best regards
Hao
Quote:
Originally Posted by calim_cfd View Post
Hello again Maka!

I'll try to be clearer. I want to model a porous zone. To achieve so, i need to set resistences in space, at all directions. Thing is, how do i "pick" a radial direction? i need a cylindrical coordinate system so that i can, lets say, set a vector similar to (100 1e7 1e7) representing a small resistence in direction e1=er and infinite in the others, i. e., e2(etheta)=e3(ez). I guess, It all boils down to: can i, and if so, how do i build a Cylindrical CS in OP.... so that i can set sth like:

coordinateSystem cylindrical_123;
// porosity 0.781;
Darcy
{
d d [0 -2 0 0 0] (x y z);
f f [0 -1 0 0 0] ( a b c);
}

http://imageshack.us/photo/my-images...radeteladn.png


thanks a lot :P


***********************

i found this post with the link and i guess i may have found what i needed.. it has a whole lof of other things..it's worth checking it out

Porous Case
hchen 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
RasInterFoam and LRR turbulence model dmoroian OpenFOAM Bugs 7 March 28, 2011 17:28
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
LES and combustion model Margherita Cadorin CFX 0 October 29, 2008 06:24
Kato-Launder model sam Main CFD Forum 13 September 21, 2006 10:15
Two-Phase Buoyant Flow Issue Miguel Baritto CFX 4 August 31, 2006 12:02


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