scholarly journals Numerical Modeling of Tsunami Waves Interaction with Porous and Impermeable Vertical Barriers

2012 ◽  
Vol 2012 ◽  
pp. 1-27 ◽  
Author(s):  
Manuel del Jesus ◽  
Javier L. Lara ◽  
Inigo J. Losada

Tsunami wave interaction with coastal regions is responsible for very important human and economic losses. In order to properly design coastal defenses against these natural catastrophes, new numerical models need to be developed that complement existing laboratory measurements and field data. The use of numerical models based on the Navier-Stokes equations appears as a reasonable approach due to their ability to evaluate complex flow patterns around coastal structures without the inherent limitations of the classical depth-averaged models. In the present study, a Navier-Stokes-based model, IH-3VOF, is applied to study the interaction of tsunami waves with porous and impermeable structures. IH-3VOF is able to simulate wave flow within the porous structures by means of the volume-averaged Reynolds-averaged Navier-Stokes (VARANS) equations. The equations solved by the model and their numerical implementation are presented here. A numerical analysis of the interaction of a tsunami wave with both an impermeable and porous vertical breakwater is carried out. The wave-induced three-dimensional wave pattern is analysed from the simulations. The role paid by the porous media is also investigated. Finally, flow around the breakwater is analyzed identifying different flow behaviors in the vicinity of the breakwater and in the far field of the structure.

2021 ◽  
Vol 929 ◽  
Author(s):  
Ravi Sudam Jadhav ◽  
Abhimanyu Gavasane ◽  
Amit Agrawal

The main goal of the present study is to thoroughly test the recently derived OBurnett equations for the normal shock wave flow problem for a wide range of Mach number ( $3 \leq Ma \leq 9$ ). A dilute gas system composed of hard-sphere molecules is considered and the numerical results of the OBurnett equations are validated against in-house results from the direct simulation Monte Carlo method. The primary focus is to study the orbital structures in the phase space (velocity–temperature plane) and the variation of hydrodynamic fields across the shock. From the orbital structures, we observe that the heteroclinic trajectory exists for the OBurnett equations for all the Mach numbers considered, unlike the conventional Burnett equations. The thermodynamic consistency of the equations is also established by showing positive entropy generation across the shock. Further, the equations give smooth shock structures at all Mach numbers and significantly improve upon the results of the Navier–Stokes equations. With no tweaking of the equations in any way, the present work makes two important contributions by putting forward an improved theory of shock waves and establishing the validity of the OBurnett equations for solving complex flow problems.


Author(s):  
Pierre Ferrant ◽  
Lionel Gentaz ◽  
Bertrand Alessandrini ◽  
Romain Luquet ◽  
Charles Monroy ◽  
...  

This paper documents recent advances of the SWENSE (Spectral Wave Explicit Navier-Stokes Equations) approach, a method for simulating fully nonlinear wave-body interactions including viscous effects. The methods efficiently combines a fully nonlinear potential flow description of undisturbed wave systems with a modified set of RANS with free surface equations accounting for the interaction with a ship or marine structure. Arbitrary incident wave systems may be described, including regular, irregular waves, multidirectional waves, focused wave events, etc. The model may be fixed or moving with arbitrary speed and 6 degrees of freedom motion. The extension of the SWENSE method to 6 DOF simulations in irregular waves as well as to manoeuvring simulations in waves are discussed in this paper. Different illlustative simulations are presented and discussed. Results of the present approach compare favorably with available reference results.


1990 ◽  
Vol 17 (2) ◽  
pp. 252-261 ◽  
Author(s):  
Kevin R. Hall

The interaction of a wave with a rubblemound breakwater results in a complex flow field which is both nonlinear and turbulent, particularly within a region close to the surface of the structure. Numerical models describing internal flow in a rubblemound breakwater are becoming increasingly important, particularly as the influence of scale effects on internal flow in physical hydraulic models are becoming understood as important. A number of numerical models to predict the internal breakwater flow kinematics have been produced in the past two decades. This paper provides a review of the state-of-the-art of numerical modelling of wave interaction with rubblemound breakwaters. Details of the theoretical development and the resulting numerical solution techniques are presented. Methods for incorporating secondary effects such as two-phase (air–water) flow, inertia, and unbalanced boundary conditions are discussed. Limitations of the models resulting from the validity of the assumptions made in order to effect a numerical solution are discussed. Key words: breakwaters, internal flow, porous media flow, numerical modelling, rubblemound breakwaters.


