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

ReactingFoam specie transport equation

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

Like Tree3Likes
  • 2 Post By alberto
  • 1 Post By smehdi609

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 28, 2007, 20:15
Default During a discussion on IRC. a
  #1
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
During a discussion on IRC. a doubt about reactingFoam implementation came out.

In reactingFoam, the transport equation for the specie mass fration contains the diffusional term in the form:

fvm::laplacian(turbulence->muEff(), Yi)

which means that the effective specie diffusivity is assumed to be equal to the effective viscosity:

Deff = D_lam + D_turb = nu_lam/Sc_lam + nu_turb/Sc_turb = nu_lam + nu_turb

which means that

Sc_lam = Sc_turb = 1

Why this assumption is done?

For example, in a laminar calculation, if we consider water, the viscosity of the fluid and the diffusivity of a specie, for example the oxygen, can be very different. At 25°C we have:

D_02 ~ 2.0 * 10^-10 m^2/s
nu_H20 ~ 1.0 * 10^-4 m^2/s

In turbulent flows, the turbulent Schmidt number is set to 0.7 in other codes, but however it can change too. For example, in a pipe with passive transport it's between 0.6 and 0.9.

Thanks in advance for any clarification.

With kind regards,
Alberto
TommyM and labyrinth01 like this.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   March 27, 2009, 08:44
Default
  #2
New Member
 
Lu,Y.
Join Date: Mar 2009
Posts: 2
Rep Power: 0
ly.foam is on a distinguished road
I guess the reason is that when dealing with turbulent combustion, the molecular transport are often small compared to their tubulent counterparts and so are neglected.
But I think it is very important of you to have pointed it out that when dealing with laminar combustion with OpenFoam, special attention may have to be paid to the issue. But I am not very sure about this.
ly.foam is offline   Reply With Quote

Old   March 27, 2009, 13:22
Default
  #3
Member
 
M. Mahdi Salehi
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 50
Rep Power: 17
smehdi609 is on a distinguished road
You are right alberto, I usually use fvm::laplacian(turbulence->muEff()/Sc, Yi) and set Sc=0.7.
labyrinth01 likes this.
smehdi609 is offline   Reply With Quote

Old   December 4, 2010, 13:03
Default
  #4
New Member
 
Join Date: Nov 2010
Posts: 10
Rep Power: 15
siavash_abadeh is on a distinguished road
Quote:
Originally Posted by alberto View Post
During a discussion on IRC. a doubt about reactingFoam implementation came out.

In reactingFoam, the transport equation for the specie mass fration contains the diffusional term in the form:

fvm::laplacian(turbulence->muEff(), Yi)

which means that the effective specie diffusivity is assumed to be equal to the effective viscosity:

Deff = D_lam + D_turb = nu_lam/Sc_lam + nu_turb/Sc_turb = nu_lam + nu_turb

which means that

Sc_lam = Sc_turb = 1

Why this assumption is done?

For example, in a laminar calculation, if we consider water, the viscosity of the fluid and the diffusivity of a specie, for example the oxygen, can be very different. At 25°C we have:

D_02 ~ 2.0 * 10^-10 m^2/s
nu_H20 ~ 1.0 * 10^-4 m^2/s

In turbulent flows, the turbulent Schmidt number is set to 0.7 in other codes, but however it can change too. For example, in a pipe with passive transport it's between 0.6 and 0.9.

Thanks in advance for any clarification.

With kind regards,
Alberto
dear alberto
i have a problem and i think you can help me.
my problem is mass transfer from solid to liquid with electrochemical method and i want simulate this problem with openfoam and i cant find solver for it.
how resolvent your problem?
use a solver or write a code for it.
if you write a code : can i have your code?
very very thank in advance
siavash_abadeh 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
ReactingFoam Laminar flame simulation transport properties gbansal OpenFOAM Running, Solving & CFD 3 November 23, 2012 17:03
Please explain the implementation of species transport Eqn in reactingFoam kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 08:29
Laminar diffusion flames and reactingFoam Transport properties jgaricano OpenFOAM Running, Solving & CFD 0 June 4, 2008 17:58
Poisson equation combined with transport equation Se-Hee CFX 0 December 27, 2007 02:00
transport equation CMB Siemens 2 July 9, 2004 06:19


All times are GMT -4. The time now is 14:25.