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

about programming a new B.C.

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 28, 2012, 05:08
Default about programming a new B.C.
  #1
New Member
 
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 5
garrison is on a distinguished road
Hi, I'm trying to programming a new LES wall condition, It's formula is exactly the same as RAS wall function omegaWallFunction.

you may find it src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions

This wall functions needs y+(yplus), but LESModel has not defined y+, I have to comply a new LESmodel using the y+ from RASModel (let's call it myLESModel, and I got a library called myLESModel.so), and after that, when complying the omegaWallfunction for LES, how can I including the head file?

for example, in RAS-omegaWallFunction, it begins with:
#include "omegaWallFunctionFvPatchScalarField.H"
#include "RASModel.H"
#include "fvPatchFieldMapper.H"
#include "volFields.H"
...
now I want to change RASModel.H to myLESModel.H, but it doesn't work..

and another question is, when i compile the new wall function, what should I write in make/file?
garrison is offline   Reply With Quote

Old   February 28, 2012, 06:22
Default
  #2
PeS
New Member
 
Petr Simanek
Join Date: Jan 2011
Location: Prague
Posts: 2
Rep Power: 0
PeS is on a distinguished road
There is an easy (but uncool) way how to do it. Download OpenFOAM 1.5, go to kOmegaSST source and you will find several *.H files. These files are included in kOmegaSST.C and work as omegaWallFunction so you can include them in the very same (there may be some changes in the code, depending on your OF version) way in your code. It is not nice... but easy.

Are you working on the SST-DDES model? Which version? With F1, F2 or the new rd blending? Do you intend to use some nuSgs wall function as well? And what is your application?

Good luck.
Petr
PeS is offline   Reply With Quote

Old   February 28, 2012, 08:34
Default
  #3
New Member
 
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 5
garrison is on a distinguished road
Hi Petr, thanks for your suggestion. it sounds like kOmegaSST in openFoam 1.5 will act in the wall boundary like what omegawallfunction does, am I right? I have already change the omegawallfunction(RANS) to LES, but the problem is the function uses ypls which is not included in LESModel, so anyway I have to compile a new LESModel(say myLESModel).

Indeed I am working on SST-DDES with wing tip vortex, with omegafunction i think maybe i can use some relatively coarse grid to get the correct results.

I have no idea about F1,F2 or rd version, I am using the standard MSST model and I found the code in RANS-kOmegaSST and changed it into a DES one. the Delay blending function is with F2, i find the function in F.R. Menter's paper, 2004.

would you pls tell me what's the new rd bending? is it the same rd from SA-DDES?

Regards
Garriosn
garrison is offline   Reply With Quote

Old   February 28, 2012, 08:57
Default
  #4
PeS
New Member
 
Petr Simanek
Join Date: Jan 2011
Location: Prague
Posts: 2
Rep Power: 0
PeS is on a distinguished road
Hi, the rd function is very similar to SA - check this article:

Development of DDES and IDDES Formulations for the k-ω Shear Stress Transport Model
Mikhail S. Gritskevich,Andrey V. Garbaruk,Jochen Schütze, Florian R. Menter


If you include the .H files in your DES model you don't need to compile new LESModel, only your DES model. The included files will work as omegaWallFunction and the ypls is computed in the included file kOmegaWallFunctionsI.H .
PeS is offline   Reply With Quote

Old   February 28, 2012, 09:39
Default
  #5
New Member
 
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 5
garrison is on a distinguished road
Thanks Petr, it sounds very attractive! just one more question, with your method, how should i set the B.C. for omega? should it be "caldulated"?

and do you have those files you mentioned above? would you please be so kind to send them to me ? it can really do me a huge favour! but anyway, thanks a million.

p.s. my email adress: garrison.liang@gmail.com
garrison is offline   Reply With Quote

Old   February 29, 2012, 08:45
Default
  #6
New Member
 
Garrison, Liang
Join Date: Feb 2012
Posts: 15
Rep Power: 5
garrison is on a distinguished road
Hi Petr, I already found it. And i think i know how to deal with it like you told me.
once again, thanks a million!
garrison is offline   Reply With Quote

Reply

Tags
compile, wall function

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
Confusion between 'non-slip' and 'slip' B.C jinwon park Main CFD Forum 7 March 6, 2008 05:56
how to set such B.C. in Fluent? Sherri FLUENT 1 January 12, 2006 13:32
Urgent! Help on UDF to set B.C. of 3rd type Ray Hong FLUENT 0 December 28, 2005 20:35
How to set B.C. of the 2nd or 3rd type in UDS? Ray Hong FLUENT 0 December 28, 2005 07:03
Non-Reflecting B.C. in NSC2KE Zou Chu Main CFD Forum 2 May 27, 1999 21:26


All times are GMT -4. The time now is 18:44.