Skip to content

Commit

Permalink
ADI: update py_ad_5MW_OC4Semi_WSt_WavesWN case
Browse files Browse the repository at this point in the history 
- The wind data was incorrect in previous version of the test case.
- Updated notes on how to test
- Reduced number of outputs for faster run time
  • Loading branch information
@andrew-platt
andrew-platt committed Apr 24, 2023
1 parent 8421cfb commit cddc9fd
Show file tree
Hide file tree
Showing 68 changed files with 3,183 additions and 3,127 deletions.
19 changes: 4 additions & 15 deletions modules/aerodyn/py_ad_5MW_OC4Semi_WSt_WavesWN/README.md
View file
Original file line number Diff line number Diff line change
@@ -1,18 +1,7 @@
Simple rigid motion test of 3 blade rotor through python interface with turbulent inflow
Based on 5MW\_OC4Semi\_WSt\_WavesWN regression test case
#py\_ad\_5MW\_OC4Semi\_WSt\_WavesWN

- Rigid rotor, root node defined only
- Steady rotation speed
- Turbulent inflow
This test is based on the fully coupled OpenFAST regression test case 5MW\_OC4Semi\_WSt\_WavesWN. The purpose of this test is to verify that if the same input motions from an OpenFAST test were passed through the c-bindings interface to AeroDyn and InflowWind, the same results would be achieved.

