A new model of shoaling and breaking waves. Part 2. Run-up and two-dimensional waves

2019 ◽  
Vol 867 ◽  
pp. 146-194 ◽  
Author(s):  
G. L. Richard ◽  
A. Duran ◽  
B. Fabrèges

We derive a two-dimensional depth-averaged model for coastal waves with both dispersive and dissipative effects. A tensor quantity called enstrophy models the subdepth large-scale turbulence, including its anisotropic character, and is a source of vorticity of the average flow. The small-scale turbulence is modelled through a turbulent-viscosity hypothesis. This fully nonlinear model has equivalent dispersive properties to the Green–Naghdi equations and is treated, both for the optimization of these properties and for the numerical resolution, with the same techniques which are used for the Green–Naghdi system. The model equations are solved with a discontinuous Galerkin discretization based on a decoupling between the hyperbolic and non-hydrostatic parts of the system. The predictions of the model are compared to experimental data in a wide range of physical conditions. Simulations were run in one-dimensional and two-dimensional cases, including run-up and run-down on beaches, non-trivial topographies, wave trains over a bar or propagation around an island or a reef. A very good agreement is reached in every cases, validating the predictive empirical laws for the parameters of the model. These comparisons confirm the efficiency of the present strategy, highlighting the enstrophy as a robust and reliable tool to describe wave breaking even in a two-dimensional context. Compared with existing depth-averaged models, this approach is numerically robust and adds more physical effects without significant increase in numerical complexity.

2013 ◽  
Vol 732 ◽  
pp. 47-76 ◽  
Author(s):  
Simo A. Mäkiharju ◽  
Brian R. Elbing ◽  
Andrew Wiggins ◽  
Sarah Schinasi ◽  
Jean-Marc Vanden-Broeck ◽  
...  

AbstractThe behaviour of a nominally two-dimensional ventilated partial cavity was examined over a wide range of size scales and flow speeds to determine the influence of Froude, Reynolds, and Weber number on the cavity shape, dynamics, and gas entrainment rate. Two geometrically similar experiments were conducted with a 14:1 length scale ratio. The results were compared to a two-dimensional semi-analytical model of the cavity flow, and Froude scaling was found to be sufficient to match basic cavity shapes. However, the air flux required to maintain a stable cavity did not scale with Froude number alone, as the dynamics of the cavity closure changed with increasing Reynolds number. The required air flux differed over one order of magnitude between the lowest and highest Reynolds number flows. But, for sufficiently high Reynolds numbers, the rate of scaled entrainment appeared to approach Reynolds number independence. Modest changes in surface tension of the small-scale experiment suggested that the Weber number was important only at the lowest speeds and smaller length scale. Otherwise, the Weber numbers of the flows were sufficiently high to make the effects of interfacial tension negligible. We also observed that modest unsteadiness in the inflow to the large-scale cavity led to a significant increase in the required air flux needed to maintain a stable cavity, with the required excess gas flux nominally proportional to the flow’s perturbation amplitude. Finally, discussion is provided on how these results relate to model testing of partial cavity drag reduction (PCDR) systems for surface ships.


2002 ◽  
Vol 459 ◽  
pp. 67-92 ◽  
Author(s):  
R. A. ANTONIA ◽  
T. ZHOU ◽  
G. P. ROMANO

Measurements have been made in nominally two-dimensional turbulent wakes generated by five different bluff bodies. Each wake has a different level of large-scale organization which is reflected in different amounts of large-scale anisotropy. Structure functions of streamwise (u) and lateral (v) velocity fluctuations at approximately the same value of Rλ, the Taylor microscale Reynolds number, indicate that inertial-range scales are significantly affected by the large-scale anisotropy. The effect is greater on v than u and more pronounced for the porous-body wakes than the solid-body wakes. In particular, ‘relative’ values of the scaling (or power-law) exponents indicate that the magnitude of the transverse exponents can exceed that of the longitudinal ones in the porous-body wakes. This is supported by the inertial-range behaviour of the spectra of u and v. The difference between the transverse and longitudinal exponents appears to depend on the large-scale anisotropy of the flow, as measured by the ratio of the variances of v and u and ratio of the integral length scales of v and u. The spanwise vorticity spectra are much less affected by the anisotropy than the spectra of u and v.


