[dickcruz]

Hello All,

I've been trying to use the y+ calculation option on a transient external flow case using pimpleFoam -postProcess to make sure that my mesh resolution is sufficient closer to the flow surfaces.

Unfortunately, I get zeros for all my surfaces as my min, max and average. I've found that if I run an LES on the same case, I'm able to get some numbers. Similarly, if after I run the RAS simulation, I change up my turbulence model to LES before I run the yPlus utility, I'm able to get non-zero numbers for the RAS. I'm uncertain how trustworthy this number is. But, more fundamentally, I don't understand what I'm doing wrong with the model.

Please let me know what data needs to be analyzed to figure this out.

[vdeschuy]

Hey,

I think I'm dealing with a similar issue!
Are you also using OF version 5?
Do you mean with 'zeros' a very low number or actually zero?

I'm having the problem where even for very coarse meshes my y+ is always smaller than 1. As in your case, when I change to LES before I postprocess yPlus, I get numbers above 1, but I'm also not sure these are trustworthy.

Maybe you or someone else has found a solution already?

[dickcruz]

Originally Posted by vdeschuy

Hey,

I think I'm dealing with a similar issue!
Are you also using OF version 5?
Do you mean with 'zeros' a very low number or actually zero?

I'm having the problem where even for very coarse meshes my y+ is always smaller than 1. As in your case, when I change to LES before I postprocess yPlus, I get numbers above 1, but I'm also not sure these are trustworthy.

Maybe you or someone else has found a solution already?

Hello, could you paste the contents of your nut file here, please?

[vdeschuy]

Here you go!

This is for a case were I used a very coarse mesh, wall functions, and here I expected to see a high y+ value.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            calculated;
        value           uniform 0;
    }
    outlet
    {
        type            calculated;
        value           uniform 0;
    }
    top
    {
        type            slip;
    }
    bottom
    {
        type            slip;
    }
    cylinder
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
    side1
    {
        type            empty;
    }
    side2
    {
        type            empty;
    }
}

// ************************************************************************* //

[cnsidero]

The nutkWallFunction requires k to compute y+. Checkout out the source code for that BC:

https://github.com/OpenFOAM/OpenFOAM...arField.C#L157

The SA turbulence model has no ability to compute k so it returns a value of zero for k. Checkout the source code:

https://github.com/OpenFOAM/OpenFOAM...llmaras.C#L302

Instead use the nutUWallFunction which uses a velocity instead of k to compute y+.

-Chris