Skip to content

Add python test case for ADI #76

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
andrew-platt merged 19 commits into from
Oct 25, 2022
Merged

Add python test case for ADI #76

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
19 commits
Select commit Hold shift + click to select a range
bd3d157
adi_c_binding: rough in a regression library test of interface [not f…
andrew-platt Jan 6, 2022
a09c4dd
adi_c_binding: pass hub and root orientation info through
andrew-platt Jan 13, 2022
6b61323
adi_c_binding: pass nacelle info to init
andrew-platt Jan 13, 2022
3e0744a
adi_c_binding: split hub and nacelle out from blades, add TransposeDC…
andrew-platt Jan 14, 2022
9d0d059
adi_c_lib: read mesh info from vtk
andrew-platt Jan 25, 2022
c0764d8
adi_c_lib: add reading of vtk motion data, add CalcOutput hub motion …
andrew-platt Feb 3, 2022
a5cb2d3
adi_c_lib: add mesh to calcoutput call
andrew-platt Feb 8, 2022
8b6d8cf
adi_c_lib: add calcoutput calls to adi_c_lib interface
andrew-platt Feb 22, 2022
e00fe7d
adi_c_lib: add updatestates calls to test case
andrew-platt Feb 24, 2022
400076c
adi_c_lib: renamed test case to reflect what it actually is
andrew-platt Feb 24, 2022
a684448
adi_c_lib: match the 5MW_OC3Semi wind conditions, and add outputs for…
andrew-platt Feb 24, 2022
55e0ca7
adi_c_lib: adjust indexing for calls (readability mostly)
andrew-platt Feb 24, 2022
587477c
ADI: add WrVTK_Type option
andrew-platt Jul 25, 2022
1a87ee5
ADI: rename python case
andrew-platt Oct 24, 2022
228cb64
ADI: reduce time of python ADI case to 10 timesteps
andrew-platt Oct 24, 2022
520bc5d
ADI: update py_ad aerodyn input file
andrew-platt Oct 24, 2022
ad617b1
ADI: further simplify py_ad driver case
andrew-platt Oct 24, 2022
314c433
ADI: add AeroProjMod as input variable
andrew-platt Oct 24, 2022
574b929
ADI: update py_ad_ case
andrew-platt Oct 24, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
28 changes: 28 additions & 0 deletions modules/aerodyn/py_ad_5MW_OC4Semi_WSt_WavesWN/NRELOffshrBsline5MW_AeroDyn_blade.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@
------- AERODYN v15.00.* BLADE DEFINITION INPUT FILE -------------------------------------
NREL 5.0 MW offshore baseline aerodynamic blade input properties; note that we need to add the aerodynamic center to this file
====== Blade Properties =================================================================
19 NumBlNds - Number of blade nodes used in the analysis (-)
BlSpn BlCrvAC BlSwpAC BlCrvAng BlTwist BlChord BlAFID
(m) (m) (m) (deg) (deg) (m) (-)
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 1.3308000E+01 3.5420000E+00 1
1.3667000E+00 -8.1531745E-04 -3.4468858E-03 0.0000000E+00 1.3308000E+01 3.5420000E+00 1
4.1000000E+00 -2.4839790E-02 -1.0501421E-01 0.0000000E+00 1.3308000E+01 3.8540000E+00 1
6.8333000E+00 -5.9469375E-02 -2.5141635E-01 0.0000000E+00 1.3308000E+01 4.1670000E+00 2
1.0250000E+01 -1.0909141E-01 -4.6120149E-01 0.0000000E+00 1.3308000E+01 4.5570000E+00 3
1.4350000E+01 -1.1573354E-01 -5.6986665E-01 0.0000000E+00 1.1480000E+01 4.6520000E+00 4
1.8450000E+01 -9.8316709E-02 -5.4850833E-01 0.0000000E+00 1.0162000E+01 4.4580000E+00 4
2.2550000E+01 -8.3186967E-02 -5.2457001E-01 0.0000000E+00 9.0110000E+00 4.2490000E+00 5
2.6650000E+01 -6.7933232E-02 -4.9624675E-01 0.0000000E+00 7.7950000E+00 4.0070000E+00 6
3.0750000E+01 -5.3393159E-02 -4.6544755E-01 0.0000000E+00 6.5440000E+00 3.7480000E+00 6
3.4850000E+01 -4.0899260E-02 -4.3583519E-01 0.0000000E+00 5.3610000E+00 3.5020000E+00 7
3.8950000E+01 -2.9722933E-02 -4.0591323E-01 0.0000000E+00 4.1880000E+00 3.2560000E+00 7
4.3050000E+01 -2.0511081E-02 -3.7569051E-01 0.0000000E+00 3.1250000E+00 3.0100000E+00 8
4.7150000E+01 -1.3980013E-02 -3.4521705E-01 0.0000000E+00 2.3190000E+00 2.7640000E+00 8
5.1250000E+01 -8.3819737E-03 -3.1463837E-01 0.0000000E+00 1.5260000E+00 2.5180000E+00 8
5.4666700E+01 -4.3546914E-03 -2.8909220E-01 0.0000000E+00 8.6300000E-01 2.3130000E+00 8
5.7400000E+01 -1.6838383E-03 -2.6074456E-01 0.0000000E+00 3.7000000E-01 2.0860000E+00 8
6.0133300E+01 -3.2815226E-04 -1.7737470E-01 0.0000000E+00 1.0600000E-01 1.4190000E+00 8
6.1499900E+01 -3.2815226E-04 -1.7737470E-01 0.0000000E+00 1.0600000E-01 1.4190000E+00 8

!bjj: because of precision in the BD-AD coupling, 61.5m didn't work, so I changed it to 61.4999m
6.1500000E+01 -3.2815226E-04 -1.7737470E-01 0.0000000E+00 1.0600000E-01 1.4190000E+00 8
83 changes: 83 additions & 0 deletions modules/aerodyn/py_ad_5MW_OC4Semi_WSt_WavesWN/NRELOffshrBsline5MW_Blade.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,83 @@
------- ELASTODYN V1.00.* INDIVIDUAL BLADE INPUT FILE --------------------------
NREL 5.0 MW offshore baseline blade input properties.
---------------------- BLADE PARAMETERS ----------------------------------------
49 NBlInpSt - Number of blade input stations (-)
0.477465 BldFlDmp(1) - Blade flap mode #1 structural damping in percent of critical (%)
0.477465 BldFlDmp(2) - Blade flap mode #2 structural damping in percent of critical (%)
0.477465 BldEdDmp(1) - Blade edge mode #1 structural damping in percent of critical (%)
---------------------- BLADE ADJUSTMENT FACTORS --------------------------------
1 FlStTunr(1) - Blade flapwise modal stiffness tuner, 1st mode (-)
1 FlStTunr(2) - Blade flapwise modal stiffness tuner, 2nd mode (-)
1.04536 AdjBlMs - Factor to adjust blade mass density (-) !bjj: value for AD14=1.04536; value for AD15=1.057344 (it would be nice to enter the requested blade mass instead of a factor here)
1 AdjFlSt - Factor to adjust blade flap stiffness (-)
1 AdjEdSt - Factor to adjust blade edge stiffness (-)
---------------------- DISTRIBUTED BLADE PROPERTIES ----------------------------
BlFract PitchAxis StrcTwst BMassDen FlpStff EdgStff
(-) (-) (deg) (kg/m) (Nm^2) (Nm^2)
0.0000000E+00 2.5000000E-01 1.3308000E+01 6.7893500E+02 1.8110000E+10 1.8113600E+10
3.2500000E-03 2.5000000E-01 1.3308000E+01 6.7893500E+02 1.8110000E+10 1.8113600E+10
1.9510000E-02 2.5049000E-01 1.3308000E+01 7.7336300E+02 1.9424900E+10 1.9558600E+10
3.5770000E-02 2.5490000E-01 1.3308000E+01 7.4055000E+02 1.7455900E+10 1.9497800E+10
5.2030000E-02 2.6716000E-01 1.3308000E+01 7.4004200E+02 1.5287400E+10 1.9788800E+10
6.8290000E-02 2.7941000E-01 1.3308000E+01 5.9249600E+02 1.0782400E+10 1.4858500E+10
8.4550000E-02 2.9167000E-01 1.3308000E+01 4.5027500E+02 7.2297200E+09 1.0220600E+10
1.0081000E-01 3.0392000E-01 1.3308000E+01 4.2405400E+02 6.3095400E+09 9.1447000E+09
1.1707000E-01 3.1618000E-01 1.3308000E+01 4.0063800E+02 5.5283600E+09 8.0631600E+09
1.3335000E-01 3.2844000E-01 1.3308000E+01 3.8206200E+02 4.9800600E+09 6.8844400E+09
1.4959000E-01 3.4069000E-01 1.3308000E+01 3.9965500E+02 4.9368400E+09 7.0091800E+09
1.6585000E-01 3.5294000E-01 1.3308000E+01 4.2632100E+02 4.6916600E+09 7.1676800E+09
1.8211000E-01 3.6519000E-01 1.3181000E+01 4.1682000E+02 3.9494600E+09 7.2716600E+09
1.9837000E-01 3.7500000E-01 1.2848000E+01 4.0618600E+02 3.3865200E+09 7.0817000E+09
2.1465000E-01 3.7500000E-01 1.2192000E+01 3.8142000E+02 2.9337400E+09 6.2445300E+09
2.3089000E-01 3.7500000E-01 1.1561000E+01 3.5282200E+02 2.5689600E+09 5.0489600E+09
2.4715000E-01 3.7500000E-01 1.1072000E+01 3.4947700E+02 2.3886500E+09 4.9484900E+09
2.6341000E-01 3.7500000E-01 1.0792000E+01 3.4653800E+02 2.2719900E+09 4.8080200E+09
2.9595000E-01 3.7500000E-01 1.0232000E+01 3.3933300E+02 2.0500500E+09 4.5014000E+09
3.2846000E-01 3.7500000E-01 9.6720000E+00 3.3000400E+02 1.8282500E+09 4.2440700E+09
3.6098000E-01 3.7500000E-01 9.1100000E+00 3.2199000E+02 1.5887100E+09 3.9952800E+09
3.9350000E-01 3.7500000E-01 8.5340000E+00 3.1382000E+02 1.3619300E+09 3.7507600E+09
4.2602000E-01 3.7500000E-01 7.9320000E+00 2.9473400E+02 1.1023800E+09 3.4471400E+09
4.5855000E-01 3.7500000E-01 7.3210000E+00 2.8712000E+02 8.7580000E+08 3.1390700E+09
4.9106000E-01 3.7500000E-01 6.7110000E+00 2.6334300E+02 6.8130000E+08 2.7342400E+09
5.2358000E-01 3.7500000E-01 6.1220000E+00 2.5320700E+02 5.3472000E+08 2.5548700E+09
5.5610000E-01 3.7500000E-01 5.5460000E+00 2.4166600E+02 4.0890000E+08 2.3340300E+09
5.8862000E-01 3.7500000E-01 4.9710000E+00 2.2063800E+02 3.1454000E+08 1.8287300E+09
6.2115000E-01 3.7500000E-01 4.4010000E+00 2.0029300E+02 2.3863000E+08 1.5841000E+09
6.5366000E-01 3.7500000E-01 3.8340000E+00 1.7940400E+02 1.7588000E+08 1.3233600E+09
6.8618000E-01 3.7500000E-01 3.3320000E+00 1.6509400E+02 1.2601000E+08 1.1836800E+09
7.1870000E-01 3.7500000E-01 2.8900000E+00 1.5441100E+02 1.0726000E+08 1.0201600E+09
7.5122000E-01 3.7500000E-01 2.5030000E+00 1.3893500E+02 9.0880000E+07 7.9781000E+08
7.8376000E-01 3.7500000E-01 2.1160000E+00 1.2955500E+02 7.6310000E+07 7.0961000E+08
8.1626000E-01 3.7500000E-01 1.7300000E+00 1.0726400E+02 6.1050000E+07 5.1819000E+08
8.4878000E-01 3.7500000E-01 1.3420000E+00 9.8776000E+01 4.9480000E+07 4.5487000E+08
8.8130000E-01 3.7500000E-01 9.5400000E-01 9.0248000E+01 3.9360000E+07 3.9512000E+08
8.9756000E-01 3.7500000E-01 7.6000000E-01 8.3001000E+01 3.4670000E+07 3.5372000E+08
9.1382000E-01 3.7500000E-01 5.7400000E-01 7.2906000E+01 3.0410000E+07 3.0473000E+08
9.3008000E-01 3.7500000E-01 4.0400000E-01 6.8772000E+01 2.6520000E+07 2.8142000E+08
9.3821000E-01 3.7500000E-01 3.1900000E-01 6.6264000E+01 2.3840000E+07 2.6171000E+08
9.4636000E-01 3.7500000E-01 2.5300000E-01 5.9340000E+01 1.9630000E+07 1.5881000E+08
9.5447000E-01 3.7500000E-01 2.1600000E-01 5.5914000E+01 1.6000000E+07 1.3788000E+08
9.6260000E-01 3.7500000E-01 1.7800000E-01 5.2484000E+01 1.2830000E+07 1.1879000E+08
9.7073000E-01 3.7500000E-01 1.4000000E-01 4.9114000E+01 1.0080000E+07 1.0163000E+08
9.7886000E-01 3.7500000E-01 1.0100000E-01 4.5818000E+01 7.5500000E+06 8.5070000E+07
9.8699000E-01 3.7500000E-01 6.2000000E-02 4.1669000E+01 4.6000000E+06 6.4260000E+07
9.9512000E-01 3.7500000E-01 2.3000000E-02 1.1453000E+01 2.5000000E+05 6.6100000E+06
1.0000000E+00 3.7500000E-01 0.0000000E+00 1.0319000E+01 1.7000000E+05 5.0100000E+06
---------------------- BLADE MODE SHAPES ---------------------------------------
0.0622 BldFl1Sh(2) - Flap mode 1, coeff of x^2
1.7254 BldFl1Sh(3) - , coeff of x^3
-3.2452 BldFl1Sh(4) - , coeff of x^4
4.7131 BldFl1Sh(5) - , coeff of x^5
-2.2555 BldFl1Sh(6) - , coeff of x^6
-0.5809 BldFl2Sh(2) - Flap mode 2, coeff of x^2
1.2067 BldFl2Sh(3) - , coeff of x^3
-15.5349 BldFl2Sh(4) - , coeff of x^4
29.7347 BldFl2Sh(5) - , coeff of x^5
-13.8255 BldFl2Sh(6) - , coeff of x^6
0.3627 BldEdgSh(2) - Edge mode 1, coeff of x^2
2.5337 BldEdgSh(3) - , coeff of x^3
-3.5772 BldEdgSh(4) - , coeff of x^4
2.376 BldEdgSh(5) - , coeff of x^5
-0.6952 BldEdgSh(6) - , coeff of x^6


18 changes: 18 additions & 0 deletions modules/aerodyn/py_ad_5MW_OC4Semi_WSt_WavesWN/README.md
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,18 @@
Simple rigid motion test of 3 blade rotor through python interface with turbulent inflow
Based on 5MW\_OC4Semi\_WSt\_WavesWN regression test case

- Rigid rotor, root node defined only
- Steady rotation speed
- Turbulent inflow

NodePositions.txt -- order of points [x,y,z]
- hub point
- blade 1 -- 19 nodes
- blade 2 -- 19 nodes
- blade 3 -- 19 nodes

NodeOrientations.txt -- order of points [r11,r12,r13,r21,r22,r23,r31,r32,r33]
- hub point
- blade 1 -- 19 nodes
- blade 2 -- 19 nodes
- blade 3 -- 19 nodes
Loading