[bathgooner]

Hi Guys,

I'm very new to CFD and OpenFOAM, but I'm interested in learning more about it! I've got a terrain that I'd like to model wind flow over, much like in the tutorial 'turbinesiting' in windSimpleFoam. I've used GPS coordinated to create an .STL file with my terrain on. I've got SnappyHexMesh working correctly in my blockMesh.

I just want some guidance on how to fill out the boundary conditions, mostly U and P in the '0' folder (using the 'turbinesiting' as a template). It would be good to just test my model with some random values for this to see how it runs.

Any help would be much appriciated.

Thanks

[AlanR]

The turbinesiting tutorial sets pressure boundary conditions in the /0 p file - use exactly the same pressure setup for all boundaries. The wind velocity is in the /include ABLCondidtions file (Uref). If you aren't familiar with the other variable names, the parameters in this file set up an inlet velocity profile (Ustar - friction velocity, z0 - atmospheric roughness height, Uref is the wind speed (m/s) at Href - meters above the zGround elevation).

To start, I would copy the entire turbinesiting tutorial, rename it, and use your own mesh. You will need to have z as the vertical axis in the domain, and set zGround as the minimum elevation of your terrain in meters. If you already have a good mesh with your own terrain, you've made a lot of progress! To move your mesh around, go to the directory where you ran snappyHexMesh and copy the /3 polymesh file to the replace the turbinesiting tutorial polymesh directory. You can experiment by changing Uref to increase or decrease inlet wind speed. Good luck,
Alan

[bathgooner]

Hello again. Thanks for the previous help, unfortunately I've been unable to look at it for a while due to other work!

I had a bit more of a play around with it a few days ago. I've managed to run windSimpleFoam on my 2D terrain in my mesh. This is however where my problems lie. The wind, which should flow from left to right over the terrain definitely doesn't look correct when viewing the glyphs on ParaView. I thought it must be something to do with my boundary conditions, but I am as yet, unable to find the source of the problem.

I have attached a few screen shots, one showing the wind I'm getting through the mesh, another of the present boundary conditions and also of the ABL conditions. As I said previously I'm very new to this, hence why my mesh is currently very simple - will be refined once I can get flow running through it!

Any help would be hugely appreciated! Thanks

[hfs]

Originally Posted by bathgooner

Hello again. Thanks for the previous help, unfortunately I've been unable to look at it for a while due to other work!

I had a bit more of a play around with it a few days ago. I've managed to run windSimpleFoam on my 2D terrain in my mesh. This is however where my problems lie. The wind, which should flow from left to right over the terrain definitely doesn't look correct when viewing the glyphs on ParaView. I thought it must be something to do with my boundary conditions, but I am as yet, unable to find the source of the problem.

I have attached a few screen shots, one showing the wind I'm getting through the mesh, another of the present boundary conditions and also of the ABL conditions. As I said previously I'm very new to this, hence why my mesh is currently very simple - will be refined once I can get flow running through it!

Any help would be hugely appreciated! Thanks

Hello!

I am using the same tutorial to help me understand how to implement a rough wall boundary conditions

I tried to change the Ks value in the nut file in 0 folder but it is giving me the same velocity field every time ... is this nutkRoughWallFunction really working??

Please see the attached photos!

Thanks for your explanation!

[JR22]

There is a master's thesis by Benjamin Martinez from Rizo-DTU in 2011. It contains some good information which is relevant to the Turbine Sitting Tutorial. He talks about the version in OpenFOAM 1.7, but I don't think the tutorial has changed much. He gives good information about nutRoughWallFunction.

http://windenergyresearch.org/downlo..._MScthesis.pdf

Best.

[rob3rt 0ng]

Hi all,

Could you please provide a soft copy or a valid URL link of the thesis? I can't seem to find it in google.

Thanks,
Robert

[marco-vs]

Originally Posted by hfs

Hello!

I am using the same tutorial to help me understand how to implement a rough wall boundary conditions

I tried to change the Ks value in the nut file in 0 folder but it is giving me the same velocity field every time ... is this nutkRoughWallFunction really working??

Please see the attached photos!

Thanks for your explanation!

Helo hfs,

I'm pretty new using openfoam, y would like to simulate de ABL using turbinesiting tutorial but, after running "Allrun" script and opening in paraview, it doesn't work, velocity an pressure is always cero. Perhaps you can provide me some tips to run the case.

thanks in advance!
Screenshot from 2019-10-20 14-20-16.jpg

[TianJoshua]

Originally Posted by marco-vs

Helo hfs,

I'm pretty new using openfoam, y would like to simulate de ABL using turbinesiting tutorial but, after running "Allrun" script and opening in paraview, it doesn't work, velocity an pressure is always cero. Perhaps you can provide me some tips to run the case.

thanks in advance!
Attachment 72878

Hi

Do you know how to run the turbine siting now? I am in the trouble for a few days.
Thank you!