Determination of the Response Amplitude Operator(s) of an FPSO

The driving necessities of cost reduction and the need to develop fields at ever increasing water depths have led to the use of floating structures. Among these structures are the Floating Production Storage and Offloading (FPSO) units whose motion analysis is considered in this paper. In actual environmental condition, it is required to accurately determine or predict large amplitude motion of the FPSO before any offshore operation. This paper seeks to present a detailed method of computing the Response Amplitude Operator(s) (RAOs) for the six (6) degrees of freedom using ANSYS AQWA. The results indicate for Heave motion a tendency for the heave peak value to move slightly higher dimensionless encounter-frequency as the wave moves from Head sea to Beam sea direction. A MATLAB source code was developed to validate the result for heave motion at head sea. Although a small difference in predicted heave motion occurred, it is pertinent to note that the comparisons between results gen-erated in the MATLAB program and ANSYS AQWA demonstrate generally good agreement, and the roll response of the FPSO is noted to be critical.


Introduction
Floating Production Storage and Offloading (FPSO) units as the name implies are used in the offshore industry for the storage and processing of oil and gas.
Her storage capacity is large with enough deck area space to accommodate the production equipment as shown in Figure 1. This floating structure has the capability to produce, store and offload the oil, but does not have the capability to drill. In actual condition, the FPSO will respond to environmental loads due to wind, waves, and currents, which make the FPSO be on the motion. The interest in the FPSO motion characteristics stems from the need to predict and control large-amplitude motions that can reduce the operational level and affects the performance of the process.
A generic procedure for predicting the FPSO motions is to have an input-output system, where the ship response is predicted from the input, usually, of an irregular seaway. Such a system is based on a transfer function, also known as "response amplitude operator" (RAO). The validity of this operator is based on the assumption of linearity between the wave excitation and the system's response [1]. As such, the RAO(s) is usually determined in a regular seaway for various wave frequencies and directions either analytical, experimentally, or via numerical simulations to predict the ship motion response, such as the Surge, Sway, Heave, Roll, Pitch and Yaw motion as described in Figure 2.
Hence, the RAOs serve to linearly relate the frequency components of one or more of the ship responses to the same components in the wave spectrum and describe how the response of the FPSO varies with the frequency.

Methods of Analyzing Hydrodynamic Response
Hydrodynamic responses of a floating structure are analyzed using the best available methods. The common method is by conducting a model test in an ocean basin, which gives more accurate prediction of the ship behavior. Over time with the introduction of computer programs and computing technology, numerical simulation became another method used to analyze the motion responses of ocean structures in waves. The role of numerical simulations in the engineering design process is constantly increasing since the virtual test is conducted in controlled environmental conditions and the amount of information available is in orders of magnitude higher than any complex physical testing. The increasing availability of powerful computational resources has recently attracted a large amount of research in the field of the simulation of waves and wave-structures interaction processes.
Computer programs for predicting the motion response of a floating structure are developed based on potential flow theory. According to [2], a popular tool  for computing RAOs of offshore structure is ANSYS AQWA, which implements a three-dimensional panel code to compute the linear wave forcing and motion characteristics of offshore structures in the frequency domain. Other predictors of motion responses are for time domain [3], strip theory [4], and 3-D frequency domain [3]. Most of the sevirtual tools are based on different theories, and most often their predictions have some discrepancies. These seeming in consistencies in numerical predictors [5] are not due to the hydrodynamic theory applied, but due to over-simplifying assumptions during modeling. Nonetheless, numerical simulation is still used independently or in conjunction with veritable test data to characterize most engineering systems.

FPSO Characteristics
The FPSO surface geometry and mass characteristics are input to ANSYS AQWA.
These parameters are given in Table 1.   The ship is assumed to be on even keel in the calculation of the products of inertia;  Moments of inertia and radii of gyration are taken with respect to the CG of the loaded ship.

Analysis Flow
The analysis to determine the response behavior of the FPSO begins with modeling of the FPSO body hull. In general, the body surface geometry and mass distribution of the FPSO model are first described using a panel model, and are followed by the description of the environment, as shown in Figure 3. The FPSO has been modelled accurately using the design modeler in the ANSYS workbench [6].   In ANSYS AQWA, the thickness of the model is 0 m and it is sliced through the waterline (draft). The origin of the global coordinate is located at LCG in X-direction, CL in Y-direction and the distance of KG from the water surface in Z-direction. Consequently, all the hydrodynamic quantities are evaluated at the center of gravity of the vessel in ANSYS AQWA. The phase is defined relative to the phase of an incident wave at the origin of the global coordinates system in X, Y and the free surface in Z. The radiation and diffraction velocity potentials on the body wetted surface are determined from the solution of an integral equation obtained by using Green's theorem with the free-surface source-potential as the DOI: 10.4236/eng.2019.119038 546 Engineering

Coordinate System Definition
Green function. All of the quantities, like the exciting forces, motions, hydrodynamic pressure and field velocity induced by an incident wave, are defined relative to X = Y = 0 in the global coordinate of ANSYS AQWA and Z = 0 which is actually located at KG level and may be located above or below the free surface. Figure 5 and Figure 6 show the graphical interface for AQWA analysis platform. The interface has a number of key areas, but it is tree driven. The tree and the details are used to define the modelling requirements for the FPSO. Depending on the selection in the tree, a detailed pane, normally located below the tree shows the details of the selected object.

