
[Sponsors] 
June 20, 2006, 03:49 
Hello,
I'm doing a DNS in a c

#1 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
Hello,
I'm doing a DNS in a channel flow. I modified the channelOodles solver removing the LES model and set up the case using the same numerical schemes adopted in dnsFoam (cubic method for U divergence term, linear for all other terms). The only difference is that I adopted the CrankNicholson scheme for time integration, while dnsFoam adopts the backward Euler. Results are quite in good agreement with experimental data, if I consider the mean velocity. But I've a significative underestimation of the u_rms (flow direction) and of the Reynolds stress u'v', if compared with the results of Kim et al. (1987). Could be this related to the choice of numerical schemes? Thanks in advance, Alberto
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

June 20, 2006, 04:45 
Hello Alberto,
What is the

#2 
Member
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 10 
Hello Alberto,
What is the Reynolds number you are using and the grid resolution in terms of wall coordinates? I am also doing some channel flow but I still need some statistical convergence. The schemes I have applied are second order in time and space, and I also modified the LES channel to make it as a DNS. The Reyn I compute is Retau=180 and the grid is almost a DNS: 100by100by64 for Lx=8h, Ly=2h and Lz=4h. Once I have some results I will communicate them to you. Anne 

June 20, 2006, 06:07 
Hello Anne,
Re_tau = 180 al

#3 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
Hello Anne,
Re_tau = 180 also for me, so Ubar = 15.66. I'm using a channel with these dimensions: Lx = 2pi*delta, Ly = 2*delta, Lz = pi*delta, where delta is half the channel height, and is taken equal to 1. The grid is: Nx = 96, Ny = 129, Nz = 64, uniform along x and z. Along y I adopted an hyperbolic tangent distribution of the nodes. I'll tell you my results too. How many times are you considering for statistics? With kind regards, Alberto
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

June 20, 2006, 08:19 
Hello Alberto,
Once the flo

#4 
Member
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 10 
Hello Alberto,
Once the flow is well established (this depends from the initial condition you use) normally about 10 "throughflow" units time is enough. But check by simple visualization of the velocity field that you flow looks well turbulent before starting statistics. The best way to check your statistic convergence is to draw the profile of the total mean shear and ensure that it is linear. Anne Apparently you grid resolutiojn should be enough. I will let you know my results. 

June 20, 2006, 08:29 
Yes, I'm using the same criter

#5 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
Yes, I'm using the same criterion for the statistics convergence.
My initial condition is a fully developed flow obtained by a previous DNS performed with spectral methods, so I've neglected the first two dimensionless times and then I started to average. Alberto
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

June 28, 2006, 10:17 
Hello Alberto,
I have now r

#6 
Member
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 10 
Hello Alberto,
I have now results of the channel flow I simulated with OpenFoam. I also added a passive scalar transport equation in a channel with fixed temperature at the bottom (T=0) and top (T=1) walls of the channel. The DNS (or better say quasiDNS because of my resolution) was performed with second order time scheme (backward) and second order spatial schemes (gauss linear for all the gradients: Give me your email address so that I send you the plots where you will also find more details about the grid resolution. Anne 

June 29, 2006, 03:06 
Hello Anne,
I sent an email

#7 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
Hello Anne,
I sent an email to you with my address. Thanks, Alberto
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

February 20, 2007, 09:27 
Hi Anne + Alberto,
I also w

#8 
Senior Member
Jens Klostermann
Join Date: Mar 2009
Posts: 117
Rep Power: 10 
Hi Anne + Alberto,
I also want to simulate the dnsChannel. How did you B and epsilon declare? (in createAverage.h) Could you please sent your results with the dnsChannel over (email adress in the profile)? regards, Jens 

February 20, 2007, 09:51 
I changed the solver to remove

#9 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
I changed the solver to remove the LES model, so I haven't B and epsilon, which would not be used in a DNS.
I'll look for the code and publish it. Regards, Alberto
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

February 20, 2007, 10:08 
Hi Alberto,
Ok than I guess

#10 
Senior Member
Jens Klostermann
Join Date: Mar 2009
Posts: 117
Rep Power: 10 
Hi Alberto,
Ok than I guess I know what to do: Just change + sgsModel>divB(U) with fvm::laplacian(nu, U) and delete sgsModel>correct() in channelOodles and some changes in the headerfiles. I thought you use B and epsilon vor postprocessing. Can you send or post your calculation results? Regards, Jens 

February 20, 2007, 14:45 
Exactly. Remove the sgs model,

#11 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
Exactly. Remove the sgs model, B and epsilon costructors, and averages calculation. That's all :)
Regards, A.
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

February 20, 2007, 16:50 
Here's the code:
Hope i

#12 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
Here's the code:
Hope it helps :)
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

