CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

2d turbulent flow parallel plates matlab

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 8, 2015, 19:07
Default 2d turbulent flow parallel plates matlab
  #1
P_P
New Member
 
Join Date: Mar 2014
Posts: 11
Rep Power: 12
P_P is on a distinguished road
Hello,

I am currently trying to determine the turbulent velocity profile with matlab for a flow between two parallel plates. In order to do so, I am using the universal velocity profiles of Karman and Taylor/Prandtl. The veloityprofile is considered to be symmetric.
In both cases the velocity gradients close to the walls show realistic behaviour. Still I have my doubts concerning the values on the centerline where the velocity profile shows a peak. Does anyone know models which represent a smoother transition on the centerline or rather a better approximation for high y+ numbers?
A plot of the velocity profile and the corresponding Matlabcode can be found in the appendix.
Every input is welcome

Greetings

P_P
Attached Images
File Type: jpg Turbulent_Velocity_profile.jpg (13.7 KB, 36 views)
Attached Files
File Type: zip Turbulenz.zip (2.2 KB, 35 views)
P_P is offline   Reply With Quote

Old   June 9, 2015, 03:57
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,810
Rep Power: 72
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by P_P View Post
Hello,

I am currently trying to determine the turbulent velocity profile with matlab for a flow between two parallel plates. In order to do so, I am using the universal velocity profiles of Karman and Taylor/Prandtl. The veloityprofile is considered to be symmetric.
In both cases the velocity gradients close to the walls show realistic behaviour. Still I have my doubts concerning the values on the centerline where the velocity profile shows a peak. Does anyone know models which represent a smoother transition on the centerline or rather a better approximation for high y+ numbers?
A plot of the velocity profile and the corresponding Matlabcode can be found in the appendix.
Every input is welcome

Greetings

P_P

Hello,
I suggest plotting the velocity u+ versus y+ so that you can better check for viscous sub-layer and log region
FMDenaro is offline   Reply With Quote

Old   June 9, 2015, 06:57
Default
  #3
P_P
New Member
 
Join Date: Mar 2014
Posts: 11
Rep Power: 12
P_P is on a distinguished road
Hello FMDenaro,

thanks a lot for your fast reply. I have plotted u+ versus y+. The linear behaviour of the laminar sublayer can be seen for y+<11.5.
You can see the the plot in the appendix.

Greetings P_P
Attached Images
File Type: jpg u+_y+.jpg (12.7 KB, 21 views)
P_P is offline   Reply With Quote

Old   June 9, 2015, 07:00
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,810
Rep Power: 72
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by P_P View Post
Hello FMDenaro,

thanks a lot for your fast reply. I have plotted u+ versus y+. The linear behaviour of the laminar sublayer can be seen for y+<11.5.
You can see the the plot in the appendix.

Greetings P_P

to check further, you can superimpose the law u+=y+ and the log law using a logarithmic scale
FMDenaro is offline   Reply With Quote

Old   June 9, 2015, 09:42
Default
  #5
P_P
New Member
 
Join Date: Mar 2014
Posts: 11
Rep Power: 12
P_P is on a distinguished road
I don't know if I unterstood you well. I have created a semilogarithmic plot u+ versus y+. What is the idea behind checking the laminar sublayer for my problem?
Attached Images
File Type: jpg u_plus_log(y_plus).jpg (13.9 KB, 20 views)
P_P is offline   Reply With Quote

Old   June 9, 2015, 11:30
Default
  #6
Super Moderator
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,412
Rep Power: 49
flotus1 has a spectacular aura aboutflotus1 has a spectacular aura about
Looks reasonable.
The "peak" in your initial figure is due to the fact that the equations you used are for the flow over a flat plate with an infinite fluid domain. Not for a fluid between parallel plates.
flotus1 is offline   Reply With Quote

Old   June 9, 2015, 12:11
Default
  #7
P_P
New Member
 
Join Date: Mar 2014
Posts: 11
Rep Power: 12
P_P is on a distinguished road
Hello flotus1,
thanks for your reply. Do you know any models which offer a better desciription of the velocity profile between two parallel plates? Or is it unevitable to excute a CFD simulation with common CFD software in order to export the velocity profile.


P_P
P_P is offline   Reply With Quote

Old   June 9, 2015, 14:28
Default
  #8
Super Moderator
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 3,412
Rep Power: 49
flotus1 has a spectacular aura aboutflotus1 has a spectacular aura about
I havent come across any analytical formula for this case.
However there is plenty of DNS data available for turbulent flow between 2 parallel flat plates. Here is just one example...
flotus1 is offline   Reply With Quote

Reply

Tags
log law, matlab, turbulent flow

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
CFX Treatment of Laminar and Turbulent Flows Jade M CFX 18 September 15, 2022 07:08
Nusselt Number for flow between flat parallel plates vijud25 FLUENT 2 November 23, 2015 12:29
turbulent or laminar flow student2008 Main CFD Forum 5 August 24, 2013 05:33
Turbulent Flow in Straight Square Duct - k-ep Validation mcclud STAR-CCM+ 2 August 7, 2012 16:56
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 12:19


All times are GMT -4. The time now is 16:06.