scholarly journals Vortex-induced vibration of a sphere close to or piercing a free surface

2021 ◽  
Vol 929 ◽  
Author(s):  
Methma M. Rajamuni ◽  
Kerry Hourigan ◽  
Mark C. Thompson
Keyword(s):  
Free Surface ◽  
Response Amplitude ◽  
Surface Boundary ◽  
Vortex Induced Vibration ◽  
Vortical Structures ◽  
Submerged Sphere ◽  
Flow Through ◽  
Lift And Drag ◽  
Experimental Findings ◽  
Submergence Depth

Vortex-induced vibration (VIV) of an elastically mounted sphere placed close to or piercing a free surface (FS) was investigated numerically. The submergence depth ( $h$ ) was systematically varied between $1$ and $-$ 0.75 sphere diameters ( $D$ ) and the response simulated over the reduced velocity range $U^*\in [3.5,14]$ . The incompressible flow was coupled with the sphere motion modelled by a spring–mass–damper system, treating the free-surface boundary as a slip wall. In line with the previous experimental findings, as the submergence depth was decreased from $h^* = h/D =1$ , the maximum response amplitude of the fully submerged sphere decreased; however, as the sphere pierced the FS, the amplitude increased until $h^* = -0.375$ , and then decreased beyond that point. The fluctuating components of the lift and drag coefficients also followed the same pattern. The variation of the near-wake vortex dynamics over this submergence range was examined in detail to understand the effects of $h^*$ on the VIV response. It was found that $h^* = 1$ is a critical submergence depth, beyond which, as $h^*$ is decreased, the vortical structures in the wake vary significantly. For a fully submerged sphere, the influence of the stress-free condition on the VIV response was dominant over the kinematic constraint preventing flow through the surface. For piercing sphere cases, two previously unseen vortical recirculations were formed behind the sphere near times of maximal displacement, enhancing the VIV response. These were strongest at $h^* = -0.375$ , and much weaker for small submergence depths, explaining the observed response-amplitude variation.

Hydrodynamic Forces on a Submerged Sphere Undergoing Large Amplitude Motion

1994 ◽  
Vol 38 (04) ◽  
pp. 272-277
Author(s):  
G. X. Wu
Keyword(s):  
Boundary Condition ◽  
Free Surface ◽  
Large Amplitude ◽  
Multipole Expansion ◽  
The Body ◽  
Surface Boundary ◽  
Free Surface Boundary Condition ◽  
Submerged Sphere ◽  
Instantaneous Position ◽  
Surface Boundary Condition

The hydrodynamic problem of a sphere submerged below a free surface and undergoing large amplitude oscillation is investigated based on the velocity potential theory. The body surface boundary condition is satisfied on its instantaneous position while the free-surface boundary condition is linearized. The solution is obtained by writing the potential in terms of the multipole expansion.

Numerical Solution of Body-Exact Problem in the Time Domain

2013 ◽  
Vol 57 (01) ◽  
pp. 13-23
Author(s):  
Wei Qiu ◽  
Hongxuan (Heather) Peng
Keyword(s):  
Free Surface ◽  
Time Domain ◽  
The Body ◽  
Floating Body ◽  
Entire Body ◽  
Surface Boundary ◽  
Time Step ◽  
Submerged Sphere ◽  
Wigley Hull ◽  
The Time Domain

Motions of a floating body in waves are computed in the time domain by solving the body-exact problem with the panel-free method and exact geometry. In the present study, the body boundary condition is imposed on the instantaneous wetted surface exactly at each time step. The free surface boundary is assumed linear so that the time-domain Green function can be applied. The body geometry is represented by NonUniform Rational B-Spline surfaces. At each time step, the instantaneous wetted surface is obtained by trimming the entire body surface. With the panel-free method, the body-exact problems are solved without involving repanelization of the wetted hull surface at each time step. Validation studies have been carried out for a submerged sphere, a flared body, and a Wigley hull. The hydrodynamic forces on the submerged sphere undergoing large-amplitude motion were computed and compared with analytical solutions. For the flared body oscillating in a free surface and the Wigley hull in waves, numerical results were compared with experimental data and solutions by other numerical methods.

