A phase-field model for non-small cell lung cancer under the effects of immunotherapy

1.   Introduction

1.1.   Rational of hybriding marine technologies

With the increasing electricity demand and global warming attributed to the combustion of fossil fuels, exploitation of offshore energy sources in the form of wind, sea currents, and waves has attracted the global pursuit of renewable energy ^[1].

While a life cycle analysis (LCA) of offshore wind generation against onshore shows a 48% improvement in the project's sustainability, the cost of electricity using offshore wind is still high ^[2,3]. The generation of tidal energy is also expensive. This obstacle makes investment in tidal turbines difficult as harsh water currents decrease the design life, and at deep-sea levels, the foundation design is complicated ^[4]. Using the same foundation for wind and tidal turbines can reduce the cost of electricity by increasing power generation enabled by two different sources of energy ^[5]. However, if a support structure is not correctly designed for offshore conditions, failure will be disastrous, and the structure will be costly to replace. For this reason, foundations are an essential design consideration. Aspects to be considered while choosing and designing the foundation for a particular site include weather conditions, seabed geology, vessels and equipment required to install, and protecting the environment in accordance with local regulations.

Deployed on floating bodies or along with cables, offshore energy harvesters can convert wave, solar, tidal, ocean current, and other renewable energy sources to stable electrical energy ^[6]. Hybrid them with wind electricity generation would reduce the currently significant operations & maintenance (O & M) of a wind turbine (WT) which is around 10-25% of the total cost of the electricity, and a lower transmission cost ^[7,8].

By bringing together two marine renewable technologies with considerable synergies, the combined harnessing of offshore energies presents an excellent potential for development. This is corroborated by some recent European Union (EU) funded projects: MARINA, ORECCA, TROPOS, MERMAID and H2OCEAN ^[9]. MARINA classifies combined wave-wind systems according to the technology, water depth (shallow, transition, or deep water), or location relative to the shoreline (shoreline, nearshore, offshore). ORECCA analyses the offshore renewable energies (ORE) combined resources in Europe. Looking particularly at the combined wave-wind resource in Europe, this can be divided into three main sea basins: the Mediterranean Sea, North and Baltic Seas, and the Atlantic Ocean. TROPOS is aimed at developing a floating multi-purpose platform system for deep water ^[9]. The MERMAID project seeks to develop concepts for the next generation of offshore activities for multi-use of ocean space. It proposes new design concepts for combining offshore activities, like energy extraction, aquaculture, and platform-related transport at various ocean areas ^[10]. H2OCEAN is developing a wind-wave power open-sea platform equipped for hydrogen generation with support for multiple energy users ^[11].

1.2.   Literature review of floating offshore systems

Floating systems are the best solution for mass production in deep water sites where bottom fixed turbines are uneconomical and can be installed easier in offshore environments. Depending on how floating offshore turbines achieve their stability, they can be categorized into three main concepts, spar buoy, Tension Leg Platform (TLP), and semi-submersibles. A spar buoy foundation enhanced stability by adding ballast water. It is based on the concepts of spar buoy and Tension Leg Platform (TLP), which maintains its stability through a mooring system and buoyancy ^[12,13]. The mooring system has six degrees of freedom inhabiting motion while connecting to the seabed, making it possible to operate under different water depths ^[14,15]. This flexibility makes its market broader and more economical for investment in any range of water height ^[15]. In recent decades, there has been a great tendency to move from standing-driven foundations in the seabed to floating moored foundations. This is illustrated by shifts from fixed foundations working less than 60 m of water (e.g., Vindey Denmark 2017) to floating moored foundations in 120 m of water (e.g., Hywind Scotland 2018) and even to 1000 m of water (e.g., 2021 call for proposals from Bureau of Ocean Energy Management (BOEM) in California, USA) ^[16]. Currently, the key companies of spar buoys are Hywind/Statoil (Norway), Toda (Japan), Sway (Norway), NauticrAft (Australia) and SeaTwirl (Sweden) ^[17].