2020 ◽  
Vol 8 (6) ◽  
pp. 419 ◽  
Author(s):  
Yun-Ta Wu ◽  
Shih-Chun Hsiao

In this article, the interaction of solitary waves and a submerged slotted barrier is investigated in which the slotted barrier consists of three impermeable elements and its porosity can be determined by the distance between the two neighboring elements. A new experiment is conducted to measure free surface elevation, velocity, and turbulent kinetic energy. Numerical simulation is performed using a two-dimensional model based on the Reynolds-Averaged Navier-Stokes equations and the non-linear k-ɛ turbulence model. A detailed flow pattern is illustrated by a flow visualization technique. A laboratory observation indicates that flow separations occur at each element of the slotted barrier and the vortex shedding process is then triggered due to the complicated interaction of those induced vortices that further create a complex flow pattern. During the vortex shedding process, seeding particles that are initially accumulated near the seafloor are suspended by an upward jet formed by vortices interacting. Model-data comparisons are carried out to examine the accuracy of the model. Overall model-data comparisons are in satisfactory agreement, but modeled results sometimes fail to predict the positions of the induced vortices. Since the measured data is unique in terms of velocity and turbulence, the dataset can be used for further improvement of numerical modeling.


1969 ◽  
Vol 37 (4) ◽  
pp. 727-750 ◽  
Author(s):  
Gareth P. Williams

A method of numerically integrating the Navier-Stokes equations for certain three-dimensional incompressible flows is described. The technique is presented through application to the particular problem of describing thermal convection in a rotating annulus. The equations, in cylindrical polar co-ordinate form, are integrated with respect to time by a marching process, together with the solving of a Poisson equation for the pressure. A suitable form of the finite difference equations gives a computationally-stable long-term integration with reasonably faithful representation of the spatial and temporal characteristics of the flow.Trigonometric interpolation techniques provide accurate (discretely exact) solutions to the Poisson equation. By using an auxiliary algorithm for rapid evaluation of trigonometric transforms, the proportion of computation needed to solve the Poisson equation can be reduced to less than 25% of the total time needed to’ advance one time step. Computing on a UNIVAC 1108 machine, the flow can be advanced one time-step in 2 sec for a 14 × 14 × 14 grid upward to 96 sec for a 60 × 34 × 34 grid.As an example of the method, some features of a solution for steady wave flow in annulus convection are presented. The resemblance of this flow to the classical Eady wave is noted.


2020 ◽  
Vol 8 (11) ◽  
pp. 903
Author(s):  
Sixtine Neuvéglise ◽  
Gaële Perret ◽  
Hassan Smaoui ◽  
François Marin ◽  
Philippe Sergent

This paper studies the behaviour of a quayside floater oscillating in front of a vertical dike. In order to study the floater motion and the impact of the dike on the floater, a linear analytical model based on 2D potential flow theory in intermediate water depth conditions and a numerical model resolving 2D Navier–Stokes equations are developed. Physical tests performed for different floater dimensions in a wave tank are used as references for the analytical and numerical models. The comparison of the results obtained analytically, numerically and experimentally leads to the validity domain of the potential model. A correction of this model is proposed, based on the optimization of the radiated coefficients, and a quadratic drag term is added according to Morison equation. The impact of the different parameters of the system on floater behaviour is considered. Results show that the draft has the most important impact on floater motion.


2013 ◽  
Vol 721 ◽  
pp. 58-85 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall ◽  
Andrew Walton

