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

Flow Rate from Turbulent Boundary Conditions

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 18, 2020, 12:06
Default Flow Rate from Turbulent Boundary Conditions
  #1
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
I did a Civil Engineering course some years ago and from my textbook I had a question on Turbulent boundary Flow.

In case you are wondering where my calculations are, I am actually blocked on this question and not sure where to start.

May someone please explain to me what things I need to understand, so that I can work towards the answer and get a better understanding on this ?



Take air density = 1.18 kg/m^3 and viscosity = 1.3 * 10^-5 Pa.s

Book Answer 1.31 kg/s


They say calculate mass flow rate entering air intake under zero pressure gradient conditions, how does zero pressure gradient come into this?

From my textbook it says
mass flow rate m = Rho * Uo (Y - boundary Layer Thickness) + Integral between limits y = boundary layer thickness to y = 0 for Rho * u dy

Rho = Density , Uo = free stream velocity u = air velocity
Rob Wilk is offline   Reply With Quote

Old   May 18, 2020, 12:17
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rob Wilk View Post
I did a Civil Engineering course some years ago and from my textbook I had a question on Turbulent boundary Flow.

In case you are wondering where my calculations are, I am actually blocked on this question and not sure where to start.

May someone please explain to me what things I need to understand, so that I can work towards the answer and get a better understanding on this ?

https://firebasestorage.googleapis.c...d-f2332b19c09c

[img]Question 9.jpg[/Iimg]

Take air density = 1.18 kg/m^3 and viscosity = 1.3 * 10^-5 Pa.s

Book Answer 1.31 kg/s


They say calculate mass flow rate entering air intake under zero pressure gradient conditions, how does zero pressure gradient come into this?

From my textbook it says
mass flow rate m = Rho * Uo (Y - boundary Layer Thickness) + Integral between limits y = boundary layer thickness to y = 0 for Rho * u dy

Rho = Density , Uo = free stream velocity u = air velocity



Zero pressure gradient is typical of flat plate flow, as illustrated in the classical Blasius solution for laminar flow. Just think of the external inviscid solution where dp/dx=0 from Bernoulli since u is constant.
FMDenaro is offline   Reply With Quote

Old   May 18, 2020, 12:44
Default Flow Rate from Turbulent Boundary Conditions
  #3
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Zero pressure gradient is typical of flat plate flow, as illustrated in the classical Blasius solution for laminar flow. Just think of the external inviscid solution where dp/dx=0 from Bernoulli since u is constant.

That equation I mentioned from my textbook is that what is used here for this question

mass flow rate m = Rho * Uo (Y - boundary Layer Thickness) + Integral between limits y = boundary layer thickness to y = 0 for Rho * u dy
Rob Wilk is offline   Reply With Quote

Old   May 18, 2020, 12:57
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
If you prescribe a velocity profile u=u(y), the mass flow rate across an unlimited vertical section is



rho0* Int[0,+Inf] u(y) dy


but this integral clearly diverges. For that reason one works in terms of mass defect.

If you consider a function delta(x) as sup extrema of the integral it depends on the arbitrary definition of the boundary layer function.
FMDenaro is offline   Reply With Quote

Old   May 19, 2020, 10:42
Default Flow Rate from Turbulent Boundary Conditions
  #5
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
If you prescribe a velocity profile u=u(y), the mass flow rate across an unlimited vertical section is



rho0* Int[0,+Inf] u(y) dy


but this integral clearly diverges. For that reason one works in terms of mass defect.

If you consider a function delta(x) as sup extrema of the integral it depends on the arbitrary definition of the boundary layer function.
Is the boundary layer profile given in the image for my question used for calculating boundary layer thickness ?
This is u/uo = (y / boundary layer thickness)^1/7
and Shear Stress Tau w = 0.025 and the rest of that formula.

It seems to me that you may have to calculate boundary layer thickness first and then the mass flow rate after that.

They do mention in the question that a turbulent boundary layer grows for a distance of 40 m, where does this come into it ?
Rob Wilk is offline   Reply With Quote

Old   May 19, 2020, 11:39
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rob Wilk View Post
Is the boundary layer profile given in the image for my question used for calculating boundary layer thickness ?
This is u/uo = (y / boundary layer thickness)^1/7
and Shear Stress Tau w = 0.025 and the rest of that formula.

It seems to me that you may have to calculate boundary layer thickness first and then the mass flow rate after that.

They do mention in the question that a turbulent boundary layer grows for a distance of 40 m, where does this come into it ?



