Numerical study of breaking waves by a two‐phase flow model

A two‐phase flow model, which solves the flow in the air and water simultaneously, is presented for modelling breaking waves in deep and shallow water, including wave pre‐breaking, overturning and post‐breaking processes. The model is based on the Reynolds‐averaged Navier–Stokes equations with the k ?ε turbulence model. The governing equations are solved by the finite volume method in a Cartesian staggered grid and the partial cell treatment is implemented to deal with complex geometries. The SIMPLE algorithm is utilised for the pressure‐velocity coupling and the air‐water interface is modelled by the interface capturing method via a high resolution volume of fluid scheme. The numerical model is validated by simulating overturning waves on a sloping beach and over a reef, and deep‐water breaking waves in a periodic domain, in which good agreement between numerical results and available experimental measurements for the water surface profiles during wave overturning is obtained. The overturning jet, air entrainment and splash‐up during wave breaking have been captured by the two‐phase flow model, which demonstrates the capability of the model to simulate free surface flow and wave breaking problems.Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

A marker-and-cell approach to viscoelastic free surface flows using the PTT model

This work is concerned with the development of a numerical method capable of simulating two-dimensional viscoelastic free surface flows governed by the non-linear constitutive equation PTT (Phan-Thien–Tanner). In particular, we are interested in flows possessing moving free surfaces. The fluid is modelled by a marker-and-cell type method and employs an accurate representation of the fluid surface. Boundary conditions are described in detail and the full free surface stress conditions are considered. The PTT equation is solved by a high order method which requires the calculation of the extra-stress tensor on the mesh contour. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. In order to validate the numerical method fully developed flow in a two-dimensional channel was simulated and the numerical solutions were compared with known analytic solutions. Convergence results were obtained throughout by using mesh refinement. To demonstrate that complex free surface flows using the PTT model can be computed, extrudate swell and a jet flowing onto a rigid plate were simulated.     

用E-L方法计算圆平头物体垂直入水空泡的面闭合和深闭合

本文根据非线性自由面理论,用E-L方法计算了圆平头物体垂直入水空泡发生、发展、面闭合、拉离和深闭合等一系列连续变化过程。空泡中气体流动分为三个阶段加以处理。采用Riabouchinsky映象板模型,初步解决了高F_r数下数值模拟空泡深闭合的困难。计算结果表明,对于给定的物体头型,F_r数、空气出度数、物体质量数是影响入水空泡发展过程的主要因素。     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical studies on non‐linear free surface flow using generalized Schwarz–Christoffel transformation

This paper deals with a technique to transform a free surface flow problem in the physical domain with an unknown boundary to a standard domain that has a fixed boundary. All the difficulties in the physical domain are reduced to finding an unknown mapping function that can be solved iteratively in a standard domain. A derivation is first presented to express an analytic function in terms of the boundary value of its imaginary part. Using a relationship between boundaries of the standard and the physical domains, a formula for the generalized Schwarz–Christoffel transformation is then developed. Based on the generalized Schwarz–Christoffel integral and the Hilbert transform, a pair of non‐linear boundary integro‐differential equations in an infinite strip is formulated for solving fully non‐linear free surface flow problems. The boundary integral equations are then discretized with quadratic elements in an untruncated standard domain and solved by the Levenberg–Marquardt algorithm. Several examples of supercritical flow past obstructions are provided to demonstrate the flexibility and the accuracy of the proposed numerical scheme. Copyright © 2000 John Wiley & Sons, Ltd.     

Direct Numerical Simulation of Flow in a Ribbed Channel

The issue of turbine lifetime is an important one, particularly for modern turbines operating at high temperature regimes. A cooling design such as ribs may achieve an improved lifetime and complex mechanisms of heat transfer need to be well studied. In this paper, a Direct Numerical Simulation (DNS) is presented for a 3-D channel flow with two square ribs on the lower wall. The full unsteady compressible Navier-Stokes equations are solved with an original hybrid finite difference/finite element scheme. The Reynolds number of the simulation is 7 000 based on the bulk velocity at the inlet and the channel height. The present study is mainly devoted to understand the mechanism of heat transfer at the wall through the topological analysis of the flow and the temperature flux. Results show that the large-scale structures generated by obstacles splash onto the lower surface and induce longitudinal vortices which enhance heat transfer at the wall. A comprehensive data base including 56 correlations was set up for testing and improving turbulence models for this complex, separated flow.     

Oscillatory line source for water waves in shear flow

