A Model for the Simulation of Steady Spilling Breaking Waves

2003 ◽  
Vol 47 (01) ◽  
pp. 13-23
Author(s):  
R. Muscari ◽  
A. Di Mascio

A numerical model for the simulation of two-dimensional spilling breaking waves is described. The model is derived from Cointe and Tulin's theory of steady breakers (Cointe & Tulin 1994), although some important changes have been introduced in order to obtain a stable algorithm when coupled with steady-state Reynolds averaged Navier-Stokes equations (RANSE) solvers. In particular, the shape of the breaker and its relation with the following wave height differ from the original model, and moreover, additional conditions for the tangential stress and the turbulent viscosity are proposed. The model has been implemented in a RANSE code, developed for the study of ship flows, through a modification in the free-surface boundary conditions below the breaker. This yields a simple but effective way to reproduce the breaker influence on the underlying flow. The algorithm was used for the simulation of the flow past a submerged hydrofoil. The numerical results are compared with the experimental data by Duncan (1983).

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.


1989 ◽  
Author(s):  
Francesco Martelli ◽  
Vittorio Michelassi

An implicit procedure based on the artificial compressibility formulation is presented for the numerical solution of the two-dimensional incompressible steady Navier-Stokes equations in the presence of large separated regions. Turbulence effects are accounted for by the Chien low Reynolds number form of the K-ε turbulence model and the Baldwin-Lomax algebraic expression for turbulent viscosity. The governing equations are written in conservative form and implicitly solved in fully coupled form using the approximate factorization technique. Preliminary tests were carried out in a laminar flow regime to check the accuracy and stability of the method in two-dimensional and cylindrical axisymmetric flow configurations. After testing in laminar and turbulent flow regimes and comparing the two turbulence models, the code was successfully applied to an actual gas turbine diffuser at low Mach numbers.


2010 ◽  
Vol 658 ◽  
pp. 33-62 ◽  
Author(s):  
XIN GUO ◽  
LIAN SHEN

Direct numerical simulation is performed for the interaction between a deformable free surface and the homogeneous isotropic turbulent flow underneath. The Navier–Stokes equations subject to fully nonlinear free-surface boundary conditions are simulated by using a pseudospectral method in the horizontal directions and a finite-difference method in the vertical direction. Statistically, steady turbulence is generated by using a linear forcing method in the bulk flow below. Through investigation of cases of different Froude and Weber numbers, the present study focuses on the effect of surface deformation of finite amplitude. It is found that the motion of the free surface is characterized by propagating waves and turbulence-generated surface roughness. Statistics of the turbulence field near the free surface are analysed in detail in terms of fluctuations of velocity, fluctuations of velocity gradients and strain rates and the energy budget for horizontal and vertical turbulent motions. Our results illustrate the effects of surface blockage and vanishing shear stress on the anisotropy of the flow field. Using conditional averaging analysis, it is shown that splats and antisplats play an essential role in energy inter-component exchange and vertical transport.


1998 ◽  
Vol 371 ◽  
pp. 207-232 ◽  
Author(s):  
G. VITTORI ◽  
R. VERZICCO

Numerical simulations of Navier–Stokes equations are performed to study the flow originated by an oscillating pressure gradient close to a wall characterized by small imperfections. The scenario of transition from the laminar to the turbulent regime is investigated and the results are interpreted in the light of existing analytical theories. The ‘disturbed-laminar’ and the ‘intermittently turbulent’ regimes detected experimentally are reproduced by the present simulations. Moreover it is found that imperfections of the wall are of fundamental importance in causing the growth of two-dimensional disturbances which in turn trigger turbulence in the Stokes boundary layer. Finally, in the intermittently turbulent regime, a description is given of the temporal development of turbulence characteristics.


Sign in / Sign up

Export Citation Format

Share Document