Gaussian Pulses over Random Topographies for the Linear Euler Equations

. This study investigates numerically the interaction between a Gaussian pulse and variable topography using the linear Euler equations. The impact of topography variation on the amplitude and behavior of the wave pulse is examined through numerical simulations and statistical analysis. On one hand, we show that for slowly varying topographies, the incoming pulse almost retains its shape, and little energy is transferred to the small reflected waves. On the other hand, we demonstrate that for rapidly varying topographies, the shape of the pulse is destroyed, which is different from previous studies.


INTRODUCTION
The study of the ocean holds significant interest in its own right.The main set of equations used to model the dynamic of water waves was first derived by Isaac Newton (1687) with limited mathematical rigor.Later, Leonhard Euler (1757) derived a set of equations for the governing water wave dynamics without incorporating any dissipative effects, these set of equations become known as the Euler equations [7].
The Euler equations are used to model several problems such as ship wakes, flow of water over rocks, the formation of storms in the ocean and even atmospheric problems such as atmospheric flows encountering obstacles [2,14].These problems are typically approached as deterministic, where the initial state is precisely known, and the evolution of the phenomenon over time can be computed using the Euler equations.However, in many cases, the scarcity of data for initializing deterministic models presents a challenge.In such instances, assuming random initial states becomes preferable, allowing for the study and formulation of statistically reliable predictions [8,9,21,22,25].Analogously, the medium in which a wave propagates can also be considered random [1,12,13,20].In particular, we mention the work of Nachbin and Solna [20] which investigates the effect of a random topography on the Gaussian pulse on the free surface.Among their results, they showed that the pulse experience an apparent diffusion, which is featured by a loss in amplitude of the incoming pulse after passing over a rapidly varying topography.
The Euler equations have been widely utilized in various studies.For example, Choi and Camassa [5] numerically computed solitary waves over a flat bottom , while Choi [4] investigated these waves in the presence of a vertically sheared current with constant vorticity.Ribeiro-Jr et al. [23] analyzed particle trajectories and pressure within the fluid bulk.Nachbin [16] introduced a conformal mapping technique to flatten uneven topography and derived an asymptotic Boussinesq-type equation with variable coefficients.Flamarion et al. [10] examined particle trajectories for the linear Euler equations considering a linear sheared current and variable topography.Moreover, the Euler equations have found applications in the study of free boundary problems at small scales, such as the motion of a fluid droplet (silicon) walking on the surface of a vibrating bath.This phenomenon exhibits similarities to the wave-particle duality observed in quantum mechanics [3,6,15].Nachbin and colleagues extensively explored the role of free surface waves as a spontaneous mechanism for synchronizing oscillating droplets confined within various types of cavities, through the Schwarz-Christoffel mapping [17,18,19].
The objective of this study is to investigate the behaviour of the interaction between a Gaussian pulse with a random topography.As our model is deterministic, we conduct numerous Monte Carlo simulations with randomly generated topographies sharing identical statistical characteristics.The primary aim of this study is to investigate statistically the effective behavior resulting from the interaction between an initial Gaussian pulse and randomly smooth topographies, an approach not fully explored in the literature.
Our study indicates that when the topography exhibits slow variation (small wavenumbers), the resulting wave shape undergoes minimal changes, with the amplitude remaining nearly constant.As the wavenumber increases, fluctuations in the resulting wave amplitude emerge, yet its shape still resembles the initial Gaussian bell, albeit with additional tails caused by the pulsetopography interaction.Ultimately, as the wavenumber further increases, the resulting wave experiences a change in shape, tending to split into two pulses.This division results in a substantial difference in amplitude between the leading wave and the initial pulse.This article is organized as follows: In Section 2, we present the mathematical formulation of the Euler equations.In Section 3, we describe the conformal mapping technique and rewrite the Euler equations in the canonical domain, which is a uniform strip.The numerical methods are presented in Section 4. In Section 5, we present the main results, and the final considerations are discussed in Section 6.

