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

problem with DNS turbulent channel flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 4, 2014, 11:18
Default problem with DNS turbulent channel flow
  #1
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 263
Rep Power: 4
huangxianbei is on a distinguished road
Hi, all
I'm doing a research of turbulent channel flow with DNS. I referred to this thread " Everything you need to compute DNS in channel vs OF 2.1.0" by adding a source term in the UEqn representing pressure gradient.
I've uploaded my solver based on icoFoam. The problem I'm facing is the unsteady solving process. That is, as time goes on, the velocity at the centerline will decrease from my desired value, and this exist consistently along the solving process, I don't know how to deal with it.

Oh, my mesh is 64*128*64,the fluid is water, the channel is 12.56*2*6.28m, Ret=194, Re=2900, according to this , the ut=1.95358e-4m/s and the dp/(dx*rho)=3.82e-8m/s^2

The result is big so I can't upload, all the things except it are uploaded here, hope anyone can give me some suggestions.

Thankssolver.tar.gz

0.tar.gz

constant.tar.gz

system.tar.gz
huangxianbei is offline   Reply With Quote

Old   March 4, 2014, 13:48
Default
  #2
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 6
ArathoN is on a distinguished road
I'm also a newbie in OF but I know that for a DNS simulation you need a finer mesh otherwise you'll not compute all the turbulent scale, hence the profile near the wall will be far from the real one. Did you try to refine the mesh and see if the results would improve or not?

From what i understood about OF on channel he will maintains the mass flow constant so the velocity profile is quite important, and depending on how the velocity develops near the wall the value at the center-line will change accordingly.

Unfortunately my knowledge on OF is quite scarce, if you found the culprit of the error please write it.
ArathoN is offline   Reply With Quote

Old   March 4, 2014, 20:05
Default
  #3
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 263
Rep Power: 4
huangxianbei is on a distinguished road
Quote:
Originally Posted by ArathoN View Post
I'm also a newbie in OF but I know that for a DNS simulation you need a finer mesh otherwise you'll not compute all the turbulent scale, hence the profile near the wall will be far from the real one. Did you try to refine the mesh and see if the results would improve or not?

From what i understood about OF on channel he will maintains the mass flow constant so the velocity profile is quite important, and depending on how the velocity develops near the wall the value at the center-line will change accordingly.

Unfortunately my knowledge on OF is quite scarce, if you found the culprit of the error please write it.
Thank you for your reply. I just adopt the same mesh as Kristofferson in his paper about rotating channel which is 128*128*128, while in order to make a quick watch at the calculation ,I chose a coarse mesh in the 2 periodic directions. So I think in the spanwise direction, the mesh is enough. While I will still try to refine the mesh in the spanwise to see if there is any improvement.

Also, I tried to use RAS to do a faster calculation, the residuals are also very large along the calculation as DNS. That's quite annoying~~
huangxianbei is offline   Reply With Quote

Old   March 4, 2014, 21:30
Default
  #4
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 6
ArathoN is on a distinguished road
I'm having problem with the residuals too when I try simulating the backward facing step, in my case they will decrease like normal then suddenly they will increase and will have an oscillatory profile. I couldn't find yet the culprit.

Another thing that can influence your case is the schemes adopted try increasing the order and see if the computed solution change or not. In my case thank to changing schemes and the solvers in fvsolution i could increase the rate of convergence and eliminate the oscillatory behavior at the start of the simulation.
ArathoN is offline   Reply With Quote

Old   March 5, 2014, 00:39
Default
  #5
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: CAU,China
Posts: 263
Rep Power: 4
huangxianbei is on a distinguished road
Quote:
Originally Posted by ArathoN View Post
I'm having problem with the residuals too when I try simulating the backward facing step, in my case they will decrease like normal then suddenly they will increase and will have an oscillatory profile. I couldn't find yet the culprit.

Another thing that can influence your case is the schemes adopted try increasing the order and see if the computed solution change or not. In my case thank to changing schemes and the solvers in fvsolution i could increase the rate of convergence and eliminate the oscillatory behavior at the start of the simulation.
yes, the schemes using is important, there are so many schemes to chose. Now I'm installing the openfoam-1.6-ext in order to use the channelfoam. I want to see if it can do a different solution, hoping a improvement then
huangxianbei 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
How do I prescribe Average Velocity /Total Inlet Flow for a simple channel problem? skuznet OpenFOAM Pre-Processing 3 December 1, 2013 23:54
How to get fully developed turbulent flow in a channel TQIM STAR-CCM+ 2 November 25, 2012 23:37
STAR-CD tutorial 16(LES of turbulent channel flow) yss STAR-CD 1 October 15, 2012 14:19
references for how to maintain a constant flow rate in turbulent channel flow amirrstg Main CFD Forum 0 October 25, 2011 03:17
turbulent channel flow profile otto mierka Main CFD Forum 3 August 23, 2005 11:56


All times are GMT -4. The time now is 14:51.