Skip to content

Commit 06d2505

Browse files
deslaughterdeslaughter
authored
Merge pull request OpenFAST#2213 from hkross/feature/MHK_AddedMass_Inertia
MHK wave-current superposition, added mass, inertia
2 parents 2b3d36c + 570c0ec commit 06d2505

47 files changed

Lines changed: 1951 additions & 310 deletions

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

‎docs/OtherSupporting/AeroDyn/Fluid_Inertia_Added_Mass_Implementation_Plan.docx‎

46.5 KB
Binary file not shown.

‎docs/OtherSupporting/AeroDyn/UMERC24-Poster.pdf‎

2.68 MB
Binary file not shown.

‎docs/source/user/aerodyn-aeroacoustics/example/AeroDyn.ipt‎

Lines changed: 8 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -3,18 +3,15 @@ IEA Wind Task 37 land-based reference wind turbine
33
====== General Options ===================================================
44
False Echo - Echo the input to "<rootname>.AD.ech"? (flag)
55
"default" DTAero - Time interval for aerodynamic calculations {or "default"} (s)
6-
1 WakeMod - Type of wake/induction model (switch) {0=none, 1=BEMT}
7-
2 AFAeroMod - Type of blade airfoil aerodynamics model (switch
8-
0 TwrPotent - Type of tower influence on wind around the tower (switch)
9-
0 TwrShadow - Type of tower influence on wind based on downstream tower shadow (switch) {0=none, 1=Powles model, 2=Eames model}
6+
1 Wake_Mod - Wake/induction model (switch) {0=none, 1=BEMT, 3=OLAF} [Wake_Mod cannot be 2 or 3 when linearizing]
7+
0 TwrPotent - Type tower influence on wind based on potential flow around the tower (switch) {0=none, 1=baseline potential flow, 2=potential flow with Bak correction}
8+
0 TwrShadow - Calculate tower influence on wind based on downstream tower shadow (switch) {0=none, 1=Powles model, 2=Eames model}
109
False TwrAero - Calculate tower aerodynamic loads? (flag)
11-
False FrozenWake - Assume frozen wake during linearization? (flag
12-
False CavitCheck - Perform cavitation check? (flag)
13-
False Buoyancy - Include buoyancy effects? (flag)
10+
False CavitCheck - Perform cavitation check? (flag) [UA_Mod must be 0 when CavitCheck=true]
1411
False NacelleDrag - Include Nacelle Drag effects? (flag)
15-
True CompAA - Flag to compute AeroAcoustics calculation
16-
"AeroAcousticsInput.dat" AA_InputFile
17-
====== Environmental Conditions ==========================================
18-
"default" AirDens - Air density (kg/m^3)
12+
True CompAA - Flag to compute AeroAcoustics calculation [used only when Wake_Mod = 1 or 2]
13+
"AeroAcousticsInput.dat" AA_InputFile - AeroAcoustics input file [used only when CompAA=true]
14+
====== Environmental Conditions ===================================================================
15+
"default" AirDens - Air density (kg/m^3)
1916

2017
File continues...

‎docs/source/user/aerodyn/appendix.rst‎

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -21,11 +21,11 @@ to vary with time. This feature can be useful for debugging the aerodynamic resp
2121
outside of OpenFAST.
2222

2323
2) Multi-rotor AeroDyn Driver Input File
24-
:download:`(driver input file example) <examples/ad_driver_multiple.dvr>`:
24+
:download:`(driver input file example) <examples/ad_driver_multiple.dvr>`
2525

2626

2727
3) AeroDyn Primary Input File
28-
:download:`(primary input file example) <examples/ad_primary_example.dat>`:
28+
:download:`(primary input file example) <examples/ad_primary_example.dat>`
2929

3030
The primary AeroDyn input file defines modeling options, environmental conditions (except freestream flow), airfoils, tower nodal discretization and properties, tower, hub, and nacelle properties, as well as output file specifications.
3131

@@ -35,14 +35,14 @@ The input file begins with two lines of header information which is for your use
3535

3636
4) Airfoil Data Input File
3737

38-
:download:`(profile data) <examples/ad_polar_example.dat>`:
38+
:download:`(profile data) <examples/ad_polar_example.dat>`
3939

40-
:download:`(profile coordinates) <examples/ad_airfoil_example.dat>`:
40+
:download:`(profile coordinates) <examples/ad_airfoil_example.dat>`
4141

4242
The airfoil data input files themselves (one for each airfoil) include tables containing coefficients of lift force, drag force, and pitching moment versus AoA, as well as UA model parameters. In these files, any line whose first non-blank character is an exclamation point (!) is ignored (for inserting comment lines). The non-comment lines should appear within the file in order, but comment lines may be intermixed as desired for reading clarity.
4343