AbstractThe recently understood relationship between high-Reynolds-number vortex–wave interaction theory and computationally generated self-sustaining processes provides a possible route to an understanding of some of the underlying structures of fully turbulent flows. Here vortex–wave interaction (VWI) theory is used in the long streamwise wavelength limit to continue the development found at order-one wavelengths by Hall & Sherwin (J. Fluid Mech., vol. 661, 2010, pp. 178–205). The asymptotic description given reduces the Navier–Stokes equations to the so-called boundary-region equations, for which we find equilibrium states describing the change in the VWI as the wavelength of the wave increases from $O(h)$ to $O(Rh)$, where $R$ is the Reynolds number and $2h$ is the depth of the channel. The reduced equations do not include the streamwise pressure gradient of the perturbation or the effect of streamwise diffusion of the wave–vortex states. The solutions we calculate have an asymptotic error proportional to ${R}^{- 2} $ when compared to the full Navier–Stokes equations. The results found correspond to the minimum drag configuration for VWI states and might therefore be of relevance to the control of turbulent flows. The key feature of the new states discussed here is the thickening of the critical layer structure associated with the wave part of the flow to completely fill the channel, so that the roll part of the flow is driven throughout the flow rather than as in Hall & Sherwin as a stress discontinuity across the critical layer. We identify a critical streamwise wavenumber scaling, which, when approached, causes the flow to localize and take on similarities with computationally generated or experimentally observed turbulent spots. In effect, the identification of this critical wavenumber for a given value of the assumed high Reynolds number fixes a minimum box length necessary for the emergence of localized structures. Whereas nonlinear equilibrium states of the Navier–Stokes equations are thought to form a backbone on which turbulent flows hang, our results suggest that the localized states found here might play a related role for turbulent spots.


2008 ◽  
Vol 136 (11) ◽  
pp. 4045-4062 ◽  
Author(s):  
Robert L. Walko ◽  
Roni Avissar

Abstract The dynamic core of the Ocean–Land–Atmosphere Model (OLAM), which is a new global model that is partly based on the Regional Atmospheric Modeling System (RAMS), is described and tested. OLAM adopts many features of its predecessor, but its dynamic core is new and incorporates a global geodesic grid with triangular mesh cells and a finite-volume discretization of the nonhydrostatic compressible Navier–Stokes equations. The spatial discretization of horizontal momentum is based on a C-staggered grid and uses a method that has not been previously applied in atmospheric modeling. The temporal discretization uses a unique form of time splitting that enforces consistency of advecting mass flux among all conservation equations. OLAM grid levels are horizontal, and topography is represented by the shaved-cell method. Aspects of the shaved-cell method that pertain to the OLAM discretization on the triangular mesh are described, and a method of conserving momentum in shaved cells on a C-staggered grid is presented. The dynamic core was tested in simulations with multiple vertical model levels and significant vertical motion. The tests include an idealized global circulation simulation, a cold density current, and mountain-wave flow over an orographic barrier, all of which are well-known standard benchmark experiments. OLAM gave acceptable results in all tests, demonstrating that its dynamic core produces accurate and robust solutions.


2015 ◽  
Vol 772 ◽  
Author(s):  
L. Kahouadji ◽  
N. Périnet ◽  
L. S. Tuckerman ◽  
S. Shin ◽  
J. Chergui ◽  
...  

We report the first simulations of the Faraday instability using the full three-dimensional Navier–Stokes equations in domains much larger than the characteristic wavelength of the pattern. We use a massively parallel code based on a hybrid front-tracking/level-set algorithm for Lagrangian tracking of arbitrarily deformable phase interfaces. Simulations performed in square and cylindrical domains yield complex patterns. In particular, a superlattice-like pattern similar to those of Douady & Fauve (Europhys. Lett., vol. 6, 1988, pp. 221–226) and Douady (J. Fluid Mech., vol. 221, 1990, pp. 383–409) is observed. The pattern consists of the superposition of two square superlattices. We conjecture that such patterns are widespread if the square container is large compared with the critical wavelength. In the cylinder, pentagonal cells near the outer wall allow a square-wave pattern to be accommodated in the centre.


Author(s):  
Hadi Karrabi ◽  
Mohsen Rezasoltani

An investigation to understand the impact of twisted, leaned and bowed blades on the performance of axial turbine was undertaken. A CFD code, which solves the Reynolds-averaged Navier–Stokes equations, was used to compute the complex flow field of axial turbine. The code was validated against existing Hannover turbine experimental data. Numerical data showed good agreement with measured data. Finally, the effect of geometry changes, focusing on blade lean, twist and bow, on the Avon turbine blade performance, was analyzed. Results show that twisted blade affects performance significantly. Leaned and bowed blade has minor effect on performance.


Sign in / Sign up

Export Citation Format

Share Document