MATHEMATICAL FORMULATION
We examine the behavior of a two-dimensional, irrotational flow in a finite depth channel filled with an inviscid and incompressible fluid of constant density.The channel includes a variable topography and features a free surface wave (ζ (x,t)) affected by gravity.The governing equations, Trends Comput.Appl.Math., 25 (2024), e01766 known as the Euler equations and referenced in [26], describe the dynamics of the free surface and the potential velocity ( φ (x, y,t)).These equations can be expressed in dimensionless form as where ε is the nonlinearity parameter, which measures the ratio between a typical wave amplitude and the depth channel, and µ is the shallow-water/long-wave parameter which measures the ratio between the depth channel and a characteristic wavelenght.For computational purposes, it is convenient to consider the new scaling y → µy.In this new scaling, the potential velocity is a harmonic function, thus the Euler equations become φxx + φyy = 0 for µh(x) < y < ε µζ (x,t), We are interested in investigating solutions of (2.1) in the linear regime (ε = 0) and weakly dispersive regime (µ ≈ 0).Under these assumptions we have the linear system of equations φxx + φyy = 0 for µh(x) < y < 0, φy = µ φx h x at y = µh(x), When the bottom is flat, the phase speed of the linear waves are given by [26] We recall that this is for example the speed of the crest of the wave as it propagates.Notice that in the limit µ → 0, all linear waves travel with the same speed c = 1.

