scholarly journals Variational modelling of extreme waves through oblique interaction of solitary waves

2016 ◽  
Author(s):  
Floriane Gidel ◽  
Onno Bokhove

Abstract. In this work, we model extreme waves that occur due to Mach reflection through the intersection of two obliquely incident solitary waves. For a given range of incident angles and amplitudes, the Mach stem wave grows linearly in length and amplitude, reaching up to four times the amplitude of the incident waves. A variational approach is used to derive the bidirectional Benney-Luke equations, an asymptotic equivalent of the three-dimensional potential-flow equations modelling water waves. This nonlinear and dispersive model has the advantage of allowing wave propagation in two horizontal directions, which is not the case with the unidirectional Kadomtsev-Petviashvili (KP) equation used in most previous studies. A variational Galerkin finite element method is applied to solve the system numerically in Firedrake with a second-order Stormer-Verlet temporal integration scheme in order to obtain stable simulations that conserve the overall mass and energy of the system. Using this approach, we are able to get close to the fourfold amplitude amplification predicted by Miles.

2017 ◽  
Vol 24 (1) ◽  
pp. 43-60 ◽  
Author(s):  
Floriane Gidel ◽  
Onno Bokhove ◽  
Anna Kalogirou

Abstract. In this work, we model extreme waves that occur due to Mach reflection through the intersection of two obliquely incident solitary waves. For a given range of incident angles and amplitudes, the Mach stem wave grows linearly in length and amplitude, reaching up to 4 times the amplitude of the incident waves. A variational approach is used to derive the bidirectional Benney–Luke equations, an asymptotic equivalent of the three-dimensional potential-flow equations modelling water waves. This nonlinear and weakly dispersive model has the advantage of allowing wave propagation in two horizontal directions, which is not the case with the unidirectional Kadomtsev–Petviashvili (KP) equation used in most previous studies. A variational Galerkin finite-element method is applied to solve the system numerically in Firedrake with a second-order Störmer–Verlet temporal integration scheme, in order to obtain stable simulations that conserve the overall mass and energy of the system. Using this approach, we are able to get close to the 4-fold amplitude amplification predicted by Miles.


Author(s):  
S. V. Subramanian ◽  
R. Bozzola ◽  
Louis A. Povinelli

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.


2012 ◽  
Vol 708 ◽  
pp. 480-501 ◽  
Author(s):  
Zhan Wang ◽  
Paul A. Milewski

AbstractThe dynamics of solitary gravity–capillary water waves propagating on the surface of a three-dimensional fluid domain is studied numerically. In order to accurately compute complex time-dependent solutions, we simplify the full potential flow problem by using surface variables and taking a particular cubic truncation possessing a Hamiltonian with desirable properties. This approximation agrees remarkably well with the full equations for the bifurcation curves, wave profiles and the dynamics of solitary waves for a two-dimensional fluid domain, and with higher-order truncations in three dimensions. Fully localized solitary waves are then computed in the three-dimensional problem and the stability and interaction of both line and localized solitary waves are investigated via numerical time integration of the equations. There are many solitary wave branches, indexed by their finite energy as their amplitude tends to zero. The dynamics of the solitary waves is complex, involving nonlinear focusing of wavepackets, quasi-elastic collisions, and the generation of propagating, spatially localized, time-periodic structures akin to breathers.


2001 ◽  
Vol 434 ◽  
pp. 301-335 ◽  
Author(s):  
R. PORTER ◽  
D. PORTER

The scattering and trapping of water waves by three-dimensional submerged topography, infinite and periodic in one horizontal coordinate and of finite extent in the other, is considered under the assumptions of linearized theory. The mild-slope approximation is used to reduce the governing boundary value problem to one involving a form of the Helmholtz equation in which the coefficient depends on the topography and is therefore spatially varying.Two problems are considered: the scattering by the topography of parallel-crested obliquely incident waves and the propagation of trapping modes along the periodic topography. Both problems are formulated in terms of ‘domain’ integral equations which are solved numerically.Trapped waves are found to exist over any periodic topography which is ‘sufficiently’ elevated above the unperturbed bed level. In particular, every periodic topography wholly elevated above that level supports trapped waves. Fundamental differences are shown to exist between these trapped waves and the analogous Rayleigh–Bloch waves which exist on periodic gratings in acoustic theory.Results computed for the scattering problem show that, remarkably, there exist zeros of transmission at discrete wavenumbers for any periodic bed elevation and for all incident wave angles. One implication of this property is that total reflection of an incident wave of a particular frequency will occur in a channel with a single symmetric elevation on the bed. The zeros of transmission in the scattering problem are shown to be related to the presence of a ‘nearly trapped’ mode in the corresponding homogeneous problem.The scattering of waves by multiple rows of periodic topography is also considered and it is shown how Bragg resonance – well-established in scattering of waves by two-dimensional ripple beds – occurs in modes other than the input mode.


