Influence of cross channel depth variation on ship wave patterns
|
|
- Chad McKenzie
- 6 years ago
- Views:
Transcription
1 Dept. of Math. Universit of Oslo Mechanics and Applied Mathematics No. ISSN 9 3 August Influence of cross channel depth variation on ship wave patterns T. Torsvik 3, G. Pedersen and K. Dsthe 1 1 Department of Mathematics, Universit of Bergen, Johannes Brunsg. 1, N-5 Bergen, Norwa Department of Mathematics, Universit of Oslo, PO box 153, 31 Oslo, Norwa 3 Bergen Center for Computational Science, Universit of Bergen, Høteknologisenteret, Thormøhlensgate 55, N-5 Bergen, Norwa August 5, Abstract The present work is a numerical investigation on waves generated b a pressure disturbance moving at constant speed in a channel with a variable cross-channel depth profile. The channel profile, which is uniform in the along-channel direction, has a deep trench located in the vicinit of the center line of the channel, and shallow banks near the channel walls. Wave fields generated upstream and downstream of the moving pressure disturbance are described, and the characteristic features present in the wave patterns are related to the parameters governing the speed of the pressure disturbance and the shape of the cross-channel profile. Our numerical investigation is based on the COULWAVE long wave model, which solves a set of Boussinesq tpe equations in two horizontal dimensions. 1 Introduction In coastal waters, ships are often required to navigate in natural or dredged trenches. Long waves are known to be generated b high speed vessels traveling at speeds close to the critical depth Froude number F = U gh = 1 (1) where U is the speed of the vessel, g is the acceleration of gravit and h is the depth. In order to examine how ship waves behave under these conditions, we consider the idealized case where waves, generated b a moving pressure disturbance, propagate in a channel with a deep trench at the centre line and shallow banks along the channel walls. We are particularl interested in cases where the speed of the pressure disturbance is subcritical with respect to the Froude number in the trench, but near critical or supercritical with respect to the shallow banks. Because the ship wave pattern changes substantiall in the transition between the subcritical and supercritical regime, these test ma provide us with a better understanding on how the different features in the bathmetr influences the wave pattern. Several researchers have studied waves propagating in channels with arbitrar cross-section profiles[19, 1] and the wave patterns, in two horizontal dimensions (HD), generated b a 1
2 disturbance moving at speeds close to the critical Froude number in channels with a rectangular cross-section profile[1, 3, 17]. Mathew and Aklas[1] brought these elements together in their stud of waves propagating in channels with a trapezoidal cross-section profile. Recentl, Teng and Wu[], Jiang et al.[], and Liu and Wu[] have made contributions to this field of research. Most of these studies deal with channels with trapezoidal cross-section profiles, but Jiang et al.[] also include results for a channel with a deep trench along the center line, which is similar to the cases that will be discussed in this paper. In our stud, we simulate the wave generation and propagation using COULWAVE, a computer model for long waves, developed at Cornell Universit b Lnett, Wu and Liu[, 1, 7]. The model is based on a set of weakl dispersive, full nonlinear Boussinesq equations, first developed b Liu[5] and Wei et al.[5]. A version of the model which includes a pressure disturbance has been implemented b Liu and Wu[]. Mathematical model equations The numerical model is based on a set of full nonlinear and weakl dispersive Boussinesq equations, which are nondimensionalized b introducing the characteristic depth h, as the length scale, g/h as the time scale, and the hdrostatic pressure ρgh, as the pressure scale. In dimensionless form, the Boussinesq equations consist of the continuit equation ζ t + [(ζ + h)u α] {( ) z + α h { + ζ ( h ( u α ) + z α + h [(z α 1 ζ ) ( (hu α )) + 1 ) } h [ (hu α )] (z α 13 ζ ) ( u α ) ]} () =, and the momentum equation u α t + (u α )u α + ζ + p { ( ( + [z α h u ))] α + t { 1 + z αu α ( u α ) + z α u α ( hu α ) + 1 ( [ (hu α)] ζ h u )} α t { + ζ( hu α )( u α ) 1 ( ) uα ζ t { 1 + ζ ( u α ) 1 } ζ u α ( u α ) =, [ ( 1 z α u )]} α t } ζu α ( hu α ) (3) where ζ(x,, t) is the surface displacement, u α (x,, t) = (u α (x,, t), v α (x,, t)) is the velocit at the depth z = z α, = ( / x, / ) is the horizontal gradient operator, ρ is the densit of water and p(x,, t) is the pressure disturbance of magnitude p a at the free surface. The reference
3 depth z α =.531h is applied, following the recommendation b Nwogu[1]. This optimizes the dispersive properties of the equations with respect to the linear dispersion relation over a range of long and intermediate waves. With the model equations formulated as above, wave propagation in terms of phase velocit and group velocit is accuratel described for waves with dimensionless wave number in the range < kh π, as shown b Madsen and Schäffer[11]. When deriving () and (3), it is assumed that the characteristic water depth h, is small relative to the horizontal length scale λ, i.e. µ = h /λ 1. The characteristic wave amplitude a ma be of the same order of magnitude as h, and the nonlinear effects, related to the parameter ǫ = a/h = O(1) are not required to be weak. The model equations () and (3) are accurate up to O(µ ). In our simulations we represent a monohull vessel b a localized pressure disturbance (see e.g. [1, 3, 17, ]). Waves generated b a thin ship representation was also examined in [17] and [], where the ship is represented b a slender bod extending from the free surface down to the sea bed. These studies showed that the different wave generating mechanisms generated qualitativel similar wave patterns, although the waves generated b a pressure disturbance were generall smoother than waves generated using the thin ship representation. In our stud we are not concerned about waves generated b a ship with a particular hull shape, and in this case the pressure disturbance representation is adequate. The forcing pressure disturbance is defined according to p(x + Ft, ) = p a f(x + Ft)q(), where the Froude number F is defined according to eq. (1). Functions f and q are defined b [ ] π(x + Ft) f(x + Ft) = cos, L x + Ft L, L ( q() = cos π ), R R, R on L x + Ft L, R R, and is zero outside this rectangle. Tpical values for the pressure disturbance is p a =.1, L =.h, and R =.h. In the following figures, the pressure disturbance moves towards decreasing values of x, and the along-channel coordinate is redefined as = x + Ft, in which the location of the pressure disturbance remains fixed at all times. In this frame of reference, upstream waves are located to the left of the pressure disturbance. At the start of the simulation, the pressure disturbance is abruptl set in motion at a particular speed, which is maintained constant thereafter. Numerical method The COULWAVE model applies an algorithm which is similar to the method of Wei et al.[5]. The algorithm is formall accurate to t in time. The first-order spatial derivatives in the equations are discretized to fourth order, whereas the spatial derivatives in the dispersive terms are discretized to second order. This ensures that the numerical dispersion related to the spatial derivatives of first order is of higher order than the phsical dispersive terms. The algorithm 3
4 is marched forward in time using the Adams-Bashforth-Moulton predictor-corrector scheme to update the dependent variables. Although the algorithm is formall of second order in space, this convergence rate has not been achieved in the simulations. The numerical results converge when the grid resolution is increased, but the convergence is slow for the steep, large amplitude waves often found immediatel downstream of the pressure disturbance. A five point filter has been applied to avoid saw-tooth noise which tends to emerge in the cross-channel direction. A small damping is also inherent in the time stepping scheme. No other dissipative terms have been added to smooth the variable fields. The model has been tested b running the same simulation with spatial grid resolutions of.1 h,. h, and. h. Using the finest resolution as reference, we calculated the relative errors E u and E d, for the amplitudes of the leading upstream wave and the first downstream wave, respectivel. For x = =. h, the errors E u =. 1 and E d = were found, while increasing the grid resolution to x = =. h, reduced the errors to E u = and E d = A spatial grid resolution of x = =. h has been used for all subsequent simulations in this paper. The pressure disturbance is often required to travel a length of order O(1 3 h ) for the wave pattern to full develop. In order to speed up the computation, the wave field is onl computed in a window including the parts of the channel immediatel upstream and downstream of the pressure disturbance. This window is shifted upstream when the free surface displacement grows to δ =.5 p a somewhere along a cross-channel line located 3 h downstream of the upstream boundar. A potential problem with this method is contamination from waves generated at the downstream boundar, which is of particular concern for simulations where the Froude number is subcritical. B running the same simulation with different sizes of computational windows, we found that b the end of the simulation, the error was O(1 ) at 5h downstream of the pressure disturbance, and O(1 3 ) near the downstream boundar of the smallest computational window, located.5 h downstream of the pressure disturbance. We also tested for the sensitivit to the shifting criteria b running the model with δ =.1 p a. The errors were found to be O(1 5 ) everwhere, except near the downstream boundar. 3 Results z h 1 h W R x L W 1 W W 3 Figure 1: Channel cross-sectional profile, and sketch of the pressure disturbance.
5 In this section we present results based on numerical simulations where waves are generated b a moving pressure disturbance in a channel with a variable cross-channel topograph. The channel has a deep section along the center line and shallow banks near the solid boundaries. The parameters for the cross-channel variation are defined according to Fig. 1. We will refer to the depth at the center line as h 1, and the depth at the channel wall as h. The total width of the channel is W, but because the wave pattern is smmetrical across the line =, we run the simulations onl on the half plane W. The widths of the different regions are given b the triplet W = (W 1, W, W 3 ), where W 1 is the half-width of the trench, W is the width of the sloping region, and W 3 is the width of the shallow bank adjacent to the channel wall. This cross-channel profile is similar to one of the examples discussed b Jiang et al.[]. A direct comparison between our results and the result presented in Jiang et al.[] has not been attempted, because Jiang et al. do not quantif the free surface displacement attained in their simulation. 3.1 Influence of the Froude number on the wave pattern In this part we examine how the wave pattern is influenced b changing the velocit of the advancing pressure disturbance. The channel is defined b W = (3,, 5), h 1 = 1. and h =., and the velocities used in the simulations corresponds to Froude numbers at the trench of =.7,.,.9 and 1., which give Froude numbers F =.937, 1.3, and 1.91, respectivel, at the shallow banks. Figure shows contour plots of the generated waves at t = 153, figure 3 shows the leading downstream waves in greater detail, and figure shows along-channel wave profiles at =, 5 and 1. For =.7 and F =.937, the upstream wave generation, shown if Figs. and, seems to be transient in nature, and depends on the startup of the simulation (see discussion in section 3.5). In the region immediatel upstream of the pressure disturbance, both the wave amplitude and the mean surface elevation tend to zero with time. A stead state is attained downstream of the pressure disturbance, shown in figure 3, where the wave pattern consists of small amplitude waves with wave lengths of λ 3. Since hk is within the range where the wave dispersion is accuratel described, we can expect the results to be reasonable even for these short waves. This is a case which clearl demonstrates the advantage of using Boussinesq models with improved dispersion properties instead of the standard form of the Boussinesq equations. The leading part of the upstream wave in figure resembles the Air function Ai(), which is the traveling wave solution for the linearized equations. These waves displa onl a slight cross-channel variation in amplitude and no cross-channel variation in phase speed despite the significant cross-channel variation in depth. It ma seem strange that an essentiall linear waves should not be influenced b the depth variation. In this case the leading waves are long not onl with respect to the depth, but also with respect to the width of the channel, and is therefore onl influenced b the mean depth in the channel. This is also consistent with the fact that the cross-channel variation increases as the wave length decreases in the wave train. When the Froude number is increased to =. and F = 1.3, an undular bore is generated upstream of the moving pressure disturbance, as seen in Figs. and. We note that the wave crests of the upstream waves span the channel at an angle which is nearl perpendicular to the channel wall, and that the amplitude, although higher at the wall than at the center of the channel, varies onl slightl in the cross-channel direction. Waves with large amplitudes are generated downstream of the pressure field for =.. These high amplitudes 5
6 1 =.7, F =.9, H = =., F = 1.33, H = =.9, F = 1.1, H = (c) Figure : Contour plot of waves generated at different Froude numbers. Contour levels given b H for each plot. occur near the channel wall ( = 1) due to wave reflection at the solid boundar, but also in the shallow part of the channel (5 < < ) due to wave-wave interaction. While the waves generated for =.7 and =. clearl conforms to subcritical wave patterns, the wave pattern generated at =.9 and F = 1.119, seen in Figs. (c) and (c), shows what appears to be a persistent generation of upstream waves which is characteristic for the near critical wave pattern. It is notable that this wave pattern occurs when the local Froude number at the trench is still subcritical. The leading upstream wave attains the bellshaped form often associated with solitar waves, but is trailed b small amplitude disturbances which interact with the following upstream waves. Diffraction of the leading upstream wave due to depth variation has also been reported b [1], who studied the propagation of an undular bore in a trapezoidal channel. The wave diffraction due to depth variation is more prominent in the middle of the channel than at the channel walls, as seen in Fig. (c). Even though wave diffraction occurs, causing significant variation in amplitude and wave profiles in the cross-
7 1 =.7, F =.9, H = =., F = 1.33, H = =.9, F = 1.1, H = (c) = 1., F = 1.91, H = (d) Figure 3: Details of the downstream wave pattern. Contour levels given b H for each plot. channel direction, each of the upstream waves can be clearl identified as distinct waves with wave crests oriented nearl perpendicular to the channel wall. The amplitudes of the downstream waves for =.9, seen in Fig. 3(c), are of the same order of magnitude as for =., and the largest peaks are located near the center line of the channel and near the channel walls. 7
8 .15.1 =.7, F =.9. =., F = =..1 = 5. = =. = 5. = =.9, F = = 1., F = =. = 5. = (c).1..3 =. = 5. = (d) Figure : Wave cuts along the channel of the wave patterns shown in Fig.. Waves do not propagate upstream of the pressure field for = 1. (F = 1.91), so the wave pattern consists onl of the downstream waves, as seen in Fig. 3(d). Although the reflection of the leading wave at the wall resembles the stem seen in Mach reflection (see e.g. [15]), this is not a feature of Mach reflection because the stem does not grow with time and the wave does not propagate upstream relative to the pressure disturbance. The lack of upstream waves for = 1. ma seem to contradict results from previous studies. We return to this point in the discussions at the end of the paper. Due to the lack of upstream wave generation, all the wave energ is transferred downstream of the pressure field. Hence the wave amplitude of downstream waves are seen to be significantl larger than for =.9, particularl where the waves amplif near the channel wall and the center line of the channel. 3. Waves radiated upstream of the pressure field For =., the upstream waves propagate as an undular bore. The position of the leading wave crest w (, t) is plotted in figure 5, for =, 5, and 1. Earl in the simulation, the position of the leading wave is located further upstream at the center of the channel than at the channel walls. While the wave crest at the center line initiall remains nearl stationar compared to the position of the pressure disturbance, with w (, t) = 3.3 ±.1, the wave crest at the channel wall is seen to propagate upstream. The straightening of the wave crest at the channel wall and the upstream propagation is clearl a feature of Mach reflection, which was previousl found b Pedersen[17] to occur for waves generated in channels with a rectangular cross-channel profile.
9 w Position of leading wave =. = 5. =1. w (,t) w (,t) Position relative to w () = 5. = time 1 3 time.5 Leading wave amplitude. Cross channel position and amplitude =. = 5. = time (c) w (,t) w (,t) ( ) (d) ( ) Figure 5: Position and amplitude of the leading upstream wave for =.. The wave propagation at the wall slows down at t 75, when the wave crest has become nearl a straight line, as seen in figure 5, where the position of the wave crest w (, t) is plotted relative to w (, t). This coincides with a temporar slump in the wave amplitude at the wall, as seen in figure 5(c). Once the wave amplitude at the wall increases, the waves crest over the shallow bank propagate a considerable distance upstream of the crest at the center line, before an abrupt adjustment occurs at t 135, which brings the wave crest at the center line into alignment with the wave crest near the channel wall. From figure 5 it is apparent that the adjustment shifts the central part of the upstream wave forward to the position of the wave near the channel wall. After the adjustment, the leading wave propagate with a constant speed c w of approximatel c w /U = 1.1, which is subcritical (c w / h 1 =.9) at = and supercritical (c w / h = 1.15) at = 1. The cross-channel amplitude increases from the center line to the channel wall, as shown b the stapled line in figure 5(d), and the mean wave amplitude is increasing slowl with time, as seen in figure 5(c). This ma however be a transient phenomenon, because the upstream waves form an undular bore which is not full developed within the time shown, which means that wave energ is still transfered to the leading wave from the rear. The wave crest is located slightl further upstream near the channel wall than the center line, as seen b the solid line in figure 5(d). The generation of the leading upstream wave for =.9 follow much the same pattern 9
10 1 5 Position of leading wave =. = 5. =1. Position relative to w () = 5. =1. w w (,t) w (,t) time time.3 Leading wave amplitude.5 Cross channel position and amplitude =. = 5. = time (c) w (,t) w (,t) ( ) (d) ( ) Figure : Position and amplitude of the leading upstream wave for =.9. as for =.. At the start of the simulation, the upstream crest at = is stationar at w (, t) = 3. ±.3, but adjusts to the wave crest in the far field at t, as seen in Fig.. The wave crest is located nearl the same distance upstream near the center line and near the channel wall, but lags behind b.h at =, as seen in figure (d). The wave amplitude decreases with time, as seen in figure (c), which is to be expected because the wave is subject to wave diffraction. The wave propagates at approximatel c w /U = 1.7, which is subcritical (c w / h 1 =.9) at = and supercritical (c w / h = 1.5) at = Comparison with cross-channel averaged theor For the channel configurations used in the examples above, we often find that the leading upstream wave is long compared to the width of the channel. These results ma therefor be comparable to results using cross-channel averaged models, which were first developed b Peters[19] and Peregrine[1], and have more recentl been applied b Teng and Wu[]. The models are based on weakl nonlinear and dispersive equations, such as KdV and standard Boussinesq equations, and are developed under the assumption that the cross-channel variation of the velocit and wave amplitude is small compared to the along-channel variation. Under these conditions, it is possible to determine the solitar wave solution analticall, given a mean wave amplitude and a specific cross-channel depth profile. The procedure is briefl outlined in 1
11 1 COULWAVE model: =.7, H =.1 x TW model: =.7, H =.1 x 1 3 ζ ζ 1 COULWAVE model: =., H =.5. 1 TW model: =., H = ζ ζ.1 (c) (d) 1 COULWAVE model: =.9, H =.5 1 TW model: =.9, H = ζ ζ.5 (e) (f) Figure 7: Comparison between TW and COULWAVE models. the appendix. Figure 7 shows a comparison between the leading upstream wave at ζ = w, found in the simulations with the COULWAVE model, and wave solutions from the Teng-Wu (TW) cross-channel averaged model. For =.7, it is appropriate to compare the simulated result in figure 7 with the Air function, which is the solution in the linearized TW model. In the simulated result, the leading wave has a steeper wave front and the trailing waves are smaller than for the TW model result based on the Air function. This ma be due to nonlinear effects, which are small but not negligible for these waves. There is however a good agreement between the models on the cross-channel amplitude variation, which deviates b approximatel ±1.% at the centre line and at the channel wall. The results for =. and =.9 are compared to the corresponding solitar wave solution of TW with the same mean amplitude. Results for =., shown in figures 7(c) and 7(d), are in reasonable agreement. The cross-channel 11
12 variation in amplitude is smaller than for the COULWAVE model, resulting in a deviation between the amplitudes of approximatel ±.5% at the center line and at the channel wall. The wave crest of the TW solitar wave is also broader than the wave crest in the COULWAVE model. For =.9, the TW solitar wave solution is not in agreement with the COULWAVE result. The leading wave in the COULWAVE simulation has a wave length that is comparable with the channel width, and is therefor influenced b the cross-channel depth variation. Because of the short wave length, it ma also be influenced b the waves following in the wave train. These results show that cross-channel averaged theor ma be applicable for small amplitude solitar waves propagating in channels with a variable cross-section topograph, but onl if the wave length is large compared to the channel width. Figure 7 shows the comparison between the solitar wave solution from the Teng-Wu (TW) cross-channel averaged model and the leading upstream wave at ζ = w, found in the simulations with the COULWAVE model. The best match is achieved for =., shown in figures 7(c) and 7(d). The cross-channel variation in amplitude is smaller than for the COULWAVE model, resulting in a deviation between the amplitudes of approximatel ±.5% at the center line and at the channel wall. The wave crest of the TW solitar wave is also broader than the wave crest in the COULWAVE model. For =.7, the amplitude deviation is approximatel ±.1% at the center line and channel wall, but now the wave length of the TW solitar wave is significantl longer than the leading wave in the COULWAVE model. This result is consistent with our previous assertion that the upstream wave behaves like an Air function, and does not develop independent solitar waves. For =.9, the TW solitar wave solution is not in agreement with the COULWAVE result. The leading wave in the COULWAVE simulation has a wave length that is comparable with the channel width, and is therefor influenced b the cross-channel depth variation. Because of the short wave length, it ma also be influenced b the waves following in the wave train. These results show that cross-channel averaged theor ma be applicable for small amplitude solitar waves propagating in channels with a variable cross-section topograph, but onl if the wave length is large compared to the channel width. 3.3 Amplification in the wake wave pattern Wave amplification due to wave-wave interaction can sometimes occur in the downstream wave pattern, as seen in figure 3. The strongest amplification is located at approximatel = 1 and = 7, where the first reflected wave in the downstream wave packet interacts with the second downstream wave. A similar feature is present at approximatel = 35. A time sequence of the downstream wave pattern is shown in figure, and the temporal development of the amplitude of the leading downstream wave at = 7 and the two leading waves at = 1 is shown in figure 9. At first, the wave amplification grows with increasing amplitude of the reflected leading wave, as seen in figure. The amplitude of the leading wave peaks at approximatel t = 3, attaining an amplitude of =.157 at the channel wall, and steadil declines after this time. The amplitude of the wave at = 7 continues to increase as the second downstream wave grows large. Due to the decreasing amplitude of the leading wave, the amplified wave at = 7 can not be sustained indefinitel, and the amplitude decreases significantl after t = 15. Amplification due to wave-wave interaction has been found in test cases with a slightl different depth profile (W = (,, )), and with different velocit of the pressure disturbance ( =.9, W = (,, )), as seen in figure 1. It is not clear at this point whether these amplified waves will alwas be transient, or if the can occur in the wave pattern after a stead state 1
13 1 =., F =1.33, H=.5, T= =., F =1.33, H=.5, T= =., F =1.33, H=.5, T= =., F =1.33, H=.5, T= (c) (d) Figure : Wave amplification on the shallow bank due to wave-wave interaction has been attained. Wave amplification due to wave-wave interaction has also been observed in channels with a rectangular cross-section profile, where the waves were generated b two pressure patches of equal size in a configuration resembling the hull of a catamaran, indicating that a crosschannel depth variation is not a necessar condition for this phenomenon. A similar phenomenon has been discussed b Peterson et al.[], who studied extreme waves occurring in the intersection of two solitar wave groups. Peterson et al.[] suggested that such a phenomenon could occur as a result of wave-wave interaction between waves from different high speed vessels, or waves from a single high speed vessel changing its direction of travel. As mentioned earlier, the model equations reproduce the linear dispersive properties of waves in the range < kh π. However, [11] showed that the accurac of these equations with respect to nonlinear properties degraded sooner than the linear properties. It could therefore be questioned whether the fairl short, steep waves in the downstream wave pattern was well represented in the model. [9] explored the possibilit of using vertical laers to represent the flow, and found that the resulting set of equations faithfull reproduced the linear and secondorder nonlinear behavior up to kh. Fig. 11 shows a comparison between results for the 1-laer and -laer versions of the COULWAVE model corresponding to the situation in Fig.. There is a slight phase difference between the results, which becomes more apparent far downstream of the pressure patch, and the difference in maximum amplitude is approximatel %, but the downstream wave amplification is clearl maintained also in the -laer model. 13
14 Figure 9: Time development of wave amplification on the shallow bank. 1 W=(,,), h =., =., F =1.33, H=.5, T= W=(,,), h =., =.9, F =1.1, H=.5, T= Figure 1: Wave amplification due to wave-wave interaction for different Froude numbers and channel profiles. Time of output at t = 9, with contour levels of H = Results for different channel parameters In this section we consider how changes in the channel width parameters W = (W 1, W, W 3 ), and channel depth parameter h (keeping h 1 = 1), influences the wave pattern. In all these simulations we have maintained the Froude number at the trench constant at =.. Figure 1 show wave cuts along the channel for five test cases. Figure 1 corresponds to the case described earlier, with W = (3,, 5), and h =.. In figures 1 and 1(d) the depth h is changed to.5 and., respectivel, while maintaining the same values for W 1, W and W 3 as in figure 1. In figure 1(c) the widths are W = (1,, 7), and in figure 1(e) the widths are 1
15 1 ψ Figure 11: Comparison between simulations using 1 and vertical laers: 1 laer laer. Results corresponds to Fig. for t = 7. W = (,, ), and the depths h 1 and h are the same as in figure 1. The variation in the wave pattern displaed in these results clearl depend both on the width mean Froude number F = 1 W W F()d, and on the channel configuration. Table 1 shows the width mean Froude number and the width mean wave amplitude of the leading upstream wave for the five plots in figure 1. Although the Constant W = (3,, 5) Constant h =. h F max W F max (,, ) (1,, 7) 1(c) (,, ) (3,, 5) (,, ) (5,, 3) (,, ) 1(e) (d).5.17 (7,, 1) (,, ) (1,, )..97 Table 1: Table of channel parameters used in Fig. 13, the corresponding width-mean Froude number, and the amplitude of the leading upstream wave. When wave profiles are shown for a particular set of channel parameters, superscripts indicate the figure number. increase in wave amplitude near critical F resemble the results from similar tests for channels with rectangular cross-channel profiles, the magnitude of the wave amplitude and the wave period is influenced b the depth variation, as seen in figures 1 and 1(c). The influence of the channel configuration is seen more clearl in figure 13, where the amplitude of the leading upstream wave is plotted as a function of the width mean Froude number. These results suggest that at a given value of F, the largest upstream waves will be generated in the channel with the most narrow trench. 15
16 =., F = =. = 5. = W=(3,,5), h =.5, =., F =1.131, T = 153 W=(1,,7), h =., =., F =1.33, T = =. = 5. = =. = 5. = (c) W=(3,,5), h =., =., F =.9, T = 153 W=(,,), h =., =., F =1.33, T = =. = 5. = (d).1. =. = 5. = (e) Figure 1: Waves generated for various channel width and depth parameters given in Tab. 1. Wave cuts along the channel at t = Effect of an initial acceleration Jiang et al. discuss the sensitivit of the results to the starting condition, and conclude that different initial accelerations ield different results. In our simulations the pressure disturbance accelerates from rest to the designated speed over a single time step. In figure 1 we compare results for the abruptl started pressure propagation with results where the pressure disturbance accelerates from rest over a time T. Results are shown for a later time in the acceleration case than in the abrupt start case, to compensate for the dela in the generation of upstream waves in the former case. Differences between the results are clearl visible for =.7, where the upstream wave generation is transient in nature. A slight deviation is also visible in the upstream waves for 1
17 .1.1 depth variation width variation Froude number Figure 13: Amplitude of the leading upstream wave for the various for various channel width and depth parameters given in Tab. 1. Open circles signif Froude number variation due to variation of h. Dots signif Froude number variation due to variation of the width parameter W. =., but the deviations clearl diminish with increasing. We therefore conclude that the deviations due to the starting condition mainl influences the transient wave phenomena, and to lesser extent the wave pattern attained after long times. Concluding remarks In the preceding section we presented results for wave generation and propagation in a channel with a variable cross-channel topograph. A case stud for a similar channel configuration have previousl been presented b [] for = 1.. In figure 3(c) in [] we clearl see a wave propagating upstream of the wave-generating disturbance. Our result for = 1. and p a =.1, seen in fig. 3(d), does not indicate upstream wave propagation. We do however find waves propagating upstream of the pressure disturbance if we run the simulation with p a =.15, as seen in Fig. 15. Wave generation obviousl depends on both the Froude number and the magnitude of the pressure patch. In channels with a uniform cross-section topograph, the latter effect is often represented in terms of the blockage coefficient; the ratio of the maximum cross-section ship area to the cross-section channel area (see e.g. [3]). In wide channels where the wave pattern displas significant cross-channel variation, Mach reflection is seen to be crucial for the generation of straight crested upstream solitar waves. Mach reflection was analzed b [13, 1], who found that a Mach stem could be formed if ψ i / 3a i 1, where a i is the amplitude of the incident wave and ψ i is the angle of incidence. [1] showed in a numerical stud that conditions for Mach reflection was more restricted for large amplitude waves, i.e. the tended to reflect normall, but also showed b example that the limit of Mach reflection Mach reflection for a i =.3 was ψ i = 37.. Table shows a comparison between the result of [1], and results for p a =.1 (Fig. 3(d)) and p a =.15 (Fig. 15) based on measurements from the figures. These results show wh Mach reflection, and subsequent generation of upstream solitar waves, occur for p a =.15 and not for p a =.1 for the = 1. case. Although the p a =.15 result (Fig. 15) corresponds well with the result of [], this particular configuration ma not be phsicall sound. The Boussinesq equations tend to allow waves of 17
18 .15 =.7, t 1 = 7, t = 533, =..15 =.7, t 1 = 7, t = 533, = T = T = T = T = =., t 1 = 7, t = 7, =..15 =., t 1 = 7, t = 7, = T = T = 15. (c).5.1 T = T = 15. (d).3 =.9, t 1 = 5, t = 7, =..3 =.9, t 1 = 5, t = 7, = T = T = 15.. (e).1 T = T = 15.. (f) Figure 1: Comparison between results with and without acceleration from rest at the start of the simulation. 1
19 1 =1., F =1.91, H=.1, T= W=(3,,5), h =., =1., F =1.91, T = ψ=. ψ= 5. ψ= Figure 15: Wave pattern generated for p a =.15, at Froude numbers = 1., F = Time of output at t = 153, with contour levels of H =.1. a i ψ i ψ i / 3a i Tanaka p a = p a = Table : Parameters for Mach reflection for = 1., comparing the limit value found b [1] for a i =.3 with results for p a =.1 and p a =
20 unreasonabl large amplitudes and wave speeds, when in realit instabilities will cause these waves to undergo wave breaking. The upstream wave in Fig. 15 propagates at a speed of c w / h = 1.3 relative to the local shallow water wave speed near the channel wall, which exceeds the largest possible wave speed, c w /c = 1.9, for the solitar wave solution (see e.g. [15]). It is therefore doubtful if the simulated result, with respect to the maximum wave height, corresponds to a phsicall reasonable solution. In this paper we have studied the wave generation and propagation in a channel with a variable cross-channel topograph. We have found that waves ma propagate upstream of the wave generating disturbance, and that the crest of these waves will span across the channel in a nearl straight line despite the variation in depth across the channel. The wave amplitude and generation time varies with the width-averaged Froude number F, but is not governed b this parameter alone, as different configurations of widths and depths ma result in different wave patterns for the same value of F. We have also found that wave amplification ma occur due to wave-wave interaction in the downstream wave packet,and can be maintained over long times. Although the variation in cross-channel topograph seems to facilitate wave amplification, this does not seem to be a necessar condition for the occurrence of this phenomenon. Acknowledgements This work has received support from the Research Council of Norwa through a grant to the program Modelling of currents and waves for sea structures. We are indebted to Patrick Lnett for providing us with the latest version of the COULWAVE model, and to Philip L.-F. Liu and Tso-Ren Wu for guidance on how to use the model. Appendix In the TW model, the along channel wave propagation is computed b using a width averaged long wave model, such as the channel-boussinesq model or the channel-kdv model ( 1 ζ c t ζ h b ζ t + [ b( h + x ζ)ū ] =, ū t + ū ū x + ζ x 1 3 κ h 3 ū x t = p a x, ) ( ζ x ζ κ h x h ) h x + 1 b b x = 1 p a x, where c = (g h) 1 is a characteristic speed, and the shape factor κ is determined b the channel geometr (see e.g. Teng and Wu[, 3, ]). Width averages and cross section averages are computed b ( ) = 1 b b b ( )d, ( ) = 1 A A ( )d dz,
21 where b is the channel width at z = and A is the equilibrium area of the cross section. The general formulation allows b to var in the along-channel direction, but we have onl used a constant value for b in our applications. The cross channel perturbation is found b solving the Poisson equation with boundar conditions Ψ + Ψ zz = 1, Ψ z z= = h, (at the frees surface), Ψ n =, (at channel walls below the unperturbed water surface), where n is the unit normal vector at the solid boundaries. This determines the shape factor κ = 3 h( Ψ Ψ). The width averaged solitar wave solution with amplitude α is 3α ζ = α sech κ(x ct), and the width averaged solution to the linearized problem can be written in terms of the Air function Ai as ( ) 1 [ ) ζ = β 1 κ h 3 1] t Ai x( 1 κ h 3 t, where β determines the wave amplitude. The wave profile including cross channel effects is ζ(x,, t) = ζ(x, t) (Ψ z= Ψ) ζ xx The problem outlined above has been solved using the Poisson solver included in MATLAB. 1
22 References [1] R. C. Ertekin, W. C. Webster, and J. V. Wehausen. Waves caused b a moving disturbance in a shallow channel of finite width. J. Fluid Mech., 19:75 9, 19. [] T. Jiang, R. Henn, and S. D. Sharma. Wash waves generated b ships moving on fairwas of varing topograph. In Proceedings for the th Smposium on Naval Hdrodnamics, Fukuoka,. [3] C. Katsis and T. R. Aklas. On the exitation of long nonlinear water waves b a moving pressure distribution. part. three-dimensional effects. J. Fluid Mech., 177:9 5, 197. [] Y. Li and P. D. Sclavounos. Three-dimensional nonlinear solitar waves in shallow water generated b an advancing disturbance. J. Fluid Mech., 7:33 1,. [5] P. L.-F. Liu. Model equations for wave propagation from deep to shallow water. Adv. Coastal Engng., 1:15 157, 199. [] P. L.-F. Liu and T.-R. Wu. Waves generated b moving pressure disturbances in rectangular and trapezoidal channels. Journal of Hdraulic Research, ():13 171,. [7] P. Lnett and P. L.-F. Liu. Coulwave model page: Modeling wave generation, evolution, and interaction with depth-integrated, dispersive wave equations. [] P. Lnett and P. L.-F. Liu. A numerical stud of submarine landslide generated waves and runup. Roal Societ of London A, 5():5 91,. [9] P. Lnett and P. L.-F. Liu. A two-laer approach to wave modelling. Proc. R. Soc. Lond. A, (9):37 9,. [1] P. Lnett, T.-R. Wu, and P. L.-F. Liu. Modeling wave runup with depth-integrated equations. Coast. Engrg., ():9 17,. [11] P. A. Madsen and H. A. Schäffer. Higher-order boussinesq-tpe equations for surface gravit waves: Derivation and analsis. Phil. Trans. R. Soc. Lond. A, 35:313 31, 199. [1] J. Mathew and T. R. Aklas. On three-dimensional long water waves in a channel with sloping sidewalls. J. Fluid Mech., 15:9 37, 199. [13] J. W. Miles. Obliquel interacting solitar waves. J. Fluid Mech., 79(1):157 19, [1] J. W. Miles. Resonantl interacting solitar waves. J. Fluid Mech., 79(1): , [15] J. W. Miles. Solitar waves. Ann. Rev. Fluid Mech., 1:11 3, 19. [1] O. Nwogu. Alternative form of boussinesq equations for nearshore wave propagation. Journal of waterwa, port, coastal, and ocean engineering, 119():1 3, [17] G. Pedersen. Three-dimensional wave patterns generated b moving disturbances at transcritical speeds. J. Fluid Mech., 19:39 3, 19.
23 [1] D. H. Peregrine. Long waves in a uniform channel of arbitrar cross-section. J. Fluid Mech., 3():353 35, 19. [19] A. S. Peters. Rotational and irrotational solitar waves in a channel with arbitrar crosssection. Commun. pure appl. Math, 19:5 71, 19. [] P Peterson, T Soomere, J Engelbrecht, and E. van Groesen. Soliton interaction as a possible model for extreme waves in shallow water. Nonlinear Processes in Geophsics, 1:53 51, 3. [1] M. Tanaka. Mach reflection of a large-amplitude solitar wave. J. Fluid Mech., :37 1, [] M. H. Teng and T. Y. Wu. Nonlinear water waves in channels of arbitrear shape. J. Fluid Mech., :11 33, 199. [3] M. H. Teng and T. Y. Wu. Evolution of long water waves in variable channels. J. Fluid Mech., :33 317, 199. [] M. H. Teng and T. Y. Wu. Effects of channel cross-sectional geometr on long wave generation and propagation. Phs. Fluids, 9(11): , [5] G. Wei, J. T. Kirb, S. T. Grilli, and R. Subramana. A full nonlinear boussinesq model for surface waves. part 1. highl nonlinear unstead waves. Journal of Fluid Mechanics, 9:71 9,
Experiment 7 Energy Loss in a Hydraulic Jump
Experiment 7 Energ Loss in a Hdraulic Jump n Purpose: The purpose of this experiment is to examine the transition from supercritical (rapid) flow to subcritical (slow) flow in an open channel and to analze
More informationLecture 12: Transcritical flow over an obstacle
Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem
More informationHydraulic jump (lab. scale)
Hdraulic jump (lab. scale) http://www.outube.com/watch?v=c x6-_ejhdxy Hdraulic jump (field scale) http://www.outube.com/watch?v=h 45FauagdXw Downstream elevation (or Tail water) conditions are controlling
More information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationBOUSSINESQ-TYPE EQUATIONS WITH VARIABLE COEFFICIENTS FOR NARROW-BANDED WAVE PROPAGATION FROM ARBITRARY DEPTHS TO SHALLOW WATERS
BOUSSINESQ-TYPE EQUATIONS WITH VARIABLE COEFFICIENTS FOR NARROW-BANDED WAVE PROPAGATION FROM ARBITRARY DEPTHS TO SHALLOW WATERS Gonzalo Simarro 1, Alvaro Galan, Alejandro Orfila 3 A fully nonlinear Boussinessq-type
More informationJoule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of
More informationLocal energy losses at positive and negative steps in subcritical open channel flows
Local energ losses at positive and negative steps in subcritical open channel flows Nura Denli Toka and A Burcu Altan-Sakara* Middle East Technical niversit, Department of Civil Engineering, 63, Ankara,
More informationA weakly dispersive edge wave model
Coastal Engineering 38 1999 47 5 wwwelseviercomrlocatercoastaleng Technical Note A weakl dispersive edge wave model A Sheremet ), RT Guza Center for Coastal Studies, Scripps Institution of Oceanograph,
More informationViolent Sloshing H.Bredmose M.Brocchini D.H.Peregrine L.Thais
Violent Sloshing H.Bredmose Technical University of Denmark, M.Brocchini University of Genoa, D.H.Peregrine School of Mathematics, University of Bristol, d.h.peregrine@bris.ac.uk L.Thais Université des
More informationStrongly nonlinear long gravity waves in uniform shear flows
Strongly nonlinear long gravity waves in uniform shear flows Wooyoung Choi Department of Naval Architecture and Marine Engineering, University of Michigan, Ann Arbor, Michigan 48109, USA Received 14 January
More informationConsideration of Shock Waves in Airbag Deployment Simulations
Consideration of Shock Waves in Airbag Deploment Simulations Doris Rieger BMW Group ABSTRACT When the inflation process of a simple flat airbag was simulated with the MADYMO gas flow module, the resulting
More informationImproved Performance in Boussinesq-type Equations
Improved Performance in Boussinesq-type Equations Andrew B. Kennedy, James T. Kirby 1 & Mauricio F. Gobbi 2 Abstract In this paper, simple but effective techniques are used to improve the performance of
More informationOscillatory heating in a microchannel at arbitrary oscillation frequency in the whole range of the Knudsen number
Journal of Phsics: Conference Series Oscillator heating in a microchannel at arbitrar oscillation freuenc in the whole range of the Knudsen number To cite this article: O Buchina and D Valougeorgis J.
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationShoaling of Solitary Waves
Shoaling of Solitary Waves by Harry Yeh & Jeffrey Knowles School of Civil & Construction Engineering Oregon State University Water Waves, ICERM, Brown U., April 2017 Motivation The 2011 Heisei Tsunami
More informationTranscritical flow over a hole
Transcritical flow over a hole R. H. J. Grimshaw 1, D.-H. Zhang and K. W. Chow 1 Department of Mathematical Sciences, Loughborough University, UK Department of Mechanical Engineering, University of Hong
More informationHydraulic adjustment to an obstacle in a rotating channel
J. Fluid Mech. (), vol. 44, pp. 7 49. Printed in the United Kingdom c Cambridge Universit Press 7 Hdraulic adjustment to an obstacle in a rotating channel B L. J. PRATT, K. R. HELFRICH AND E. P. CHASSIGNET
More informationThe effect of disturbances on the flows under a sluice gate and past an inclined plate
J. Fluid Mech. (7), vol. 576, pp. 475 49. c 7 Cambridge University Press doi:.7/s7486 Printed in the United Kingdom 475 The effect of disturbances on the flows under a sluice gate and past an inclined
More informationInhomogeneity and the Post-Inflationary Universe
Inhomogeneit and the Post-Inflationar Universe Richard Easther and Matthew Parr Department of Phsics, Brown Universit, Providence, RI9 Abstract. We discuss the interaction between perturbations in the
More informationL O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S
L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S Yoël FORTERRE, Olivier POULIQUEN Institut Universitaire des Sstèmes Thermiques et Industriels (UMR 6595 CNRS) 5, rue Enrico Fermi 13453
More informationSaturation of inertial instability in rotating planar shear flows
J. Fluid Mech. (27), vol. 583, pp. 413 422. c 27 Cambridge Universit Press doi:1.117/s2211276593 Printed in the United Kingdom 413 Saturation of inertial instabilit in rotating planar shear flows R. C.
More information2017 年環境流體力學短期講座 Short Course on Environmental Flows
2017 年環境流體力學短期講座 Short Course on Environmental Flows 數學 海浪 與沿海動態過程 Mathematics, ocean waves and coastal dynamic processes Philip L-F. Liu National University of Singapore Cornell University September 2017
More informationThe effect of a background shear current on large amplitude internal solitary waves
The effect of a background shear current on large amplitude internal solitary waves Wooyoung Choi Dept. of Mathematical Sciences New Jersey Institute of Technology CAMS Report 0506-4, Fall 005/Spring 006
More informationChapter 4 Transport of Pollutants
4- Introduction Phs. 645: Environmental Phsics Phsics Department Yarmouk Universit hapter 4 Transport of Pollutants - e cannot avoid the production of pollutants. hat can we do? - Transform pollutants
More informationFrom bell-shaped solitary wave to W/M-shaped solitary wave solutions in an integrable nonlinear wave equation
PRAMANA c Indian Academ of Sciences Vol. 74, No. journal of Januar 00 phsics pp. 9 6 From bell-shaped solitar wave to W/M-shaped solitar wave solutions in an integrable nonlinear wave equation AIYONG CHEN,,,
More informationPan Pearl River Delta Physics Olympiad 2005
1 Jan. 29, 25 Morning Session (9 am 12 pm) Q1 (5 Two identical worms of length L are ling on a smooth and horizontal surface. The mass of the worms is evenl distributed along their bod length. The starting
More informationLévy stable distribution and [0,2] power law dependence of. acoustic absorption on frequency
Lév stable distribution and [,] power law dependence of acoustic absorption on frequenc W. Chen Institute of Applied Phsics and Computational Mathematics, P.O. Box 89, Division Box 6, Beijing 88, China
More informationCONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 742-761, 21. c Association for Scientific Research CONSERVATION LAWS AND CONSERVED QUANTITIES FOR LAMINAR RADIAL JETS WITH SWIRL R. Naz 1,
More informationA two-layer approach to wave modelling
1.198/rspa.4.135 A two-layer approach to wave modelling By Patrick Lynett and Philip L.-F. Liu School of Civil and Environmental Engineering, Cornell University, Ithaca, NY 14853, USA (plynett@civil.tamu.edu)
More informationWell-balanced shock-capturing hybrid finite volume-finite difference schemes for Boussinesq-type models
NUMAN 2010 Well-balanced shock-capturing hybrid finite volume-finite difference schemes for Boussinesq-type models Maria Kazolea 1 Argiris I. Delis 2 1 Environmental Engineering Department, TUC, Greece
More informationWater Flow in Open Channels
The Islamic Universit of Gaza Facult of Engineering Civil Engineering Department Hdraulics - ECIV 33 Chapter 6 Water Flow in Open Channels Introduction An open channel is a duct in which the liquid flows
More information* Ho h h (3) D where H o is the water depth of undisturbed flow, D is the thickness of the bridge deck, and h is the distance from the channel floor t
The Seventh International Colloquium on Bluff Body Aerodynamics and Applications (BBAA7) Shanghai, China; September -6, 01 Numerical simulation of hydrodynamic loading on submerged rectangular bridge decks
More informationExperimental study of a round jet in cross-flow by means of PIV
Experimental stud of a round jet in cross-flow b means of PIV Gennaro Cardone, Francesco G. Nese and T. Astarita Dipartimento di Energetica TErmofluidodinamica applicata e Condizionamenti Ambientali DETEC,
More informationChapter 4: Non uniform flow in open channels
Chapter 4: Non uniform flow in open channels Learning outcomes By the end of this lesson, students should be able to: Relate the concept of specific energy and momentum equations in the effect of change
More informationCombined tidal and wind driven flows and residual currents
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Combined tidal and wind driven flows and residual currents 1 Lars Erik Holmedal, Hong Wang Abstract : The effect of a residual current on the combined
More informationFast waves during transient flow in an asymmetric channel
Int. Conference on Boundar and Interior Laers BAIL 6 G. Lube, G. Rapin (Eds) c Universit of Göttingen, German, 6. Introduction Fast waves during transient flow in an asmmetric channel Dick Kachuma & Ian
More informationPhysics of Tsunamis. This is an article from my home page: Ole Witt-Hansen 2011 (2016)
Phsics of Tsunamis This is an article from m home page: www.olewitthansen.dk Ole Witt-Hansen 11 (16) Contents 1. Constructing a differential equation for the Tsunami wave...1. The velocit of propagation
More informationNumerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder
Bulletin of Environment, Pharmacolog and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue 1] 215: 318-323 214 Academ for Environment and Life Sciences, India Online ISSN 2277-188 Journal
More informationSurvey of Wave Types and Characteristics
Seminar: Vibrations and Structure-Borne Sound in Civil Engineering Theor and Applications Surve of Wave Tpes and Characteristics Xiuu Gao April 1 st, 2006 Abstract Mechanical waves are waves which propagate
More informationInternational Journal of Scientific Research and Reviews
Research article Available online www.ijsrr.org ISSN: 2279 543 International Journal of Scientific Research and Reviews Soret effect on Magneto Hdro Dnamic convective immiscible Fluid flow in a Horizontal
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationChapter 2 Basic Conservation Equations for Laminar Convection
Chapter Basic Conservation Equations for Laminar Convection Abstract In this chapter, the basic conservation equations related to laminar fluid flow conservation equations are introduced. On this basis,
More informationSubmarine sand ripples formation in a viscous fluid: 2D and 3D linear stability analysis
Marine Sandwave and River Dune Dnamics 1 & April 4 - Enschede, the Netherlands Submarine sand ripples formation in a viscous fluid: D and 3D linear stabilit analsis V. Langlois (1) and A. Valance (1) Groupe
More informationPrimary Results on Simulation of a New Kind Traveling Wave. of Permanent Form
Primary Results on Simulation of a New Kind Traveling Wave of Permanent Form Min CHEN y Abstract This talk concerns the nonlinear dispersive waves in a water channel. Results from the study of a new version
More informationNOTES ON OPEN CHANNEL FLOW
NOTES ON OPEN CHANNEL FLOW Prof. Marco Pilotti DICATAM, niversità degli Studi di Brescia Profili di moto permanente in un canale e in una serie di due canali - Boudine, 86 OPEN CHANNEL FLOW: uniform motion
More informationChapter 16 Mechanical Waves
Chapter 6 Mechanical Waves A wave is a disturbance that travels, or propagates, without the transport of matter. Examples: sound/ultrasonic wave, EM waves, and earthquake wave. Mechanical waves, such as
More informationSHADOW AND BACKGROUND ORIENTED SCHLIEREN INVESTIGATION OF SHOCK WAVES IN GAS-DISCHARGE MEDIUM
SHADOW AND BACKGROUND ORIENTED SCHLIEREN INVESTIGATION OF SHOCK WAVES IN GAS-DISCHARGE MEDIUM J. JIN 1, I.V. MURSENKOVA 1,c, N. N. SYSOEV 1 and I.A. ZNAMENSKAYA 1 1 Department of Phsics, Lomonosov Moscow
More informationThe Attractors of the Rayleigh-Bénard Flowof ararefied Gas
The Attractors of the Raleigh-Bénard Flowof ararefied Gas S.Stefanov, V. Roussinov and C. Cercignani Inst. of Mech., Bulg. Acad. of Sciences, Sofia, Bulgaria Dipart. di Matematica, Polit. di Milano, Ital
More informationMooring Model for Barge Tows in Lock Chamber
Mooring Model for Barge Tows in Lock Chamber by Richard L. Stockstill BACKGROUND: Extensive research has been conducted in the area of modeling mooring systems in sea environments where the forcing function
More informationSecond-order wave interaction with a vertical plate
Second-order wave interaction with a vertical plate B. Molin, F. Rem, O. Kimmoun École Généraliste d Ingénieurs de Marseille, 3 4 Marseille cedex 2, France Abstract The second-order interaction of monochromatic
More informationInstabilities of coupled density fronts
and their nonlinear evolution in the two-laer rotating shallow water model. Influence of the lower laer and of the topograph. B. Ribstein and V. Zeitlin Laborator of Dnamical Meteorolog, Univ. P.et M.
More informationarxiv: v1 [nlin.ps] 14 Aug 2012
Nonlinear shallow ocean wave soliton interactions on flat beaches Mark J. Ablowitz and Douglas E. Baldwin Department of Applied Mathematics, Universit of Colorado, Boulder, Colorado, 80309-056, USA (Dated:
More informationInternational Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India
International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow
More informationy R T However, the calculations are easier, when carried out using the polar set of co-ordinates ϕ,r. The relations between the co-ordinates are:
Curved beams. Introduction Curved beams also called arches were invented about ears ago. he purpose was to form such a structure that would transfer loads, mainl the dead weight, to the ground b the elements
More informationOn the Effect of the Boundary Conditions and the Collision Model on Rarefied Gas Flows
On the Effect of the Boundar Conditions and the Collision Model on Rarefied Gas Flows M. Vargas a, S. Stefanov a and D. Valougeorgis b a Institute of Mechanics, Bulgarian Academ of Sciences, Acad. G. Bonchev
More informationADVANCES IN MECHAN ICS Aug125, ,, KdV. , KdV. KdV. (m KdV). (non - propagating soliton or standing soliton), (precursor soliton)
8 Vol18 No1 1998 8 5 ADVANCES IN MECHAN ICS Aug15, 1998 LNM,, 100080 80,,, 150, Russell,,,,,,, KdV,,, :, KdV,,, KdV (m KdV) Schrodinger (NL S ) KdV, KP, KP 80, KdV :, (non - propagating soliton or standing
More informationDNS CODE VALIDATION FOR COMPRESSIBLE SHEAR FLOWS USING LST
Copright c 2016 ABCM DNS CODE VALIDATION FOR COMPRESSIBLE SHEAR FLOWS USING LST Jônatas Ferreira Lacerda, jonatasflacerda@hotmail.com Leandro Franco de Souza, lefraso@icmc.usp.br Josuel Kruppa Rogenski,
More informationVISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL CHANNEL THROUGH A POROUS MEDIUM
Volume 2, No. 1, Januar 214 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info VISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL
More information5. Two-layer Flows in Rotating Channels.
5. Two-layer Flows in Rotating Channels. The exchange flow between a marginal sea or estuary and the open ocean is often approximated using two-layer stratification. Two-layer models are most valid when
More informationTIME EVOLUTION OF SCOUR AND DEPOSITION AROUND A CYLINDRICAL PIER IN STEADY FLOW
A-7 Fourth International Conference on Scour and Erosion 28 TIME EVOLUTION OF SCOUR AND DEPOSITION AROUND A CYLINDRICAL PIER IN STEADY FLOW Shina UMEDA 1,Tatsuo YAMAZAKI 2, Hajime ISHIDA 1 Member of JSCE,
More informationMagnetohydrodynamic Fluid Flow Between Two Parallel Infinite Plates Subjected To An Inclined Magnetic Field Under Pressure Gradient
Journal of Multidisciplinar Engineering Science and Technolog (JMEST) Magnetohdrodnamic Fluid Flow Between Two Parallel Infinite Plates Subjected To An Inclined Magnetic Field Under Pressure Gradient Agnes
More informationFluid Mechanics II. Newton s second law applied to a control volume
Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.
More informationGENERAL SOLUTIONS FOR THE INITIAL RUN-UP OF A BREAKING TSUNAMI FRONT
International Symposium Disaster Reduction on Coasts Scientific-Sustainable-Holistic-Accessible 14 16 November 2005 Monash University, Melbourne, Australia GENERAL SOLUTIONS FOR THE INITIAL RUN-UP OF A
More informationTwo-Dimensional Analysis of the Power Transfer between Crossed Laser Beams
Two-Dimensional Analsis of the Power Transfer between Crossed Laser Beams The indirect-drive approach to inertial confinement fusion involves laser beams that cross as the enter the hohlraum. Ionacoustic
More informationVisualizing Non-Equilibrium Flow Simulations using 3-D Velocity Distribution Functions
Visualizing Non-Equilibrium Flow Simulations using 3-D Velocit Distribution Functions A. Venkattraman and A.A. Alexeenko School of Aeronautics & Astronautics, Purdue Universit, West Lafaette IN 797 Abstract.
More informationLECTURE NOTES - VIII. Prof. Dr. Atıl BULU
LECTURE NOTES - VIII «LUID MECHNICS» Istanbul Technical Universit College of Civil Engineering Civil Engineering Department Hdraulics Division CHPTER 8 DIMENSIONL NLYSIS 8. INTRODUCTION Dimensional analsis
More informationP10.5 Water flows down a rectangular channel that is 4 ft wide and 3 ft deep. The flow rate is 15,000 gal/min. Estimate the Froude number of the flow.
P10.5 Water flows down a rectangular channel that is 4 ft wide and ft deep. The flow rate is 15,000 gal/min. Estimate the Froude number of the flow. Solution: Convert the flow rate from 15,000 gal/min
More information!y + v!d.!y = 0. (1.2.4)
L. Pratt J. Whitehead 5/23/06 1.2 The shallow water equations one-dimensional wave propagation. Traditional discussions of hydraulic effects such as those found in engineering text books are often based
More informationExperimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper
Tamkang Journal of Science and Engineering, Vol. 8, No 4, pp. 319-35 (005) 319 Experimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper Jong-Cheng Wu Department of Civil Engineering,
More informationProf. B.S. Thandaveswara. A short horizontal reach of a prismatic channel is considered. Further, the external
Hdraulics 9. Speciic Force short horizontal reach o a prismatic channel is considered. Further, the external rictional orce and the eect o weight component o water can be considered as negligible. Then
More informationGALLOPING OF SMALL ASPECT RATIO SQUARE CYLINDER
OL. NO. JANUARY 5 ISSN 89-668 6-5 Asian Research Publishing Network (ARPN). All rights reserved. GALLOPING OF SMALL ASPET RATIO SQUARE YLINDER A. N. Rabinin and. D. Lusin Facult of Mathematics and Mechanics
More informationExperimental Uncertainty Review. Abstract. References. Measurement Uncertainties and Uncertainty Propagation
Experimental Uncertaint Review Abstract This is intended as a brief summar of the basic elements of uncertaint analsis, and a hand reference for laborator use. It provides some elementar "rules-of-thumb"
More information=1. The thinner contours continue to represent constant Q 2. 1/2 w *], except that w* rather than h* is considered
5.4 Steady flow through a pure contraction. If the bottom remains horizontal (h*=constant=0, say) and the flow is choked only by contractions in the width of the rectangular channel then a new type of
More informationEP225 Note No. 4 Wave Motion
EP5 Note No. 4 Wave Motion 4. Sinusoidal Waves, Wave Number Waves propagate in space in contrast to oscillations which are con ned in limited regions. In describing wave motion, spatial coordinates enter
More informationModelling of Tsunami Waves
MATEMATIKA, 2008, Volume 24, Number 2, 211 230 c Department of Mathematics, UTM. Modelling of Tsunami Waves 1 Nazeeruddin Yaacob, Norhafizah Md Sarif & 2 Zainal Abdul Aziz Department of Mathematics, Faculty
More informationMETHOD OF WATER QUALITY IMPROVEMENT BY USE OF ICE FLOE CONTROL STRUCTURE
METHOD OF WATER QUALITY IMPROVEMENT BY USE OF ICE FLOE CONTROL STRUCTURE Shunsuke Makita 1 and Hiroshi Saeki 1 Graduate School of Engineering, Hokkaido Universit, Sapporo, Japan Hokkaido Universit, Sapporo,
More informationSound Propagation in Ducts
Sound Propagation in Ducts Hongbin Ju Department of Mathematics Florida State Universit, Tallahassee, FL.3306 www.aeroacoustics.info Please send comments to: hju@math.fsu.edu In this section we will discuss
More informationOn the Whitham Equation
On the Whitham Equation Henrik Kalisch Department of Mathematics University of Bergen, Norway Joint work with: Handan Borluk, Denys Dutykh, Mats Ehrnström, Daulet Moldabayev, David Nicholls Research partially
More informationGoals of this Chapter
Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature in the presence of positive static stability internal gravity waves Conservation of potential vorticity in the presence
More informationThe Gas-assisted Expelled Fluid Flow in the Front of a Long Bubble in a Channel
C. H. Hsu, P. C. Chen,. Y. ung, G. C. uo The Gas-assisted Epelled Fluid Flow in the Front of a Long Bubble in a Channel C.H. Hsu 1, P.C. Chen 3,.Y. ung, and G.C. uo 1 1 Department of Mechanical Engineering
More informationUnsteady free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion *
CHAPTER 4 Unstead free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion * 4. Nomenclature A : Suction parameter C : Dimensionless species
More informationDIRECT NUMERICAL SIMULATION OF AN AIRFOIL WITH UNSTEADY WAKE (ECCOMAS 2006)
European Conference on Computational Fluid Dnamics ECCOMAS CFD 26 P. Wesseling, E. Oñate, J. Périau (Eds) TU Delft, The Netherlands, 26 DIRECT NUMERICAL SIMULATION OF AN AIRFOIL WITH UNSTEADY WAKE (ECCOMAS
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More information!y + v!d.!y = 0. (1.2.4)
L. Pratt J. Whitehead 9/13/05 1.2 The shallow water equations one-dimensional wave propagation. Traditional discussions of hydraulic effects such as those found in engineering text books are often based
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1
AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the
More informationGeneration of solitary waves by transcritical flow over a step. Citation Journal Of Fluid Mechanics, 2007, v. 587, p
Title Generation of solitary waves by transcritical flow over a step Author(s) Grimshaw, RHJ; Zhang, DH; Chow, KW Citation Journal Of Fluid Mechanics, 27, v. 587, p. 235-254 Issued Date 27 URL http://hdl.handle.net/1722/57194
More informationLECTURE NOTES - III. Prof. Dr. Atıl BULU
LECTURE NOTES - III «FLUID MECHANICS» Istanbul Technical University College of Civil Engineering Civil Engineering Department Hydraulics Division CHAPTER KINEMATICS OF FLUIDS.. FLUID IN MOTION Fluid motion
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationDepartment of Mathematics, Usmanu Danfodiyo University, Sokoto, Nigeria. DOI: / Page
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. II (Jan - Feb. 5), PP 4-58 www.iosrjournals.org Finite difference solutions of magneto hdrodnamic free convective
More informationSoliton Molecules. Fedor Mitschke Universität Rostock, Institut für Physik. Benasque, October
Soliton Soliton Molecules Molecules and and Optical Optical Rogue Rogue Waves Waves Benasque, October 2014 Fedor Mitschke Universität Rostock, Institut für Physik fedor.mitschke@uni-rostock.de Part II
More informationThomas Pierro, Donald Slinn, Kraig Winters
Thomas Pierro, Donald Slinn, Kraig Winters Department of Ocean Engineering, Florida Atlantic University, Boca Raton, Florida Applied Physics Laboratory, University of Washington, Seattle, Washington Supported
More informationQUANTITATIVE ANALYSIS ON THE TSUNAMI INUNDATION AREA AND EROSION ALONG THE ABUKUMAGAWA RIVER AT THE GREAT EAST JAPAN EARTHQUAKE
UANTITATIVE ANALYSIS ON THE TSUNAMI INUNDATION AREA AND EROSION ALONG THE ABUKUMAGAWA RIVER AT THE GREAT EAST JAPAN EARTHUAKE Masauki Sato Graduate School of Science and Engineering Saitama Universit,
More informationPhase diagram of interfacial instabilities in a two-layer Couette flow and mechanism of the long-wave instability
J. Fluid Mech. (), vol. 44, pp. 95 3. Printed in the United Kingdom c Cambridge Universit Press 95 Phase diagram of interfacial instabilities in a two-laer Couette flow and mechanism of the long-wave instabilit
More informationFIRST- AND SECOND-ORDER IVPS The problem given in (1) is also called an nth-order initial-value problem. For example, Solve: Solve:
.2 INITIAL-VALUE PROBLEMS 3.2 INITIAL-VALUE PROBLEMS REVIEW MATERIAL Normal form of a DE Solution of a DE Famil of solutions INTRODUCTION We are often interested in problems in which we seek a solution
More informationACCEPTED VERSION AIP Publishing. Originally published at:
ACCEPTED VERSION D. J. Arnold, Y. M. Stokes, and J. E. F. Green Thin-film flow in helicall wound shallow channels of arbitrar cross-sectional shape Phsics of Fluids, 217; 29(1): 1312-1-1312-11 217 AIP
More informationChemical Reaction and Viscous Dissipation Effects on an Unsteady MHD Flow of Heat and Mass Transfer along a Porous Flat Plate
International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 2 Issue 2 December 24 PP 977-988 ISSN 2347-37X (Print) & ISSN 2347-342 (Online) www.arcjournals.org Chemical Reaction
More informationLateral Inflow into High-Velocity Channels
Lateral Inflow into High-Velocity Channels by Richard L. Stockstill PURPOSE: This Coastal and Hydraulics Engineering Technical Note (CHETN) investigates lateral flow discharging into a high-velocity channel.
More informationExperiments on capillary-gravity waves of solitary type on deep water
Experiments on capillary-gravity waves of solitary type on deep water Michael Longuet-Higgins Institute for Nonlinear Science, University of California San Diego, La Jolla, California 92093-0402 Xin Zhang
More informationExperimental and numerical investigation of 2D sloshing with slamming
Experimental numerical investigation of 2D sloshing with slamming A. Colagrossi C. Lugni M. Greco O. M. Faltinsen a.colagrossi@insean.it c.lugni@insean.it m.greco@insean.it oddfal@marin.ntnu.no INSEAN,
More informationCONSERVATION OF ENERGY FOR ACONTINUUM
Chapter 6 CONSERVATION OF ENERGY FOR ACONTINUUM Figure 6.1: 6.1 Conservation of Energ In order to define conservation of energ, we will follow a derivation similar to those in previous chapters, using
More information