Numerical Modeling of Breaking Waves Over a Reef With a Level-Set Based Numerical Wave Tank

2013 ◽  
Author(s):  
Mayilvahanan Alagan Chella ◽  
Hans Bihs ◽  
Arun Kamath ◽  
Michael Muskulus

Wave breaking is a highly unsteady, non-linear and extremely turbulent phenomenon. During the wave breaking process, the energy of the wave system is focused close to the crest of the wave and a spatial spread of wave energy occurs. Thus, the description of such a physical phenomenon is highly complex and it requires a deep insight into the breaking wave process. The accurate assessment of breaking wave kinematics is essential for an accurate prediction of hydrodynamic loads on structures. Besides, the understanding of the transformation of waves propagating over an artificial or natural reef is important concerning the coastal processes. The numerical model used in this study is a two-phase model, which solves the flow problem for air and water simultaneously. The Navier-Stokes equations are solved on uniform Cartesian grids in two dimensions. The complex free surface is captured by the level set method. A staggered grid is used for the computation with the velocities defined at the cell edges and the pressure at the cell centres. This avoids unphysical pressure oscillations that can occur due to the coupling of pressure and velocity in the incompressible Navier-Stokes equations. The Ghost Cell Immersed Boundary Method is employed to handle the boundary conditions for complex boundaries. Turbulence modelling is carried out using the k-ω model. Discretization of the convective terms is performed using the 5th order Weighted Essentially Non-Oscillatory (WENO) scheme. In this study, a two-dimensional numerical wave tank is used to simulate waves propagating over steep slopes and wave dissipation. The main objective of the present study is to investigate the wave breaking process over a submerged reef. This is accomplished by examining the wave profile during wave breaking and the breaker indices. Also, the numerical results are compared to data from physical experiments and the numerical results exhibit reasonable agreement with experimental data.

2016 ◽  
Author(s):  
Hans Bihs ◽  
Mayilvahanan Alagan Chella ◽  
Arun Kamath ◽  
Øivind A. Arnsten

For the stability of offshore structures, such as offshore wind foundations, extreme wave conditions need to be taken into account. Waves from extreme events can become critical from design perspective. In a numerical wave tank, extreme waves can be generated through focussed waves. Here, linear waves are generated from a wave spectrum. The wave crests of the generated waves coincide at a pre-selected location and time. In order to test the generated waves, the time series of the free surface elevation are compared with experimental benchmark cases. The numerically simulated free surface shows good agreement with the measurements from experiments. In further computations, the wave impact of the focussed waves on a vertical circular cylinder is investigated. The focussed wave generation is implemented in the numerical wave tank module of REEF3D, which has been extensively and successfully tested for various wave hydrodynamics and wave-structure interaction problems in particular and for free surface flows in general. The open-source CFD code REEF3D solves the three-dimensional Navier-Stokes equations on a staggered Cartesian grid. Solid boundaries are taken into account with the ghost cell immersed boundary method. For the discretization of the convection terms of the momentum equations, the conservative finite difference version of the fifth-order WENO (weighted essentially non-oscillatory) scheme is used. For temporal treatment, the third-order TVD (total variation diminishing) Runge-Kutta scheme is employed. For the pressure, the projection method is used. The free surface flow is solved as two-phase fluid system. For the interface capturing, the level set method is selected. The level set function can be discretized with high-order differencing schemes. This makes it the appropriate solution for wave propagation problems based on Navier-Stokes solvers, which requires high-order numerical methods to avoid artificial wave damping. The numerical model is fully parallelized based on the domain decomposition, using MPI (message passing interface) for internode communication.


2020 ◽  
Vol 14 (4) ◽  
pp. 7369-7378
Author(s):  
Ky-Quang Pham ◽  
Xuan-Truong Le ◽  
Cong-Truong Dinh

Splitter blades located between stator blades in a single-stage axial compressor were proposed and investigated in this work to find their effects on aerodynamic performance and operating stability. Aerodynamic performance of the compressor was evaluated using three-dimensional Reynolds-averaged Navier-Stokes equations using the k-e turbulence model with a scalable wall function. The numerical results for the typical performance parameters without stator splitter blades were validated in comparison with experimental data. The numerical results of a parametric study using four geometric parameters (chord length, coverage angle, height and position) of the stator splitter blades showed that the operational stability of the single-stage axial compressor enhances remarkably using the stator splitter blades. The splitters were effective in suppressing flow separation in the stator domain of the compressor at near-stall condition which affects considerably the aerodynamic performance of the compressor.


2021 ◽  
Vol 11 (2) ◽  
pp. 527
Author(s):  
Tim Tofan ◽  
Harald Kruggel-Emden ◽  
Vytautas Turla ◽  
Raimondas Jasevičius

The numerical simulation and analysis of the ejection of an ink droplet through a nozzle as well its motion through air until its contact with a surface and taking up of a stable form is performed. The fluid flow is modeled by the incompressible Navier–Stokes equations with added surface tension. The presented model can be solved using either a level set or a phase field method to track the fluid interface. Here, the level set method is used to determinate the interface between ink and air. The presented work concentrates on the demonstration how to check the suitability of ink for inkjet printhead nozzles, for instance, for the use in printers. The results such as velocity, change of size, and volume dependence on time of an ink droplet are presented. Recommendations for the use of specific inks are also given.


