CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM News & Announcements > OpenFOAM Announcements from Other Sources

New library for wall-modelled LES with OpenFOAM

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

Like Tree9Likes
  • 8 Post By tiam
  • 1 Post By amanbearpig

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 17, 2018, 13:04
Default New library for wall-modelled LES with OpenFOAM
  #1
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Gothenburg, Sweden
Posts: 114
Rep Power: 10
tiam is on a distinguished road
Dear all,


For some years now, me and my colleagues have been working on implementing a library of wall models for LES, expanding the native capabilities of OpenFOAM. In the latest release, v0.4.0, we have added support for the latest versions of OpenFOAM, both Foundation and OpenCFD releases. It is thus a good time-point to (re-)introduce the code to the community.


The library can be found here

https://bitbucket.org/lesituu/libwallmodelledles

The README file provides quite a bit of information about what the library does, I copy the Key features section below
  • Provides a number of wall models, based on both non-linear algebraic and ordinary differential equations, see the class-headers in the wallModels folder.
  • Makes it possible to specify the distance to the wall model's sampling point, h, on a per-face basis.
  • Allows the user to control all the other parameters of wall modelling, e.g. model constants, iterative solver settings etc.
  • Serves a as a convenient framework for implementing new models without a lot of code duplication.
An in-depth discussion along with simulation results for turbulent channel flow and flow over a BFS can be found here.
https://arxiv.org/abs/1807.11786
More works using the library can be found in the end of the README.


For obvious reasons, it has not be possible to thoroughly validate the code using all the supported OF versions. Therefore any issues with compilations and running are very welcome to be reported. Please also feel free to contact me regarding any possible collaborations, WMLES is becoming a hot topic, and joining forces is usually of benefit.
tiam is offline   Reply With Quote

Old   February 8, 2019, 08:26
Default
  #2
Member
 
Join Date: Apr 2011
Posts: 55
Rep Power: 11
amanbearpig is on a distinguished road
Hi Timofey,

Really great work with this! I've been playing around with it and have been pleased with my early testing. Thanks so much for making this publically available to the community. Is there a foam-extend version available or in the works?
tiam likes this.
amanbearpig is offline   Reply With Quote

Old   February 8, 2019, 09:22
Default
  #3
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Gothenburg, Sweden
Posts: 114
Rep Power: 10
tiam is on a distinguished road
Quote:
Originally Posted by amanbearpig View Post
Hi Timofey,

Really great work with this! I've been playing around with it and have been pleased with my early testing. Thanks so much for making this publically available to the community. Is there a foam-extend version available or in the works?

Hi!
Great to hear that it's working out for you. foam-extend support is not really in the making at the moment. Unfortunately, I have no experience with that fork, but I can imagine that for someone who does it would be pretty easy to make the adjustments! In any case, the fact that you request it is a push towards me considering to spend time on it . It won't be quick though, since I have other developments in the pipeline, which have higher priority.
tiam is offline   Reply With Quote

Old   April 11, 2019, 05:00
Default
  #4
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Gothenburg, Sweden
Posts: 114
Rep Power: 10
tiam is on a distinguished road
The published version of the article describing the library can be downloaded for free until may the 30th by following this link


https://authors.elsevier.com/a/1Ysl32OInZsxl
tiam is offline   Reply With Quote

Old   August 1, 2019, 09:10
Default
  #5
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Gothenburg, Sweden
Posts: 114
Rep Power: 10
tiam is on a distinguished road
I've just tagged v0.5.0 on Bitbucket. Some highlights


- Unit and integration tests.
- Stability improvements
- Code refactoring for better testebility and extensibility.
- Multicell sampling is included as a "beta version", currently not used by any wall models.
tiam is offline   Reply With Quote

Old   May 22, 2020, 01:03
Default
  #6
New Member
 
Tommy
Join Date: May 2020
Location: United States
Posts: 2
Rep Power: 0
Tj_m3 is on a distinguished road
Hi Timofey,

Many thanks for making this available. I am running into an issue with the parallel decomposition of the h value on the wall. The case is an airfoil in freestream flow. After the mesh decomposition, the processor0 folder shows an h value of "0" on the airfoil boundary, although it is set to the 2nd off the wall cell distance in the 0 folder. This seems to be triggering this warning in the log file:

"SingleCellSampler: max(h) is -1e+300 but no cell centres within distance 2*max(h) were found. Will sample from wall-adjacent cells"

h is being correctly set in the other processor folders - the ones which the airfoil wall resides on. But I am not sure those values are being used in the calculation.

Thanks!
Tommy
Tj_m3 is offline   Reply With Quote

Old   May 22, 2020, 13:54
Default
  #7
Senior Member
 
Timofey Mukha
Join Date: Mar 2012
Location: Gothenburg, Sweden
Posts: 114
Rep Power: 10
tiam is on a distinguished road
Quote:
Originally Posted by Tj_m3 View Post
Hi Timofey,

Many thanks for making this available. I am running into an issue with the parallel decomposition of the h value on the wall. The case is an airfoil in freestream flow. After the mesh decomposition, the processor0 folder shows an h value of "0" on the airfoil boundary, although it is set to the 2nd off the wall cell distance in the 0 folder. This seems to be triggering this warning in the log file:

"SingleCellSampler: max(h) is -1e+300 but no cell centres within distance 2*max(h) were found. Will sample from wall-adjacent cells"

h is being correctly set in the other processor folders - the ones which the airfoil wall resides on. But I am not sure those values are being used in the calculation.

Thanks!
Tommy

Hi Tommy,


Glad to see you pick up the library. I think you might be the first to use it for external flows, so I'm excited to see how it goes. Do I understand correctly that the wall boundary has no faces on processor0? In that case you don't have to worry about it. In any case, I would do the following. Open the case in Paraview, and under Case Type, select Decomposed Case. This way you actually open the parallelized data. If you load in the data for your wall patch you can directly inspect your h values.
To be 100% sure what cells are used for sampling, you can run the simulation for 1 time-step and save the data. There will be a samplingCells field, which marks the cells that are used for sampling. The value that they get should be equal to the index of the patch, whereas for the rest of the cells the value is set to -1.



It is highly possible that you can just ignore the warning about the h values, I am pretty sure it sometimes prints that for no reason, but I hadn't had the time to fix that---to busy with other projects.


Best,
Timofey
tiam is offline   Reply With Quote

Old   May 27, 2020, 14:44
Default
  #8
New Member
 
Tommy
Join Date: May 2020
Location: United States
Posts: 2
Rep Power: 0
Tj_m3 is on a distinguished road
Hi Timofey

Yes, there are no wall boundary faces on processor0, but when inspecting the decomposed "h" files and the samplingCells as you said I can see that the sampling cells are correct.

Thanks for clearing that up! I just wanted to be extra sure before starting the grid study with a fixed h value.

-Tommy
Tj_m3 is offline   Reply With Quote

Reply

Tags
les, wall modelling
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
Setting up LES of Airfoil in OpenFoam Jack001 OpenFOAM 4 April 12, 2020 03:09
LES models in OpenFOAM The King OpenFOAM Running, Solving & CFD 0 April 1, 2017 10:19
OpenFOAM Training, London, Chicago, Munich, Sep-Oct 2015 cfd.direct OpenFOAM Announcements from Other Sources 2 August 31, 2015 13:36
Serious bug in LES interface fs82 OpenFOAM Bugs 21 November 16, 2009 08:15


All times are GMT -4. The time now is 07:21.