Numerical Computations for a Nonlinear Free Surface Problem in Shallow Water

Author(s):  
K. J. Bai ◽  
J. H. Kyoung ◽  
J. W. Kim

This paper describes a finite element method applied to a nonlinear free surface flow problem for a ship moving in three dimensions. The physical model is taken to simulate the towing tank experimental conditions. The exact nonlinear free-surface flow problem formulated by an initial/boundary value problem is replaced by an equivalent weak formulation. The same problem was considered earlier by Bai, et. al. [1] where some irregularities were observed in the downstream waves and a transom stern ship geometry could not be treated. In the present paper, specifically, three improvements are made from the earlier work. The first improvement is the introduction of the 5-point Chebyshev filtering scheme which eliminates the irregular and saw-toothed waves in the downstream. The second improvement is that now we can treat a transom stern ship geometry. The third improvement is the introduction of a new boundary condition to simulate a dry bottom behind a transom stern ship which is stretched from the free surface to the bottom at a high Froude number. Computations are made for two models. The first model is tested for the generation of the solitons in the upstream and smooth waves in the downstream. The second model is used to compute the generation of a dry bottom behind a transom stern which is one of highly nonlinear phenomena. The results of the first model show a good agreement with previous results for the generation of the solitons. The results of the second model also show a good agreement with the preliminary experimental observation for a dry-bottom, which will be reported in near future. The numerical simulation of the second model can be applied to the local flow behind a sail of a submarine in cruise, a sloshing problem in LNG tankers, and a dam breaking problem. Both computed models are assumed to be in shallow water for simplicity. However, the present numerical method can treat arbitrary water-depth and practical ship geometries.

2016 ◽  
Vol 808 ◽  
pp. 441-468 ◽  
Author(s):  
S. L. Gavrilyuk ◽  
V. Yu. Liapidevskii ◽  
A. A. Chesnokov

A two-layer long-wave approximation of the homogeneous Euler equations for a free-surface flow evolving over mild slopes is derived. The upper layer is turbulent and is described by depth-averaged equations for the layer thickness, average fluid velocity and fluid turbulent energy. The lower layer is almost potential and can be described by Serre–Su–Gardner–Green–Naghdi equations (a second-order shallow water approximation with respect to the parameter $H/L$, where $H$ is a characteristic water depth and $L$ is a characteristic wavelength). A simple model for vertical turbulent mixing is proposed governing the interaction between these layers. Stationary supercritical solutions to this model are first constructed, containing, in particular, a local turbulent subcritical zone at the forward slope of the wave. The non-stationary model was then numerically solved and compared with experimental data for the following two problems. The first one is the study of surface waves resulting from the interaction of a uniform free-surface flow with an immobile wall (the water hammer problem with a free surface). These waves are sometimes called ‘Favre waves’ in homage to Henry Favre and his contribution to the study of this phenomenon. When the Froude number is between 1 and approximately 1.3, an undular bore appears. The characteristics of the leading wave in an undular bore are in good agreement with experimental data by Favre (Ondes de Translation dans les Canaux Découverts, 1935, Dunod) and Treske (J. Hydraul Res., vol. 32 (3), 1994, pp. 355–370). When the Froude number is between 1.3 and 1.4, the transition from an undular bore to a breaking (monotone) bore occurs. The shoaling and breaking of a solitary wave propagating in a long channel (300 m) of mild slope (1/60) was then studied. Good agreement with experimental data by Hsiao et al. (Coast. Engng, vol. 55, 2008, pp. 975–988) for the wave profile evolution was found.


2007 ◽  
Vol 24 (12) ◽  
pp. 1699-1722
Author(s):  
Meng-Chi Hung ◽  
Te-Yung Hsieh ◽  
Tung-Lin Tsai ◽  
Jinn-Chuang Yang

2007 ◽  
Vol 51 (01) ◽  
pp. 47-64
Author(s):  
James C. Huan ◽  
Thomas T. Huang

