Numerical Simulations of Tsunami Wave Generation by Submarine and Aerial Landslides Using RANS and SPH Models

2009 ◽  
Author(s):  
Kaushik Das ◽  
Ron Janetzke ◽  
Debashis Basu ◽  
Steve Green ◽  
John Stamatakos
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Numerical Simulations ◽  
Tsunami Wave ◽  
Wave Generation ◽  
Numerical Results ◽  
Navier Stokes ◽  
Two Dimensional ◽  
Tsunami Waves ◽  
Time Histories

Tsunami wave generation by submarine and aerial landslides is examined in this paper. Two different two-dimensional numerical methods have been used to simulate the time histories of fluid motion, free surface deformation, shoreline movement, and wave runup from tsunami waves generated by aerial and submarine landslides. The first approach is based on the Navier-Stokes equation and the volume of fluid (VOF) method: the Reynolds Averaged Navier-Stokes (RANS)-based turbulence model simulates turbulence, and the VOF method tracks the free surface locations. The second method uses Smoothed Particle Hydrodynamics (SPH)—a numerical model based on a fully Lagrangian approach. In the current work, two-dimensional numerical simulations are carried out for a freely falling wedge representing the landslide and subsequent wave generations. Numerical simulations for the landslide-driven tsunami waves have been performed with different values of landslide material densities. Numerical results obtained from both approaches are compared with experimental data. Simulated results for both aerial and submerged landslides show the complex flow patterns in terms of the velocity field, shoreline evolution, and free-surface profiles. Flows are found to be strongly transient, rotational, and turbulent. Predicted numerical results for time histories of free-surface fluctuations and the runup/rundown at various locations are in good agreement with the available experimental data. The similarity and discrepancy between the solutions obtained by the two approaches are explored and discussed.

Numerical Simulations of Regular and Irregular Wave Forces on a Horizontal Semi-Submerged Cylinder

2017 ◽  
Author(s):  
Shengnan Liu ◽  
Muk Chen Ong ◽  
Charlotte Obhrai ◽  
Sopheak Seng
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Numerical Simulations ◽  
Stokes Equations ◽  
Wave Length ◽  
Wave Structure ◽  
Numerical Results ◽  
Irregular Waves ◽  
Wave Forces ◽  
Submerged Cylinder

Two-dimensional (2D) numerical simulations have been performed using OpenFOAM (an open source CFD software package [1]) and waves2Foam (an OpenFOAM based add-on library for wave generations and absorption [2]) to investigate free surface waves past one fixed horizontally semi-submerged cylinder. The 2-D simulations are carried out by solving Navier-Stokes equations which are discretized based on finite volume method (FVM). Volume of Fluid (VOF) method is employed to capture the free surface in the numerical wave tank. Validation studies have been performed by comparing the numerical results of Stokes first-order wave past a semi-submerged circular cylinder with the published experimental data at different incident wave properties. The numerical results are in good agreement with the experimental data. Subsequently, regular and irregular waves past semi-submerged cylinder at different wave heights and the wave lengths are computed numerically to investigate the effect of the wave height and wave length on wave-structure interaction. The numerical results for irregular waves are compared with those induced by regular waves.

Numerical Simulation of Irregular Wave Forces on a Horizontal Cylinder

2016 ◽  
Author(s):  
Ankit Aggarwal ◽  
Mayilvahanan Alagan Chella ◽  
Arun Kamath ◽  
Hans Bihs ◽  
Øivind Asgeir Arnsten
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Numerical Model ◽  
Circular Cylinder ◽  
Wave Generation ◽  
Numerical Results ◽  
Irregular Waves ◽  
Wave Forces ◽  
Ghost Cell ◽  
Irregular Wave

