# On the Study of Second-Order Wave Theory and Its Convergence for a Two-Fluid System

Authors: Chi-Min Liu; Hwung-Hweng Hwung; Ray-Yeng Yang
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Second-order solutions of internal and surface waves in a two-fluid system are theoretically analyzed in this study. Using the perturbation technique, the derivation of second-order solutions for internal waves is revisited, and the results are expressed in one-by-one forms instead of a matrix form. Second-order solutions arising from the interactions of two arbitrary linear waves of different frequencies contain the sum-frequency (superharmonic) and the difference-frequency (subharmonic) components, which are separately examined. Internal Stokes wave being a special case of present solutions is firstly investigated. Next, the convergence of second-order theory and the second-order effects on wave profiles are analyzed. For general cases, the effects of the thickness ratio of two fluids and the ratio of wavenumbers of two first-order waves on second-order wave characteristics, which include transfer functions and particle velocities, are also examined. Moreover, most existing theories for the one-fluid and two-fluid systems can be deduced from present solutions.


---

## Body

## 1. Introduction

A natural wave train is usually composed of thousands of waves with different frequencies. Second-order waves are generated due to the interaction between two arbitrary linear waves. Hence, there are two types of second-order components, the sum-frequency and the difference-frequency components. These two types of second-order terms are the so-called superharmonic and subharmonic components. For second-order wave theories, most of the past studies focused on waves propagating over a single-layer fluid. For sinusoidal waves at first order, Stokes [1] provided second-order solutions describing the self-interaction of first-order waves using wave steepness as a perturbation parameter. An important conclusion is that the wave profile has a sharper crest and a flatter trough rather than a simple sine wave due to second-order effects. Stokes’ theory not only plays an important role in high-order analyses of sinusoidal waves but also virtually forms a foundation for further nonlinear wave studies. It is, however, merely a special case of superharmonic interactions. For random waves with different frequencies and wavelengths, super- and subharmonic waves are naturally generated to form more complicated wave profiles rather than Stokes waves. With the restriction of slightly different frequencies (the narrow-band restriction), Longuet-Higgins and Stewart [2, 3] provided solutions for the subharmonic interactions. Ottesen-Hansen [4] overcame the narrow-band restriction of frequencies and represented the expressions of subharmonic waves in a more adequate form with the introduction of transfer functions, which are defined as the ratio of the second-order amplitude to the first-order amplitudes of waves. Hence, the application of transfer functions enables one to easily measure the second-order effects on wave profiles. Dean and Sharma [5] gave a unified and more compact theory for both superharmonic and subharmonic waves. In practical applications, the theory of second-order surface waves is frequently applied to design the signals of wavemakers [6, 7].Though the results for second-order surface waves are abundant, studies on second-order internal waves are comparatively few. Hunt [8] provided the solutions for internal Stokes waves in a two-fluid system with a rigid lid. Thorpe [9] further examined internal Stokes waves and discussed a number of cases with various thickness ratios and relative wave lengths. He also carried out a comparison between experimental measurements and theoretical results. Tsuji and Nagata [10] applied Stokes theory to study the fifth-order superharmonic interfacial waves between two-layer fluids of infinite thicknesses. This was later checked by Holyer [11] who numerically calculated the series of internal waves to the 31st order using the Fourier expansion procedure. Umeyama [12] theoretically and experimentally studied third-order Stokes waves for slow-mode motion. His derived solution and experimental results are quite coincident. Though internal Stokes wave theory has been widely developed in above studies, it seems that other superharmonic and subharmonic interactions are rarely stressed. Later, Song [13] derived the superharmonic and subharmonic solutions for random internal waves in a two-fluid system bounded by two rigid plates. In oceans, a very slight difference in densities across the pycnocline leads to an infinitesimal displacement at the free surface; the rigid upper boundary is reasonably assumed to be rigid in order to simplify the theoretical analysis without a great loss of accuracy. Liu [14] removed this rigid restriction on the free surface and theoretically derived the second-order super- and subharmonic solutions. The solutions are expressed in matrix form. In practice, his results can be applied to a considerably wide range of density ratios making his solutions superior to Song’s solutions. Preliminary investigations of transfer functions were also addressed in his study.The main purpose of the present paper is to theoretically examine the convergence of second-order wave theory and the second-order effects on wave properties. In Section2, using the perturbation technique, the derivation of second-order solutions for internal waves is revisited and the results are expressed in a one-by-one form instead of a matrix form [14]. Superharmonic and subharmonic interactions are, respectively, examined. In addition, most existing theories for the two-fluid and one-fluid systems are included in the present solutions. Wave characteristics of the slow-mode motion are analyzed. In Section 3, the internal Stokes wave being a special case of superharmonic interactions is investigated. The convergence range for the second-order theory is studied. All second-order wave properties should be subject to this convergence range to ensure the perturbation valid. Wave profiles influenced by second-order effects are subsequently examined. Interactions arising from two arbitrary linear waves are investigated in Section 4. Wave characteristics including the convergence range, transfer functions, and particle velocities are studied. Conclusions are presented in Section 5.


