The Focusing Waves Induced by Bragg Resonance with V-Shaped Undulating Bottom

: Intensive wave reﬂection occurs when the wavelengths of the incident waves and bottom undulations are in a 2:1 ratio. Existing studies have included the Bragg resonance phenomenon of waves passing over a continuous undulating bottom parallel to and oblique to the shoreline. More generally, the Bragg resonance mechanism is used as a means of coastal protection, rather than wave power generation. To focus the wave energy in a speciﬁc area, here, we propose sinusoidal sandbars of a horizontal V-shaped pattern, which is formed by two continuous undulating bottoms inclined at an angle to each other and the center axis of the angle is perpendicular to the shoreline. Based on the high-order spectral (HOS) numerical model, both the characteristics of Bragg resonance induced by the regular waves and random waves are investigated. In the scenario of regular waves, it shows that the wave-focusing effect is related to the angle of the V-shaped undulating bottom, and the optimal angle of inclination for the V-shaped undulating bottom is 162.24 ◦ . On that basis, considering the interactions between the random waves and the V-shaped undulating bottom of 162.24 ◦ , the Bragg resonance characteristics of random waves are studied. The BFI factor combining wave steepness and spectrum width can evaluate the focusing intensity of the Bragg resonance of the random waves. For BFI , in the range of 0.15–1.0, the values of H smax / H s0 linearly increase with the increase of BFI .


