Skip to content

Regression tests for the MoorDyn SeaState coupling #151

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 5 commits into from
Apr 28, 2025
Merged

Regression tests for the MoorDyn SeaState coupling #151

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
4059be6
Update WaterKin file comments in md_VIV
RyanDavies19 Apr 25, 2025
6aa89e7
Add waterkin2 test for hybrid MD SS coupling
RyanDavies19 Apr 17, 2025
dcf9ecd
Update md_WaterKin2 test for testing wavekin
RyanDavies19 Apr 18, 2025
b5c3bcb
md_waterkin2 run in double precision
RyanDavies19 Apr 23, 2025
40df91a
Adds waterkin3 test for MoorDyn
RyanDavies19 Apr 24, 2025
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
14 changes: 7 additions & 7 deletions modules/moordyn/md_VIV/WaterKin.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,20 +1,20 @@
MoorDyn Waves and Currents input file
...any notes here...
--------------------------- WAVES -------------------------------------
0 WaveKinMod - type of wave input {0 no waves; 3 set up grid of wave data based on time series}
"" WaveKinFile - file containing wave elevation time series at 0,0,0
0.5 dtWave - time step to use in setting up wave kinematics grid (s)
0 WaveDir - wave heading (deg)
0 WaveKinMod - type of wave input {0 no waves; 1 use the old method; 2 use the hybrid method}
"" WaveKinFile - file containing wave elevation time series at 0,0,0 # Ignored if WaveKinMod = 2
0 dtWave - time step to use in setting up wave kinematics grid (s) # Ignored if WaveKinMod = 2
0 WaveDir - wave heading (deg) # Ignored if WaveKinMod = 2
2 - X wave input type (0: not used; 1: list values in ascending order; 2: uniform specified by -xlim, xlim, num)
-24, 150, 100 - X wave grid point data
2 - Y wave input type (0: not used; 1: list values in ascending order; 2: uniform specified by -xlim, xlim, num)
-100, 100, 5 - Y wave grid point data
2 - Z wave input type (0: not used; 1: list values in ascending order; 2: uniform specified by -xlim, xlim, num)
-600, 0, 60 - Z wave grid point data
--------------------------- CURRENT -------------------------------------
1 CurrentMod - type of current input {0 no current; 1 steady current profile described below}
z-depth x-current y-current
(m) (m/s) (m/s)
1 CurrentMod - type of current input {0 no current; 1 steady current profile described below; 2 hybrid method}
z-depth x-current y-current # Ignored if CurrentMod = 2
(m) (m/s) (m/s) # Ignored if CurrentMod = 2
-11.0 0.2 0
-10.775510204081632 0.2 0
-10.551020408163266 0.2 0
Expand Down
71 changes: 71 additions & 0 deletions modules/moordyn/md_waterkin2/SeaState.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,71 @@
------- SeaState Input File ----------------------------------------------------
For the hybrid configuration with MoorDyn, the Spatial Discretization values are not used, and the Outputs are not necessarily reflective of what is run in MoorDyn.
False Echo - Echo the input file data (flag)
---------------------- ENVIRONMENTAL CONDITIONS --------------------------------
1025.0 WtrDens - Water density (kg/m^3)
600.0 WtrDpth - Water depth (meters) relative to MSL
"DEFAULT" 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]
---------------------- SPATIAL DISCRETIZATION ---------------------------------------------------
10 X_HalfWidth – Half-width of the domain in the X direction (m) [>0, NOTE: X[nX] = nX*dX, where nX = {-NX+1,-NX+2,…,NX-1} and dX = X_HalfWidth/(NX-1)]
10 Y_HalfWidth – Half-width of the domain in the Y direction (m) [>0, NOTE: Y[nY] = nY*dY, where nY = {-NY+1,-NY+2,…,NY-1} and dY = Y_HalfWidth/(NY-1)]
600 Z_Depth – Depth of the domain the Z direction (m) relative to SWL [0 < Z_Depth <= WtrDpth+MSL2SWL; "default": Z_Depth = WtrDpth+MSL2SWL; Z[nZ] = ( COS( nZ*dthetaZ ) – 1 )*Z_Depth, where nZ = {0,1,…NZ-1} and dthetaZ = pi/( 2*(NZ-1) )]
3 NX – Number of nodes in half of the X-direction domain (-) [>=2]
3 NY – Number of nodes in half of the Y-direction domain (-) [>=2]
50 NZ – Number of nodes in the Z direction (-) [>=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, 7: wave frequency components [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]
600 WaveTMax - Analysis time for incident wave calculations (sec) [unused when WaveMod=0; determines WaveDOmega=2Pi/WaveTMax in the IFFT]
0.1 WaveDT - Time step for incident wave calculations (sec) [unused when WaveMod=0 or 7; 0.1<=WaveDT<=1.0 recommended; determines WaveOmegaMax=Pi/WaveDT in the IFFT]
5.0 WaveHs - Significant wave height of incident waves (meters) [used only when WaveMod=1, 2, or 3]
6 WaveTp - Peak-spectral period of incident waves (sec) [used only when WaveMod=1 or 2]
"DEFAULT" 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, 6, or 7]
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, 6, or 7]
0 WaveDir - Incident wave propagation heading direction (degrees) [unused when WaveMod=0, 6 or 7]
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]
"RANLUX" WaveSeed(2) - Second random seed of incident waves [-2147483648 to 2147483647] for intrinsic pRNG, or an alternative pRNG: "RanLux" (-) [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, 6 or 7]
---------------------- 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]
---------------------- CONSTRAINED WAVES ---------------------------------------
0 ConstWaveMod - Constrained wave model: 0=none; 1=Constrained wave with specified crest elevation, alpha; 2=Constrained wave with guaranteed peak-to-trough crest height, HCrest (flag)
1 CrestHmax - Crest height (2*alpha for ConstWaveMod=1 or HCrest for ConstWaveMod=2), must be larger than WaveHs (m) [unused when ConstWaveMod=0]
60 CrestTime - Time at which the crest appears (s) [unused when ConstWaveMod=0]
0 CrestXi - X-position of the crest (m) [unused when ConstWaveMod=0]
0 CrestYi - Y-position of the crest (m) [unused when ConstWaveMod=0]
---------------------- CURRENT ------------------------------------------------- [unused with WaveMod=6]
1 CurrMod - Current profile model {0: none=no current, 1: standard, 2: user-defined from routine UserCurrent} (switch)
10 CurrSSV0 - Sub-surface current velocity at still water level (m/s) [used only when CurrMod=1]
0 CurrSSDir - Sub-surface current heading direction (degrees) or DEFAULT (string) [used only when CurrMod=1]
10 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]
---------------------- MacCamy-Fuchs diffraction model -------------------------
0 MCFD - MacCamy-Fuchs member radius (ignored if radius <= 0) [must be 0 when WaveMod 0 or 6]
---------------------- OUTPUT --------------------------------------------------
False SeaStSum - Output a summary file [flag]
2 OutSwtch - Output requested channels to: [1=SeaState.out, 2=GlueCode.out, 3=both files]
"E15.7e2" 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!]
0 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]
0 NWaveKin - Number of points where the wave kinematics can be output (-) [maximum of 9 output locations]
0 WaveKinxi - List of xi-coordinates for points where the wave kinematics can be output (meters) [NWaveKin points, separated by commas or white space; usused if NWaveKin = 0]
0 WaveKinyi - List of yi-coordinates for points where the wave kinematics can be output (meters) [NWaveKin points, separated by commas or white space; usused if NWaveKin = 0]
0 WaveKinzi - List of zi-coordinates for points where the wave kinematics can be output (meters) [NWaveKin points, separated by commas or white space; usused if NWaveKin = 0]
---------------------- OUTPUT CHANNELS -----------------------------------------
END of output channels and end of file. (the word "END" must appear in the first 3 columns of this line)
18 changes: 18 additions & 0 deletions modules/moordyn/md_waterkin2/WaterKin.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,18 @@
MoorDyn Waves and Currents input file
This uses the hybrid approach, where wave and current init information is set by SeaState, and the grid is set by MoorDyn
--------------------------- WAVES -------------------------------------
2 WaveKinMod - type of wave input {0 no waves; 1 use the old method; 2 use the hybrid method}
"" WaveKinFile - file containing wave elevation time series at 0,0,0 # Ignored if WaveKinMod = 2
0 dtWave - time step to use in setting up wave kinematics grid (s) # Ignored if WaveKinMod = 2
0 WaveDir - wave heading (deg) # Ignored if WaveKinMod = 2
2 - X wave input type (0: not used; 1: list values in ascending order; 2: uniform specified by -xlim, xlim, num)
-800, 10, 101 - X wave grid point data
2 - Y wave input type (0: not used; 1: list values in ascending order; 2: uniform specified by -xlim, xlim, num)
-5, 5, 3 - Y wave grid point data
2 - Z wave input type (0: not used; 1: list values in ascending order; 2: uniform specified by -xlim, xlim, num)
-600, 0, 61 - Z wave grid point data
--------------------------- CURRENT -------------------------------------
2 CurrentMod - type of current input {0 no current; 1 steady current profile described below; 2 hybrid method}
2 - Z wave input type (0: not used; 1: list values in ascending order; 2: uniform specified by -xlim, xlim, num) # Ignored if CurrentMod = 1
-600, 0, 50 - Z wave grid point data # Ignored if CurrentMod = 1
--------------------- need this line ------------------
Loading