4444
5) Blade Data Input File
45-
:download:`(blade data input file example) <examples/ad_blade_example.dat>`:
45+
:download:`(blade data input file example) <examples/ad_blade_example.dat>`
4646

4747
The blade data input file contains the nodal discretization, geometry, twist, chord, airfoil identifier, and buoyancy properties for a blade. Separate files are used for each blade, which permits modeling of aerodynamic imbalances.
4848

‎docs/source/user/aerodyn/driver.rst‎

Lines changed: 32 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -118,6 +118,21 @@ An example of inputs is given below:
118118
119119
120120
121+
**SeaState data**
122+
123+
SeaState can be called by the AeroDyn driver to define a wave field as part of the inflow information.
124+
For MHK turbines with waves and currents, SeaState will query InflowWind and sum the velocities and
125+
accelerations from the wave and current fields. If SeaState is activated, InflowWind must also be
126+
activated, though the current can be set to 0 if desired. An example of inputs for this section is given below:
127+
128+
.. code::
129+
130+
----- SeaState Data ---------------------------------------------------------------------
131+
1 CompSeaSt - Compute wave velocities (switch) {0=No Waves; 1=SeaState}
132+
"MHK_RM1_Floating_SeaState.dat" SeaStFile - Name of the SeaState input file [used only when CompSeaSt=1]
133+
134+
135+
121136
**Turbine data**
122137

123138
The user specifies the number of turbines as follows:
@@ -381,19 +396,19 @@ An example is given below for two turbines:
381396

382397
.. code::
383398
384-
====== Hub Properties ============================================================================== [used only when Buoyancy=True]
399+
====== Hub Properties ============================================================================== [used only when MHK=1 or 2]
385400
7.0 VolHub - Hub volume (m^3)
386401
0.0 HubCenBx - Hub center of buoyancy x direction offset (m)
387-
====== Hub Properties ============================================================================== [used only when Buoyancy=True]
402+
====== Hub Properties ============================================================================== [used only when MHK=1 or 2]
388403
5.0 VolHub - Hub volume (m^3)
389404
0.2 HubCenBx - Hub center of buoyancy x direction offset (m)
390-
====== Nacelle Properties ========================================================================== [used only when Buoyancy=True or NacelleDrag=True]
405+
====== Nacelle Properties ========================================================================== [used only when MHK=1 or 2 or when NacelleDrag=True]
391406
32.0 VolNac - Nacelle volume (m^3)
392407
0.3, 0.0, 0.05 NacCenB - Position of nacelle center of buoyancy from yaw bearing in nacelle coordinates (m)
393408
4.67, 20.15, 20.15 NacArea - Projected area of the nacelle in X, Y, Z in the nacelle coordinate system (m^2)
394409
0.5, 0.5, 0.5 NacCd - Drag coefficient for the nacelle areas defined above (-)
395410
0.43, 0, 0 NacDragAC - Position of aerodynamic center of nacelle drag in nacelle coordinates (m)
396-
====== Nacelle Properties ========================================================================== [used only when Buoyancy=True or NacelleDrag=True]
411+
====== Nacelle Properties ========================================================================== [used only when MHK=1 or 2 or when NacelleDrag=True]
397412
32.0 VolNac - Nacelle volume (m^3)
398413
0.3, 0.0, 0.05 NacCenB - Position of nacelle center of buoyancy from yaw bearing in nacelle coordinates (m)
399414
4.67, 20.15, 20.15 NacArea - Projected area of the nacelle in X, Y, Z in the nacelle coordinate system (m^2)
@@ -412,14 +427,14 @@ An example is given below for two turbines:
412427

