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.

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 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.

Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.

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.

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.

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.

MHD Effect of Liquid Metal Film Flows as Plasma-Facing Components

Stability of liquid metal film flow under gradient magnetic field is investigated. Three dimensional numerical simulations on magnetohydrodynamics （MHD） effect of free surface film flow were carried out, with emphasis on the film thickness variation and its surface stability. Three different MHD phenomena of film flow were observed in the experiment, namely, retardant, rivulet and flat film flow. From our experiment and numerical simulation it can be concluded that flat film flow is a good choice for plasma-facing components （PFCs）

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.

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.

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.

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.

Kinetic-spectrophotometric determination of free cyanide by stopped-flow reversed flow injection analysis

A stopped-flow reversed flow injection method for the determination of free cyanide is proposed. Pyridine-barbituric acid mixture is injected in the flow system as reagent to form the colour species with cyanide. The flow is stopped when the reagent zone comes in the flow cell, where absorbance-time data are collected at 580nm wavelength. The linear range of the determination is 0.1 -10μg/ml CN-. The sampling rate is 60h-1 and the relative standard deviation is 1.6% (n=16) at 5.0 μg/ml CN-1 level. With satisfactory results, the proposed method was applied to the determination of free cyanide in wastewater without sample pretreatment.

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.

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.

A Meshless and Matrix-Free Approach to Modeling Turbulent Fluid Flow Dedicated to Professor Karl Stark Pister for his 95th birthday

A meshless and matrix-free fluid dynamics solver(SOMA)is introduced that avoids the need for user generated and/or analyzed grids,volumes,and meshes.Incremental building of the approximation avoids creation and inversion of possibly dense block diagonal matrices and significantly reduces user interaction.Validation results are presented from the application of SOMA to subsonic,compressible,and turbulent flow over an adiabatic flat plate.

Radiation Effects on Free Convection MHD Couette Flow Started Exponentially with Variable Wall Temperature in Presence of Heat Generation

Radiation effects on free convection MHD Couette flow started exponentially with variable wall temperature in the presence of heat generation have been studied. The governing equations are solved analytically using the Laplace transform technique. The variations of velocity and fluid temperature are presented graphically. It is observed that the velocity decreases with an increase in either magnetic parameter or radiation parameter or Prandtl number. It is also observed that the velocity increases with an increase in either heat generation parameter or Grashof number or accelerated parameter or time. An increase in either radiation parameter or Prandtl number leads to fall in the fluid temperature. It is seen that the fluid temperature increases with an increase in either heat generation parameter or time. Further, it is seen that the shear stress at the moving plate decreases with an increase in either magnetic parameter or radiation parameter while it increases with an increase in either heat generation parameter or Prandtl number. The rate of heat transfer increases with an increase in either Prandtl number or time whereas it decreases with an increase in heat generation parameter.