A Fourier Series Approximation for Deep-water Waves

Article information

J. Ocean Eng. Technol. 2022;36(2):101-107
Publication date (electronic) : 2022 April 12
doi : https://doi.org/10.26748/KSOE.2021.092
1Engineer, Offshore structure design department, Daewoo Shipbuilding & Marine Engineering co., LTD, Geoje, Korea
Corresponding author JangRyong Shin: +82-55-735-5117, jrshin@dsme.co.kr
Received 2021 December 1; Revised 2022 February 15; Accepted 2022 March 21.

Abstract

Dean (1965) proposed the use of the root mean square error (RMSE) in the dynamic free surface boundary condition (DFSBC) and kinematic free-surface boundary condition (KFSBC) as an error evaluation criterion for wave theories. There are well known wave theories with RMSE more than 1%, such as Airy theory, Stokes theory, Dean’s stream function theory, Fenton’s theory, and trochodial theory for deep-water waves. However, none of them can be applied for deep-water breaking waves. The purpose of this study is to provide a closed-form solution for deep-water waves with RMSE less than 1% even for breaking waves. This study is based on a previous study (Shin, 2016), and all flow fields were simplified for deep-water waves. For a closed-form solution, all Fourier series coefficients and all related parameters are presented with Newton’s polynomials, which were determined by curve fitting data (Shin, 2016). For verification, a wave in Miche’s limit was calculated, and, the profiles, velocities, and the accelerations were compared with those of 5th-order Stokes theory. The results give greater velocities and acceleration than 5th-order Stokes theory, and the wavelength depends on the wave height. The results satisfy the Laplace equation, bottom boundary condition (BBC), and KFSBC, while Stokes theory satisfies only the Laplace equation and BBC. RMSE in DFSBC less than 7.25×10−2% was obtained. The series order of the proposed method is three, but the series order of 5th-order Stokes theory is five. Nevertheless, this study provides less RMSE than 5th-order Stokes theory. As a result, the method is suitable for offshore structural design.

1. Introduction

The Brazilian multinational petroleum corporation Petrobras and partners are developing the Buzios oil field, which is approximately 210 km offshore of Brazil. The Floating production storage and offloading (FPSO) is spread-moored in a maximum water depth of 2,030 m and has facilities to receive oil from sub-sea wells. It also has production plant facilities to process fluids, stabilize them, and separate produced water and natural gas, which is re-injected into a dedicated reservoir. Processed liquids are metered, stored in the FPSO cargo storage tanks, and offloaded to export tankers. The design life of the FPSO is 30 years.

The relative water depth is defined as kh = ω2h/g, where k = 2π/L is the wave number, ω = 2π/T is the angular frequency, L is the wavelength, T is the wave period, h is water depth, and g is gravity. According to DNV (2010a), it is normally not necessary to investigate wave periods longer than 18 s. Therefore, the relative water depth is greater than 25 for the project. When the relative water depth is greater than 2, deep-water wave theories are applicable (Chakrabarti, 1987; DNV, 2010b; Shin, 2019).

Well-known wave theories include Airy theory, Stokes theory, Dean’s stream function theory, Fenton’s theory, and trochodial (Gerstner) theory for deep-water waves in offshore structure design. Trochodial theory is an exact solution of the Euler equation with vorticity. The first rotational solution was described by Gerstner in 1802 and was independently rediscovered later by Rankine (1863). A mathematical analysis of trochodial theory was performed by Constantin (Henry, 2008). The wavelength is independent of the trochoidal wave’s height, unlike in Stokes’ wave theory and observations. The trajectories of a water particle are closed circles, in contrast with the usual experimental observation of Stokes drift associated with wave motion. Therefore, trochodial theory is of limited use for offshore structure design.

Airy theory, Stokes theory, Dean’s stream function theory, and Fenton’s theory are irrotational wave theories, unlike Trochodial theory. The wavelength is also independent of Airy wave’s height, and Airy theory is applicable for kH ≤1/25 (Chakrabarti, 1987), where H is the wave height. Therefore, Airy theory is unsuitable for describing waves near the Miche limit (DNV, 2010b); i.e., H/Lo = 0.14, where Lo is the wavelength calculated by Airy theory.

