AeroDyn and FAST

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

Moderators: Bonnie.Jonkman, Jason.Jonkman

Qiang.Liu
Posts: 3
Joined: Mon Sep 26, 2011 9:04 pm
Organization: Institute of Engineering Thermophysics
Location: Beijing, China

AeroDyn and FAST

Postby Qiang.Liu » Mon Oct 17, 2011 11:08 pm

Hello all,
I am new to this forum and a beginner to FAST. Recently, to learn the use of FAST, I downloaded FAST_v7.00.00a-bjj_AeroDyn_v13.00.00a-bjj_BladedDLLInterface.exe and NRELOffshrBsline5MW_TLP.zip from Jason's public folder, and carried out some simulations. I encountered two problems:

1. In the input file of AeroDyn NRELOffshrBsline5MW_AeroDyn.ipt, when the IndModel is set to be SWIRL, the TanInd is always 0.00000E+00, which is inconsistent with AeroDyn User’s Guide, Version 12.50. I downloaded AeroDyn source files (AD_v130000a-bjj.exe) from NREL website, and I found a few lines of code seems wrong, starting from line 424 to 441 in AeroSubs.f90:

Code: Select all

   SELECT CASE ( TRIM(Line) )
      CASE ('NONE')
         Wake = .FALSE.
         Swirl = .FALSE.

         CALL ProgWarn( ' All wake calculations are turned off! This option is recommended only '// &
                        'in high winds or for debugging.' )
      CASE ('WAKE')
         Wake = .TRUE.
         Swirl = .FALSE.
      CASE ('SWIRL')
         Wake = .TRUE.
         Swirl = .TRUE.
      CASE DEFAULT
         CALL ProgWarn( ' Error: Expecting "NONE", "WAKE", or "SWIRL" wake model option.')
         ErrStat = 1
         RETURN
   END SELECT

CASE ('SWIRL')
Wake = .TRUE.
Swirl = .TRUE.

Here, I think Wake = .TRUE. should be Wake = .FALSE.
2. In the platform file NRELOffshrBsline5MW_Platform_TLP.dat, there is a line:
"200.0 WtrDpth - Water depth (meters) [USE THE SAME VALUE SPECIFIED IN THE WAMIT .POT FILE!]"
It seems that I can not change the value of water depth with the HydroData unchanged. However, I do not have a licensed WAMIT.
So, I have two questions:
Can I change the WaveHs and WaveTp with the HydroData unchanged?
Is there other ways for me to change the water depth?
Confucius said: Pig have dream, so can you.

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

Re: AeroDyn and FAST

Postby Jason.Jonkman » Tue Oct 18, 2011 6:24 am

Dear Qiang,

When IndModel is set to SWIRL, AeroDyn enables both axial and tangential induction. In the AeroDyn source code, variable Wake is a flag to enable axial induction and variable Swirl is a flag to enable tangential induction. So, the source looks correct to me. I'm not sure why your model shows zero-valued TanInd. What is InfModel set to? Is TanInd zero for all elements for all time, or just some elements some time?

The radiation and diffraction solutions of WAMIT are typically found with a finite-depth boundary condition. The statement in the HydroDyn input file that you refer to is reminding you of this. For the greatest accuracy, it is recommended that you set the water depth in WAMIT and HydroDyn to be identical. Of course, the finite-depth boundary condition may have only minor influence on the WAMIT hydrodynamic coefficients, depending on the geometry of the floater, sea state, etc. When the dependency is minor, it reasonable to alter the HydroDyn water depth without updating the WAMIT hydrodynamic coefficients. I haven't performed a sensitivity study on water depth for the WAMIT hydrodynamic coefficients of the MIT/NREL TLP, but it is probably reasonable to assume that the coefficients are reasonable for deviations of the water depth around 200 m.

The WAMIT hydrodynamic coefficients are derived independent of the sea state, so, it is possible to set WaveHs and WaveTp in HydroDyn without modifying the WAMIT hydrodynamic coefficients.

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

Qiang.Liu
Posts: 3
Joined: Mon Sep 26, 2011 9:04 pm
Organization: Institute of Engineering Thermophysics
Location: Beijing, China

Re: AeroDyn and FAST

Postby Qiang.Liu » Tue Oct 18, 2011 7:03 pm

Dear Jason,
Thank you so much for you detailed response. It is very helpful for me.
I'm sorry I have misunderstood the code.
What is InfModel set to? Is TanInd zero for all elements for all time, or just some elements some time?

The InfModel is set to DYNIN, and TanInd is zero for all elements for all time. However, I have only tested three cases when the inflow wind speed is 5 m/s, 8 m/s and 12 m/s.
I have just found that when the InfModel is set to EQUIL, and IndModel is set to SWIRL, TanInd will not always be zero.
Confucius said: Pig have dream, so can you.

Bonnie.Jonkman
Posts: 575
Joined: Thu Nov 10, 2005 10:51 am
Organization: Envision Energy USA
Location: Colorado
Location: Boulder, CO
Contact:

Re: AeroDyn and FAST

Postby Bonnie.Jonkman » Thu Oct 20, 2011 11:10 am

Dear Qiang,

It looks like you have found a bug in AeroDyn. After a little checking into the code, we have discovered that while TanInd is calculated, it is not written in the output file when SWIRL is TRUE.

To fix this, you can can add a line in AeroSubs.f90's vindinf() subroutine. The following code starts at line 7284 in AeroDyn v13.00.00a-bjj, with the additional line surrounded by comments:

Code: Select all

   IF ( SWRLARG > 0.0 ) THEN
      A2P = 0.5 * ( -1.0 + SQRT( SWRLARG ) )
      VT  = VT * ( 1.0 + A2P)
