CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Shear stress from velocity profile

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

LinkBack Thread Tools Display Modes
Old   November 12, 2009, 07:46
Default Shear stress from velocity profile
New Member
Join Date: Sep 2009
Posts: 3
Rep Power: 7
baranii2000 is on a distinguished road
Dear all,

I was trying to calculate the shear stress due to fluid flow inbetween two parallel plates from the velocity profile taken from the middle of the flow channel (along vertical axis 'y') using comsol. The solution i got is two times the analytical solution i have calculated. The analytical equation i have used here is

where tou is shear stress, umax is the max velocity
h is half of the distance between plates.

The calculation with velocity profile from comsol is done by

Tou= viscosity*(du/dy)
du/dy, i calculated with excel from the u,y values in the velocity profile.

the velocity profile is parabolic and looks laminar.
is that i am doing something wrong?


baranii2000 is offline   Reply With Quote

Old   November 12, 2009, 08:04
Senior Member
sbaffini's Avatar
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 534
Blog Entries: 14
Rep Power: 17
sbaffini will become famous soon enough
I don't remember the formulas but, if i understand well your problem, it is impossible that for the parabolic velocity profile the derivative at wall is equal to the one from the linear profile. Hence you should check your analytical formula for the stress (just take by yourself the derivative of the analytic velocity profile)
sbaffini is offline   Reply With Quote

Old   November 16, 2009, 21:16
New Member
Charles Rundle
Join Date: Oct 2009
Posts: 9
Rep Power: 7
Charles Rundle is on a distinguished road

I agree with the previous poster. I think your analytical solution is incorrect. Using a back of the envelope calculation I get a significantly different solution for the shear stress. I am assuming a laminar parabolic flow so depending on conditions my solution may not apply. Also since I mean the back of envelope quite literally I may have made an error.

Assuming your analytical solution is correct I would consider using a finer mesh. I would suggest that you plot your velocity profile against the analytical solution. The profile can look parabolic without truly being parabolic. The problem may simply be that grid resolution or insufficient convergence.

Lastly, your analytical solution may not valid for your Re.

I hope this is helpful. Have a good day.

Charles Rundle is offline   Reply With Quote

Old   November 18, 2009, 10:20
New Member
Join Date: Sep 2009
Posts: 3
Rep Power: 7
baranii2000 is on a distinguished road
Thanks very much Paolo & Charles, i think i did a mistake in my analytical is supposed to be Tou=2*Umax/h*viscosity i dont know in many books they just refer to Umax/h...

baranii2000 is offline   Reply With Quote


comsol, fluid flow, shear stress

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
Velocity Profile Jeff FLUENT 1 November 24, 2008 09:21
LES Near-Wall Velocity Profile Alex May FLUENT 2 June 12, 2008 12:36
Macro to set Wall Shear Stress Satish FLUENT 4 November 26, 2003 15:46
Re: About wall shear stress Mike FLUENT 9 November 17, 2003 15:41
How to get shear stress all over the flow domain lingo FLUENT 2 November 3, 2003 22:16

All times are GMT -4. The time now is 09:34.