February 20, 2007, 16:50 
Here's the code:
http://ww

#13 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
__________________
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 opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

March 18, 2010, 05:36 
channelDNS for OF1.6

#14  
Member
Heng Xiao
Join Date: Mar 2009
Location: Zurich, Switzerland
Posts: 58
Rep Power: 10 
Thanks, alberto, for sharing the code. I modified a bit (just removed the averaging and probing part, which are not necessary anymore for OF1.6), and it works perfectly!
Since my modification part is so trivial, I won't post the code here. If anyone needs it, let me know. best, Heng Quote:


June 15, 2011, 02:30 

#15  
New Member
Join Date: Oct 2010
Posts: 24
Rep Power: 8 
Quote:


March 26, 2012, 04:17 

#16  
Member
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 8 
Hello,
guys i need this DNSchannel solver for OF 2.1.0. Regards Lev Quote:


May 6, 2012, 02:41 
Everything you need to compute DNS

#17 
Member
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 8 
Have a look here:"Everything you need to compute DNS"
http://www.cfdonline.com/Forums/ope...tml#post359522 Regards, Lev 

December 5, 2012, 16:42 
Channel flow DNS with constant pressure gradient

#18  
New Member
subhendu
Join Date: May 2012
Posts: 10
Rep Power: 7 
Hello
I am trying to do channel flow simulation with constant pressure gradient. I am giving an initial streamwise vortices and a sinous pertubation in W along with U=(1y^2) to generate my initial velocity field in 0/ folder using funkysetfield. I have done the same simulation with exactly the same parameters in another spectral code written in fortran. My initial conditions are excatly the same in both the cases .But in OPENFOAM dns case no matter how small amplitude of the initial pertubation I give the pertubation energy is always increasing and the solution is becoming turbulent . Can any help to solve this problem with constant pressure gradient case . I want to keep constant pressure gradient because I want to validate OpenFoam result with my Fortan code results. Any kind of help would be highly appreciated thanks Quote:


December 6, 2012, 04:43 

#19  
Member
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 8 
Quote:
1dp/dx 2initial perturbations If dp/dx small enough (when Re<<5000)then any perturbations will vanish in time. If Re>>5000 then no matter value of perturbations the flow will develop to turbulent state If you deal with transient Re around 5000 (that you have defined by dp/dx) then perturbation will play major role. And the flow can reach either turb/ or lam. states depending on initial perturbations. dx/dx/rho is defined in transportProperties. The value you can estimate in Excel (dp/dx/rho VIA Re_tau) file that is included in the package. 

December 6, 2012, 05:05 

#20  
New Member
subhendu
Join Date: May 2012
Posts: 10
Rep Power: 7 
Hi
Thanks for your reply. Yes you are right my Re =4000 based on Uc*h/nu . Where Uc is the centreline velocity =1 , h=1 so it means that my Re=1/nu. I have defined my gradP =2/Re . Now my problem as I said before is that even when my pertubations are very small my solution is never becoming laminar. Even when the initial pertubatiojn energy is very small the solution is still becoming turbulent and these results are not matching with my previous findings. Quote:


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Channel flow DNS  Abhinav  Main CFD Forum  1  April 11, 2013 06:37 
maintain a constant flow rate in channel flow ?  Lewis  Main CFD Forum  2  September 28, 2010 11:35 
Channel flow set up  Jung  Main CFD Forum  2  November 21, 2007 07:51 
channel flow  Khan  Main CFD Forum  0  February 3, 2007 11:39 
channel flow help  Tajul  FLUENT  4  February 22, 2006 21:23 