413428
.. code::
414429
415-
====== Turbine(1) Tower Influence and Aerodynamics ================================================ [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True]
416-
2 NumTwrNds - Number of tower nodes used in the analysis (-) [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True]
430+
====== Turbine(1) Tower Influence and Aerodynamics ================================================ [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or MHK=1 or 2]
431+
2 NumTwrNds - Number of tower nodes used in the analysis (-) [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or MHK=1 or 2]
417432
TwrElev TwrDiam TwrCd TwrTI TwrCb
418433
(m) (m) (-) (-) (-)
419434
0.0 2.0 1.0 0.1 0.0
420435
10.0 1.0 1.0 0.1 0.0
421-
====== Turbine(2) Tower Influence and Aerodynamics ================================================ [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True]
422-
3 NumTwrNds - Number of tower nodes used in the analysis (-) [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True]
436+
====== Turbine(2) Tower Influence and Aerodynamics ================================================ [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or MHK=1 or 2]
437+
3 NumTwrNds - Number of tower nodes used in the analysis (-) [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or MHK=1 or 2]
423438
TwrElev TwrDiam TwrCd TwrTI TwrCb
424439
(m) (m) (-) (-) (-)
425440
0.0 4.0 1.0 0.1 0.0
@@ -466,13 +481,16 @@ An example of an AeroDyn driver for a basic inflow, basic HAWT, and combined cas
466481
1.700000000000000e+03 Pvap - Vapour pressure of working fluid (Pa) [used only for an MHK turbine cavitation check]
467482
0 WtrDpth - Water depth (m)
468483
----- Inflow Data -----------------------------------------------------------------------
469-
0 CompInflow - Compute inflow wind velocities (switch) {0=Steady Wind; 1=InflowWind}
470-
"unused" InflowFile - Name of the InflowWind input file [used only when CompInflow=1]
471-
9.0 HWindSpeed - Horizontal wind speed [used only when CompInflow=0 and AnalysisType=1] (m/s)
472-
140 RefHt - Reference height for horizontal wind speed [used only when CompInflow=0] (m)
473-
0.10 PLExp - Power law exponent [used only when CompInflow=0 and AnalysisType=1] (-)
484+
0 CompInflow - Compute inflow wind velocities (switch) {0=Steady Wind; 1=InflowWind}
485+
"unused" InflowFile - Name of the InflowWind input file [used only when CompInflow=1]
486+
9.0 HWindSpeed - Horizontal wind speed [used only when CompInflow=0 and AnalysisType=1] (m/s)
487+
140 RefHt - Reference height for horizontal wind speed [used only when CompInflow=0] (m)
488+
0.10 PLExp - Power law exponent [used only when CompInflow=0 and AnalysisType=1] (-)
489+
----- SeaState Data ---------------------------------------------------------------------
490+
0 CompSeaSt - Compute wave velocities (switch) {0=No Waves; 1=SeaState}
491+
"unused" SeaStFile - Name of the SeaState input file [used only when CompSeaSt=1]
474492
----- Turbine Data ----------------------------------------------------------------------
475-
1 NumTurbines - Number of turbines
493+
1 NumTurbines - Number of turbines
476494
----- Turbine(1) Geometry ---------------------------------------------------------------
477495
True BasicHAWTFormat(1) - Flag to switch between basic or generic input format {True: next 7 lines are basic inputs, False: Base/Twr/Nac/Hub/Bld geometry and motion must follow}
478496
0,0,0 BaseOriginInit(1) - Coordinate of tower base in base coordinates (m)

‎docs/source/user/aerodyn/examples/ad_blade_example.dat‎

