Learning Forum

[ashwin Gurjar]

Hello 

I am using Ansys aqwa for hydrodynamic analysis of offhore wind turbines and I want to couple it with another tool like OpenFast for wind load calculation.

or

I want to add externally calcualted force in to the simulation via user_force.dll but I have no idea how to do it. How dll files will as a input will work.

Can you please explain this to me?

Thank You

[Yu Bi]

Hi, from 2022R2, Aqwa-Rigid Dynamics-AeroDyn co-simulation system for floating wind turbines is available. This workflow needs Ansys Twin Builder and Ansys Mechanical. More details about it can be found in Aqwa User's Manual Section 8.3.  

Regarding external force, you can find a python example in the installation folder under aqwa/utils/ExternalForceCalculation while the FORTRAN and C examples can be found under aqwa/doc. If Fortran or C language is used, then compile your own routine and then put the *.dll file (Windows) into the same directory as the Aqwa executables. More information can be found in Aqwa Reference Manual Chapter 30.

Hope this helps.

[mohsen ab.1999]

hello .....i want to simulate floating offshore wind turbine with the new add on pakage (co-simulation) between ansys aqwa and rigid dynamic and AeroDyn. in wizard option when i upload primary input file and inflow wind input file .... the software give me erorre which it cannot read the file ..... the file which i upload to this has the (.text) format ..... and it the example file which i download from AeroDyn user guid .... can you help me.... what is the problem ?

[Yu Bi]

It could be due to the AeroDyn file version is not what the wizard designed for. Please use the file format from this manual: https://raf-openfast.readthedocs.io/en/docs-legacy/source/user/index.html

[mohsen ab.1999]

thank you for fast answering ....i read the file format manual but and upload (blade data input file + Airfoil data input file) for AeroDyn primary input file and upload (inflow wind input file ) for inflow wind input file in co simulation wizard and again it does not work and this massage appear in extension log file:AqwaCosimulation: Failed to find the AeroDyn blade data input file defined in the AeroDyn primary input file (C:\Users\fatemeh.km\Downloads\ad_primary_example.dat)
AqwaCosimulation: Failed to prepare the AeroDyn marker data.

can you send me primery input file and inflow wind input file sample which can be upload in co simulation wizared

that's my gmail: abnmohaen@gmail.com

[Yu Bi]

Hi, the blade data input file and Airfoil data input file should follow the directory and file name written in the AeroDyn primary input file, and the same for the inflow wind file. So please check if those files are in the correct directory and with the correct name. For example in the ad_primary_example.dat file downloaded from the AeroDyn manual I sent, the blade data input file should be in the same directory as the primary file and named Test01_UAE_AeroDyn_blade.dat.

[mohsen ab.1999]

i checked all over the points which you mentioned in above massages for me but this give me erorre yet like before :

AqwaCosimulation: Inadequate blade data files (ADBlFile) given in the AeroDyn primary input file (E:\hydraulic pitch system of wind turbin\new reasearch\open fast\r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338\r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338\modules\aerodyn\py_ad_5MW_OC4Semi_WSt_WavesWN\ad_primary.dat)
AqwaCosimulation: Failed to prepare the AeroDyn marker data.

now i want to shear with you me Aero Dyn primary input file which actualle relate to NREL 5 mw and please help me to correct this because there is no toturial exist about that yet ..

this is my Aero Dyn primary input file :

------- AERODYN v15 for OpenFAST INPUT FILE -----------------------------------------------
NREL 5.0 MW offshore baseline aerodynamic input properties.
======  General Options  ============================================================================
True          Echo               - Echo the input to ".AD.ech"?  (flag)
          1    DTAero             - Time interval for aerodynamic calculations {or "default"} (s)
          1   WakeMod            - Type of wake/induction model (switch) {0=none, 1=BEMT, 2=DBEMT, 3=OLAF} [WakeMod cannot be 2 or 3 when linearizing]
          1   AFAeroMod          - Type of blade airfoil aerodynamics model (switch) {1=steady model, 2=Beddoes-Leishman unsteady model} [AFAeroMod must be 1 when linearizing]
          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}
          0   TwrShadow          - Calculate tower influence on wind based on downstream tower shadow (switch) {0=none, 1=Powles model, 2=Eames model}
