1. Introduction
Internal solitary waves (ISWs), which are a kind of nonlinear internal wave, are ubiquitous in marginal sea areas [
1,
2,
3]. Based on the balance between nonlinearity and dispersion, ISWs can maintain a stable shape for long distances in the horizontal direction [
4]. Due to the opposite flow direction induced by the upper and lower layers of the wave, a strong shear flow can be vertically induced, which can have a significant impact on the safety of marine structures and underwater platforms [
5,
6,
7,
8,
9,
10]. In addition, ISWs with large amplitudes, which are often accompanied by relatively large vertical velocities, can have an important impact on seawater mixing [
11], nutrient transport [
12] and sediment resuspension [
13] in the ocean. Therefore, studying the vertical structure characteristics of ISWs is greatly significant.
Vertical structures of ISWs usually refer to the distribution of amplitudes at different depths [
14]. A large amount of information can be obtained from the vertical structure of ISWs, including the depth of the maximum amplitude and the vertical scale of ISWs. Additionally, vertical structures can reflect the evolution of ISWs to a certain extent [
15]. Many theoretical studies on vertical structures of ISWs have been carried out thus far [
16,
17,
18,
19,
20,
21,
22]. Fliegel and Hunkins [
18] first calculated the dispersion relation of internal waves by using the Thomson–Haskell method and obtained vertical structures of internal waves in different modes. Based on the Thomson–Haskell method, Cai and Gan [
19] calculated the internal wave dispersion relation without using the Boussinesq approximation and theoretically analyzed the vertical structure characteristics of different modes at two stations in the South China Sea. Zhang and Gao [
20] obtained vertical structures and the dispersion relation of internal waves by using the transformation method of Russian scholars. Based on a numerical method, Deng et al. [
21] analyzed the vertical structures of different modes of internal waves in background currents. Kurkina et al. [
22] solved the theoretical vertical structures of first-mode and second-mode ISWs with nonlinear corrections in two-layer and three-layer fluids.
The propagation and evolution characteristics of ISWs have been investigated in several studies by using in situ observational data [
23,
24,
25,
26,
27]. Vlasenko et al. [
28] obtained measured data of ISWs in the Strait of Messina and found that the observed vertical structures of ISWs are consistent with numerical simulation results. Geng et al. [
29] obtained the vertical structures of ISWs near the Dongsha Atoll of the South China Sea using seismic technology and found that the observational vertical structures of ISWs are generally consistent with the theoretical results except when ISWs interact with the topography. Similarly, Gong et al. [
14] compared the vertical structures of ISWs obtained by the seismic method with two theories. The results show that ISWs with low nonlinearity can be described by a linear vertical structure function, while ISWs with high nonlinearity are in agreement with the first-order nonlinear vertical structure function. Because obtaining refined vertical structure data from field observations is difficult, the verification of the theoretical vertical structure of ISWs based on observational data is relatively limited. In terms of laboratory experiments, many experimental methods have been used to study the characteristics of ISWs [
30,
31,
32,
33,
34,
35]. At present, experimental research objects of ISWs are not limited to the first mode. Many experimental studies on the generation mechanism and evolution characteristics of second-mode ISWs have been carried out [
36,
37,
38]. However, most experimental studies focus on the propagation and evolution characteristics of ISWs at the interface (position of maximum amplitude), while there are few studies on the vertical structure of ISWs in stratified fluids.
Korteweg de Vries’ (KdV) theory is the most commonly used theoretical model to describe the propagation and evolution characteristics of ISWs in laboratory experiments [
39,
40,
41]. Since weak nonlinearity is an important condition for validating the KdV model, KdV theory is not applicable to describe ISWs with large amplitudes [
42,
43,
44]. Michallet and Barthelemy [
45] compared the theoretical and experimental results of the waveform, phase velocity and frequency–amplitude relationship of ISWs and found that KdV theory can describe small-amplitude solitary waves in different stratified fluids well. Grue et al. [
46] analyzed the propagation characteristics of ISWs in two-layer fluids and found that KdV theory can provide predictions that are consistent with the experimental results of small-amplitude ISWs. Du et al. [
47] compared the experimental results of the elevation- and depression-type ISWs with the results from the following three different theoretical models: KdV, extended KdV (eKdV) and Miyata–Choi–Camassa (MCC). The results show that the experimental results conform to those of KdV theory when the ratio of the amplitude of ISWs to the total water depth is less than 0.04. However, KdV theory is presently used to describe the characteristics of ISWs at the interface in laboratory experiments, and the validity of KdV theory describing the properties of ISWs at different depths in stratified fluids is not clear.
In view of the above discussion, an experimental method is used in this paper to study the vertical structure of the amplitude of ISWs in a large, gravitationally stratified fluid flume and to analyze the domain of validity of KdV theory by comparing the characteristics of ISWs at different depths with theoretical models. The influence of environmental factors on the vertical structure characteristics of ISWs is further analyzed by setting different stratifications and amplitudes. This paper is organized in the following manner: The KdV model is introduced in
Section 2, including the two-layer model and the continuously stratified model. The experimental techniques and methods are presented in
Section 3. The results and discussion are provided in
Section 4, and the conclusion is presented in
Section 5.
2. Theoretical Models
Considering a stratified fluid of a finite depth with rigid boundaries where buoyancy frequency is a continuous function of the depth, the eigenfunction
is determined by the solution of the boundary value problem:
where
is the buoyancy frequency, which is determined by the formula
;
is the linear phase speed of ISWs; and
is the total water depth. Equation (1) can be solved by using the Thomson–Haskell method [
18], and
can be obtained. Then,
can be calculated under the normalization condition of
.
For ISWs with a maximum amplitude
in a stratified fluid, their amplitude at different depths
can be calculated by the following expression (with the nonlinear terms removed):
Under the Boussinesq approximation, the vertical isopycnal displacement
of ISWs propagating in the positive direction of axis
can be described by the classical KdV equations:
where
is the time, and
is the horizontal coordinate along the propagation direction of ISWs. The parameters
and
are called the nonlinear coefficient and the dispersion coefficient, respectively.
The well-known solitary wave solution to Equation (3) is
where
is the amplitude of ISWs, and the nonlinear velocity
and the characteristic width
of this soliton can be calculated by:
The characteristic frequency of the soliton can be obtained by . Since and are both relative to , the characteristic frequency will change with . Moreover, the obtained frequency–amplitude relationship varies for the two-layer fluid model and the continuously stratified fluid model due to the different calculation methods of , and in these two models.
In the continuously stratified model,
can also be obtained when the eigenfunction
is calculated by the Thomson–Haskell method. The terms
,
can be determined by the following formula:
In the two-layer model,
,
and
can be determined by the following expressions:
where
and
are the thicknesses of the upper and lower layers, respectively, and
and
are the densities of the upper and lower layers, respectively.
3. Experimental Device and Method
The experiments were performed in a gravitationally stratified fluid flume that is 1200 cm long, 120 cm wide and 100 cm high. At the right end of the flume, the wave maker aimed to achieve gravity collapse in the stratified fluid, which is presented in
Figure 1. The step length and depth are
and
, respectively, in the experimental generation principle of ISWs. This device can produce high-quality elevation- or depression-type ISWs at arbitrary pycnocline depths. At the left end of the flume, a wedge-shaped wave absorber was installed. According to the amplitude of ISWs and interface depth, the wedge angle can be adjusted to effectively eliminate the reflected waves. A stratified fluid was created by using the classic “double-tank” method. At first, the flume was filled with fresh water of density
to the depth
. Then a brine solution of density
was slowly injected beneath the fresh water through two small openings along the bottom of the flume until the brine depth reached
, and the total depth was
. Due to the mixing of fresh water in the upper layer and brine solution in the lower layer, a continuous density distribution structure with a certain thickness was formed at the interface, namely the pycnocline, which can simulate the actual vertical density distribution of the ocean with the characteristics of a typical pycnocline.
The vertical distribution of density in stratified fluids and vertical structures of ISWs was measured by an internal wave dynamic measurement system. This system consists of a conductivity probe array, signal control circuit and data processing system. A vertically displaceable probe was used to measure the vertical density distribution of the fluid. A conductivity probe array consisting of seven probes vertically arranged at equal intervals was used to measure the vertical structure of ISWs. The horizontal distance between the probes and the wave maker is 6 m, so the fully developed ISWs can be detected (
Figure 1). Adjusting the position of the conductivity probe array to cover the whole pycnocline according to the stratification conditions and amplitude of ISWs, the fluctuation of ISWs at different depths can be completely captured. The origin of the axis
is set at the surface with the positive axis
running downward.
Figure 2a shows the vertical distribution of isopycnic lines of ISWs in a certain stratification. The back portion of the wave was contaminated by the small-amplitude disturbances in the tail. Therefore, as shown in
Figure 2b, only the forward portion of the wave was used to calculate the amplitude and characteristic length of ISWs [
40,
41]. The characteristic length and frequency of ISWs can be experimentally calculated by the following expressions (Equations (10) and (11), respectively):
4. Results and Discussion
To fully study the vertical structure characteristics of ISWs, four stratification conditions were set in the experiment, and two types of ISWs (elevation- and depression-type) were generated during operation. The specific experimental conditions are summarized in
Table 1. S1 is a representative of the stratification in the deep sea. S2–S4 are used to simulate the different stratifications in the continental shelves. Both elevation-type (
) and depression-type (
) ISWs exist in real oceans. Compared with the depression-type ISWs, the amplitudes of the elevation-type ISWs are generally smaller. During the experiment, the total water depth
and the densities of upper and lower layers (
and
, respectively) remained unchanged. Additionally, the dimensionless maximum amplitude
ranged from 0.008 to 0.120.
4.1. Vertical Structure of the Amplitude
According to the theoretical model, the distribution of the buoyancy frequency
can be obtained from the experimental density profile
, and the vertical modal function of the internal wave
can be subsequently calculated.
Figure 3 illustrates the vertical distribution of ISWs in S3, with
Figure 3a showing the density profile and buoyancy frequency distribution and
Figure 3b presenting the vertical modal function of ISWs obtained by the Thomson–Haskell method.
Figure 3 shows that under the condition of this density stratification, the variation in the buoyancy frequency and internal wave vertical modal function with depth are consistent, which show a trend of first increasing and then decreasing with an increasing depth. The maximum value of the buoyancy frequency profile is at a depth of
, and the maximum value of the corresponding vertical modal function is at a depth of approximately
(red line in
Figure 3b). There is a small deviation between the depth of the maximum value of the buoyancy frequency profile and the depth of the maximum value of the vertical modal function that is obtained by the Thomson–Haskell method, which can also be found from the theoretical results of Kurkina et al. [
22].
Figure 4 demonstrates a comparison between the theoretical and experimental results of the vertical structure of the amplitude of ISWs in four stratification conditions. Under different stratification conditions, whether elevation-type or depression-type, the theoretical vertical structure of the amplitude is basically consistent with most of the experimental results. In addition, the amplitude initially increases and then decreases with an increasing depth both theoretically and experimentally. For S1, the depths of the maximum amplitudes of ISWs under different cases are basically the same, while the depth of the maximum amplitude of ISW with large amplitude (
) moves up. For S2, most of the vertical distributions of the amplitude conform to the theory, while the vertical amplitude structure of ISWs with relatively large amplitudes (
and
) moves up. For S3, the depths of the maximum amplitudes under different cases and the vertical structure of the amplitude are consistent with the theory. For S4, similar to S3, the theory can describe the vertical structure of ISWs under different cases well.
Figure 5 shows the variation in root mean square error
with increasing
under different stratification conditions. Root mean square errors in different stratifications are relatively small when
. However, once
, the root mean square errors start to increase significantly and exceed the allowable range.
For S1 and S2, the vertical structures of ISWs basically conform to the theory in most cases, while for ISWs with large amplitudes (), the experimental vertical structures deviate from the theoretical model. These results occur because the effect of the nonlinear terms will be more significant with an increasing amplitude. In addition, with an increasing , the difference between the theoretical and experimental vertical structures of ISWs with large amplitudes is more obvious. However, for S3 and S4, the vertical structures of ISWs are consistent with the theoretical value due to the relatively small amplitudes ().
4.2. Vertical Distribution of the Characteristic Frequency
As shown in
Figure 6, the variation in the theoretical characteristic frequency of ISWs with amplitudes in continuously stratified fluid and two-layer fluid can be obtained according to the expression of the characteristic frequency
. The results show that the characteristic frequency of ISWs increases with an increasing amplitude whether it is for the continuously stratified fluid or the two-layer fluid. The theoretical value of the characteristic frequency in the two-layer fluid is larger than that in the continuously stratified fluid at the same amplitude. Generally, this is mainly due to smaller absolute values of
and larger values of
in the continuously stratified model for different stratifications. In addition, the variation in characteristic frequency is closely related to the stratification conditions. For the depression-type ISWs, both the characteristic frequency under the same amplitude and the difference in the theoretical values between the two-layer model and the continuously stratified model will correspondingly decrease with an increasing
. The variation in the characteristic frequency of the elevation-type ISWs with the amplitude in S4 is basically consistent with that in S3. These trends are also closely linked to values of
and
in different stratification conditions.
The theoretical and experimental results of the frequency–amplitude relationship of ISWs in S1–S4 are presented in
Figure 7. For S1, the experimental characteristic frequencies at small amplitudes (
) are slightly larger than their theoretical values. The characteristic frequencies of ISWs are close to the theoretical values of the two-layer model under the condition of
. When
is satisfied, the experimental characteristic frequencies approach the theoretical values of the continuously stratified model. The characteristic frequencies of ISWs are smaller than the theoretical frequencies for relatively large amplitudes (
). For S2, which is similar to S1, the experimental characteristic frequencies at small amplitudes (
) are larger than their theoretical values. Regarding the amplitude that satisfies
, the characteristic frequencies of ISWs are basically between the theoretical values of the two models but closer to the continuously stratified model. The characteristic frequencies of ISWs are still smaller than the theoretical frequencies for relatively large amplitudes (
). For S3, the theoretical values of the two models are coinciding, but the characteristic frequencies of most ISWs obtained by experiments are larger. For S4, the theoretical values of the two models are also close, and the characteristic frequencies obtained experimentally are dispersed with the amplitude. Different from other stratification conditions, the experimental values of the characteristic frequencies can be larger or smaller than the theoretical ones under the same amplitude.
Comparing different stratification conditions, the characteristic frequency of ISWs increases with an increasing amplitude for the four stratifications, which is consistent with the trend of the theoretical model. Among these stratifications, S1 and S2 can describe the variation in the characteristic frequency with the amplitude to some extent by the two-layer model or the continuously stratified model. For S3 and S4, these two models have difficulty describing the frequency–amplitude relationship. For the depression-type ISWs (S1–S3), the characteristic frequency distribution is more concentrated with the amplitude, while the characteristic frequency distribution of elevation-type ISWs (S4) is relatively dispersed with the amplitude. One possible reason may be that the characteristic frequencies do not change much with the amplitude due to the interactive relationship of the waveforms at different depths. Therefore, the experimental characteristic frequency does not increase with an increasing amplitude as fast as that in models. In addition, we also consider that the waveforms of elevation-type ISWs can be easily affected by some disturbances due to the small amplitudes, which results in the dispersive amplitude–frequency relationship in S4. Another explanation for the discrepancy is that the density profile of the fluid may change with the extension of time, leading to the difference between the experimental characteristic frequency and the theoretical values.
Figure 8 shows the vertical variation in the characteristic frequency under different stratifications and cases, where the solid line is the fitted curve reflecting the variation in the characteristic frequency with depth.
Figure 8a–c show that the variation in the characteristic frequency with depth can be divided into two stages. In the first stage, the characteristic frequency increases with an increasing depth and then begins to decrease when it reaches a certain value. Then, there is a local minimum of the characteristic frequency. In the second stage, the characteristic frequency initially increases with an increasing depth and then decreases when it reaches a certain value. Some studies have shown that the characteristic frequency initially increases with an increasing amplitude. When the amplitude reaches a certain value, the characteristic frequency will stop increasing and begin to decrease [
46]. Since the amplitude first increases and then decreases with an increasing depth, such variation in the characteristic frequency with the amplitude is formed. The red line in
Figure 7 represents the depth where the maximum amplitude of ISWs lies (position of the maximum buoyancy frequency in the stratified fluid), which corresponds to the depth where the local minimum of the characteristic frequency lies.
Figure 8d shows that the characteristic frequency first increases and then decreases with an increasing depth under different cases in S4, and the characteristic frequency reaches the maximum near the depth of the maximum buoyancy frequency. This is mainly because the amplitudes in this stratification are relatively small. Therefore, the characteristic frequency only increases with an increasing amplitude in S4. In summary, for a single small-amplitude ISW, the characteristic frequency initially increases and then decreases with an increasing depth. The characteristic frequency is the largest at the depth where the maximum buoyancy frequency lies. For an ISW with relatively large amplitude, there is likely to be a local minimum of the characteristic frequency near the depth where the maximum buoyancy frequency is located.
4.3. Vertical Distribution of the Waveform
The ISW solution obtained by KdV theory is
. This ISW solution contains three basic elements (amplitude, characteristic frequency and
function). The amplitude and characteristic frequency of ISWs have been discussed above. To investigate whether the
function in KdV theory can describe waveforms of ISWs at different vertical depths, the amplitude and characteristic frequency obtained by experiments are used to calculate the theoretical waveform distribution of ISWs. Then, the theoretical waveform distribution is compared with the experimental results to verify the validity of the
function.
Figure 9 illustrates the theoretical and experimental results of the vertical distribution of waveforms in four different stratification conditions.
Figure 9 shows that in four different stratification conditions, the experimental distribution of solitary wave waveforms at different depths is consistent with the theory. Since some errors may occur when the conductivity probes are used to measure the density in the flume, the experimental waveforms of the ISWs at some depths deviate from the theory to some extent. Generally, the deviations are within the allowable range. Thus, the
function in KdV theory can describe the waveforms of ISWs at different depths well.
5. Conclusions
In this paper, four stratifications are established by using a gravitationally stratified fluid flume. The experimental results of the vertical structure characteristics of ISWs in each stratification are compared with the theoretical model, and the influence of different stratifications and different maximum amplitudes on the agreement between the theoretical and experimental results is analyzed. The conclusions are as follows:
Under different stratification conditions, the vertical amplitude structures of most ISWs with relatively small amplitudes are in good agreement with the theoretical model, whether it is elevation-type or depression-type. Additionally, with the increase in the maximum amplitude of ISWs, the influence of the nonlinear terms of the vertical structure function will be more significant, which leads to deviations between the experimental results and the theoretical model.
For the same stratification and the same wave amplitude, the theoretical value of the characteristic frequency in a two-layer fluid is larger than that in a continuously stratified fluid. For the depression-type ISWs, with the increase in , the theoretical values of the characteristic frequency in the two models are both larger, and the difference between the theoretical values of the two models is more significant under the same amplitude. The frequency–amplitude relationship in experiments indicates that the frequency–amplitude relationship can be described to a certain extent by the two-layer model or the continuously stratified model only when is relatively small.
The vertical distribution of the characteristic frequency of a single ISW in the experiment indicates that the characteristic frequency of a single small-amplitude ISW initially increases and then decreases with an increasing depth. The characteristic frequency is largest at the depth where the maximum buoyancy frequency lies. For an ISW with a relatively large amplitude, there is likely to be a local minimum of the characteristic frequency near the depth where the maximum buoyancy frequency is located.
Under different stratification conditions, the function of KdV theory can be in good agreement with the waveforms of ISWs at different depths regardless of whether the ISWs are depression-type or elevation-type.