The first nonlinear solution was reported by Stokes (1847), who calculated waves in infinite/finite depth. De (1955) and Skjelbreia et al. (1960) obtained the 5th-order solutions for waves in finite depth. Schwartz (1974) and Cokelet (1977) obtained solutions of higher-order expansions. Stokes theory is based on a perturbation method and provides a closed-form solution, but it does not satisfy two boundary conditions on the free surface. The 5th-order theory (De, 1955; Skjelbreia and Hendrickson, 1960) is applicable for all deep waves (Chakrabarti, 1987), but the error of the theory is more than 1% for waves near the Miche limit (DNV, 2010b). Stokes theory gives non-conservative results for offshore structure design.

Based on the use of truncated Fourier expansions for field quantities, Chappelear (1961) and Dean (1965) developed a Fourier approximation method. By choosing expansions to satisfy the Laplace equation and the BBC, the problem is reduced to solving a set of nonlinear equations for each of the Fourier coefficients. Dean’s stream function theory (Dean, 1965) is a purely numerical procedure (DNV, 2010b). It requires order higher than 11 for describing waves near the Miche limit (DNV, 2010b). Nevertheless, its error is more than 1% (DNV, 2010b).

Rienecker and Fenton (1981) adopted the stream function expansion (Dean, 1965) and directly solved for the coefficients by Newton’s method. This method was further simplified by Fenton (1988). The major simplification is that all the partial derivatives are obtained numerically. The most error results from the reason that the wave profile and the Fourier coefficients are simultaneously calculated in Fenton’s theory (Fenton, 1988). To reduce the error, Fenton’s theory requires more series order than 64 (Rienecker and Fenton, 1981). However, since the stream function expansions contain hyperbolic functions, neither of these stream function approaches can be applied for deep-water waves.

Shin (2016) used a Fourier approximation and directly solved the coefficients by Newton’s method, but there are some differences between Shin’s (2016) and Fenton’s theories (Rienecker et al., 1981; Fenton 1988). While Fenton’s theory (Rienecker and Fenton, 1981; Fenton 1988) adopted the stream function (Dean, 1965), Shin (2016) adopted the velocity potential. Fenton’s theory adopted the moving coordinate system introduced by Dean (1965), but Shin (2016) adopted a dimensionless coordinate system. The moving coordinate system requires partial derivatives with regard to the wavelength since the wavelength is coupled to the Fourier coefficients, but the dimensionless coordinate system does not. Therefore, Newton’s method was further simplified (Shin, 2016). While the wave profile and the Fourier coefficients are simultaneously calculated in Fenton’s theory, they are independently calculated by Shin (2016). As a result, the required series order and total error are dramatically reduced (Shin, 2016).

These problems were discussed with Petrobas for the project. It was agreed that Shin’s (2016) wave theory would be applied for calculating wave loads on hard piping riser structures because the theory has less error and provides greater wave loads than 5th-order Stokes theory (Chakrabarti, 1987; De, 1955; Stokes, 1847; Stokes, 1880; DNV, 2010b; Skjelbreia and Hendrickson, 1960). For the application, the results from Shin (2016) were further simplified in this study. The Fourier series is represented for deep-water waves, and all Fourier coefficients and all related parameters are represented with Newton’s polynomials. As a result, this study provides a closed-form solution for deep-water waves. A wave in Miche’s limit was calculated and compared to that of 5th-order Stokes theory. The profiles, velocities, and accelerations were also compared to each other for verification.

2. Coordinate Systems

For the description of a progressive water wave, two coordinate systems are adopted. One is the conventional coordinate system (t,x,y) shown in Fig. 1. The origin is located on the still water line. The x-axis is in the direction of wave propagation, the y-axis points upwards, and t is time. The fluid domain is bounded by a free surface y = η(t, x).

Fig. 1.

A conventional coordinate system for a progressive wave (SWL: still water line)

