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

Ustar and Atmospheric Boundary Layer

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 9, 2012, 11:51
Default Ustar and Atmospheric Boundary Layer
  #1
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 277
Rep Power: 9
vinz is on a distinguished road
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?

Thanks

Vincent
vinz is offline   Reply With Quote

Old   July 10, 2012, 03:18
Default
  #2
Senior Member
 
Vincent RIVOLA
Join Date: Mar 2009
Location: France
Posts: 277
Rep Power: 9
vinz is on a distinguished road
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:

Code:
    Ustar_(0),
    n_(dict.lookup("n")),
    z_(dict.lookup("z")),
    z0_(readScalar(dict.lookup("z0"))),
    kappa_(dict.lookupOrDefault<scalar>("kappa", 0.41)),
    Uref_(readScalar(dict.lookup("Uref"))),
    Href_(readScalar(dict.lookup("Href"))),
    zGround_("zGround", dict, p.size())
so the Ustart imposed into the ABLconditions file is not taken into account?

and apparently it is recomputed afterwards:
Code:
    Ustar_ = kappa_*Uref_/(log((Href_  + z0_)/max(z0_ , 0.001)));
vinz is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On



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