True          TwrAero            - Calculate tower aerodynamic loads? (flag)
False         FrozenWake         - Assume frozen wake during linearization? (flag) [used only when WakeMod=1 and when linearizing]
False         CavitCheck         - Perform cavitation check? (flag) [AFAeroMod must be 1 when CavitCheck=true]
True         Buoyancy           - Include buoyancy effects? (flag)
False         CompAA             - Flag to compute AeroAcoustics calculation [used only when WakeMod = 1 or 2]
"unused"      AA_InputFile       - AeroAcoustics input file [used only when CompAA=true]
======  Environmental Conditions  ===================================================================
      1.225   AirDens            - Air density (kg/m^3)
    0.00001   KinVisc            - Kinematic viscosity of working fluid (m^2/s)
      340.3   SpdSound           - Speed of sound in working fluid (m/s)
     101325   Patm               - Atmospheric pressure (Pa) [used only when CavitCheck=True]
       2000   Pvap               - Vapour pressure of working fluid (Pa) [used only when CavitCheck=True]
======  Blade-Element/Momentum Theory Options  ====================================================== [unused when WakeMod=0 or 3]
          1   SkewMod            - Type of skewed-wake correction model (switch) {1=uncoupled, 2=Pitt/Peters, 3=coupled} [unused when WakeMod=0 or 3]
          2    SkewModFactor      - Constant used in Pitt/Peters skewed wake model {or "default" is 15/32*pi} (-) [used only when SkewMod=2; unused when WakeMod=0 or 3]
True          TipLoss            - Use the Prandtl tip-loss model? (flag) [unused when WakeMod=0 or 3]
True          HubLoss            - Use the Prandtl hub-loss model? (flag) [unused when WakeMod=0 or 3]
True          TanInd             - Include tangential induction in BEMT calculations? (flag) [unused when WakeMod=0 or 3]
False         AIDrag             - Include the drag term in the axial-induction calculation? (flag) [unused when WakeMod=0 or 3]
False         TIDrag             - Include the drag term in the tangential-induction calculation? (flag) [unused when WakeMod=0,3 or TanInd=FALSE]
    0.00001   IndToler           - Convergence tolerance for BEMT nonlinear solve residual equation {or "default"} (-) [unused when WakeMod=0 or 3]
        100   MaxIter            - Maximum number of iteration steps (-) [unused when WakeMod=0]
======  Dynamic Blade-Element/Momentum Theory Options  ============================================== [used only when WakeMod=2]
          2   DBEMT_Mod          - Type of dynamic BEMT (DBEMT) model {1=constant tau1, 2=time-dependent tau1, 3=constant tau1 with continuous formulation} (-) [used only when WakeMod=2]
          4   tau1_const         - Time constant for DBEMT (s) [used only when WakeMod=2 and DBEMT_Mod=1 or 3] 
======  OLAF -- cOnvecting LAgrangian Filaments (Free Vortex Wake) Theory Options  ================== [used only when WakeMod=3]
"unused"      OLAFInputFileName  - Input file for OLAF [used only when WakeMod=3]
======  Beddoes-Leishman Unsteady Airfoil Aerodynamics Options  ===================================== [used only when AFAeroMod=2]
          3   UAMod              - Unsteady Aero Model Switch (switch) {2=B-L Gonzalez, 3=B-L Minnema/Pierce, 4=B-L HGM 4-states, 5=B-L 5 states, 6=Oye, 7=Boeing-Vertol} [used only when AFAeroMod=2]
True          FLookup            - Flag to indicate whether a lookup for f' will be calculated (TRUE) or whether best-fit exponential equations will be used (FALSE); if FALSE S1-S4 must be provided in airfoil input files (flag) [used only when AFAeroMod=2]
======  Airfoil Information =========================================================================
          1   AFTabMod           - Interpolation method for multiple airfoil tables {1=1D interpolation on AoA (first table only); 2=2D interpolation on AoA and Re; 3=2D interpolation on AoA and UserProp} (-)
          1   InCol_Alfa         - The column in the airfoil tables that contains the angle of attack (-)
          2   InCol_Cl           - The column in the airfoil tables that contains the lift coefficient (-)
          3   InCol_Cd           - The column in the airfoil tables that contains the drag coefficient (-)
          4   InCol_Cm           - The column in the airfoil tables that contains the pitching-moment coefficient; use zero if there is no Cm column (-)
          0   InCol_Cpmin        - The column in the airfoil tables that contains the Cpmin coefficient; use zero if there is no Cpmin column (-)
          8   NumAFfiles         - Number of airfoil files used (-)
          1    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/Cylinder1.dat"    AFNames            - Airfoil file names (NumAFfiles lines) (quoted strings)
          2    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/Cylinder2.dat"
          3    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/DU40_A17.dat"
          4    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/DU35_A17.dat"
          5    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/DU30_A17.dat"
          6    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/DU25_A17.dat"
          7    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/DU21_A17.dat"
          8    "E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/Airfoils/NACA64_A17.dat"