Introduction
Ocean waves contain large untapped renewable energy sources, which could be used to fulfil the energy demand [1]. Wave energy has the characteristics of high power density and wide distribution, and is renewable [2]. Hitherto, about one thousand wave energy converter (WEC) inventions have been patented [3], and more than 200 of them have entered the model testing stage [2]. However, present research mainly focuses on promoting the conversion efficiency of wave energy converters by improved mechanical design [4]. There are few research results to improve the wave energy generation efficiency from the perspective of increasing the wave energy density of the target sea area. Taking the characteristics of Bragg resonance into account is a good idea for wave energy focusing.
Bragg reflection or resonance is initially referred to as a special physical optical phenomenon of X-rays [5]. In the 1980s, Davies [6] first studied that Bragg resonance occurs if the bottom wavelength is an integral multiple of half wavelength of the incident waves. A followed idea originated for coastal protection was put forward. The strong reflection of water waves by seabed bars was investigated experimentally [7,8] as well as theoretically [9][10][11]. Afterwards, Bailard et al. [12] have extensively studied theoretically, computationally and experimentally to learn explicitly the feasibility of using longshore seabed bars to protect the shore and summarized its effectiveness and limitations. Some issues, including the partially standing-wave pattern [13] and shoreward increase 2 of 16 of wave amplitude [14][15][16] by Bragg reflection, were exposed. Howard & Yu [17] and Weidman et al. [18] have demonstrated that a specific phase relationship of wave-bottom interactions is significant, and sometimes it can cause exponentially varying standing-wave patterns. It is illustrated that Bragg resonance leads to offward intense reflection, even shoreward wave amplitude increase. Considering the strong effect of Bragg resonance, more studies focus on coastal protection or beach erosion from the aspect of the backward of the seabed bars.
Recently, Couston et al. [22] have proposed to revamp the Bragg resonance mechanism as a means of coastal protection by considering oblique seabed bars that divert, rather than reflect, shore-normal incident waves to the shore-parallel direction. Indeed, the novelty that with two superposed sets of oblique seabed bars reflecting waves along the shore-parallel direction, becomes efficiently deflected far to the sides, leaving a wake of decreased wave activity downstream of the patch. However, the wave characteristics of the forward of the oblique seabed bars are not studied by Couston et al. By contrast, few studies concern the wave reflection of the forward. From another perspective of the seabed bars reflecting waves, the wave-bottom interactions due to Bragg resonance may be of interest and be more in accordance with wave power generation need.
In the cases where closed-form solutions cannot be obtained, a wide range of numerical models are available to study water-wave scattering by seabed topographies [22]. These include the extended versions, e.g., the mild-slope equations [11,25], the coupledmode approach [26,27], the integral matching/discretized bottom method [28,29], the fully nonlinear Boussinesq equation adjusted for rapid bottom undulations [30], an asymptotic linear analytical solution (ALAS) [24], and the high-order spectral method [31]. Numerical investigations of Bragg scattering have helped explain several discrepancies between theory and experiments, including the difference between the observed and predicted Class II Bragg resonance frequency due to evanescent modes [32] and the resonant frequency downshift/upshift for the subharmonic/superharmonic Class III Bragg condition due to high-order nonlinearity [30].
Simultaneously, Bragg resonance theory has made some achievements in application fields. A series of submerged breakwaters have been developed to achieve coastal protection based on the Bragg resonance characteristic of reflecting waves to the sea in recent years. For the two-dimensional numerical cases, the parameters' influence on the Bragg resonance reflection coefficient has been investigated by scholars from different aspects, e.g., the shape, height, width, the number of the submerged breakwaters, the slope of the seabed, and wave conditions [33][34][35][36][37]. Generally speaking, these numerical simulations are based on various theoretical methods to study normally incident waves over a series of rectangular or trapezoidal breakwaters. Ning et al. [38] investigated numerically and experimentally the effect of the breakwater shape on hydrodynamic behavior in a two-dimensional flume. Shih & Weng [39] studied the interactions between waves and various combinations of undulating breakwaters in a three-dimensional basin. The reflection coefficient, transmission coefficient, and attenuation of wave energy were analyzed. Liu et al. [40] examined the Bragg reflection of waves by multiple submerged semi-circular breakwaters, by considering obliquely and normally incident waves independently. Although two-dimensional flume or three-dimensional basin experimental studies have been carried out, the study based on Bragg resonance for improving wave energy generation has not been widely considered. At present, research on the use of wave energy is gradually carried out. Water flow and wave propagation can be altered by artificial terrains or control devices. Zheng et al. [41] carried out a laboratory study on wave-induced setup and wave-driven current in a 2DH reef-lagoon-channel system. Elandt et al. [42] found that gravity waves can focus at a specific location by a concave mirror or a convex lens of gravity waves. Tao et al. [43] investigated experimentally strong wave reflection induced by Bragg resonance. They found that wave reflection is effective to amplify the free-surface oscillation amplitude and focus the wave energy in front of the undulating bottom. Zhang & Ning [44] confirmed that reflected waves from the parabolic opening of the breakwater can travel towards a fixed focus position. For a specified wave environment, the wave heights at the focus positions can reach over several times of the incident wave heights, indicating wave energy multiplication. They all proposed the concept of the focal point. Instead of putting a multitude of wave energy harvesting devices over a large area, a large wave energy absorber can be placed at the focal point. However, the influence of peculiar spatial bottom layouts is ignored, which also affect the focal point and focusing area. Previous physical experiment research [43] shows that waves are strongly reflected in front of the undulating bottom. Based on that, this study considers the wave-focusing effect in a specific area by optimizing the undulating bottom pattern.
The objective of this study is to investigate both the characteristics of Bragg resonance induced by the regular waves and random waves. Three main research questions need to be answered: (1) How does the angle of inclination for the V-shaped undulating bottom affect the wave-focusing characteristics induced by the regular waves? (2) What analytical methods are used to determine the optimal angle of the V-shaped undulating bottom? (3) How do the steepness and spectrum width of random waves affect the wave-focusing characteristics with the V-shaped undulating bottom, and what is the law of their influence on wave height? By the high-order spectral (HOS) numerical simulation, these questions have been examined. The concept of V-shaped undulating bottom is first proposed, i.e., a horizontal V-shaped pattern with two continuous undulating bottoms inclined at the same angle perpendicular to the shoreline. The V-shaped undulating bottom can not only play the role of reflecting waves in front of the undulations, but also exploit the advantage of symmetrical "V-shaped", which can focus waves on the central axis and increase the wave heights in the focusing areas. However, waves in the ocean are stochastic, and the random waves contain wave components of different frequencies, each of which has different Bragg resonance effect due to interactions with the undulating bottom of fixed frequency. Most of the previous studies on wave Bragg resonance only consider the regular waves at the dominant resonance frequency. In addition, there are few studies on Bragg resonance involving the interactions between random waves and undulating bottom. In this work, we have further studied the Bragg resonance characteristics of the interactions between random waves and V-shaped undulating bottom by the numerical simulation of HOS. This study has been organized as follows. In Section 2, the generalized Bragg conditions, the high-order spectral (HOS) method, and the model establishment are introduced. In Section 3, Bragg resonance reflection coefficients of regular waves are studied from the perspectives of different angles θ and different values f. In Section 4, the Bragg resonance focusing characteristics of regular waves and random waves for V-shaped undulating bottom are expounded. In Section 5, the conclusions are provided.

