scholarly journals Assessment of Numerical Methods for Plunging Breaking Wave Predictions

2021 ◽  
Vol 9 (3) ◽  
pp. 264
Author(s):  
Shanti Bhushan ◽  
Oumnia El Fajri ◽  
Graham Hubbard ◽  
Bradley Chambers ◽  
Christopher Kees

This study evaluates the capability of Navier–Stokes solvers in predicting forward and backward plunging breaking, including assessment of the effect of grid resolution, turbulence model, and VoF, CLSVoF interface models on predictions. For this purpose, 2D simulations are performed for four test cases: dam break, solitary wave run up on a slope, flow over a submerged bump, and solitary wave over a submerged rectangular obstacle. Plunging wave breaking involves high wave crest, plunger formation, and splash up, followed by second plunger, and chaotic water motions. Coarser grids reasonably predict the wave breaking features, but finer grids are required for accurate prediction of the splash up events. However, instabilities are triggered at the air–water interface (primarily for the air flow) on very fine grids, which induces surface peel-off or kinks and roll-up of the plunger tips. Reynolds averaged Navier–Stokes (RANS) turbulence models result in high eddy-viscosity in the air–water region which decays the fluid momentum and adversely affects the predictions. Both VoF and CLSVoF methods predict the large-scale plunging breaking characteristics well; however, they vary in the prediction of the finer details. The CLSVoF solver predicts the splash-up event and secondary plunger better than the VoF solver; however, the latter predicts the plunger shape better than the former for the solitary wave run-up on a slope case.

2002 ◽  
Vol 456 ◽  
pp. 295-318 ◽  
Author(s):  
YING LI ◽  
FREDRIC RAICHLEN

The run-up of non-breaking and breaking solitary waves on a uniform plane beach connected to a constant-depth wave tank was investigated experimentally and numerically. If only the general characteristics of the run-up process and the maximum run-up are of interest, for the case of a breaking wave the post-breaking condition can be simplified and represented as a propagating bore. A numerical model using this bore structure to treat the process of wave breaking and subsequent shoreward propagation was developed. The nonlinear shallow water equations (NLSW) were solved using the weighted essentially non-oscillatory (WENO) shock capturing scheme employed in gas dynamics. Wave breaking and post-breaking propagation are handled automatically by this scheme and ad hoc terms are not required. A computational domain mapping technique was used to model the shoreline movement. This numerical scheme was found to provide a relatively simple and reasonably good prediction of various aspects of the run-up process. The energy dissipation associated with wave breaking of solitary wave run-up (excluding the effects of bottom friction) was also estimated using the results from the numerical model.


2016 ◽  
Vol 139 (2) ◽  
Author(s):  
Luca Marocco ◽  
Andrea Franco

A turbulent convective flow of an incompressible fluid inside a staggered ribbed channel with high blockage at ReH ≈ 4200 is simulated with direct numerical simulation (DNS) and Reynolds-averaged Navier–Stokes (RANS) techniques. The DNS results provide the reference solution for comparison of the RANS turbulence models. The k–ε realizable, k–ω SST, and v2¯–f model are accurately analyzed for their strengths and weaknesses in predicting the flow and temperature field for this geometry. These three models have been extensively used in literature to simulate this configuration and boundary conditions but with discordant conclusions upon their performance. The v2¯–f model performs much better than the k–ε realizable while the k–ω SST model results to be inadequate.


2013 ◽  
Vol 135 (3) ◽  
Author(s):  
C. Chin ◽  
M. Li ◽  
C. Harkin ◽  
T. Rochwerger ◽  
L. Chan ◽  
...  

A numerical study of compressible jet flows is carried out using Reynolds averaged Navier–Stokes (RANS) turbulence models such as k-ɛ and k-ω-SST. An experimental investigation is performed concurrently using high-speed optical methods such as Schlieren photography and shadowgraphy. Numerical and experimental studies are carried out for the compressible impinging at various impinging angles and nozzle-to-wall distances. The results from both investigations converge remarkably well and agree with experimental data from the open literature. From the flow visualizations of the velocity fields, the RANS simulations accurately model the shock structures within the core jet region. The first shock cell is found to be constraint due to the interaction with the bow-shock structure for nozzle-to-wall distance less than 1.5 nozzle diameter. The results from the current study show that the RANS models utilized are suitable to simulate compressible free jets and impinging jet flows with varying impinging angles.


1997 ◽  
Vol 342 ◽  
pp. 141-158 ◽  
Author(s):  
M. J. COOKER ◽  
P. D. WEIDMAN ◽  
D. S. BALE

The collision of a solitary wave, travelling over a horizontal bed, with a vertical wall is investigated using a boundary-integral method to compute the potential fluid flow described by the Euler equations. We concentrate on reporting new results for that part of the motion when the wave is near the wall. The wall residence time, i.e. the time the wave crest remains attached to the wall, is introduced. It is shown that the wall residence time provides an unambiguous characterization of the phase shift incurred during reflection for waves of both small and large amplitude. Numerically computed attachment and detachment times and amplitudes are compared with asymptotic formulae developed using the perturbation results of Su & Mirie (1980). Other features of the flow, including the maximum run-up and the instantaneous wall force, are also presented. The numerically determined residence times are in good agreement with measurements taken from a cine film of solitary wave reflection experiments conducted by Maxworthy (1976).


2007 ◽  
Vol 01 (03) ◽  
pp. 257-272 ◽  
Author(s):  
KAO-SHU HWANG ◽  
YU-HSUAN CHANG ◽  
HWUNG-HWENG HWUNG ◽  
YI-SYUAN LI