CONFORMAL MAPPING AND NUMERICAL METHODS
We remember that one of the primary characteristics of conformal mappings is the conformal invariance of the Laplacian operator.Consequently, to solve equation (2.2), we can employ a conformal mapping and address it within a simpler domain.For this purpose, we construct a conformal mapping denoted as which flattens the bottom topography and maps a strip of width D onto the fluid domain, where the value of D is determined afterwards.This conformal mapping satisfies the boundary conditions y(ξ , 0) = 0 and y(ξ where The formulation presented here is identical to the one described in [11], hence we will briefly summarize the main steps. Let φ (ξ , η,t) = φ (x(ξ , η), y(ξ , η),t) represent the potential velocity in the canonical domain.
By utilizing the Dirichlet-to-Neumann condition, we obtain Furthermore, at η = 0, we have where in the shallow-water regime (µ ≈ 0), we have the limit Therefore, all the topography information is carried in the metric coefficient M(ξ ).

NUMERICAL METHODS
We solve (3.2) and (3.3) in a computational domain [−L, L], with a uniform grid with N ξ points and step ∆ξ = 2L/N ξ .The spatial derivatives are computed using the Fast Fourier Transform (FFT) [24].In addition, the time advance of the system (3.3) is computed through the Runge-Kutta fourth order method (RK4) with time step ∆t.The numerical method proposed here is implemented in MATLAB and it consists of a variation of the code available in [11].
The initial wave profile is taken as where ξ 0 is its initial location.The initial potential velocity is taken as where c(k) is given by (2.3) and F −1 denotes the inverse Fourier transform.This choice of initial data guarantees that the solution travels to right when the bottom is flat.A zero-mean random topography is modeled as the Fourier series and contains N harmonics where χ [−L 0 ,L 0 ] (ξ ) is the characteristic function in the interval [−L 0 , L 0 ], S(k) is an initial power spectrum, k i = i∆k and ∆k is a sampling wave number, and a phase ϕ i is a random variable, uniformly distributed in the interval (0, 2π).The length of initial realization is L = 2π/∆k.An initial power spectrum is assumed to have a Gaussian distribuition In what follows we set the following parameters: ∆k = 0.013, L = 2π/∆k, N = 256, N ξ = 2 12 , ∆t = 0.005, σ = 0.27, Q 0 = 0.01, ξ 0 = −60 and L 0 = 50, µ = 0.001.

RESULTS
We examine various values for the peak wavenumber, denoted as k 0 , which play a crucial role in determining the rate of topography variation.When k 0 ≫ 1, the topography is said to be rapidly varying, while for k 0 ≪ 1, the topography is said to be slowly varying.Consequently, we focus our simulations on investigating three distinct regimes of k 0 , namely, intermediate varying, rapidly and slowly.We begin by considering the case where k 0 = 0.5.In this scenario, the Gaussian-shaped wave pulse propagates without significant reflections or dispersion effects in the positive direction.However, when the wave encounters the intermediate varying topography, its amplitude undergoes variations, leading to the generation of reflected waves.This phenomenon is illustrated in detail in Figure 1, which demonstrates the interaction between the wave pulse and the topography.The interaction with the variable topography causes the bell-shaped pulse to lose energy, resulting in a decrease in its amplitude.This energy is then transferred to the reflected waves propagating in the opposite direction.To gain a deeper understanding of this behavior, we conducted numerous Monte Carlo simulations to draw statistical conclusions regarding the pulse-  topography interaction.Figures 2 and 3 provide detailed information on the maximum amplitude variation of the pulse after interacting with the topography and throughout the entire simulation, respectively.These figures highlight the statistical nature of the amplitude fluctuations, with the standard deviation of the amplitude on the order of O(0.1) following the interaction.The influence of the topography on the amplitude-wave behavior becomes evident.To further explore the statistical aspects of the pulse-topography interaction, we examine the probability distribution of the output wave amplitude during the interaction.Figures 2 (right) and 3 (right) display the probability distribution of the output wave maximum amplitude and the maximum amplitude during the interaction, respectively.These probability distributions provide valuable insights into the behavior of the pulse as it interacts with the variable topography.For instance, in Figure 2 (left) we see that the probability of obtaining an output wave with amplitude greater than 0.95 is approximately 60%.
Next, we investigate the case where k 0 = 2.In this regime, the resulting pulse shape undergoes a significant transformation upon interacting with the topography, as depicted in Figure 4.The interaction between the Gaussian pulse and the topography becomes more pronounced, resulting in wave reflection.In order to give more details, we monitor the changes in the maximum amplitude of the output wave and the wave during propagation, as shown in Figures 5 and 6, respectively.It is important to note that compared to the previous case, the fluctuation of the output wave amplitude is more prominent.This increase in fluctuation is a direct consequence of the stronger interaction between the pulse and the topography.The reflection of waves caused by this interaction introduces additional complexity to the system, resulting in larger amplitude variations in the output wave.Notice that differently from the previous case we see that the probability of obtaining an output wave with amplitude greater than 0.90 is zero, see Figure 5 (left).Lastly, we consider the case where the topography varies slowly (k 0 is small).We set k 0 = 0.05. Figure 7 displays a typical wave solution at different times for slowly varying topography.Observe that its amplitude remains almost unchanged after the Gaussian pulse interacts with the topography.Additionally, we note that the reflected waves have less energy compared to the previous cases.This behavior is also evident from the statistical perspective, as shown in Figure 8 for the output wave.The fluctuation in the wave amplitude is very small, indicating that its initial shape is almost preserved.The same holds true for the fluctuation of the Gaussian pulse as it Trends Comput.Appl.Math., 25 (2024), e01766 passes over the topography, as seen in Figure 9.It is worth to mention that the probability of the output pulse retains 99% of its initial amplitude is about 80%.In summary, our investigation demonstrates that different values of the peak wavenumber k 0 , which results in different topographies, significantly affecting the amplitude and behavior of the wave pulse.These results ties with the study of Nachbin and Solna [20].We notice that the changing in the dynamic presented in this work is different from the ones reported by Flamarion and Ribeiro-Jr [11].In this article, the authors showed that a monochromatic rapidly varying topography barely affects the shape of the pulse, while a slowly varying topography causes a considerable variation in amplitude and an increasing of reflected waves.This occurs because only monochromatic topographies are considered in [11].Therefore, our study shows that when the topography is modeled by a series of Fourier modes the wave field changes considerably.The statistical analysis of the pulse-topography interaction sheds light on the energy transfer and amplitude fluctuations involved in the process.

CONCLUSIONS
In conclusion, our investigation into the interaction between wave pulses and variable topography has provided valuable insights into the behavior of the system.On one hand, we show that for slowly varying topographies, the incoming pulse almost retains its shape, and little energy is transferred to the small reflected waves.On the other hand, we demonstrate that for rapidly varying topographies, the shape of the pulse is destroyed.This article contributes to a deeper understanding of wave propagation in variable topographic environments.Further studies can be built upon this work to explore additional aspects of wave-topography interactions for more complex models.

Figure 1 :
Figure 1: The Gaussian pulse evolution at different times for the wavenumber peak k 0 = 0.5.

Figure 2 :Figure 3 :
Figure 2: Left: The maximum amplitude of the pulse after passing over the topography (time t = 150).Right: The distribution of the amplitudes over 1000 realizations for the Gaussian pulse at t = 150.Parameters: k 0 = 0.5.

Figure 4 :
Figure 4: The Gaussian pulse evolution at different times for the wavenumber peak k 0 = 2.

Figure 5 :Figure 6 :
Figure 5: Left: The maximum amplitude of the pulse after passing over the topography (time t = 150).Right: The distribution of the amplitudes over 1000 realizations for the Gaussian pulse at t = 150.Parameters: k 0 = 2.

Figure 7 :Figure 8 :
Figure 7: The Gaussian pulse evolution at different times for the wavenumber peak k 0 = 0.05.

Figure 9 :
Figure 9: Left: The maximum amplitude of the pulse over time.Right: The distribution of the amplitudes over 1000 realizations for the Gaussian pulse over time.Parameters: k 0 = 0.05.