Effects of Submergence on Low and Moderate Reynolds Number Free-Surface Flow Around a Matrix of Cubes

2016 ◽  
Vol 138 (5) ◽  
Author(s):  
Z. Ikram ◽  
E. J. Avital ◽  
J. J. R. Williams
Keyword(s):  
Reynolds Number ◽  
Free Surface ◽  
Free Surface Flow ◽  
Surface Flow ◽  
Immersed Boundary ◽  
Surface Boundary ◽  
Reynolds Number Flow ◽  
Eddy Simulation ◽  
Moderate Reynolds Number ◽  
Submergence Depth

The effect of reducing submergence depth at a low and moderate Reynolds number flow is investigated using large eddy simulation (LES) around a matrix of cubes. The submerged body is modeled using an immersed boundary method, while the free-surface is accounted for using a moving mesh. Results show that for reducing the submergence depth, the forces acting on the cube reduce as the force variation increased. Variation in depth is also found to influence the level of damping and redistribution of turbulence near the free-surface boundary. Both submergence depth and Reynolds number are also found to influence the dominant free-surface signature and shedding frequencies from the cube. In the interobstacle region (IOR), the variation of Reynolds number and submergence depth is found to have little effect.

scholarly journals Experiments on generation of surface waves by an underwater moving bottom

2015 ◽  
Vol 471 (2178) ◽  
pp. 20150069 ◽  
Author(s):  
Timothée Jamin ◽  
Leonardo Gordillo ◽  
Gerardo Ruiz-Chavarría ◽  
Michael Berhanu ◽  
Eric Falcon
Keyword(s):  
Free Surface ◽  
Surface Waves ◽  
Surface Deformation ◽  
Fluid Layer ◽  
Fluid Velocity ◽  
Pass Filter ◽  
Low Pass Filter ◽  
Low Pass ◽  
Spatio Temporal ◽  
Experimental Findings

We report laboratory experiments on surface waves generated in a uniform fluid layer whose bottom undergoes an upward motion. Simultaneous measurements of the free-surface deformation and the fluid velocity field are focused on the role of the bottom kinematics (i.e. its spatio-temporal features) in wave generation. We observe that the fluid layer transfers bottom motion to the free surface as a temporal high-pass filter coupled with a spatial low-pass filter. Both filter effects are often neglected in tsunami warning systems, particularly in real-time forecast. Our results display good agreement with a prevailing linear theory without any parameter fitting. Based on our experimental findings, we provide a simple theoretical approach for modelling the rapid kinematics limit that is applicable even for initially non-flat bottoms: this may be a key step for more realistic varying bathymetry in tsunami scenarios.

Marangoni instability of a thin liquid film heated from below by a local heat source

2003 ◽  
Vol 475 ◽  
pp. 377-408 ◽  
Author(s):  
SERAFIM KALLIADASIS ◽  
ALLA KIYASHKO ◽  
E. A. DEMEKHIN
Keyword(s):  
Free Surface ◽  
Liquid Film ◽  
Discrete Spectrum ◽  
Essential Spectrum ◽  
Marangoni Number ◽  
Thin Liquid Film ◽  
Nonlinear Partial Differential Equations ◽  
Spanwise Direction ◽  
Surface Boundary ◽  
Integral Boundary

We consider the motion of a liquid film falling down a heated planar substrate. Using the integral-boundary-layer approximation of the Navier–Stokes/energy equations and free-surface boundary conditions, it is shown that the problem is governed by two coupled nonlinear partial differential equations for the evolution of the local film height and temperature distribution in time and space. Two-dimensional steady-state solutions of these equations are reported for different values of the governing dimensionless groups. Our computations demonstrate that the free surface develops a bump in the region where the wall temperature gradient is positive. We analyse the linear stability of this bump with respect to disturbances in the spanwise direction. We show that the operator of the linearized system has both a discrete and an essential spectrum. The discrete spectrum bifurcates from resonance poles at certain values of the wavenumber for the disturbances in the transverse direction. The essential spectrum is always stable while part of the discrete spectrum becomes unstable for values of the Marangoni number larger than a critical value. Above this critical Marangoni number the growth rate curve as a function of wavenumber has a finite band of unstable modes which increases as the Marangoni number increases.

