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

Everything you need to compute DNS in channel vs OF 2.1.0

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

Like Tree14Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   May 6, 2012, 02:39
Default Everything you need to compute DNS in channel vs OF 2.1.0
  #1
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 7
levka is on a distinguished road
Hello everyone.
Here are everything you need to compute DNS in channel vs OF 2.1.
Attached pdf is comparison this solver with DNS data of Kawamura and "Kim and Moin, 1987". Max divergence by mean "u" vs Moin is 2%, vs Kawamura 0%. Also in pdf there is view of full turbulent velocity profile that was taken during average procedure(10,000-30,000).
Following reference is to download zip file:
http://www.sendspace.com/file/xw4myd

Zip contains:
- solver (modified ico_Foam);
- utility to create tangential mesh spacing;
- test case with perturbed initial data;
- computed averaged data+log file;
- manipulation with data Excel file;
- paper Kim and Moin 1987;
- DNS data from paper above;
- DNS data Kawamura;
- master thesis of Steven van Haren"Testing DNS capability of OpenFOAM
and STAR-CCM+";

In case zip is not available i will e-mail it .


Have fun guys!
Attached Files
File Type: pdf Compare.pdf (51.7 KB, 509 views)
ehsan, Bernhard, eRzBeNgEl and 7 others like this.
levka is offline   Reply With Quote

Old   May 7, 2012, 11:22
Default
  #2
Senior Member
 
Goutam Saha
Join Date: Dec 2011
Location: UK
Posts: 131
Rep Power: 6
Goutam is on a distinguished road
Thanks for your contribution. What about LES for a pipe flow?
Goutam is offline   Reply With Quote

Old   May 7, 2012, 16:18
Default
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,907
Rep Power: 27
alberto will become famous soon enoughalberto will become famous soon enough
You might want to consider to submit this validation study to the Foundation, once they will start accepting code and contributions. It is surely valuable.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats.
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   May 8, 2012, 02:40
Default nice
  #4
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 7
levka is on a distinguished road
Nice idea! But how to do it ?
levka is offline   Reply With Quote

Old   May 8, 2012, 02:46
Default
  #5
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,907
Rep Power: 27
alberto will become famous soon enoughalberto will become famous soon enough
There is no defined procedure yet. When I tried to submit some code I used OpenCFD mantis, and they said they are defining the rules for contributors, and will contact me when ready. Maybe you could ask them if they are interested also in validation studies.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats.
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   May 23, 2012, 14:54
Default DNS and Airfoil
  #6
Senior Member
 
ehsan
Join Date: Mar 2009
Posts: 106
Rep Power: 9
ehsan is on a distinguished road
Dear all

Could someone help me in taking a test of DNS over airfoil?

Thanks
ehsan is offline   Reply With Quote

Old   May 24, 2012, 02:02
Default ...
  #7
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 7
levka is on a distinguished road
Hello ehsan,

I suggest you to download solver solver from the first post and utilize that for your purpose. Just change mesh and BC according to your case, not forget to use initial conditions with perturbations( you can map them from my case to yours using mapField utility). And enjoy. Be in touch. I am interesting in the results.
songwukong likes this.
levka is offline   Reply With Quote

Old   May 24, 2012, 02:32
Default
  #8
Senior Member
 
ehsan
Join Date: Mar 2009
Posts: 106
Rep Power: 9
ehsan is on a distinguished road
Dear Levka

Thanks a lot, we will work on this topic and will be in touch. However, does 2D DNS solution is very time consuming?

Regards
Ehsan
ehsan is offline   Reply With Quote

Old   May 24, 2012, 03:54
Default ...
  #9
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 7
levka is on a distinguished road
Quote:
Originally Posted by ehsan View Post
Dear Levka

Thanks a lot, we will work on this topic and will be in touch. However, does 2D DNS solution is very time consuming?

Regards
Ehsan
That depends on your mesh and mesh it self should be defined according to comparison some integral parameter (for example U+ profile or Cf, Cd ) with some closest benchmark problem. Got it?
Cheers
levka is offline   Reply With Quote

