Implementation of True Linear Hydrodynamic Model in Hydrodyn

Dear @Jason.Jonkman ,
Thanks for you paper entitled “Dynamics of Offshore Floating Wind Turbines—Model Development and Verification” The detailed content improves my understanding of theory implemented in Hydrodyn.
I want to share my understanding here. Could you please correct if there are something misundertood?

In the following, omega denoted angular frequency of incident wave, and beta denotes incident wave direction.

  1. Hydrodyn adopted true linear hydrodynamic model, rather than frequency-domain hydrodynamic model. But the excitation force Xi(omega, beta) used in frequency-domain hydrodynamic model is also used in true linear hydrodynamic model (diffraction term).

  2. For added mass matrix Aij(omega) and damping Bij(omega), true linear hydrodynamic model only uses infinite-frequency limit of Aij(omega). Bij(omega) is actually used to integrate wave radiation-retardation kernel Kij(t).

  3. For Hydrodyn input files, .1 file contains Aij(omega), Bij(omega). Based on above understanding, I would guess only Aij(infinity) is directly used by Hydrodyn, and other Aij(omega) will be ignored. A range of Bij(omega) are used to derive wave radiation-retardation kernel Kij(t) through FFT.

I still have a question irrelevant to Hydrodyn:
When I use BEMRosetta (an open-source software to animate hydrodynamic coefficients) to animate WAMIT results of OC4 platform offered in OpenFAST baseline folder. I found a variable called “Kirf Impulse Response Function”. I checked the WAMIT documentation, .hst, .1 and .3 files do not contain this variable.
4. Do you know what “Kirf Impulse Response Function” is used for? It seems that Hydrodyn does not use this variable.


Dear @Ran.Tu,

I agree with your points 1-3.

Regarding 4, I’m not too familiar with BEMRosetta, but the radiation retardation kernel is sometimes referred to as the radiation impulse response function. The frequency-to-time (F2T) utility of WAMIT provides something similar, although in that case the impulse-response function is called L, where K = d/dt(L) (equation (10) in my Wind Energy journal article that you referenced). So, you should check the documentation of BEMRosetta to figure out which form is being computed.

Best regards,