Author(s):  
Hans Bihs ◽  
Mayilvahanan Alagan Chella ◽  
Arun Kamath ◽  
Øivind Asgeir Arntsen

For the stability of offshore structures, such as offshore wind foundations, extreme wave conditions need to be taken into account. Waves from extreme events are critical from the design perspective. In a numerical wave tank, extreme waves can be modeled using focused waves. Here, linear waves are generated from a wave spectrum. The wave crests of the generated waves coincide at a preselected location and time. Focused wave generation is implemented in the numerical wave tank module of REEF3D, which has been extensively and successfully tested for various wave hydrodynamics and wave–structure interaction problems in particular and for free surface flows in general. The open-source computational fluid dynamics (CFD) code REEF3D solves the three-dimensional Navier–Stokes equations on a staggered Cartesian grid. Higher order numerical schemes are used for time and spatial discretization. For the interface capturing, the level set method is selected. In order to test the generated waves, the time series of the free surface elevation are compared with experimental benchmark cases. The numerically simulated free surface elevation shows good agreement with experimental data. In further computations, the impact of the focused waves on a vertical circular cylinder is investigated. A breaking focused wave is simulated and the associated kinematics is investigated. Free surface flow features during the interaction of nonbreaking focused waves with a cylinder and during the breaking process of a focused wave are also investigated along with the numerically captured free surface.


2013 ◽  
Vol 397-400 ◽  
pp. 783-788
Author(s):  
Xing Wei Zhang ◽  
Chao Wang ◽  
Hang Liu

This paper investigates the aerodynamic forces of several plunging wing models by means of computational fluid dynamics. A finite volume method was used to solve the two-dimensional unsteady incompressible Navier-Stokes equations. The forces and power efficiency have been calculated and compared between sets of different models. Current work found that the nonsymmetrical moving can enhance the lift and thrust forces. The numerical results also prove that the flexible wing model can be use to improve the efficiency and reduce the input. Additionally, a new conceptual model for flapping wing mechanism with active deformation and adaptive nonsymmetrical driving motion is proposed base on the numerical results.


Author(s):  
Fayçal Hammami ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya

This paper deals with the numerical study of bifurcations in a two-sided lid driven cavity flow. The flow is generated by moving the upper wall to the right while moving the left wall downwards. Numerical simulations are performed by solving the unsteady two dimensional Navier-Stokes equations using the finite volume method and multigrid acceleration. In this problem, the ratio of the height to the width of the cavity are ranged from H/L = 0.25 to 1.5. The code for this cavity is presented using rectangular cavity with the grids 144 × 36, 144 × 72, 144 × 104, 144 × 136, 144 × 176 and 144 × 216. Numerous comparisons with the results available in the literature are given. Very good agreements are found between current numerical results and published numerical results. Various velocity ratios ranged in 0.01≤ α ≤ 0.99 at a fixed aspect ratios (A = 0.5, 0.75, 1.25 and 1.5) were considered. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. The stability analysis depending on the aspect ratio, velocity ratios α and the Reynolds number when transition phenomenon occurs is considered in this paper.


Author(s):  
Rahma Yanti

Penelitian ini bertujuan untuk mengetahui efek dari posisi dari bukaan yang ada di rumah Balai Padang. Metode yang digunakan adalah eksperimental dengan menggunakan bantuan software CFD (computational Fluid Dimension) berdasarkan persamaan Navier-Stoke, menggunakan K-Epsilon RNG. Eksperimen dilakukan dengan validasi hasil pengukuran lapangan This study aims to investigate the effect of position openings in Balai Padang house. The numerical methodology is based on solution of the Navier-Stokes equations, using K-epsilon RNG. Numerical results are validated with available field measurement data.


2014 ◽  
Author(s):  
P. Bigay ◽  
A. Bardin ◽  
G. Oger ◽  
D. Le Touzé

In order to efficiently address complex problems in hydrodynamics, the advances in the development of a new method are presented here. This method aims at finding a good compromise between computational efficiency, accuracy, and easy handling of complex geometries. The chosen method is an Explicit Cartesian Finite Volume method for Hydrodynamics (ECFVH) based on a compressible (hyperbolic) solver, with a ghost-cell method for geometry handling and a Level-set method for the treatment of biphase-flows. The explicit nature of the solver is obtained through a weakly-compressible approach chosen to simulate nearly-incompressible flows. The explicit cell-centered resolution allows for an efficient solving of very large simulations together with a straightforward handling of multi-physics. A characteristic flux method for solving the hyperbolic part of the Navier-Stokes equations is used. The treatment of arbitrary geometries is addressed in the hyperbolic and viscous framework. Viscous effects are computed via a finite difference computation of viscous fluxes and turbulent effects are addressed via a Large-Eddy Simulation method (LES). The Level-Set solver used to handle biphase flows is also presented. The solver is validated on 2-D test cases (flow past a cylinder, 2-D dam break) and future improvements are discussed.


1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES

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.


Sign in / Sign up

Export Citation Format

Share Document