Old   May 24, 2012, 04:58
Default
  #10
Senior Member
 
ehsan
Join Date: Mar 2009
Posts: 106
Rep Power: 9
ehsan is on a distinguished road
Dear Levka

In fact, I recently saw this paper:
http://www.sciencedirect.com/science...45793004001203

where the flow simulation over two dimensional and three dimensional airfoil had been reported, in Fig. 7 the U velocity had been reported and we will try to compare our results with this paper.

In fact, I am not familiar with DNS at all while I had worked on LES before, so it took me some time to understand DNS. I downloaded the file you prepared and work on it now.

Regards and thanks
Ehsan
ehsan is offline   Reply With Quote

Old   May 26, 2012, 15:34
Default Run b y using dnsFoam
  #11
New Member
 
Farshad Rezaei
Join Date: Apr 2012
Posts: 10
Rep Power: 6
Farshad585 is on a distinguished road
Hi everybody
I want to simulate an airfoil and analyze it by dnsFoam. I've used boundary condition as it is in incompressible/simplefoam/2Dairfoil
I've started to run by dnsFoam, but there are the below errors:
--> FOAM FATAL ERROR:
calculated number of cells is incorrect

From function Kmesh::Kmesh(const fvMesh& mesh)
in file Kmesh/Kmesh.C at line 84.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Foam::Kmesh::Kmesh(Foam::fvMesh const&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/librandomProcesses.so"
#3
in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/dnsFoam"
#4 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#5
in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/dnsFoam"
Aborted


Could anybody help me what should I do?
Farshad585 is offline   Reply With Quote

Old   May 28, 2012, 10:08
Default ....
  #12
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 7
levka is on a distinguished road
Hello,Farshad585
that error maybe coursed by either you did not run blockmesh, either your initial field is generated not for your mesh.
Best

Quote:
Originally Posted by Farshad585 View Post
Hi everybody
I want to simulate an airfoil and analyze it by dnsFoam. I've used boundary condition as it is in incompressible/simplefoam/2Dairfoil
I've started to run by dnsFoam, but there are the below errors:
--> FOAM FATAL ERROR:
calculated number of cells is incorrect

From function Kmesh::Kmesh(const fvMesh& mesh)
in file Kmesh/Kmesh.C at line 84.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Foam::Kmesh::Kmesh(Foam::fvMesh const&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/librandomProcesses.so"
#3
in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/dnsFoam"
#4 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#5
in "/opt/openfoam210/platforms/linuxGccDPOpt/bin/dnsFoam"
Aborted


Could anybody help me what should I do?
Farshad585 likes this.
levka is offline   Reply With Quote

Old   May 28, 2012, 18:38
Default
  #13
New Member
 
Farshad Rezaei
Join Date: Apr 2012
Posts: 10
Rep Power: 6
Farshad585 is on a distinguished road
Quote:
Originally Posted by levka View Post
Hello,Farshad585
that error maybe coursed by either you did not run blockmesh, either your initial field is generated not for your mesh.
Best
Thanks a lot
Farshad585 is offline   Reply With Quote

Old   November 13, 2012, 10:30
Default
  #14
Senior Member
 
aerospain
Join Date: Sep 2009
Location: Madrid, Spain
Posts: 125
Rep Power: 9
aerospain is on a distinguished road
Send a message via Skype™ to aerospain
Hi levka,

Could you explain how to proceed with the zip file you made available in order to build the ico_DNS and run it, please?

thanks!
aerospain is offline   Reply With Quote

Old   December 26, 2012, 06:56
Default least sqaures & fourth order
  #15
New Member
 
Farshad Rezaei
Join Date: Apr 2012
Posts: 10
Rep Power: 6
Farshad585 is on a distinguished road
Hi everybody
I've used least squares and fourth order for discretization for pressure gradient (grad(p)), but I don't know the performance of least sqaures and fourth order. Could somebody give me some information about it?
regards
Farshad585 is offline   Reply With Quote

Old   June 21, 2013, 06:07
Default
  #16
New Member
 
Join Date: Jun 2013
Posts: 1
Rep Power: 0
student_07 is on a distinguished road
Hallo Levka,



I have one question to your solver. I'm still a beginner to OpenFoam, so don't wonder, if it is obvious to you.



I ran the test file and it works fine. Now I want to modify it and want to use a constant inlet velocity as boundry condition. In "constant/transportProperties" you set a pressure gradient "gradP". But there is also a "transportProperties_channel" file in which U_bar (probably an average velocity) can be set instead of the pressure gradient.


How do I have to change the solver to use a constant inlet velocity? When I just use Ubar, OpenFoam still requires a pressure gradient?


Regards,


Klaus
student_07 is offline   Reply With Quote

Old   June 21, 2013, 14:09
Default
  #17
New Member
 
Farshad Rezaei
Join Date: Apr 2012
Posts: 10
Rep Power: 6
Farshad585 is on a distinguished road
Quote:
Originally Posted by student_07 View Post
Hallo Levka,



I have one question to your solver. I'm still a beginner to OpenFoam, so don't wonder, if it is obvious to you.



I ran the test file and it works fine. Now I want to modify it and want to use a constant inlet velocity as boundry condition. In "constant/transportProperties" you set a pressure gradient "gradP". But there is also a "transportProperties_channel" file in which U_bar (probably an average velocity) can be set instead of the pressure gradient.


How do I have to change the solver to use a constant inlet velocity? When I just use Ubar, OpenFoam still requires a pressure gradient?


Regards,


Klaus
hi
you shoul go to system/fvscheme and set pressure gradient as it show gradP
Farshad585 is offline   Reply With Quote

Old   June 21, 2013, 14:10
Default
  #18
New Member
 
Farshad Rezaei
Join Date: Apr 2012
Posts: 10
Rep Power: 6
Farshad585 is on a distinguished road
hi
you shoul go to system/fvscheme and set pressure gradient as it show gradP
Farshad585 is offline   Reply With Quote

Old   June 23, 2013, 07:16
Default
  #19
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 7
levka is on a distinguished road
I suggest you to look at regular channel solver and it's tutorial. "Ubar" is used over there without any gradP.
levka is offline   Reply With Quote

Old   July 9, 2013, 05:18
Default
  #20
New Member
 
Join Date: Nov 2012
Posts: 13
Rep Power: 5
lixx is on a distinguished road
Nice work!

In the test folder, I can not find the 0/ folder. So you just started from a precursor run or else?

Quote:
Originally Posted by levka View Post
Hello everyone.
Here are everything you need to compute DNS in channel vs OF 2.1.
Attached pdf is comparison this solver with DNS data of Kawamura and "Kim and Moin, 1987". Max divergence by mean "u" vs Moin is 2%, vs Kawamura 0%. Also in pdf there is view of full turbulent velocity profile that was taken during average procedure(10,000-30,000).
Following reference is to download zip file:
http://www.sendspace.com/file/xw4myd

Zip contains:
- solver (modified ico_Foam);
- utility to create tangential mesh spacing;
- test case with perturbed initial data;
- computed averaged data+log file;
- manipulation with data Excel file;
- paper Kim and Moin 1987;
- DNS data from paper above;
- DNS data Kawamura;
- master thesis of Steven van Haren"Testing DNS capability of OpenFOAM
and STAR-CCM+";

In case zip is not available i will e-mail it .


Have fun guys!
Farshad585 likes this.
lixx 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
Channel flow DNS alberto OpenFOAM Running, Solving & CFD 20 February 17, 2016 18:26
DNS data of Jimenez et al. for turbulent channel flow @ Ret = 2000 sbaffini Main CFD Forum 2 July 29, 2013 06:03
Channel flow DNS Abhinav Main CFD Forum 1 April 11, 2013 06:37
DNS: how to compute nonlinear term in the spectral kinetic energy equation magicsquirrel Main CFD Forum 4 December 13, 2011 18:45
DNS of turbulence in a channel flow YANGLIANG OpenFOAM 0 March 4, 2010 09:40


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