!bjj start of proposed change v13.00.01a-bjj
      AP(iradius,iblade) = A2P
!bjj end of proposed change v13.00.01a-bjj
   ENDIF


We have modified our code and will include it in the next release of AeroDyn. Thanks for finding the error.

Regards,
Bonnie
Bonnie Jonkman

Envision Energy USA, 2016-
National Renewable Energy Laboratory, 2003-2016

Qiang.Liu
Posts: 3
Joined: Mon Sep 26, 2011 9:04 pm
Organization: Institute of Engineering Thermophysics
Location: Beijing, China

Re: AeroDyn and FAST

Postby Qiang.Liu » Thu Oct 20, 2011 6:08 pm

Dear Bonnie,

Thank you very much! It's so exciting that I actually found a bug! :D
Confucius said: Pig have dream, so can you.

Yiqing.Xia
Posts: 9
Joined: Mon Mar 11, 2013 9:52 pm
Organization: Shanghai Jiao Tong University
Location: China

Re: AeroDyn and FAST

Postby Yiqing.Xia » Wed Mar 19, 2014 7:01 am

Dear
Can I just change the water depth to 150m? Thank YOU.
200.0 WtrDpth - Water depth (meters) [USE THE SAME VALUE SPECIFIED IN THE WAMIT .POT FILE!]


Jason.Jonkman wrote:Dear Qiang,

When IndModel is set to SWIRL, AeroDyn enables both axial and tangential induction. In the AeroDyn source code, variable Wake is a flag to enable axial induction and variable Swirl is a flag to enable tangential induction. So, the source looks correct to me. I'm not sure why your model shows zero-valued TanInd. What is InfModel set to? Is TanInd zero for all elements for all time, or just some elements some time?

The radiation and diffraction solutions of WAMIT are typically found with a finite-depth boundary condition. The statement in the HydroDyn input file that you refer to is reminding you of this. For the greatest accuracy, it is recommended that you set the water depth in WAMIT and HydroDyn to be identical. Of course, the finite-depth boundary condition may have only minor influence on the WAMIT hydrodynamic coefficients, depending on the geometry of the floater, sea state, etc. When the dependency is minor, it reasonable to alter the HydroDyn water depth without updating the WAMIT hydrodynamic coefficients. I haven't performed a sensitivity study on water depth for the WAMIT hydrodynamic coefficients of the MIT/NREL TLP, but it is probably reasonable to assume that the coefficients are reasonable for deviations of the water depth around 200 m.

The WAMIT hydrodynamic coefficients are derived independent of the sea state, so, it is possible to set WaveHs and WaveTp in HydroDyn without modifying the WAMIT hydrodynamic coefficients.

I hope that helps.

Best regards,

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

Re: AeroDyn and FAST

Postby Jason.Jonkman » Fri Mar 21, 2014 1:00 pm

Dear Yiqing,

You can, but without updating the hydrodynamic coefficients through a WAMIT calculation you run the risk of reduced accuracy of the solution. The hydrodynamic coefficients computed by WAMIT are only independent of depth for deep water. My guess is that the hydrodynamic coefficients wont change much between 200 m and 150 m, but I wouldn't go much shallower than that without updating the coefficients.

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

Mustafa.Vardaroglu
Posts: 25
Joined: Sat Jan 06, 2018 11:45 am
Organization: University of Campania
Location: Italy

Re: AeroDyn and FAST

Postby Mustafa.Vardaroglu » Thu Sep 17, 2020 2:22 pm

[size=100]Dear All,

I am trying to obtain the hydrodynamic coefficients (added mass, radiation damping), first order wave excitation force (magnitude and the phase angle) of MIT/NREL TLP by using NEMOH panel method solver. I am attaching my results where I compared :

- NEMOH solution (Blue data series)
- HydroD solution of SESAM (HydroD data series: red)
- HydroD solution of the same platform except with the draft= 50m. (EB data series: pink)
- I have digitized the values of Denis Matha , where he published in: NREL/SR-500-45891 Figure 16, 17,18 (Matha 2010 data series: dark dashed)
- tlpmit.1 & tlpmit.3 WAMIT outputs that are available in the FAST distribution folder (NREL-WAMIT data series: dark-solid)

With the comparisons, except B33 all looks fine (B11 and B55 is acceptable)
All the analyses considers the same platform geometry, draft height. All the geometries are modeled without he pontoons as Matha (2010) also stated.
I am confused with the comparison between Matha (2010) and WAMIT data series which I except them to be the same. What I am missing? Have any idea?
Thanks in advance for your time, suggestions and help.

Sincerely,
Mustafa
[/size]
Attachments
Radiation_damping.jpg
Radiation_damping.jpg (127.22 KiB) Viewed 1060 times
first_order_wave_excitation_mag.jpg
first_order_wave_excitation_mag.jpg (112.25 KiB) Viewed 1060 times
Added_mass.jpg
Added_mass.jpg (114 KiB) Viewed 1060 times

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

Re: AeroDyn and FAST

Postby Jason.Jonkman » Fri Sep 18, 2020 10:21 am

Dear Mustafa,

I would expect the WAMIT output files provided in the FAST/OpenFAST archives for the MIT/NREL TLP (tlpmit.1, .3, and .hst) would be the same as what is in Matha (2010).

The WAMIT output data is nondimensional, so, I would first check that you've redimensionalized the data properly. I'm questioning particularly the B33 plot, where the B33 term does not tend toward zero at low frequency, even though it should be redimensionalized by multiplying by rho*omega. See the following forum post for information on how to dimensionalize the WAMIT output: viewtopic.php?f=4&t=487&p=1764.

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


Return to “Computer-Aided Engineering Software Tools”

Who is online

Users browsing this forum: No registered users and 1 guest