[maddy11]

Hi all,

I'm an OpenFoam beginner. In the simulation of supersonic cavity, I'm trying to provide the boundary layer using Blasius 1/7th Power law at the walls ( 0 folder --> U)

So I have instead provided the following command in the Terminal

applyBoundaryLayer -ybl 0.0025

where 0.0025 is my boundary layer thickness. So, this calculates the velocity in the internalfield in the U file. Is my understanding correct? Could anyone kindly help please?

[maddy11]

Anyone? Kindly share information pls!

[RobertHB]

Originally Posted by maddy11

Anyone? Kindly share information pls!

I gave the source code a quick read, and if i understand it correctly you need a velocity field beforehand. Not sure if only applying it to a empty dict with an internalField will work.
And you can not use wall functions for your turbulence parameters.
applyBoundaryLayer.C @ github

[maddy11]

Hi Robert,

Thank you for your kind reply. I'm a beginner. I have placed my file for velocity field below ( from '0' folder ) . I write the above mentioned command after this. And values are initialized in the 'internal field' in the velocity file subsequently for all the cells.

Can you kindly have a look at it? I have also placed the turbulence parameter 'k'. Kindly correct my mistakes.


FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (516 0 0);

boundaryField
{
outlet
{
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform (516 0 0);
}
wall
{
type fixedValue;
value uniform (0 0 0);
}
top
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}
}


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


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

internalField uniform 10;

boundaryField
{
outlet
{
type zeroGradient;
}
inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.005;
value uniform 10;
}
wall
{
type kqRWallFunction;
value uniform 10;
}
top
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;
}
}

[RobertHB]

Originally Posted by maddy11

Hi Robert,

Thank you for your kind reply. I'm a beginner. I have placed my file for velocity field below ( from '0' folder ) . I write the above mentioned command after this. And values are initialized in the 'internal field' in the velocity file subsequently for all the cells.

O.k., than forget my comment about needing a vel. field beforehand.

Can you kindly have a look at it? I have also placed the turbulence parameter 'k'. Kindly correct my mistakes.

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


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

internalField   uniform 10;

boundaryField
{
    outlet
    {
        type            zeroGradient;
    }
    inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
    intensity       0.005;
        value           uniform 10;
    }
    wall
    {
        type            kqRWallFunction;
        value           uniform 10;
    }
    top
    {
        type            zeroGradient;
    }
    frontAndBackPlanes
    {
        type            empty;
    }
 }

From what i can figure from the source code, you are not supposed to use a wallfunction for k (and nut, epsilon, omega aswell). Might be a good idea to go with zeroGradient instead. But i'm guessing here.

[maddy11]

Okay Robert, I will change as per your suggestion. Thank you

[Malkocoglu]

Originally Posted by maddy11

Okay Robert, I will change as per your suggestion. Thank you

Hello @maddy11

I am facing with similar situation. I wonder that if it worked to convert wall function into zeroGradient or fixedValue?

Thanks in advance, stay healthy,
Malkocoglu