The other coordinate system (β,α) is the dimensionless stationary frame shown in Fig. 2. The origin is located at the point under the crest on the reference line, which is the horizontal line passing through two points at ηo = η(±90°) on the free surface and two phases of β = ±90°. Therefore, the wave profile is a fixed, periodic, even function in the system. The horizontal axis is the phase, β = kxωt (πβπ). The vertical axis is the dimensionless elevation from the reference line, α = k(yη0) in αζ, where ζ = k(ηη0) is the dimensionless free surface elevation from the reference line.

Fig. 2.

A dimensionless coordinate system for a progressive wave

The dimensionless coordinate system provides several advantages besides the one mentioned in previous section. One independent variable is removed in the coordinate system. Therefore, the Korteweg-de Vries equation and the Benjamin-Bona equation are presented with ordinary differential equations in the coordinate system. The KFSBC is presented with an exact differential equation in the coordinate system.

One parameter is removed in this coordinate system. While deep-water waves are described with two parameters (i.e., wave height and wave period) in a conventional coordinate system or in a moving coordinate system (Dean, 1965), they are described with only one parameter, the linear steepness θ = (ω2H)/g, in the dimensionless coordinate system. Furthermore, the linear steepness is in the range of 0≤θ<1 since there is a wave-height limitation. Therefore, it is possible to convert the result from Shin (2016) to a closed-form solution.

3. Fourier Series Approximation

For deep-water waves, the solution is represented as follows. By using Fourier series approximation, the velocity potential is presented as follows:

(1) ϕ=ωk2n=1Nanenαsinnβ
where an is a Fourier coefficient, and N is the required Fourier series order. The velocity potential satisfies the Laplace equation and the BBC. After differentiating Eq. (1) with respect to x and y, the horizontal water-particle velocity is:
(2) u(β,α)=cn=1Nnanenαcosnβ
where c=defω/kis the wave celerity. The vertical water-particle velocity is:
(3) v(β,α)=cn=1Nnanenαsinnβ

The water-particle accelerations in the horizontal direction and the vertical direction are represented by:

(4) ut=ω2kn=1Nn2anenαsinnβ
(5) vt=-ω2kn=1Nn2anenαcosnβ

From the KFSBC, the wave profile is calculated as follows:

(6) ζ=n=1Nan{enζcosnβ-cosnπ2}

Using Bernoulli’s equation, the pressure is calculated as follows:

(7) p(β,α)ρc2=U(β,α)-U0-12{U2(β,α)+V2(β,α)}+12(U02+V02)-Sαθ
where ρ is the water density, and U and V are dimensionless horizontal and vertical velocities defined by U=defu/c and V=defv/c. U0 =U(π/2,0) and V0 =V(π/2,0) are velocities at the phase of β90° on the free surface. By applying the DFSBC in Eq. (7), the other wave profile is calculated as follows:
(8) ζSθ=U(β,ζ)-U0-12{U2(β,ζ)+V2(β,ζ)}+12(U02+V02)

The linear steepness θ = ω2H/g is a known constant for a particular wave. Dimensionless wave height (steepness) is defined by S = kH (H: wave height), which provides the dispersion relation because the wave number is calculated as k =S/H.

4. Wave Profile

An implicit function f(β,ζ)= 0 can be considered as an equation with respect to the dependent variable ζ because the independent variable β is known. Therefore, using Newton’s method, the explicit representation of the function is presented:

(9) ζ(β)=limnζn+1(β)
where n stands for a step of Newton’s method, and ζn + 1 is:
(10) ζn+1(β)=ζn(β)-f(β,ζn)f(β,ζn)
where f(β,ζ) is the error in the KFSBC defined by Eq. (6) as follows:
(11) f(β,ζ)=-ζ+n=1Nan{enζcosnβ-cosnπ2}
f′(β,ζ) is the first-order partial derivative with respect to the variable ζ. Therefore, we have:
(12) f(β,ζ)=-1+n=1Nnanenζcosnβ