Lines changed: 25 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -2,29 +2,29 @@
22
Description line for this file -- file corresponds to inputs in Test01_UAE_AeroDyn.dat
33
====== Blade Properties =================================================================
44
23 NumBlNds - Number of blade nodes used in the analysis (-)
5-
BlSpn BlCrvAC BlSwpAC BlCrvAng BlTwist BlChord BlAFID BlCb BlCenBn BlCenBt
6-
(m) (m) (m) (deg) (deg) (m) (-) (-) (m) (m)
7-
0.0 0.0 0.0 0.0 0.000 0.219 1 0.0 0.0 0.0
8-
0.1360 0.0 0.0 0.0 0.000 0.219 1 0.0 0.0 0.0
9-
0.4481 0.0 0.0 0.0 -0.098 0.181 1 0.0 0.0 0.0
10-
0.8001 0.0 0.0 0.0 19.423 0.714 3 0.0 0.0 0.0
11-
1.0767 0.0 0.0 0.0 14.318 0.711 4 0.0 0.0 0.0
12-
1.2779 0.0 0.0 0.0 10.971 0.691 5 0.0 0.0 0.0
13-
1.4958 0.0 0.0 0.0 8.244 0.668 6 0.0 0.0 0.0
14-
1.7137 0.0 0.0 0.0 6.164 0.647 7 0.0 0.0 0.0
15-
1.9149 0.0 0.0 0.0 4.689 0.627 7 0.0 0.0 0.0
16-
2.1160 0.0 0.0 0.0 3.499 0.606 8 0.0 0.0 0.0
17-
2.3340 0.0 0.0 0.0 2.478 0.584 8 0.0 0.0 0.0
18-
2.5520 0.0 0.0 0.0 1.686 0.561 8 0.0 0.0 0.0
19-
2.7530 0.0 0.0 0.0 1.115 0.542 8 0.0 0.0 0.0
20-
2.9542 0.0 0.0 0.0 0.666 0.522 8 0.0 0.0 0.0
21-
3.1721 0.0 0.0 0.0 0.267 0.499 8 0.0 0.0 0.0
22-
3.3900 0.0 0.0 0.0 -0.079 0.478 8 0.0 0.0 0.0
23-
3.5912 0.0 0.0 0.0 -0.381 0.457 9 0.0 0.0 0.0
24-
3.7924 0.0 0.0 0.0 -0.679 0.437 9 0.0 0.0 0.0
25-
3.9684 0.0 0.0 0.0 -0.933 0.419 9 0.0 0.0 0.0
26-
4.1444 0.0 0.0 0.0 -1.184 0.401 10 0.0 0.0 0.0
27-
4.3456 0.0 0.0 0.0 -1.466 0.381 10 0.0 0.0 0.0
28-
4.5216 0.0 0.0 0.0 -1.711 0.363 10 0.0 0.0 0.0
29-
4.5970 0.0 0.0 0.0 -1.711 0.363 10 0.0 0.0 0.0
5+
BlSpn BlCrvAC BlSwpAC BlCrvAng BlTwist BlChord BlAFID t_c BlCb BlCenBn BlCenBt BlCpn BlCpt BlCan BlCat BlCam
6+
(m) (m) (m) (deg) (deg) (m) (-) (-) (-) (m) (m) (-) (-) (-) (-) (-)
7+
0.0 0.0 0.0 0.0 0.000 0.219 1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
8+
0.1360 0.0 0.0 0.0 0.000 0.219 1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
9+
0.4481 0.0 0.0 0.0 -0.098 0.181 1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
10+
0.8001 0.0 0.0 0.0 19.423 0.714 3 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
11+
1.0767 0.0 0.0 0.0 14.318 0.711 4 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
12+
1.2779 0.0 0.0 0.0 10.971 0.691 5 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
13+
1.4958 0.0 0.0 0.0 8.244 0.668 6 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
14+
1.7137 0.0 0.0 0.0 6.164 0.647 7 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
15+
1.9149 0.0 0.0 0.0 4.689 0.627 7 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
16+
2.1160 0.0 0.0 0.0 3.499 0.606 8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
17+
2.3340 0.0 0.0 0.0 2.478 0.584 8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
18+
2.5520 0.0 0.0 0.0 1.686 0.561 8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
19+
2.7530 0.0 0.0 0.0 1.115 0.542 8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
20+
2.9542 0.0 0.0 0.0 0.666 0.522 8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
21+
3.1721 0.0 0.0 0.0 0.267 0.499 8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
22+
3.3900 0.0 0.0 0.0 -0.079 0.478 8 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
23+
3.5912 0.0 0.0 0.0 -0.381 0.457 9 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
24+
3.7924 0.0 0.0 0.0 -0.679 0.437 9 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
25+
3.9684 0.0 0.0 0.0 -0.933 0.419 9 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
26+
4.1444 0.0 0.0 0.0 -1.184 0.401 10 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
27+
4.3456 0.0 0.0 0.0 -1.466 0.381 10 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
28+
4.5216 0.0 0.0 0.0 -1.711 0.363 10 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
29+
4.5970 0.0 0.0 0.0 -1.711 0.363 10 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
3030

‎docs/source/user/aerodyn/examples/ad_driver_example.dvr‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -20,6 +20,9 @@ False Echo - Echo input parameters to "<rootname>.ech"?
2020
9.0 HWindSpeed - Horizontal wind speed [used only when CompInflow=0 and AnalysisType=1] (m/s)
2121
140 RefHt - Reference height for horizontal wind speed [used only when CompInflow=0] (m)
2222
0.10 PLExp - Power law exponent [used only when CompInflow=0 and AnalysisType=1] (-)
23+
----- SeaState Data [used only when MHK = 1 or 2] ---------------------------------------
24+
0 CompSeaSt - Compute wave velocities (switch) {0=No Waves; 1=SeaState}
25+
"unused" SeaStFile - Name of the SeaState input file [used only when CompSeaSt=1]
2326
----- Turbine Data ----------------------------------------------------------------------
2427
1 NumTurbines - Number of turbines
2528
----- Turbine(1) Geometry ---------------------------------------------------------------

‎docs/source/user/aerodyn/examples/ad_driver_multiple.dvr‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -20,6 +20,9 @@ False Echo - Echo input parameters to "<rootname>.ech"?
2020
10 HWindSpeed - Horizontal wind speed [used only when CompInflow=0 and AnalysisType=1] (m/s)
2121
200 RefHt - Reference height for horizontal wind speed [used only when CompInflow=0] (m)
2222
0 PLExp - Power law exponent [used only when CompInflow=0 and AnalysisType=1] (-)
23+
----- SeaState Data [used only when MHK = 1 or 2] ---------------------------------------
24+
0 CompSeaSt - Compute wave velocities (switch) {0=No Waves; 1=SeaState}
25+
"unused" SeaStFile - Name of the SeaState input file [used only when CompSeaSt=1]
2326
----- Turbine Data ----------------------------------------------------------------------
2427
2 NumTurbines - Number of turbines
2528
----- Turbine(1) ------------------------------------------------------------------------

0 commit comments

Comments
 (0)