NodePositions.txt -- order of points [x,y,z]
- hub point
- blade 1 -- 19 nodes
- blade 2 -- 19 nodes
- blade 3 -- 19 nodes
For this test, the 5MW\_OC4Semi\_WSt\_WavesWN case is run with no tower effects, and slightly simplified aero (see the `ad_primary.dat` in this case). All motions are output for all timesteps from OpenFAST using the `WrVTK` option 3 with all fields.

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
The resulting vtk motion files are then read in by the simple driver with this test case and used to set all mesh motions passed to AD15 through the c-bindings interface from Python. The resulting aero loads should be exactly identical to the full test case.
131 changes: 72 additions & 59 deletions modules/aerodyn/py_ad_5MW_OC4Semi_WSt_WavesWN/ad_primary.dat
View file
Original file line number Diff line number Diff line change
Expand Up @@ -84,70 +84,83 @@ TwrElev TwrDiam TwrCd TwrTI TwrCb !TwrTI used on
8.7600000E+01 3.8700000E+00 1.0000000E+00 1.0000000E-01 0.0
====== Outputs ====================================================================================
True SumPrint - Generate a summary file listing input options and interpolated properties to "<rootname>.AD.sum"? (flag)
1 NBlOuts - Number of blade node outputs [0 - 9] (-)
19 BlOutNd - Blade nodes whose values will be output (-)
3 NBlOuts - Number of blade node outputs [0 - 9] (-)
1, 9, 19 BlOutNd - Blade nodes whose values will be output (-)
0 NTwOuts - Number of tower node outputs [0 - 9] (-)
1, 2, 6 TwOutNd - Tower nodes whose values will be output (-)
OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx for a listing of available output channels, (-)
"RtAeroCp"
"RtAeroCt"
"RtAeroCq"
"RtTSR"
"RtSpeed, RtSkew, B1Azimuth"
"B1N1VDisx, B1N1VDisy, B1N1VDisz" ! disturbed wind velocity at Blade 1, Node 1
"B1N2VDisx, B1N2VDisy, B1N2VDisz" ! disturbed wind velocity at Blade 1, Node 2
"B1N3VDisx, B1N3VDisy, B1N3VDisz" ! disturbed wind velocity at Blade 1, Node 3
"B1N1STVx, B1N1STVy, B1N1STVz" ! structural translational velocity at Blade 1, Node 1
"B1N2STVx, B1N2STVy, B1N2STVz" ! structural translational velocity at Blade 1, Node 2
"B1N3STVx, B1N3STVy, B1N3STVz" ! structural translational velocity at Blade 1, Node 2
"B1N1AxInd, B1N2AxInd, B1N3AxInd"
"B1N1Alpha, B1N2Alpha, B1N3Alpha"
"B1N1Theta, B1N2Theta, B1N3Theta"
END of input file (the word "END" must appear in the first 3 columns of this last OutList line)
---------------------- NODE OUTPUTS --------------------------------------------
3 BldNd_BladesOut - Blades to output
99 BldNd_BlOutNd - Blade nodes on each blade (currently unused)
OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx, AeroDyn_Nodes tab for a listing of available output channels, (-)
"VUndx" - x-component of undisturbed wind velocity at each node
"VUndy" - y-component of undisturbed wind velocity at each node
"VUndz" - z-component of undisturbed wind velocity at each node
"VDisx" - x-component of disturbed wind velocity at each node
"VDisy" - y-component of disturbed wind velocity at each node
"VDisz" - z-component of disturbed wind velocity at each node
"STVx" - x-component of structural translational velocity at each node
"STVy" - y-component of structural translational velocity at each node
"STVz" - z-component of structural translational velocity at each node
"VRel" - Relvative wind speed at each node
"DynP" - Dynamic pressure at each node
"Re" - Reynolds number (in millions) at each node
"M" - Mach number at each node
"Vindx" - Axial induced wind velocity at each node
"Vindy" - Tangential induced wind velocity at each node
"AxInd" - Axial induction factor at each node
"TnInd" - Tangential induction factor at each node
"Alpha" - Angle of attack at each node
"Theta" - Pitch+Twist angle at each node
"Phi" - Inflow angle at each node
"Curve" - Curvature angle at each node
"Cl" - Lift force coefficient at each node
"Cd" - Drag force coefficient at each node
"Cm" - Pitching moment coefficient at each node
"Cx" - Normal force (to plane) coefficient at each node
"Cy" - Tangential force (to plane) coefficient at each node
"Cn" - Normal force (to chord) coefficient at each node
"Ct" - Tangential force (to chord) coefficient at each node
"Fl" - Lift force per unit length at each node
"Fd" - Drag force per unit length at each node
"Mm" - Pitching moment per unit length at each node
"Fx" - Normal force (to plane) per unit length at each node
"Fy" - Tangential force (to plane) per unit length at each node
"Fn" - Normal force (to chord) per unit length at each node
"Ft" - Tangential force (to chord) per unit length at each node
"Clrnc" - Tower clearance at each node (based on the absolute distance to the nearest point in the tower from blade node B#N# minus the local tower radius, in the deflected configuration); please note that this clearance is only approximate because the calculation assumes that the blade is a line with no volume (however, the calculation does use the local tower radius); when blade node B#N# is above the tower top (or below the tower base), the absolute distance to the tower top (or base) minus the local tower radius, in the deflected configuration, is output
"Vx" - Local axial velocity
"Vy" - Local tangential velocity
"GeomPhi" - Geometric phi? If phi was solved using normal BEMT equations, GeomPhi = 1; otherwise, if it was solved geometrically, GeomPhi = 0.
"Chi" - Skew angle (used in skewed wake correction) -- not available for OLAF
"UA_Flag" - Flag indicating if UA is turned on for this node. -- not available for OLAF
"CpMin" - Pressure coefficient
"SgCav" - Cavitation number
"SigCr" - Critical cavitation number
"Gam" - Gamma -- circulation on blade
"Cl_Static" - Static portion of lift force coefficient at each node, without unsteady effects -- not available for BEMT/DBEMT
"Cd_Static" - Static portion of drag force coefficient at each node, without unsteady effects -- not available for BEMT/DBEMT
"Cm_Static" - Static portion of pitching moment coefficient at each node, without unsteady effects -- not available for BEMT/DBEMT
"Uin" - Axial induced velocity in rotating hub coordinates. Axial aligned with hub axis. rotor plane polar hub rotating coordinates
"Uit" - Tangential induced velocity in rotating hub coordinates. Tangential to the rotation plane. Perpendicular to blade aziumth. rotor plane polar hub rotating coordinates
"Uir" - Radial induced velocity in rotating hub coordinates. Radial outwards in rotation plane. Aligned with blade azimuth. rotor plane polar hub rotating coordinates
END of input file (the word "END" must appear in the first 3 columns of this last OutList line)
---------------------------------------------------------------------------------------


#### The following was for testing, but bloats the regression test case too much
#---------------------- NODE OUTPUTS --------------------------------------------
# 3 BldNd_BladesOut - Blades to output
# 99 BldNd_BlOutNd - Blade nodes on each blade (currently unused)
# OutList - The next line(s) contains a list of output parameters. See OutListParameters.xlsx, AeroDyn_Nodes tab for a listing of available output channels, (-)
#"VUndx" - x-component of undisturbed wind velocity at each node
#"VUndy" - y-component of undisturbed wind velocity at each node
#"VUndz" - z-component of undisturbed wind velocity at each node
#"VDisx" - x-component of disturbed wind velocity at each node
#"VDisy" - y-component of disturbed wind velocity at each node
#"VDisz" - z-component of disturbed wind velocity at each node
#"STVx" - x-component of structural translational velocity at each node
#"STVy" - y-component of structural translational velocity at each node
#"STVz" - z-component of structural translational velocity at each node
#"VRel" - Relvative wind speed at each node
#"DynP" - Dynamic pressure at each node
#"Re" - Reynolds number (in millions) at each node
#"M" - Mach number at each node
#"Vindx" - Axial induced wind velocity at each node
#"Vindy" - Tangential induced wind velocity at each node
#"AxInd" - Axial induction factor at each node
#"TnInd" - Tangential induction factor at each node
#"Alpha" - Angle of attack at each node
#"Theta" - Pitch+Twist angle at each node
#"Phi" - Inflow angle at each node
#"Curve" - Curvature angle at each node
#"Cl" - Lift force coefficient at each node
#"Cd" - Drag force coefficient at each node
#"Cm" - Pitching moment coefficient at each node
#"Cx" - Normal force (to plane) coefficient at each node
#"Cy" - Tangential force (to plane) coefficient at each node
#"Cn" - Normal force (to chord) coefficient at each node
#"Ct" - Tangential force (to chord) coefficient at each node
#"Fl" - Lift force per unit length at each node
#"Fd" - Drag force per unit length at each node
#"Mm" - Pitching moment per unit length at each node
#"Fx" - Normal force (to plane) per unit length at each node
#"Fy" - Tangential force (to plane) per unit length at each node
#"Fn" - Normal force (to chord) per unit length at each node
#"Ft" - Tangential force (to chord) per unit length at each node
#"Clrnc" - Tower clearance at each node (based on the absolute distance to the nearest point in the tower from blade node B#N# minus the local tower radius, in the deflected configuration); please note that this clearance is only approximate because the calculation assumes that the blade is a line with no volume (however, the calculation does use the local tower radius); when blade node B#N# is above the tower top (or below the tower base), the absolute distance to the tower top (or base) minus the local tower radius, in the deflected configuration, is output
#"Vx" - Local axial velocity
#"Vy" - Local tangential velocity
#"GeomPhi" - Geometric phi? If phi was solved using normal BEMT equations, GeomPhi = 1; otherwise, if it was solved geometrically, GeomPhi = 0.
#"Chi" - Skew angle (used in skewed wake correction) -- not available for OLAF
#"UA_Flag" - Flag indicating if UA is turned on for this node. -- not available for OLAF
#"CpMin" - Pressure coefficient
#"SgCav" - Cavitation number
#"SigCr" - Critical cavitation number
#"Gam" - Gamma -- circulation on blade
#"Cl_Static" - Static portion of lift force coefficient at each node, without unsteady effects -- not available for BEMT/DBEMT
#"Cd_Static" - Static portion of drag force coefficient at each node, without unsteady effects -- not available for BEMT/DBEMT
#"Cm_Static" - Static portion of pitching moment coefficient at each node, without unsteady effects -- not available for BEMT/DBEMT
#"Uin" - Axial induced velocity in rotating hub coordinates. Axial aligned with hub axis. rotor plane polar hub rotating coordinates
#"Uit" - Tangential induced velocity in rotating hub coordinates. Tangential to the rotation plane. Perpendicular to blade aziumth. rotor plane polar hub rotating coordinates
#"Uir" - Radial induced velocity in rotating hub coordinates. Radial outwards in rotation plane. Aligned with blade azimuth. rotor plane polar hub rotating coordinates
#END of input file (the word "END" must appear in the first 3 columns of this last OutList line)
#---------------------------------------------------------------------------------------
16 changes: 8 additions & 8 deletions modules/aerodyn/py_ad_5MW_OC4Semi_WSt_WavesWN/py_ad_driver.out
View file

Large diffs are not rendered by default.

4 changes: 2 additions & 2 deletions ...y_ad_5MW_OC4Semi_WSt_WavesWN/vtkRef/5MW_OC4Semi_WSt_WavesWN.AD_BladeRootMotion1.00000.vtp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -16,10 +16,10 @@
0.93565 -0.00000 -0.11321
</DataArray>
<DataArray type="Float32" Name="TranslationalAcceleration" NumberOfComponents="3" format="ascii">
-0.09341 -0.01365 -1.60564
-0.08747 -0.01414 -1.57571
</DataArray>
<DataArray type="Float32" Name="RotationalAcceleration" NumberOfComponents="3" format="ascii">
0.03947 -0.06181 -0.01507
0.03877 -0.05620 -0.01474
</DataArray>
<DataArray type="Float32" Name="OrientationX" NumberOfComponents="3" format="ascii">
0.99705 0.00000 -0.07671
Expand Down
8 changes: 4 additions & 4 deletions ...y_ad_5MW_OC4Semi_WSt_WavesWN/vtkRef/5MW_OC4Semi_WSt_WavesWN.AD_BladeRootMotion1.00001.vtp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -10,16 +10,16 @@
</Points>
<PointData>
<DataArray type="Float32" Name="TranslationalVelocity" NumberOfComponents="3" format="ascii">
-0.00111 -1.41247 -0.02000
-0.00104 -1.41248 -0.01963
</DataArray>
<DataArray type="Float32" Name="RotationalVelocity" NumberOfComponents="3" format="ascii">
0.93608 -0.00076 -0.11339
0.93607 -0.00069 -0.11338
</DataArray>
<DataArray type="Float32" Name="TranslationalAcceleration" NumberOfComponents="3" format="ascii">
-0.07992 -0.01449 -1.58931
-0.07520 -0.01555 -1.56336
</DataArray>
<DataArray type="Float32" Name="RotationalAcceleration" NumberOfComponents="3" format="ascii">
0.03099 -0.05816 -0.01094
0.03032 -0.05286 -0.01059
</DataArray>
<DataArray type="Float32" Name="OrientationX" NumberOfComponents="3" format="ascii">
0.99705 -0.00051 -0.07671
Expand Down
10 changes: 5 additions & 5 deletions ...y_ad_5MW_OC4Semi_WSt_WavesWN/vtkRef/5MW_OC4Semi_WSt_WavesWN.AD_BladeRootMotion1.00002.vtp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,21 +5,21 @@
NumberOfStrips=" 0" NumberOfPolys=" 0">
<Points>
<DataArray type="Float32" NumberOfComponents="3" format="ascii">
3.10065 -0.03531 91.61136
3.10066 -0.03531 91.61137
</DataArray>
</Points>
<PointData>
<DataArray type="Float32" Name="TranslationalVelocity" NumberOfComponents="3" format="ascii">
-0.00190 -1.41268 -0.03960
-0.00179 -1.41270 -0.03896
</DataArray>
<DataArray type="Float32" Name="RotationalVelocity" NumberOfComponents="3" format="ascii">
0.93638 -0.00143 -0.11349
0.93636 -0.00130 -0.11348
</DataArray>
<DataArray type="Float32" Name="TranslationalAcceleration" NumberOfComponents="3" format="ascii">
-0.04551 -0.00079 -1.54326
-0.03993 -0.01322 -1.52015
</DataArray>
<DataArray type="Float32" Name="RotationalAcceleration" NumberOfComponents="3" format="ascii">
0.00561 -0.04859 -0.00429
0.01009 -0.04383 -0.00390
</DataArray>
<DataArray type="Float32" Name="OrientationX" NumberOfComponents="3" format="ascii">
0.99705 -0.00103 -0.07670
Expand Down
14 changes: 7 additions & 7 deletions ...y_ad_5MW_OC4Semi_WSt_WavesWN/vtkRef/5MW_OC4Semi_WSt_WavesWN.AD_BladeRootMotion1.00003.vtp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,30 +5,30 @@
NumberOfStrips=" 0" NumberOfPolys=" 0">
<Points>
<DataArray type="Float32" NumberOfComponents="3" format="ascii">
3.10063 -0.05297 91.61074
3.10063 -0.05297 91.61076
</DataArray>
</Points>
<PointData>
<DataArray type="Float32" Name="TranslationalVelocity" NumberOfComponents="3" format="ascii">
-0.00216 -1.41253 -0.05846
-0.00196 -1.41283 -0.05755
</DataArray>
<DataArray type="Float32" Name="RotationalVelocity" NumberOfComponents="3" format="ascii">
0.93622 -0.00195 -0.11349
0.93631 -0.00176 -0.11347
</DataArray>
<DataArray type="Float32" Name="TranslationalAcceleration" NumberOfComponents="3" format="ascii">
0.01024 0.00964 -1.46960
0.01511 -0.00444 -1.45176
</DataArray>
<DataArray type="Float32" Name="RotationalAcceleration" NumberOfComponents="3" format="ascii">
-0.02478 -0.03379 0.00416
-0.01893 -0.02998 0.00418
</DataArray>
<DataArray type="Float32" Name="OrientationX" NumberOfComponents="3" format="ascii">
0.99705 -0.00155 -0.07669
0.99705 -0.00155 -0.07670
</DataArray>
<DataArray type="Float32" Name="OrientationY" NumberOfComponents="3" format="ascii">
0.00425 0.99937 0.03510
</DataArray>
<DataArray type="Float32" Name="OrientationZ" NumberOfComponents="3" format="ascii">
0.07659 -0.03532 0.99644
0.07660 -0.03532 0.99644
</DataArray>
</PointData>
</Piece>
Expand Down
Loading

0 comments on commit cddc9fd

Please sign in to comment.