In the present study, the irregular wave forces on a fully submerged circular cylinder are investigated using the open-source computational fluid dynamics (CFD) model REEF3D. A complete three dimensional representation of the ocean waves requires the consideration of the sea surface as an irregular wave train with the random characteristics. The numerical model uses the incompressible Reynolds-averaged Navier-Stokes (RANS) equations together with the continuity equation to solve the fluid flow problem. Turbulence modeling is carried out using the two equation k-ω model. Spatial discretization is done using an uniform Cartesian grid. The level set method is used for computing the free surface. For time discretization, third-order total variation diminishing (TVD) Runge Kutta scheme is used. Ghost cell boundary method is used for implementing the complex geometries in the numerical model. MPI is used for the exchange of the value of a ghost cell. Relaxation method is used for the wave generation. The numerical model is validated for the irregular waves for a wave tank without any structure. Further, the numerical model is validated by comparing the numerical results with the experimental data for a fully submerged circular cylinder under regular waves and irregular waves. The numerical results are in a good agreement with the experimental data for the regular and irregular wave forces. The JONSWAP spectrum is used for the wave generation. The free surface features and kinematics around the cylinder is also presented and discussed.

Numerical Simulation of the Oscillating Water Column inside a Vertical Cylinder in Regular Waves Using IHFOAM

2017 ◽  
Vol 372 ◽  
pp. 40-49
Author(s):  
J.M.P. Conde ◽  
J.M.S.M. Cardoso
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Water Column ◽  
Vertical Cylinder ◽  
Cylindrical Tube ◽  
Numerical Results ◽  
Navier Stokes ◽  
Oscillating Water Column ◽  
Good Correspondence ◽  
Ocean Wave Energy

The motion of the free-surface inside a surface-piercing vertical cylindrical tube can be seen as a simplified approximation of an oscillating-water-column ocean-wave energy-converter (OWC-OWEC). In the present work the IHFOAM code, which solves the Reynolds-averaged Navier-Stokes (RANS) equations by a finite-volume method using the volume-of-fluid (VoF) technique, was used to simulate the action of a regular wave on the free-surface inside a vertical cylinder open to the atmosphere. In this paper the results obtained by IHFOAM are compared with other numerical results and with experimental data, showing a good correspondence between these results. In this way the IHFOAM code was verified (by comparison with other numerical results) and validated (by comparison with experimental data).

Simulation of the Gap Resonance Between Two Rectangular Barges in Regular Waves by a Free Surface Viscous Flow Solver

2013 ◽  
Author(s):  
B. Elie ◽  
G. Reliquet ◽  
P.-E. Guillerm ◽  
O. Thilleul ◽  
P. Ferrant ◽  
...  
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Viscous Flow ◽  
Stokes Equations ◽  
Resonance Phenomenon ◽  
Navier Stokes ◽  
Regular Waves ◽  
Navier Stokes Equations ◽  
Flow Solver ◽  
Gap Resonance

This paper compares numerical and experimental results in the study of the resonance phenomenon which appears between two side-by-side fixed barges for different sea-states. Simulations were performed using SWENSE (Spectral Wave Explicit Navier-Stokes Equations) approach and results are compared with experimental data on two fixed barges with different headings and bilges. Numerical results, obtained using the SWENSE approach, are able to predict both the frequency and the magnitude of the RAO functions.

Nonlinear Wave Crest Distribution on a Vertical Breakwater

2018 ◽  
Author(s):  
Valentina Laface ◽  
Giovanni Malara ◽  
Felice Arena ◽  
Ioannis A. Kougioumtzoglou ◽  
Alessandra Romolo
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Vertical Wall ◽  
Surface Elevation ◽  
Wave Crest ◽  
Height Distribution ◽  
Free Surface Elevation ◽  
Pressure Transducers ◽  
Pressure Time ◽  
Time Histories

The paper addresses the problem of deriving the nonlinear, up to the second order, crest wave height probability distribution in front of a vertical wall under the assumption of finite spectral bandwidth, finite water depth and long-crested waves. The distribution is derived by relying on the Quasi-Deterministic representation of the free surface elevation in front of the vertical wall. The theoretical results are compared against experimental data obtained by utilizing a compressive sensing algorithm for reconstructing the free surface elevation in front of the wall. The reconstruction is pursued by starting from recorded wave pressure time histories obtained by utilizing a row of pressure transducers located at various levels. The comparison shows that there is an excellent agreement between the proposed distribution and the experimental data and confirm the deviation of the crest height distribution from the Rayleigh one.

scholarly journals Linear pulse structure and signalling in a film flow on an inclined plane

1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES
Keyword(s):  
Free Surface ◽  
Saddle Point ◽  
Convective Instability ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Numerical Results ◽  
Navier Stokes ◽  
Surface Boundary ◽  
Navier Stokes Equations ◽  
Surface Boundary Conditions

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.