The prescribed power-law is one approximation of the mean velocity profile in the turbulent BL. I strongly suggest to read Sec.7.3 in the textbook of Pope. You will see that the definition of delta(x), as the boundary at which the 99% of the mean velocity is obtained, is arbitrary and different lenghts can be deduced (7.126-127).

Your power-law velocity profile is depicted in Exercise 7.22.
FMDenaro is offline   Reply With Quote

Old   May 19, 2020, 11:42
Default Flow Rate from Turbulent Boundary Conditions
  #7
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The prescribed power-law is one approximation of the mean velocity profile in the turbulent BL. I strongly suggest to read Sec.7.3 in the textbook of Pope. You will see that the definition of delta(x), as the boundary at which the 99% of the mean velocity is obtained, is arbitrary and different lenghts can be deduced (7.126-127).

Your power-law velocity profile is depicted in Exercise 7.22.
What is the textbook actually called, I might be able to find it online maybe through google books ?
Rob Wilk is offline   Reply With Quote

Old   May 19, 2020, 11:50
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rob Wilk View Post
What is the textbook actually called, I might be able to find it online maybe through google books ?

https://www.cambridge.org/core/books...FAC9ED16486B3A
FMDenaro is offline   Reply With Quote

Old   May 20, 2020, 11:55
Default Flow Rate from Turbulent Boundary Conditions
  #9
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Hi Denaro

I've had a look at that Turbulent Flows book with the section 7 you suggested
on page 299 and 313 in the book.

I'm still confused on how this can be used to work towards the answer.

Are you able to please show me the process of how these calculations are meant to work to work towards the answer for this question ?
Rob Wilk is offline   Reply With Quote

Old   May 20, 2020, 11:59
Default Flow Rate from Turbulent Boundary Conditions
  #10
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Hi Denaro

I've had a look at that Turbulent Flows book with the section 7 you suggested
on page 299 and 313 in the book.

I'm still confused on how this can be used to work towards the answer.

Are you able to please show or explain to me, the process of how these calculations are meant to work to get to the answer for this question ?
Rob Wilk is offline   Reply With Quote

Old   May 20, 2020, 17:04
Default
  #11
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rob Wilk View Post
Hi Denaro

I've had a look at that Turbulent Flows book with the section 7 you suggested
on page 299 and 313 in the book.

I'm still confused on how this can be used to work towards the answer.

Are you able to please show or explain to me, the process of how these calculations are meant to work to get to the answer for this question ?