A fast turnaround and an accurate computational fluid dynamics (CFD) approach for ship total resistance prediction is developed. The approach consists of a nonlinear free surface potential flow solver (PShip code) with a wet-or-dry transom stern model, and a Reynolds-averaged Navier-Stokes (RANS) equation solver that solves viscous free surface flow with a prescribed free surface given from the PShip. The prescribed free surface RANS predicts a viscous correction to the pressure resistance (viscous form) and viscous flow field around the hull. The viscous free surface flow solved this way avoids the time-consuming RANS iterations to resolve the free surface profile. The method, however, requires employing a flow characteristic-based nonreflecting boundary condition at the free surface. The approach can predict the components of ship resistance, the associated wave profile around the hull, and the sinkage and trim of the ship. Validation of the approach is presented with Wigley, Series 60 (CB = 0.6), and NSWCCD Model 5415 hulls. An overall accuracy of ±2% for ship total resistance prediction is achieved. The approach is applied to evaluating the effects of a stern flap on a DD 968 model on ship performance. An empirical viscous form resistance formula is also devised for a quick ship total resistance estimate.


2005 ◽  
Vol 127 (3) ◽  
pp. 572-582 ◽  
Author(s):  
Shin Hyung Rhee

The present study is concerned with liquid tank sloshing at low filling level conditions. The volume of fluid method implemented in a Navier–Stokes computational fluid dynamics code is employed to handle the free-surface flow of liquid sloshing. The geometric reconstruction scheme for the interface representation is employed to ensure sharpness at the free surface. The governing equations are discretized by second order accurate schemes on unstructured grids. Several different computational approaches are verified and numerical uncertainties are assessed. The computational results are validated against existing experimental data, showing good agreement. The capability is demonstrated for a generic membrane-type liquefied natural gas carrier tank with a simplified pump tower inside. The validation results suggest that the present computational approach is both easy to apply and accurate enough for more realistic problems.


2017 ◽  
Vol 97 (10) ◽  
pp. 1701-1716
Author(s):  
N. Foukroun ◽  
R. Ait-Yahia-Djouadi ◽  
D. Hernane-Boukari

1998 ◽  
Vol 364 ◽  
pp. 325-347 ◽  
Author(s):  
MING XUE ◽  
DICK K. P. YUE

The unsteady fully nonlinear free-surface flow due to an impulsively started submerged point sink is studied in the context of incompressible potential flow. For a fixed (initial) submergence h of the point sink in otherwise unbounded fluid, the problem is governed by a single non-dimensional physical parameter, the Froude number, [Fscr ]≡Q/4π(gh5)1/2, where Q is the (constant) volume flux rate and g the gravitational acceleration. We assume axisymmetry and perform a numerical study using a mixed-Eulerian–Lagrangian boundary-integral-equation scheme. We conduct systematic simulations varying the parameter [Fscr ] to obtain a complete quantification of the solution of the problem. Depending on [Fscr ], there are three distinct flow regimes: (i) [Fscr ]<[Fscr ]1≈0.1924 – a ‘sub-critical’ regime marked by a damped wave-like behaviour of the free surface which reaches an asymptotic steady state; (ii) [Fscr ]1<[Fscr ]<[Fscr ]2≈0.1930 – the ‘trans-critical’ regime characterized by a reversal of the downward motion of the free surface above the sink, eventually developing into a sharp upward jet; (iii) [Fscr ]>[Fscr ]2 – a ‘super-critical’ regime marked by the cusp-like collapse of the free surface towards the sink. Mechanisms behind such flow behaviour are discussed and hydrodynamic quantities such as pressure, power and force are obtained in each case. This investigation resolves the question of validity of a steady-state assumption for this problem and also shows that a small-time expansion may be inadequate for predicting the eventual behaviour of the flow.


Sign in / Sign up

Export Citation Format

Share Document