In recent years, many researchers have worked in the application of floating offshore domain in oil and gas platforms. Musial et al. ^[18] have observed that individual mooring and anchor costs are significant for single-turbine systems compared to a shared system. Fontana et al. ^[19] investigated the hydrodynamic performance and loading analysis of shared anchors for various FOWT configurations and found that shared anchors must be structurally strong enough to handle loading from unexpected directions. Anchors with a directional preference in their holding position and capacity are generally not suited for multi-line moorings, but they can be adapted by extra structural outfitting for handling mooring loads in various directions ^[20]. The type of floater is also found to have a significant influence on the anchor forces. Balakrishnan et al. ^[21] analyzed and compared the anchor forces on a semi-submersible floater system and a spar-type floater system and found that the anchor forces on the latter were less. This is because spar-type platforms have less surface area interacting with waves compared to that of a semi-submersible. Goldschmidt and Muskulus ^[22] investigated the performance of coupled-mooring systems involving 1, 5 and 10 floaters in various configurations. Semi-submersible floaters were arranged in a row and triangular and rectangular configurations separately, and the system dynamics were studied. It was found that mooring system cost reductions up to 60% and total system cost reductions up to 8% were achievable using shared moorings. However, it was noted that the displacements of the floaters were higher when the number of floaters in the system increased, which would be a problem for large wind farms. Further investigation is required to improve the behavior of floaters in larger wind farms.

2.   Concept description

In contrast to grounded foundations, aerodynamic and hydrodynamic loads create considerable motion for floating ones ^[15]. So, it is essential to define a straightforward methodology to estimate the loads on the anchor to establish a feasibility analysis ^[13]. The paper aims to provide a simplified approach for finding the maximum allowable loads necessary for a foundation which can support the selected wind and tidal turbines with spar buoy foundation. An example of offshore wind and tidal turbines supported on a floating spar is then considered based on the environmental data from NIWA. Finally, a simplified anchor sizing procedure is presented, demonstrating conservative upper bound estimates for the required suction caisson for Cook Strait. What is novel described in this paper is that it defines a method of foundation design of integrated wind and tidal structures in deep water using spar buoy floating in terms of structural modelling. The methodology used in this research for hybrid foundation design is an improved method for current method which was used before for foundation design of wind turbine.

The design will be carried out by developing spreadsheets relevant to different aspects of foundation design based on standard IEC 61400, IEA Wind TCP Task 37-May 2019, DNV 2014 and IEC 2019 ^{[23,24,25,26]}.

A single long column extends far below the water to dampen movement using water and sand as the ballast shown in Figure 1. It is loosely anchored to the sea floor. It is also loosely moored to the sea floor ^[17]. The total ballast is approximately 8, 000 tons, split between a solid part and water. The percentage split and type of ballast material are subject to further design optimization. The ballast water will typically be dosed with lye, and NaOH, giving a pH value above 10.5 ^[27].

The base case is secured to the seabed using suction anchors. The suction anchors are likely to have a maximum diameter of 7 m, corresponding to an estimated footprint of 40 m² per anchor. Suction anchors are designed such that 1 m of sand erosion/scouring is acceptable. The base case is that there is there is no need for rock dump around anchors to prevent scour. However, due to the likely presence of mobile sediments in the area, scour protection around the anchors may be required to some extent (e.g., rock dumping, mattresses). The footprint of such scour protection is expected to extend no more than 15 m out from the anchor perimeter ^[27].

The schematic structure of turbines installed with a floating foundation is shown in Figure 2.

Table 1 summarizes the features of the proposed foundation considering the total weight of wind and tidal turbines used for the hybrid system is approximately equal to the Statoil turbine. Key design parameters will receive further discussion in section 5 ^[28].

3.   Environmental parameters for design

Severe environmental conditions where the wind speed and the significant wave height have a large value are more common in deep water seas ^[14]. Therefore, tidal current modelling conducted by MetOcean Solutions Limited (MSL) was used for site optimization ^[29]. The resolution of this model in the NZ-wide grid is 0.06° (5.6 × 6.6 km). MetOcean used The Princeton Ocean Model (POM) model with a high-resolution domain over Cook Strait (0.002°; 170 × 230 m) as the best area for offshore supply, as shown in Figure 3 ^[30]. POM was used in a vertically integrated two-dimensional mode with boundaries provided by the global TPX0 solution ^[31].