Methods
According to the high-order spectral (HOS) method [20], an effective computational method for general nonlinear wave-bottom (resonant) interactions and the boundary-value problem are stated.

Generalized Bragg Conditions
The mechanism for Bragg resonances is analogous to that for nonlinear (surface) wave-wave resonant interactions in the absence of bottom undulations. Thus, general Bragg conditions can be deduced from the well-known resonance condition for nonlinear wave-wave interactions [45]. For a wave field over uniform depth h, interactions among different wave components become resonant at order m (in wave steepness) if the wave numbers k j and the corresponding frequencies ω j satisfy: where the same combination of signs is to be taken in both equations, and k j and ω j satisfy the linear dispersion relation: Generalized Bragg resonance conditions in the presence of bottom ripples are obtained by replacing one or more of the free-surface wave components in Equation (1) by periodic bottom ripple components of corresponding wavenumbers k bj but with zero frequencies (since the ripples are fixed) [20]. Thus, by combining wavenumbers and frequencies of surface waves and bottom ripples, we obtain general conditions for Bragg resonances at each order, m = 2, 3 . . .
Here, the Bragg resonance between free wave surface and V-shaped undulating bottom originated satisfies generalized Bragg conditions. In addition, we propose that the first case considers two surface wave components of an incident wavenumber k 1 and a reflected wavenumber k 2 , propagating over a V-shaped undulating horizontal bottom containing a single wavenumber k b (analogously to surface waves, which refers to a fixed sinusoidally varying bottom with crest lines normal to k b ). The first case belongs to Class I Bragg resonance. The second case considers random surface wave generated by the Gaussian spectrum, propagating over a V-shaped undulating horizontal bottom containing a single wavenumber k b (analogously to surface waves, which refers to a fixed sinusoidally varying bottom with crest lines normal to k b ). In the paper, Class I Bragg resonance is considered, and numerical simulations all adopt the high-order spectral (HOS) method.

High-Order Spectral (HOS) Method
Combined with the advantages of Zakharov equations and mode-coupling ideas based on potential theory (with the application of the Fast Fourier Transform), the HOS method has been proved as a suitable, robust, and highly efficient numerical method for direct phase-resolved simulation of nonlinear ocean wave field evolution [46][47][48], wave nonlinear mechanism analysis [49][50][51][52] and nonlinear wave-wave interactions and wavebody interactions [31,53,54]. In this paper, the HOS numerical model for three-dimensional Bragg resonance is considered based on the HOS method for general nonlinear wavebottom interactions developed by Liu [20]. In addition, the focusing properties of the three-dimensional wave field caused by Bragg resonance with waves passing through V-shaped undulating bottom are studied.
The difficulty of establishing a numerical model is to deal with free-surface boundary conditions. To simplify the calculation, the high-order spectral (HOS) method introduces the velocity potential of the free surface in reference to Zakharov theory [55], so the freesurface boundary condition can be written in the following form:

Model Establishment
The V-shaped undulating bottom is shown in Figure 1. It shows that the two continuous undulating bottoms of a horizontal V-shaped pattern, a set of oblique undulations k b1 (the upper part of the V-shaped pattern). To match with k b1 , the other set of oblique undulations k b2 (the lower part of the V-shaped pattern matched) is given. The values of k b1 and k b2 are the same for a given V-shaped undulating bottom, while the difference of the k b1 and k b2 is in the directions of the two vectors. The vectors of k b1 and k b2 are symmetric about the x-axis. The color vectors in Figure 1 describe simply the V-shaped undulating bottom configuration, where the number of the continuous sinusoidal undulations N b = 10, an incident wave wavenumber k 1 (the direction of k 1 along the x-axis); two bottom undulations wavenumber k b1 and k b2 are perpendicular to crest lines of the each side of V-shaped undulating bottom respectively; the reflected waves k r1 and k r2 due to the two parts of bottom superpose each other ahead of the bottom. In addition, the relationship of the angle θ between k b1 /k b2 and x-axis and the angle α between the two sets of oblique crest lines satisfy: of oblique undulations kb2 (the lower part of the V-shaped pattern matched) is given. The values of kb1 and kb2 are the same for a given V-shaped undulating bottom, while the difference of the kb1 and kb2 is in the directions of the two vectors. The vectors of kb1 and kb2 are symmetric about the x-axis. The color vectors in Figure 1 describe simply the V-shaped undulating bottom configuration, where the number of the continuous sinusoidal undulations Nb = 10, an incident wave wavenumber k1 (the direction of k1 along the x-axis); two bottom undulations wavenumber kb1 and kb2 are perpendicular to crest lines of the each side of V-shaped undulating bottom respectively; the reflected waves kr1 and kr2 due to the two parts of bottom superpose each other ahead of the bottom. In addition, the relationship of the angle θ between kb1/kb2 and x-axis and the angle α between the two sets of oblique crest lines satisfy: 2. Page 7 of 15 in the jmse-09-00708-v2

Bragg Resonance Reflection Coefficients of Regular Waves under Different Values f
The wavenumber ratios f is determined as 0.9-1.1 due to the study of Liu and Yue [20]. For f, in the range of 0.9-1.1, 17 different wavenumber ratios f are uniformly selected to investigate the regular wave reflection coefficients with different values of f under the forward (θ = 0°) and oblique (θ = 19.5°) incident conditions. The results of the threedimensional HOS model and perturbation theory [19] are compared.  Table 1.

Results
In this study, the three-dimensional Bragg resonance numerical model is validated from two aspects: the first one is to verify the regular wave reflection coefficients under different incident angles θ; the other one is to verify the regular wave reflection coefficients under different wavenumber ratios f (f is the ratio of 2 times the incident wavenumber (2k 1 ) to the undulating bottom wavenumber (k b ), i.e., f = 2k 1 /k b ). This will prove the accuracy of the model in studying the focusing characteristics of Bragg resonance on V-shaped undulating bottom under different incident angles and different resonance wave numbers.

Bragg Resonance Reflection Coefficients of Regular Waves under Different Angles θ
The Bragg resonance reflection coefficient is the most important parameter to measure the magnitude of Bragg resonance. The Bragg resonance reflection coefficient is calculated according to the method proposed by Liu and Yue [20]. The Bragg reflection resonance coefficients R(θ) of the Class I Bragg resonance at different incident angles θ are simulated by the three-dimensional HOS method, as shown in Table 2, which are compared with the perturbation theory solutions of Mei [19] (oblique incidence on a finite strip of bars) in Figure 2.

Bragg Resonance Reflection Coefficients of Regular Waves under Different Values f
The wavenumber ratios f is determined as 0.9-1.1 due to the study of Liu and Yue [20]. For f, in the range of 0.9-1.1, 17 different wavenumber ratios f are uniformly selected to investigate the regular wave reflection coefficients with different values of f under the forward (θ = 0 • ) and oblique (θ = 19.5 • ) incident conditions. The results of the threedimensional HOS model and perturbation theory [19] are compared.
According to Figure 3, when f = 0.9-1.1, the Bragg resonance reflection coefficients increase first and then decrease with the increase of f. When f is slightly less than 1.0, the Bragg resonance reflection coefficients reach the peak, which also accords with the phenomenon of Bragg resonance frequency descending. The verification results confirm the reliability of the three-dimensional HOS numerical model. Figure 3, when f = 0.9-1.1, the Bragg resonance reflection coefficients increase first and then decrease with the increase of f. When f is slightly less than 1.0, the Bragg resonance reflection coefficients reach the peak, which also accords with the phenomenon of Bragg resonance frequency descending. The verification results confirm the reliability of the three-dimensional HOS numerical model.

