https://linen.dev logo
Join Discord
Powered by
# les-ras
  • k

    kandelabr

    05/25/2023, 5:27 AM
    ...and the dreaded source code... https://www.openfoam.com/documentation/guides/latest/api/dir_033281ae6918016c4a70b6356445b423.html
  • t

    t_bo99

    05/25/2023, 11:30 PM
    Thank you very much for answering. It will be nice for future reference since I plan to use kOmegaSSTIDDES also, but I forgot to mention that its was for SpalartAlamarasIDDES....sorry my bad....but I figured out how to find it for SAIDDES from your link...thanks again! I guess I'll need to use the same approach in order to know how the different section (e.i. IDDESdelta, IDDESDeltaCoeffs, smoothCoeffs, etc.) work all together. And about the
    smoothCoeffs
    subdictionary, does someone has reference to suggest about the topic? Thank again for the help!
  • g

    godoftheuniverse

    05/27/2023, 7:11 AM
    has anyone used kw sst on rhoCentralFoam?
  • t

    t_bo99

    06/30/2023, 2:29 AM
    Im wondering if someone know, in LES, how y+ is calculated? Do we use nu, nuEff or nuSgs? Since its unsteady, is it preferable to get a mean value or y+min and y+max over every recorded time step? Thanks for the help.
  • s

    slopezcastano

    06/30/2023, 12:07 PM
    $$ y^+ = \frac{u_\tau y}{\nu} $$
  • t

    TeXit

    06/30/2023, 12:07 PM
    slopezcastano

    https://cdn.discordapp.com/attachments/598785402391887872/1124310272119685170/725986814787059794.png▾

  • s

    slopezcastano

    06/30/2023, 12:09 PM
    $$ u_\tau = \sqrt{\frac{\tau_w}{\rho}} $$
  • t

    TeXit

    06/30/2023, 12:10 PM
    slopezcastano

    https://cdn.discordapp.com/attachments/598785402391887872/1124311018584146030/725986814787059794.png▾

  • t

    t_bo99

    06/30/2023, 12:30 PM
    That's the function object code. So, nu, nuEff or nuSgs? y+mean or y+max for instance over all time step?

    https://cdn.discordapp.com/attachments/598785402391887872/1124315938251091988/image.png▾

  • s

    slopezcastano

    06/30/2023, 12:36 PM
    I understand LES as wall-resolved LES, so that function object is not correct for wall-resolved LES
  • s

    slopezcastano

    06/30/2023, 12:36 PM
    WMLES is another story
  • t

    t_bo99

    06/30/2023, 12:44 PM
    Thank for answering and for your help. I'm trying to implement the Schaumann BC (related to atmospheric boundary layer simulation) and some authors say that the BC perform better when y+>20-30, so I'm trying to find out how it must be calculated..
  • s

    slopezcastano

    06/30/2023, 12:51 PM
    they use the formula above, but for u_tau they calculate the following:
  • s

    slopezcastano

    06/30/2023, 12:52 PM
    $$ u_\tau = \sqrt{\nu_t \frac{\partial U}{\partial y}} $$
  • t

    TeXit

    06/30/2023, 12:52 PM
    slopezcastano

    https://cdn.discordapp.com/attachments/598785402391887872/1124321601345830972/725986814787059794.png▾

  • s

    slopezcastano

    06/30/2023, 12:53 PM
    IT should be nu, not nu_t, if LES
  • s

    slopezcastano

    06/30/2023, 12:54 PM
    So, if you want to guarantee such in LES you must set your nuSgs to be zero at the wall boundary
  • t

    t_bo99

    06/30/2023, 12:55 PM
    Ah.....I have to take a look....I maybe set this to zeroGradient....
  • t

    t_bo99

    06/30/2023, 12:57 PM
    But again, regarding y+, the Shaumann BC involve applying a shear stress at the wall and it is function of the roughness lenght (z0), which mean its basically using this no? $$ u_\tau = \sqrt{\frac{\tau_w}{\rho}} $$
  • t

    TeXit

    06/30/2023, 12:57 PM
    t\_bo99

    https://cdn.discordapp.com/attachments/598785402391887872/1124322869921456170/574993721712443415.png▾

  • s

    slopezcastano

    06/30/2023, 12:59 PM
    Careful, here you are talking about WMLES
  • s

    slopezcastano

    06/30/2023, 1:00 PM
    $$ u_\tau = \sqrt{\frac{\tau_w}{\rho}} = \sqrt{\nu \frac{\partial U}{\partial y}} $$
  • t

    TeXit

    06/30/2023, 1:00 PM
    slopezcastano

    https://cdn.discordapp.com/attachments/598785402391887872/1124323598497235036/725986814787059794.png▾

  • t

    t_bo99

    06/30/2023, 1:02 PM
    So, if I understand well, then using this should be approriate?
  • t

    t_bo99

    06/30/2023, 1:02 PM
    This is what it look like:
  • t

    t_bo99

    06/30/2023, 1:03 PM
    Copy code
    scalar uStar = (kappa * U) / Foam::log(z1 / z0);
    
    //     Info << "planarAverage" << endl;
            Rw[facei].xx() = 0.0;
            Rw[facei].xy() = 0.0;
            Rw[facei].xz() = -sqr(uStar[facei]) * (UParallel[facei].x() / max(UParallelMeanMag, 1.0E-5));
            Rw[facei].yy() = 0.0;
            Rw[facei].yz() = -sqr(uStar[facei]) * (UParallel[facei].y() / max(UParallelMeanMag, 1.0E-5));
            Rw[facei].zz() = 0.0;
  • q

    qr

    07/11/2023, 6:03 AM
    What is an appropriate window when calcualting UPrime2Mean for the purpose of characterizing LES sims?
  • u

    .noble6ix

    07/22/2023, 9:19 PM
    Hello! I am trying to implement the kOmegaSST turbulence model into my compressible simulations of flow around a wedge. My solver is rhoCentralFoam. I believe I have set up my initial conditions correctly (0/ k omega nut alphat p U T), however, I keep receiving a floating point error when running rhoCentralFoam, relating to the nutkWallFunction. Is there more I had to do besides adding the necessary 0/ files?
  • u

    .noble6ix

    07/22/2023, 9:19 PM
    looking for general guidance!
  • u

    ⵣAryazⵣ

    07/22/2023, 10:10 PM
    > I have set up my initial conditions correctly
1...4849505152Latest