Linear and non-linear potential-flow calculations of free-surface waves with lift and induced drag

2000 ◽  
Vol 214 (6) ◽  
pp. 801-812
Author(s):  
C-E Janson
Keyword(s):  
Free Surface ◽  
Potential Flow ◽  
Lift Force ◽  
Radiation Condition ◽  
Linear Potential ◽  
Surface Boundary ◽  
Model Approximation ◽  
Non Linear ◽  
Lifting Surfaces ◽  
The Waves

A potential-flow panel method is used to compute the waves and the lift force from surface-piercing and submerged bodies. In particular the interaction between the waves and the lift produced close to the free surface is studied. Both linear and non-linear free-surface boundary conditions are considered. The potential-flow method is of Rankine-source type using raised source panels on the free surface and a four-point upwind operator to compute the velocity derivatives and to enforce the radiation condition. The lift force is introduced as a dipole distribution on the lifting surfaces and on the trailing wake, together with a flow tangency condition at the trailing edge of the lifting surface. Different approximations for the spanwise circulation distribution at the free surface were tested for a surface-piercing wing and it was concluded that a double-model approximation should be used for low speeds while a single-model, which allows for a vortex at the free surface, was preferred at higher speeds. The lift force and waves from three surface-piercing wings, a hydrofoil and a sailing yacht were computed and compared with measurements and good agreement was obtained.

Scalability studies and large grid computations for surface combatant using CFDShip-Iowa

2011 ◽  
Vol 25 (4) ◽  
pp. 466-487 ◽  
Author(s):  
Shanti Bhushan ◽  
Pablo Carrica ◽  
Jianming Yang ◽  
Frederick Stern
Keyword(s):  
Free Surface ◽  
Message Passing Interface ◽  
Free Surface Flows ◽  
Processing Unit ◽  
Mean Flow ◽  
Time Step ◽  
Vortical Structures ◽  
Cpu Time ◽  
Surface Combatant ◽  
Straight Ahead

Scalability studies and computations using the largest grids to date for free-surface flows are performed using message-passing interface (MPI)-based CFDShip-Iowa toolbox curvilinear (V4) and Cartesian (V6) grid solvers on Navy high-performance computing systems. Both solvers show good strong scalability up to 2048 processors, with V6 showing somewhat better performance than V4. V6 also outperforms V4 in terms of the memory requirements and central processing unit (CPU) time per time-step per grid point. The explicit solvers show better scalability than the implicit solvers, but the latter allows larger time-step sizes, resulting in a lower total CPU time. The multi-grid HYPRE solver shows better scalability than the portable, extensible toolkit for scientific computation solver. The main scalability bottleneck is identified to be the pressure Poisson solver. The memory bandwidth test suggests that further scalability improvements could be obtained by using hybrid MPI/open multi-processing (OpenMP) parallelization. V4-detached eddy simulation (DES) on a 300 M grid for the surface combatant model DTMB 5415 in the straight-ahead condition provides a plausible description of the vortical structures and mean flow patterns observed in the experiments. However, the vortex strengths are over predicted and the turbulence is not resolved. V4-DESs on up to 250 M grids for DTMB 5415 at 20° static drift angle significantly improve the forces and moment predictions compared to the coarse grid unsteady Reynolds averaged Navier–Stokes, due to the improved resolved turbulence predictions. The simulations provide detailed resolution of the free-surface and breaking pattern and vortical and turbulent structures, which will guide planned experiments. V6 simulations on up to 276 M grids for DTMB 5415 in the straight-ahead condition predict diffused vortical structures due to poor wall-layer predictions. This could be due to the limitations of the wall-function implementation for the immersed boundary method.

