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

Linear Stability Analysis of Blasius Boundary Layer

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 14, 2015, 16:11
Default Linear Stability Analysis of Blasius Boundary Layer
  #1
Member
 
Obad
Join Date: Sep 2013
Posts: 42
Rep Power: 6
Obad is on a distinguished road
Hi folks,

I'm currently trying to solve the Orr-Sommerfeld equation (OSE) for Blasius flow. I want to use an approach where the equation is discretizes with chebychev polinomials (chebychev collocation) and then calculate the eigenvalues. I do all of that with Matlab.

However, I'm already stuck with it for a couple of days...
I have a Matlab code to solve the OSE for Poiseulle flow, which works fine. But it don't know how to adapt the code to the Blasius case.

It seems to me that no one is using Chebychev Collocation to solve Blasius flow, at least Google doesn't tell me.


Does anyone have a clue how to handle Blasius flow?

Cheers!
Obad is offline   Reply With Quote

Old   June 15, 2015, 00:58
Default
  #2
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 189
Rep Power: 10
harishg is on a distinguished road
Look into www.channelflow.org.
harishg is offline   Reply With Quote

Old   July 21, 2015, 12:18
Default
  #3
Member
 
Jingchang.Shi
Join Date: Aug 2012
Location: Xi'an, China
Posts: 69
Rep Power: 7
aerosjc is on a distinguished road
I'm also stuck with blasius linear stability calculation!
I've already written f90 program to solve poiseuille flow successfully. so did my spectral collocation method. however, blasius case is wrong for both codes.
I suspect that the key is some algebra mapping on the coordinates between blasius similarly variable and collocation points.
I'm not successful for the moment and I hope to discuss with you!
sincerely
aerosjc is offline   Reply With Quote

Old   July 22, 2015, 18:06
Default
  #4
Member
 
Obad
Join Date: Sep 2013
Posts: 42
Rep Power: 6
Obad is on a distinguished road
Hi,
well I made it

If you worked it out for the Poiseuille flow then you are very close to solving the Blasius flow. The code is exactly the same.

Concerning the coordinate system you might have look at this master thesis: http://tuprints.ulb.tu-darmstadt.de/3173/ There is a chapter about the linear stability analysis of laminar boundary layers with a spectral collocation method using Chebychev polynomials.
Furthermore check out this website: http://www.lmm.jussieu.fr/~hoepffner/codes.php

The important thing about solving the stability of a boundary layer is, that you use a base velocity profile that was calculated across a couple of times the boundary layer thickness, so that a substantial amount of the velocity profile lies withing the free stream. This has to do with the farfield boundary condition that you impose on the upper edge where you say that the perturbation goes to zero as the distance from the wall goes towards large values.

If you are using a temporal theory approach then try to solve the equations for a range of alpha = 0.1 - 0.5.

Cheers!
Obad is offline   Reply With Quote

Old   July 23, 2015, 11:19
Default
  #5
Member
 
Jingchang.Shi
Join Date: Aug 2012
Location: Xi'an, China
Posts: 69
Rep Power: 7
aerosjc is on a distinguished road
Quote:
Originally Posted by Obad View Post
Hi,
well I made it

If you worked it out for the Poiseuille flow then you are very close to solving the Blasius flow. The code is exactly the same.

Concerning the coordinate system you might have look at this master thesis: http://tuprints.ulb.tu-darmstadt.de/3173/ There is a chapter about the linear stability analysis of laminar boundary layers with a spectral collocation method using Chebychev polynomials.
Furthermore check out this website: http://www.lmm.jussieu.fr/~hoepffner/codes.php

The important thing about solving the stability of a boundary layer is, that you use a base velocity profile that was calculated across a couple of times the boundary layer thickness, so that a substantial amount of the velocity profile lies withing the free stream. This has to do with the farfield boundary condition that you impose on the upper edge where you say that the perturbation goes to zero as the distance from the wall goes towards large values.

If you are using a temporal theory approach then try to solve the equations for a range of alpha = 0.1 - 0.5.

Cheers!
Thanks for your help!
I will read the paper in the link. I hope I can also make it.
And I have another question. Do you have any experience on shooting method to solve OS equation? I mean the compound matrix method. I've no idea why my shooting method does not work on Blasius case though it works on Poiseuille case. I want to show you my codes implementing shooting method and hope you could spare some time finding why it does not work. However, it's not an appropriate request. And you can ignore it.
aerosjc is offline   Reply With Quote

Old   July 25, 2015, 18:39
Default
  #6
Member
 
Obad
Join Date: Sep 2013
Posts: 42
Rep Power: 6
Obad is on a distinguished road
Hey,

yeah at first I also wanted to solve the spatial theory by using the CMM. I tried it out, but left it unfinished since the Spectral method worked.

Concerning your code I'm sorry, but I'm not familiar with Fortran coding.
Obad is offline   Reply With Quote

Old   December 2, 2017, 16:18
Default
  #7
Member
 
sajad
Join Date: Aug 2013
Posts: 65
Rep Power: 6
seju is on a distinguished road
Quote:
Originally Posted by Obad View Post
Hi,
well I made it

If you worked it out for the Poiseuille flow then you are very close to solving the Blasius flow. The code is exactly the same.

Concerning the coordinate system you might have look at this master thesis: http://tuprints.ulb.tu-darmstadt.de/3173/ There is a chapter about the linear stability analysis of laminar boundary layers with a spectral collocation method using Chebychev polynomials.
Furthermore check out this website: http://www.lmm.jussieu.fr/~hoepffner/codes.php

The important thing about solving the stability of a boundary layer is, that you use a base velocity profile that was calculated across a couple of times the boundary layer thickness, so that a substantial amount of the velocity profile lies withing the free stream. This has to do with the farfield boundary condition that you impose on the upper edge where you say that the perturbation goes to zero as the distance from the wall goes towards large values.

If you are using a temporal theory approach then try to solve the equations for a range of alpha = 0.1 - 0.5.

Cheers!
Hi, Good For you. I wanna do stability analysis for boundary layer. I confuse For doing that. I just can drive Orr Sommerfeld but I dont know How I should solve It. can You help me please?
best regards
Sajad
seju 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
sliding mesh problem in CFX Saima CFX 45 September 22, 2015 10:53
Divide Prismatic Boundary Layer Mesh causes overlapping faces SilentRunner42 enGrid 4 May 4, 2015 03:37
Radiation interface hinca CFX 15 January 26, 2014 18:11
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
Boundary layer in a pipe Clementhuon OpenFOAM Native Meshers: snappyHexMesh and Others 6 March 12, 2012 13:41


All times are GMT -4. The time now is 22:39.