Because |ζ|<1 for all waves, the first-step solution ζ1 is calculated with the power series expansion of f in ζ, i.e., f(β,ζ)=m=0f(m)(β,0)m!ζm while ignoring the higher-order terms after the first order. Because f(β,ζ)= 0 for all phases, the power series expansion is a linear equation with respect to ζ, i.e., f(β,0)+ζf’(β,0)= 0. Then we have:

(13) ζ1=n=1Nan{cosnβ-cosnπ2}1-n=1Nnancosnβ

When θ≤0.98992, Newton’s method rapidly and absolutely converges to a complete solution. Therefore, θmax=0.98992 is the limit in this study. Since the limit is greater than the Miche limit (DNV, 2010b; Shin, 2019; Chakrabarti, 1987), which is 0.28π, the proposed method can be used for all deep-water waves.

5. Verification

Fourier coefficients were determined so that the two profiles meet at four phases: β = 0°, 30°, 90°, and 180°. Therefore, Eq. (6) satisfies the KFSBC for all phases and satisfies the DFSBC at only the four phases. Eq. (8) satisfies the DFSBC for all phases and satisfies the KFSBC at only the four phases, but the difference of the two profiles is very small.

The other wave profile can be calculated by applying the same method presented in section 4 to Eq. (8) or by substituting the result of Eq. (9) into the right side of Eq. (8). A wave with a period of 6 s and height of 7.87 m was calculated by the proposed method and by 5th-order Stokes theory. The wave is in the Miche limit. Wave profiles are compared in Fig. 3. The red solid line is the profile from KFSBC and was calculated with Eq. (6). The green dotted line is the profile from DFSBC and was calculated with Eq. (8). It is extremely difficult to distinguish the two curves. The purple dashed line is the profile calculated with 5th-order Stokes theory. Bernoulli’s constants were also calculated on the free surface. The blue solid line is Bernoulli’s constant from the proposed method, and the blue dashed line is Bernoulli’s constant calculated with 5th-order Stokes theory.

Fig. 3.

Wave profiles and Bernoulli’s constants on the free surface

The constants represent error in the DFSBC. If there is no error in the DFSBC, the constant is presented with a horizontal straight line like the blue solid line from Bernoulli’s principle. Therefore, the blue dashed line shows that Stokes theory has some error in the DFSBC. The blue solid line has a total RMSE (Chakrabarti, 1987) of 7.25×10−2%, but the blue dashed line has a total RMSE of 0.37%. In addition to the error, Stokes theory has 1.42% error in KFSBC, unlike the proposed method. Therefore, the total error of Stokes theory is 1.79%.

The wavelength is 64.8 m in the proposed method and 64.52 m in Stokes theory. L0= 56.21 m is obtained by Airy theory. The wavelength calculated by the proposed method depends on the wave height, unlike Airy theory or trochodial theory.

Water particle velocities on the free surface are presented in Fig. 4. The two solid lines were calculated by the proposed method. The two dashed lines were calculated by 5th-order Stokes theory. The two symmetric curves represent horizontal velocities. There is a big difference in the horizontal velocity under the crest, unlike the wave profiles shown in Fig. 3. As a result, there is more difference in the drag force since the force is proportional to the square of the velocities.

Fig. 4.

Water particle velocities on the free surface

Water particle accelerations on the free surface are presented in Fig. 5. The two solid lines were calculated by the proposed method. The two dashed lines were calculated by 5th-order Stokes theory. The two symmetric curves represent vertical accelerations. There are big differences in the acceleration. As a result, there are big differences in the inertia force since the force is proportional to the acceleration.

Fig. 5.

Water particle accelerations on the free surface

Figs. 45 show that Stokes theory gives non-conservative results for offshore structure design. The analysis results are summarized in Table 1.

Analysis results of the wave with period of 6 s and height of 7.87 m.

Even though the series order of the proposed method is less than in Stokes theory’s, the error of the proposed method is lower.

6. Solution Method and Curve Fitting Results

This study is based on a Fourier approximation. There are differences in numerical methods to determine the Fourier coefficients and the related parameters between Shin’s (2016) method and the other Fourier approximations (Dean, 1965; Chappeler, 1961; Chaplin, 1979; Rienecker and Fenton, 1981; Fenton, 1988). The result from Shin (2016) is summarized for a closed-form solution in the following.