Bragg Resonance Focusing Characteristics of Regular Waves for V-Shaped Undulating Bottom
In the section, the wavelength of the incident wave is 1 m, which need 16 nodes to represent. The other parameters of the incident waves and undulating bottom can be calculated according to the section of the Model Establishment.

Spatial Distributions of Wave Amplitudes in the Focusing Areas for Regular Waves
The V-shaped layout can be regarded as the combination of two parts of oblique bottom undulations relative to the incident waves. The wave-focusing effects, therefore, are the effects of two parts of oblique bottom undulations relative to the incident waves. When the wave incident angle (between k 1 and k b ) θ is greater than 45 • , R(θ) has a secondary maximum. There is a relationship between α and θ, as shown in Figure 1.
As θ increases, the α of the V-shaped layout becomes smaller. As a result, the wavefocusing area becomes smaller, so the wave-focusing effect is weakened. The values of α parameter for the specific simulation cases are listed in Table 1.
Bragg resonance occurs when regular waves pass over the V-shaped undulating bottom and satisfy the conditions that the wavelengths of the incident waves and bottom undulations are in a 2:1 ratio. Incident waves and reflected waves from two directions superimposed in front of the bottom causes wave energy focusing effect. Under the condition of regular wave incidence with the same wave amplitude, the Bragg resonance effects of wave-bottom interactions at different V-shaped undulating bottom angles α are simulated by the three-dimensional HOS numerical model. The maximum wave amplitudes calculated (A max ) are compared with the incident wave amplitude (A 0 ), and the spatial distribution characteristics of wave amplitudes are analyzed by A max /A 0 . The maximum wave amplitudes calculated (A max ) are different at different node positions. So, the values of A max /A 0 is also different. Table 3 shows that the calculated wave amplitude increases change at the focal points corresponding to different α (90 • ≤ α ≤ 180 • ). The node position of the maximum value of A max /A 0 for each simulation case is called the focal point, and then A maxp /A 0 is used to denote the value of A max /A 0 at the focal point. When α = 180 • (the incident wave direction k 1 is parallel to the bottom k b ), the calculated wave amplitude at the focal point increases to 1.81 times the initial incident wave amplitude. As θ increases, the α of the V-shaped layout becomes smaller. As a result, the wave-focusing area becomes smaller, so the wave-focusing effect is weakened. In conclusion, 90 • ≤ α ≤ 180 • is mainly considered in this study. To better analyze the spatial distribution characteristics of wave amplitudes, α = 178.   To better analyze the spatial distribution characteristics of wave amplitudes, α = 178.21°, 168.40°, 162.24°, 155.32°, 148.58°, and 140.47° are selected for comparison. The spatial distributions of wave amplitudes in the focusing areas due to Bragg resonance at different V-shaped undulating bottom angles α are analyzed, the values of Amax/A0 are shown in Figure 4.  In addition, the wave-focusing areas are generally expand from right to left in front of the bottom in a "V" shape that is slightly smaller than the "V" shape of the bottom. Figure 4 also shows that through the 6 selected representative α angles of the Vshaped undulating bottom, the spatial distributions of the wave amplitudes in the focusing areas owning to Bragg resonance are symmetric about the center of the V-shaped undulating bottom, and the largest wave amplitudes are along the axis of symmetry. By comparison, it is found that the wave-focusing intensity first increases and then weakens with the decreasing V-shaped undulating bottom angles α = 178.21°, 168.40°, 162.24°, 155.32°, 148.58°, and 140.47°. When α is about 162.24°, there are significant focusing areas and focal points, and the focusing intensity at this angle is the strongest.

