CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Ustar and Atmospheric Boundary Layer (

vinz July 9, 2012 11:51

Ustar and Atmospheric Boundary Layer
Hi everybody,

I was trying to use the ABL condition present in the turbineSitting tutorial in other simulation.

However I am wondering how is fixed the value of Ustar in this tutorial?
From what I know: Ustar = sqrt(Tau / rho) with Tau being the wall shear stress. The problem is this Tau is normally computed as: Tau = mu*dU/dz and I can only know the du/dz once I have done the computation.

So I have the felling that I am turning around and missing something.

Could somebody tell me how to know this value before doing the computation or advise me on how to proceed?



vinz July 10, 2012 03:18

Ok, so I did some more testing and this value of U* does not seem to be used.

I tried to do the computation by commenting the line, and by changing U* of 0.82, 5 and even 10, and I get the same result at the end.

Then I looked in the code "atmBoundaryLayerInletVelocityFvPatchVectorField.C " and I saw the following lines:


    kappa_(dict.lookupOrDefault<scalar>("kappa", 0.41)),
    zGround_("zGround", dict, p.size())

so the Ustart imposed into the ABLconditions file is not taken into account?

and apparently it is recomputed afterwards:

    Ustar_ = kappa_*Uref_/(log((Href_  + z0_)/max(z0_ , 0.001)));

All times are GMT -4. The time now is 09:12.