CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Implementation of trip in Boundary Layer flow

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

Like Tree1Likes
  • 1 Post By mAlletto

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 29, 2018, 18:01
Exclamation Implementation of trip in Boundary Layer flow
  #1
New Member
 
Pedro
Join Date: Feb 2016
Location: United States
Posts: 23
Rep Power: 8
pedramtx is on a distinguished road
Hi All,
I am trying to implement a numerical tripping by applying a Gaussian force in the wall-normal direction in OpenFOAM. The reason for this implementation is that I'd like to force the transition to turbulence in flat plate boundary layer flow. This idea can be done by adding a term to the momentum equation of the solver. The volume force has been defined in Schlatter and Orlu (2012) "Schlatter, P., & Örlü, R. (2012). Turbulent boundary layers at moderate Reynolds numbers: inflow length and tripping effects. Journal of Fluid Mechanics, 710, 5-34."
I am seeking help regarding implementation of this force defined in spectral domain (the ineverse Fourier) into the momentum equation in physical domain. Any help is greatly appreciated.

Thanks,
Pedram
Attached Images
File Type: jpg Tripping.jpg (156.5 KB, 60 views)

Last edited by pedramtx; August 30, 2018 at 19:43.
pedramtx is offline   Reply With Quote

Old   September 5, 2018, 13:41
Default
  #2
Senior Member
 
Michael Alletto
Join Date: Jun 2018
Location: Bremen
Posts: 603
Rep Power: 14
mAlletto will become famous soon enough
Just from the picture you showed it seems that the force is a function of x,y,z,t. The random numbers h are function of z and t. Can you provide some more detail?
mAlletto is offline   Reply With Quote

Old   September 5, 2018, 13:55
Default
  #3
New Member
 
Pedro
Join Date: Feb 2016
Location: United States
Posts: 23
Rep Power: 8
pedramtx is on a distinguished road
Quote:
Originally Posted by mAlletto View Post
Just from the picture you showed it seems that the force is a function of x,y,z,t. The random numbers h are function of z and t. Can you provide some more detail?
Hi Michael,

You are right. The random force g is a function of z and t. I want to apply a 1D wavenumber in the spanwise direction. Then the whole term F is expressed by a local Gaussian distribution in z plane.
pedramtx is offline   Reply With Quote

Old   September 6, 2018, 03:01
Default
  #4
Senior Member
 
Michael Alletto
Join Date: Jun 2018
Location: Bremen
Posts: 603
Rep Power: 14
mAlletto will become famous soon enough
Maybe I did not und understand correctly the question.
mAlletto is offline   Reply With Quote

Old   September 13, 2018, 04:18
Default
  #5
Senior Member
 
Michael Alletto
Join Date: Jun 2018
Location: Bremen
Posts: 603
Rep Power: 14
mAlletto will become famous soon enough
I found by accident this post of Schlatter where he explains how he generated the random harmonic signal: https://lists.mcs.anl.gov/pipermail/nek5000-users/2018-February/005087.html

in order to generate the random source term in order to simulate the boundary layer tipping.

pedramtx likes this.
mAlletto is offline   Reply With Quote

Old   September 13, 2018, 13:42
Default
  #6
New Member
 
Pedro
Join Date: Feb 2016
Location: United States
Posts: 23
Rep Power: 8
pedramtx is on a distinguished road
Quote:
Originally Posted by mAlletto View Post
I found by accident this post of Schlatter where he explains how he generated the random harmonic signal: https://lists.mcs.anl.gov/pipermail/nek5000-users/2018-February/005087.html

in order to generate the random source term in order to simulate the boundary layer tipping.
Thank you, Michael. This is so helpful.
pedramtx is offline   Reply With Quote

Old   September 14, 2018, 12:08
Default
  #7
Senior Member
 
Michael Alletto
Join Date: Jun 2018
Location: Bremen
Posts: 603
Rep Power: 14
mAlletto will become famous soon enough
You're welcome

Let me know how it goes. I find the topic interesting
mAlletto is offline   Reply With Quote

Old   September 18, 2018, 04:37
Default
  #8
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 307
Rep Power: 13
agustinvo is on a distinguished road
Very interesting topic! I will follow it. It is also interesting for me. For the moment, what I do is to use swak4Foam and introduce some disturbances in the wall.
agustinvo is offline   Reply With Quote

Old   September 18, 2018, 23:26
Default
  #9
New Member
 
