scholarly journals CFD Analysis of Water Solitary Wave Reflection

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. 

2017 ◽  
Vol 24 (4) ◽  
pp. 751-762 ◽  
Author(s):  
Vladimir Maderich ◽  
Kyung Tae Jung ◽  
Kateryna Terletska ◽  
Kyeong Ok Kim

Abstract. The dynamics and energetics of a head-on collision of internal solitary waves (ISWs) with trapped cores propagating in a thin pycnocline were studied numerically within the framework of the Navier–Stokes equations for a stratified fluid. The peculiarity of this collision is that it involves trapped masses of a fluid. The interaction of ISWs differs for three classes of ISWs: (i) weakly non-linear waves without trapped cores, (ii) stable strongly non-linear waves with trapped cores, and (iii) shear unstable strongly non-linear waves. The wave phase shift of the colliding waves with equal amplitude grows as the amplitudes increase for colliding waves of classes (i) and (ii) and remains almost constant for those of class (iii). The excess of the maximum run-up amplitude, normalized by the amplitude of the waves, over the sum of the amplitudes of the equal colliding waves increases almost linearly with increasing amplitude of the interacting waves belonging to classes (i) and (ii); however, it decreases somewhat for those of class (iii). The colliding waves of class (ii) lose fluid trapped by the wave cores when amplitudes normalized by the thickness of the pycnocline are in the range of approximately between 1 and 1.75. The interacting stable waves of higher amplitude capture cores and carry trapped fluid in opposite directions with little mass loss. The collision of locally shear unstable waves of class (iii) is accompanied by the development of instability. The dependence of loss of energy on the wave amplitude is not monotonic. Initially, the energy loss due to the interaction increases as the wave amplitude increases. Then, the energy losses reach a maximum due to the loss of potential energy of the cores upon collision and then start to decrease. With further amplitude growth, collision is accompanied by the development of instability and an increase in the loss of energy. The collision process is modified for waves of different amplitudes because of the exchange of trapped fluid between colliding waves due to the conservation of momentum.


Author(s):  
Huan-huan Wang ◽  
Xian-long Jin

The arbitrary Lagrangian–Eulerian method is used to simulate solitary wave interacting with wave wall of a rubble mound breakwater. The Navier–Stokes equations in arbitrary Lagrangian–Eulerian description are solved using operator splitting technique. The penalty-based coupling method is introduced to analyze the coupling of structure and fluid. A piston-type wave maker generating the incident solitary wave is set up in the computational domain. In order to evaluate the numerical model’s performance, a set of experimental studies are carried out in a wave flume using solitary waves at a 1:25 scale. The accuracy of simulation of the solitary wave is verified by comparing with the theoretical value and testing data of the wave surface. Wave transformation, impact and overtopping on the breakwater are simulated in this numerical flume, and the time history of pressure on the wave wall is analyzed. The distributions of maximum impact pressure for solitary waves with different wave heights from the numerical simulations show a good agreement with the experimental results. Based on the dynamic response of the breakwater, the stress distribution and the deformation of the wave wall are discussed. The numerical model can be used as a complementary tool for the design of this kind of structures.


2009 ◽  
Vol 12 (18) ◽  
pp. 59-68
Author(s):  
Thao Danh Nguyen ◽  
Duy The Nguyen

This paper applies and develops a numerical model based on the two-dimensional vertical Navier-Stokes equations to simulate the temporal and spatial variations of wave parameters in front of vertical walls. A non-uniform grids system is performed in the numerical solution of the model by transforming a variable physical domain to a fixed computational domain. Through present model, beside some basic hydrodynamic problems of water waves such as wave profile and water particle velocities, standing wave pressures at the wall are examined. Numerical results of the present model are compared with laboratory data and with existing empirical and theoretical models. The comparisons show that the model can simulate reasonably the wave processes of the waves in front of vertical walls as well as the wave pressures on the wall.


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.


2001 ◽  
Vol 11 (04) ◽  
pp. 565-599 ◽  
Author(s):  
CRISTIAN A. COCLICI ◽  
WOLFGANG L. WENDLAND

We analyze a nonoverlapping domain decomposition method for the treatment of two-dimensional compressible viscous flows around airfoils. Since at some distance to the given profile the inertial forces are strongly dominant, there the viscosity effects are neglected and the flow is assumed to be inviscid. Accordingly, we consider a decomposition of the original flow field into a bounded computational domain (near field) and a complementary outer region (far field). The compressible Navier–Stokes equations are used close to the profile and are coupled with the linearized Euler equations in the far field by appropriate transmission conditions, according to the physical properties and the mathematical type of the corresponding partial differential equations. We present some results of flow around the NACA0012 airfoil and develop an a posteriori analysis of the approximate solution, showing that conservation of mass, momentum and energy are asymptotically attained with the linear model in the far field.


2021 ◽  
Vol 17 ◽  
Author(s):  
B. Kanimozhi ◽  
M. Muthtamilselvan ◽  
Qasem M. Al-Mdallal ◽  
Bahaaeldin Abdalla