The solution contains 5 unknown constants: 3 Fourier coefficients an, the steepness S, and the reference line − 0. Two wave profiles from the KFSBC and the DFSBC were presented with implicit functions in Eqs. (6) and (8). Considering the dimensionless coordinate system, the two wave profiles should be even functions. The Fourier series of a periodic even function is presented as f(x)=n=0ancosnxin-πxπ. The coefficients are determined using the orthogonality of trigonometric functions in general.

When f(x) is an implicit function like Eqs. (6) and (8), the method is impossible. A method to solve the problem is to convert the series to a set of algebraic equations that are obtained when we calculate the series at some phases instead of all phases by replacing the infinite series with a truncated series: n=0Nancosnxm=f(xm) for m=1, ... , N.

Because cosnxm and f(xm) are known, we have N algebraic equations for calculating the coefficients, an. Therefore, we have a set of N equations such that Eqs. (6) and (8) are equal to each other at N phases. In this study, N=3 was considered. Additionally, we have two equations: the wave height condition and water depth condition. As a result, we have 5 equations to determine the unknown constants.

Assuming that ζc = ζ(0o), ζ30 = ζ(30o), and ζt = ζ(180o) are known values and applying the idea above to Eq. (6), we have a set of linear algebraic equations with regard to the three coefficients. From the set of equations, the Fourier coefficient a2 is determined as follows:

(14) a2=A22A21
where
(15) A22=defζce-ζc-2ζ30e-ζ303+2ζ30e-(ζ30+2ζt-2ζc)3+ζte-(3ζt-2ζc)
and
(16) A21=def13e(ζ30-2ζt+2ζc)+23e-(ζ30+2ζt-2ζc)+2e-2(ζt-ζc)coshζt-13eζ30-23e-ζ30+2coshζc

The Fourier coefficient a1 is determined as follows:

(17) a1=23{ζ30e-ζ30-a2(12eζ30+e-ζ30)}

The Fourier coefficient a3 is determined as follows:

(18) a3=-2ζ30e-(ζ30+2ζt)3-ζte-3ζt+(13e(ζ30-2ζt)+23e-(ζ30+2ζt)+2e-2ζtcoshζt)a2

Applying the same idea to Eq. (8) and substituting Eqs. (14), (17), and (18) in Eq. (8), we have a set of three nonlinear equations with regard to ζc, ζ30, and ζt. From the wave height condition, we have the following:

(19) S=ζc-ζt

We also have an equation from the water depth condition, in which the still water line is the average of the wave profile during a wave period or over a wavelength. Because the horizontal axis of the conventional coordinate system is on the still water line, the reference line is represented as:

(20) kη0=-12π-ππζdβ

Therefore, we have five nonlinear equations, which were solved by Newton’s method by Shin (2016). All the unknown constants are functions of the linear steepness θ, which is a known constant for a particular wave in the range of 0≤θ<1. They were calculated with a highly dense interval of the linear steepness, and the results were shown by Shin (2016). The constants are represented with Newton’s polynomials in this study. For the curve fitting, the divided differences method (Gerald and Wheatley, 2006) is considered. For θ≤0.9, the steepness is presented as follows:

(21) S=0.99751F0-0.0729F1-0.21233333F2+0.11125F3+0.07333333F4-0.14722222F5+0.176587302F6-0.12400794F7+0.270061728F8
where Fi is a Newton basis polynomial defined as Fi=defn=0i(θ-0.1n). For 0.9<θ<1, the steepness is:
(22) S=0.776933+0.6819G0-0.11G1+1.6666667G2-50G3+2333.33333G4-75000G5+2242063.492G6-52579365.1G7+1948302469G8
where Gi=n=0i{θ-(0.9+0.01n)}. The wave number is determined by Eq. (21) or (22) and k =S/H. For θ≤0.9, the ratio of the trough depth to crest height is presented as follows:
(23) λ=def-ζtζc=1-0.95155346F0+0.449837676F1-0.12765248F2-0.09413857F3+0.312588656F4-0.90252346F5+1.622263968F6-3.25526214F7+1.865838796F8