Pedro
Join Date: Feb 2016
Location: United States
Posts: 23
Rep Power: 8
pedramtx is on a distinguished road
Quote:
Originally Posted by agustinvo View Post
Very interesting topic! I will follow it. It is also interesting for me. For the moment, what I do is to use swak4Foam and introduce some disturbances in the wall.
Agustín,
Thank you for your suggestion. Still trying to make it work.
pedramtx is offline   Reply With Quote

Old   September 20, 2018, 04:06
Default
  #10
Senior Member
 
Ehsan Asgari
Join Date: Apr 2010
Posts: 471
Rep Power: 16
syavash is on a distinguished road
Quote:
Originally Posted by agustinvo View Post
Very interesting topic! I will follow it. It is also interesting for me. For the moment, what I do is to use swak4Foam and introduce some disturbances in the wall.



Why not use a more sophisticated method like LeMOS or SEM at the inlet? Or a step further ahead, use the method introduced by Lund et al. (JCP 1998) in your simulation.


Regards,
Syavash
syavash is offline   Reply With Quote

Old   September 20, 2018, 04:14
Default
  #11
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 307
Rep Power: 13
agustinvo is on a distinguished road
Quote:
Originally Posted by syavash View Post
Why not use a more sophisticated method like LeMOS or SEM at the inlet? Or a step further ahead, use the method introduced by Lund et al. (JCP 1998) in your simulation.


Regards,
Syavash

Hi Syavash,
my case is a natural convection boundary layer, not a forced convection one. The boundary layer starts developing when the wall is heated. If I introduce an inlet flow, I'd have a mixed convection case.


I will check later the different approaches you said, I have some colleagues who may be interested.


Agustín
agustinvo is offline   Reply With Quote

Old   September 20, 2018, 04:20
Default
  #12
Senior Member
 
Ehsan Asgari
Join Date: Apr 2010
Posts: 471
Rep Power: 16
syavash is on a distinguished road
Quote:
Originally Posted by agustinvo View Post
Hi Syavash,
my case is a natural convection boundary layer, not a forced convection one. The boundary layer starts developing when the wall is heated. If I introduce an inlet flow, I'd have a mixed convection case.


I will check later the different approaches you said, I have some colleagues who may be interested.


Agustín



Agustin,


Having a naturally developing boundary-layer to turbulence should be very interesting, and challenging! Actually I got very interested in this problem. Could you please provide me a snapshot of your simulation? Also would you introduce me a paper dealing with the same problem using LES?


Regards,
Syavash
syavash is offline   Reply With Quote

Old   September 20, 2018, 11:07
Default
  #13
New Member
 
Pedro
Join Date: Feb 2016
Location: United States
Posts: 23
Rep Power: 8
pedramtx is on a distinguished road
Quote:
Originally Posted by syavash View Post
Why not use a more sophisticated method like LeMOS or SEM at the inlet? Or a step further ahead, use the method introduced by Lund et al. (JCP 1998) in your simulation.


Regards,
Syavash
Syavash,

Thank you for your suggestion. We already captured the turbulent BL using Lund method, but the current simulation aims at resolving the transition behavior of BL. As you may know, rescaling procedure in Lund method necessitates equilibrium state (developed turbulent field) at recycling station.


Thanks,
Pedram

Last edited by pedramtx; September 20, 2018 at 13:58.
pedramtx is offline   Reply With Quote

Old   September 15, 2022, 11:20
Default
  #14
New Member
 
Jin Zhang
Join Date: May 2018
Location: Germany
Posts: 15
Rep Power: 6
sjlouie91 is on a distinguished road
Quote:
Originally Posted by pedramtx View Post
Agustín,
Thank you for your suggestion. Still trying to make it work.
Hi Pedram,

I also find this numerical tripping quite interesting. Although it's been a long time since you last posted, I am still very curious about how this method is implemented in OpenFOAM. Have you already solved this problem? I would appreciate it if you could share your code or give me some clues on how to implement it. Thanks.

kind regards,
Jin
sjlouie91 is offline   Reply With Quote

Reply

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
sliding mesh problem in CFX Saima CFX 46 September 11, 2021 08:38
Wind turbine simulation Saturn CFX 58 July 3, 2020 02:13
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 08:30
CFD analaysis of Pelton turbine amodpanthee CFX 31 April 19, 2018 19:02
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 09:44


All times are GMT -4. The time now is 01:49.