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

reactingFoam non unity Lewis number

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

Like Tree2Likes
  • 1 Post By clapointe
  • 1 Post By clapointe

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 22, 2021, 12:45
Default reactingFoam non unity Lewis number
  #1
New Member
 
Lorenzo
Join Date: Apr 2021
Posts: 9
Rep Power: 3
lollo15 is on a distinguished road
Hello everyone !

I'm new with reactingFoam ( OpenFoam v8) and i'm looking for the option or where set a Lewis number different from 1.

I know that exist the file nonUnityLewisEddyDiffusivity.H that define Sc and Prt so i'm trying to use it instead of the conventional method.

Someone knows how to do it ??

Best regards !
Thank you!

Lorenzo
lollo15 is offline   Reply With Quote

Old   April 22, 2021, 13:33
Default
  #2
Senior Member
 
Join Date: Aug 2015
Posts: 494
Rep Power: 13
clapointe is on a distinguished road
Checking out the relevant .H file (nonUnityLewis) you'll see the following :

Code:
    LES
    {
        model         nonUnityLewisEddyDiffusivity;
        Prt             0.85;
        Sct             0.7;
    }
This would go into a thermophysicalTransport file (similarly to how you have a momentumTransport file). The description mentions for RAS or LES, so either option should be available. There are a few tutorials that include a thermophysicalTransport file (e.g., https://github.com/OpenFOAM/OpenFOAM...sicalTransport) that could also be helpful.

Caelan
lollo15 likes this.
__________________
Public git repository : https://github.com/clapointe2011/public
clapointe is offline   Reply With Quote

Old   April 23, 2021, 04:52
Default
  #3
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 517
Blog Entries: 1
Rep Power: 14
dlahaye is on a distinguished road
Wonderful, thx! Can you point to references in which the non-unit Lewis model has been validated?
dlahaye is offline   Reply With Quote

Old   April 23, 2021, 17:20
Default
  #4
Senior Member
 
Join Date: Aug 2015
Posts: 494
Rep Power: 13
clapointe is on a distinguished road
I've unfortunately not used that functionality, so I'm not familiar if/where it's been used since its release. I imagine, though, that updating the reactingFoam tutorial for the Sandia or DLR flames would be a good place to start, given that the experimental data is readily available.

Caelan
dlahaye likes this.
__________________
Public git repository : https://github.com/clapointe2011/public
clapointe is offline   Reply With Quote

Old   April 25, 2021, 17:13
Default
  #5
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 517
Blog Entries: 1
Rep Power: 14
dlahaye is on a distinguished road
Thx for your advice on this matter. Best wishes, Domenico.
dlahaye is offline   Reply With Quote

Old   April 28, 2021, 06:36
Default
  #6
New Member
 
Lorenzo
Join Date: Apr 2021
Posts: 9
Rep Power: 3
lollo15 is on a distinguished road
Thanks you for the response,

I have another question, the unity of lewis number means also that the Prt is equal to the Sct ?

Best wishes, Lorenzo.
lollo15 is offline   Reply With Quote

Old   April 28, 2021, 08:15
Default
  #7
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 517
Blog Entries: 1
Rep Power: 14
dlahaye is on a distinguished road
Yes, affermative, it does. See e.g. the book by Versteeg and Malalasekera, the master thesis of Ali Kadar at https://repository.tudelft.nl/island...tion=education or https://en.wikipedia.org/wiki/Lewis_number among other references.
dlahaye is offline   Reply With Quote

Old   April 30, 2021, 08:31
Default
  #8
New Member
 
Lorenzo
Join Date: Apr 2021
Posts: 9
Rep Power: 3
lollo15 is on a distinguished road
Hello

I serch in the code where the Sct is taken equal to the Prt but if the unity Lewis number is imposed the Sct is not considered. Is it mean that Sct is 1?
And if I do not specify the Prt in the momentumTransport file in the KEpsilonCoeffs section, is it take as 1 also ?

Best wishes, Lorenzo.
lollo15 is offline   Reply With Quote

Reply

Tags
edc model, reacting flow, reactingfoam

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
GenerateVolumeMesh Error - Surface Wrapper Self Interacting (?) AndreP STAR-CCM+ 10 August 2, 2018 08:48
Inconsistencies in reading .dat file during run time in new injection model Scram_1 OpenFOAM 0 March 23, 2018 23:29
[mesh manipulation] Mesh Refinement Luiz Eduardo Bittencourt Sampaio (Sampaio) OpenFOAM Meshing & Mesh Conversion 42 January 8, 2017 13:55
Unaligned accesses on IA64 andre OpenFOAM 5 June 23, 2008 11:37
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Meshing & Mesh Conversion 2 July 15, 2005 05:15


All times are GMT -4. The time now is 07:22.