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

ReactingFoam specie transport equation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 28, 2007, 19:15
Default During a discussion on IRC. a
  #1
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,880
Rep Power: 25
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 25C 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
__________________
Alberto

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
GeekoCFD 32bit - The 32bit edition of GeekoCFD.
GeekoCFD text mode - A smaller version of GeekoCFD, text-mode only, with only OpenFOAM. Available in a variety of virtual formats.
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: 7
smehdi609 is on a distinguished road
You are right alberto, I usually use fvm::laplacian(turbulence->muEff()/Sc, Yi) and set Sc=0.7.
smehdi609 is offline   Reply With Quote

Old   December 4, 2010, 12:03
Default
  #4
New Member
 
Join Date: Nov 2010
Posts: 10
Rep Power: 5
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 25C 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
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
ReactingFoam Laminar flame simulation transport properties gbansal OpenFOAM Running, Solving & CFD 3 November 23, 2012 16:03
Please explain the implementation of species transport Eqn in reactingFoam kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 07:29
Laminar diffusion flames and reactingFoam Transport properties jgaricano OpenFOAM Running, Solving & CFD 0 June 4, 2008 16:58
Poisson equation combined with transport equation Se-Hee CFX 0 December 27, 2007 01:00
transport equation CMB CD-adapco 2 July 9, 2004 05:19


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