Effect of an internal nonlinear rotational dissipative element on vortex shedding and vortex-induced vibration of a sprung circular cylinder

2017 ◽  
Vol 828 ◽  
pp. 196-235 ◽  
Author(s):  
Ravi Kumar R. Tumkur ◽  
Arne J. Pearlstein ◽  
Arif Masud ◽  
Oleg V. Gendelman ◽  
Antoine B. Blanchard ◽  
...  
Keyword(s):  
Circular Cylinder ◽  
Nonlinear Energy Sink ◽  
Stream Velocity ◽  
Rectilinear Motion ◽  
Mean Flow ◽  
Number Range ◽  
Vortex Induced Vibration ◽  
Concomitant Reduction ◽  
Fixed Radius ◽  
Lift And Drag

We computationally investigate coupling of a nonlinear rotational dissipative element to a sprung circular cylinder allowed to undergo transverse vortex-induced vibration (VIV) in an incompressible flow. The dissipative element is a ‘nonlinear energy sink’ (NES), consisting of a mass rotating at fixed radius about the cylinder axis and a linear viscous damper that dissipates energy from the motion of the rotating mass. We consider the Reynolds number range $20\leqslant Re\leqslant 120$, with $Re$ based on cylinder diameter and free-stream velocity, and the cylinder restricted to rectilinear motion transverse to the mean flow. Interaction of this NES with the flow is mediated by the cylinder, whose rectilinear motion is mechanically linked to rotational motion of the NES mass through nonlinear inertial coupling. The rotational NES provides significant ‘passive’ suppression of VIV. Beyond suppression however, the rotational NES gives rise to a range of qualitatively new behaviours not found in transverse VIV of a sprung cylinder without an NES, or one with a ‘rectilinear NES’, considered previously. Specifically, the NES can either stabilize or destabilize the steady, symmetric, motionless-cylinder solution and can induce conditions under which suppression of VIV (and concomitant reduction in lift and drag) is accompanied by a greatly elongated region of attached vorticity in the wake, as well as conditions in which the cylinder motion and flow are temporally chaotic at relatively low $Re$.

Fourth‐order finite‐difference P-SV seismograms

Geophysics ◽  
1988 ◽  
Vol 53 (11) ◽  
pp. 1425-1436 ◽  
Author(s):  
Alan R. Levander
Keyword(s):  
Free Surface ◽  
Finite Difference ◽  
Fourth Order ◽  
Equations Of Motion ◽  
Dispersion Analysis ◽  
Numerical Dispersion ◽  
Staggered Grid ◽  
Order System ◽  
Surface Boundary ◽  
Lamb’S Problem

I describe the properties of a fourth‐order accurate space, second‐order accurate time, two‐dimensional P-SV finite‐difference scheme based on the Madariaga‐Virieux staggered‐grid formulation. The numerical scheme is developed from the first‐order system of hyperbolic elastic equations of motion and constitutive laws expressed in particle velocities and stresses. The Madariaga‐Virieux staggered‐grid scheme has the desirable quality that it can correctly model any variation in material properties, including both large and small Poisson’s ratio materials, with minimal numerical dispersion and numerical anisotropy. Dispersion analysis indicates that the shortest wavelengths in the model need to be sampled at 5 gridpoints/wavelength. The scheme can be used to accurately simulate wave propagation in mixed acoustic‐elastic media, making it ideal for modeling marine problems. Explicitly calculating both velocities and stresses makes it relatively simple to initiate a source at the free‐surface or within a layer and to satisfy free‐surface boundary conditions. Benchmark comparisons of finite‐difference and analytical solutions to Lamb’s problem are almost identical, as are comparisons of finite‐difference and reflectivity solutions for elastic‐elastic and acoustic‐elastic layered models.

Sign in / Sign up

Export Citation Format

Share Document