Use the law delta(x) for turbulent BL you can find in many textbooks (for example 0.38*x/Re_x^1/5 and integrate u(y)/U0 from zero to delta(x).
FMDenaro is offline   Reply With Quote

Old   May 23, 2020, 11:23
Default Flow Rate from Turbulent Boundary Conditions
  #12
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Use the law delta(x) for turbulent BL you can find in many textbooks (for example 0.38*x/Re_x^1/5 and integrate u(y)/U0 from zero to delta(x).
What is law delta(x) ?
Rob Wilk is offline   Reply With Quote

Old   May 23, 2020, 11:26
Default
  #13
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
The approximate expression for the boundary layer thickness along x
FMDenaro is offline   Reply With Quote

Old   May 27, 2020, 14:08
Default Flow Rate from Turbulent Boundary Conditions
  #14
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Use the law delta(x) for turbulent BL you can find in many textbooks (for example 0.38*x/Re_x^1/5 and integrate u(y)/U0 from zero to delta(x).
I found this Boundary Layer 01 youtube video https://www.youtube.com/watch?v=35qN3jjyIL0.
5 minute in it talks about mass flow rate.

It looks like it does the integration u(y)/U0 from zero to delta(x)

From the velocity profile and shear stress distribution I get boundary layer thickness= 0.37*x/Re_x^1/5 for this question.

x maximum for this question is 40 m isn't it.
in the video they talk about b, that seems like it is the same as x
Rob Wilk is offline   Reply With Quote

Old   May 27, 2020, 14:15
Default
  #15
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rob Wilk View Post
I found this Boundary Layer 01 youtube video https://www.youtube.com/watch?v=35qN3jjyIL0.
5 minute in it talks about mass flow rate.

It looks like it does the integration u(y)/U0 from zero to delta(x)

From the velocity profile and shear stress distribution I get boundary layer thickness= 0.37*x/Re_x^1/5 for this question.

x maximum for this question is 40 m isn't it.
in the video they talk about b, that seems like it is the same as x



For each position x you compute the thickness delta(x), then the flow rate through the section y=0 .. delta(x). Change the position to a different x and repeat. Of course you have a flow rate also across delta from two different positions.

That's all.
FMDenaro is offline   Reply With Quote

Old   May 28, 2020, 05:02
Default Flow Rate from Turbulent Boundary Conditions
  #16
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
For each position x you compute the thickness delta(x), then the flow rate through the section y=0 .. delta(x). Change the position to a different x and repeat. Of course you have a flow rate also across delta from two different positions.

That's all.

I think I am on the right track with getting closer to solving this, just want to check a few things.

When you say for each position x, do you mean first for position x = 0 at the start of the air intake and for position x = 40 m (end of boundary layer)

I know that flow rate = area * velocity and we are given 300mm * 300mm
and for mass flow rate you have to take into account density as well.

From velocity profile u / uo = (y / delta(x))^1/7

We can say that the velocity u = uo * (y / delta(x))^1/7

so from this equation mass flow rate I think needs to be integrated for y between limits of 0 and maximum boundary layer thickness as taken from the video I showed you.

so mass flow rate = Integral between limits of 0 to boundary layer thickness for uo * (y / delta(x))^1/7 * area * dy

and y gets integrated between the limits I mentioned

Is this the final equation ?

Last edited by Rob Wilk; May 28, 2020 at 06:12.
Rob Wilk is offline   Reply With Quote

Old   May 28, 2020, 05:17
Default
  #17
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Yes, at each position between x=0 and x=end you compute the BL thickness and integrate in y. The assumption is for 2d flow so that the third dimension is simple assumed to be unity in lenght (1 meter dimensionally). However for the mass flow rate you need to multiply for the constant density otherwise it is the volumetric flow rate.
FMDenaro is offline   Reply With Quote

Old   May 28, 2020, 09:09
Default Flow Rate from Turbulent Boundary Conditions
  #18
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Yes, at each position between x=0 and x=end you compute the BL thickness and integrate in y. The assumption is for 2d flow so that the third dimension is simple assumed to be unity in lenght (1 meter dimensionally). However for the mass flow rate you need to multiply for the constant density otherwise it is the volumetric flow rate.
It looks like i am on the right track then with the
equation mass flow rate = Integral between limits of 0 to boundary layer thickness for Rho * uo * (y / delta(x))^1/7 * area * dy

I added Rho because you have to take density into account as well.

In the video https://www.youtube.com/watch?v=35qN3jjyIL0 5 minutes 55 seconds in, they have used elemental area of A = b * dy.
For this question though, they say that the intake is square and 300 mm * 300 mm, so how do I deal with this ?

When I get the chance I'll go through the calculations and see if I can get the answer of 1.31 kg/s
Rob Wilk is offline   Reply With Quote

Old   May 29, 2020, 10:02
Default Flow Rate from Turbulent Boundary Conditions
  #19
Member
 
Rob Wilk's Avatar
 
Rob Wilkinson
Join Date: Apr 2020
Location: Wellington, New Zealand
Posts: 69
Rep Power: 5
Rob Wilk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Yes, at each position between x=0 and x=end you compute the BL thickness and integrate in y. The assumption is for 2d flow so that the third dimension is simple assumed to be unity in lenght (1 meter dimensionally). However for the mass flow rate you need to multiply for the constant density otherwise it is the volumetric flow rate.
I'm nearly there with this, I think.
I have one important part to clarify.

In the Boundary Layer 01 video https://www.youtube.com/watch?v=35qN3jjyIL0 5 minutes 55 seconds in, they have used elemental area of A = b * dy.

For this question though, they have given both the width and height of 300 mm * 300 mm for the square air intake.

So how do I deal with the elemental area for this question or do I not need to, can I take the 300m * 300mm as the area ?
Rob Wilk is offline   Reply With Quote

Old   May 29, 2020, 10:35
Default
  #20
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rob Wilk View Post
I'm nearly there with this, I think.
I have one important part to clarify.

In the Boundary Layer 01 video https://www.youtube.com/watch?v=35qN3jjyIL0 5 minutes 55 seconds in, they have used elemental area of A = b * dy.

For this question though, they have given both the width and height of 300 mm * 300 mm for the square air intake.

So how do I deal with the elemental area for this question or do I not need to, can I take the 300m * 300mm as the area ?
You need of a vertical surface that is delta(x)*(spanwise length). In the 2D theory, the spanwise lenght is assumed to be 1
FMDenaro is offline   Reply With Quote

Reply


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
Centrifugal fan j0hnny CFX 13 October 1, 2019 13:55
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 06:20
Difficulty In Setting Boundary Conditions Moinul Haque CFX 4 November 25, 2014 17:30
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
[Help] Boundary conditions for turbulent flow (kEpsilon) schaumy OpenFOAM 1 May 24, 2011 09:48


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