## turbine geometry

Provide feedback, request enhancements, and get help with wind-turbine computer-aided engineering tools.

Moderators: Bonnie.Jonkman, Jason.Jonkman

Jordi Serret
Posts: 11
Joined: Fri Dec 30, 2016 4:36 am
Organization: University of Edinburgh
Location: United Kingdom

### turbine geometry

Dear Bonnie/Jason,

I have a question regarding how to define the geometry of the turbine.

Please, have a look at the attached drawing.

levenmouth.png (165.35 KiB) Viewed 3371 times

TowerHt would be B? (B is 107m)
TowerBsHt would be (C=9m)+18m= 27m?

If so, should I provide the tower distributed properties from 0 to 107m? Or only for the 80m (107-27) of 'real' tower?
If I give values only for the 80m tower, I will violate the inequality TipRad<TowerHt+Twr2Shft+(OverHang*sin(ShftTilt). My blades are extraordinarily large and flexible.
If I give values for the 107m tower (let's imagine I can figure out how), is there any problem later in SubDyn when I give the full description of the substructure?

My turbine is a kind of nearshore or quasi-offshore configuration (see the photography). It has water (0.5m or so, and there are include into the 18m in the drawing) only during the high tide. In this case, and having in account that the wave loads are going to be small, would be better to model it as an onshore turbine?
levenmouth2.png (283.43 KiB) Viewed 3371 times

If so, how that will affect the previous questions?

Kind Regards,
Jordi

Jason.Jonkman
Posts: 5771
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

### Re: turbine geometry

Dear Jordi,

Based on your drawing, TowerHt in ElastoDyn should be 107 m, TowerBsHt in ElastoDyn should be 27 m, the distributed tower properties should be specified along 80 m (that is the tower base is at 27 m above mean sea level (MSL), and the tower top is 107 m above MSL), and the SubDyn model should extend up to 27 m above MSL. HydroDyn can be enabled/disabled independent from SubDyn, so, you can choose to enable HydroDyn for only cases where hydrodynamic loads are important.

I hope that helps.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Jordi Serret
Posts: 11
Joined: Fri Dec 30, 2016 4:36 am
Organization: University of Edinburgh
Location: United Kingdom

### Re: turbine geometry

Dear Jason,

Giving distributed values for the 80m, and setting up TowerHt=107m gives me the next error message:

Code: Select all

`FAST_InitializeAll:InitModuleMappings:ED_L_2_AD_L_T:MeshMapCreate:CreateMotionMap_L2_to_L2:CreateMapping_ProjectToLine2:Node 1 does not project onto any line2 element. FAST encountered an error during module initialization. Simulation error level: FATAL ERROR Aborting FAST.`

Could be because I didn't set up SubDyn yet?
In AeroDyn15, I have provided distributed tower properties from 0 to 80m. Is it correct?

Kind regards,
Jordi

Jason.Jonkman
Posts: 5771
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

### Re: turbine geometry

Dear Jordi,

That mesh-mapping error has been discussed in the following forum topic: viewtopic.php?f=13&t=1747&p=8378. Basically, AeroDyn can only provide aerodynamic loads on the tower; it cannot map aerodynamic loads below the tower or onto a substructure modeled in SubDyn.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Achim.Struve
Posts: 12
Joined: Thu Aug 06, 2015 2:10 am
Organization: Wind Energy Technology Institute of Fh Flensburg
Location: Germany - Schleswig-Holstein

### Re: turbine geometry

Dear Jason,

in case I would like to model a land-based full lattice tower and based on Jordi's drawing in the beginning of this topic, TowerHt should be equal to TowerBsHt, shouldn't it?

If so, I get the following error:

Code: Select all

`FAST_InitializeAll:ED_Init:ED_ValidateInput:ValidatePrimaryData:TowerBsHt must be less than TowerHt.`

Furthermore I tried to overcome this by setting TowerHt to a slightly greater value, such as TowerHt = TowerHt +0.001, but then the next error arose:

Code: Select all

`CalcOutputs_And_SolveForInputs:SolveOption1:FullOpt1_InputOutputSolve:LAPACK_SGETRF: U(7,      7)=0. Facut U is exactly singular.`

As prescribed in the SubDyn Manual on pages 29-30, all 6 platform DOF's have been enabled and all tower bending DOF's have been disabled. Here I show you the relevant part in my ElastoDyn file, where the full lattice tower should have the same hight as the reference 5MW NREL tower:

Code: Select all

`     87.601   TowerHt     - Height of tower above ground level [onshore] or MSL [offshore] (meters)       87.6   TowerBsHt   - Height of tower base above ground level [onshore] or MSL [offshore] (meters)          0   PtfmCMxt    - Downwind distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters)          0   PtfmCMyt    - Lateral distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters)       87.6   PtfmCMzt    - Vertical distance from the ground level [onshore] or MSL [offshore] to the platform CM (meters)       87.6   PtfmRefzt   - Vertical distance from the ground level [onshore] or MSL [offshore] to the platform reference point (meters)`

Probably I misunderstood or forgot something. Can you help me, please?

Thank you and best regards,
Achim

Jason.Jonkman
Posts: 5771
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

### Re: turbine geometry

Dear Achim,

I'm not sure, but it sounds like the Jacobian relating the loads and accelerations between SubDyn, ElastoDyn (and perhaps HydroDyn if you have that enabled) calculated within FAST is singular for some reason e.g. because a mass or inertia is specified as zero. You can look at the Jacobian directly by recompiling FAST with compiler directive OUTPUT_JACOBIAN defined and then rerunning the simulation. I suggest starting there to debug.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Achim.Struve
Posts: 12
Joined: Thu Aug 06, 2015 2:10 am
Organization: Wind Energy Technology Institute of Fh Flensburg
Location: Germany - Schleswig-Holstein

### Re: turbine geometry

Dear Jason,

thank you so much for your helpful reply. I did what you said and found out that the Jacobian has been singular, first. Therefore I introduced artificial inertias to the platform:

Code: Select all

`2.60789E+05   PtfmRIner   - Platform inertia for roll tilt rotation about the platform CM (kg m^2)2.60789E+05   PtfmPIner   - Platform inertia for pitch tilt rotation about the platform CM (kg m^2)2.60789E+05   PtfmYIner   - Platform inertia for yaw rotation about the platform CM (kg m^2)`

This solved the problem and the Jacobian is not singular anymore. Furthermore the loads seem to be consistent in terms of this full lattice configuration. What I mean by "consistent" is that the Yaw loads ≈ Tower Base Loads ≈ SD Interface Loads. In this case I arbitrary chose the artificial platform inertias to be 10% of the nacelle yaw inertia and I am wondering, if these inertias would have a significant influence to the simulation results. This questions comes up especially because I run into a "small angle violation error due to a large platform displacement" for smaller artificial platform inertia values. Maybe there is another way, do you have an idea?

Best regards,
Achim

Jason.Jonkman
Posts: 5771
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

### Re: turbine geometry

Dear Achim,

My understanding is that you've placed the "platform" at the tower-top. So, effectively, the platform is now coincident with the nacelle.

Must you apply the same inertia in all directions (roll, pitch, yaw)? In yaw, I would guess that you'd only to specify a nonzero platform inertia if you've also enabled the nacelle-yaw DOF; otherwise, I would think the nacelle-yaw inertia that you've specified would suffice. In roll and pitch, neither of these inertias are specified for the nacelle (and they are assumed to be zero within ElastoDyn); perhaps this is why you need nonzero values for those.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Achim.Struve
Posts: 12
Joined: Thu Aug 06, 2015 2:10 am
Organization: Wind Energy Technology Institute of Fh Flensburg
Location: Germany - Schleswig-Holstein

### Re: turbine geometry

Dear Jason,

Jason.Jonkman wrote:My understanding is that you've placed the "platform" at the tower-top. So, effectively, the platform is now coincident with the nacelle.

Exactly.

Jason.Jonkman wrote:Must you apply the same inertia in all directions (roll, pitch, yaw)? In yaw, I would guess that you'd only to specify a nonzero platform inertia if you've also enabled the nacelle-yaw DOF

Yes, if I disable the yaw-DOF, only roll, and pitch platform inertias are required. Enabling the nacelle-yaw DOF requires additional artificial platform inertia around the yaw axis.

Jason.Jonkman wrote:; otherwise, I would think the nacelle-yaw inertia that you've specified would suffice. In roll and pitch, neither of these inertias are specified for the nacelle (and they are assumed to be zero within ElastoDyn); perhaps this is why you need nonzero values for those.

Okay, but for me, artificial platform inertias of >10e4 kg m² seem to have no negliable influence to the whole system. I assume that they would reduce the global bending and torsional natural frequencies of the system integrated tower, which would be an assumption on the safe side in terms of following a stiff-stiff tower design philosophy, wouldn't it?

Best regards,
Achim

Jason.Jonkman
Posts: 5771
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

### Re: turbine geometry

Dear Achim,

Did you mean to say that the artificial platform inertia has negligible influence on the whole system?

Yes, reducing the natural frequency of a stiff-stiff tower would be conservative.

Regardless, I was implying that you don't specify an "artificial" inertia, but instead apply the "real" inertia of the nacelle about the roll and pitch axes that are currently neglected by ElastoDyn.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Achim.Struve
Posts: 12
Joined: Thu Aug 06, 2015 2:10 am
Organization: Wind Energy Technology Institute of Fh Flensburg
Location: Germany - Schleswig-Holstein

### Re: turbine geometry

Dear Jason,

Jason.Jonkman wrote:Did you mean to say that the artificial platform inertia has negligible influence on the whole system?

I am not sure - that was my question. Furthermore I assumed that artificial platform inertias would decrease the bending and torsional natural frequencies. Therefore it would be conservative to introduce those artificial platfrom inertias for a stiff-stiff tower.

Jason.Jonkman wrote:Regardless, I was implying that you don't specify an "artificial" inertia, but instead apply the "real" inertia of the nacelle about the roll and pitch axes that are currently neglected by ElastoDyn.

I would like to do so, but do we have those pitch and roll inertias of the nacelle? Neither the 5MW NREL definition document (https://www.nrel.gov/docs/fy09osti/38060.pdf page 12-13) nor the pre-design of the Dowec 6 MW (https://www.ecn.nl/fileadmin/ecn/units/wind/docs/dowec/10046_009.pdf page 26-27) proposes them.

Thanks and best regards,
Achim

Jason.Jonkman
Posts: 5771
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

### Re: turbine geometry

Dear Achim,

No, the roll and pitch inertias of the nacelle where never specified for the NREL 5-MW turbine. Using artificial values as you propose is probably best.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Jharna.Pokhrel
Posts: 15
Joined: Mon Jul 11, 2016 2:49 pm
Organization: South Dakota State University
Location: South Dakota, USA

### Re: turbine geometry

Hi Jason,
I came across some problem relating to hub height. I have noticed that my simulation is aborted when hub-height is around 80s.
I have reduced the grid-height in TurbSim to a small value(=2) which took a good amount of time to terminate TurbSim input file. But when I run the FAST input file, I get the following error. I know such kind of error has already been addressed in another forum, but I couldn't resolve this issue by any means.
Do you think taking such a small hub height causes such kind of error?
I have attached error I am receiving when hub height is 85m and rotor diameter is 137m.

Regards,
Jharna
Attachments
Capture.PNG (18.71 KiB) Viewed 2959 times

Jason.Jonkman
Posts: 5771
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

### Re: turbine geometry

Dear Jharna,

The grid height in TurbSim should be set larger than the rotor diameter (in your case larger than 137 m); I'm not sure why you have set it too a small value.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Jharna.Pokhrel
Posts: 15
Joined: Mon Jul 11, 2016 2:49 pm
Organization: South Dakota State University
Location: South Dakota, USA

### Re: turbine geometry

Thank you so much. The model is working fine when GridHeight is increased to 140 and NumGrid_Z reduced to 10.
I might have misinterpreted to what has been discussed.

Thanks,
Jharna