CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Epsilon Boundary condition

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

Like Tree3Likes
  • 1 Post By Bernhard
  • 1 Post By MaLa
  • 1 Post By MaLa

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 3, 2014, 11:05
Default Epsilon Boundary condition
  #1
Member
 
Manan
Join Date: Oct 2013
Location: Göteborg
Posts: 37
Rep Power: 12
MaLa is on a distinguished road
Hello

I am trying to validate a new turbulence model (zeta f) and am solving for epsilon amongst the other variables. The epsilon equation that will be used is the one used in standard k-epsilon model.

I want to make use of the formulation, epsilon @Y=0 is specified as 2*nu* sqr(d/dy(k^0.5)). (Turbulence models for Near Wall and Low Reynolds Number Flows: A Review, Patel et. al 1985)

I don't understand though how I should specify this in OF? Epsilon @ the wall i.e. at Y=0, is a derivative of k, that is: it is using the value of k @wall and the the value of k @the first node.
How do I instruct the solver to use these two values?

Please let me know if I haven't made myself clear. Thanks.
MaLa is offline   Reply With Quote

Old   February 3, 2014, 11:43
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21
Bernhard is on a distinguished road
Are you aware that the zeta-f0 model was already implemented in OpenFOAM by Mirza Popovac? See: http://openfoamwiki.net/index.php/Si...bulence_models

Also, if you use the epsilon equation from the standard k-epsilon, you don't include any Low-Re modifications, which you would need. Maybe you want to start from the

To answer your question: your best change is to start from the epsilonWallFunction, see how the loop over the faces is done and access of the field data, and adjust to your needs
Ramzy1990 likes this.
Bernhard is offline   Reply With Quote

Old   February 11, 2014, 07:48
Default
  #3
Member
 
Manan
Join Date: Oct 2013
Location: Göteborg
Posts: 37
Rep Power: 12
MaLa is on a distinguished road
Hi Bernhard

Thanks so much for your reply. I wasn't aware of the existence of the zeta-f formulation and had written my own code. Fortunately I have Popovac's code to compare it with now.
Also, the tip about modifying the epsilonWallFunction is helpful. It did serve my need.
Ramzy1990 likes this.
MaLa is offline   Reply With Quote

Old   February 11, 2014, 07:50
Default
  #4
Member
 
Manan
Join Date: Oct 2013
Location: Göteborg
Posts: 37
Rep Power: 12
MaLa is on a distinguished road
One question though, when I try compiling Popovac's formulation of zeta-f, using wmake libso, I get the following error message:

/chalmers/sw/unsup64/OpenFOAM/OpenFOAM-1.6-ext/wmake/MakefileOptions:37: /chalmers/sw/unsup64/OpenFOAM/OpenFOAM-1.6-ext/wmake/rules/linux64Gcc47/general: No such file or directory
make: *** No rule to make target `/chalmers/sw/unsup64/OpenFOAM/OpenFOAM-1.6-ext/wmake/rules/linux64Gcc47/general'. Stop.
/chalmers/sw/unsup64/OpenFOAM/OpenFOAM-1.6-ext/wmake/MakefileFiles:39: /chalmers/sw/unsup64/OpenFOAM/OpenFOAM-1.6-ext/wmake/rules/linux64Gcc47/general: No such file or directory
/chalmers/sw/unsup64/OpenFOAM/OpenFOAM-1.6-ext/wmake/MakefileFiles:40: linux64Gcc47DPOpt/options: No such file or directory
make: *** No rule to make target `linux64Gcc47DPOpt/options'. Stop.
wmake error: file 'Make/linux64Gcc47DPOpt/objectFiles' could not be created


I figured out that there are only a few compilers available in the OF-1.6-ext folder (Gcc45, Gcc46) but not Gcc47. Is there any way to direct $WM_COMPILER to seek Gcc46 (which is available), instead of adding the folder for Gcc 47 in the OF1.6-ext wmake folder explicitly?
MaLa is offline   Reply With Quote

Old   February 11, 2014, 08:12
Default
  #5
Member
 
Manan
Join Date: Oct 2013
Location: Göteborg
Posts: 37
Rep Power: 12
MaLa is on a distinguished road
Figured it out. It was as simple as saying:

export WM_COMPILER = Gcc46

Thanks!
Ramzy1990 likes this.
MaLa 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
Mathematical representation of fixedDisplacementZeroShear boundary condition Sargam05 OpenFOAM 14 January 11, 2022 07:55
Wind turbine simulation Saturn CFX 58 July 3, 2020 02:13
No-slip condition for non-resolved boundary layer in open channel banks Lupocci Main CFD Forum 1 January 17, 2013 04:11
External Radiation Boundary Condition for Grid Interface CFD XUE FLUENT 0 July 9, 2010 03:53
External Radiation Boundary Condition (Two sided wall), Grid Interface CFD XUE FLUENT 0 July 8, 2010 07:49


All times are GMT -4. The time now is 03:02.