Current and Wave simulation

Hi everybody,

I was performing OC4 load case 2.4 (current and regular waves) and a similar one but with the wave and current in the opposite direction with the semisubmersible. The idea was to evaluate if there is some kind of symmetrical behavior between the two. The only difference between the hydrodyn files is, thus, CurrSSDir and WaveDir. In one case they are both 0, in the other, they are both 180. The issue is that the results seem to be odd. I’ve attached some of the results below.

What do you think of them? I would say, for instance that that surge in the 180º case would be negative and the sway null… I don’t understand if it is correct and unexpected or if there’s some modelling mistake.

Best regards,

Tomás
sway.JPG
roll.JPG
surge.JPG

Dear Tomas,

I agree that these results really don’t make sense to me, but I’m not sure what is set incorrectly. From the surge response, it looks like the 180-degree waves and current, causes a reduction in mean surge, but there must still be some force pushing the semi in positive surge. Do you have aerodynamic loads enabled? And I would not expect any sway and roll motion in a simulation with waves and current only.

Best regards,

Dear Jason,

Thanks for your reply. I do not have any wind in the simulation: CompInFlow=0.
In order to try to understand better what’s going on, I’ve run other two cases splitting the previous one: only current and only waves, both with CurrSSDir and WaveDirr at 180º. The results are below. Somehow it makes more sense (at least as current is concerned), though the wave effect is still not expected. But maybe that’s the way the platform reacts to the waves…
sway.JPG
surge.JPG
roll.JPG

Dear Tomas,

Are you using a floating platform model NREL has provided as an example (e.g., the OC4-DeepCwind semisubmersible) or a model you made yourself? Do you have second-order difference-frequency effects enabled (mean drift), which will induce a mean load for regular waves? I’d probably need to know more about your model set up to comment more.

Best regards,

Dear Jason,

Yes, I’m using OC4-DeepCwind semisubmersible and without second-order difference-frequency effects enabled.

Here’s the hydrodyn file used:

------- HydroDyn v2.03.* Input File --------------------------------------------
NREL 5.0 MW offshore baseline floating platform HydroDyn input properties for the OC4 Semi-submersible.
False Echo - Echo the input file data (flag)
---------------------- ENVIRONMENTAL CONDITIONS --------------------------------
1025 WtrDens - Water density (kg/m^3)
200 WtrDpth - Water depth (meters)
0 MSL2SWL - Offset between still-water level and mean sea level (meters) [positive upward; unused when WaveMod = 6; must be zero if PotMod=1 or 2]
---------------------- WAVES ---------------------------------------------------
1 WaveMod - Incident wave kinematics model {0: none=still water, 1: regular (periodic), 1P#: regular with user-specified phase, 2: JONSWAP/Pierson-Moskowitz spectrum (irregular), 3: White noise spectrum (irregular), 4: user-defined spectrum from routine UserWaveSpctrm (irregular), 5: Externally generated wave-elevation time series, 6: Externally generated full wave-kinematics time series [option 6 is invalid for PotMod/=0]} (switch)
0 WaveStMod - Model for stretching incident wave kinematics to instantaneous free surface {0: none=no stretching, 1: vertical stretching, 2: extrapolation stretching, 3: Wheeler stretching} (switch) [unused when WaveMod=0 or when PotMod/=0]
4600 WaveTMax - Analysis time for incident wave calculations (sec) [unused when WaveMod=0; determines WaveDOmega=2Pi/WaveTMax in the IFFT]
0.2 WaveDT - Time step for incident wave calculations (sec) [unused when WaveMod=0; 0.1<=WaveDT<=1.0 recommended; determines WaveOmegaMax=Pi/WaveDT in the IFFT]
6 WaveHs - Significant wave height of incident waves (meters) [used only when WaveMod=1, 2, or 3]
10 WaveTp - Peak-spectral period of incident waves (sec) [used only when WaveMod=1 or 2]
1.05 WavePkShp - Peak-shape parameter of incident wave spectrum (-) or DEFAULT (string) [used only when WaveMod=2; use 1.0 for Pierson-Moskowitz]
0.314159 WvLowCOff - Low cut-off frequency or lower frequency limit of the wave spectrum beyond which the wave spectrum is zeroed (rad/s) [unused when WaveMod=0, 1, or 6]
1.570796 WvHiCOff - High cut-off frequency or upper frequency limit of the wave spectrum beyond which the wave spectrum is zeroed (rad/s) [unused when WaveMod=0, 1, or 6]
180 WaveDir - Incident wave propagation heading direction (degrees) [unused when WaveMod=0 or 6]
0 WaveDirMod - Directional spreading function {0: none, 1: COS2S} (-) [only used when WaveMod=2,3, or 4]
1 WaveDirSpread - Wave direction spreading coefficient ( > 0 ) (-) [only used when WaveMod=2,3, or 4 and WaveDirMod=1]
1 WaveNDir - Number of wave directions (-) [only used when WaveMod=2,3, or 4 and WaveDirMod=1; odd number only]
0 WaveDirRange - Range of wave directions (full range: WaveDir +/- 1/2*WaveDirRange) (degrees) [only used when WaveMod=2,3,or 4 and WaveDirMod=1]
123456789 WaveSeed(1) - First random seed of incident waves [-2147483648 to 2147483647] (-) [unused when WaveMod=0, 5, or 6]
1011121314 WaveSeed(2) - Second random seed of incident waves [-2147483648 to 2147483647] (-) [unused when WaveMod=0, 5, or 6]
FALSE WaveNDAmp - Flag for normally distributed amplitudes (flag) [only used when WaveMod=2, 3, or 4]
“” WvKinFile - Root name of externally generated wave data file(s) (quoted string) [used only when WaveMod=5 or 6]
1 NWaveElev - Number of points where the incident wave elevations can be computed (-) [maximum of 9 output locations]
0 WaveElevxi - List of xi-coordinates for points where the incident wave elevations can be output (meters) [NWaveElev points, separated by commas or white space; usused if NWaveElev = 0]
0 WaveElevyi - List of yi-coordinates for points where the incident wave elevations can be output (meters) [NWaveElev points, separated by commas or white space; usused if NWaveElev = 0]
---------------------- 2ND-ORDER WAVES ----------------------------------------- [unused with WaveMod=0 or 6]
FALSE WvDiffQTF - Full difference-frequency 2nd-order wave kinematics (flag)
FALSE WvSumQTF - Full summation-frequency 2nd-order wave kinematics (flag)
0 WvLowCOffD - Low frequency cutoff used in the difference-frequencies (rad/s) [Only used with a difference-frequency method]
1.256637 WvHiCOffD - High frequency cutoff used in the difference-frequencies (rad/s) [Only used with a difference-frequency method]
0.618319 WvLowCOffS - Low frequency cutoff used in the summation-frequencies (rad/s) [Only used with a summation-frequency method]
3.141593 WvHiCOffS - High frequency cutoff used in the summation-frequencies (rad/s) [Only used with a summation-frequency method]
---------------------- CURRENT ------------------------------------------------- [unused with WaveMod=6]
0 CurrMod - Current profile model {0: none=no current, 1: standard, 2: user-defined from routine UserCurrent} (switch)
0.5 CurrSSV0 - Sub-surface current velocity at still water level (m/s) [used only when CurrMod=1]
30 CurrSSDir - Sub-surface current heading direction (degrees) or DEFAULT (string) [used only when CurrMod=1]
20 CurrNSRef - Near-surface current reference depth (meters) [used only when CurrMod=1]
0 CurrNSV0 - Near-surface current velocity at still water level (m/s) [used only when CurrMod=1]
0 CurrNSDir - Near-surface current heading direction (degrees) [used only when CurrMod=1]
0 CurrDIV - Depth-independent current velocity (m/s) [used only when CurrMod=1]
0 CurrDIDir - Depth-independent current heading direction (degrees) [used only when CurrMod=1]
---------------------- FLOATING PLATFORM --------------------------------------- [unused with WaveMod=6]
1 PotMod - Potential-flow model {0: none=no potential flow, 1: frequency-to-time-domain transforms based on WAMIT output, 2: fluid-impulse theory (FIT)} (switch)
“…/5MW_Baseline/HydroData/marin_semi” PotFile - Root name of potential-flow model data; WAMIT output files containing the linear, nondimensionalized, hydrostatic restoring matrix (.hst), frequency-dependent hydrodynamic added mass matrix and damping matrix (.1), and frequency- and direction-dependent wave excitation force vector per unit wave amplitude (.3) (quoted string) [MAKE SURE THE FREQUENCIES INHERENT IN THESE WAMIT FILES SPAN THE PHYSICALLY-SIGNIFICANT RANGE OF FREQUENCIES FOR THE GIVEN PLATFORM; THEY MUST CONTAIN THE ZERO- AND INFINITE-FREQUENCY LIMITS!]
1 WAMITULEN - Characteristic body length scale used to redimensionalize WAMIT output (meters) [only used when PotMod=1]
13917 PtfmVol0 - Displaced volume of water when the platform is in its undisplaced position (m^3) [only used when PotMod=1; USE THE SAME VALUE COMPUTED BY WAMIT AS OUTPUT IN THE .OUT FILE!]
0 PtfmCOBxt - The xt offset of the center of buoyancy (COB) from the platform reference point (meters) [only used when PotMod=1]
0 PtfmCOByt - The yt offset of the center of buoyancy (COB) from the platform reference point (meters) [only used when PotMod=1]
1 RdtnMod - Radiation memory-effect model {0: no memory-effect calculation, 1: convolution, 2: state-space} (switch) [only used when PotMod=1; STATE-SPACE REQUIRES *.ss INPUT FILE]
60 RdtnTMax - Analysis time for wave radiation kernel calculations (sec) [only used when PotMod=1; determines RdtnDOmega=Pi/RdtnTMax in the cosine transform; MAKE SURE THIS IS LONG ENOUGH FOR THE RADIATION IMPULSE RESPONSE FUNCTIONS TO DECAY TO NEAR-ZERO FOR THE GIVEN PLATFORM!]
0.0125 RdtnDT - Time step for wave radiation kernel calculations (sec) [only used when PotMod=1; DT<=RdtnDT<=0.1 recommended; determines RdtnOmegaMax=Pi/RdtnDT in the cosine transform]
---------------------- 2ND-ORDER FLOATING PLATFORM FORCES ---------------------- [unused with WaveMod=0 or 6, or PotMod=0 or 2]
0 MnDrift - Mean-drift 2nd-order forces computed {0: None; [7, 8, 9, 10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero]
0 NewmanApp - Mean- and slow-drift 2nd-order forces computed with Newman’s approximation {0: None; [7, 8, 9, 10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero. Used only when WaveDirMod=0]
12 DiffQTF - Full difference-frequency 2nd-order forces computed with full QTF {0: None; [10, 11, or 12]: WAMIT file to use} [Only one of MnDrift, NewmanApp, or DiffQTF can be non-zero]
12 SumQTF - Full summation -frequency 2nd-order forces computed with full QTF {0: None; [10, 11, or 12]: WAMIT file to use}
---------------------- FLOATING PLATFORM FORCE FLAGS -------------------------- [unused with WaveMod=6]
True PtfmSgF - Platform horizontal surge translation force (flag) or DEFAULT
True PtfmSwF - Platform horizontal sway translation force (flag) or DEFAULT
True PtfmHvF - Platform vertical heave translation force (flag) or DEFAULT
True PtfmRF - Platform roll tilt rotation force (flag) or DEFAULT
True PtfmPF - Platform pitch tilt rotation force (flag) or DEFAULT
True PtfmYF - Platform yaw rotation force (flag) or DEFAULT
---------------------- PLATFORM ADDITIONAL STIFFNESS AND DAMPING --------------
0 0 0 0 0 0 AddF0 - Additional preload (N, N-m)
0 0 0 0 0 0 AddCLin - Additional linear stiffness (N/m, N/rad, N-m/m, N-m/rad)
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1451298897 0 0
0 0 0 0 1451298897 0
0 0 0 0 0 0
0 0 0 0 0 0 AddBLin - Additional linear damping(N/(m/s), N/(rad/s), N-m/(m/s), N-m/(rad/s))
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0 AddBQuad - Additional quadratic drag(N/(m/s)^2, N/(rad/s)^2, N-m(m/s)^2, N-m/(rad/s)^2)
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
---------------------- AXIAL COEFFICIENTS --------------------------------------
2 NAxCoef - Number of axial coefficients (-)
AxCoefID AxCd AxCa AxCp
(-) (-) (-) (-)
1 0.00 0.00 1.00
2 9.60 0.00 1.00
---------------------- MEMBER JOINTS -------------------------------------------
44 NJoints - Number of joints (-) [must be exactly 0 or at least 2]
JointID Jointxi Jointyi Jointzi JointAxID JointOvrlp [JointOvrlp= 0: do nothing at joint, 1: eliminate overlaps by calculating super member]
(-) (m) (m) (m) (-) (switch)
1 0.00000 0.00000 -20.00000 1 0
2 0.00000 0.00000 10.00000 1 0
3 14.43376 25.00000 -14.00000 1 0
4 14.43376 25.00000 12.00000 1 0
5 -28.86751 0.00000 -14.00000 1 0
6 -28.86751 0.00000 12.00000 1 0
7 14.43376 -25.00000 -14.00000 1 0
8 14.43376 -25.00000 12.00000 1 0
9 14.43375 25.00000 -20.00000 2 0
10 -28.86750 0.00000 -20.00000 2 0
11 14.43375 -25.00000 -20.00000 2 0
12 9.23760 22.00000 10.00000 1 0
13 -23.67130 3.00000 10.00000 1 0
14 -23.67130 -3.00000 10.00000 1 0
15 9.23760 -22.00000 10.00000 1 0
16 14.43375 -19.00000 10.00000 1 0
17 14.43375 19.00000 10.00000 1 0
18 4.04145 19.00000 -17.00000 1 0
19 -18.47520 6.00000 -17.00000 1 0
20 -18.47520 -6.00000 -17.00000 1 0
21 4.04145 -19.00000 -17.00000 1 0
22 14.43375 -13.00000 -17.00000 1 0
23 14.43375 13.00000 -17.00000 1 0
24 1.62500 2.81500 10.00000 1 0
25 11.43376 19.80385 10.00000 1 0
26 -3.25000 0.00000 10.00000 1 0
27 -22.87000 0.00000 10.00000 1 0
28 1.62500 -2.81500 10.00000 1 0
29 11.43376 -19.80385 10.00000 1 0
30 1.62500 2.81500 -17.00000 1 0
31 8.43376 14.60770 -17.00000 1 0
32 -3.25000 0.00000 -17.00000 1 0
33 -16.87000 0.00000 -17.00000 1 0
34 1.62500 -2.81500 -17.00000 1 0
35 8.43376 -14.60770 -17.00000 1 0
36 1.62500 2.81500 -16.20000 1 0
37 11.43376 19.80385 9.13000 1 0
38 -3.25000 0.00000 -16.20000 1 0
39 -22.87000 0.00000 9.13000 1 0
40 1.62500 -2.81500 -16.20000 1 0
41 11.43376 -19.80385 9.13000 1 0
42 14.43376 25.00000 -19.94000 1 0
43 -28.86751 0.00000 -19.94000 1 0
44 14.43376 -25.00000 -19.94000 1 0
---------------------- MEMBER CROSS-SECTION PROPERTIES -------------------------
4 NPropSets - Number of member property sets (-)
PropSetID PropD PropThck
(-) (m) (m)
1 6.50000 0.03000 ! Main Column
2 12.00000 0.06000 ! Upper Columns
3 24.00000 0.06000 ! Base Columns
4 1.60000 0.01750 ! Pontoons
---------------------- SIMPLE HYDRODYNAMIC COEFFICIENTS (model 1) --------------
SimplCd SimplCdMG SimplCa SimplCaMG SimplCp SimplCpMG SimplAxCa SimplAxCaMG SimplAxCp SimplAxCpMG
(-) (-) (-) (-) (-) (-) (-) (-) (-) (-)
0.00 0.00 0.00 0.00 1.00 1.00 0.00 0.00 1.00 1.00
---------------------- DEPTH-BASED HYDRODYNAMIC COEFFICIENTS (model 2) ---------
0 NCoefDpth - Number of depth-dependent coefficients (-)
Dpth DpthCd DpthCdMG DpthCa DpthCaMG DpthCp DpthCpMG DpthAxCa DpthAxCaMG DpthAxCp DpthAxCpMG
(m) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-)
---------------------- MEMBER-BASED HYDRODYNAMIC COEFFICIENTS (model 3) --------
25 NCoefMembers - Number of member-based coefficients (-)
MemberID MemberCd1 MemberCd2 MemberCdMG1 MemberCdMG2 MemberCa1 MemberCa2 MemberCaMG1 MemberCaMG2 MemberCp1 MemberCp2 MemberCpMG1 MemberCpMG2 MemberAxCa1 MemberAxCa2 MemberAxCaMG1 MemberAxCaMG2 MemberAxCp1 MemberAxCp2 MemberAxCpMG1 MemberAxCpMG2
(-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) ! Main Column
1 0.56 0.56 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Upper Column 1
2 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Upper Column 2
3 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Upper Column 3
4 0.61 0.61 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Base Column 1
5 0.68 0.68 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Base Column 2
6 0.68 0.68 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Base Column 3
7 0.68 0.68 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Base column cap 1
23 0.68 0.68 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Base column cap 2
24 0.68 0.68 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Base column cap 3
25 0.68 0.68 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Delta Pontoon, Upper 1
8 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Delta Pontoon, Upper 2
9 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Delta Pontoon, Upper 3
10 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Delta Pontoon, Lower 1
11 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Delta Pontoon, Lower 2
12 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Delta Pontoon, Lower 3
13 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Y Pontoon, Upper 1
14 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Y Pontoon, Upper 2
15 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Y Pontoon, Upper 3
16 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Y Pontoon, Lower 1
17 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Y Pontoon, Lower 2
18 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Y Pontoon, Lower 3
19 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Cross Brace 1
20 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Cross Brace 2
21 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 ! Cross Brace 3
22 0.63 0.63 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
-------------------- MEMBERS -------------------------------------------------
25 NMembers - Number of members (-)
MemberID MJointID1 MJointID2 MPropSetID1 MPropSetID2 MDivSize MCoefMod PropPot [MCoefMod=1: use simple coeff table, 2: use depth-based coeff table, 3: use member-based coeff table] [ PropPot/=0 if member is modeled with potential-flow theory]
(-) (-) (-) (-) (-) (m) (switch) (flag)
1 1 2 1 1 1.0000 3 TRUE ! Main Column
2 3 4 2 2 1.0000 3 TRUE ! Upper Column 1
3 5 6 2 2 1.0000 3 TRUE ! Upper Column 2
4 7 8 2 2 1.0000 3 TRUE ! Upper Column 3
5 42 3 3 3 1.0000 3 TRUE ! Base Column 1
6 43 5 3 3 1.0000 3 TRUE ! Base Column 2
7 44 7 3 3 1.0000 3 TRUE ! Base Column 3
23 9 42 3 3 1.0000 3 TRUE ! Base column cap 1
24 10 43 3 3 1.0000 3 TRUE ! Base column cap 2
25 11 44 3 3 1.0000 3 TRUE ! Base column cap 3
8 12 13 4 4 1.0000 3 TRUE ! Delta Pontoon, Upper 1
9 14 15 4 4 1.0000 3 TRUE ! Delta Pontoon, Upper 2
10 16 17 4 4 1.0000 3 TRUE ! Delta Pontoon, Upper 3
11 18 19 4 4 1.0000 3 TRUE ! Delta Pontoon, Lower 1
12 20 21 4 4 1.0000 3 TRUE ! Delta Pontoon, Lower 2
13 22 23 4 4 1.0000 3 TRUE ! Delta Pontoon, Lower 3
14 24 25 4 4 1.0000 3 TRUE ! Y Pontoon, Upper 1
15 26 27 4 4 1.0000 3 TRUE ! Y Pontoon, Upper 2
16 28 29 4 4 1.0000 3 TRUE ! Y Pontoon, Upper 3
17 30 31 4 4 1.0000 3 TRUE ! Y Pontoon, Lower 1
18 32 33 4 4 1.0000 3 TRUE ! Y Pontoon, Lower 2
19 34 35 4 4 1.0000 3 TRUE ! Y Pontoon, Lower 3
20 36 37 4 4 1.0000 3 TRUE ! Cross Brace 1
21 38 39 4 4 1.0000 3 TRUE ! Cross Brace 2
22 40 41 4 4 1.0000 3 TRUE ! Cross Brace 3
---------------------- FILLED MEMBERS ------------------------------------------
2 NFillGroups - Number of filled member groups (-) [If FillDens = DEFAULT, then FillDens = WtrDens; FillFSLoc is related to MSL2SWL]
FillNumM FillMList FillFSLoc FillDens
(-) (-) (m) (kg/m^3)
3 2 3 4 -6.17 1025
3 5 6 7 -14.89 1025
---------------------- MARINE GROWTH -------------------------------------------
0 NMGDepths - Number of marine-growth depths specified (-)
MGDpth MGThck MGDens
(m) (m) (kg/m^3)
---------------------- MEMBER OUTPUT LIST --------------------------------------
0 NMOutputs - Number of member outputs (-) [must be < 10]
MemberID NOutLoc NodeLocs [NOutLoc < 10; node locations are normalized distance from the start of the member, and must be >=0 and <= 1] [unused if NMOutputs=0]
(-) (-) (-)
---------------------- JOINT OUTPUT LIST ---------------------------------------
0 NJOutputs - Number of joint outputs [Must be < 10]
0 JOutLst - List of JointIDs which are to be output (-)[unused if NJOutputs=0]
---------------------- OUTPUT --------------------------------------------------
True HDSum - Output a summary file [flag]
False OutAll - Output all user-specified member and joint loads (only at each member end, not interior locations) [flag]
2 OutSwtch - Output requested channels to: [1=Hydrodyn.out, 2=GlueCode.out, 3=both files]
“ES11.4e2” OutFmt - Output format for numerical results (quoted string) [not checked for validity!]
“A11” OutSFmt - Output format for header strings (quoted string) [not checked for validity!]
---------------------- OUTPUT CHANNELS -----------------------------------------
“Wave1Elev” - Wave elevation at the platform reference point (0, 0)
END of output channels and end of file. (the word “END” must appear in the first 3 columns of this line)

Best regards,

Tomás

Dear Tomas,

Actually, I see that you do have the second-order difference-frequency and sum-frequency potential-flow effects enabled (DiffQTF = SumQTF = 12). This likely explains the mean offset in the surge direction.

I’m still not sure about sway and roll, but perhaps this is related to something else in your model, not set in HydroDyn. To start debugging, I would suggest to disable all platform structural degrees of freedom in ElastoDyn and output the total hydrodynamic loads from HydroDyn (HydroFxi - HydroMzi) to verify that the hydrodynamic loads are only in the surge, pitch, and heave directions.

Best regards,

Dear Jason,

Indeed, I was just looking at the 2ND-ORDER WAVES part, so I missed 2ND-ORDER FLOATING PLATFORM FORCES…

So, I turned off both DiffQTF and SumQTF and the results are now as expected. I just have one doubt regarding this issue. When I have the waves heading at 0º, there’s almost no difference between having DiffQTF and SumQTF turned on or off. On the other hand, if the waves are heading at 180º, the difference is huge. Any hint on why this is happening?

Best regards,

Tomás

Dear Tomás,

Did you change the second-order WAMIT data for the OC4-DeepCwind semi in any way (marin_semi.12d and marin_semi.12s)? The second-order WAMIT data that NREL has provided is only applicable to the zero-degree wave heading. HydroDyn should abort if you select a wave direction outside of the range of the WAMIT data, but there is a known bug in HydroDyn that has not yet been fixed related to the wave directions for the second-order potential-flow solution–see: github.com/OpenFAST/openfast/issues/9. My guess is you have not changed the second-order WAMIT data and–because of the bug–HydroDyn is not functioning properly when the wave direction is not zero, when it should have aborted in this case.

Best regards,

Sir,

I have simulated wave loading by modeling a sea state with significant wave height 5m and peak spectral period of 12.4 sec using the JONSWAP spectrum for bottom fixed monopile. And I observed that inertia force was more than drag force. But my structure is small compared to the wave.

  1. if so Froude-Krylov theory should be used for calculating wave loads. Am I correct? But in hydrodyn manual relative form of Morison equation is there. I am not getting a clear view of the formulation of wave force calculation on bodies

  2. What are the equations FAST is using for calculating water particle velocities and accelerations and wave forces on a bottom fixed and flexible monopile for a random wave?

Thank you.

Dear Dhaneesh,

Here are my answers to your questions:

  1. The strip-theory formulation of FAST’s HydroDyn hydrodynamics module uses the relative form of Morison equation, with extensions for end effects, tapered members, buoyancy, etc. The Froude-Krylov force is the integration of the undisturbed pressure field over the surface of the body, which is captured in Morison’s equation, based on the long wavelength approximation (i.e, the “1” in the (1+Ca) term, but HydroDyn uses a user-specified pressure coefficient, Cp, in place of “1”).

  2. The first-order (linear) terms use Airy wave theory. If the second-order wave kinematics are enabled (added to the first-order terms using perturbation theory), these are implemented using analytical expressions developed by Sharma and Dean.

Best regards,

Dear Jason,

Does hydrodyn consider/model wave-current interaction ? Are there any wave frequency limitations while using current and waves together in regular wave simulations?

Best regards,
Rahul.

Dear @Rahul.Ramachandran,

HydroDyn currently superimposes wave kinematics and current through a simple vector summation. More sophisticated modeling of waves plus current, e.g., to account for the encounter frequency, is not considered. The superimposed wave plus current velocity is used in the strip-theory solution’s calculation of viscous-drag loads.

Best regards,

Dear Jason,

Thank you for your answer.

Rahul.