The linearized water-wave radiation problem for the oscillating 2D submerged source in an inviscid shear flow with a free surface is investigated analytically. The vorticity is uniform, with zero velocity at the free surface. Then there will be at most two emitted waves, and no Doppler effects. Exact far-field waves are derived, with radiation conditions applied at infinity. An upstream wave will always exist, whereas the downstream wave exists only when the angular frequency of oscillation exceeds the vorticity. The wave radiation problem is solved also for oscillating vortex and dipoles. The amplitudes and energy fluxes are calculated.     

Symmetry groups and similarity solutions for a free convective boundary-layer problem

The free convective boundary-layer problem due to the motion of an elastic surface into an electrically conducting fluid is studied with group-theoretical methods. The symmetry groups admitted by the corresponding boundary value problem are obtained. Particular attention is paid on the group of scaling which provides the similarity solution of the problem. Also, the admissible form of the data, in order to be conformed to the obtained symmetries, is provided. Finally, with the use of the entailed similarity solution the problem is transformed into a boundary value problem of ODEs and is solved numerically.     

近边界三维水下爆炸气泡动态特性研究

模拟了近壁、近自由面的水下爆炸气泡的非线性动态特性,假定水下爆炸气泡脉动阶段的流场是无旋、不可压缩的,采用高阶曲面三角形单元离散三维气泡表面,用边界积分法求解气泡的运动,在计算奇异积分时通过重新构造双层位势的主值积分消除双层奇异积分的奇异性,得到更精确的结果,并通过合理的加权方法精确的求解边界面上各节点的真实速度,结合弹性网格技术（elastic mesh technique, EMT）得到优化速度,在整个模拟过程中不需要采用数值光顺。将本文的三维模型与轴对称模型进行的对比分析表明,两种模型计算结果吻合很好,并用三维模型模拟了气泡与自由表面及圆筒的相互作用,水下爆炸气泡在自由表面及圆筒的联合作用下呈现出强非线性。     

Mixing layer under a free surface

In the long-wave approximation, the flow of a homogeneous fluid with a free surface in the gravity field is considered. Mathematical models of the surface turbulent layer in shear flows are derived. Steady solutions of the problem of evolution of the mixing layer under the free surface and formation of a surface turbulent jet are constructed. In particular, the problem of the structure of a turbulent bore in a supercritical flow is solved, and the conditions for the formation of a local subcritical zone ahead of the obstacle are studied.     

A net inflow method for incompressible viscous flow with moving free surface

A new finite element procedure called the net inflow method has been developed to simulate time-dependent incompressible viscous flow including moving free surfaces and inertial effects. As a fixed mesh approach with triangular element, the net inflow method can be used to analyse the free surface flow in both regular and irregular domains. Most of the empty elements are excluded from the computational domain, which is adjusted successively to cover the entire region occupied by the liquid. The volume of liquid in a control volume is updated by integrating the net inflow of liquid during each iteration. No additional kinetic equation or material marker needs to be considered. The pressure on the free surface and in the liquid region can be solved explicitly with the continuity equation or implicitly by using the penalty function method. The radial planar free surface flow near a 2D point source and the dam-breaking problem on either a dry bed or a still liquid have been analysed and presented in this paper. The predictions agree very well with available analytical solutions, experimental measurements and/or other numerical results.     

Flow Past a Submerged Hydrofoil

The problem of translational heavy-fluid flow past a hydrofoil with nonlinear boundary conditions is solved by means of a numerical-analytic method based on the singularity distribution over the free surface using the perturbation method. The calculation results are compared with the experimental and calculated data of other authors.     

有限深水中骑行波的显式Hamilton描述

小尺度波（扰动波）迭加在大尺度波（未受扰动波）上形成的波动一般之为“骑行波”。研究了有限可变深度的理想不可压缩流体中的骑行波的显式Hamliltn表示，考虑了自由面上流体与空气之间的表面张力。采用自由面高度和自由面上速度势构成的Hamilton正则变量表示骑行波的动能密度，并在未受扰动波的自由面上作一阶展开。运用复变函数论方法处理了二维流动。先用保角变换将物理平面上的流动区域变换到复势平面上的无限长带形区域，然后在复势平面上用Fourier变换解出Laplace方程，最后经Fourier逆变换求出了扰动波速度热所满足的积分方程。作为特例考虑了平坦底部的流动，导出了动能密度的显式表达式。这里给出的积分方程可以替代相当难解的Hamilton正则方程。通过求解积分方程可得出agrange密度的显式表达式。本文提出的方法约研究骑行波的Hamilton描述以及波的相互作用问题提供了新的途径，有助于了解海面的小尺度波的精细结构。     

Conjugate mixed convection with surface radiation from a vertical electronic board with multiple discrete heat sources