The main geographical parameters of the optimized site in Cook Strait are shown in Table 2.

The most important data from NIWA's ^[34] analysis is summarized in Table 3.

Equivalent wind data will be sourced from the meteorological recording site closest to the Terawhiti taken from NASA. This data will be presented in Table 4 and Figure 4. These data are essential to estimate the wind stresses transmitted through the turbine's support structure to its foundation.

Naturally, the wind's speed constantly varies. In order to be able to predict a wind turbine 's production it is necessary to know exactly how often the wind blows how strongly. Normally, the wind is measured with an anemometer and the mean wind speed is recorded every 10 minutes. This data can be sorted into wind speed classes of 1 m/s each. The energy contained in the wind at a certain site may then be expressed by this frequency distribution. The Weibull distribution is often a good approximation for the wind speed distribution. The Weibull two-parameter distribution function often describes wind speed variability. It is considered a standard approach for evaluating local wind load probabilities because it has been found to fit a wide collection of wind data ^[39]. The Weibull shape and scale parameters are denoted by s and K, respectively. It is dimensionless and indicates how peak the site under consideration is, while K has a unit of wind speed (m/s), and it shows how windy the site is ^[40].

To calculate the Weibull parameters, the frequency percentage of Terawhiti results by Homer is used in the Weibull calculator ^[37], which results in s = 1.98 and K = 7.99 m/s.

The turbulence intensity varies with mean wind speed, which for Terawhiti is 7.10 m/s, and quantifies how much the wind varies typically within 10 minutes ^[36]. The reference turbulence intensity value (16%) may be obtained from the IEC 61400 standard ^[23].

Based on the DNV code, for height above sea level (z) less than 60 metres, Lk is 5.67z, and for z above 60 meters, Lk is 340.2 meters ^[36]. As the height above sea level is 87m, the turbulence integral length scale is 340.2 m in this case.

4.   Methods

The steps given in Figure 5 are used to evaluate the feasibility of an offshore hybrid structure in Cook Strait.

The design criteria which will be checked for the possibility of a floating design are presented in Table 5.

Design criteria specific to the selected turbines are presented in Tables 6 and 7.

There are five scenarios for an offshore foundation design exposed to wind and wave loads ^[25], summarized in Table 8.

In Table 8, the Normal Turbulence Model (NTM) relates to the normal working conditions of the turbine. The Extreme Turbulence model (ETM) is for extreme turbulence conditions. The Extreme Operating Gust (EOG) is the highest single occurrence wind load by a sudden change in the wind speed. ESS and EWH denote Extreme Sea State and Extreme Wave Height, respectively. The significant wave height H_S, used in ESS scenarios, is the average of the maximum one-third of all waves in the three hours, while the maximum wave height Hm, used in EWH scenarios, is the maximum wave height for three hours ^[36].

5.   Results

5.1.   Wind loads

In this section, the main parameters of the wind scenarios are calculated. The loads and moments will be used later in section 5.3 to identify driving combined scenarios, either E1 or E5.

Wind scenario U-1: Normal Turbulence (NTM) at Rated Wind Speed (U_R):

The mean wind speed is the rated wind speed (U_R) where the highest thrust force (Th) is expected, and the wind turbulence is modelled by the Normal Turbulence Model (NTM). This scenario is typical for the normal operation of the turbine. The standard deviation of wind speed in normal turbulence following IEC (2005) can be written as

where I_ref is the reference turbulence intensity (expected value at u = 15[m/s]), ${\sigma }_{u, NTM}$ is the standard deviation of wind speed in the normal turbulence model (NTM), and b is the model parameter defined as 5.6 according to IEC 2005. For the calculation of the maximum turbulent wind speed component (u_NTM), the time constant of the pitch control is assumed to be the same as the time period of the rotation of the rotor. In other words, it is assumed that the pitch control can follow changes in the wind speed that occur at a lower frequency than the turbine's rotational speed. Then maximum turbulent wind speed component (u_NTM) may be determined by calculating the contribution of variations in the wind speed with a higher frequency than f_1P, max to the total standard deviation of wind speed. From the Kaimal spectrum used for the wind turbulence process, its standard derivation can be calculated using Eq (2).