Quantitative Analysis of Wave Amplitudes in the Focusing Areas for Regular Waves
To further quantitatively analyze and compare wave amplitude focusing characteristics at different angles α, the wave amplitude increases A max /A 0 are classified and counted. In addition, under the four levels of A max /A 0 > 1, A max /A 0 ≥ 1.5, A max /A 0 ≥ 2, and A max /A 0 ≥ 2.5, the occurrence frequencies of the values of A max /A 0 in the whole simulation range are calculated and expressed respectively in the form of P(A max /A 0 > 1), P(A max /A 0 ≥ 1.5), P(A max /A 0 ≥ 2) and P(A max /A 0 ≥ 2.5), respectively. The frequency formula of each level is as follows: where C is the specified level constant; N (A max /A 0 ≥ C) is the number of nodes corresponding to A max /A 0 ≥ C; N x and N y are the numbers of nodes on the x and y-axis in the whole simulation range respectively. In this section, C is assigned 1, 1.5, 2, and 2.5 respectively, N x = N y = 512.  Figure 5 shows the level-frequency statistical diagrams of wave amplitude changes at different angles α.
In Figure 5a, it shows that a series of nodes appear in the simulation range at all selected angles, but P(A max /A 0 > 1) and P(A max /A 0 ≥ 1.5) gradually decrease with the decrease of angles α. P(A max /A 0 > 1) decreases from 21.49% to 15.74%, while P(A max /A 0 ≥ 1.5) decreases from 5.32% to 0. It can be seen that the focusing area is decreasing. That is because the area of Bragg resonance decreases with the decrease of the angle α, and then the wave height focusing areas also reduce. In addition, that is why the range of angles α of V-shaped undulating bottom is 90 However, the wave-focusing effects are not only judged by the wave height in the focusing areas but also based on the overall amplitude increases of nodes in the wave height focusing areas, i.e., the values of A max /A 0 . Therefore, Figure 5b further shows the frequency statistic results of P(A max /A 0 ≥ 2) and P(A max /A 0 ≥ 2.5). It shows that both P(A max /A 0 ≥ 2) and P(A max /A 0 ≥ 2.5) increase first and then decrease with the decrease of angles α. Moreover, when 160 • < α < 168 • , the overall value of P(A max /A 0 ≥ 2) is larger. That is to say, the total number of nodes of A max /A 0 ≥ 2 is the largest. Meanwhile, when α = 162.24 • , the values of P(A max /A 0 ≥ 2) and P(A max /A 0 ≥ 2.5) are the largest in the simulation range, which implies that the number of nodes corresponding to A max /A 0 ≥ 2 and A max /A 0 ≥ 2.5 is the most, and the wave-focusing effect is the strongest. The spatial distribution of wave amplitude change shows that when α = 162.24 • , the scope of wave height focusing areas are relatively concentrated, the wave amplitude increases are the largest as a whole, and the wave-focusing effect is the best.
In conclusion, α = 162.24 • is considered to be the optimal angle for V-shaped undulating bottom.

The Gaussian Spectrum
The Gaussian spectrum is a symmetric spectrum, and its spectrum pattern is simple and regular. The expression of the Gaussian spectrum requires only 3 parameters, compared with the Jonswap formula, which requires 5 parameters. Therefore, it is more convenient to use the Gaussian spectrum to calculate. In this paper, the two-dimensional Gaussian spectrum pattern is selected, and the expression is shown in Equation (6).
where: Gaussian spectrum expression S(k i ), the standard deviation of the height of the wave surface η = ( S(k)dk) 1/2 = √ m 0 , zero-order spectrum moment m 0 , spectral peak wave number k 0 , spectral width parameter σ, value range of the i-th wave number 0 ≤ k i /dk ≤ N/4, the node number in the simulation range N.
The initial wave steepness and spectral width of the random wave are changed by controlling η and σ respectively. The initial wave steepness ε 0 , initial amplitude a 0 , and spectral width B 0 are calculated according to Equations (7)-(9) respectively.
The amplitude of the i-th wave can be calculated by Equation (10).
The initial wave surface and potential function can be obtained from Equation (11).
where: the initial wave surface η(x, 0), the initial potential function φ s (x, 0), the initial phase of the i-th wave generated randomly through the program θ i , the gravitational acceleration g, the circular frequency of the i-th wave ω i , according to the dispersion relation ω 2 = gktanhkh (water depth h) in the finite water depth to calculate the k i . Assuming that the initial wave surface and the initial potential function along the y-axis are equal, the two-dimension random wave field is extended into a three-dimension random wave field along the y-axis.