For 0.9<θ<1, the ratio of the trough depth to crest height is presented as follows:

(24) λ=0.385658-0.5997714G0-1.12845501G1-10.1573878G2-177.266276G3+470.1400869G4-123978.831G4-1305636.873G6-130109452G7-19577000000G8

For θ≤0.9, the ratio of the elevation ζ30 at a phase of 30° to the crest height is presented as follows:

(25) λ30=def-ζ30ζc=32-0.0639043F0-0.06427581F1-0.04834541F2-0.06425683F3-0.0302724F4-0.19240911F5+0.003370865F6-1.18380101F7-3.46302865F8

For 0.9<θ<1, the ratio of the elevation ζ30 at phase of 30° to the crest height is presented as follows:

(26) λ30=0.697103-0.67038918G0-2.80988095G1-21.3795326G2-191.85558G3-4414.7284G4-40908.9389G5-2822028.06G6-99414564.7G7-37945000000G8

By using the definition in Eq. (23) and Eq. (19), we have:

(27) ζc=S1+λ

Therefore, by substituting Eqs. (21) or (22) and (23) or (24) in Eq. (27), the crest elevation is determined. By using the definition in Eq. (23), we have:

(28) ζt=-λζc

By substituting Eqs. (23) or (24) and (27) in Eq. (28), the trough elevation is determined. By using the definition in Eq. (25), we have:

(29) ζ30=λ30ζc

By substituting Eqs. (25) or (26) and (27) in Eq. (29), the wave elevation at a phase of β = 30° is determined. The reference line is numerically calculated by applying Eq. (9) in Eq. (20). The reference line was also calculated with a highly dense interval of the linear steepness. The result is represented with Newton’s polynomials as follows:

(30) kη0=-(ζt+Ikπ)

For θ≤0.9, the parameter Ik is presented as:

(31) Ik=1.5277031F0-0.502451F1-0.3040983F2+0.208770833F3-0.03169167F4-0.214652778F5-0.102063492F6-0.6900298F7-4.066964286F8

For 0.9<θ<1, the parameter Ik is presented as:

(32) Ik=0.89942883+0.123682G0-2.9221G1-19.328333G2-216.75G3-3270G4-67625G5-2072023.81G6-284632937G7-21853395062G8

The dimensionless elevation is determined by Eq. (9), the reference line 0 is determined by Eq. (30), wave number is determined by Eq. (21) or Eq. (22) and then, substituting the three results in the definition of the dimensionless elevation, ζ = k(ηη0), the wave profile is calculated as:

(33) η=ζ+kη0k

ζc, ζ30, and ζt are merely parameters to calculate the Fourier coefficients. The wave profile is calculated with Eqs. (9) and (33). This is one of the major differences between Shin’s (2016) theory and Fenton’s theory (Rienecker et al., 1981; Fenton 1988).

7. Conclusions

The purpose of this study was to provide a closed-form solution for deep-water waves with error less than 1%. For this purpose, the result from Shin (2016) was simplified. The major simplification is that all Fourier coefficients and all parameters are presented with Newton’s polynomials for a closed-form solution and engineering application. Therefore, a numerical procedure is not necessary for calculating the coefficients and the parameters any more.

The required Fourier series order is N = 3 in this study, while more than 11 is required by Dean’s stream function theory (DNV, 2010b), and more than N = 64 is required by Fenton’s theory (Rienecker and Fenton, 1981). The series order of 5th-order Stokes theory is 5. Therefore, the required series order is dramatically reduced. The result has less error, is simpler, and gives greater water particle velocities and water particle accelerations than the 5th-order Stokes wave theory. Therefore, the Stokes wave theory is not conservative, and the proposed method is more suitable for offshore structure design.

Acknowledgements

We would like to thank Petrobras for applying the proposed method in the FPSO project.

Notes

No potential conflict of interest relevant to this article was reported.

References

