[Daniele111]

Hi
I try to simulate an atmospheric boundary layer, on bottom I use wall function, but if i use nutRoughWallFunction, epsilon has a overshooting in last cells why?

Best Regards

[farbfilm]

Hi,

the overshooting is a known problem! Look for example here:
http://www.openfoamworkshop.org/2009...ner_slides.pdf

What Version of OpenFOAM are you using? I think, there have been some problems with the nutRoughWall-Condition in OpenFOAM 1.6! So, you should use OpenFOAM 1.6.x!!

Hope, that helps!

[Daniele111]

Hi
I read this paper, useful. But the problem is openfoam 1.6! I use now wall function of 1.6.x and I haven't overshooting problem,thanks you

[Daniele111]

using openfoam 1.6.x the overshooting there isn't; but I Have another problem: in atmospheric boundary layer simulation, th wall value of epsilon and k would be costant? In my simulation these value are constant only after some cells along x axis Why?

[Daniele111]

Hi
What are the correct boundary conditions for p field in these simulation?

Thanks

[farbfilm]

What is x?? Is it the height of your domain or the main flow direction??

[Daniele111]

Main flow direction, change inlet p condition:

[farbfilm]

Sorry, I don't understand what Wss means...
I've chosen a "zeroGradient"-condition for my ABL-Simulations at the inlet for p!

To achieve a homogenous ABL (constant values in flow direction), you have to calculate the right Ks (and Cs) in subject to your z0 (conversion of z0 to Ks due to your roughWallBC)!

Hope, that helps...

[Daniele111]

Wss, is wall shear stress, my z0 is 0.001, Ks=0.02, Cs=0.5 it's wrong?

Thanks you!

[Daniele111]

My Bc. Thanks you again

[farbfilm]

I'm getting a value of 0.019586! Almost the same...

I've used the formula:
Ks = E*z0/Cs where E=9.793 is a constant of the wall function!

[Daniele111]

My Ks=20*y0, but where is th mistake? Your Cs=0.5? E in your simulation epsilon, k are constant in the first cell along main flow direction, my x-axis?

[farbfilm]

I think, the velocity at the bottom should be zero (no slip condition)!

Unfortunately, I don't have any experiences with the timeVaryingMappedFixedValue-BC! So, I can't help with this! I used the GroovyBC for the analytical profiles!

Maybe, you should define the constants of the wall function in the nut-file! Just to be sure, you use the same like in your conversion from z0 to Ks! Actually, Cs is 0.327 in your BC's and not 0.5 like you mentioned before!

[Daniele111]

I just use groovyBc, timeVaryingMappedFixedValue-BC is not the problem, it's the same thing. I don't understand when you say that I would define the constants of the wall function in the nut-file. What do i do this? The velocity at the bottom has a value becouse I use like reference a paper that imposed it, but I pose it 0 and the problem there is. The p bc is the same that you have?

[farbfilm]

Yes, p BC is the same!!

[Daniele111]

Hi
Bah I try to use absurd values of Ks and Cs, but the result is always the same; I copied the turbulenceModels folder of OpenFOAM 1.6.x and I copy it in src folder of OpenFOAM 1.6 and the I compile it. It's wrong? Because if I change ks and cs not change the result?

[Daniele111]

Can i see your bc?

[farbfilm]

Hi,

sure! Here are my BC's! I definitely get a change when I use different values for Ks (see the picture)!

Hope that helps! Why do you want to simulate an ABL? What is the topic of your project??

[Daniele111]

Thank you! I must simulate sand dune evolution. So I would calculate shear stress at bottom. Thank you again
What are these fields in your groovy bc?

refValue uniform 0;
refGradient uniform 0;
valueFraction uniform 1;
value uniform 0.1;
valueExpression "para";
gradientExpression "0";
fractionExpression "1";

[Daniele111]

I try to use your bc setup, but the velocity field along main flow is these. It's correct? What is your first cell higth? My domain is higth 100m