Evolution Characteristics of Random Waves on Flat Bottom
The Gaussian spectrum is used to generate the initial wave field. The initial phases of each simulated wave field are randomly generated. Different initial phases, dispersion relation, and wave modulation instability will cause uneven wave height distributions along the x-axis. Therefore, the evolution characteristics of random waves on the flat bottom are first studied in this study. The distributions of significant wave heights H s along the x-axis under three random initial phases in the wave order M = 3 are shown in Figure 6. It shows that the distributions of significant wave heights along the x-axis are greatly different with different initial phases. Here, the node number along the x-axis is 512, and 16 nodes represent one wavelength of incident free surface.
As the number of simulation increases, the significant wave heights H s under different random initial phases are averaged. The results are shown in Figure 7.
With the increase of the simulation groups, the distributions of significant wave heights H s along the x-axis become more uniform. When the number of simulation groups is more than 10, the significant wave heights are basically stable in the form of a horizontal line, which meets the requirements of analysis. Figure 6. It shows that the distributions of significant wave heights along the x-a greatly different with different initial phases. Here, the node number along the x-512, and 16 nodes represent one wavelength of incident free surface.  With the increase of the simulation groups, the distributions heights Hs along the x-axis become more uniform. When the number is more than 10, the significant wave heights are basically stable in the line, which meets the requirements of analysis.

Evolution Characteristics of Random Waves on V-Shaped Undulating
BFI was proposed by Janssen [56], which is determined by wave trum width. The initial definition is shown in Equation (12), which definition of wavenumber spectrum (14) by combining with Equatio where: the circular frequency  , the initial circular frequency 0  , incident random waves k , the wavenumber of initial incident random trum width B , the initial spectrum width 0 B , the initial wave steep Janssen proposed that when BFI is greater than or equal to 1.0, tions for generation of modulation instability, the possibility of freak To eliminate the influence of wave modulation instability on the si BFI is set at less than 1.0 in this study. Under different combinations ness and spectrum width, the values of BFI are listed in Table 4.

Evolution Characteristics of Random Waves on V-Shaped Undulating Bottom
BFI was proposed by Janssen [56], which is determined by wave steepness and spectrum width. The initial definition is shown in Equation (12), which is converted into the definition of wavenumber spectrum (14) by combining with Equation (13): 2∆ω/ω 0 = ∆k/k 0 (13) where: the circular frequency ω, the initial circular frequency ω 0 , the wavenumber of incident random waves k, the wavenumber of initial incident random waves k 0 , the spectrum width B, the initial spectrum width B 0 , the initial wave steepness ε 0 . Janssen proposed that when BFI is greater than or equal to 1.0, satisfying the conditions for generation of modulation instability, the possibility of freak waves is increased. To eliminate the influence of wave modulation instability on the simulation results, the BFI is set at less than 1.0 in this study. Under different combinations of initial wave steepness and spectrum width, the values of BFI are listed in Table 4. The calculated wave heights (H smax ) are compared with the incident wave height (H s0 ). To examine the interactions between random waves and undulating bottom, the spatial distribution characteristics of wave heights due to Bragg resonance are analyzed. Analyzing the combinations of initial wave steepness and initial spectrum width, it is found that there is a good linear relationship between H smax /H s0 and BFI, as shown in Figure 8. For BFI, in the range of 0.15-1.0, the values of H smax /H s0 linearly increase with the increase of BFI, so the fitting formula is shown in Equation (15): H smax /H s0 = 0.39 × BFI + 1.09 (15) found that there is a good linear relationship between Hsmax/Hs0 and ure 8. For BFI, in the range of 0.15-1.0, the values of Hsmax/Hs0 linea increase of BFI, so the fitting formula is shown in Equation (15) The goodness of fit (R 2 ) in Equation (15) can reach 0.94, indica fits well. It is speculated that the calculation results will be more co with the increase of the simulation groups. When the value of BFI is is a linear relationship between Hsmax/Hs0 and BFI; whereas the value Hsmax/Hs0 and BFI are weakly linear relation. When the BFI factor is s the initial wave steepness of a given spectrum width is small, wh nonlinear wave interaction is weak. However, a given initial wave s spectrum width, indicating that the frequency bandwidth of wave larger. In summary, the smaller the BFI factor is, the smaller Hsmax/H are becoming more independent, when the incident wave nonlinear

