Question about the ABL results of the LES precursor

[wcha533]

Hi,
First, thank you NREL for making this SOWFA-6 publicly available.

I have a question about the setup for the LES precursor.
I am working on the LES precursor to produce an idealized neutral ABL.
I tested with a coarse mesh (20 x 20 x 20m) in a 3km x 3km x 1km domain. I ran this test for 18000 s physical time and the profiles were averaged during the last 2000 s.

The horizontal velocity, resolved tke and potential temperature profiles from the current dev branch look different from SOWFA-2.x and the previous version of SOWFA-6 as shown in the plots below.

I understand that there are some changes in UEqn and TEqn of the current solver. Is it something about the buoyancy term?

The temperature profile of the current dev version show that it diffuses down very quickly from the strong capping inversion zone and that caused the velocity profile at the strong temperature gradient to speed up.

So I would like to ask if there is any suggestion to setup the parameters properly in the setUp file for this current solver version?

Thank you very much

[mchurchf]

Hi. Could you please send your input files (not the polymesh and anything large) to me at [email protected] and I will help figure out what is going on.