2019 ◽  
Vol 862 ◽  
pp. 552-591 ◽  
Author(s):  
M. Kazakova ◽  
G. L. Richard

We present a new approach to model coastal waves in the shoaling and surf zones. The model can be described as a depth-averaged large-eddy simulation model with a cutoff in the inertial subrange. The large-scale turbulence is explicitly resolved through an extra variable called enstrophy while the small-scale turbulence is modelled with a turbulent-viscosity hypothesis. The equations are derived by averaging the mass, momentum and kinetic energy equations assuming a shallow-water flow, a negligible bottom shear stress and a weakly turbulent flow assumption which is not restrictive in practice. The model is fully nonlinear and has the same dispersive properties as the Green–Naghdi equations. It is validated by numerical tests and by comparison with experimental results of the literature on the propagation of a one-dimensional solitary wave over a mild sloping beach. The wave breaking is characterized by a sudden increase of the enstrophy which allows us to propose a breaking criterion based on the new concept of virtual enstrophy. The model features three empirical parameters. The first one governs the turbulent dissipation and was found to be a constant. The eddy viscosity is determined by a turbulent Reynolds number depending only on the bottom slope. The third parameter defines the breaking criterion and depends only on the wave initial nonlinearity. These dependences give a predictive character to the model which is suitable for further developments.


2019 ◽  
Vol 4 (12) ◽  
Author(s):  
C. Marchioli ◽  
H. Bhatia ◽  
G. Sardina ◽  
L. Brandt ◽  
A. Soldati

Author(s):  
Scott Martin ◽  
Aleksandar Jemcov ◽  
Björn de Ruijter

Here the premixed Conditional Moment Closure (CMC) method is used to model the recent PIV and Raman turbulent, enclosed reacting methane jet data from DLR Stuttgart [1]. The experimental data has a rectangular test section at atmospheric pressure and temperature with a single inlet jet. A jet velocity of 90 m/s is used with an adiabatic flame temperature of 2,064 K. Contours of major species, temperature and velocities along with velocity rms values are provided. The conditional moment closure model has been shown to provide the capability to model turbulent, premixed methane flames with detailed chemistry and reasonable runtimes [2]. The simplified CMC model used here falls into the class of table lookup turbulent combustion models where the chemical kinetics are solved offline over a range of conditions and stored in a table that is accessed by the CFD code. Most table lookup models are based on the laminar 1-D flamelet equations, which assume the small scale turbulence does not affect the reaction rates, only the large scale turbulence has an effect on the reaction rates. The CMC model is derived from first principles to account for the effects of small scale turbulence on the reaction rates, as well as the effects of the large scale mixing, making it more versatile than other models. This is accomplished by conditioning the scalars with the reaction progress variable. By conditioning the scalars and accounting for the small scale mixing, the effects of turbulent fluctuations of the temperature on the reaction rates are more accurately modeled. The scalar dissipation is used to account for the effects of the small scale mixing on the reaction rates. The original premixed CMC model used a constant value of scalar dissipation, here the scalar dissipation is conditioned by the reaction progress variable. The steady RANS 3-D version of the open source CFD code OpenFOAM is used. Velocity, temperature and species are compared to the experimental data. Once validated, this CFD turbulent combustion model will have great utility for designing lean premixed gas turbine combustors.


2014 ◽  
Vol 955-959 ◽  
pp. 2425-2429 ◽  
Author(s):  
Yun Fei Li ◽  
Jian Guo Yang ◽  
Yan Yan Wang ◽  
Xiao Guo Wang

