P. Naaijen
Please Note
17 records found
1
Data Driven Ship Motion Prediction
Forecasting Vessel Response from Radar and Ship Motion History
Industry-standard systems reconstruct the surrounding wave field from X-band navigation radar images, propagate it forward using linear wave theory, and convert it to motion through pre-computed Response Amplitude Operators (RAOs).
Existing machine learning approaches to vessel motion forecasting typically rely on past motion alone, sometimes augmented with scalar wave parameters or sparse wave measurements, and therefore do not exploit the full spatial structure of the surrounding wave field.
No prior study has demonstrated direct prediction of 6-DoF vessel motions from raw X-band radar image sequences within a single end-to-end learnable model.
This thesis addresses that gap by proposing the Spatiotemporal Cross-Attention Transformer (SCAT), a neural architecture that takes sequences of X-band radar images together with 45~seconds of past 6-DoF motion and predicts the next 60~seconds of motion.
A convolutional radar encoder produces a tokenised spatiotemporal representation of the radar history, which is queried through cross-attention by learned future motion tokens.
Both radar and motion representations are conditioned on the relative wave--vessel heading through a learned embedding.
Four model configurations are studied by varying the radar input representation and training loss in order to identify the most effective design.
All variants are pre-trained on synthetic data spanning five JONSWAP sea states and eight wave directions, and then fine-tuned on 20~minutes of operational data recorded aboard a semi-submersible rig in the Norwegian Sea.
Their performance is compared against state-of-the-art phase-resolved predictions provided alongside the same recordings and evaluated on four five-hour blocks covering representative wave headings.
On synthetic data, all four variants achieve heave correlations between 0.89 and 0.94 and 6-DoF mean correlations between 0.87 and 0.91 at the 60-second horizon, showing that the architecture can learn wave-to-motion relationships directly from radar imagery.
After fine-tuning, the selected variant reaches heave correlations of 0.69 to 0.76 across all four real-data blocks and reproduces per-DOF directional patterns consistent with the underlying wave physics.
Joint training across all blocks further improves correlation in three of the four blocks.
The state-of-the-art baseline still leads the 6-DoF mean correlation at 60~seconds by 0.14 to 0.31, indicating that synthetic-to-real domain mismatch remains the main limitation.
SCAT is therefore presented as a proof of concept showing that an attention-based model operating on raw radar imagery and motion history can learn physically meaningful 6-DoF motion predictions without hard-coded wave physics, dispersion relations, or RAOs.
These results establish a viable and complementary path toward data-driven maritime motion prediction and suggest that reducing the synthetic-to-real gap through higher-fidelity pre-training data is the most direct route toward operational deployment.
...
Industry-standard systems reconstruct the surrounding wave field from X-band navigation radar images, propagate it forward using linear wave theory, and convert it to motion through pre-computed Response Amplitude Operators (RAOs).
Existing machine learning approaches to vessel motion forecasting typically rely on past motion alone, sometimes augmented with scalar wave parameters or sparse wave measurements, and therefore do not exploit the full spatial structure of the surrounding wave field.
No prior study has demonstrated direct prediction of 6-DoF vessel motions from raw X-band radar image sequences within a single end-to-end learnable model.
This thesis addresses that gap by proposing the Spatiotemporal Cross-Attention Transformer (SCAT), a neural architecture that takes sequences of X-band radar images together with 45~seconds of past 6-DoF motion and predicts the next 60~seconds of motion.
A convolutional radar encoder produces a tokenised spatiotemporal representation of the radar history, which is queried through cross-attention by learned future motion tokens.
Both radar and motion representations are conditioned on the relative wave--vessel heading through a learned embedding.
Four model configurations are studied by varying the radar input representation and training loss in order to identify the most effective design.
All variants are pre-trained on synthetic data spanning five JONSWAP sea states and eight wave directions, and then fine-tuned on 20~minutes of operational data recorded aboard a semi-submersible rig in the Norwegian Sea.
Their performance is compared against state-of-the-art phase-resolved predictions provided alongside the same recordings and evaluated on four five-hour blocks covering representative wave headings.
On synthetic data, all four variants achieve heave correlations between 0.89 and 0.94 and 6-DoF mean correlations between 0.87 and 0.91 at the 60-second horizon, showing that the architecture can learn wave-to-motion relationships directly from radar imagery.
After fine-tuning, the selected variant reaches heave correlations of 0.69 to 0.76 across all four real-data blocks and reproduces per-DOF directional patterns consistent with the underlying wave physics.
Joint training across all blocks further improves correlation in three of the four blocks.
The state-of-the-art baseline still leads the 6-DoF mean correlation at 60~seconds by 0.14 to 0.31, indicating that synthetic-to-real domain mismatch remains the main limitation.
SCAT is therefore presented as a proof of concept showing that an attention-based model operating on raw radar imagery and motion history can learn physically meaningful 6-DoF motion predictions without hard-coded wave physics, dispersion relations, or RAOs.
These results establish a viable and complementary path toward data-driven maritime motion prediction and suggest that reducing the synthetic-to-real gap through higher-fidelity pre-training data is the most direct route toward operational deployment.
Hydrodynamic Response of a Subsea Piling Rig During Offshore Lowering
Studying dynamic loads and slack formation
The objective of this thesis was to quantify how the rig's specific shape affects its vertical oscillatory motion and to identify slack-inducing conditions in a given sea state. The hydrodynamic response of the rig and resulting line tensions were investigated using time-domain simulations in OrcaFlex, focusing on the deep submergence and positioning/landing phases. A representative Service Operation Vessel (SOV) was modeled using Response Amplitude Operators (RAOs), and the mudmat was treated using a heave-plate analogy with coefficients derived for low oscillation movements.
Three distinct operational scenarios were analyzed to capture the varying physics of the descent:
•
Unbounded Deep-Water Region (Mid-water Transit): In this regime, where the mudmat is sufficiently far from both the surface and the seabed, the rig behaved with stable hydrodynamic forces. Simulations showed that slack events were present but infrequent.
•
Near-Seabed Region: The rig's response when positioned close to the seabed was drastically different. Flow confinement amplified both added mass and damping significantly. This caused the slack probability to increase drastically compared to the unbounded case, consistent with the confinement-dominated regime.
•
Passive Suction Release: Uncontrolled suction release at the maximum upward vessel movement (heave crest) was found to trigger an immediate, full slack event. The sudden release of stored elastic energy in the lifting line caused the mudmat to shoot upward, leading to a complete loss of line tension even before the vessel began its downward motion.
Based on these findings, operational conclusions emphasize the need for caution. The increased hydrodynamic resistance near the seabed makes the system prone to slack, indicating that the rig should be lifted far away from the seabed before horizontal repositioning. Furthermore, controlled or gradual suction venting is essential, preferably timed near a point of neutral vertical vessel movement, to prevent the severe dynamic snap loads associated with sudden suction failure. Future modeling should incorporate porous boundary conditions using suitable open source potential flow solvers and dedicated experiments to better capture complex flow interactions. ...
The objective of this thesis was to quantify how the rig's specific shape affects its vertical oscillatory motion and to identify slack-inducing conditions in a given sea state. The hydrodynamic response of the rig and resulting line tensions were investigated using time-domain simulations in OrcaFlex, focusing on the deep submergence and positioning/landing phases. A representative Service Operation Vessel (SOV) was modeled using Response Amplitude Operators (RAOs), and the mudmat was treated using a heave-plate analogy with coefficients derived for low oscillation movements.
Three distinct operational scenarios were analyzed to capture the varying physics of the descent:
•
Unbounded Deep-Water Region (Mid-water Transit): In this regime, where the mudmat is sufficiently far from both the surface and the seabed, the rig behaved with stable hydrodynamic forces. Simulations showed that slack events were present but infrequent.
•
Near-Seabed Region: The rig's response when positioned close to the seabed was drastically different. Flow confinement amplified both added mass and damping significantly. This caused the slack probability to increase drastically compared to the unbounded case, consistent with the confinement-dominated regime.
•
Passive Suction Release: Uncontrolled suction release at the maximum upward vessel movement (heave crest) was found to trigger an immediate, full slack event. The sudden release of stored elastic energy in the lifting line caused the mudmat to shoot upward, leading to a complete loss of line tension even before the vessel began its downward motion.
Based on these findings, operational conclusions emphasize the need for caution. The increased hydrodynamic resistance near the seabed makes the system prone to slack, indicating that the rig should be lifted far away from the seabed before horizontal repositioning. Furthermore, controlled or gradual suction venting is essential, preferably timed near a point of neutral vertical vessel movement, to prevent the severe dynamic snap loads associated with sudden suction failure. Future modeling should incorporate porous boundary conditions using suitable open source potential flow solvers and dedicated experiments to better capture complex flow interactions.
Offshore Monopile Installation Optimization
A Parametric Study and Optimization Framework to Address Resonance-Induced Operational Constraints
A planar, linearized Lagrangian model was developed to simulate the dynamic response of the rigging system suspended from a vessel-mounted crane. The model includes multiple pendulum bodies and uses a frequency-domain formulation to calculate the transfer function from crane tip motion to resulting sideloads. This transfer function is used to compute the most probable maximum (MPM) sideloads across a range of sea states. A grid-based workability analysis links these loads to operational thresholds, and a numerical optimization was implemented to tune the rigging geometry for maximum uptime.
The analysis identified a dominant resonance near T_p = 8.16 s — the second natural mode of the suspended system — as the primary driver of limit exceedance. By modifying the lengths of the rigging elements, the natural frequencies could be shifted away from this critical range. The optimized configuration improved workability by 10–20 percentage points, depending on the metocean conditions.
Although the simplified model underestimates absolute load levels compared to a detailed OrcaFlex simulation, the relative response behavior is well captured. A conservative threshold correction was applied to enable consistent limit checking. The study shows that dynamic rigging optimization is a viable, low-effort strategy for reducing downtime in offshore monopile installation, and should be considered during project engineering. ...
A planar, linearized Lagrangian model was developed to simulate the dynamic response of the rigging system suspended from a vessel-mounted crane. The model includes multiple pendulum bodies and uses a frequency-domain formulation to calculate the transfer function from crane tip motion to resulting sideloads. This transfer function is used to compute the most probable maximum (MPM) sideloads across a range of sea states. A grid-based workability analysis links these loads to operational thresholds, and a numerical optimization was implemented to tune the rigging geometry for maximum uptime.
The analysis identified a dominant resonance near T_p = 8.16 s — the second natural mode of the suspended system — as the primary driver of limit exceedance. By modifying the lengths of the rigging elements, the natural frequencies could be shifted away from this critical range. The optimized configuration improved workability by 10–20 percentage points, depending on the metocean conditions.
Although the simplified model underestimates absolute load levels compared to a detailed OrcaFlex simulation, the relative response behavior is well captured. A conservative threshold correction was applied to enable consistent limit checking. The study shows that dynamic rigging optimization is a viable, low-effort strategy for reducing downtime in offshore monopile installation, and should be considered during project engineering.
We investigate two MIMO systems, the force-to-motion system and wave-to-motion system, employing Spectral Analysis method and Subspace method, respectively. Both synthetic and real-field data are utilized to study the methods. While Spectral Analysis method demon- strates better accuracy, it relies on the assumption of perfect accuracy of the pre-computed force Response Amplitude Operators (RAOs), which is a limitation of this approach. Sub- space method is less precise, it also requires real-time re-evaluation of the system order due to its sensitivity to wave conditions, making this method less versatile.
Comparing these methods with SISO system method and conventional RAO-based method, Spectral Analysis method emerges as the most accurate, followed by Subspace method. These findings strongly support the potential of MIMO system identification to significantly improve ship motion prediction accuracy.
...
We investigate two MIMO systems, the force-to-motion system and wave-to-motion system, employing Spectral Analysis method and Subspace method, respectively. Both synthetic and real-field data are utilized to study the methods. While Spectral Analysis method demon- strates better accuracy, it relies on the assumption of perfect accuracy of the pre-computed force Response Amplitude Operators (RAOs), which is a limitation of this approach. Sub- space method is less precise, it also requires real-time re-evaluation of the system order due to its sensitivity to wave conditions, making this method less versatile.
Comparing these methods with SISO system method and conventional RAO-based method, Spectral Analysis method emerges as the most accurate, followed by Subspace method. These findings strongly support the potential of MIMO system identification to significantly improve ship motion prediction accuracy.
The frequency-domain method (FD) is chosen for further development. These methods are considered as efficient tools for analyzing first design iterations. Given the ongoing development of new concepts for FOWTs, the industry could benefit from such a tool. Currently, only one such open-sourced method is developed (RAFT), but the method still has room for improvements.
After enhancing RAFT's control modeling, the frequency-dependent aerodynamic added mass and damping coefficients are determined following RAFT's methodology. These coefficients are then implemented in a wave diffraction program to predict the surge, pitch, and mooring line tension response of FOWTs for four different design load cases (DLCs). These results are compared against time-domain (TD) simulations.
The method proved capable of predicting the surge and pitch motions well, but the mooring line tension prediction still needs improvements. ...
The frequency-domain method (FD) is chosen for further development. These methods are considered as efficient tools for analyzing first design iterations. Given the ongoing development of new concepts for FOWTs, the industry could benefit from such a tool. Currently, only one such open-sourced method is developed (RAFT), but the method still has room for improvements.
After enhancing RAFT's control modeling, the frequency-dependent aerodynamic added mass and damping coefficients are determined following RAFT's methodology. These coefficients are then implemented in a wave diffraction program to predict the surge, pitch, and mooring line tension response of FOWTs for four different design load cases (DLCs). These results are compared against time-domain (TD) simulations.
The method proved capable of predicting the surge and pitch motions well, but the mooring line tension prediction still needs improvements.
This study presents an identification procedure to handle the inherent uncertainties of vessel model parameters, aiming to improve vessel motion prediction. The identification procedure identifies the vessel's RAO by the measured response spectrum and nowcast wave spectrum, with the goal of finding the heave and roll natural frequencies. The natural frequencies provide information on the vessel’s parameters. This is used to identify the parameters related to the mass distribution and damping of the vessel. These were found by minimizing a cost function, that quantified the difference between the measured and predicted response spectrum, using an optimization method. Identifiability analyses of the parameters were performed on two case studies.
For the first case study, a synthetic data set is created with the vessel response model to simulate the vessel motions. Tests were conducted with five different wave spectra and several vessel headings, constituting diversified scenarios. The RAO was identified by the measured response, the wave spectrum, and a sinusoidal function to describe the directional dependency of the RAO. Using the synthetic data set, the identification algorithm successfully identified the parameters with good agreement to their actual values. The second case study involved the examination of parameter identification on real onboard vessel motion measurements. In most of the cases, the RAO could be identified from the measurements and the natural heave and roll frequency was found. The identified parameters resulted from the identification procedure and improved the vessel motion prediction compared to the initial prediction, but still, deviations remained. The identified parameters are verified against a different measured data set. The results show that the identified response spectra approach the measured responses, indicating that the identified parameters are reusable.
In summary, it was found that the parameters have a great influence on the output of the vessel response model. Therefore, it is essential to have a thorough understanding of the correct operational parameters for accurate motion prediction. The established identification procedure shows to be a good addition to existing vessel motion models to identify input parameters at relatively low computational cost. ...
This study presents an identification procedure to handle the inherent uncertainties of vessel model parameters, aiming to improve vessel motion prediction. The identification procedure identifies the vessel's RAO by the measured response spectrum and nowcast wave spectrum, with the goal of finding the heave and roll natural frequencies. The natural frequencies provide information on the vessel’s parameters. This is used to identify the parameters related to the mass distribution and damping of the vessel. These were found by minimizing a cost function, that quantified the difference between the measured and predicted response spectrum, using an optimization method. Identifiability analyses of the parameters were performed on two case studies.
For the first case study, a synthetic data set is created with the vessel response model to simulate the vessel motions. Tests were conducted with five different wave spectra and several vessel headings, constituting diversified scenarios. The RAO was identified by the measured response, the wave spectrum, and a sinusoidal function to describe the directional dependency of the RAO. Using the synthetic data set, the identification algorithm successfully identified the parameters with good agreement to their actual values. The second case study involved the examination of parameter identification on real onboard vessel motion measurements. In most of the cases, the RAO could be identified from the measurements and the natural heave and roll frequency was found. The identified parameters resulted from the identification procedure and improved the vessel motion prediction compared to the initial prediction, but still, deviations remained. The identified parameters are verified against a different measured data set. The results show that the identified response spectra approach the measured responses, indicating that the identified parameters are reusable.
In summary, it was found that the parameters have a great influence on the output of the vessel response model. Therefore, it is essential to have a thorough understanding of the correct operational parameters for accurate motion prediction. The established identification procedure shows to be a good addition to existing vessel motion models to identify input parameters at relatively low computational cost.
Ship Motion Predictions in the Time Domain
Using identified linear and non-linear force coefficients
To enable the use of nonlinear forces, the equations of motions are first implemented in the time domain using Cummins’ equations. The linear coefficients in this equation are determined from the frequency dependent coefficients in Acta Auriga’s hydrodynamic database. To ensure the correctness and investigate the limitations of the implemented Cummins equation, the results from the time and frequency domain are extensively compared. Both monochromatic excitations, as well as spectrum (JONSWAP) excitations are considered. In contrast to linear force coefficients, nonlinear force coefficients are generally not known for a given vessel. Implementing these forces into the equations of motion imposes the need for an approximation of these coefficients. Estimates for the linear, quadratic and cubic viscous damping forces are made, using both the measured motions of a vessel operating at sea and the predictions of the corresponding excitation force, made by Next Ocean. A multivariate regression algorithm is used for the identification.
The Cummins equation was successfully implemented. However, the translation from frequency to time domain was more error prone than anticipated; the frequency dependent coefficients need to meet requirements that are typically not met by databases meant for frequency domain calculations only. Furthermore, degrees of freedom without a restoring force showed running away behaviour that could not be negated without adding extra damping. The roll motion was susceptible to instabilities. The exact origin of this instability was looked for, but could not be found. Adding damping resolved the instability. Furthermore, the interpolation in the roll response amplitude operator introduced errors in the initialisation of time domain calculations, which led to inaccurate results when spectrum excitations corresponding to more severe sea states were used. Only under specific conditions, these high energy spectra led to accurate roll predictions. Adding damping made for a close match between frequency and time domain calculation for all degrees of freedom. The complexities mentioned made that an easily scalable algorithm could not be obtained using time domain calculations. Easy scalability be- ing important to Next Ocean means that running time domain simulations in Next Ocean’s product is deemed unrealistic. This also means that no nonlinear forces can be used in real time applications.
In an attempt to improve linear predictions, coefficients in the linear seakeeping model, including the linear damping coefficient, are identified and vessel motions are re-predicted with the added damping and updated linear force coefficients, using the Cummins equation. None of the identified parameters led to better predictions than were obtained by Next Ocean. Identifying and updating parameters was therefore concluded to be not beneficial to the quality of the motion predictions. Adding a fixed amount of linear roll damping, which was not identified from the field data, did lead to improved prediction quality; correlations between predictions and measurements increased by 9.6%.
While the motion prediction could not be improved using the Cummins equation, valuable information on the time domain simulations was obtained. The finding that better predictions were consistently obtained by adding a fixed amount of damping, sparks an opportunity for further research into the ideal amount of damping. This, and more elaborate identification schemes could provide meaningful insight to improve motion predictions. ...
To enable the use of nonlinear forces, the equations of motions are first implemented in the time domain using Cummins’ equations. The linear coefficients in this equation are determined from the frequency dependent coefficients in Acta Auriga’s hydrodynamic database. To ensure the correctness and investigate the limitations of the implemented Cummins equation, the results from the time and frequency domain are extensively compared. Both monochromatic excitations, as well as spectrum (JONSWAP) excitations are considered. In contrast to linear force coefficients, nonlinear force coefficients are generally not known for a given vessel. Implementing these forces into the equations of motion imposes the need for an approximation of these coefficients. Estimates for the linear, quadratic and cubic viscous damping forces are made, using both the measured motions of a vessel operating at sea and the predictions of the corresponding excitation force, made by Next Ocean. A multivariate regression algorithm is used for the identification.
The Cummins equation was successfully implemented. However, the translation from frequency to time domain was more error prone than anticipated; the frequency dependent coefficients need to meet requirements that are typically not met by databases meant for frequency domain calculations only. Furthermore, degrees of freedom without a restoring force showed running away behaviour that could not be negated without adding extra damping. The roll motion was susceptible to instabilities. The exact origin of this instability was looked for, but could not be found. Adding damping resolved the instability. Furthermore, the interpolation in the roll response amplitude operator introduced errors in the initialisation of time domain calculations, which led to inaccurate results when spectrum excitations corresponding to more severe sea states were used. Only under specific conditions, these high energy spectra led to accurate roll predictions. Adding damping made for a close match between frequency and time domain calculation for all degrees of freedom. The complexities mentioned made that an easily scalable algorithm could not be obtained using time domain calculations. Easy scalability be- ing important to Next Ocean means that running time domain simulations in Next Ocean’s product is deemed unrealistic. This also means that no nonlinear forces can be used in real time applications.
In an attempt to improve linear predictions, coefficients in the linear seakeeping model, including the linear damping coefficient, are identified and vessel motions are re-predicted with the added damping and updated linear force coefficients, using the Cummins equation. None of the identified parameters led to better predictions than were obtained by Next Ocean. Identifying and updating parameters was therefore concluded to be not beneficial to the quality of the motion predictions. Adding a fixed amount of linear roll damping, which was not identified from the field data, did lead to improved prediction quality; correlations between predictions and measurements increased by 9.6%.
While the motion prediction could not be improved using the Cummins equation, valuable information on the time domain simulations was obtained. The finding that better predictions were consistently obtained by adding a fixed amount of damping, sparks an opportunity for further research into the ideal amount of damping. This, and more elaborate identification schemes could provide meaningful insight to improve motion predictions.
Single blade installation with a floating monohull crane vessel
Establishing the operational limits while using dynamic controlled taglines
Larger wind turbines are being developed to fill the growing demand for green energy. Wind turbines that need to be installed in the next decade will reach a capacity of 15 MW, with an approximate rotor diameter of 240 meters. The currently used wind turbine installation vessels cannot install these larger wind turbines; hence, investments need to be made. An opportunity arises from the oil and gas industry: investments in this industry are decreasing. Thus, the floating monohull fleet working in this industry will need repurposing in the upcoming decade.
In this thesis, the operational limit of a single blade installation with a floating monohull crane vessel is sought. An installation method is proposed where taglines are used to compensate vessel motions and wind loads on the blade. Motion compensation with taglines is a technology that is already being used in the industry, which makes it attractive to contractors who are hesitant to use new technologies. The proposed installation concept uses three different compensation systems, allowing the blade motions to be decoupled and controlled separately. One of the compensation systems is selected and investigated in more detail in this research. The selected compensation system compensates for the blade motion in y- and yaw-direction with two dynamically controlled taglines horizontally attached from the blade to the crane boom.
The operational limit is determined for the final installation stage, the mating phase between the blade and the hub. The alignment of the blade and the hub is modeled in frequency- and time-domain. From these models, it can be concluded that the control system can limit the y-motions at the blade root below typically accepted motion limits, when the wave direction is limited to an interval of 150-210\deg, for a sea state up to a wave height of $H_s$= 5.6 m and wave peak period $T_p$ of 10 s. However, this is only true when the predisplacement of the blade towards the crane boom is larger than the maximum crane tip displacement.
Based on these preliminary findings, it seems that floating installation of turbine blades could be feasible with proper tagline control. Further investigation, involving dynamic control in all degrees of freedom, is therefore recommended. ...
Larger wind turbines are being developed to fill the growing demand for green energy. Wind turbines that need to be installed in the next decade will reach a capacity of 15 MW, with an approximate rotor diameter of 240 meters. The currently used wind turbine installation vessels cannot install these larger wind turbines; hence, investments need to be made. An opportunity arises from the oil and gas industry: investments in this industry are decreasing. Thus, the floating monohull fleet working in this industry will need repurposing in the upcoming decade.
In this thesis, the operational limit of a single blade installation with a floating monohull crane vessel is sought. An installation method is proposed where taglines are used to compensate vessel motions and wind loads on the blade. Motion compensation with taglines is a technology that is already being used in the industry, which makes it attractive to contractors who are hesitant to use new technologies. The proposed installation concept uses three different compensation systems, allowing the blade motions to be decoupled and controlled separately. One of the compensation systems is selected and investigated in more detail in this research. The selected compensation system compensates for the blade motion in y- and yaw-direction with two dynamically controlled taglines horizontally attached from the blade to the crane boom.
The operational limit is determined for the final installation stage, the mating phase between the blade and the hub. The alignment of the blade and the hub is modeled in frequency- and time-domain. From these models, it can be concluded that the control system can limit the y-motions at the blade root below typically accepted motion limits, when the wave direction is limited to an interval of 150-210\deg, for a sea state up to a wave height of $H_s$= 5.6 m and wave peak period $T_p$ of 10 s. However, this is only true when the predisplacement of the blade towards the crane boom is larger than the maximum crane tip displacement.
Based on these preliminary findings, it seems that floating installation of turbine blades could be feasible with proper tagline control. Further investigation, involving dynamic control in all degrees of freedom, is therefore recommended.
Currently operations are planned based on DP capability plots and experience of captain and DPO. DP capability plots have little operational value as this is a static calculation and only provide information for average station keeping capability. During operations, the displacements made by the vessel around the DP set-point, also referred to as DP offset, are of great importance to determine the operability of an operation. Currently, the only way of calculating the DP offset is by conducting extensive time domain simulations, which are hard to integrate in the operational workflow of a DP vessel involved in walk-to-work operations. Therefore, a new approach is developed which predicts the vessel’s DP offset in the frequency domain, which enables a quick and robust calculation of the DP offset which is suited to merge into the on-board workflow. A frequency domain model is per definition a linear model. This leads to the main challenge of this research. A vessel operating on DP is non-linear. Currently there is no insight in what the effect is of non-linear components present in a DP system, on the linear approximation of a frequency domain model.
To investigate the effect of non-linear components onto the DP frequency domain model, a time domain model is developed that is capable of systematically enabling/disabling different non-linear components. The time domain model will serve as the ’truth’ in this research as no actual vessel data is available. Furthermore, this helps identify the effects more easily, as the input for both models are identical. From the time domain model transfer functions can be derived that serve as the basis for the frequency domain model. The transfer function is a linear relation between two variables. In this case, between second order wave drift forces and displacement of the vessel in surge, sway and yaw direction. The following non-linear components are investigated in this research: Thruster ramp up, thruster turning rate, forbidden zones, saturation and thruster allocation. Thruster allocation is present in each model that will be tested, as this is an essential part of a DP system.
Using two methods of determining transfer functions the model and the effects of all non-linear components are tested. The model is subjected to a variety sea-state, with different wave directions. Both methods offer similar results even though different approaches to determine the transfer functions are used. The selected method is capable of accurately predicting vessel offsets, although some extreme offsets are not captured.
It is concluded that the presence of non-linear components have little to no effect on the DP offset as calculated by the time domain model. Because natural frequencies characteristic to these non-linear components are expected to exist at much higher frequencies that naturally present in second order wave drift forces. Thus, making a linear frequency domain model suitable for DP offset forecasting. It is advised to investigate the effect of including 2D input spectra as this is expected to improve the current model.
...
Currently operations are planned based on DP capability plots and experience of captain and DPO. DP capability plots have little operational value as this is a static calculation and only provide information for average station keeping capability. During operations, the displacements made by the vessel around the DP set-point, also referred to as DP offset, are of great importance to determine the operability of an operation. Currently, the only way of calculating the DP offset is by conducting extensive time domain simulations, which are hard to integrate in the operational workflow of a DP vessel involved in walk-to-work operations. Therefore, a new approach is developed which predicts the vessel’s DP offset in the frequency domain, which enables a quick and robust calculation of the DP offset which is suited to merge into the on-board workflow. A frequency domain model is per definition a linear model. This leads to the main challenge of this research. A vessel operating on DP is non-linear. Currently there is no insight in what the effect is of non-linear components present in a DP system, on the linear approximation of a frequency domain model.
To investigate the effect of non-linear components onto the DP frequency domain model, a time domain model is developed that is capable of systematically enabling/disabling different non-linear components. The time domain model will serve as the ’truth’ in this research as no actual vessel data is available. Furthermore, this helps identify the effects more easily, as the input for both models are identical. From the time domain model transfer functions can be derived that serve as the basis for the frequency domain model. The transfer function is a linear relation between two variables. In this case, between second order wave drift forces and displacement of the vessel in surge, sway and yaw direction. The following non-linear components are investigated in this research: Thruster ramp up, thruster turning rate, forbidden zones, saturation and thruster allocation. Thruster allocation is present in each model that will be tested, as this is an essential part of a DP system.
Using two methods of determining transfer functions the model and the effects of all non-linear components are tested. The model is subjected to a variety sea-state, with different wave directions. Both methods offer similar results even though different approaches to determine the transfer functions are used. The selected method is capable of accurately predicting vessel offsets, although some extreme offsets are not captured.
It is concluded that the presence of non-linear components have little to no effect on the DP offset as calculated by the time domain model. Because natural frequencies characteristic to these non-linear components are expected to exist at much higher frequencies that naturally present in second order wave drift forces. Thus, making a linear frequency domain model suitable for DP offset forecasting. It is advised to investigate the effect of including 2D input spectra as this is expected to improve the current model.
Feed forward control of U anti-roll tanks
Research on the effect of using a wave prediction system for the control of an active U anti-roll tank on the workability of an SOV operating at zero speed
Floating Installation of Windturbine Towers
A Conceptual Design Study
This novel concept requires a multibody modeling approach to perform a dynamic loads and response analysis, as the stiffness between the floating platform and the counter weight is provided by chains. Additional design criteria are required for the counter weight system dependent on a combination of chain capacity and maintaining positive tension in all of the suspension lines. To satisfy these design criteria a global hydrodynamic load and response analysis of the floating system is performed. In this concept, the counter weight depth contributes significantly to the dynamic properties of the system and therefore a parametric study is conducted. The global response parameters of the rigid-body motion, natural frequencies, nacelle accelerations, counter weight chain tensions, and maximum platform-pitch angles are compared. Following the parametric study, an ultimate limit state analysis is conducted on the original and alternative designs. Design recommendations are made for the counter weight depth and configuration of the suspension system layout.
...
This novel concept requires a multibody modeling approach to perform a dynamic loads and response analysis, as the stiffness between the floating platform and the counter weight is provided by chains. Additional design criteria are required for the counter weight system dependent on a combination of chain capacity and maintaining positive tension in all of the suspension lines. To satisfy these design criteria a global hydrodynamic load and response analysis of the floating system is performed. In this concept, the counter weight depth contributes significantly to the dynamic properties of the system and therefore a parametric study is conducted. The global response parameters of the rigid-body motion, natural frequencies, nacelle accelerations, counter weight chain tensions, and maximum platform-pitch angles are compared. Following the parametric study, an ultimate limit state analysis is conducted on the original and alternative designs. Design recommendations are made for the counter weight depth and configuration of the suspension system layout.
Rudder loads of moored FPSO's in waves and current
A study by experimental means
Deterministic Wave Prediction as Support for Launch & Recovery Operations
A study into applications of the waveradar in the Feadship Comfort System
For these operations, experience-based decision making does however not always guarantee that the operations are carried out in the most safe and effective manner possible. Also, the level of comfort as experienced by a yacht-owner or guests can be directly affected by this type of decision making. With this in mind, Feadship aims to develop a Feadship Comfort System (FCS), which can contribute to the decision making process. One of the elements of this system is ought to be a waveradar. The aim of this thesis is to describe how the waveradar can be included effectively in the FCS, and how it can be used for the above mentioned operations.
The waveradar is in basis a common navigation radar that can perform measurements of the surface elevation in the direct surrounding of a vessel. These measurements can be used to predict the ship's motions up to approximately two minutes ahead in time. With this, windows of opportunity can be distinguished when operations can be carried out best. These windows depend on the corresponding limiting criteria that are set on the ship's motions. A major issue is however that sea trials have shown that the roll prediction is inaccurate.
Due to this problem, two methods are evaluated in this thesis to improve the accuracy of the roll prediction. These methods are set up such that they are practical in use and can be applied directly to the deterministic wave predictions from the waveradar. The first method is scaling the response based on the roll motion history, whereas the second method is based on scaling the roll damping. The latter is effectively linearising the viscous contribution to the roll damping.
As no data from sea-trials is present to compare these methods, a benchmark calculation is carried out based on Cummins approach. In these time-domain calculations, a linear and non-linear roll damping coefficient are taken into account, which are based on a decay test that is available at De Voogt Naval Architects (DVNA). This follows from model tests performed by MARIN.
Comparing both methods to this benchmark calculation showed that the method of scaled damping resulted in the best approximation, with a Pearson correlation coefficient of (only) 33 %. Also, in terms of the behaviour of the motion envelope, this method showed the closest approximation of the two methods. It is found that the way in which both methods influence the RAO of the roll motion has a great influence on the results found. More suitable methods that should provide a better approximation are suggested from this, such as a frequency dependent scaling of the roll damping. This is however not discussed further.
With the method of scaled damping, the practical application of the waveradar within the FCS is evaluated. For the common operations of yachts and naval vessels, limiting criteria on ship motions are obtained from literature or stated by the author. The majority of these criteria are RMS values, for which a method is suggested to use them real-time. In this method, the most probable maxima are obtained from the RMS criteria, which are then applied as real-time criteria. This showed impractically large results, which showed that this method does not suffice.
From this conclusion it follows that the analysis of the practical application of the waveradar is only carried out for the (dis)embarkment of small craft. From this, design considerations for the FCS are obtained, from which the main result is that heading suggestion is one of the most important tasks that the system needs to be capable of.
Finally, a performance monitor has been evaluated for yachts at anchor. This performance monitor displays the predicted Comfort Rating (CR) real time, which can be used by the yachts crew to evaluate the current anchor location. This is a specific desire of Feadship and is developed next to the other operations. For this application, heading suggestion is also one of the most important tasks that the FCS needs to be capable of.
...
For these operations, experience-based decision making does however not always guarantee that the operations are carried out in the most safe and effective manner possible. Also, the level of comfort as experienced by a yacht-owner or guests can be directly affected by this type of decision making. With this in mind, Feadship aims to develop a Feadship Comfort System (FCS), which can contribute to the decision making process. One of the elements of this system is ought to be a waveradar. The aim of this thesis is to describe how the waveradar can be included effectively in the FCS, and how it can be used for the above mentioned operations.
The waveradar is in basis a common navigation radar that can perform measurements of the surface elevation in the direct surrounding of a vessel. These measurements can be used to predict the ship's motions up to approximately two minutes ahead in time. With this, windows of opportunity can be distinguished when operations can be carried out best. These windows depend on the corresponding limiting criteria that are set on the ship's motions. A major issue is however that sea trials have shown that the roll prediction is inaccurate.
Due to this problem, two methods are evaluated in this thesis to improve the accuracy of the roll prediction. These methods are set up such that they are practical in use and can be applied directly to the deterministic wave predictions from the waveradar. The first method is scaling the response based on the roll motion history, whereas the second method is based on scaling the roll damping. The latter is effectively linearising the viscous contribution to the roll damping.
As no data from sea-trials is present to compare these methods, a benchmark calculation is carried out based on Cummins approach. In these time-domain calculations, a linear and non-linear roll damping coefficient are taken into account, which are based on a decay test that is available at De Voogt Naval Architects (DVNA). This follows from model tests performed by MARIN.
Comparing both methods to this benchmark calculation showed that the method of scaled damping resulted in the best approximation, with a Pearson correlation coefficient of (only) 33 %. Also, in terms of the behaviour of the motion envelope, this method showed the closest approximation of the two methods. It is found that the way in which both methods influence the RAO of the roll motion has a great influence on the results found. More suitable methods that should provide a better approximation are suggested from this, such as a frequency dependent scaling of the roll damping. This is however not discussed further.
With the method of scaled damping, the practical application of the waveradar within the FCS is evaluated. For the common operations of yachts and naval vessels, limiting criteria on ship motions are obtained from literature or stated by the author. The majority of these criteria are RMS values, for which a method is suggested to use them real-time. In this method, the most probable maxima are obtained from the RMS criteria, which are then applied as real-time criteria. This showed impractically large results, which showed that this method does not suffice.
From this conclusion it follows that the analysis of the practical application of the waveradar is only carried out for the (dis)embarkment of small craft. From this, design considerations for the FCS are obtained, from which the main result is that heading suggestion is one of the most important tasks that the system needs to be capable of.
Finally, a performance monitor has been evaluated for yachts at anchor. This performance monitor displays the predicted Comfort Rating (CR) real time, which can be used by the yachts crew to evaluate the current anchor location. This is a specific desire of Feadship and is developed next to the other operations. For this application, heading suggestion is also one of the most important tasks that the FCS needs to be capable of.