A High-Order Conservative Semi-Lagrangian Solver for 3D Free Surface Flows with Sediment Transport on Voronoi Meshes

In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.

Vertical 2D Modeling of Free Surface Flow with Hydrodynamic Pressure Using SIMPLE Arithmetic in σ Coordinates

A numerical model for shallow water flow has been developed based on the unsteady Reynolds-averaged Navier-Stokes equations with the hydrodynamic pressure instead of hydrostatic pressure assumption. The equations are transformed into the σ-coordinate system and the eddy viscosity is calculated with the standard k-ε turbulence model. The control volume method is used to discrete the equations, and the boundary conditions at the bed for shallow water models only include vertical diffusion terms expressed with wall functions. And the semi-implicit method for pressure linked equation arithmetic is adopted to solve the equations. The model is applied to the 2D vertical plane flow of a current over two steep-sided trenches for which experiment data are available for comparison and good agreement is obtained. And the model is used to predicting the flow in a channel with a steep-sided submerged breakwater at the bottom, and the streamline is drawn.

Validation of a CIP-based tank for numerical simulation of free surface flows

A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the present paper.The tangent of hyperbola for interface capturing（THINC） scheme is applied for capturing complex free surfaces.The new model is capable of simulating a flow with violently varied free surface.A series of computations are conducted to assess the developed algorithm and its versatility.These tests include the collapse of water column with and without an obstacle,sloshing in a fixed tank,the generation of regular waves in a tank,the generation of extreme waves in a tank.Excellent agreements are obtained when numerical results are compared with available analytical,experimental,and other numerical results.

A Higher-Efficient Non-Hydrostatic Finite Volume Model for Strong Three-Dimensional Free Surface Flows and Sediment Transport

In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional （2-D） grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.

Effect of Internal Flow on Vortex-Induced Vibration of Submarine Free Spanning Pipelines

At present, most researches on the vortex-induced vibration of submarine free spanning pipelines ignore the effect of internal flowing fluid; furthermore, there are no research reports considering the coupling effect of internal and external fluid with the free span. In this paper, combining Iwan's wake oscillator model with the differential equation derived for the dynamic response of submarine free spanning pipelines with inclusion of internal flow, the pipe-fluid coupling equations are developed to investigate the effect of internal flow on the vortex-induced vibration of the free spans. The finite element approximation is implemented to derive the matrix equations of equilibrium. The Newmark method combined with simple iteration is used to solve the system of equations. The results indicate that the internal fluid flow may cause the shift of resonance band to the lower frequency and a slight decrease in the peak value; the effect will be more pronounced with the increase of the span length and can be weakened in the presence of the axial tension.

Computational modeling of free-surface slurry flow problems using particle simulation method

The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed.

Phase transition in evolution of traffic flow with scale-free property

This paper investigates the behaviour of traffic flow in traffic systems with a new model based on the NaSch model and cluster approximation of mean-field theory. The proposed model aims at constructing a mapping relationship between the microcosmic behaviour and the macroscopic property of traffic flow. Results demonstrate that scale-free phenomenon of the evolution network becomes obvious when the density value of traffic flow reaches at the critical point of phase transition from free flow to traffic congestion, and jamming is limited in this scale-free structure.

Effect of Viscosity on Free-Surface Waves in Oseen Flows

Based on the complex dispersion relation for the two-dimensional free-surface waves generated by a moving body in the steady Oseen flows, the effect of viscosity on wavelength and wave amplitude was investigated by means of an asymptotic method and a numerical analysis. A comparison between the asymptotic and numerical analysis for the viscous decay factor demonstrates the validity of the perturbation expansions for the wave profile. The numerical result shows that the wavelength of viscous wave is slightly elongated in comparison with that of inviscid wave.

A Higher-Efficient Three-Dimensional Numerical Model for Small Amplitude Free Surface Flows

A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to solve the Navier-Stokes equations in two major steps. A top-layer pressure method is proposed to minimize the number of vertical layers and subsequently the computational cost. Three classical examples of small amplitude free surface flows are used to demonstrate the capability and efficiency of the model. The satisfactory results demonstrated the capability and efficiency of modelling a range of small amplitude free surface flows with only a small number of vertical layers.

A New Efficient Finite Volume Modeling of Small Amplitude Free Surface Flows with Unstructured Grid

A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.

Model for seawater fouling and effects of temperature,flow velocity and surface free energy on seawater fouling

A kinetic model was proposed to predict the seawater fouling process in the seawater heat exchangers.The new model adopted an expression combining depositional and removal behaviors for seawater fouling based on the Kern–Seaton model.The present model parameters include the integrated kinetic rate of deposition(k d)and the integrated kinetic rate of removal(k r),which have clear physical signi ficance.A seawater-fouling monitoring device was established to validate the model.The experimental data were well fitted to the model,and the parameters were obtained in different conditions.SEM and EDX analyses were performed after the experiments,and the results show that the main components of seawater fouling are magnesium hydroxide and aluminum hydroxide.The effects of surface temperature,flow velocity and surface free energy were assessed by the model and the experimental data.The results indicate that the seawater fouling becomes aggravated as the surface temperature increased in a certain range,and the seawater fouling resistance reduced as the flow velocity of seawater increased.Furthermore,the effect of the surface free energy of metals was analyzed,showing that the lower surface free energy mitigates the seawater fouling accumulation.

Analysis on MHD Stability of Free Surface Jet flow in a Gradient Magnetic Fields

The simplified modeling for analysis on MHD stability of free surface jet flow in a gradient magnetic fields is based on the theoretical and experimental results on channel liquid metal MHD flow, especially, the results of MHD flow velocity distribution in cross-section of channels (rectangular duct and circular pipe), and the expected results from the modeling are well agreed with the recent experimental data obtained. It is the first modeling which can efficiently explain the experimental results of liquid-metal free surface jet flow.