Thruster and FPSO Hull Interaction

2014 ◽  
Author(s):  
Ye Tian ◽  
Spyros A. Kinnas
Keyword(s):  
Experimental Data ◽  
Hybrid Method ◽  
Computational Efficiency ◽  
Vortex Lattice ◽  
Numerical Results ◽  
Navier Stokes ◽  
Dynamic Positioning ◽  
Lattice Method ◽  
Vortex Lattice Method ◽  
Number Of Cells

A hybrid method which couples a Vortex-Lattice Method (VLM) solver and a Reynolds-Averaged Navier-Stokes (RANS) solver is applied to simulate the interaction between a Dynamic Positioning (DP) thruster and an FPSO hull. The hybrid method could significantly reduce the number of cells to fifth of that in a full blown RANS simulation and thus greatly enhance the computational efficiency. The numerical results are first validated with available experimental data, and then used to assess the significance of the thruster/hull interaction in DP systems.

On the Hydrodynamic Interaction Between Ship and Free-Surface Motions on Vessels With Moonpools

2019 ◽  
Author(s):  
Senthuran Ravinthrakumar ◽  
Trygve Kristiansen ◽  
Babak Ommani
Keyword(s):  
Free Surface ◽  
Flow Separation ◽  
Coupling Effect ◽  
Three Dimensional ◽  
Dimensional Case ◽  
Navier Stokes ◽  
Linear Potential ◽  
Two Dimensional ◽  
Sloshing Mode ◽  
Vessel Motion

Abstract Coupling between moonpool resonance and vessel motion is investigated in two-dimensional and quasi three-dimensional settings, where the models are studied in forced heave and in freely floating conditions. The two-dimensional setups are with a recess, while the quasi three-dimensional setups are without recess. One configuration with recess is presented for the two-dimensional case, while three different moonpool sizes (without recess) are tested for the quasi three-dimensional setup. A large number of forcing periods, and three wave steepnesses are tested. Boundary Element Method (BEM) and Viscous BEM (VBEM) time-domain codes based on linear potential flow theory, and a Navier–Stokes solver with linear free-surface and body-boundary conditions, are implemented to investigate resonant motion of the free-surface and the model. Damping due to flow separation from the sharp corners of the moonpool inlets is shown to matter for both vessel motions and moonpool response around the piston mode. In general, the CFD simulations compare well with the experimental results. BEM over-predicts the response significantly at resonance. VBEM provides improved results compared to the BEM, but still over-predicts the response. In the two-dimensional study there are significant coupling effects between heave, pitch and moonpool responses. In the quasi three-dimensional tests, the coupling effect is reduced significantly as the moonpool dimensions relative to the displaced volume of the ship is reduced. The first sloshing mode is investigated in the two-dimensional case. The studies show that damping due to flow separation is dominant. The vessel motions are unaffected by the moonpool response around the first sloshing mode.

scholarly journals Numerical modeling of fluid resonance in narrow gaps of fixed rectangular structures

2019 ◽  
Vol 11 (8) ◽  
pp. 168781401987230
Author(s):  
Ming-ming Liu ◽  
Rui-jia Jin ◽  
Zhen-dong Cui
Keyword(s):  
Free Surface ◽  
Numerical Model ◽  
Stokes Equations ◽  
Theoretical Data ◽  
Numerical Results ◽  
Two Dimensional ◽  
Arbitrary Lagrangian Eulerian ◽  
Fluid Resonance ◽  
Resonance Wave ◽  
Narrow Gaps

A two-dimensional numerical model is developed to investigate the phenomenon of resonance in narrow gaps. Instead of using commonly used Volume of Fluid method to capture the free surface which is sometimes difficult to capture the geometric properties of the geometrically complicated interface, the free surface is traced by using Arbitrary Lagrangian–Eulerian method. The numerical model is based on the two-dimensional Reynolds-Averaged Navier–Stokes equations. The numerical model is validated against wave propagation in wave flume. Comparisons between the numerical results and available theoretical data show satisfactory agreements. Fluid resonance in narrow gaps of fixed rectangular structures are simulated. Numerical results show that resonance wave height and wave frequency for rectangle boxes with sphenoid corners is larger than for rectangle boxes.

Sign in / Sign up

Export Citation Format

Share Document