Author(s):  
Paul A. Milewski ◽  
Zhan Wang

When both gravity and surface tension effects are present, surface solitary water waves are known to exist in both two- and three-dimensional infinitely deep fluids. We describe here solutions bridging these two cases: travelling waves which are localized in the propagation direction and periodic in the transverse direction. These transversally periodic gravity–capillary solitary waves are found to be of either elevation or depression type, tend to plane waves below a critical transverse period and tend to solitary lumps as the transverse period tends to infinity. The waves are found numerically in a Hamiltonian system for water waves simplified by a cubic truncation of the Dirichlet-to-Neumann operator. This approximation has been proved to be very accurate for both two- and three-dimensional computations of fully localized gravity–capillary solitary waves. The stability properties of these waves are then investigated via the time evolution of perturbed wave profiles.


2021 ◽  
Vol 2021 ◽  
pp. 1-18
Author(s):  
Ghazala Akram ◽  
Naila Sajid ◽  
Muhammad Abbas ◽  
Y. S. Hamed ◽  
Khadijah M. Abualnaja

In this study, the solutions of 2 + 1 -dimensional nonlinear Date–Jimbo–Kashiwara–Miwa (DJKM) equation are characterized, which can be used in mathematical physics to model water waves with low surface tension and long wavelengths. The integration scheme, namely, the extended direct algebraic method, is used to extract complex trigonometric, rational and hyperbolic functions. The complex-valued solutions represent traveling waves in different structures, such as bell-, V-, and W-shaped multiwaves. The results obtained in this article are novel and more general than those contained in the literature (Wang et al., 2014, Yuan et al., 2017, Pu and Hu 2019, Singh and Gupta 2018). Furthermore, the mechanical features and dynamical characteristics of the obtained solutions are demonstrated by three-dimensional graphics.


Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


2020 ◽  
Vol 2020 (1) ◽  
Author(s):  
Marin Marin ◽  
M. M. Bhatti

AbstractThe present study deals with the head-on collision process between capillary–gravity solitary waves in a finite channel. The present mathematical modeling is based on Nwogu’s Boussinesq model. This model is suitable for both shallow and deep water waves. We have considered the surface tension effects. To examine the asymptotic behavior, we employed the Poincaré–Lighthill–Kuo method. The resulting series solutions are given up to third-order approximation. The physical features are discussed for wave speed, head-on collision profile, maximum run-up, distortion profile, the velocity at the bottom, and phase shift profile, etc. A comparison is also given as a particular case in our study. According to the results, it is noticed that the free parameter and the surface tension tend to decline the solitary-wave profile significantly. However, the maximum run-up amplitude was affected in great measure due to the surface tension and the free parameter.


2003 ◽  
Vol 87 (510) ◽  
pp. 468-476
Author(s):  
Paul R. Turner

Water waves are familiar to all of us and we encounter them in a variety of guises in many places, be it crashing to shore at the beach, rippling concentrically outward where a pebble lands in a pond or simply splashing at the sides of the bath. The study of waves can be simplified by idealising them as graphs, each graph being thought of as a cross-section of a physical wave at an instant in time. A sequence of such graphs can represent the progress of the wave as time passes.


Author(s):  
W. T. Tiow ◽  
M. Zangeneh

The development and application of a three-dimensional inverse methodology is presented for the design of turbomachinery blades. The method is based on the mass-averaged swirl, rV~θ distribution and computes the necessary blade changes directly from the discrepancies between the target and initial distributions. The flow solution and blade modification converge simultaneously giving the final blade geometry and the corresponding steady state flow solution. The flow analysis is performed using a cell-vertex finite volume time-marching algorithm employing the multistage Runge-Kutta integrator in conjunction with accelerating techniques (local time stepping and grid sequencing). To account for viscous effects, dissipative forces are included in the Euler solver using the log-law and mixing length models. The design method can be used with any existing solver solving the same flow equations without any modifications to the blade surface wall boundary condition. Validation of the method has been carried out using a transonic annular turbine nozzle and NASA rotor 67. Finally, the method is demonstrated on the re-design of the blades.


Sign in / Sign up

Export Citation Format

Share Document