Time-Domain Nonlinear Wave Making Simulation of the Catamaran Based on Velocity Potential Boundary Element Method

2010 ◽  
Author(s):  
Dakui Feng ◽  
Xianzhou Wang ◽  
Zhiguo Zhang ◽  
Yanming Guan
Keyword(s):  
Free Surface ◽  
Radiation Condition ◽  
Flow Fields ◽  
Design Stage ◽  
Hydrodynamic Characteristics ◽  
Outer Side ◽  
Surface Boundary ◽  
Computational Domain ◽  
Time Step ◽  
Surface Boundary Conditions

The catamaran is composed of two monohulls, the flow fields between the inner and outer side of each monohull are different, the bodies must be considered as lifting bodies. So it is very important to know the lifting effect on hydrodynamic characteristics of catamaran hull at the preliminary design stage of its hull form. The pressure Kutta condition is imposed on the trailing-surface of the lifting body by determining the dipole distribution, which generates required circulation on the lifting part. The method is based on Green’s second theorem. Rankine Sources and dipoles are placed on boundary surfaces. Time-stepping scheme is adopted to simulate the wave generated by the catamaran with a uniform speed in deep water. The values of the potential and position of the free surface are updated by integrating the nonlinear Lagrangian free surface boundary conditions for every time. A moving computational window is used in the computations by truncating the fluid domain (the free surface) into a computational domain. The grid regeneration scheme is developed to determine the approximate position of the free surface for the next time step. An implicit implement of far field condition is enforced automatically at the truncation boundary of the computational window, Radiation condition is satisfied automatically. The influences on the wave making resistance of the distance between the twin hulls of the Wigley catamaran on the hydrodynamic characteristics are discussed. The numerical results are presented compared with the existing simulation result. The method can be used to simulate the flow fields around the foil near free surface.

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.

scholarly journals A BEM for the Propagation of Nonlinear Planar Free-surface Waves

2007 ◽  
Vol 5 (1) ◽  
Author(s):  
V. Vinayan ◽  
S. A. Kinnas
Keyword(s):  
Free Surface ◽  
Surface Waves ◽  
Mixed Boundary ◽  
Hydrodynamic Characteristics ◽  
Roll Motion ◽  
Surface Boundary ◽  
Free Surface Waves ◽  
Surface Boundary Conditions ◽  
Stokes Waves ◽  
Fifth Order

A Boundary Element Method (BEM) model for the propagation of non- linear free-surface waves is described and its application to the study of the hydrodynamic characteristics associated with the roll-motion of 2-D hull sec- tions is presented. The roll-motion of the hull section is modeled as a mixed boundary value problem and solved using a higher-order (linear strength dis- tribution) BEM coupled with a Mixed-Eulerian-Lagrangian (MEL) scheme for the time-dependent free-surface boundary conditions. Applications, that in- clude the propagation of fifth-order Stokes waves and waves generated by a piston wave-maker, used to validate the BEM scheme prior to its application to the hull roll-motion are also described.

Time-Domain Simulations of Radiation and Diffraction Forces

2010 ◽  
Vol 54 (02) ◽  
pp. 79-94 ◽  
Author(s):  
Xinshu Zhang ◽  
Piotr Bandyk ◽  
Robert F. Beck
Keyword(s):  
Free Surface ◽  
Time Domain ◽  
Three Dimensional ◽  
The Body ◽  
Surface Boundary ◽  
Two Dimensional ◽  
Forward Speed ◽  
Time Step ◽  
Surface Boundary Conditions ◽  
Strip Theory

Large-amplitude, time-domain, wave-body interactions are studied in this paper for problems with forward speed. Both two-dimensional strip theory and three-dimensional computation methods are shown and compared by a number of numerical simulations. 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 flat panels on the exact body surface, the boundary integral equations are solved numerically at each time step. The strip theory method implements Radial Basis Functions to approximate the longitudinal derivatives of the velocity potential on the body. 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 wetted body geometry. Extensive results are presented to validate the efficiency of the present methods. These results include the added mass and damping computations for a Wigley III hull and an S-175 hull with forward speed using both two-dimensional and three-dimensional approaches. Exciting forces acting on a Wigley III hull due to regular head seas are obtained and compared using both the fully three-dimensional method and the two-dimensional strip theory. All the computational results are compared with experiments or other numerical solutions.

Three-Dimensional Large Amplitude Body Motions in Waves

2008 ◽  
Vol 130 (4) ◽  
Author(s):  
Xinshu Zhang ◽  
Robert F. Beck
Keyword(s):  
Free Surface ◽  
Time Domain ◽  
Three Dimensional ◽  
Surface Boundary ◽  
Forward Speed ◽  
Time Step ◽  
Surface Boundary Conditions ◽  
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 flat panels on the exact wetted body surface, the boundary integral equations are numerically solved 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 wetted body geometry. The desingularized method applied on the free surface produces nonsingular kernels in the integral equations by moving the fundamental singularities a small distance outside of the fluid domain. Constant-strength flat 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 proceeded 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 to the experiments for both linear computations and body-exact computations.

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.

Enhancing Signal-to-Noise Ratios of High-Frequency Rayleigh Waves Extracted from Ambient Seismic Noises in Topographic Region

2020 ◽  
Vol 110 (2) ◽  
pp. 793-802
Author(s):  
Ping Ping ◽  
Risheng Chu ◽  
Yu Zhang ◽  
Jun Xie
Keyword(s):  
Free Surface ◽  
Correlation Functions ◽  
High Frequency ◽  
Rayleigh Waves ◽  
Vertical Component ◽  
Noise Correlation ◽  
Surface Boundary ◽  
Elastic Wedge ◽  
Signal To Noise ◽  
Surface Boundary Conditions

ABSTRACT High-frequency Rayleigh waves can be extracted from ambient seismic noises through noise correlation functions (NCFs), which provides a useful tool to image shallow structures in topographic regions, for example, landslides. Topography may affect signal-to-noise ratios (SNRs) of extracted Rayleigh waves. It is necessary to investigate the propagation features of Rayleigh waves passing a 3D topography. Based on the incident and scattered waves satisfying the free surface boundary conditions, we first derive the displacement responses of Rayleigh waves across a 3D elastic wedge. The results show that the particle motions of Rayleigh waves are an ellipse whose longer axis is always perpendicular to the topographic free surface. Therefore, the Qg component, perpendicular to the topographic free surface, is a better choice to extract high-frequency Rayleigh waves than the conventional vertical component. To verify the choice, we carry out numerical simulations to extract high-frequency NCFs for a typical 3D massif model. Finally, we apply this approach to extract high-frequency Rayleigh-wave NCFs on the Xishancun landslide in southwestern China. The NCFs obtained using the Qg component have more coherent waveforms and higher SNRs than those using the vertical component. We conclude that the Qg component has advantages in extracting high-frequency Rayleigh waves over the conventional vertical component.

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

scholarly journals Linear pulse structure and signalling in a film flow on an inclined plane

1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES
Keyword(s):  
Free Surface ◽  
Saddle Point ◽  
Convective Instability ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Numerical Results ◽  
Navier Stokes ◽  
Surface Boundary ◽  
Navier Stokes Equations ◽  
Surface Boundary Conditions

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