Three dimensional free convection couette flow with transpiration cooling

Free convection flow between two vertical parallel plates with transverse sinusoidal injection of the fluid at the stationary plate and its corresponding removal by constant suction through the plate in uniform motion has been analyzed. Due to this type of injection velocity, the flow becomes three-dimensional. Analytical expressions for the velocity, temperature, skin friction and rate of heat transfer were obtained. The important characteristics of the problem, namely the skin friction and the rate of heat transfer are discussed in detail with the help of graphs.

A Study of Two Dimensional Unsteady MHD Free Convection Flow over a Vertical Plate

In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined.

Effect of Variable Viscosity, Dufour, Soret and Thermal Conductivity on Free Convective Heat and Mass Transfer of Non-Darcian Flow past Porous Flat Surface

The motion of incompressible fluid of a variable fluid viscosity and variable thermal conductivity with thermal radiation, Dufour, Soret with heat and mass transfer over a linearly moving porous vertical semi-infinite plate with suction is investigated. The governing equations are transformed into a system of coupled nonlinear ordinary differential equations using similarity transformations with dimensionless variables and solved numerically using shooting method with Runge-Kutta fourth-order method and Newton-Raphson’s interpolation scheme implemented in MATLAB. The result showed that with increase in Dufour and Soret parameter, fluid velocity increases and temperature increases with increase in variation of Dufour while, temperature decreases with increase in Soret. The effects of variable fluid viscosity, variable thermal conductivity, thermal radiation, Soret, Dufour, Prandtl and Schmidt parameters on the dimensionless velocity, temperature and concentration profiles are shown graphically.

Free Surface Flow Generated by Submerged Twin-cylinders in Forward Motion Using a Fully Nonlinear Method

免费表面流动产生了由在强迫的运动成双柱体在免费表面下面沉没了基于边界元素方法被学习。二个相对地点，也就是水平、垂直，为成双的柱体被检验。在两个盒子中，成双的柱体从休息正在开始并且最终通过一个加速的过程以一样的经常的速度移动。假设液体是 inviscid 并且不可压缩并且无旋的流动，不可分的 Laplace 方程能基于边界元素方法是 discretized。充分非线性的边界条件在未知免费表面和动人的身体表面上满足。免费表面被一种 Lagrangian 技术跟踪。Regriding 并且重新协调被使用，它对数字结果的质量关键。单个圆形的柱体和椭圆的柱体被线性方法和充分非线性的方法为精确性检查计算然后充分非线性的方法分别地在成双的柱体盒子上被进行。产生波浪举起和结果的力量被分析象水深度一样讨论在二柱体之间的差距的影响。当动人的柱体接近对方时，不管分发的类型，这被发现当动人的柱体接近对方时，他们在方向与大绝对值承受水动力学力量打手势。随差距的增加变化的力量的趋势能清楚地从数字分析被看见。当当他们接近对方时，水平地分布式的成双的柱体是相反的时，垂直地分布式的成双的柱体似乎与对方一起吸引。

Dissipative free-surface solver for potential flow around hydrofoil distributed with doublets

A doublet integral equation is formulated for the two-dimensional dissipative potential flow around a hydrofoil submerged below a free-water surface. The free-water surface is assumed to involve energy dissipation, and thus it is the source of damping. A doublet panel method is developed from incorporation of the dissipative Green function approach and the doublet distributions on the hydrofoil surface. Numerical computations are implemented, and the derived numerical results are in good agreement with analytic solutions and experimental measurements.

Transient free convection flow of a viscoelastic fluid over vertical surface

In this paper, the viscoelsatic boundary layer flow and the heat transfer near a vertical isothermal impermeable surface and in a quiescent fluid are examined. The governing equations are formulated and solved numerically using MackCormak＇s technique. The results show excellent agreement with previously published results by a comparision. Representative results for the velocity and temperature profiles, boundary layer thicknesses, Nusselt numbers, and local skin friction coefficients are shown graphically for different values of viscoelsatic parameters. In general, it is found that the velocities increase inside the hydrodynamic boundary layers and the temperatures decrease inside the thermal boundary layers for the viscoelsatic fluid as compared with the Newtonian fluid due to favorable tensile stresses. Consequently, the coefficients of friction and heat transfer enhance for higher viscoelsatic parameters.

Development and Flow Modes of Vertical Taylor-Couette System with Free Surface

We have numerically and experimentally investigated the flow modes of Taylor-Couette system consisting of coaxial two cylinders with vertical axes. The inner cylinder rotates and the outer cylinder and the bottom end of the cylinders remain stationary. The upper top boundary is the free surface of the working liquid between the inner and outer cylinders and it contacts with the air. While this flow appears in fluid machinery and chemical reactors and includes industrial interests, it also contains problems of fluid mechanics, which is about the behavior of the free surface in the rotating field. In this paper, we concretely show the developments of the one cell mode flow and the three cell mode flow at a small aspect ratio. We also represent the bifurcation diagram of the flow at the moderate aspect ratio about 5.5. In the numerical simulation, the flow is rest in the initial state, and the inner cylinder is linearly or suddenly accelerated to attain a flow with a prescribed Reynolds number. When the acceleration of the inner cylinder is high, an imperfect bifurcation occurs and the flows of the secondary modes emerge. At high Reynolds numbers, the flow first has many vortices and then some of the vortices collapse and the final stable flow arises. The loci of the normal five cell mode, the anomalous six cell mode and the secondary seven cell mode are determined.