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

channel flow, too large reynolds stress <u'u'> near the wall

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 5, 2020, 03:12
Default channel flow, too large reynolds stress <u'u'> near the wall
  #1
Member
 
niewiemnic
Join Date: Jan 2018
Location: Niemcy
Posts: 80
Rep Power: 8
spalartallmaras is on a distinguished road
Hallo,

I simulate channel flow for Re_tau=395 for long time, but I have always really large <u'u'> near the wall. I simulate with pimpleFoam.
Inlet, Outlet, left and right are cyclic and walls are wall.

My setup is,

1. nu=2,53e-03 (so 1/nu=395)
2. bulk velocity is about 22-23 m/s (patchMeanVelocityForce on inlet), so wallShearStressMean and u_tau are also about 1. Density is supposed to be 1 kg/m^3.
3. turbulence model is WALE-Modell
4. because Re_bulk is low, I use createBoxTurb before the pimpleFoam.
5. Mesh is really coarse, 32*32*32, because I want to make the mesh fine later.
6. After the flow is turbulent enough, I switch on the fieldAverage1.

but I have always <u'u'> about 20-22 m^2/s^2 near the wall. But it should be about 7-8 m^2/s^2 from Reference.

What will be the problem?

Last edited by spalartallmaras; January 5, 2020 at 22:15.
spalartallmaras is offline   Reply With Quote

Old   January 6, 2020, 03:30
Default
  #2
Member
 
niewiemnic
Join Date: Jan 2018
Location: Niemcy
Posts: 80
Rep Power: 8
spalartallmaras is on a distinguished road
Nobody has idea?
spalartallmaras is offline   Reply With Quote

Old   January 6, 2020, 03:30
Default
  #3
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 16
Santiago is on a distinguished road
Quote:
Originally Posted by spalartallmaras View Post
Hallo,

I simulate channel flow for Re_tau=395 for long time, but I have always really large <u'u'> near the wall. I simulate with pimpleFoam.
Inlet, Outlet, left and right are cyclic and walls are wall.

My setup is,

1. nu=2,53e-03 (so 1/nu=395)
2. bulk velocity is about 22-23 m/s (patchMeanVelocityForce on inlet), so wallShearStressMean and u_tau are also about 1. Density is supposed to be 1 kg/m^3.
3. turbulence model is WALE-Modell
4. because Re_bulk is low, I use createBoxTurb before the pimpleFoam.
5. Mesh is really coarse, 32*32*32, because I want to make the mesh fine later.
6. After the flow is turbulent enough, I switch on the fieldAverage1.

but I have always <u'u'> about 20-22 m^2/s^2 near the wall. But it should be about 7-8 m^2/s^2 from Reference.

What will be the problem?
Which schemes are you using? time step?
Santiago is offline   Reply With Quote

Old   January 6, 2020, 03:31
Default
  #4
Member
 
niewiemnic
Join Date: Jan 2018
Location: Niemcy
Posts: 80
Rep Power: 8
spalartallmaras is on a distinguished road
Quote:
Originally Posted by Santiago View Post
Which schemes are you using? time step?
Thank you for your answer, I will post it in 5-6h.

Time step is maxCo = 0.99.
spalartallmaras is offline   Reply With Quote

Old   January 6, 2020, 06:17
Default fvSchemes
  #5
Member
 
niewiemnic
Join Date: Jan 2018
Location: Niemcy
Posts: 80
Rep Power: 8
spalartallmaras is on a distinguished road
Here I post my fvSchemes file,

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         backward;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss linear;
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,B)      Gauss limitedLinear 1;
    div(B)          Gauss linear;
    div(U)          Gauss linear;
    div(phi,nuTilda) Gauss limitedLinear 1;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}



// ************************************************************************* //
spalartallmaras is offline   Reply With Quote

Old   October 22, 2024, 17:30
Default
  #6
Member
 
Shravan
Join Date: Mar 2017
Posts: 75
Rep Power: 9
Severus is on a distinguished road
Hello,
This issue has been discussed and solved in another post.
I am adding the link below so that it could be useful for others with the same problem.
channel flow, too large reynolds stress <u'u'> near the wall

Thanks
Severus is offline   Reply With Quote

Reply


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
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 37 July 20, 2020 06:38
Wall shear stress in channel flow sreeyuth OpenFOAM Running, Solving & CFD 2 September 5, 2014 11:02
UDF for wall slipping HFLUENT Fluent UDF and Scheme Programming 0 April 27, 2011 13:03
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13
X-Y plot of Yplus in Fluent 5.3 Luo FLUENT 24 April 11, 2000 07:07


All times are GMT -4. The time now is 13:02.