where ${\sigma }_{u, NTM, f > {f}_{1P}}$ is the variations in the wind speed with a higher frequency than f_1P, ${\sigma }_{u, NTM}$ is the standard deviation of wind speed in normal turbulence, L_k is the integral length scale, u_R is rated wind speed, and f_1P is the rotor's rotation frequency (1P).

The turbulent wind speed encountered in normal operation in normal turbulence conditions is found by assuming a normal distribution of the turbulent wind speed component and taking the 90% confidence level value. This is substituted into the quasi-static equation used in Eq (4). Equation (5) shows expressions for the corresponding wind moment at the mudline ^[36].

where u_NTM is the maximum turbulent wind speed component, ${\sigma }_{u, NTM, f > {f}_{1P}}$ is the variations in the wind speed with a higher frequency than f_1P, ${\rho }_{a}$ is the air density, A_R is rotor swept area, C_T is the thrust coefficient, u_R is rated wind speed, ${F}_{wind, NTM}$ and ${M}_{wind, NTM}$ are mudline load and moment in normal turbulence, S is water depth, and z_hub is hub height.

I_ref is assumed to be 16% from Table 4. Thus, ${\sigma }_{u, NTM} = 0.16(0.75\times 15+5.6) = 2.69m/s.$ 1.1 f_{1P, max} is assumed to be 0.24 Hz or f_{1P, max} = 0.218 Hz from Table 5. L_K is assumed to be 340.2 m from Table 4 and u_R = 16.5 m/s from Table 6. So, ${\sigma }_{u, NTM, f > {f}_{1P}}$ (Eq (2)) will be 0.89 m/s. This value in Eq (3) results, ${u}_{NTM} = 1.28\times 0.89 = 1.13m/s.$ Using the total wind load and moment for the Terawhiti are estimated as F_{wind, NTM} $=$ 1.68 [MN], ${M}_{wind, NTM} =$ 374.64 [MNm]

Wind scenario U-2: Extreme Turbulence (ETM) at Rated Wind Speed (U_R):

The ETM calculates the standard deviation of wind speed at the rated wind speed and, from that, the maximum wind load under normal operation in extreme turbulence conditions. The standard deviation of wind speed in ETM is given in IEC (2005) as

where Eq (6) terms are already defined (same as Eq (1)), but U_avg is the long-term average wind speed at the site. The maximum turbulent wind speed component u_ETM is determined similarly to the previous case.

The turbulent wind speed encountered in normal operation in extreme turbulence conditions, which is used for cyclic/dynamic load analysis, is found by assuming a normal distribution of the turbulent wind speed component. Unlike normal turbulence situations, the 95% confidence level value is taken. This is substituted into the quasi-static equation used in Eq (4). Equation (10) shows the expression for wind moments at the mudline ^[36].

where Eqs (8) to (10) terms are already defined (same as Eqs (3) to (5)), but U_avg from Table 4 is 7.1 m/s for Terawhiti. So, the results are estimated as:

Wind scenario U-3: Extreme operating Gust (EOG) at Rated Wind Speed (UR):

The EOG is a sudden change in the wind speed, and if it hits the rotor when it is operating at the turbine's rated wind speed, the wind turbine's pitch control has no time to alleviate the loading. IEC (2005) suggests that for simplified foundation design, 50-year extreme operating gust (EOG) is assumed to estimate the highest single occurrence wind load. The 50-year extreme wind speed, which is typically used in wind turbine design for extreme wind conditions, can be determined by standard Eqs (11) to (14) sourced from Bhattacharya ^[36].

where the number 52596 = 365.25[days/year] × 24[hours/day] × 6[10 min intervals/hour] represents the number of 10-minutes intervals in a year. The 10-year extreme wind speed is represented by

The characteristic standard deviation of wind speed can be calculated by

which c is the return period of EOG.

The extreme gust speed is then calculated at the rated wind speed from

where D is the rotor diameter = 107 m, Λ₁ = L_k∕8, with L_k being the integral length scale. Using this, the total wind load is estimated as

where Eq (15) terms are already defined (same as Eqs (4)), for wind speeds between cut-in and cut-out, thrust coefficient is estimated as recommended by Arany et al. ^[42]:

Using Eq (16), C_T = 7∕ $\stackrel{-}{u}$ = 7∕7.1 = 0.98 and rated wind speed u_R = 16.5 m/s. Using the water depth S and the hub height above sea level z_hub, the mudline bending moment (without the load factor 𝛾_L) can be estimated using Eq (17) sourced from ^[36].

The EOG wind speed is calculated using data from Table 5 and Eqs (11) to (14)

Using this, the total wind load is estimated as F_{wind, EOG}

and using the water depth S = 110 [m] and the hub height above mean sea level z_hub = 103 [m]

Wind scenario U-4: Extreme operating Gust (EOG) at the Cut-Out Wind Speed (U_out):

If the 50-year EOG hits the rotor slightly below the cut-out speed of the turbine due to the sudden change in wind speed, the turbine cannot shut down. Note that the EOG calculated at the cut-out wind speed differs from that evaluated at the rated wind speed according to IEC 2005.

Wind load caused by the EOG at the cut-out wind speed U_out is calculated taking into consideration that according to thrust coefficient expression of Frohboese and Schmuck 2010 (Eq (16)) ^[43] is no longer valid. Instead, it should be determined by assuming that the pitch control keeps the power constant. This means that the thrust force is inversely proportional to the wind speed above the rated wind speed U_R and the thrust coefficient is inversely proportional to the cube of the wind speed. So, between cut-in (U_in) and rated wind speed (U_R), Eq (16) is used for calculating the thrust coefficient and after rated wind speed, Eq (18) is used for calculating the thrust coefficient.

The EOG speed at cut-out wind speed u_{EOG, Uout} is determined from Eqs (11) to (14) by substituting a value for U_out (which is not the same at U_R). The thrust force and moment can be estimated using Eqs (19) and (20) sourced from ^[36].

where Eqs (19) & (20) terms are already defined (same as Eqs (4) & (5)). So, using Eq (14) ${u}_{{EOG, u}_{out}} =$ 7.1 m/s and as a result:

Table 9 summarizes the important parameters and presents the wind loads for the different wind scenarios to determine the turbulent wind speed component and, through that, the thrust force and overturning moment. The wind loads on the structure are independent of the substructure diameter; therefore, the wind loads can be evaluated before the pile and substructure design are available.

Applying an environmental load factor of 𝛾_L = 1.35 as recommended by DNV 2014 ^[26] and IEC 2019 ^[24], the total wind moment is 862.65 [MNm] for U3. It is found that the EOG at U_R (U-3) gives the highest load and overturning moment.

5.2.   Wave loads

The wave conditions recommended by Bhattacharya ^[36] for calculating critical wave loads acting on the substructure are:

(W-1) One-year extreme sea state (ESS). A wave with height equal to the 1-year significant wave height H_{S, 1}.

(W-2) One-year extreme wave height (EWH). A wave with a height equal to the 1-year maximum wave height H_{m, 1}.

(W-3) 50-year ESS. A wave with a height equal to the 50-year significant wave height H_{S, 50}.

(W-4) 50-year EWH. A wave with a height equal to the 50-year maximum wave height H_{m, 50}.

Bhattacharya ^[36] states that the one-year ESS and EWH are used as a conservative overestimation of the normal wave height (_NWH) prescribed in IEC 2019 ^[24]. In relation to the ESS, it is important to note that the significant wave height and the maximum wave height have different meanings. The significant wave height H_S is the average of the highest one-third of all waves in the three-hour sea state, while the maximum wave height Hm is the single highest wave in the same three-hour sea state.

The highest wind load is expected to be caused by scenario U-3, and the highest wave load is due to scenario W-4. In practice, the 50-year extreme wind load and the 50-year extreme wave load have a negligible probability of occurring at the same time and the DNV.

2014 code ^[26] also doesn't require these extreme load cases to be evaluated together. The designer has to find the most severe event with a 50-year return period based on the joint probability of wind and wave loading. Therefore, for the ULS analysis, two combinations are suggested by Arany et al. 2017 ^[42]:

(1) The ETM wind load at rated wind speed combined with the 50-year EWH-the

combination of wind scenario U-2 and wave scenario W-4. This will provide higher

loads in deeper water with higher waves.

