Dear Dr. Jonkman,
I am currently attempting to reproduce the visous damping coefficients of the UMaine VolturnUS reference platform in the NREL Report (NREL/TP-5000-76773, Allen et al., 2020) using OpenFOAM. Based on the brief description in the report, my goal is to apply the same methodology (or approach) in OpenFOAM to obtain values consistent with those listed in the attached table.
I am using simpleFoam, a steady-state incompressible solver, to apply unit translational (1 m/s) and rotational velocities (1 degee/sec) to calculate the damping coefficients for each degree-of-freedom (DOF).
However, with the limited information provided in the report, it has been challenging to reproduce the consistent results shown in the table above. So I would be greatly appreciate if you could share more details on the methods used to calculate the damping coefficient in OpenFOAM. Since the OpenFOAM results can vary substantially with domain boundary conditions, meshing quality, convergence behavior, and other input settings, additional modeling details or information would be extremely appreciated.
Especially,
-
the specific convergence criteria and procedure used to determine the damping coefficients,
-
critical simulation settings or conditions (e.g., domain boundary conditions for velocity(U) and pressure (p), how to impose the rotational unit flow (angular velocity), meshing or key input parameters, etc.).
-
How was the rotational unit flow (angular velocity) imposed? When using a steady-state solver (simpleFoam), I understand that the unit flow imposed at the “inlet” is used to calculate the forces and moments acting on the surface of the stationary platform body, and the surge motion can be induced by applying a flow through the inlet in the direction normal to the surface. However, to induce the heave or sway motion, under the assumption that the flow is applied normal to the inlet surface, does this mean that the platform inside the domain must be rotated? Also, how was the rotational (angular) velocity imposed from the inlet to induce the pitch, roll, or yaw motions? I was wondering if methods like “codedFixedValue” or “MRF” should be used for simulating the rotational behaviors in the steady-state solver.
-
Could you please share the procedure used to convert the force/ moment from OpenFOAM for each DOF into the damping coefficients shown in the table? Based on the units in the table, it appears that the forces or moments were divided by the square of the velocity (or angular velocity (omega)), but I would like to confirm whether a particular governing equation was used for this calculation.
I apologize for the questions and any requests I have made, and any suggestions or comments you may be willing to share would be greatly appreciated.
Best regards,
Jeongjoo Kim