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

Caracteristic length used for Peclet number

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

Like Tree1Likes
  • 1 Post By alexeym

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 2, 2015, 11:00
Default Caracteristic length used for Peclet number
  #1
Senior Member
 
Laurent DASTUGUE
Join Date: May 2014
Location: Chartres, France
Posts: 122
Rep Power: 11
laurentD is on a distinguished road
Hi foamers,
Does anybody know what is the caracteristic length used when OpenFOAM calculate the Peclet number ?
I don't have access these days to OpenFOAM and i have to answer this question.
Thanks a lot.
Laurent
laurentD is offline   Reply With Quote

Old   June 2, 2015, 11:38
Default
  #2
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

It seems you have access to the Internet.

https://github.com/OpenFOAM/OpenFOAM...d/Pe/Pe.C#L118

https://github.com/OpenFOAM/OpenFOAM...olation.C#L136

So, Pe is calculated in assumption Sct = 1 and length scale is cell size.
alexeym is offline   Reply With Quote

Old   June 2, 2015, 11:51
Default Sct ?
  #3
Senior Member
 
Laurent DASTUGUE
Join Date: May 2014
Location: Chartres, France
Posts: 122
Rep Power: 11
laurentD is on a distinguished road
Hi Alexey,
thank you for the links, and for your answer concerning the length scale.
Just one question : What is Sct ?
Laurent
laurentD is offline   Reply With Quote

Old   June 2, 2015, 12:00
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Sct is turbulent Schmidt number (https://en.wikipedia.org/wiki/Schmidt_number). As you can see, in Pe.C everything is divided by viscosity, yet for Peclet one has to divide by diffusivity. If we assume that Sc = 1, we can just divide by viscosity (for turbulent case we have to assume that Sc = 1 and Sct = 1, so nuEff( = nu + nut) == Deff( = nu/Sc + nut/Sct)).
alexeym is offline   Reply With Quote

Old   June 3, 2015, 05:56
Default
  #5
Senior Member
 
Laurent DASTUGUE
Join Date: May 2014
Location: Chartres, France
Posts: 122
Rep Power: 11
laurentD is on a distinguished road
Ok i understand now.
Is the Peclet number calculated here a massic one ? (Pe=Re*Sc) or a thermic one (Pe=Re*Pr).
If this is a massic one, can i say that Pe = Re, since Sc = 1 ?
Laurent
laurentD is offline   Reply With Quote

Old   June 3, 2015, 06:08
Default
  #6
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

If you assume Pr = Prt (turbulent Prandtl number) = 1, then it becomes "termic" And yes, you can call the value calculated by Pe utility Reynolds number.
alexeym is offline   Reply With Quote

Old   June 3, 2015, 06:18
Default
  #7
Senior Member
 
Laurent DASTUGUE
Join Date: May 2014
Location: Chartres, France
Posts: 122
Rep Power: 11
laurentD is on a distinguished road
In fact what i want to do is the calculation of the convection coefficient h.
I know that there is a relation between h and the Nusselt number, and i know a relation letting me to have Nusselt, by the use of Re and Pr values.
Since i have the value of Prandtl number easily, knowing the fluid caracteristics, i just need the value of Re to have Nu and to have finally the value of h.
So reading your posts, i understand the following :
using Pe utility, i can do the operation Re = Pe/Pr to have Reynolds number.
Am i right ? Or must i apply some operations on Pe utility results before having the Peclet number value which will let me calculate Re ?
Laurent
laurentD is offline   Reply With Quote

Old   June 3, 2015, 08:54
Default
  #8
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

Short answer: it depends.

Pe field produced by Pe utility is surfaceScalarField and its value is

\text{Pe}_{\text{calculated by Pe utility}} = \frac{Lu}{\nu}

where L is of order of cell size, u is in fact

u = \frac{\vec{u}\times\vec{S}}{|\vec{S}|}

where S is face normal vector with module equal to the surface of the face.

So, in fact Pe is Re.

If you are happy with h as a surfaceScalarField, then yes, you just manipulate result of Pe utility; if, for example, you need convection coefficient as volume field, you need additional steps.
wayne14 likes this.
alexeym is offline   Reply With Quote

Old   June 3, 2015, 09:07
Default
  #9
Senior Member
 
Laurent DASTUGUE
Join Date: May 2014
Location: Chartres, France
Posts: 122
Rep Power: 11
laurentD is on a distinguished road
Hi,
thank you very much for your help.
So what we have is : Pe = Pe / Pr, isn't it ?
Why is this utility not called Re ???
Laurent
laurentD is offline   Reply With Quote

Old   June 3, 2015, 10:37
Default
  #10
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

Quote:
Originally Posted by laurentD View Post
So what we have is : Pe = Pe / Pr, isn't it ?
Yes, though not sure if I get the notation right.

Quote:
Why is this utility not called Re ???
I don't know. Due to certain historical reasons?
alexeym is offline   Reply With Quote

Old   June 3, 2015, 11:14
Default
  #11
Senior Member
 
Laurent DASTUGUE
Join Date: May 2014
Location: Chartres, France
Posts: 122
Rep Power: 11
laurentD is on a distinguished road
Thank you very much Alexey for your time. Have a good day.
Maybe this utility will be called Re in a next future thanks to your explanations ;-)
Laurent
laurentD 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
[snappyHexMesh] snappyHexMesh sticking point natty_king OpenFOAM Meshing & Mesh Conversion 11 February 20, 2024 10:12
decomposePar no field transfert Jeanp OpenFOAM Pre-Processing 3 June 18, 2022 13:01
AMI speed performance danny123 OpenFOAM 21 October 24, 2020 05:13
simpleFoam parallel AndrewMortimer OpenFOAM Running, Solving & CFD 12 August 7, 2015 19:45
Stable boundaries marcoymarc CFX 33 March 13, 2013 07:39


All times are GMT -4. The time now is 12:24.