(2) The 50-year EOG wind load combined with the one-year maximum wave height.

This will provide higher loads in shallow water in sheltered locations where wind load dominates. These scenarios are somewhat more conservative than those required by standards and can be adapted for simplified analysis. From the point of view of SLS and FLS, the single largest loading on the foundation is not representative because the structure is expected to experience this level of loading only once throughout its lifetime.

In simplified load calculation methodologies, simple linear waves are assumed to determine the wave loading. Higher-order theories like Stokes waves or Dean's stream function theory would provide better estimates, especially in shallow waters. However, the linear theory allows for simpler load calculation, and its application is justified for foundation design loads.

The circular substructure area A_S is also calculated from this diameter. The methodology used here builds on linear (Airy) wave theory, which gives the surface elevation 𝜂, horizontal particle velocity w, and the horizontal particle acceleration $ẇ$ as

where H_m is maximum wave height, T_s is peak wave period, S is mean water depth, t is time, and x is the horizontal coordinate in the along-wind direction. The wave number k is obtained from the dispersion relation

The force on a unit length strip of the substructure is the sum of the drag force F_D and the inertia force F_I:

where w horizontal particle velocity is a function of z coordinate and t time, the C_D is the drag coefficient, C_m is the inertia coefficient, 𝜌_w is the density of seawater, D_S and A_S are the diameter and circular substructure area. The total horizontal force and bending moment at the mudline are then given by integration as

The peak load of the drag and inertia loads occur at different time instants; therefore, the maxima are evaluated separately. The maximum inertia load occurs at the time instant t = 0 when 𝜂 = 0, and the maximum of the drag load occurs when t = T_S/4 and 𝜂 = H_m∕2.

The maximum drag load & moment ${F}_{D, max}$ & ${M}_{D, max}$ is then obtained by carrying out the integrations:

The maximum inertia load & moment ${F}_{I, max}$ & ${M}_{I, max}$ are obtained by:

In the simplified method for obtaining foundation loads, it can be conservatively assumed that the sum of the maxima of drag and inertia loads is the design wave load. This assumption is conservative because the drag and inertia load maxima occur at different time instants. All wave scenarios (W-1)-(W-4) are evaluated with the same procedure, using different values of wave height H and wave period T ^[36].

In this step, a simplified estimation of the extreme wave height and the corresponding wave period for a given site is explained, which involves the following sub-steps:

1. Obtain the relevant significant wave height H_S from a reliable source. The 1-year equivalents are calculated following DNV 2014 ^[26] from the 50-year significant wave height according to

Calculate the peak wave period corresponding to the significant wave height, T_S using:

2. Calculate the number of waves in a 3-hour period, N: Typically, significant wave heights are given for a 3-hour period. In other words, this means that the significant wave height is calculated as the mean of the highest 1/3 of all waves. Therefore, many different wave heights occur within this 3-hour period, and the highest occurring wave height is called the maximum wave height H_m. To find this, one needs to know the number of waves in the 3-hour period because the more waves there are, the higher the chance of higher waves occurring.

3. Calculate the maximum wave height, H_m by

4. Calculate the peak wave period corresponding to the maximum wave height, Tm. The same formulae can be used as in step 2.

The 50-year significant wave height, ${H}_{S, 50}$ at Terawhiti can be found from NIWA ^[44]. ERA-40 Wave Atlas ^[35] provides an approximate value of 15 m shown in Table 11. So, according to Eq (36), H_{S, 1} is 12 m. Going to step 2, T_{S, 1} can be 12.2s. Equation (39) gives N = 885.24, and as a result, H_{m, 1} = 22.10 m and T_{m, 1} = 16.66 s. Also, H_{m, 50} is 27.62m according to Eq (39). The wave periods for 4 scenarios are:

(W-1): ${T}_{S, 1} = 11.1\sqrt{\frac{{H}_{S, 1}}{g}}$ = $11.1\sqrt{\frac{12}{9.8}}$ = 12.2 s

(W-2): ${T}_{m, 1} = 11.1\sqrt{\frac{{H}_{m, 1}}{g}}$ = $11.1\sqrt{\frac{22.10}{9.8}}$ = 16.66 s

