UMEP-dev / UMEP-processing

7 stars 9 forks source link

Urock Analyzer output #41

Open jessikalonn opened 8 months ago

jessikalonn commented 8 months ago

Hi, I get a wierd output when i runt URock. I have a dummy area with a building running N-S and wind coming from SW. The building is 6m high: MicrosoftTeams-image

I draw a line near the corner and the output shows wierd wind changes far above the building and the building divided in two: _line1

When I run with vegetation i also dosn't seem to get the vegetation in the vertical wind profile.

Output folder: Urock.zip

Data used:

URockDataIssue.zip

j3r3m1 commented 8 months ago

Thank you for reporting. The wind change is unfortunately more a URock calc problem than a URock analyze problem, coming from the rooftop corner scheme. image image

The building divided in two is a bug in the postprocessor. You should consider that there is a continuity of building (and I should consider fixing this in the future).

Concerning vegetation, it is not yet included in the vertical wind profile.

We might choose to turn off the postprocessing since it is not fully operational but in the mean time it is useful to identify bugs in the vertical extend...

biglimp commented 8 months ago

Thanks Jeremy. So the reason why the distortion to the profile is extended to unreasonable high height is because the building is "long". If we out the profile more in the middle it would look ok?

Maybe we should just state that the URock Postprocessor is "under development"...

Off topic: Did you submit the article poofs?

j3r3m1 commented 8 months ago

Thanks Jeremy. So the reason why the distortion to the profile is extended to unreasonable high height is because the building is "long". If we out the profile more in the middle it would look ok?

Yes. But in my opinion it shouldn't be a problem to have long building. The denominator of the $H{ccp}$ equation should counter balance the $y{L_p}$ on the numerator. But I have not made a test for that (neither manually, neither in the code). First I would check the Equation and second the code if the equation is right. But no time right now...

Maybe we should just state that the URock Postprocessor is "under development"...

Right ! Possible adding a flag in the QGIS plug-in ?

Off topic: Did you submit the article poofs?

I should resent, GMD cannot read my comments. Can you ?

biglimp commented 8 months ago

Thanks Jeremy. So the reason why the distortion to the profile is extended to unreasonable high height is because the building is "long". If we out the profile more in the middle it would look ok?

Yes. But in my opinion it shouldn't be a problem to have long building. The denominator of the Hccp equation should counter balance the yLp on the numerator. But I have not made a test for that (neither manually, neither in the code). First I would check the Equation and second the code if the equation is right. But no time right now...

Maybe we should just state that the URock Postprocessor is "under development"...

Right ! Possible adding a flag in the QGIS plug-in ? No, but some text to inform users. An experimental flag would be set fro the whole plugin and we dont want that.

Off topic: Did you submit the article poofs?

I should resent, GMD cannot read my comments. Can you ? No but I didnt check.

biglimp commented 8 months ago

Jessika, try to put line in the middle of the building.

biglimp commented 8 months ago

We have students that put the line in the middle of the building and it also produce strange output with a vertex that is many times higher than it should be. linje 270 till fredrik 270_line1

j3r3m1 commented 8 months ago

What is strange is the first layer of air above the roof. Can you do the simulation with a much smaller vertical resolution ? It seems you used something like 4 or five while your buildings are less than 5...

biglimp commented 8 months ago

More tests with a simpler setup: image

Modelinput.zip

1 meter resolution image

2 meter resolution 11_line1

4 meter resolution _line1

j3r3m1 commented 8 months ago

Thanks for reporting. In my opinion results are consistent with what is expected by the models whenever the vertical resolution is lower than your building high. The rooftop zone is high but consistent with the Equation: image

where $W_{eff}$ the building width (cross-wind) and $H_F$ the building height

There is clearly no limit of height using this Equation while it should probably be a threshold above a certain width. For the building shown I am not that shocked by the amplitude. Now if we have measurement evidence that the zone should be lower for such large buildings we can of course update the scheme but first we need measurement data for that.. I know the QUIC-URB developpers were working on updating their models dedicated to large buildings. We may adress this issue to them to know what their finding about it ?

biglimp commented 8 months ago

Ok, so the width and height of the building is giving this result? I increased the height of the building to 20 and there is still a rather large vortex above the building. Are you saying that if I had a building that is not so wide, the vortex would be lower?

image

j3r3m1 commented 8 months ago

There is no real vortex in this case above the buildings but a clear decrease of the wind speed for sure. And indeed if the building is less wide the vortex will be lower. Here is what we got for a square based building: image

biglimp commented 8 months ago

aha... then we have an explanation. @jessikalonn and @nilswallenberg , note to us. See the dicussion above. Now we know why the votrex above our kindergarten is so high. MAybe we should make ti more squared for next year... I will close this now.

j3r3m1 commented 8 months ago

Still it might be interesting to know more about the real wind for such large building since it seems the Röckle scheme is not quite appropriate for them and that it concerns quite de lot of buildings (at least in Sweden).