model Validation

Dear Satish,

Your waves are not sinusoidal because you’ve enabled second-order wave kinematics. Disabling second-order wave kinematics (WvDiffQTF = WvSumQTF = False) will result in a sinusoidal wave elevation.

Best regards,

Dear Jason,

Thanks for your reply.

Okay, got that.

  • Does this result attached is correct? (I have attached plots resulted my simulation and the referred simulation results)

  • Do I get any reference where I can validate this OpenFast model?

Thanks,
Satish J


Dear Satish,

Your results make sense. You are getting a lower amplitude of mudline fore-aft bending moment because version of HyrdoDyn you are using in OpenFAST does not inherently consider wave stretching, which are included in the results you are comparing to.

If you need to match the OC3 results exactly, until wave stretching is added to OpenFAST*, you could apply wave stretching yourself and use the corresponding wave kinematics within OpenFAST using the WaveMod = 6 option of HydroDyn. This would require that you have wave kinematics at all hydrodynamic analysis nodes along the monopile.

Best regards,

Footnote:
*NREL did have a project to add wave stretching into OpenFAST, but the project ran out of funding before a successful completion, and the solution was not numerically stable for structural flexible substructures; and we have not yet been funded to restart and finish the project.

Dear Jason,

Thanks for your reply.

Okay, Then Is the solution I obtained is acceptable without wheeler stretching as it is getting difference of around 20000 KNm? Any recommendation on going ahead with this solution or the software version itself (Openfast)?

I can use all these in Fastv8 right? If so, could you provide me the link where I can download the Fastv8?

Thanks,
Satish J

Dear Satish,

Not including Wheeler stretching will result in a lower amplitude of wave loading. But Wheeler stretching is known to overpredict the wave loading amplitude anyway. I already suggested the use of WaveMod = 6 in OpenFAST as a way to apply stretching, if needed.

I don’t see any reason for you to run FAST v8 since you are already running a newer version–OpenFAST. And downgrading to FAST v8 would require some small input file changes. FAST v8 does not support Wheeler stretching either (but FAST v7 did). Regardless, FAST v8 is available here: drive.google.com/drive/folders/ … sp=sharing.

Best regards,

Dear Jason,

Thanks for your suggestions.

Okay, Could you please give suggestions on how to generate wave kinematics at all hydrodynamic analysis nodes along the monopile?

Thanks,
Satish J

Dear Satish,

The hydrodynamic analysis node locations are written to the HydroDyn summary file. You’d presumably need access to a software that can calculate wave kinematics (including stretching) at these locations. NREL has not provided such a software.

The HydroDyn module of OpenFAST can output wave kinematics at the various hydrodynamic analysis nodes. So, with a bit a scripting, you could probably output these, apply a stretching theory, and write the stretched wave data to the wave kinematics input file formats the WaveMod = 6 option requires.

Best regards,

Dear Jason,

Thanks for your suggestions.

Thanks,
Satish J

Dear Jason,

Could I know where can I get the wave kinematics input file format when WaveMod=6.

Thanks,
Satish J

Dear Satish,

You’ll need 8 files to specify user-generated wave kinematics via WaveMod = 6: 3 for the wave particle velocities (*.Vxi, *.Vyi, .Vyi), 3 for the wave particle accelerations (.Axi, .Ayi, .Azi), 1 for dynamic pressure (.DynP), and 1 for the wave elevation (.Elev). Examples of these files are provided in Test_004 in the standalone HydroDyn CertTest from FAST v8. The draft HydroDyn User’s Guide and Theory Manual provides some more information on these input files. The HydroDyn archive and documentation from FAST v8 are available on my Google drive: drive.google.com/drive/folders/ … sp=sharing.

Best regards,

Dear Jason,

Thanks for your reply.

Could you please share the Test_004 in the standalone HydroDyn CertTest from FAST v8 as I couldn’t find the file in which you provided the link below.

Thanks,
Satish J

Dear Satish,

HD_v2.05.00 in the link below is a self-extracting ZIP archive. Once you extract the files, the Test_004 files should be found in the CertTest directory.

Best regards,

Dear Jason,

Thanks for your reply.

Could I know what does ########## indicate in the input files?

Thanks,
Satish J

Dear Satish,

As explained in the HydroDyn documentation, a nonnumeric string in the WaveMod = 6 wave-kinematics files designates that the node is outside of the water at that time step (above the instantaneous water elevation or below the seabed).

Best regards,

Dear Jason,

Thanks for your reply.
Could you please suggest me on which software to use to generate all wave kinematics.

Thanks,
Satish J

Dear Satish,

As I mentioned before, NREL has not provided such a software. And we are not in a position to recommend one from elsewhere.

Best regards,

Dear Jason,

Thanks for your reply.

  1. I have generated wave kinematics externally. I have velocities (x,y,z) and acceleration (x,y,z) at each node at each time step. So my doubt is how do I need to consider these values in .Vxi, .Vyi, .Vzi, .Axi, .Ayi, .Azi as each row represents each time step.

  2. Can I consider dynamic pressure file (.DynP) as single column of zero as there is no tapered monopile section (no change in diameter of MP)?

Thanks,
Satish J

Dear Satish,

It sounds like you have all of the data you need (except dynamic pressure). So, all you should need to do is structure the data into the file format readable by HydroDyn. A monopile with no taper would not need the dynamic pressure, so I agree, you can set those to zero.

Best regards,

Dear Jason,

Thanks for your reply.

Regarding the structuring these files, I have considered columns as nodes and rows as time step. Here, in my case I have 41 nodes and total time 60 sec time step of 0.005sec, so 41 columns and 12001 rows. Is that correct way of structuring?

Thanks,
Satish J

Dear Satish,

Yes, that is correct except that there need only be WaveTMax/WaveDT = 12000 rows in your case because no wave kinematics definition is only needed from time = 0 to time = WaveTMax - WaveDT.

Best regards,