Chakrabarti SK. 1987. Hydrodynamics of Offshore Structures London, UK: Computational Mechanics Publications.
Chaplin JR. 1979;Developments of Stream-Function Wave Theory. Journal of Coastal Engineering 3:179–205. https://doi.org/10.1016/0378-3839(79)90020-6 .
Chappelear JE. 1961;Direct Numerical Calculation of Wave Properties. Journal of Geophysical Research 66(2):501–508. https://doi.org/10.1029/JZ066i002p00501 .
Cokelet ED. 1977;Steep Gravity Waves in Water of Arbitrary Uniform Depth. Philosophical Transaction of the Royal Society of London, Series A, Mathematical and Physical Sciences 286(1335):183–230. https://doi.org/10.1098/rsta.1977.0113 .
De SC. 1955;Contributions to the Theory of Stokes Waves. Proceedings of the Cambridge Philosophical Society 51:713–736.
Dean RG. 1965;Stream Function Representation of Nonlinear Ocean Waves. Journal of Geophysical Research 70(18):4561–4572. https://doi.org/10.1029/JZ070i018p04561 .
Det Norske Veritas. (DNV). 2010a. Column-Stabilized Units. Recommended Practice DNV-RP-C103.
Det Norske Veritas. (DNV). 2010b. Environmental Conditions and Environmental Loads. Recommended Practice DNV-RP-C205.
Fenton JD. 1988;The Numerical Solution of Steady Water Wave Problems. Computers & Geosciences 14(3):357–368. https://doi.org/10.1016/0098-3004(88)90066-0 .
Henry D. 2008;On Gerstner’s Water Wave. Journal of Nonlinear Mathematical Physics 15(2):87–95. https://doi.org/10.2991/jnmp.2008.15.S2.7 .
Gerald CF, Wheatley PO. 2006. Applied Numerical Analysis 7th edth ed. Pearson Education Inc.
Rankine WJM. 1863;On the Exact form of Waves Near the Surface of Deep Water. Philosophical Transactions 153:127–138. https://doi.org/10.1098/rstl.1863.0006 .
Rienecker MM, Fenton JD. 1981;A Fourier Approximation Method for Steady Water Waves. Journal of Fluid Mechanics 104:119–137. https://doi.org/10.1017/S0022112081002851 .
Schwartz LW. 1974;Computer Extension and Analytic Continuation of Stokes’ Expansion for Gravity Waves. Journal of Fluid Mechanics 62(3):553–578. https://doi.org/10.1017/S0022112074000802 .
Skjelbreia L, Hendrickson J. 1960. Fifth Order Gravity Wave Theory. In : Proceedings of 7th Conference on Coastal Engineering. The Hague, Netherlands. p. 184–196. https://doi.org/10.9753/icce.v7.10 .
Shin J. 2016;Analytical Approximation in Deep Water Waves. Journal of Advanced Research in Ocean Engineering 2(1):1–11. https://doi.org/10.5574/JAROE.2016.2.1.001 .
Shin J. 2019;A Regression Analysis Result for Water Waves on Irrotational Flow Over a Horizontal Bed. International Journal of Offshore Polar Engineering. IJOPE 29(4):461–466. https://doi.org/10.17736/ijope.2019.hc17 .
Stokes GG. 1847;On the Theory of Oscillatory Waves. Transactions of the Cambridge Philosophical Society 8:441–473.
Stokes GG. 1880;Supplement to a Paper on the Theory of Oscillatory Wave. Mathematical and Physical Papers 1:314–326.

Article information Continued

Fig. 1.

A conventional coordinate system for a progressive wave (SWL: still water line)

Fig. 2.

A dimensionless coordinate system for a progressive wave

Fig. 3.

Wave profiles and Bernoulli’s constants on the free surface

Fig. 4.

Water particle velocities on the free surface

Fig. 5.

Water particle accelerations on the free surface

Table 1.

Analysis results of the wave with period of 6 s and height of 7.87 m.

Stokes theory This study
Error in KFSBC 1.42% 0
Error in DFSBC 0.37% 7.25×10−2%
Wavelength 64.52 m 64.8 m
Series order 5 3