## 2.    Review of the Second-Order Theory for a Two-Fluid System

In this section, the second-order solutions of internal waves in a two-fluid system are briefly revisited. The schematic diagram is shown in Figure1. General assumptions and definitions of symbols are described as follows. The flow in each layer is assumed to be irrotational and inviscid. Φ1 and Φ2 are the velocity potentials of the upper- and lower-layer fluids, respectively. ρr describes the density ratio of the upper fluid ρ1 relative to the lower fluid ρ2. h1 and h2 stand for the undisturbed thicknesses of the upper and the lower layers, respectively. The displacements of the free surface and the interface are represented as ηS and ηI, respectively. The horizontal coordinate is denoted by x, and z represents the vertical coordinate originating at the undisturbed interface and pointing upward. Hence, the governing equations and boundary conditions are expressed as
(1)∂2Φ1∂x2+∂2Φ1∂z2=0,(2)∂2Φ2∂x2+∂2Φ2∂z2=0,(3)∂Φ1∂z=∂ηI∂t+∂Φ1∂x∂ηI∂xatz=ηI,(4)∂Φ2∂z=∂ηI∂t+∂Φ2∂x∂ηI∂xatz=ηI,(5)ρr[gηI+∂Φ1∂t+12(∂Φ1∂x)2+12(∂Φ1∂z)2]=[gηI+∂Φ2∂t+12(∂Φ2∂x)2+12(∂Φ2∂z)2]atz=ηI,(6)∂Φ1∂z=∂ηS∂t+∂Φ1∂x∂ηS∂xatz=h1+ηS,(7)gηS+∂Φ1∂t+12(∂Φ1∂x)2+12(∂Φ1∂z)2=0atz=h1+ηS,(8)∂Φ2∂z=0atz=-h2.Figure 1
Definition sketch.Next, velocity potentials and displacements of the free surface and the interface are expanded in powers of the wave slopeε as shown below
(9)Φ1=Φ1(1)+Φ1(2)+⋯,Φ2=Φ2(1)+Φ2(2)+⋯,ηI=ηI(1)+ηI(2)+⋯,ηS=ηS(1)+ηS(2)+⋯,
where the superscripts indicate the order of magnitude in powers of ε, that is, Φ1(1)=O(ε), Φ1(2)=O(ε2), and so forth. Using the Taylor-series expansion, the boundary conditions on z=ηI and z=h1+ηS ((3) to (7)) are, respectively, expanded about z=0 and z=h1. Accordingly, the amplitudes of internal and surface waves must be small enough to make this expansion valid. Hence, governing equations and boundary conditions for first-order waves are
(10)∂2Φ1(1)∂x2+∂2Φ1(1)∂z2=0,(11)∂2Φ2(1)∂x2+∂2Φ2(1)∂z2=0,(12)∂Φ1(1)∂z-∂ηI(1)∂t=0atz=0,(13)∂Φ2(1)∂z-∂ηI(1)∂t=0atz=0,(14)ρr[gηI(1)+∂Φ1(1)∂t]-[gηI(1)+∂Φ2(1)∂t]=0atz=0,(15)∂Φ1(1)∂z-∂ηS(1)∂t=0atz=h1,(16)gηS(1)+∂Φ1(1)∂t=0atz=h1,(17)∂Φ2(1)∂z=0atz=-h2.
The first-order solutions are assumed to be
(18)ηI(1)=Re[∑j=m,najeiθj],(19)ηS(1)=Re[∑j=m,ndjeiθj],(20)Φ1(1)=Im{∑j=m,n[Ajcosh(kjz)+Bjsinh(kjz)]·eiθj},(21)Φ2(1)=Im{∑j=m,nCjcosh[kj(z+h2)]·eiθj},
where θj=kjx-σjt+ψj (ψj is the phase) and Re and Im indicate the real and imaginary parts of the complex, respectively. The assumed velocity potentials in (20) and (21) have automatically satisfied boundary conditions (10), (11), and (17). Inserting (18) to (21) into (12) to (16) results in the following relations between the assumed coefficients:
(22)Aj=ajσj[(1-1ρr)g+σj2ρrkjtanh(kjh2)],(23)Bj=σjajkj,(24)Cj=σjajkjsinh(kjh2),(25)dj=σj2σj2cosh(kjh1)-gkjsinh(kjh1)·aj,(26)σj4[coth(kjh1)coth(kjh2)+ρr]-σj2gkj[coth(kjh1)+coth(kjh2)]+(1-ρr)g2kj2=0.
Equation (26) is the dispersion relation in which two solutions of σj2 with a given wavenumber kj exist. It usually refers to the fast-mode and slow-mode solutions. For a case with a very small density difference (ρ2-ρ1≪ρ2), the two solutions of σj2 are shown below:
(27)σj2=gkjtanh[kj(h1+h2)],forthefast-modemotion,σj2=gkj·1-ρrcoth(kjh1)+coth(kjh2),fortheslow-modemotion.
The slow-mode motion will be excited when the density difference tends to zero. In this mode, the free surface displacement will be very small in comparison with the internal wave displacement. Therefore, the free surface boundary was usually replaced by a rigid lid in many studies for simplification. Subsequently, the governing equations and boundary conditions at the second order are shown as
(28)∂2Φ1(2)∂x2+∂2Φ1(2)∂z2=0,(29)∂2Φ2(2)∂x2+∂2Φ2(2)∂z2=0,(30)∂Φ1(2)∂z-∂ηI(2)∂t=-∂2Φ1(1)∂z2ηI(1)+∂Φ1(1)∂x∂ηI(1)∂xatz=0,(31)∂Φ2(2)∂z-∂ηI(2)∂t=-∂2Φ2(1)∂z2ηI(1)+∂Φ2(1)∂x∂ηI(1)∂xatz=0,(32)ρr[gηI(2)+∂Φ1(2)∂t]-[gηI(2)+∂Φ2(2)∂t]=-ρr[∂2Φ1(1)∂t∂zηI(1)+12(∂Φ1(1)∂x)2+12(∂Φ1(1)∂z)2]+[∂2Φ2(1)∂t∂zηI(1)+12(∂Φ2(1)∂x)2+12(∂Φ2(1)∂z)2]atz=0,(33)∂Φ1(2)∂z-∂ηS(2)∂t=-∂2Φ1(1)∂z2ηS(1)+∂Φ1(1)∂x∂ηS(1)∂xatz=h1,(34)gηS(2)+∂Φ1(2)∂t=-∂2Φ1(1)∂t∂zηS(1)-12(∂Φ1(1)∂x)2-12(∂Φ1(1)∂z)2atz=h1,(35)∂Φ2(2)∂z=0atz=-h2.
The second-order components generated by two first-order m- and n-waves are assumed as
(36)ηI(2)=Re[∑m,namn±eiθmn±],ηS(2)=Re[∑m,ndmn±eiθmn±],Φ1(2)=Im{∑m,n[Amn±cosh(kmn±z)+Bmn±sinh(kmn±z)]·eiθmn±∑m,n[Amn±cosh(kmn±z)+Bmn±sinh(kmn±z)]},Φ2(2)=Im{∑m,n[Cmn±cosh[kmn±(z+h2)]·eiθmn±]},
where θmn±=kmn±x-σmn±t+ψmn±, kmn±=km±kn, σmn±=σm±σn, ψmn±=ψm±ψn, and the superscripts plus and minus, respectively, denote the superharmonic and subharmonic components, respectively. By inserting (36) into (30) to (34), amn±, dmn±, Amn±, Bmn±, and Cmn± are solved to be
(37)Amn±=1Δ±{[(σmn±)2-gkmn±Ct1±]·[(σmn±)2(kmn±)-1Ct2±(χ2-χ1)+(1-ρr)gχ1-σmn±χ3±(σmn±)2(kmn±)-1]+(Sh1±)-1[(1-ρr)gkmn±-(σmn±)2Ct2±]·[gχ4+σmn±χ5±]{[(σmn±)2-gkmn±Ct1±]},(38)Bmn±=1Δ±{[(σmn±)2Ct1±-gkmn±]·[(σmn±)2(kmn±)-1Ct2±(χ1±-χ2±)-(1-ρr)gχ1±+σmn±χ3±(σmn±)2(kmn±)-1]-ρr(σmn±)2(Sh1±)-1[gχ4±+σmn±χ5±]},(39)Cmn±=(Sh2±)-1Δ±·{ρr(σmn±)2(χ2±-χ1±)·[(σmn±)2(kmn±)-1-gCt1±]+[(σmn±)2Ct1±-g(kmn±)]·[σmn±χ3±-(1-ρr)gχ2±]-ρr(σmn±)2(Sh1±)-1[gχ4±+σmn±χ5±]{ρr(σmn±)2(χ2±-χ1±)},(40)dmn±=1Δ±{{((σmn±)2Ct2±-(1-ρr)gσmn±kmn±)}(1-ρr)gkmn±(σmn±Ct1±χ4±+kmn±χ5±)+(Ch1±Ct1±-Sh1±)·[((σmn±)2Ct2±-(1-ρr)gσmn±kmn±)χ1±-(σmn±)3Ct2±χ2±+(σmn±)2kmn±χ3±]-(σmn±)2[σmn±(ρr+Ct1±Ct2±)χ4±+kmn±(ρrCt1±+Ct2±)χ5±]{((σmn±)2Ct2±-(1-ρr)gσmn±kmn±)}},(41)amn±=1Δ±{ρrσmn±[gkmn±Ct1±-(σmn±)2]χ1±+(σmn±Ct2±χ2±-kmn±χ3±)·[gkmn±-(σmn±)2Ct1±]-ρrσmn±kmn±(Sh1±)-1(gχ4±+σmn±χ5±)},
in which
(42)χ1±=-aman2kmn±·{[(1-ρr-1)gkmσm+σmρrcoth(kmh2)]+[(1-ρr-1)gknσn+σnρrcoth(knh2)]},χ2±=-aman2kmn±[σmcoth(kmh2)+σncoth(knh2)],χ3±=aman2ρr·{·[(1-ρr-1)gknσn+σnρrcoth(knh2)]}σm2+σn2±σmσn+[(1-ρr-1)gkmσm+σmρrcoth(kmh2)]·[(1-ρr-1)gknσn+σnρrcoth(knh2)]}-aman2{σm2+σn2-σmσnsinh(kmh2)sinh(knh2)·cosh(kmn∓h2)·[(1-ρr-1)gknσn+σnρrcoth(knh2)]}},χ4±=-aman2kmn±·{σn2σn2cosh(knh1)-gknsinh(knh1)·[((1-ρr-1)gkmσm+σmρrcoth(kmh2))·cosh(kmh1)+σmsinh(kmh1)((1-ρr-1)gkmσm+σmρrcoth(kmh2))]+σm2σm2cosh(kmh1)-gkmsinh(kmh1)·[((1-ρr-1)gknσn+σnρrcoth(knh2))·cosh(knh1)+σnsinh(knh1)((1-ρr-1)gknσn+σnρrcoth(knh2))]{σn2σn2cosh(knh1)-gknsinh(knh1)},χ5±=aman2·{σn2σn2cosh(knh1)-gknsinh(knh1)·[((1-ρr-1)gkm+σm2ρrcoth(kmh2))·sinh(kmh1)+σm2cosh(kmh1)((1-ρr-1)gkm+σm2ρrcoth(kmh2))]+σm2σm2cosh(kmh1)-gkmsinh(kmh1)·[((1-ρr-1)gkn+σn2ρrcoth(knh2))·sinh(knh1)+σn2cosh(knh1)((1-ρr-1)gkn+σn2ρrcoth(knh2))]-[((1-ρr-1)gkmσm+σmρrcoth(kmh2))·((1-ρr-1)gknσn+σnρrcoth(knh2))∓σmσn((1-ρr-1)gkmσm+σmρrcoth(kmh2))]Ch1∓-[((1-ρr-1)gknσn+σnρrcoth(knh2))σm∓σn·((1-ρr-1)gkmσm+σmρrcoth(kmh2))]·Sh1∓((1-ρr-1)gknσn+σnρrcoth(knh2))},
where
(43)Δ±=(σmn±)4(ρr+Ct1±Ct2±)Δ±=-(σmn±)2gkmn±(Ct1±+Ct2±)+(1-ρr)g2(kmn±)2,Shi±≡sinh(kmn±hi),Chi±≡cosh(kmn±hi),Cti±≡coth(kmn±hi).
Present two-dimensional solutions can be readily converted to three-dimensional forms by using the vectors of wavenumbers instead of using the scalar forms that are presented. It is also noted that the above solutions can be reduced to most existing results including surface waves in a one-fluid system (ρr=0) and internal waves in a two-fluid system [8, 9, 12, 13]. In the following analyses, the slow-mode motion is considered, and the density ratio ρr is assumed to be 0.99 for the purpose of approaching the real circumstance in ocean.


## 3. Internal Stokes Wave and Its Convergence

In this section internal Stokes wave and its convergence are examined. First, the first-order solution is assumed as(44)ηI(1)=aS(1)cos(kx-σt).
Based on (41), the second-order solution is solved to be
(45)ηI(2)=aS(2)cos2(kx-σt)=aS(1)(kaS(1))·HΔS·{{[EC^t1+1-EH-HC^t1(HC^t1-1)+(E2-1)2]}ρrH(Ct12-1)0.5·[EC^t1+1-EH-HC^t1(HC^t1-1)+(E2-1)2]+ρrE(2H-Ct1)+(2HCt1-14)·(4Ct2C^t2+3ρr+ρrE2-2+C^t22-1){[EC^t1+1-EH-HC^t1(HC^t1-1)+(E2-1)2]}}·cos2(kx-σt)≡aS(1)(kaS(1))·ℱ(kh1,kh2)·cos2(kx-σt),
where
(46)ℱ(kh1,kh2)≡HΔS·{{[EC^t1+1-EH-HC^t1(HC^t1-1)+(E2-1)2]}ρrH(Ct12-1)0.5·[EC^t1+1-EH-HC^t1(HC^t1-1)+(E2-1)2]+ρrE(2H-Ct1)+(2HCt1-14)·(4Ct2C^t2+3ρr+ρrE2-2+C^t22-1){[EC^t1+1-EH-HC^t1(HC^t1-1)+(E2-1)2]}},Cti=coth(2khi),C^ti=coth(khi),H=σ2gk,E=ρr-1C^t2+H-1(1-ρr-1),ΔS=4H2(ρr+Ct1Ct2)-2H(Ct1+Ct2)+(1-ρr).
The above solution agrees with the theoretical solutions recently given by Umeyama [12] and that obtained earlier by Thorpe [9] (some misprints in (A1.1) of Thorpe’s study should be corrected). It is also remarked that the dispersion relation remains unchanged until the third order. Above second-order solution is found to depend only on kh1, kh2, and the wave slope kaS(1). Before further investigations are carried out, the convergence for validating the perturbation has to be determined. Though the convergence for surface Stokes waves had been examined by Levi-Civita [15], such an examination for convergence has not yet been well studied for internal Stokes wave. To meet the convergence, the amplitude of second-order wave is required to be less than that of first-order wave, that is,
(47)|kaS(1)·ℱ(kh1,kh2)|<1.
In general, this convergence requires that the wave slope is sufficiently small. As for the detailed observation, the convergence zone determined by (47) for the slow-mode motion is shown in Figure 2. It is seen that larger values of wave slope result in a more rigorous condition between kh1 and kh2. All of the following analyses in this section should be subject to this convergence condition. Ratios of the second-order amplitude relative to the first-order amplitude are displayed in Figures 3 and 4. Figure 3 shows the results under different wave slopes. Solid and dash curves represent the results of internal-wave amplitude ratios aS(2)/aS(1) and surface-wave amplitude ratios dS(2)/dS(1), respectively. Waves with larger wave slopes are found to have stronger second-order effects on amplitudes. The effects of the thickness ratio of two fluid layers are shown in Figure 4. It is found that as values of h1/h2 become small enough, the amplitude ratios become negative. Positive values of amplitude ratios imply that the wave crests (with respect to the below fluid) will sharpen and the wave troughs will flatten. Such a phenomenon is analogous to that appearing in a single-fluid system. On the contrary, negative values suggest the occurrence of a sharper trough and a flatter crest (also with respect to the below fluid). Other properties of Stokes waves will be discussed in the next section.Figure 2
Convergence zone for internal Stokes waves with various wave slopes.Figure 3
Ratios of the second-order amplitude to the first-order amplitude for Stokes waves with various wave slopes. Solid lines for internal waves and dash lines for surface waves.Figure 4
Ratios of the second-order amplitude to the first-order amplitude for Stokes waves with various thickness ratios. Solid lines for internal waves and dash lines for surface waves.


## 4. Superharmonic and Subharmonic Waves and the Convergence

Superharmonic and subharmonic interactions of two arbitrary first-order waves (m-wave and n-wave) are examined in this section. Because six length parameters, h1, h2, km-1, kn-1, am, and an, appear in the second-order solutions, the condition am=an, which will cause the strongest interaction, is assumed to investigate the effects of two important parameters, kn/km and h1/h2, on second-order wave properties. The convergence is firstly examined. It may require the ratio of the second-order amplitude relative to the summation of two first-order amplitudes to be less than unity:
(48)|amn±am+an|<1.
For Stokes waves, the convergences defined by (47) and (48) are identical. As for arbitrary wave interactions, Figure 5 shows the convergences under various wavenumber ratios kn/km for superharmonic (solid lines) and subharmonic (dash line) interactions. It is found that the smaller wavenumber ratio kn/km results in a more rigorous convergence condition between kmh1 and kmh2. It implies that when one wave (m-wave) interacts with another shorter wave (n-wave), it yields a wider convergence zone. In addition, the convergence for the subharmonic case is more severe than that for the superharmonic case. As mentioned in Section 3, the larger wave slope will narrow the convergence range.Figure 5
Convergence zone for the super- and the subharmonic interactions with various wavenumber ratios. Solid line for superharmonic interactions and dash lines for subharmonic interactions.After examining the convergence condition, transfer functions defined to measure the second-order effects on wave profiles are subsequently investigated. The superharmonic and subharmonic transfer functions for surface (with the subscriptS) and internal (with the subscript I) waves are, respectively, defined as
(49)Gmn,S±≡dmn±dmdn,Gmn,I±≡amn±aman.
Note that the transfer functions are independent of the amplitudes of both first-order waves. Figures 6, 7, and 8 show the behaviors of the superharmonic transfer functions for internal and surface waves over the depth range 0<kmh2<5. All of these figures are plotted under the same condition h1=h2. The transfer functions for internal waves with various wavenumber ratios are plotted in Figure 6. A larger value of the wavenumber ratio results in stronger second-order effects on the wave profile. In addition, the tendency of Gmn,I+ is quite similar to that for the single-fluid case (see Sand and Mansard [6, Figure 2]). Results for surface waves are plotted in Figure 7. The magnitudes of the transfer functions for surface waves are much greater than those for internal waves. It is also found that all transfer functions with different wavenumber ratios are always negative over the whole range of kmh2. Figure 8 displays the amplitude ratio of surface waves to internal waves at the first order (dash line) and the second order (solid line). In virtue of only a very small undulation at the free surface for the slow-mode motion, the values of these ratios are multiplied by one thousand for a clearer understanding. It reveals that internal and surface waves are 180∘ out of phase for the slow-mode motion. Moreover, a larger amplitude ratio of wavenumbers of the first-order waves causes a faster decay of second-order effects while kmh2 grows.Figure 6
Superharmonic transfer functions for internal waves with various wavenumber ratios.Figure 7
Superharmonic transfer functions for surface waves with various wavenumber ratios.Figure 8
Ratios of the amplitude of surface wave to that of internal wave for superharmonic interactions (solid lines for second-order waves with various wavenumbers and dash lines for first-order waves).Results for the subharmonic cases are plotted in Figures9, 10, 11, 12, and 13. For the condition h1=h2, Figure 9 shows that values of transfer functions for internal waves remain negative over the whole range of kmh2. The behaviors are analogous to those of the single-fluid case (see Sand [16, Figure 2]). Figure 10 indicates that the transfer functions for surface waves are, respectively, positive and negative in the long-wave and the short-wave regimes. Similar to results for the superharmonic case, the thickness ratios of two fluid layers strongly dominate the signs and magnitudes of the transfer functions, as shown in Figures 11 and 12. Motions of the second-order internal and surface waves could be either in phase or 180∘out of phase depending on the values of kmh2, as shown in Figure 13.Figure 9
Subharmonic transfer functions for internal waves with various wavenumber ratios.Figure 10
Subharmonic transfer functions for surface waves with various wavenumber ratios.Figure 11
Subharmonic transfer functions for internal waves with various thickness ratios.Figure 12
Subharmonic transfer functions for surface waves with various thickness ratios.Figure 13
Ratios of the amplitude of surface wave to that of internal wave for subharmonic interactions (solid lines for second-order waves with various wavenumbers and dash lines for first-order waves).Another important wave property, the horizontal and vertical velocities, is also investigated. They can be directly obtained by differentiating the velocity potentials with respect tox and z:
(50)ui=ui(1)+ui(2)=∂∂x(Φi(1)+Φi(2)),vi=vi(1)+vi(2)=∂∂z(Φi(1)+Φi(2)).
Since the phase velocities of all first-order waves and the corresponding second-order wave are different, the phases of m-wave, n-wave, and the associated second-order wave in (50) are set to be zero (θm=θn=θmn±=0) for examination. For the superharmonic case, ratios of the second-order velocity to the first-order velocity with different wavenumber ratios for the shallower-water (kmh2=0.5) and the deeper-water (kmh2=1.0) regimes are, respectively, depicted in Figures 14 and 15. Solid lines and dash lines, respectively, stand for the horizontal and vertical velocity ratios along the normalized z-axis. In the shallower-water regime, it is also found that the second-order effects are stronger in the upper layer. As for the deeper-water regime, larger wavenumber ratios lead to greater second-order effects near the interface. All second-order effects near the free surface and the bottom are comparatively slight.Figure 14
Velocity ratios of the second-order waves to the first-order waves along thez-axis for the superharmonic interactions with kmh2=0.5. Solid lines for horizontal velocity ratios and dash lines for vertical velocity ratios.Figure 15
Velocity ratios of the second-order waves to the first-order waves along thez-axis for the superharmonic interactions with kmh2=1.0. Solid lines for horizontal velocity ratios and dash lines for vertical velocity ratios.Results for the subharmonic case are plotted in Figures16 and 17 for the conditions kmh2=0.5 and kmh2=1.0, respectively. The results show that second-order effects are stronger in the lower fluid. It is also observed that stronger second-order effects occur near the interface in the case of larger wavenumber ratios. Moreover, the effects of the thickness ratio h1/h2 are not as important as those of the wavenumber ratio kn/km. The reason is that the particle velocities far from the interface are usually smaller in comparison with those near the interface. Hence, the corresponding diagrams for various values of h1/h2 are ignored herein.Figure 16
Velocity ratios of the second-order waves to the first-order waves along thez-axis for the subharmonic interactions with kmh2=0.5. Solid lines for horizontal velocity ratios and dash lines for vertical velocity ratios.Figure 17
Velocity ratios of the second-order waves to the first-order waves along thez-axis for the subharmonic interactions with kmh2=1.0. Solid lines for horizontal velocity ratios and dash lines for vertical velocity ratios.Finally, we note that particle horizontal velocities of two fluids are identical at the interface due to the boundary conditions, (3) and (4). As viscous effects are not taken into account, particle horizontal velocities of two fluids do not satisfy the continuity requirement at the interface.


## 5. Conclusions

Second-order theory for internal waves in a two-fluid system and the corresponding convergence are presented in this paper. Internal Stokes waves, the second-order interactions of two arbitrary linear waves, the associated convergence condition, and the particle velocities are examined. The crucial conclusions are summarized in the following.(1)
For internal Stokes waves, our second-order solutions are in agreement with those of past studies, and the convergence of the second-order theory is found to depend only onkh1, kh2, and the wave slope kaS(1). A larger value of wave slope will require a more rigorous convergence condition between kh1 and kh2. Larger wave slopes also result in stronger second-order effects either on internal or surface waves. The depth ratio h1/h2 is an important parameter for measuring wave profiles influenced by second-order effects. As h1/h2 becomes small enough, the internal wave profile will have a sharp trough and a flattened crest with respect to the below fluid. On the contrary, internal waves with sharp crests and flattened troughs will occur as h1/h2 becomes large enough.(2)
For second-order interactions arising from two arbitrary waves, a smaller wavenumber ratiokn/km results in a more rigorous convergence condition between kmh1 and kmh2. The convergence for the subharmonic interaction is more severe than that for the superharmonic case. In addition, larger values of wave slopes will narrow the convergence range.(3)
Transfer functions measure the magnitude of second-order effects on wave profiles. For internal waves, transfer functions of all wavenumber ratios respectively remain positive and negative for the super- and subharmonic cases. The thickness ratio of the two fluid layers is found to influence not only the magnitude but also the sign of transfer functions. In the long-wave regime, transfer functions usually have opposite signs for the thicker-upper-layer case and the thicker-lower-layer case. Moreover, the first-order surface and internal waves are180∘ out of phase for the slow-mode motion. The second-order surface and internal waves are also 180∘ out of phase for the superharmonic case and are in phase over most of the range of kmh2 for the subharmonic case.(4)
As for particle velocities, larger wave number ratios cause stronger second-order effects near the interface. All second-order effects near the free surface and the bottom are much slighter than those near the interface.Finally, it is noted again that most existing theories for two-fluid and one-fluid systems are included in the present solutions.



---

*Source: 253401-2013-05-02.xml*