petebachant / ntnu-hawt-turbinesfoam

Simulating the NTNU horizontal axis wind turbine(s) with the turbinesFoam actuator line model.

NTNU-HAWT-turbinesFoam

Case files for simulating the NTNU horizontal axis wind turbine(s) with the turbinesFoam actuator line model.

The purpose of this simulation is to test the ability to predict:

1. Performance and near-wake of a standalone turbine as measured in Krogstad and Adaramola (2012) "Performance and near wake measurements of a model horizontal axis wind turbine".
2. Performance curves of two in-line turbines as described in Pierella et al. (2014) "Blind Test 2 calculations for two in-line model wind turbines where the downstream turbine operates at various rotational speeds".
3. Wake profiles behind the downwind turbine (in fig. 6) for the case above.
4. Tasks 2 and 3 for the offset turbines described in Krogstad et al. (2015) "' Blind Test 3' calculations of the performance and wake development behind two in-line and offset model wind turbines".

Citation

This case can be cited via DOI. Follow the link above for more information.

License

Code is MIT licensed. See LICENSE for details.

All other materials licensed under a Creative Commons Attribution 4.0 International License.

ntnu-hawt-turbinesfoam's Issues

nut and nuSgs files

Hi,
AFAIK, nuSgs is no longer used for LES, instead nut is used for both, starting from a certain version of OF. Can anyone confirm that?
https://www.cfd-online.com/Forums/openfoam/200655-about-nut-nusgs.html

Plotting functionality not working

From turbinesFoam/turbinesFoam#297 cc @chandukec

• wake -- needs sampling
• spanwise -- file names not formed correctly
• Remove strut-perf -- only relevant to RM2
• recovery -- might want to remove that since it may not be relevant here

Probably happened when switching over to a single run script, and focusing on sampling nacelle values. Basically we need to do similar plotting at https://github.com/petebachant/UNH-RVAT-turbinesFoam.

Plot wake profiles

Should be able to specify behind which turbine. Quantities should be the mean streamwise velocity and the streamwise velocity variance normalized by U_infty**2.

Recheck blade radii

Looks like since they've been offset, the blade roots of the upstream turbine are exposed and creating vortices. May need to extend the cylindrical part inward towards the axis a bit.

Move all run functionality into run.py script

Merge Allrun and paramsweep.py. Similar to https://github.com/UNH-CORE/RM2-CACTUS. Should be able to run parameter sweeps, turn on/off turbines, VGs, set TSRs, etc.

Set height of wind tunnel for zero pressure gradient

According to Pierella et al. (2012) this was set to:

X (m)	Height (m)
0.000	1.801
2.810	1.801
5.621	1.813
8.435	1.842
11.150	1.851

Add args to plotting script

Check LES fvSchemes

ddtSchemes should probably be backward.

Blade data for the turbine

I've searched but could not find the blade geometry and pitching vs free stream velocity data for the NTNU turbine. Can you please share the source, if possible?

Confirm zero pressure gradient

Perform simulation with turbines deactivated. Must finish #1 first.

Check free stream turbulence levels

From Pierella et al. (2012):

the turbulence intensity has been measured to be 0.3%.

Adjust rotor radii

The table in Pierella et al. (2012) seems to be incorrect in that the max radius does not correspond to the rotor radius.

Has anyone been able to validate this? Power coefficient

The experimental results are in blue from Krogstad and Adaramola (2012) "Performance and near wake measurements of a model horizontal axis wind turbine". The numerical results I obtained are in green.

I ran the single turbine case without any changes and was not able to reproduce the experimental results. I also tried refining the mesh later, but it has no influence in this TSR shift. Am I missing something I should be setting before replicating this?

Get S826 foil data at lower Re

For downstream turbine.

Verification

Should be done for both RANS and LES.

Set time step proportional to TSR

Keep steps per rev constant when running, e.g., performance curves.

Problem in Solving

I tried to run the case but I got this Error in pimpleFoam log file. what could be the reason?

((` Creating finite volume options from "system/fvOptions"

Selecting finite volume options model type axialFlowTurbineALSource
Source: turbine1
- selecting cells using cellSet turbine1

--> FOAM FATAL ERROR:
Cannot find directory "polyMesh/sets" in times 0 down to constant

From function Time::findInstance(const fileName&, const word&, const IOobject::readOption, const word&)
in file db/Time/findInstance.C at line 142.

FOAM exiting `))
I have changed the Debug switches in controlDict . does it have any effects on this issu?

Calculate mean fields via function object

Also adjust wake sampling and loading accordingly. See https://github.com/petebachant/UNH-RVAT-turbinesFoam.