(W-3): ${T}_{S, 50} = 11.1\sqrt{\frac{{H}_{S, 50}}{g}}$ = $11.1\sqrt{\frac{15}{9.8}}$ = 13.73 s

(W-4): ${T}_{m, 50} = 11.1\sqrt{\frac{{H}_{m, 50}}{g}}$ = $11.1\sqrt{\frac{27.62}{9.8}}$ = 18.63 s

The wave heights and wave periods are summarized for all wave scenarios (W-1) to (W-4) in Table 10 according to the explanations presented.

The wave loads (Eqs (26) and (27)) are calculated only for the most severe wave scenarios used for Load Cases E-2 and E-3, i.e., wave scenario (W-2) and (W-4), the 1-year and 50-year extreme wave heights (EWHs), using the definitions presented in Table 8 where ${\rho }_{w} = 1030$ kg/m³ and S = 110 m (from Table 2)

This section can be modelled using an equivalent diameter of D_D = 11.33 m for drag load calculations and D_I = 12.89 m for inertia load calculations.

C_D = 0.5 and C_m = 2 (based on DNV ^[26])

Hm, Ts, 𝜂, k for W-1 scenario: H_m = 12 m, T_s = 12.2 s, 𝜂 = 0 for inertia loads and 𝜂 = 6 m for drag loads, k = 0.034 m^-1 (based on Eq (24))

Using Eqs (28) to (35) for the W-1 scenario:

The results are:

Hm, Ts, 𝜂, k for W-2 scenario: H_m = 22.10 m, T_s = 12.2 s, 𝜂 = 0 for inertia loads and 𝜂 = 11.05 m for drag loads, k = 0.034m^-1(based on Eq (24))

Using Eqs (28) to (35) for the W-2 scenario:

The results are:

F_{wave, W-2} = 2.8+1.05 = 3.85 MN, M_{wave, W-2} = 643.2+86.7 = 729.9 MNm

Hm, Ts, 𝜂, k for W-3 scenario: H_m = 15 m, T_s = 13.73 s, 𝜂 = 0 for inertia loads and 𝜂 = 7.5 m for drag loads, k = 0.029 m^-1(based on Eq (24))

Using Eqs (28) to (35) for the W-3 scenario:

The results are:

F_{wave, W-3} = 1.2 MN, M_{wave, W-3} = 15279.4 MNm

Hm, Ts, 𝜂, k for W-4 scenario: H_m = 27.62 m, T_s = 13.73 s, 𝜂 = 0 for inertia loads and 𝜂 = 13.81 m for drag loads, k = 0.029 m^-1(based on Eq (24))

Using Eqs (28) to (35) for the W-4 scenario:

The results are:

5.3.   Anchor load combinations for Ultimate Limit State (ULS)

The current load is calculated as ^[36]

The loads under the combined actions of wind and waves were calculated using different scenarios in Table 8 and adding the current load as below:

As expected, the wave load dominates, and the scenario with the combination of the 1-year EWH and the 50-year extreme operating gust (EOG) or E3 produces the ULS load. It should be noted here that this load is conservative for anchor design, as the load that acts on the anchor is reduced by the weight of the suspended section of the mooring line, the friction on the horizontal section (Touch Down Zone) of the mooring line, the soil reaction on the inverse catenary-shaped forerunner in the soil, and the weight of the forerunner. The vertical load acts on the spar at the instant when the surface elevation at the spar is at its highest point (wave crest). In contrast, the horizontal load is dominated by the inertia load, highest when the surface elevation is at the mean water level. Therefore, the ultimate load is taken as the horizontal load as calculated above ^[36].

5.4.   Minimum caisson dimensions for Terawhiti with sandy soils

The diameter of the caisson D and the embedment depth L are the two main independent parameters that govern the holding capacity of the caisson for a given soil profile. Soil type is Soft/medium sand with the angle of internal friction of ϕ´ = 30°, effective unit weight of 𝛾′ = 9 and mooring chain friction on sand $\mu = 0.25$ ^[45].

The holding capacity of suction caissons is typically determined in terms of an envelope based on the horizontal and vertical load components at the anchor as:

where