The purpose of this study is to construct a turbulent aggregation device which has specific performance for fine particle aggregation in flue gas. The device consists of two cylindrical pipes and an array of vanes. The pipes extending fully and normal to the gas stream induce large scale turbulence in the form of vortices, while the vanes downstream a certain distance from the pipes induce small one. The process of turbulent aggregation was numerically simulated by coupling the Eulerian multiphase model and population balance model together with a proposed aggregation kernel function taking the size and inertia of particles into account, and based on data of particles’ size distribution measured from the flue of one power plant. The results show that the large scale turbulence generated by pipes favours the aggregation of smaller particles (smaller than 1μm) notably, while the small scale turbulence benefits the aggregation of bigger particles (larger than 1μm) notably and enhances the uniformity of particle size distribution among different particle groups.


2000 ◽  
Vol 663 ◽  
Author(s):  
J. Samper ◽  
R. Juncosa ◽  
V. Navarro ◽  
J. Delgado ◽  
L. Montenegro ◽  
...  

ABSTRACTFEBEX (Full-scale Engineered Barrier EXperiment) is a demonstration and research project dealing with the bentonite engineered barrier designed for sealing and containment of waste in a high level radioactive waste repository (HLWR). It includes two main experiments: an situ full-scale test performed at Grimsel (GTS) and a mock-up test operating since February 1997 at CIEMAT facilities in Madrid (Spain) [1,2,3]. One of the objectives of FEBEX is the development and testing of conceptual and numerical models for the thermal, hydrodynamic, and geochemical (THG) processes expected to take place in engineered clay barriers. A significant improvement in coupled THG modeling of the clay barrier has been achieved both in terms of a better understanding of THG processes and more sophisticated THG computer codes. The ability of these models to reproduce the observed THG patterns in a wide range of THG conditions enhances the confidence in their prediction capabilities. Numerical THG models of heating and hydration experiments performed on small-scale lab cells provide excellent results for temperatures, water inflow and final water content in the cells [3]. Calculated concentrations at the end of the experiments reproduce most of the patterns of measured data. In general, the fit of concentrations of dissolved species is better than that of exchanged cations. These models were later used to simulate the evolution of the large-scale experiments (in situ and mock-up). Some thermo-hydrodynamic hypotheses and bentonite parameters were slightly revised during TH calibration of the mock-up test. The results of the reference model reproduce simultaneously the observed water inflows and bentonite temperatures and relative humidities. Although the model is highly sensitive to one-at-a-time variations in model parameters, the possibility of parameter combinations leading to similar fits cannot be precluded. The TH model of the “in situ” test is based on the same bentonite TH parameters and assumptions as for the “mock-up” test. Granite parameters were slightly modified during the calibration process in order to reproduce the observed thermal and hydrodynamic evolution. The reference model captures properly relative humidities and temperatures in the bentonite [3]. It also reproduces the observed spatial distribution of water pressures and temperatures in the granite. Once calibrated the TH aspects of the model, predictions of the THG evolution of both tests were performed. Data from the dismantling of the in situ test, which is planned for the summer of 2001, will provide a unique opportunity to test and validate current THG models of the EBS.


2016 ◽  
Author(s):  
Janek Meyer ◽  
Hannes Renzsch ◽  
Kai Graf ◽  
Thomas Slawig

While plain vanilla OpenFOAM has strong capabilities with regards to quite a few typical CFD-tasks, some problems actually require additional bespoke solvers and numerics for efficient computation of high-quality results. One of the fields requiring these additions is the computation of large-scale free-surface flows as found e.g. in naval architecture. This holds especially for the flow around typical modern yacht hulls, often planing, sometimes with surface-piercing appendages. Particular challenges include, but are not limited to, breaking waves, sharpness of interface, numerical ventilation (aka streaking) and a wide range of flow phenomenon scales. A new OF-based application including newly implemented discretization schemes, gradient computation and rigid body motion computation is described. In the following the new code will be validated against published experimental data; the effect on accuracy, computational time and solver stability will be shown by comparison to standard OF-solvers (interFoam / interDyMFoam) and Star CCM+. The code’s capabilities to simulate complex “real-world” flows are shown on a well-known racing yacht design.


2002 ◽  
Vol 14 (7) ◽  
pp. 2475 ◽  
Author(s):  
L. Danaila ◽  
F. Anselmet ◽  
R. A. Antonia

Sign in / Sign up

Export Citation Format

Share Document