Conclusions
In this study, the V-shaped undulating bottom can be regarde of two parts of oblique bottom undulations relative to the incident w tributions of the wave heights in the focusing areas owning to Bragg wave-bottom interactions are symmetric about the center of the V-sh tom, and the largest wave height is distributed along the x-axis of sy two aspects of spatial distributions and quantitative analysis, it is about 162.24°, obvious focusing areas and focal points appear, and t at this angle is the strongest. By optimizing α, the wave-focusing effe 145° < α < 180°, the wave amplitude increase at the focal points can times of the initial incident wave amplitude, which is higher than th when α = 180°. In particular, when α = 162.24°, the wave-focusing eff maximum amplitude increases by 2.87 times of the initial incident w α < 145°, the maximum amplitude increase of the focal points is less When α = 90°, the maximum amplitude of the focal points is unchan focusing does not occur. The goodness of fit (R 2 ) in Equation (15) can reach 0.94, indicating that the formula fits well. It is speculated that the calculation results will be more consistent with the line with the increase of the simulation groups. When the value of BFI is more than 0.4, there is a linear relationship between H smax /H s0 and BFI; whereas the value of BFI is less than 0.4, H smax /H s0 and BFI are weakly linear relation. When the BFI factor is small, in other words, the initial wave steepness of a given spectrum width is small, which indicates that the nonlinear wave interaction is weak. However, a given initial wave steepness has a larger spectrum width, indicating that the frequency bandwidth of wave height distribution is larger. In summary, the smaller the BFI factor is, the smaller H smax /H s0 is; BFI and H smax /H s0 are becoming more independent, when the incident wave nonlinearity is low.

Conclusions
In this study, the V-shaped undulating bottom can be regarded as the combination of two parts of oblique bottom undulations relative to the incident waves. The spatial distributions of the wave heights in the focusing areas owning to Bragg resonance of regular wave-bottom interactions are symmetric about the center of the V-shaped undulating bottom, and the largest wave height is distributed along the x-axis of symmetry. From these two aspects of spatial distributions and quantitative analysis, it is found that when α is about 162.24 • , obvious focusing areas and focal points appear, and the focusing intensity at this angle is the strongest. By optimizing α, the wave-focusing effect is improved. When 145 • < α < 180 • , the wave amplitude increase at the focal points can reach more than 1.81 times of the initial incident wave amplitude, which is higher than that of the focal points when α = 180 • . In particular, when α = 162.24 • , the wave-focusing effect is the best and the maximum amplitude increases by 2.87 times of the initial incident wave amplitude. When α < 145 • , the maximum amplitude increase of the focal points is less than that of α = 180 • . When α = 90 • , the maximum amplitude of the focal points is unchanged, i.e., wave energy focusing does not occur.
Therefore, the V-shaped undulating bottom at the angle of 162.24 • is selected to study the random wave-bottom interactions characteristics. It is found that there is a good linear relationship between H smax /H s0 and BFI, which combines wave steepness and spectrum width. For BFI, in the range of 0.15-1.0, the values of H smax /H s0 linearly increase with the increase of BFI.