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

Channel flow DNS

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 20, 2006, 04:49
Default 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,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
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 Crank-Nicholson 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 (current status: http://albertopassalacqua.com/?p=1541)
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   June 20, 2006, 05:45
Default Hello Alberto, What is the
  #2
Member
 
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 17
anne is on a distinguished road
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
anne is offline   Reply With Quote

Old   June 20, 2006, 07:07
Default Hello Anne, Re_tau = 180 al
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
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 (current status: http://albertopassalacqua.com/?p=1541)
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   June 20, 2006, 09:19
Default Hello Alberto, Once the flo
  #4
Member
 
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 17
anne is on a distinguished road
Hello Alberto,

Once the flow is well established (this
depends from the initial condition you use)
normally about 10 "through-flow" 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.
anne is offline   Reply With Quote

Old   June 20, 2006, 09:29
Default Yes, I'm using the same criter
  #5
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
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 (current status: http://albertopassalacqua.com/?p=1541)
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   June 28, 2006, 11:17
Default Hello Alberto, I have now r
  #6
Member
 
anne dejoan
Join Date: Mar 2009
Location: madrid, spain
Posts: 66
Rep Power: 17
anne is on a distinguished road
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
anne is offline   Reply With Quote

Old   June 29, 2006, 04:06
Default Hello Anne, I sent an e-mail
  #7
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
Hello Anne,
I sent an e-mail 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 (current status: http://albertopassalacqua.com/?p=1541)
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   February 20, 2007, 09:27
Default Hi Anne + Alberto, I also w
  #8
Senior Member
 
Jens Klostermann
Join Date: Mar 2009
Posts: 117
Rep Power: 17
jens_klostermann is on a distinguished road
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
jens_klostermann is offline   Reply With Quote

Old   February 20, 2007, 09:51
Default I changed the solver to remove
  #9
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
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 (current status: http://albertopassalacqua.com/?p=1541)
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   February 20, 2007, 10:08
Default Hi Alberto, Ok than I guess
  #10
Senior Member
 
Jens Klostermann
Join Date: Mar 2009
Posts: 117
Rep Power: 17
jens_klostermann is on a distinguished road
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 header-files.

I thought you use B and epsilon vor postprocessing.

Can you send or post your calculation results?

Regards,

Jens
jens_klostermann is offline   Reply With Quote

Old   February 20, 2007, 14:45
Default Exactly. Remove the sgs model,
  #11
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
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 (current status: http://albertopassalacqua.com/?p=1541)
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   February 20, 2007, 16:50
Default Here's the code: Hope i
  #12
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
Here's the code:



Hope it helps :-)
Honda likes this.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
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   February 20, 2007, 16:50
Default Here's the code: http://ww
  #13
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
Here's the code:

channelDNS.tar.gz

Hope it helps :-)
Honda likes this.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
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   March 18, 2010, 05:36
Default channelDNS for OF1.6
  #14
Member
 
Heng Xiao
Join Date: Mar 2009
Location: Zurich, Switzerland
Posts: 58
Rep Power: 17
xiao is on a distinguished road
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:
Originally Posted by alberto View Post
Here's the code:

channelDNS.tar.gz

Hope it helps :-)
xiao is offline   Reply With Quote

Old   June 15, 2011, 03:30
Default
  #15
New Member
 
Join Date: Oct 2010
Posts: 24
Rep Power: 15
LijieNPIC is on a distinguished road
Quote:
Originally Posted by alberto View Post
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 Crank-Nicholson 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
Hi, alberto. After doing DNS, how do you get mean field value and statistics such as Reynolds stress?
LijieNPIC is offline   Reply With Quote

Old   March 26, 2012, 05:17
Default
  #16
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 15
levka is on a distinguished road
Hello,
guys i need this DNSchannel solver for OF 2.1.0.
Regards Lev

Quote:
Originally Posted by xiao View Post
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
levka is offline   Reply With Quote

Old   May 6, 2012, 03:41
Default Everything you need to compute DNS
  #17
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 15
levka is on a distinguished road
Have a look here:"Everything you need to compute DNS"

http://www.cfd-online.com/Forums/ope...tml#post359522

Regards, Lev
levka is offline   Reply With Quote

Old   December 5, 2012, 16:42
Default Channel flow DNS with constant pressure gradient
  #18
New Member
 
subhendu
Join Date: May 2012
Posts: 10
Rep Power: 13
raw17 is on a distinguished road
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=(1-y^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:
Originally Posted by levka View Post
Have a look here:"Everything you need to compute DNS"

http://www.cfd-online.com/Forums/ope...tml#post359522

Regards, Lev
raw17 is offline   Reply With Quote

Old   December 6, 2012, 04:43
Default
  #19
Member
 
Lev
Join Date: Dec 2010
Posts: 31
Rep Power: 15
levka is on a distinguished road
Quote:
Originally Posted by raw17 View Post
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=(1-y^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
Turbulent or laminar state you can manage with two parameters:
1-dp/dx
2-initial 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.
levka is offline   Reply With Quote

Old   December 6, 2012, 05:05
Default
  #20
New Member
 
subhendu
Join Date: May 2012
Posts: 10
Rep Power: 13
raw17 is on a distinguished road
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:
Originally Posted by levka View Post
Turbulent or laminar state you can manage with two parameters:
1-dp/dx
2-initial 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.
raw17 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
Channel flow DNS Abhinav Main CFD Forum 1 April 11, 2013 07:37
maintain a constant flow rate in channel flow ? Lewis Main CFD Forum 2 September 28, 2010 12: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


All times are GMT -4. The time now is 13:42.