A Continuous Desingularized Source Distribution Method Describing Wave–Body Interactions of a Large Amplitude Oscillatory Body

2015 ◽  
Vol 137 (2) ◽  
Author(s):  
Aichun Feng ◽  
Zhi-Min Chen ◽  
W. G. Price
Keyword(s):  
Free Surface ◽  
Boundary Conditions ◽  
Water Surface ◽  
Numerical Solutions ◽  
Free Water ◽  
Source Distribution ◽  
Surface Boundary ◽  
Distribution Method ◽  
Collocation Points ◽  
Calm Water

A Rankine source method with a continuous desingularized free surface source panel distribution is developed to solve numerically a wave–body interaction problem with nonlinear boundary conditions. A body undergoes forced oscillatory motion in a free water surface and the variation of wetted body surface is captured by a regridding process. Free surface sources are placed in continuous panels, rather than points in isolation, over the calm water surface, with free surface collocation points placed on the calm water surface. Nonlinear kinematic and dynamic free surface boundary conditions along the collocation points on the calm water surface are solved in a time domain simulation based on a Lagrange time dependent formulation. Compared with isolated desingularized source points distribution methods, a significantly reduced number of free surface collocation points with sparse distribution are utilized in the present numerical computation. The numerical scheme of study is shown to be computationally efficient and the accuracy of numerical solutions is compared with traditional numerical methods as well as measurements.

A Comparative Study of Wave Breaking Models in a High-Order Spectral Model

2017 ◽  
Author(s):  
Betsy R. Seiffert ◽  
Guillaume Ducrozet
Keyword(s):  
Free Surface ◽  
Energy Dissipation ◽  
Boundary Conditions ◽  
Wave Breaking ◽  
Water Surface ◽  
Eddy Viscosity ◽  
Breaking Waves ◽  
Surface Boundary ◽  
Diffusion Term ◽  
Surface Boundary Conditions

We examine the implementation of two different wave breaking models into the nonlinear potential flow solver HOS-NWT. HOS-NWT is a computationally efficient, open source code that solves for surface elevation in a numerical wave tank using the High-Order Spectral (HOS) method [1]. The first model is a combination of a kinematic wave breaking onset criteria proposed by Barthelemey, et al. [2] and validated by Saket, et al. [3], and an energy dissipation mechanism proposed by Tian, et al. [4, 5]. The wave breaking onset parameter is based on the ratio of local energy flux velocity to the local crest velocity. Once breaking is initiated, an eddy viscosity parameter is estimated based on the pre-breaking local wave geometry, as described in [4, 5]. This eddy viscosity is then added as a diffusion term to the kinematic and dynamic free surface boundary conditions for the duration of wave breaking. Results implementing this wave breaking mechanism in HOS-NWT have shown that the model can successfully calculate the surface elevation and corresponding frequency spectra, as well as the energy dissipation associated with breaking waves [6–8]. The second model implemented to account for wave breaking in HOS-NWT is based on the method proposed by Chalikov, et al. [9–11]. This model defines wave breaking onset by the curvature of the water surface and defines the wave as broken if it exceeds a certain value. A diffusion term is added to the kinematic and dynamic free surface boundary conditions which dissipates energy based on the local curvature of the water surface, which is consequently not constant in space nor time. Calculations made using the two models are compared with large scale experimental measurements conducted at the Hydrodynamics, Energetics and Atmospheric Environment Lab (LHEEA) at Ecole Centrale de Nantes. Comparison of calculations with measurements suggest that both models are successful at predicting wave breaking onset and energy dissipation. However, the model proposed by Barthelemy, et al. [2] and Tian, et al. [4] can be applied without knowing anything about the breaking waves a priori, whereas the model proposed by Chalikov [9] requires tuning to specific conditions.

Novel free surface boundary conditions for spilling breaking waves

2020 ◽  
Vol 159 ◽  
pp. 103717
Author(s):  
Nikta Iravani ◽  
Peyman Badiei ◽  
Maurizio Brocchini
Keyword(s):  
Free Surface ◽  
Boundary Conditions ◽  
Breaking Waves ◽  
Surface Boundary ◽  
Surface Boundary Conditions

The mechanism of vortex connection at a free surface

1999 ◽  
Vol 384 ◽  
pp. 207-241 ◽  
Author(s):  
CHIONG ZHANG ◽  
LIAN SHEN ◽  
DICK K. P. YUE
Keyword(s):  
Free Surface ◽  
Boundary Conditions ◽  
Vortex Ring ◽  
Stokes Equations ◽  
Incidence Angle ◽  
Surface Boundary ◽  
Viscous Layer ◽  
Vertical Vorticity ◽  
Vorticity Component ◽  
Stress Boundary Conditions

Vortex connections at the surface are fundamental and prominent features in free-surface vortical flows. To understand the detailed mechanism of such connection, we consider, as a canonical problem, the laminar vortex connections at a free surface when an oblique vortex ring impinges upon that surface. We perform numerical simulations of the Navier–Stokes equations with viscous free-surface boundary conditions. It is found that the key to understanding the mechanism of vortex connection at a free surface is the surface layers: a viscous layer resulting from the dynamic zero-stress boundary conditions at the free surface, and a thicker blockage layer which is due to the kinematic boundary condition at the surface. In the blockage layer, the vertical vorticity component increases due to vortex stretching and vortex turning (from the transverse vorticity component). The vertical vorticity is then transported to the free surface through viscous diffusion and vortex stretching in the viscous layer leading to increased surface-normal vorticity. These mechanisms take place at the aft-shoulder regions of the vortex ring. Connection at the free surface is different from that at a free-slip wall owing to the generation of surface secondary vorticity. We study the components of this surface vorticity in detail and find that the presence of a free surface accelerates the connection process. We investigate the connection time scale and its dependence on initial incidence angle, Froude and Reynolds numbers. It is found that a criterion based on the streamline topology provides a precise definition for connection time, and may be preferred over existing definitions, e.g. those based on free-surface elevation or net circulation.