======  Rotor/Blade Properties  =====================================================================
True          UseBlCm            - Include aerodynamic pitching moment in calculations?  (flag)
"E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/NRELOffshrBsline5MW_AeroDyn_blade1.dat"    ADBlFile(1)        - Name of file containing distributed aerodynamic properties for Blade #1 (-)
"E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/NRELOffshrBsline5MW_AeroDyn_blade2.dat"    ADBlFile(2)        - Name of file containing distributed aerodynamic properties for Blade #2 (-) [unused if NumBl < 2]
"E:/hydraulic pitch system of wind turbin/new reasearch/open fast/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/r-test-ef6ab575919f1db597702f1c3f86a16bdc2b4338/glue-codes/openfast/5MW_Baseline/NRELOffshrBsline5MW_AeroDyn_blade3.dat"    ADBlFile(3)        - Name of file containing distributed aerodynamic properties for Blade #3 (-) [unused if NumBl < 3]
======  Hub Properties ============================================================================== [used only when Buoyancy=True]
22    VolHub             - Hub volume (m^3)
6.0   HubCenBx           - Hub center of buoyancy x direction offset (m)
======  Nacelle Properties ========================================================================== [used only when Buoyancy=True]
40.0  VolNac             - Nacelle volume (m^3)
1,2,3 NacCenB            - Position of nacelle center of buoyancy from yaw bearing in nacelle coordinates (m)
======  Tail fin Aerodynamics ======================================================================== 
False         TFinAero           - Calculate tail fin aerodynamics model (flag)
"unused"      TFinFile           - Input file for tail fin aerodynamics [used only when TFinAero=True]
======  Tower Influence and Aerodynamics ============================================================ [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True]
         11   NumTwrNds         - Number of tower nodes used in the analysis  (-) [used only when TwrPotent/=0, TwrShadow/=0, TwrAero=True, or Buoyancy=True]
TwrElev        TwrDiam        TwrCd          TwrTI          TwrCb !TwrTI used only with TwrShadow=2, TwrCb used only with Buoyancy=True
(m)              (m)           (-)            (-)           (-)
1.0000000E+01  6.5000000E+00  1.0000000E+00  1.0000000E-01  0.0
1.7760000E+01  6.2400000E+00  1.0000000E+00  1.0000000E-01  0.0
2.5520000E+01  5.9700000E+00  1.0000000E+00  1.0000000E-01  0.0
3.3280000E+01  5.7100000E+00  1.0000000E+00  1.0000000E-01  0.0
4.1040000E+01  5.4500000E+00  1.0000000E+00  1.0000000E-01  0.0
4.8800000E+01  5.1800000E+00  1.0000000E+00  1.0000000E-01  0.0
5.6560000E+01  4.9200000E+00  1.0000000E+00  1.0000000E-01  0.0
6.4320000E+01  4.6600000E+00  1.0000000E+00  1.0000000E-01  0.0
7.2080000E+01  4.4000000E+00  1.0000000E+00  1.0000000E-01  0.0
7.9840000E+01  4.1300000E+00  1.0000000E+00  1.0000000E-01  0.0
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 ".AD.sum"?  (flag)
          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)

#### 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)
#---------------------------------------------------------------------------------------

 

 

and in the airfoil input file we should call the airfoil cordinate file with file name or file path?

 

[mohsen ab.1999]

.

 

[Yu Bi]

Hi, I think it is because your blade data file path has spaces in it. Please exclude the spaces in the path of files given in the AeroDyn primary file and input wind file. You can use the relative path instead of the full path.