H_m is the horizontal capacity, and V_m is the vertical capacity. On the other hand, H_u and V_u are the applied load. FP is the failure criterion, and the maximum value can be 1 (limiting condition). The horizontal capacity in the sand can be calculated as follows:

where A_b is the effective unit bearing area of the forerunner (equals the diameter of the rope or wire, and 2.5-2.6 times the bar diameter for a chain), Nq is the bearing capacity factor as defined by:

with 𝜙 being the internal angle of friction of the soil.

The vertical capacity in the sand will be

where

$K~~ tan\delta :$ the factor that only appears together.

K: the effective stress factor used to calculate the effective horizontal stress as constant times the effective vertical stress (𝜎_H = K𝜎_V′), 𝛿 is the mobilized angle of friction between the caisson wall and the soil, e represents the external, and i is the internal circumference of the caisson.

Z: D∕[4K tan 𝛿] with e and i referring to external and internal values, respectively.

The anchor padeye tensions and angles can be determined by (look at Figure 6):

where

T_a tension at the anchor padeye.

T_m tension at the mudline.

𝜃_a angle of the tension at the anchor padeye to horizontal.

𝜃_m angle of the tension at the mudline to horizontal.

z_a depth of the anchor padeye below mudline.

𝜇 friction coefficient between the forerunner (chain, rope, or wire) and the soil.

Q_av average soil resistance between the mudline and the padeye is calculated as:

where N_c is the bearing capacity factor and between 9 and 14, s_u(z) is the distribution of the undrained shear strength with depth.

For sand, z_a/L = 2/3 ^[36]. Statoil 2015 ^[27] advised each anchor's maximum diameter of 7 m. As this design is a hybrid design of Hywind statoil by adding tidal turbines and considering wind and tidal conditions of New Zealand, the caisson dimensions will be according to ^[13]:

6.   Discussion and conclusions

The aim of this paper is to develop a unique method for structural design of hybrid offshore systems in deep waters using spar buoy floating. The method enables to find required dimensions tolerating against wind and wave loads. Terawhiti in Cook Strait was selected as a case study because according to MetOcean model, this location can generate maximum electricity and suitable for this type of foundation. Based on the selection of turbines and Terawhiti as facility site and by including the wind and wave loads, the design provides the stability of structure.

The wind and water produce aerodynamic and hydrodynamic loads (thrust and drag) on the structure which depends on the operational speed of turbines. But, to know the acceptability of foundation design, it is necessary to combine wind and wave loads in ULS design and calculate maximum loads and find driven scenario. Based on different scenarios of foundation design, next step will be calculating maximum loads, and find driven scenario. Then, the minimum required sizes of caisson based on the maximum load of the driven scenario for the Hybrid floating offshore structure can be estimated. These dimensions satisfy the failure criterion.

The results, as presented and discussed in Section 5, enable several conclusions:

● The highest wind load for ULS is determined from the 50-year extreme operating gust (EOG); this is wind scenario U-3 which is 3 MN. Total wind loads for U-2, U-1 and U-4 are 1.86, 1.68 and 0.69 MN respectively.

● The maximum drag and inertia loads occurs for wave scenario (W-2) with the 1-year EWH which are 2.8 and 1.05 MN respectively or totally 3.85 MN. This values for W-1, W-3 and W-4 are 0.94, 1.2 and 2.14 MN respectively.

● The maximum load occurs in extreme wind load scenario when wind and wave act in same direction or E-3. This equals to combined values of U-3 and W-2 which is 7.32 MN. It is higher than combined wind and wave loads of scenarios E-1, E-2, E-4 and E-5 (as defined in Table 8) which are 3.09, 4.47, 3.3 and 4.47 MN respectively. Based on mentioned maximum load, the required diameter of the caisson and the embedment depth are 6.9 and 22.07 meter respectively.

The future work in this aspect can be estimating the cost analysis in terms of comparing the proposed hybrid system with a wind floating and with a floating tidal system and analyzing how much using the same structure for both wind, and tidal turbines can reduce the cost of electricity generation rather than having a foundation design just for wind and also just for tidal turbines.

Funding

This research did not receive any specific grant from public, commercial, or not-for-profit funding agencies.

Data availability statement

Data is available on request from the authors.

Conflicts of interest

The authors declare no conflict of interest.