Background: This article numerically examines the effect of buoyancy and Marangoni convection in a porous enclosure formed by two concentric cylinders filled with Ag-MgO water hybrid nanofluid. The inner wall of the cavity is maintained at a hot temperature and the outer vertical wall is considered to be cold. The adiabatic condition is assumed for other two boundaries. The effect of magnetic field is considered in radial and axial directions. The Brinkman-extended Darcy model has been adopted in the governing equations. Methods: The finite difference scheme is employed to work out the governing Navier-Stokes equations. The numerically simulated outputs are deliberated in terms of isotherms, streamlines, velocityand average Nusselt number profiles for numerous governing parameters. Results: Except for a greater magnitude of axial magnetic field, our results suggest that the rate of thermal transport accelerates as the nanoparticle volume fraction grows.Also, it is observed that there is an escalation in the profile of average Nusselt numberwith an enhancement in Marangoni number. Conclusion: Furthermore, the suppression of heat and fluid flow in the tall annulus is mainly due to the radial magnetic field whereas in shallow annulus, the axial magnetic field profoundly affects the flow field and thermal transfer.


Author(s):  
Wolfgang Höhn

During the design of the compressor and turbine stages of today’s aeroengines, aerodynamically induced vibrations become increasingly important since higher blade load and better efficiency are desired. In this paper the development of a method based on the unsteady, compressible Navier-Stokes equations in two dimensions is described in order to study the physics of flutter for unsteady viscous flow around cascaded vibrating blades at stall. The governing equations are solved by a finite difference technique in boundary fitted coordinates. The numerical scheme uses the Advection Upstream Splitting Method to discretize the convective terms and central differences discretizing the viscous terms of the fully non-linear Navier-Stokes equations on a moving H-type mesh. The unsteady governing equations are explicitly and implicitly marched in time in a time-accurate way using a four stage Runge-Kutta scheme on a parallel computer or an implicit scheme of the Beam-Warming type on a single processor. Turbulence is modelled using the Baldwin-Lomax turbulence model. The blade flutter phenomenon is simulated by imposing a harmonic motion on the blade, which consists of harmonic body translation in two directions and a rotation, allowing an interblade phase angle between neighboring blades. Non-reflecting boundary conditions are used for the unsteady analysis at inlet and outlet of the computational domain. The computations are performed on multiple blade passages in order to account for nonlinear effects. A subsonic massively stalled unsteady flow case in a compressor cascade is studied. The results, compared with experiments and the predictions of other researchers, show reasonable agreement for inviscid and viscous flow cases for the investigated flow situations with respect to the Steady and unsteady pressure distribution on the blade in separated flow areas as well as the aeroelastic damping. The results show the applicability of the scheme for stalled flow around cascaded blades. As expected the viscous and inviscid computations show different results in regions where viscous effects are important, i.e. in separated flow areas. In particular, different predictions for inviscid and viscous flow for the aerodynamic damping for the investigated flow cases are found.


2016 ◽  
pp. 92-97
Author(s):  
R. E. Volkov ◽  
A. G. Obukhov

The rectangular parallelepiped explicit difference schemes for the numerical solution of the complete built system of Navier-Stokes equations. These solutions describe the three-dimensional flow of a compressible viscous heat-conducting gas in a rising swirling flows, provided the forces of gravity and Coriolis. This assumes constancy of the coefficient of viscosity and thermal conductivity. The initial conditions are the features that are the exact analytical solution of the complete Navier-Stokes equations. Propose specific boundary conditions under which the upward flow of gas is modeled by blowing through the square hole in the upper surface of the computational domain. A variant of parallelization algorithm for calculating gas dynamic and energy characteristics. The results of calculations of gasdynamic parameters dependency on the speed of the vertical blowing by the time the flow of a steady state flow.


2014 ◽  
Vol 16 (5) ◽  
pp. 901-918 ◽  

<div> <p>Three-dimensional calculations were performed to simulate the flow around a cylindrical vegetation element using the Scale Adaptive Simulation (SAS) model; commonly, this is the first step of the modeling of the flow through multiple vegetation elements. SAS solves the Reynolds Averaged Navier-Stokes equations in stable flow regions, while in regions with unstable flow it goes unsteady producing a resolved turbulent spectrum after reducing eddy viscosity according to the locally resolved vortex size represented by the von Karman length scale. A finite volume numerical code was used for the spatial discretisation of the rectangular computational domain with stream-wise, cross-flow and vertical dimensions equal to 30D, 11D and 1D, respectively, which was resolved with unstructured grids. Calculations were compared with experiments and Large Eddy Simulations (LES). Predicted overall flow parameters and mean flow velocities exhibited a very satisfactory agreement with experiments and LES, while the agreement of predicted turbulent stresses was satisfactory. Calculations showed that SAS is an efficient and relatively fast turbulence modeling approach, especially in relevant practical problems, in which the very high accuracy that can be achieved by LES at the expense of large computational times is not required.</p> </div> <p>&nbsp;</p>


Sign in / Sign up

Export Citation Format

Share Document