Input Data Tree in ANSYS AQWA
The maximum diagonal length of any element to be model is first ascertained.
Since the diagonal length of panel mesh is influenced by the maximum allowed frequency, then it should be less than 1/7 of the smallest wave length of the highest frequency analyzed. Now, considering wave length associated to a period of 4 s, this criterion leads to a maximum diagonal length slightly bigger than 4 m. The wave length for 4 s equals: where: g-gravity accelerationg; T-waveperiod = (4 s); d-waterdepth = (1000 m); In this study the maximum mesh size is 3.5 m: which leads to 3.5 m × 7 =24.5 m. The mesh size complies with the criterion since the maximum diagonal length of the mesh is slightly smaller than 4 m as shown in Figure 7.
Mass properties have to be set. The data are put to point mass and point buoyancy. These characteristics are estimated in CAD. The mass properties of the FPSO are given in Table 3. The geometry of the FPSO hull has been de-

Simplified Method of Computing Heave Response Amplitude Operator (RAO)
The dynamic pressure and acceleration of the added mass of water causing Heave motion at the bottom of the vessel when; e sin e sin That is ( )  B Therefore, the heave response amplitude operator can be written as; where 3 Q the heave dynamic amplitude factor; 3 R the frequency ratio; R ω the Natural frequency. Figure 9 shows the flowchart for the MATLAB program source code. Figure 10 displays the RAO determined for different angle of encounter-frequency. From the results, it can be seen that the maximum RAO occurs at encounter-angle of 180˚ which corresponds to Head sea as it increases significantly at wave period from zero and tends to maintain a constant peak value of one as the wave period increases. This implies that the maximum disturbance in the surge  motion is experienced in the Head sea. The Surge has natural periods within the range of periods plotted and will have resonance only in low wave periods. The value of sway is not big. Small keels on the hull will reduce the sway.

RAO for Sway Motion
Similarly, Figure 11 shows the RAO for the sway motion of the vessel. The predominant peak value for sway response occurs at beam sea as it increases significantly at wave period from zero and tends to maintenance a constant peak value of one as the wave period increases. However, the value of sway is not big. Small keels on the hull will reduce the sway. The sway motion impacts more on the ship-handling quality than the stability. The surge and sway motion RAOs show similar trends as can be seen from the comparison of Figure 10 and Figure 11. Figure 12 presents the RAO of the FPSO in Heave for different angles of encounter of the wave. The heave is reasonably well behaved and increases from zero at short periods to one in long periods. The heave has natural period within the range of periods plotted and could have resonance only in low wave periods. The value of the heave translation in 100 yr return-period is in absolute response of 1.25 m. This means that the acceleration will not cause high seasickness for the FPSO crew. However, the data from the analysis suggests a tendency for the peak value to move to slightly higher dimensionless encounter-frequency as the sea moves from the head sea to beam sea direction.

RAO for Roll Motion
The roll motion in beam sea are the most dangerous inclination. In Figure 13, Figure 11. Sway RAO for FPSO.

RAO for Pitch Motion
The pitch response amplitude operator also shows significant increase in quartering seas. In Figure 14, the inclination is less than equal to 1.

Heave RAO from the Simplified Method with MATLAB
In order to validate the results obtained in ANSYS AQWA for the Heave RAO in Head Sea a simplified method is developed in this work. The simplified virtual  predictor, in MATLAB source code, is implemented for the Heave RAO in Head Sea. The results obtained from the program is imported to excel and plotted against those obtained from ANSYS AQWA. Since different programs use different theories, the output results are slightly different. The differences in output is not due to the hydrodynamic theory but due to overly simplified assumptions and modelling errors.

Validation of Results
The MATLAB program developed in this work for the determination of the Heave RAO in Head Sea is validated by ANSYS AQWA as the FPSO responds to a sinusoidal wave. Figure 16 shows

Conclusions
The Response Amplitude Operator(s) of an FPSO unit is characterized using ANSYS AQWA hydrodynamic model. The RAOs are obtained for the six degrees From the results, the maximum RAO occurs at encounter-angle of 180˚ which corresponds to head sea. This indicates that maximum disturbance due to surge motion can be experienced in the head sea, and could have resonance only in low wave periods. The value of sway is typically small and can be further reduced by hull keels. The surge and sway motion RAOs show similar trends. Like surge, resonance due to heave occurs only in low wave periods. Also, the analysis indicates a tendency for the heave peak value to move to slightly higher dimensionless encounter-frequency as the wave moves from the head sea to beam sea direction. FPSO stability is adversely affected when Roll RAO value is well over 2.0 deg/m. In such instance, it is best to stay in mooring position and take maximum possible ballast. RAO of yaw becomes high only in quartering sea from the bow, but reduces with increasing wave period. Although a small difference in predicted heave motion occurred, it is pertinent to note that the comparison between the MATLAB program and ANSYS AQWA showed a reasonably good agreement.