The problem of combined conduction-mixed convection-surface radiation from a vertical electronic board provided with three identical flush-mounted discrete heat sources is solved numerically. The cooling medium is air that is considered to be radiatively transparent. The governing equations for fluid flow and heat transfer are converted from primitive variable form to stream function-vorticity formulation. The equations, thus obtained, are normalised and then are converted into algebraic form using a finite volume based finite difference method. The resulting algebraic equations are then solved using Gauss–Seidel iterative method. An optimum grid system comprising 151 grids along the board and 111 grids across the board is chosen. The effects of various parameters, such as modified Richardson number, surface emissivity and thermal conductivity on temperature distribution along the board, maximum board temperature and relative contributions of mixed convection and radiation to heat dissipation are studied in detail. Further, the contributions of free and forced convection components of mixed convection to board temperature distribution and peak board temperature are brought out. The exclusive roles played by surface radiation and buoyancy in the present problem are clearly elucidated.     

Finite difference method for 2-D and 3-D nonlinear free surface wave problems

A finite difference method is developed for the numerical modelling of the 2-D and 3-D unsteady potential flow generated by transient disturbances on the free surface, on which the nonlinear boundary conditions are fully satisfied. The unknown function is computed with an iteration scheme processing in a transformed time-invariant space. After the velocity is calculated, the location of the free surface is renewed and so is the value of velocity on it. The boundary-value problem of the governing equation is then solved at the next time step. The present method incorporates the FFT. Consequently, a tri-diagonal equation system is obtained which could be readily solved. The feasibility of this method has been demonstrated by 2-D and 3-D examples corresponding to different initial disturbances. This work is supported by the science foundation of Academia Sinica. The paper had been accepted by the XVIth International Congress of IUTAM, Lyngby, Denmark, August, 1984.     

A level set technique applied to unsteady free surface flows

An unsteady Navier–Stokes solver for incompressible fluid is coupled with a level set approach to describe free surface motions. The two‐phase flow of air and water is approximated by the flow of a single fluid whose properties, such as density and viscosity, change across the interface. The free surface location is captured as the zero level of a distance function convected by the flow field. To validate the numerical procedure, two classical two‐dimensional free surface problems in hydrodynamics, namely the oscillating flow in a tank and the waves generated by the flow over a bottom bump, are studied in non‐breaking conditions, and the results are compared with those obtained with other numerical approaches. To check the capability of the method in dealing with complex free surface configurations, the breaking regime produced by the flow over a high bump is analyzed. The analysis covers the successive stages of the breaking phenomenon: the steep wave evolution, the falling jet, the splash‐up and the air entrainment. In all phases, numerical results qualitatively agree with the experimental observations. Finally, to investigate a flow in which viscous effects are relevant, the numerical scheme is applied to study the wavy flow past a submerged hydrofoil. Copyright © 2001 John Wiley & Sons, Ltd.     

Parallel finite element simulation of mooring forces on floating objects

The coupling between the equations governing the free‐surface flows, the six degrees of freedom non‐linear rigid body dynamics, the linear elasticity equations for mesh‐moving and the cables has resulted in a fluid‐structure interaction technology capable of simulating mooring forces on floating objects. The finite element solution strategy is based on a combination approach derived from fixed‐mesh and moving‐mesh techniques. Here, the free‐surface flow simulations are based on the Navier–Stokes equations written for two incompressible fluids where the impact of one fluid on the other one is extremely small. An interface function with two distinct values is used to locate the position of the free‐surface. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian–Eulerian domain. This allows us to handle the motion of the time dependent geometries. Forces and momentums exerted on the floating object by both water and hawsers are calculated and used to update the position of the floating object in time. In the mesh moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The non‐linear rigid body dynamics equations are coupled with the governing equations of fluid flow and are solved simultaneously to update the position of the floating object. The numerical examples includes a 3D simulation of water waves impacting on a moored floating box and a model boat and simulation of floating object under water constrained with a cable. Copyright © 2003 John Wiley & Sons, Ltd.     

Similarity solution of free convective boundary-layer behavior at a stretching surface

In the present study, free convection and heat transfer behavior of electrically conducting fluid in the boundary layer over a vertical continuously stretching surface is investigated. The effects of free convection, magnetic field, suction/blowing at the surface and the stretching speed of the surface on the flow and heat transfer characteristics are considered. By applying one-parametric group theory to analysis of the problem, a similarity solution is found. The governing equations of continuity, momentum and energy are solved numerically by a fourth-order Runge-Kutta scheme. The numerical results, which are obtained for the flow and heat transfer characteristics, reveal the influences of the parameters. Received on 9 September 1998