KdV Equation Model in Open Cylindrical Channel under Precession

A new model for Korteweg and de-Vries equation (KdV) is derived. The system under study is an open channel consisting of two concentric cylinders, rotating about their vertical axis, which is tilted by slope {\tau} from the inertial vertical z, in uniform rate {\Omega}_1={\tau}{\Omega}, and the whole tank is elevated over other table rotating at rate {\Omega}. Under these conditions, a set of Kelvin waves is formed on the free surface depending on the angle of tilt, characterized by the slope {\tau}, volume of water, and rotation rate. The resonant mode in the system appears in the form of a single Kelvin solitary wave, whose amplitude satisfies the Korteweg-de Vries equation with forced term. The equation was derived following classical perturbation methods, the additional term made the equation a non-integrable one, that cannot be solved without the help of numerical methods. Invoking the simple finite difference scheme method, it was found that the numerical results are in a good agreement with the experiment.


Introduction.
Most wave phenomena in modern science are decribed by nonlinear equations. The most beautiful model was the one derived in 1895 by two Dutch mathematicians, which the equation takes their names: Korteweg and de-Vries (KdV) equation 2 . Although it was first dedicated for water waves in open channel, after the discovery of solitary wave by Russell (1834), later in the sixties the equation was noticed in other applications like: internal gravity waves in stratified fluid, ion-acoustic waves in plasma, axisymmetric waves in a nonuniformly rotating fluid, nonlinear waves in cold plasmas, axisymmetric magnetohydrodynamic waves and several other physical applications. (Debnath 2012).
In this paper I introduce a new model of this KdV equation with azimuthal dependence where the cylindrical geometry, incorporates both tilt and rotation effects. This is a new case, as all the previous cases that discussed the solitary wave invoked the rotation effect only. The rotation 1 (E-mail: hajaralshoufi@ymail.com) 2 The KdV equation was also derived by Boussinesq in the 1870's. He found its first three conservation laws, and its one-soliton and periodic traveling wave solutions, see Equation (30), page 77, of the paper J. Boussinesq, Theorie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond, J. Math. Pures Appl. 17 (2) (1872) and Equations (283,291) of the work J. Boussinesq, Essai sur la theorie des eaux courants, Mem. Acad. Sci. Inst. Nat. France 23 (1) (1877) 1-680. This historical observation is due to Prof. P.J. Olver (~ 2001). effect, however, is old enough as the internal solitary waves are the famous example that may occur in lakes, oceans, fjords, etc. The earliest experiment was carried out by Maxworthy (1983) to study the rotation effect on the internal solitary wave in a rectangular tank mounted over a rotating table, he noticed that the rotation affected the wave form as it was noticed to be curved backward while moving. He called this wave the internal Kelvin solitary wave, that generated when the pressure gradient normal to the lateral boundary opposes the Coriolis effect generated by the wave motion. The author finally confirmed that the wave speed is independent of the rate at which the system rotates and depends only on the stratification and maximum wave amplitude. A similar experiment was carried out later by Renouard et al. (1987) with approximately similar observations of Maxworthy (1983). Their observations about the wave showed that the crest of the wave is neither horizontal nor contained in a vertical plane perpendicular to the side due to rotation effect, but is curved backward, hence there is a spatial phase shift which increases by increasing the distance from the wall, and at a given distance from the wall increases by increasing Coriolis parameter.
Deep mathematical studies about rotation effect was thoroughly discussed by Grimshaw (1985), who proposed two different models discussing whether the rotation effect is strong or weak. His results showed that if scaling the transverse direction with where −1 the dimensionless measure of Rossby radius, if it has either unity value or of order (the shallowness parameter) i.e., strong rotation conditions, then the amplitude of the linear internal Kelvin satisfies a Korteweg-de Vries equation for shallow fluids at the wall, or its counterpart for the deep fluids. The second case, when weak rotation occurs, the scaling of the transverse direction is of order 2 , then the effects of rotation are small and comparable to those of nonlinearity and dispersion, in this case the transverse variation of wave amplitude is undetermined at the leading order, and is instead contained in the evolution equation, which is a modification of the two-dimensional Korteweg-de Vries equation for shallow fluids, or its counterpart for deep fluids. Although the solitary wave is of the stable type, that does not change its shape while moving, Grimshaw et al. (1997) studied the extinction of KdV solution in finite time due to the presence of low-frequency rotational term (additive term to KdV). The extinction time was 2√ 0 / 2 , where 0 the initial solitary wave amplitude, and = 2 2 (with long wave velocity = √ ℎ, and Coriolis parameter = 2 sin( ), the gravity acceleration, ℎ the depth of the basin, Earth rotation frequency, and geographic latitude location) is responsible for low frequency dispersion, which typically describes the effect of rotation in a physical system. They noted that this extinction time could be as important as frictional decay times for internal waves in the ocean.
It was shown that cylindrical symmetric flows can support stationary waves of finite amplitude under certain circumstances. (Leibovich 1970). This encouraged both Derzho and Grimshaw (2002) to study the structure of steady travelling solitary waves of large amplitudes in axially symmetric rotating flows in a long tube including the recirculation zone effect. There was critical amplitude where the flow reversal (recirculation) first occurred. The critical amplitude was related to the incipient rotation and the lowest mode of boundary problem solution at the lowest order. It was noticed that the small parameter that determined the critical width of recirculation area could be related to the amplitude parameter that determined the nonlinearity of the problem, and it took two extreme values, at the minimum level the wave was moving out of the recirculation zone with the normal solitary form, and the recirculation zone reaches a minimum radial extent . And at the maximum level the width of the recirculation zone became infinite, and the crest of the wave looked flattened. In fact, the flow in tubes under rotation conditions was also studied when some obstacle in the direction of the flow exists. Grimshaw (1990) proposed the theory for this case, precisely when upstream flow is nearly resonant (critical) defined with respect to long waves mode, that has a phase celerity almost zero. They used the forced KdV equation for the simulation with forcing term related to the obstacle, and advection term that contained the detuning effect. Later Grimshaw and Zhu (1994) completed the same problem for the unbounded radial swirling flow, they showed that this wave equation is similar to another one derived by Leibovich for freely propagating weakly nonlinear waves on a radially unbounded swirling flow.
The numerical methods that processed Kortweg-de Vries equation under rotation effects are varied, for instance  prepared a numerical scheme of the modified Kadomtsev-Petviashvili (KP) equation to investigate the nature of the waves generated from typical initial conditions. The KP model was solved numerically where the trapezoidal rule was used for performing the integrals, while the equation was discretized through a simple explicit second-order Lax-Wendroff finite-difference scheme. Their results suggested that weakly nonlinear inviscide theory actually revealed the main features of the observed waves, and that wave front curvature is possible because of slow attenuation of the wave amplitude. In their later work similar numerical scheme using Lax-Wendorff finite-difference scheme was used, but this time for the forced KP equation, for the sake of the dependence of soliton radiation on channel width. This equation was derived under conditions of pressure distribution travelling at a speed near the linear-long-wave speed, and it had three-dimensional formulae to express the effects of nonlinearity, linear dispersion, and three-dimensions close to critical speed. The linearized KP equation showed that no steady state is reached at critical conditions. Motivated by studying the dynamics of internal Kelvin waves, Melville et al. (1988) studied the effect of nonlinear resonance with Poincaré waves on the stability of weakly nonlinear dispersive Kelvin waves of finite bandwidth, this resonance arose from the nonlinear correction to the speed of the Kelvin wave. Their numerical analysis using Fourier transformation showed that the leading-order Kelvin wave solution in the horizontal direction was unstable, owing to resonant forcing of the transverse velocity associated with the linear Poincaré modes of the channel.
Motivated by studying the nonlinear plasma waves, a group of Chinese researchers Peng et al.
(2019) carried out the derivation of KP equation in cylindrical geometry, where the effect of ions' densities and their temperature were noticed in the coefficients of the final CKP (cylindrical KP equation). This equation was solved numerically to study three different types of waves: breather, rogue waves, and solitons. The authors invoked the Grammian solutions to solve the soliton problem, and the bilinear Hirota method to solve the breather waves. And finally the rogue equation that was treated in the limit of long wave function that is assumed an exponential polynomial. Their final result showed that the amplitude of rogue waves decreased with the absolute value of time, thus the energy of the wave was damped with the increase of time value. Another novel method used to solve the extended modified KdV equation, using the ̅ -dressing methods, is proposed by Li &Tian (2021). The authors introduced a spectral transform matrix to construct the spatial-time spectral problems of the nonlocal emKdV based on Cauchy-Green integral operator. And then to find the hierarchy of nonlinear evolution equations they used a recursive operator based on the proper spectral transformation matrix, so that the final N-soliton solution is derived and solved.
The new techniques in solving the nonlinear equations are varied, like the inverse scattering transform, the Lax Pairs, the Darboux transformation. However, Hirota bilinear method is of high interest, it was used to solve many type of KdV equation. For instance the form of soliton breaking was solved by Hossen et al. (2018), where they discussed the interaction with rogue waves, by introducing different quadratic polynomials that include either ℎ term or exponential term, and solved the problem accordingly. The ( A relatively new case of the solitary wave under rotation conditions was examined by Abderrahmane et al. (2011), this wave was observed during the liquid drainage through a centrally located circular opening on the bottom of the container, when shallow water conditions were reached. The free surface first displayed long wavelength increases as the water height decreased, then when the shallow water conditions were reached the free surface undulations transformed almost instantly into a surface bulge or solitary wave rotating around the cylindrical wall of the container. Their case here was treated under swirl conditions using full Euler equations, but by taking the wave variation dominantly along the azimuthal direction of the flow and the radial motion was neglected. In later work by Amaouche et al. (2013), the same case study of the wave in the cylindrical tank was treated in different model contained all velocity components plus Coriolis effect. In this case the linearized solution of the boundary value problem leaded to a modified Bessel equation that included rotational effects, and in both cases the material derivative was considered as a function of time and the azimuthal coordinates from which they included the rotational effect of the cylinder itself. Even if the dimensionless scaling was taken into consideration, the normal shallow water speeded up.
In open channels one may find simple linear waves, as well as other nonlinear ones, this includes the weakly and strongly nonlinear waves. Those waves can be generated by navigation, or the operation of hydraulic controls or accidentally, then all waves happened to occur because of them are forced ones, and one the "great primary wave", i.e., the solitary wave, is forced wave that accompanied a canal boat in its motion. (Darrigol 2008). Such phenomena affect the hydraulic conditions in the channel and pose engineering problems. The violent impacts of water waves on walls create velocities and pressures that are bigger in magnitude than the normal ones, which may increase the erosive action and warrant the reinforcement or redesign of the channel. Also, one should not forget the dynamics of sediment transport processes generated by the action of the flow phenomena. Thus, this work in the present paper is motivated by studying the forced oscillations and their instabilities in open channels. Precisely, studying the solitary wave which is carried out here theoretically by deriving the new KdV version, and experimentally by comparing the mathematical model with the experiment. This wave is also called the single Kelvin fundamental resonance mode. As the interested reader can see the previous work (cfr. ALSHOUFI 2021) about the forced oscillations that were studied linearly first, it turned out that when the sytem resonates the wave becomes singularly large and the linear analysis is not valid anymore. Thus, to complete the linear analysis in this work I introduce the weakly nonlinear case where this specific wave is under consideration, because it has similar features to the famous solitary wave first observed by Russell (1843). The paper is divided as follows: in §2 The governing equations where the new system of coordinates after transformations from the laboratory frame to frame fixed on the tank is introduced, in addition to flow assumptions, in §3 experiment and results where a description of the design is carried out with comparing the results experimentally, in §4 conclusion.

Governing Equations.
The motion of the wave is supposed to occupy in cylindrical annulus that is confined between the inner cylinder with radius , and the outer cylinder with radius . This distance does not exceed 10cm, and it was deduced by the work of  that in such small geometries if the width 2 smaller than specific ratio which is ≤ ℎ ̅ 2 / , where ℎ ̅ the average depth of water, the typical wave amplitude, the three-dimensional effect can be negligible. This ratio was computed in the designed flume for relatively big amounts of water inside it, those are the cases where the solitary wave was observed. Table (1) presents those results. Thus, the motion is assumd two-dimensional one, where the predominant direction is the azimuthal one, as water was the main liquid fluid used in the experiment, then the flow is assumed incompressible and inviscid one. Additional assumption was on replacing velocity components (which are , in the azimuthal and axial directions, respectively) with scalar potential function ɸ, and assume the flow irrotational one. As the surface of the wave is under the atmospheric conditions then the pressure is assumed negligible. The flume was filled with different volumes of water that are characterized with their average depth ℎ ̅ . The evolution of the free surface function is a time-space depenent function ( , ). Any solitary wave appear on the surface has an amplitude . In this work we assume that the total mechanical energy of the flow is conserved in the fully developed flow, thus, turbulence, wall friction due to boundary layers, internal friction due to viscosity, capillarity due to the surface tension and dissipation of vorticity are all connected to energy dissipation and will be neglected. As mentioned that the channel is tilted from the axial direction, this is characterized by the slope of the upper tilted table . (cfr. §3 Experimental Procedure). As the solitary wave is a long wave, thus the platform in studying this wave is the shallow water theory, this theory depends on the hydrostatic pressure distribution, and one important feature of this theory is that no assumption is made regarding the magnitude of surface elevation or the velocity components. (Stoker1957). Thus, the parameter that describes the evolution of waves in shallow water is called the shallowness parameter and it is going to appear in the system of equation when introducing the scaling process for the dimensionless analysis. This is also called the dispersion parameter : where = 2 represents the length scale or the horizontal scale, which is assumed to be the periphery of the outer cylinder. In the experiment this ratio was all the time satisfied, as it is clear in Table (1). As usual -axis is taken vertically upward, the bottom is given by: = 0, while the free surface as: = ℎ ̅ + ( , ), the flow is confined between the inner and the outer radius of the channel. The boundary conditions on rigid walls like the bottom satisfy the impermeability condition so that the normal velocity is zero: = 0 at = 0. At the free surface of water, the kinematic condition is given as: We can now write Euler equations as the following: where represents the material derivative operator: The conservation of mass or continuity equation can be written as: The condition for irrotational flow is worth to be split as follows: The treaties provided by Korteweg-de-Vries (1895) showed that to find the nonlinear type of wavy motion represented by KdV equation, there should be other important variable to measure the nonlinearity is called the amplitude parameter, and it is going to be included in the process of dimensionless analysis as well, it has the form: Thus, the unknown quantities of velocity field components, the pressure and water surface are first going to be scaled as follows to include the shallowness effect: Where the long wave speed is given by: = √ ℎ ̅ . Where we assume that = ℎ ̅ . After substitution into the major equations of continuity, momentum, and bounday conditions, we find: with the irrotational condition: 2 ̅̅̅̅ ̅ = ̅̅̅. In effect equations (7), (8), (9), (10) represent the shallow water equations at the leading order of ( ), as it is obvious from (9) that this scaling leads to the hydrostatic pressure distribution, the main assumption of shallow water theory. One can also notice at the leading order that the effect of both Coriolis and Euler forces in equations (8) and (9) are of second order of the shallowness parameter ( 2 ), thus, they vanish and only pressure and gravity act. The KdV derivation can be further pursued depending on equations (7), (8), (9), and (10) by taking higher orders of the factor this was done by Keller (1947) for instance: it leads to fully nonlinear solution. However, we can instead make it easier, especially since our equations are more complicated and here comes the usage of the amplitude parameter : We find after rescaling the following: Equation (14) is the free surface kinematic condition that is applied at: ̅ = 1 + ̅ . KdV equation has beautiful balance between quadratic nonlinear term and the linear dispersive one (Grimshaw 2007), this can be invoked by assuming: Or equivalently for simplicity: So, we find after substitution: Before expanding in terms of the amplitude parameter for the unknown quantities, it is better to integrate the equation replacing the nonlinear terms with the corresponding irrotational versions, and also by using the scalar velocity potential function for each velocity projection: ̅ = ɸ , ̅ = ɸ , so that we get the Forced Bernoulli Equation in the system, where the forcing term is the integral of azimuthal gravity force projection radially as in (19): Now the expansion of the unknown quantities (the amplitude and velocity potential) will be in terms of the amplitude parameter and takes the form: ɸ = ɸ (0) + ɸ (1) + 2 ɸ (2) + ⋯, ̅ = ̅ (0) + ̅ (1) + 2 ̅ (2) + ⋯.
Let's re-write the equations, mainly the continuity, Bernoulli, the kinematic and boundary conditions, respectively. We get: The unified form of the phase of tilt , with the azimuthal angle in the forced term can be harnessed in favor of d'Alambert wave solution, as in the first order of the problem we have: This will give the linear long wave equation: which has d'Alembert solution in the dimensional form as: where here is assumed the mean azimuthal velocity that the wave in the channel flume is assumed to rotate with about the outer periphery, thus it has the form: = .
Thus, in the dimensionless form of the problem we can write: Hence, we can introduce the following: We also assume slow time evolution for the wave by introducing: Substituting back into Bernoulli and continuity, and dropping the bars we get: Equations (22), (23), (24) are going to be expanded in terms of . It is worth to mention up to this point that the tilt angles in this system are in general small ones, and that the cases where the solitary wave were noticed corresponded to very small tilt character factor: , (cf. Table 1), thus it is assumed that if → 0 (which is the solution we are interested in for small wave amplitude) then → 0, so that at the leading order of the problem the forced term can be neglected. The most important result at the leading order then is the fact that the connection between potential velocity and the amplitude is the dimensionless rotation speed: Collecting similar terms in higher orders will lead to the final version of KdV equation in the channel under precession conditions, is given as: This equation is KdV type which unfortunately does not have integrable solution and the usage of numerical method is carried out in the next section. Of course, the form in (26) is dimensionless and using the assumed scaling the dimensional version is recovered to be compared with the experiment. The final version in the dimensional form takes the form: where the coefficients A, B, C, D, are respectively given by: If we assume periodic initial condition for as: then the forced term can be approximately assumed to be: − sin( − ) ≈ . Thus, equation (26) if it is rescaled into the dimensional form can be written as: Then the linearization of this equation by searching solution of the form solution: = ( − ) , will lead to the dispersion relation between the frequency and the wavenumber takes the form:

• Finite Difference Scheme
The finite difference method is the oldest and the simplest numerical method. The modern version is the spectral method, which is not going to be considered here. The KdV equation has three different differential terms, the third derivative, and the quadratic nonlinear term, and the variation with time of the problem. Each term can be extracted using Taylor series expansion. Starting from the third derivative, we are interested in the central difference, so we can write: By adding (27), (28) and then (29), (30), we get: +2 + −2 = 2 + 4∆ 2 2 2 , +1 + −1 = 2 + ∆ 2 2 2 .
The same when subtracting (27), (28) and then (29), (30) we get: By manipulating and re-arranging the previous equations we get: Which finally leads to: The first order differentiation also is going to be treated using central difference which is: And finally forward Euler finite difference for the time is assumed as follows: Then the final equation is given as: • Implementation The grid in the azimuthal direction is going to be discretized depending on the number of points extracted from the grid, which is different for each case depending on the extracted pixels and their corresponding millimeters, where the grid is built on the outer periphery of the cylinder and has center at the paper clip number 5 as it is clear in Figure (2) that is why we have a domain of negative and positive values. Those points are relatively big like for the case of 10000ml we have = 749 points, for the case of 12000ml we have = 899 points, and for the case of 14000ml we have = 649 points, the starting point of the domain is the first value of each angle and it is called 1 , and the ending point is 1 , the initial guess of the problem is assumed a Gaussian fit, where it takes its values based on the extracted pixles from which we determine both the average and the standard deviation . The distance between each point of the pixel coordinates causes instability problem using Courant number defined in (35) below, unless small value is introduced to ensure equation (35) is satisfied, which is a specific ratio between the time and the distance represented by Courant number : Where is assumed to be the mean velocity or the stream velocity, and is the disturbance speed. It is well known that choosing the value of Courant number bigger than one will certainly lead to instability of the assumed solution, that is why it is preferable to take value between zero and one. In the experiment, the time difference was chosen so that we satisfy the following:

Experiment and Results.
Full details of the experimental apparatus and procedure are in (ALSHOUFI 2021). In order to track the different oscillations in the channel under precession conditions, a real flume was constructed for this purpose. The Figure (1 Table (2) shows some characteristic and geometrical features of the design. This includes the three adjustable control parameters of the system: (i) The first is the amount of water in the tank, which can be quantified either by its volume or, preferably, by the mean water depth = ℎ ̅ , by which we ensure the mean water level constraint, which is one-one relationship between volume of water and the average water level: (ii) The second is the tilt angle, characterized by its tangent as mentioned above, which is also the largest slope of the tilting table. (iii) The third is the angular velocity of the rotating table.
The energy is provided to the channel by using DC electrical motor that has external electrical supply, which is battery. The whole system is placed on the laboratory floor by a metal base, noting that currently the rotating table can be turned only in the counterclockwise sense. The inner cylinder has wall length higher than the outer one, the inner one was covered from inside with white papers A4 in size to reduce the light reflection. The motion inside the channel was tracked using JAI PULNiX TM-1405 GE camera, which has Charge-Coupled Device (CCD) sensor type. It has high picture resolution 1392 × 1040, with focal length 8mm, and pixel resolution 4.65 × 4.65 .
(a) (b) Figure (1). (a). The Real Flume Channel in the Laboratory. (b). Sketch of experimental setup. The camera system in front of the flume almost at the same level of the water level inside the tank. The pictures of this camera were extracted first, and then they were calibrated using a program was designed by the current author for this purpose, it depended on Zhang's method (2000), and the algorithm was applied is the one introduced by Burger (2016), with slight modification that included more radial distortion coefficients. The camera was set at fixed position in front of the channel where it could take continuous pictures, the frequency of capturing the pictures was 1 frame per 30 second. The outer cylinder was also covered with some white papers so that only the place in front of the camera was not covered for better vision. The outer illumination for the system was two bulbs fixed on the inner cylinder's wall to cover the whole inside area, and other external big lamp which was adjustable in its height was sometimes used and fixed in the laboratory floor next to the camera.
The waves that may form in such system are called the inertial forced oscillations, and the linear problem of their solutions leads to a set of modes that are called the Kelvin modes. Such oscillations appear as long as their frequencies smaller than twice the natural one, thus in this system, they all have the same frequency, but they may differ in their wavelength. The solitary wave in the Figure (2) is a solitary Kelvin wave, which accords with the fundamental mode = 1 in the axial direction of motion, = 1 radially, and azimuthal base mode 1 = 2. This result was extracted when analyzing the linear case of this problem, from which it was found that resonance in the system occurs when the wavelength matches the outer periphery of the channel. (cfr. ALSHOUFI 2021).
The observations of single Kelvin mode in the channel was first for the smallest tilt = 0.0117, for different volumes of water varied from 2000ml to 12000ml, the crest was little flattened, and had symmetric form, stable, and of course rotating without changing i its speed, approximately like the wave in the Figure (2), but it looked as if it is going to break, the crest slightly proceded the whole wave body. During the computation of those waves for relatively small volumes 2000ml, to 6000ml, it was noticed that the nonlinearity is much bigger than the dispersion. This leads to assume that such Kelvin solitary mode can be considered for bore development (EL 2007), and that work also is being carried out by the present author. Increasing the volume of water in the channel, leads to bigger amplitudes thus bigger wave lengths, like the wave in the figure (3) for instance. Those cases look really like solitary waves, symmetric in form, not changing their speeds as propagating ahead, but unfortunately they are not stable, they preserve this smooth symmetric form for only one round about the outer periphery of the channel, then disperse into waves have smaller amplitudes and flattened crests. In connections with the forced oscillations in such systems under precession, similar observations in closed tanks were first noticed by McEwan (1970). He called this phenomenon resonance collapse, where big waves retain their smooth form for relatively short time, then they collapse. Where the original form cannot be regained even if the flow waxes and wans. And in fact, the computations for the wave in Figure (3) showed that slight priority for dispersion at the expense of nonlinearity, still in good balance to be computed using KdV equation form assumed in (26), using its numerical estimation in (34). Taking time history for the wave in such case is presented in Figure (    The finite scheme model introduced in the previous section in equation (34) is applied to the observed waves in different volumes, and comparison has been carried out. The initial condition was a Gaussian initial fit from the experimental data:  As no exact analytical solution for the derived equation, an analytical solution for the problem was derived as proposed by Brauer (2000) with little change as in equation (38). It is kind of fit solution for the results, and it was compared with the experiments at each time, as it is clear in Figure (10).
To track the wave evolution with time we construct the solution on the outer periphery of the cylinder to show the far field, Figure (9) shows how the wave evolves with time, the equation is not stable thus the need for weight coefficients is harnessed to keep it stable with time, the nonlinear term is multiplied with 0.01, and the dispersion term was multiplied with 0.001, the number of points = 70, and the time difference 0.00176. Of course, the wave starts to disperse under rotation effect, and the system is not stable at all, but the procedure of wave evolution with time shows that the wave starts to decay to the right which is similar to the real decay depicted in Figure (4). However, we cannot assume according to Figure (4) that we still have the solitary wave in essence, but this is the sample of how this system works.   The solitary wave really can be formed in such systems if the amount of water is big enough, but this wave also can suffer from breaking. It is well known in natural cases, out of rotational effects, from the gas dynamic theory that waves which travel in one direction only, their crests travel faster than troughs, so that the front of each crest continuously steepens, until it breaks down. (Ursell 1953). In Figure ( 11) we got such collapse of the inertial Solitary Kelvin mode. It is clear how the wave crest precedes the whole body of the wave in the same direction of the system propagation, which is counterclockwise. Russell (1844) stated that the wave breaks when its amplitude equals the depth of water under it, which is the case in Figure (11), ℎ ̅ = 9.33cm, = 9.17cm. In the transverse direction on the other hand, it was noticed that the wave has exponential decay form away from the wall, without curving backward in that direction. The radial wavenumber that is used in the experiment is given by: where = the ratio between the inner and the outer radiuses as in table (2). Simple fit was used assuming exponential form as follows: where is an arbitrary constant that suits the fit, and it is different in each experiment, the wave amplitude, the radial distance between the inner and the outer cylinders, 1 is dimensionless Rossby Radius, which is given as: where the shallow water velocity: = √ ℎ ̅ , and the frequency of oscillation. This exponential fit for the radial amplitude decay is presented in Figure (  To estimate the rotation and viscous effects, a computation of Reynolds, Ekman, and Rossby numbers was carried out, for each Solitary Kelvin wave mentioned above, at different volumes of water. Reynolds number determines the type of the flow, like laminar and turbulent. In open channels it can be given in terms of its characteristic length, the hydraulic radius, in the system under study it is assumed as: where is the kinematic viscosity which is has a constant value assumed during the experiment 1.03 × 10 −6 , ℎ = / the hydraulic radius, the wetted perimeter: = + 2ℎ ̅ , as the crosssection for the channel has rectangular form, with = − the channel width and ℎ ̅ is the average depth, the rotation speed. It is known that Reynolds in open channels has different ratios in comparing the turbulent and laminar effects. If ≤ 500 the flow is laminar, and ≥ 1000 the flow is turbulent. In this system, the computations showed that for resonant conditions that agree with Kelvin Solitary wave, the flow is always turbulent. To measure the rotational effect Rossby number is computed, which is a dimensionless number that represents the ratio between inertial and Coriolis forces, it is assumed in this system as in (43): where 1 = precession rate. The computations of Rossby number showed that for all the above-mentioned cases of tilt is constant and equal to 0.492, which shows moderate rotation effect on wave motion. Ekman number represents the balance between the small viscous forces and the small Coriolis ones, and it plays an important role in determining the boundary layer effect on the solid walls. It is assumed to be the inverse of Reynolds number as it is clear in (44), to ensure that the inviscid effect is satisfied. The computations of this number showed that it is of order (10 −4 ): = ℎ 2 . (44) As we deal with open channel, computations for Froude number using (45) showed that the flow all the time supercritical regardless the volume of water in the channel: The velocity of the wave was measured by taking the difference in time between two different pictures, from which the difference between the wave crest in each one is approximately = 2 , which is the crossed distance per time. Amaouche et al. (2013) proposed that the velocity in their tank theoretically was given by: In comparison with this work, it was found high match with experiment as shown in table (3):

Conclusion.
In this paper, a new KdV model equation was derived for the case of the flow in an open cylindrical channel under precession conditions. The derivation starts from assuming the inviscid irrotational conditions in two-dimensions of the problem, as the radial part is short and barely affects the dominant azimuthal flow direction. The derivation first was to the forced Bernoulli equation, then it followed the classical perturbation method to derive the final version of this eqaution. The equation has a forcing term, that makes the equation is not solvable directly, and a simple Finite Difference Scheme is used to solve it. However, as noticed from the experiment that this term is small and can be neglected from the equation and kind of analytical solution was derived, accordingly. Both the numerical and the analytical one were compared with the observed wave in the experiment, with very good agreement. It is worth to mention that during the experiment an assumption was made for Rossby radius that takes the local frequency, not its double value. The radial variation of the waves observed in this system has exponential decay, and the wave crest was not curved backward. It was noticed that the effect of rotation in this case is severe where once the wave amplitude is bigger or equal to the average depth of water it either breakes or disperses, in waxing and wanning mode from which no recovery to its original form, this is called the resonance collapse for the case of forced oscillations. Other thing is worth to mention that the breaking effect when starting from rest conditions is similar to the normal rotation effect in gas dynamics (Ursell 1953) where the observed waves have crests proceed the whole body of the wave. Although the velocity of the wave is taken the solid-body rotation one, another one that was proposed by Amaouch et al.
(2013) based on the long wave one was used. It was in a perfect match with the wave velocity noticed in the experiment.