Numerical Simulation of Nonlinear Water Wave Propagation Over Rippled Bed

2007 ◽  
Author(s):  
Gerasimos A. Kolokythas ◽  
Athanassios A. Dimas
Keyword(s):  
Free Surface ◽  
Boundary Conditions ◽  
Water Waves ◽  
Stokes Equations ◽  
Surface Flow ◽  
Surface Boundary ◽  
Computational Domain ◽  
Time Step ◽  
Outflow Boundary ◽  
Rippled Bed

In the present study, numerical simulations of the free-surface flow, developing by the propagation of nonlinear water waves over a rippled bottom, are performed assuming that the corresponding flow is two-dimensional, incompressible and viscous. The simulations are based on the numerical solution of the Navier-Stokes equations subject to the fully-nonlinear free-surface boundary conditions and the suitable bottom, inflow and outflow boundary conditions. The equations are properly transformed so that the computational domain becomes time-independent. For the spatial discretization, a hybrid scheme with finite-differences and Chebyshev polynomials is applied, while a fractional time-step scheme is used for the temporal discretization. A wave absorption zone is placed at the outflow region in order to efficiently minimize reflection of waves by the outflow boundary. The numerical model is validated by comparison to the analytical solution for the laminar, oscillatory, current flow which develops a uniform boundary layer over a horizontal bottom. For the propagation of finite-amplitude waves over a rigid rippled bed, the case with wavelength to water depth ratio λ/d0 = 6 and wave height to wavelength ratio H0/λ = 0.05 is considered. The ripples have parabolic shape, while their dimensions — length and height — are chosen accordingly to fit laboratory and field data. Results indicate that the wall shear stress over the ripples and the form drag forces on the ripples increase with increasing ripple height, while the corresponding friction force is insensitive to this increase. Therefore, the percentage of friction in the total drag force decreases with increasing ripple height.

Three-Dimensional Large Amplitude Body Motions in Waves

2007 ◽  
Author(s):  
Xinshu Zhang ◽  
Robert F. Beck
Keyword(s):  
Free Surface ◽  
Three Dimensional ◽  
Surface Boundary ◽  
Forward Speed ◽  
Time Step ◽  
Surface Boundary Conditions ◽  
Submerged Body ◽  
Calm Water ◽  
Wigley Hull ◽  
Constant Strength

Three-dimensional, time-domain, wave-body interactions are studied in this paper for cases with and without forward speed. In the present approach, an exact body boundary condition and linearized free surface boundary conditions are used. By distributing desingularized sources above the calm water surface and using constant-strength panels on the exact submerged body surface, the boundary integral equations are solved numerically at each time step. Once the fluid velocities on the free surface are computed, the free surface elevation and potential are updated by integrating the free surface boundary conditions. After each time step, the body surface and free surface are regrided due to the instantaneous changing submerged body geometry. The desingularized method applied on the free surface produces non-singular kernels in the integral equations by moving the fundamental singularities a small distance outside of the fluid domain. Constant strength panels are used for bodies with any arbitrary shape. Extensive results are presented to validate the efficiency of the present method. These results include the added mass and damping computations for a hemisphere. The calm water wave resistance for a submerged spheroid and a Wigley hull are also presented. All the computations with forward speed are started from rest and proceed until a steady state is reached. Finally, the time-domain forced motion results for a modified Wigley hull with forward speed are shown and compared with the experiments for both linear computations and body-exact computations.

REEF3D::FNPF: A Flexible Fully Nonlinear Potential Flow Solver

2019 ◽  
Author(s):  
Hans Bihs ◽  
Weizhi Wang ◽  
Tobias Martin ◽  
Arun Kamath
Keyword(s):  
Wave Propagation ◽  
Free Surface ◽  
Boundary Conditions ◽  
Potential Flow ◽  
Surface Boundary ◽  
Flow Solver ◽  
Surface Boundary Conditions ◽  
Nonlinear Potential ◽  
Numerical Wave ◽  
Computational Resources

Abstract In situations where the calculation of ocean wave propagation and impact on offshore structures is required, fast numerical solvers are desired in order to find relevant wave events in a first step. After the identification of the relevant events, Computational Fluid Dynamics (CFD) based Numerical Wave Tanks (NWT) with an interface capturing two-phase flow approach can be used to resolve the complex wave structure interaction, including breaking wave kinematics. CFD models emphasize detail of the hydrodynamic physics, which makes them not the ideal candidate for the event identification due to the large computational resources involved. In the current paper a new numerical wave model is represented that solves the Laplace equation for the flow potential and the nonlinear kinematic and dynamics free surface boundary conditions. This approach requires reduced computational resources compared to CFD based NWTs. In contrast to existing approaches, the resulting fully nonlinear potential flow solver REEF3D::FNPF uses a σ-coordinate grid for the computations. Solid boundaries are incorporated through a ghost cell immersed boundary method. The free surface boundary conditions are discretized using fifth-order WENO finite difference methods and the third-order TVD Runge-Kutta scheme for time stepping. The Laplace equation for the potential is solved with Hypres stabilized bi-conjugated gradient solver preconditioned with geometric multi-grid. REEF3D::FNPF is fully parallelized following the domain decomposition strategy and the MPI communication protocol. The model is successfully tested for wave propagation benchmark cases for shallow water conditions with variable bottom as well as deep water.

Sign in / Sign up

Export Citation Format

Share Document