The evolution and run-up of breaking solitary waves on plane beaches are investigated in this paper. A series of large-scale experiments were conducted in the SUPER TANK of Tainan Hydraulics Laboratory with three plane beaches of slope 0.05, 0.025 and 0.017 (1:20, 1:40 and 1:60). Solitary waves of which relative wave heights, H/h0, ranged from 0.03 to 0.31 were generated by two types of wave-board displacement trajectory: the ramp-trajectory and the solitary-wave trajectory proposed by Goring (1979). Experimental results show that under the same relative wave height, the waveforms produced by the two generation procedures becomes noticeably different as the waves propagate prior to the breaking point. Meanwhile, under the same relative wave height, the larger the constant water depth is, the larger the dimensionless run-up heights would be. Scale effects associated with the breaking process are discussed.


2020 ◽  
Author(s):  
Richard Marcer ◽  
Camille Journeau ◽  
Kévin Pons

<p>This work has been performed within the framework of the TANDEM project (Tsunamis in northern AtlaNtic: Definition of Effects by Modelling) which is dedicated to the appraisal of coastal effects due to tsunami waves on the French coastlines. One of the identified objectives of TANDEM consisted in designing, adapting and validating numerical codes for tsunami hazard assessment, addressing the various stages of a tsunami event: generation, propagation, run-up and coastal inundation.</p><p>PRINCIPIA has been working on the development and qualification of two in-house CFD software’s: a 2D Saint-Venant model (often called NLSW for Non-Linear Shallow Water) using an Adaptive Mesh Refinement (AMR) for simulation of large scale tsunami propagation from the source up to coastal scale, and a 3D Navier-Stokes model dedicated to tsunami coastal impact modelling.</p><p>An overview of the results obtained with both codes aiming at being applicable to tsunami modelling, is presented. The validation process has been done on several academic test cases having experimental data for comparisons, as the breaking of a solitary wave on a reef, the generation of a long wave induced by a vertical bloc (massive cliffs, ice bodies) falling down an underlying water volume, the tsunami generation due to a submarine landslide and the tsunami impact on a coastal city.</p><p>A real case simulation is concerned as well, the devastating 2011 Tohoku event which is compared with in-situ data.</p><p>The work was supported by the Tandem project in the frame of French PIA grant ANR-11-RSNR-00023.</p>


2016 ◽  
Vol 40 (3) ◽  
pp. 317-329 ◽  
Author(s):  
Mustafa Kemal Isman

The turbulent flow over backward-facing step (BFS) is numerically investigated by using FLUENT® code. Both uniform and non-uniform velocity profiles are used as inlet boundary condition. Five different Reynolds averaged Navier–Stokes (RANS) turbulence models are employed. The Std. k–ω model shows the best agreement with the experimental data among the models used under the conditions considered in this study. The results show that using a uniform velocity profile has a negative effect on predictions if the domain is not sufficiently extended upstream from the inlet. To eliminate this effect, the domain should be extended upstream by about 10Dh from the inlet. However, results show that this extension causes absorption effects of inlet parameters such as inlet turbulence intensity.


Author(s):  
Gu¨nther F. Clauss ◽  
Robert Stu¨ck ◽  
Florian Stempinski ◽  
Christian E. Schmittner

For the analysis of loads and motions of marine structures in harsh seaways precise information about the hydrodynamics of waves is required. While the surface motion of waves can easily be measured in physical wave tanks other critical characteristics such as the instantaneous particle velocity and acceleration as well as the pressure field, especially under the wave crest are difficult and time-consuming to obtain. Therefore a new method is presented to approximate the wave potential of a given instantaneous wave contour. Numerical methods — so called numerical wave tanks (NWTs) — are developed to provide the desired insight into wave hydrodynamics. A potential theory method based on the Finite Element method (Pot/FE), a RANSE (Reynolds-Averaged Navier-Stokes Equations) method applying VOF (Volume of Fluid) and a combination of both is utilized for the simulation of different model wave trains. The coupling of both CFD (computational fluid dynamics) solvers is a useful approach to benefit from the advantages of the two different methods: The Pot/FE solver WAVETUB (wave simulation code developed at Technical University Berlin) allows a very fast and accurate simulation of the propagation of nonbreaking waves while the RANSE/VOF solver has the capability of simulating breaking waves. Two different breaking criteria for the detection of wave breaking are implemented in WAVETUB for triggering the automated coupling process by data transfer at the interface. It is shown that an efficient method for the simulation of breaking wave trains including wave-structure interaction in 2D and 3D is established by the coupling of both CFD codes. All results are discussed in detail.


2011 ◽  
Vol 8 (2) ◽  
pp. 10 ◽  
Author(s):  
K. Smida ◽  
H. Lamloumi ◽  
K. Maalel ◽  
Z. Hafsia

 A new numerical wave generation method is used to investigate the head-on collision of two solitary waves. The reflection at vertical wall of a solitary wave is also presented. The originality of this model, based on the Navier-Stokes equations, is the specification of an internal inlet velocity, defined as a source line within the computational domain for the generation of these non linear waves. This model was successfully implemented in the PHOENICS (Parabolic Hyperbolic Or Elliptic Numerical Integration Code Series) code. The collision of two counter-propagating solitary waves is similar to the interaction of a soliton with a vertical wall. This wave generation method allows the saving of considerable time for this collision process since the counter-propagating wave is generated directly without reflection at vertical wall. For the collision of two solitary waves, numerical results show that the run-up phenomenon can be well explained, the solution of the maximum wave run-up is almost equal to experimental measurement. The simulated wave profiles during the collision are in good agreement with experimental results. For the reflection at vertical wall, the spatial profiles of the wave at fixed instants show that this problem is equivalent to the collision process. 


Sign in / Sign up

Export Citation Format

Share Document