A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated wate...A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.展开更多
This is a numerical study of a falling droplet surrounding by air under the electric field modeled with finite volume method by means of CFD.The VOF method has been employed to model the two-phase flow of the present ...This is a numerical study of a falling droplet surrounding by air under the electric field modeled with finite volume method by means of CFD.The VOF method has been employed to model the two-phase flow of the present study.Various capillary numbers are investigated to analyze the effects of electric field intensity on the falling droplet deformation.Also,the effects of electric potential on the heat transfer coefficient have been examined.The obtained results showed that by applying the electric field at a capillary number of 0.2 the droplet tends to retain its primitive shape as time goes by,with a subtle deformation to an oblate form.Intensifying the electric field to a capillary number of 0.8 droplet deformation is almost insignificantwith time progressing;however,further enhancement in capillary number to 2 causes the droplet to deform as a prolate shape and higher values of this number intensify the prolate form deformation of the droplet and result in pinch-off phenomenon.Ultimately,it is showed that as the electric potential augments the heat transfer coefficient increases in which for electric potential values higher than 2400 V the heat transfer coefficient enhances significantly.展开更多
This paper aims to investigate a dam break in a channel with a bend in the presence of several obstacles.To accurately determine the flood zones,it is necessary to take into account many factors such as terrain,reserv...This paper aims to investigate a dam break in a channel with a bend in the presence of several obstacles.To accurately determine the flood zones,it is necessary to take into account many factors such as terrain,reservoir volume.Numerical modeling was used to determine the flood zone.Numerical modeling based on the Navier-Stokes equations with a turbulent k-epsilon RNG model,the Volume of Fluid(VOF)method and the PISO algorithm were used to analyze the flow in a bend channel at an angle of 10 with the obstacles.To verify the numerical model,a test on dam break in the 450 channel was conducted.The simulation results were compared with experimental data and with the numerical data of existing data.Having been convinced of the correctness of the mathematical model,the authors carried out a numerical simulation of the main problem in three versions:without barriers,with one obstacle,with two obstacles.According to the obtained numerical results,it can be noted that irregular landforms held the flow,a decrease in water level and a slower time for water emergence could be seen.Thus,the water flow without an obstacle,with one obstacle and with two obstacles showed 4.2 s,4.4 s and 4.6 s of the time of water appearance,respectively.This time shift can give a certain advantage when conducting various events to evacuate people.展开更多
Sediments in the seabed hold vital clues to the study of marine geology,microbial communities and history of ocean life,and the remote operated vehicle(ROV)mounted tubular sampling is an important way to obtain sedime...Sediments in the seabed hold vital clues to the study of marine geology,microbial communities and history of ocean life,and the remote operated vehicle(ROV)mounted tubular sampling is an important way to obtain sediments.However,sampling in the seabed is a particularly difficult and complicated task due to the difficulty accessing deep water layers.The sampling is affected by the sampler’s structural parameters,operation parameters and the interaction between the sampling tube and sediments,which usually results in low volume and coring rate of sediments obtained.This paper simulated the soft viscous seabed sediments as non-Newtonian Herschel-Bulkley viscoplastic fluids and established a numerical model for the tubular sampling based on the volume of fluid(VOF)method.The influence rules of the sampling tube diameter,drainage area rate,penetration velocity,and sediments dynamic viscosity on coring rate and volume were studied.The results showed that coring volume was negatively correlated with all the parameters except the sampling tube diameter.Furthermore,coring rate decreased with increases in penetration velocity,drainage area rate,and sediments dynamic viscosity.The coring rate first increased and then decreased with increasing of the sampling tube diameter,and the peak value was also influenced by penetration velocity.Then,based on the numerical simulation results,an experimental sampling platform was set up and real-world sampling experiments were conducted.The simulation results tallied with the experimental results,with a maximum absolute error of only 4.6%,which verified that the numerical simulation model accurately reflected real-world sampling.The findings in this paper can provide a theoretical basis for facilitating the optimal design of the geometric structure of the seabed sediments samplers and the parameters in the sampling process.展开更多
A stencil-like volume of fluid （VOF） method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions...A stencil-like volume of fluid （VOF） method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions are possible in 2D cases, and the interface can be reconstructed by only requiring the known local volume fraction information. On the other hand, the fluid-occupying-length is defined on each side of the stencil, through which a unified fluid-occupying volume model and a unified algorithm can be obtained to solve the interface advection equation. The method is suitable for the arbitrary geometry of the grid cell, and is extendible to 3D cases. Typical numerical examples show that the current method can give ＂sharp＂ results for tracking free interface.展开更多
Jetting succeeded by accumulation is the characteristic of the vacuum filling,which is different from the conventional pressure-driven flow.In order to simulate this kind of flow,a three-dimensional theoretical model ...Jetting succeeded by accumulation is the characteristic of the vacuum filling,which is different from the conventional pressure-driven flow.In order to simulate this kind of flow,a three-dimensional theoretical model in terms of incompressible and viscous flow is established,and an iterative method combined with finite element method(FEM)is proposed to solve the flow problem.The Lagranian-VOF method is constructed to trace the jetting and accumulated flow fronts.Based on the proposed model and algorithm,a simulation program is developed to predict the velocity,pressure,temperature,and advancement progress.To validate the model and algorithm,a visual experimental equipment for vacuum filling is designed and construted.The vacuum filling experiments with different viscous materials and negative pressures were conducted and compared with the corresponding simulations.The results show the flow front shape closely depends on the fluid viscosity and less relates to the vacuum pressure.展开更多
CO_(2)capture and storage technology is favorable for the reduction of CO_(2)emissions.In recent years,a great number of research achievements have been obtained on CO_(2)geological storage from nano scale to oil/gas ...CO_(2)capture and storage technology is favorable for the reduction of CO_(2)emissions.In recent years,a great number of research achievements have been obtained on CO_(2)geological storage from nano scale to oil/gas reservoir scale,but most studies only focus on theflow behaviors in single-dimension porous media.Besides,the physical experiment method is influenced by many uncertain factors and consumes a lot of time and cost.In order to deeply understand theflow behaviors in the process of CO_(2)geological storage in microscopic view and increase the volume of CO_(2)geological storage,this paper established 2D and 3D models by using VOF(Volume of Fluid)method which can track the dynamic change of two-phase interface,to numerically simulate supercritical CO_(2)-brine two-phaseflow.Then,the distribution characteristics of CO_(2)clusters and the variation laws of CO_(2)saturation under different wettability,capillary number and viscosity ratio conditions were compared,and the intrinsic mechanisms of CO_(2)storage at pore scale were revealed.And the following research results were obtained.First,with the increase of rock wettability to CO_(2),the sweep range of CO_(2)enlarged,and the disconnection frequency of CO_(2)clusters deceased,and thus the volume of CO_(2)storage increased.Second,with the increase of capillary number,the displacement mode transformed from capillaryfingering to stable displacement,and thus the volume of CO_(2)storage increased.Third,as the viscosity of injected supercritical CO_(2)gradually approached that of brine,theflow resistance between two-phasefluids decreased,promoting the"lubricating effect".As a result,theflow capacity of CO_(2)phase was improved,and thus the volume of CO_(2)storage was increased.Fourth,the influence degrees of wettability,capillary number and viscosity ratio on CO_(2)saturation were different in multi-dimensional porous media models.In conclusion,the CO_(2)-brine two-phaseflow simulation based on VOF method revealed theflow mechanisms in the process of CO_(2)geological storage at pore scale,which is of guiding significance to the development of CCUS technology and provides theoretical guidance and technical support for the study of CO_(2)geological storage in a larger scale.展开更多
The icing of areo-engine inlet components during flight can affect engine operational safety.Conventional hot-air anti-icing systems require a large amount of bleed air,which compromises engine performance.Consequentl...The icing of areo-engine inlet components during flight can affect engine operational safety.Conventional hot-air anti-icing systems require a large amount of bleed air,which compromises engine performance.Consequently,low-energy anti/de-icing methods based on superhydrophobic surfaces have attracted widespread attention.Previous studies have demonstrated that for stationary components,superhydrophobic surfaces can significantly reduce anti-icing energy consumption by altering the flow behavior of runback water.However,for rotating inlet components of aero-engines,the effectiveness of superhydrophobic surfaces and the influence of surface wettability on the evolution of runback water flow remain unclear due to the effects of centrifugal and Coriolis forces.This study establishes a 3D liquid water flow simulation model using the volume of fluid(VOF)method to investigate the effects of rotational speed,airflow velocity,and surface wettability on the runback water flow behavior over the rotating spinner under dynamic rotation conditions.The results show that the rotational effects and surface wettability mutually reinforce one another.Specifically,increasing the rotational speed and contact angle can both enhance the flow velocity of liquid water and accelerate the breakup and rupture of liquid film,leading to the formation of rivulets,droplets,and subsequent detachment from the surface.A theoretical model based on force balance is proposed to describe the evolution of runback water flow,and the analysis reveals that as the rotational speed and contact angle increase,the water film is more likely to break up to form rivulets and beads,and the critical radius for droplet detachment from the surface decreases,making it easier removal from the surface.展开更多
The numerical wave channel has been developed based on the volume of fluid method (VOF) in conjunction with the Navier-Stokes equations. The absorbing wave-maker boundary on the left side of the channel is presented b...The numerical wave channel has been developed based on the volume of fluid method (VOF) in conjunction with the Navier-Stokes equations. The absorbing wave-maker boundary on the left side of the channel is presented by prescribing velocity reference to linear wave-maker theory. The principle of which is that the numerical wave-maker is designed to move in a way that generates the required incident wave and cancels out any reflected wave that reach it at the same time. On the right side of the channel, the open boundary is set to permit incident waves to be transmitted freely. The parametric studies have been carried out at a range of ratios of water depth to wave length d/ L from 0.124 to 0.219, with wave height in the front of paddle/water depth ratio (H0 / d) from 0.1 to 0.3. Wave height, wave pressure distribution along the channel and velocity field are obtained for both open boundary condition and reflective boundary condition at the other end of the channel. For a reflective case, it is shown that the absorbing wave-maker is very effective in canceling out the reflected wave that reaches the numerical paddle and highly repeatable waves can be generated.展开更多
The numerical investigation of regular waves interacting with a submerged horizontal twin-plate breakwater is pre- sented in this paper. A numerical model with an absorbing wave-maker is established based on the VOF m...The numerical investigation of regular waves interacting with a submerged horizontal twin-plate breakwater is pre- sented in this paper. A numerical model with an absorbing wave-maker is established based on the VOF method. The validity of the model is verified by experimental results. Comparisons between the numerical and experimental results show that beth the water surface profiles and the wave-induced pressures can be modeled accurately. Wave deformation over the breakwater, water particle velocities around the breakwater, and the wave-induced pressures on the structure are nu- merically investigated. Tile pressure amplitudes of the fundamental and second harmonies on the model surface are exanl- ined in various water depths. The computed and experimental results have revealed that the higher frequency components are generated at the onshore side of the breakwater. Furthermore, the computed results demonstrate a circulating flow formed at the onshore side of the breakwater.展开更多
A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational r...A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.展开更多
Comprehensive experimental and numerical studies have been undertaken to investigate wave energy dissipation performance and main influencing factors of a lower arc-plate breakwater. The numerical model, which conside...Comprehensive experimental and numerical studies have been undertaken to investigate wave energy dissipation performance and main influencing factors of a lower arc-plate breakwater. The numerical model, which considers nonlinear interactions between waves and the arc-plate breakwater, has been constructed by using the velocity wave- generating method, the volume of fluid (VOF) method and the finite volume method. The results show that the relative width, relative height and relative submergence of the breakwater are three main influencing factors and have significant influence on wave energy dissipation of the lower arc-plate open breakwater. The transmission coefficient is found to decrease with the increasing relative width, and the minimum transmission coefficient is 0.15 when the relative width is 0.45. The reflection coefficient is found to vary slightly with the relative width, and the maximum reflection coefficient is 0.53 when the relative width is 0.45. The transmission and reflection coefficients are shown to increase with the relative wave height for approximately 85% of the experimental tests when the relative width is 0.19 0.45. The transmission coefficients at relative submergences of 0.04, 0.02 and 0 are clearly shown to be greater than those at relative submergences of 0.02 and 0.04, while the reflection coefficient exhibits the opposite relationship. After the wave interacts with the lower arc-plate breakwater, the wave energy is mainly converted into transmission, reflection and dissipation energies. The wave attenuation performance is clearly weakened for waves with greater heights and longer periods.展开更多
The influence of baffle position on liquid sloshing during the braking and turning of a tank truck was studied using a volume of fluid (VOF) model. The forces,their positions and weight distribution during braking and...The influence of baffle position on liquid sloshing during the braking and turning of a tank truck was studied using a volume of fluid (VOF) model. The forces,their positions and weight distribution during braking and the forces and rolling moment during turning were calculated. The reliability of the calculation method was validated by comparisons with experimental results. The results showed that during braking,liquid splashes in the tank and the maximum forces and G (the ratio of weight acting on the front axle to the rear axle) are large when A (the ratio of the arch area above the baffle to the area of cross section)≤0.1. When A≥0.2,as the position of the baffle is lowered,the maximum of Fx (the force in direction x) first decreases then increases,and the maximum of Fy (the force in direction y) and G increase. During turning,liquid splashes in the tank and the maximum forces and M (the rolling moment) are large when D (the ratio of the arch area above the baffle to the area of cross section)≤0.2. When D≥0.3,as the position of the baffle is lowered,the maximums of Fy,Fz (the force in direction z) and M increase.展开更多
The volume of fluid (VOF) method is presented to determine the reflection coefficient of and the total horizontal wave force on perforated caisson breakwaters. The present numerical model is compared with a linear ana...The volume of fluid (VOF) method is presented to determine the reflection coefficient of and the total horizontal wave force on perforated caisson breakwaters. The present numerical model is compared with a linear analytic solution obtained by Sahoo et al. (2000). Also this model is verified with the authors′ laboratory data. It is found that the numerical model is in good agreement with the regression equations obtained from the experimental data. The present numerical method is further discussed to relate porosity, the relative wave absorbing chamber depth, the reflection coefficient of perforated caissons and the total horizontal force on them.展开更多
An enhanced numerical model for simulating two-dimensional incompressible viscous flow with distorted free surface is reported. The numerical simulation is carried out through the CIP （Constrained Interpolation Prof...An enhanced numerical model for simulating two-dimensional incompressible viscous flow with distorted free surface is reported. The numerical simulation is carried out through the CIP （Constrained Interpolation Profile）-based method, which is described in the paper. A more accurate interface capturing scheme, the VOF/WLIC scheme （VOF：Volume-of-Fluid;WLIC：weighed line interface calculation）, is adopted as the interface capturing method. To assess the developed algorithm and its versatility, a selection of test problems are examined, i.e. the square wave propagation, the Zalesak’s rigid body rotation, dam breaking problem with and without obstacles, wave sloshing in an excited wave tank and interaction between extreme waves and a floating body. Excellent agreements are obtained when numerical results are compared with available analytical, experimental, and other numerical results. These examples demonstrate that the use of the VOF/WLIC scheme in the free surface capturing makes better results and also the proposed CIP-based model is capable of predicting the freak wave-related phenomena.展开更多
A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） met...A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） method in the frame of two incompressible and immiscible viscous fluids under the action of gravity, taking into consideration of surface tension effects. A comparison of the numerical results by this method with those by other works indicates the validity of the method. In the frame of inviseid and incompressible fluids without taking into consideration of surface tension effects, the mechanisms of the generation of the liquid jet and the transition from spherical shape to toroidal shape during the bubble or drop deformation, the increase of the ring diameter of the toroidal bubble or drop and the decrease of its cross-section area during its motion, and the effects of the density ratio of the two fluids on the deformation of the bubble or drop are analysed both theoretically and numerically.展开更多
Floating structures are commonly seen in coastal and offshore engineering. They are often subjected to extreme waves and, therefore, their nonlinear dynamic behaviors are of great concern. In this paper, an in-house C...Floating structures are commonly seen in coastal and offshore engineering. They are often subjected to extreme waves and, therefore, their nonlinear dynamic behaviors are of great concern. In this paper, an in-house CFD code is developed to investigate the accurate prediction of nonlinear dynamic behaviors of a two-dimensional（2-D） box-shaped floating structure in focused waves. Computations are performed by an enhanced Constrained Interpolation Profile（CIP）-based Cartesian grid model, in which a more accurate VOF（Volume of Fluid） method, the THINC/SW scheme（THINC： tangent of hyperbola for interface capturing; SW： Slope Weighting）, is used for interface capturing. A focusing wave theory is used for the focused wave generation. The wave component of constant steepness is chosen. Comparisons between predictions and physical measurements show good agreement including body motions and free surface profiles. Although the overall agreement is good, some discrepancies are observed for impact pressure on the superstructure due to water on deck. The effect of grid resolution on the results is checked. With a fine grid, no obvious improvement is seen in the global body motions and impact pressures due to water on deck. It is concluded that highly nonlinear phenomena, such as distorted free surface, large-amplitude body motions, and violent impact flow, have been predicted successfully.展开更多
The flow pattern of supercritical flow in bend channels is complicated due to the shock wave phenomenon, which creates difficulties with regard to research and design of bend channels. Using the spillway of an actual ...The flow pattern of supercritical flow in bend channels is complicated due to the shock wave phenomenon, which creates difficulties with regard to research and design of bend channels. Using the spillway of an actual project as an example, a three-dimensional numerical investigation was conducted to simulate the flow in a steep-slope bend based on the renormalization group(RNG) k-ε turbulence flow model and the volume of fluid(VOF) method. The validity of the numerical simulation was demonstrated by comparison between the results of numerical simulation and physical model tests. An optimal scheme of setting vertical vanes in the bend channel is presented. The results of numerical simulation and physical model tests are in agreement, which demonstrates the effectiveness of optimization of vertical vanes and the validity of the three-dimensional numerical simulation. Water depths along both bend walls were analyzed numerically and theoretically. The formula for calculating supercritical water depth along either bend wall was derived, and the critical condition of flow separation from the inner wall was determined.展开更多
A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume...A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid （VOF） method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force.展开更多
In this paper, the solitary wave deformation along a gentle slope and the impact pressure, on the wall are investigated experimentally and the results are compared with numerical results obtained based on the volume o...In this paper, the solitary wave deformation along a gentle slope and the impact pressure, on the wall are investigated experimentally and the results are compared with numerical results obtained based on the volume of fluid (VOF) method. The topography used in the experiment consists of three segments. The left segment is a 1:4 slope, the middle segment 1 :SO slope and the right segment a horizontal bed. Both the wave heights and breaking points obtained from numerical simulation and experiments are in good agreement. Numerical results give reasonable pressure distributions of breaking waves on the wall.展开更多
基金The National Natural Science Foundation of China under contract No.51279023the Public Science and Technology Research Funds Projects of Ocean under contract No.201205023+1 种基金the Special Funds for Postdoctoral Innovative Projects of Liaoning Province of China under contract No.2011921018the Special Funds for Talent Projects of Dalian Ocean University under contract No.SYYJ2011004
文摘A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.
文摘This is a numerical study of a falling droplet surrounding by air under the electric field modeled with finite volume method by means of CFD.The VOF method has been employed to model the two-phase flow of the present study.Various capillary numbers are investigated to analyze the effects of electric field intensity on the falling droplet deformation.Also,the effects of electric potential on the heat transfer coefficient have been examined.The obtained results showed that by applying the electric field at a capillary number of 0.2 the droplet tends to retain its primitive shape as time goes by,with a subtle deformation to an oblate form.Intensifying the electric field to a capillary number of 0.8 droplet deformation is almost insignificantwith time progressing;however,further enhancement in capillary number to 2 causes the droplet to deform as a prolate shape and higher values of this number intensify the prolate form deformation of the droplet and result in pinch-off phenomenon.Ultimately,it is showed that as the electric potential augments the heat transfer coefficient increases in which for electric potential values higher than 2400 V the heat transfer coefficient enhances significantly.
基金supported by the grant from the Ministry of science and Higher education of the Republic of Kazakhstan(AP23489948).
文摘This paper aims to investigate a dam break in a channel with a bend in the presence of several obstacles.To accurately determine the flood zones,it is necessary to take into account many factors such as terrain,reservoir volume.Numerical modeling was used to determine the flood zone.Numerical modeling based on the Navier-Stokes equations with a turbulent k-epsilon RNG model,the Volume of Fluid(VOF)method and the PISO algorithm were used to analyze the flow in a bend channel at an angle of 10 with the obstacles.To verify the numerical model,a test on dam break in the 450 channel was conducted.The simulation results were compared with experimental data and with the numerical data of existing data.Having been convinced of the correctness of the mathematical model,the authors carried out a numerical simulation of the main problem in three versions:without barriers,with one obstacle,with two obstacles.According to the obtained numerical results,it can be noted that irregular landforms held the flow,a decrease in water level and a slower time for water emergence could be seen.Thus,the water flow without an obstacle,with one obstacle and with two obstacles showed 4.2 s,4.4 s and 4.6 s of the time of water appearance,respectively.This time shift can give a certain advantage when conducting various events to evacuate people.
基金Supported by National Key R&D Program of China(Grant No.2016YFC0300502)the National Natural Science Foundation of China(Grant Nos.51705145 and 517779092)+1 种基金Scientific Research Fund of Hunan Provincial Education Department(Grant No.18B205)Hunan Province Natural Science Foundation(Grant No.2019 JJ50182).
文摘Sediments in the seabed hold vital clues to the study of marine geology,microbial communities and history of ocean life,and the remote operated vehicle(ROV)mounted tubular sampling is an important way to obtain sediments.However,sampling in the seabed is a particularly difficult and complicated task due to the difficulty accessing deep water layers.The sampling is affected by the sampler’s structural parameters,operation parameters and the interaction between the sampling tube and sediments,which usually results in low volume and coring rate of sediments obtained.This paper simulated the soft viscous seabed sediments as non-Newtonian Herschel-Bulkley viscoplastic fluids and established a numerical model for the tubular sampling based on the volume of fluid(VOF)method.The influence rules of the sampling tube diameter,drainage area rate,penetration velocity,and sediments dynamic viscosity on coring rate and volume were studied.The results showed that coring volume was negatively correlated with all the parameters except the sampling tube diameter.Furthermore,coring rate decreased with increases in penetration velocity,drainage area rate,and sediments dynamic viscosity.The coring rate first increased and then decreased with increasing of the sampling tube diameter,and the peak value was also influenced by penetration velocity.Then,based on the numerical simulation results,an experimental sampling platform was set up and real-world sampling experiments were conducted.The simulation results tallied with the experimental results,with a maximum absolute error of only 4.6%,which verified that the numerical simulation model accurately reflected real-world sampling.The findings in this paper can provide a theoretical basis for facilitating the optimal design of the geometric structure of the seabed sediments samplers and the parameters in the sampling process.
基金Project supported by the National Natural Science Foundation of China (No.10672097)Shanghai Leading Academic Discipline Project (No.Y0103)
文摘A stencil-like volume of fluid （VOF） method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions are possible in 2D cases, and the interface can be reconstructed by only requiring the known local volume fraction information. On the other hand, the fluid-occupying-length is defined on each side of the stencil, through which a unified fluid-occupying volume model and a unified algorithm can be obtained to solve the interface advection equation. The method is suitable for the arbitrary geometry of the grid cell, and is extendible to 3D cases. Typical numerical examples show that the current method can give ＂sharp＂ results for tracking free interface.
基金the National Science Foundation of China(No.11672271)Shenzhen Zhaowei Machinery&Electronics CO.,LTD.(No.20210035A and 20210035B)for this research work are gratefully acknowledged.
文摘Jetting succeeded by accumulation is the characteristic of the vacuum filling,which is different from the conventional pressure-driven flow.In order to simulate this kind of flow,a three-dimensional theoretical model in terms of incompressible and viscous flow is established,and an iterative method combined with finite element method(FEM)is proposed to solve the flow problem.The Lagranian-VOF method is constructed to trace the jetting and accumulated flow fronts.Based on the proposed model and algorithm,a simulation program is developed to predict the velocity,pressure,temperature,and advancement progress.To validate the model and algorithm,a visual experimental equipment for vacuum filling is designed and construted.The vacuum filling experiments with different viscous materials and negative pressures were conducted and compared with the corresponding simulations.The results show the flow front shape closely depends on the fluid viscosity and less relates to the vacuum pressure.
基金supported by the National Natural Science Foun-dation of China“Research on basic theory for development of deep and ultra-deep reservoirs”(No.52034010)“Flow control during drilling/production of ultra-deep oil and gas wells”(No.52288101).
文摘CO_(2)capture and storage technology is favorable for the reduction of CO_(2)emissions.In recent years,a great number of research achievements have been obtained on CO_(2)geological storage from nano scale to oil/gas reservoir scale,but most studies only focus on theflow behaviors in single-dimension porous media.Besides,the physical experiment method is influenced by many uncertain factors and consumes a lot of time and cost.In order to deeply understand theflow behaviors in the process of CO_(2)geological storage in microscopic view and increase the volume of CO_(2)geological storage,this paper established 2D and 3D models by using VOF(Volume of Fluid)method which can track the dynamic change of two-phase interface,to numerically simulate supercritical CO_(2)-brine two-phaseflow.Then,the distribution characteristics of CO_(2)clusters and the variation laws of CO_(2)saturation under different wettability,capillary number and viscosity ratio conditions were compared,and the intrinsic mechanisms of CO_(2)storage at pore scale were revealed.And the following research results were obtained.First,with the increase of rock wettability to CO_(2),the sweep range of CO_(2)enlarged,and the disconnection frequency of CO_(2)clusters deceased,and thus the volume of CO_(2)storage increased.Second,with the increase of capillary number,the displacement mode transformed from capillaryfingering to stable displacement,and thus the volume of CO_(2)storage increased.Third,as the viscosity of injected supercritical CO_(2)gradually approached that of brine,theflow resistance between two-phasefluids decreased,promoting the"lubricating effect".As a result,theflow capacity of CO_(2)phase was improved,and thus the volume of CO_(2)storage was increased.Fourth,the influence degrees of wettability,capillary number and viscosity ratio on CO_(2)saturation were different in multi-dimensional porous media models.In conclusion,the CO_(2)-brine two-phaseflow simulation based on VOF method revealed theflow mechanisms in the process of CO_(2)geological storage at pore scale,which is of guiding significance to the development of CCUS technology and provides theoretical guidance and technical support for the study of CO_(2)geological storage in a larger scale.
基金supported by the National Natural Science Foundation of China(No.12172029)the open Resarch Project of Key Laboratory of Icing and Anti/De-icing of CARDC(No.IADL 20230101)。
文摘The icing of areo-engine inlet components during flight can affect engine operational safety.Conventional hot-air anti-icing systems require a large amount of bleed air,which compromises engine performance.Consequently,low-energy anti/de-icing methods based on superhydrophobic surfaces have attracted widespread attention.Previous studies have demonstrated that for stationary components,superhydrophobic surfaces can significantly reduce anti-icing energy consumption by altering the flow behavior of runback water.However,for rotating inlet components of aero-engines,the effectiveness of superhydrophobic surfaces and the influence of surface wettability on the evolution of runback water flow remain unclear due to the effects of centrifugal and Coriolis forces.This study establishes a 3D liquid water flow simulation model using the volume of fluid(VOF)method to investigate the effects of rotational speed,airflow velocity,and surface wettability on the runback water flow behavior over the rotating spinner under dynamic rotation conditions.The results show that the rotational effects and surface wettability mutually reinforce one another.Specifically,increasing the rotational speed and contact angle can both enhance the flow velocity of liquid water and accelerate the breakup and rupture of liquid film,leading to the formation of rivulets,droplets,and subsequent detachment from the surface.A theoretical model based on force balance is proposed to describe the evolution of runback water flow,and the analysis reveals that as the rotational speed and contact angle increase,the water film is more likely to break up to form rivulets and beads,and the critical radius for droplet detachment from the surface decreases,making it easier removal from the surface.
文摘The numerical wave channel has been developed based on the volume of fluid method (VOF) in conjunction with the Navier-Stokes equations. The absorbing wave-maker boundary on the left side of the channel is presented by prescribing velocity reference to linear wave-maker theory. The principle of which is that the numerical wave-maker is designed to move in a way that generates the required incident wave and cancels out any reflected wave that reach it at the same time. On the right side of the channel, the open boundary is set to permit incident waves to be transmitted freely. The parametric studies have been carried out at a range of ratios of water depth to wave length d/ L from 0.124 to 0.219, with wave height in the front of paddle/water depth ratio (H0 / d) from 0.1 to 0.3. Wave height, wave pressure distribution along the channel and velocity field are obtained for both open boundary condition and reflective boundary condition at the other end of the channel. For a reflective case, it is shown that the absorbing wave-maker is very effective in canceling out the reflected wave that reaches the numerical paddle and highly repeatable waves can be generated.
基金supported by the National Natural Science Foundation of China ( Grant Nos .50921001 and 51079025)
文摘The numerical investigation of regular waves interacting with a submerged horizontal twin-plate breakwater is pre- sented in this paper. A numerical model with an absorbing wave-maker is established based on the VOF method. The validity of the model is verified by experimental results. Comparisons between the numerical and experimental results show that beth the water surface profiles and the wave-induced pressures can be modeled accurately. Wave deformation over the breakwater, water particle velocities around the breakwater, and the wave-induced pressures on the structure are nu- merically investigated. Tile pressure amplitudes of the fundamental and second harmonies on the model surface are exanl- ined in various water depths. The computed and experimental results have revealed that the higher frequency components are generated at the onshore side of the breakwater. Furthermore, the computed results demonstrate a circulating flow formed at the onshore side of the breakwater.
基金Trans-Century Training program Fund for the Talent,Ministry of Education of China
文摘A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.
基金supported by the NSFC-Shandong Joint Fund(Grant Nos.U1706220 and U1806227)the National Natural Science Foundation of China(Grant Nos.51709140 and 51879019)the Key Laboratory of Coastal Disasters and Defence of Ministry of Education(Grant No.201703)
文摘Comprehensive experimental and numerical studies have been undertaken to investigate wave energy dissipation performance and main influencing factors of a lower arc-plate breakwater. The numerical model, which considers nonlinear interactions between waves and the arc-plate breakwater, has been constructed by using the velocity wave- generating method, the volume of fluid (VOF) method and the finite volume method. The results show that the relative width, relative height and relative submergence of the breakwater are three main influencing factors and have significant influence on wave energy dissipation of the lower arc-plate open breakwater. The transmission coefficient is found to decrease with the increasing relative width, and the minimum transmission coefficient is 0.15 when the relative width is 0.45. The reflection coefficient is found to vary slightly with the relative width, and the maximum reflection coefficient is 0.53 when the relative width is 0.45. The transmission and reflection coefficients are shown to increase with the relative wave height for approximately 85% of the experimental tests when the relative width is 0.19 0.45. The transmission coefficients at relative submergences of 0.04, 0.02 and 0 are clearly shown to be greater than those at relative submergences of 0.02 and 0.04, while the reflection coefficient exhibits the opposite relationship. After the wave interacts with the lower arc-plate breakwater, the wave energy is mainly converted into transmission, reflection and dissipation energies. The wave attenuation performance is clearly weakened for waves with greater heights and longer periods.
文摘The influence of baffle position on liquid sloshing during the braking and turning of a tank truck was studied using a volume of fluid (VOF) model. The forces,their positions and weight distribution during braking and the forces and rolling moment during turning were calculated. The reliability of the calculation method was validated by comparisons with experimental results. The results showed that during braking,liquid splashes in the tank and the maximum forces and G (the ratio of weight acting on the front axle to the rear axle) are large when A (the ratio of the arch area above the baffle to the area of cross section)≤0.1. When A≥0.2,as the position of the baffle is lowered,the maximum of Fx (the force in direction x) first decreases then increases,and the maximum of Fy (the force in direction y) and G increase. During turning,liquid splashes in the tank and the maximum forces and M (the rolling moment) are large when D (the ratio of the arch area above the baffle to the area of cross section)≤0.2. When D≥0.3,as the position of the baffle is lowered,the maximums of Fy,Fz (the force in direction z) and M increase.
文摘The volume of fluid (VOF) method is presented to determine the reflection coefficient of and the total horizontal wave force on perforated caisson breakwaters. The present numerical model is compared with a linear analytic solution obtained by Sahoo et al. (2000). Also this model is verified with the authors′ laboratory data. It is found that the numerical model is in good agreement with the regression equations obtained from the experimental data. The present numerical method is further discussed to relate porosity, the relative wave absorbing chamber depth, the reflection coefficient of perforated caissons and the total horizontal force on them.
基金financially supported by the National Natural Science Foundation of China(Grant No.51209184)the Fundamental Research Funds for the Central Universities(Grant No.2012QNA4020)+1 种基金the Zhejiang Open Foundation of the Most Important Subjects,the Key Laboratory of Water-Sediment Sciences and Water Disaster Prevention of Hunan Province(Grant No.2013SS03)the Educational Commission of Zhejiang Province of China(Grant No.Y201225713)
文摘An enhanced numerical model for simulating two-dimensional incompressible viscous flow with distorted free surface is reported. The numerical simulation is carried out through the CIP （Constrained Interpolation Profile）-based method, which is described in the paper. A more accurate interface capturing scheme, the VOF/WLIC scheme （VOF：Volume-of-Fluid;WLIC：weighed line interface calculation）, is adopted as the interface capturing method. To assess the developed algorithm and its versatility, a selection of test problems are examined, i.e. the square wave propagation, the Zalesak’s rigid body rotation, dam breaking problem with and without obstacles, wave sloshing in an excited wave tank and interaction between extreme waves and a floating body. Excellent agreements are obtained when numerical results are compared with available analytical, experimental, and other numerical results. These examples demonstrate that the use of the VOF/WLIC scheme in the free surface capturing makes better results and also the proposed CIP-based model is capable of predicting the freak wave-related phenomena.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10672043 and 10272032)
文摘A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） method in the frame of two incompressible and immiscible viscous fluids under the action of gravity, taking into consideration of surface tension effects. A comparison of the numerical results by this method with those by other works indicates the validity of the method. In the frame of inviseid and incompressible fluids without taking into consideration of surface tension effects, the mechanisms of the generation of the liquid jet and the transition from spherical shape to toroidal shape during the bubble or drop deformation, the increase of the ring diameter of the toroidal bubble or drop and the decrease of its cross-section area during its motion, and the effects of the density ratio of the two fluids on the deformation of the bubble or drop are analysed both theoretically and numerically.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51209184 and 51479175)Zhejiang Provincial Natural Science Foundation of China(Grant No.LR16E090002)the Open Foundation of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering(Grant No.2013490211)
文摘Floating structures are commonly seen in coastal and offshore engineering. They are often subjected to extreme waves and, therefore, their nonlinear dynamic behaviors are of great concern. In this paper, an in-house CFD code is developed to investigate the accurate prediction of nonlinear dynamic behaviors of a two-dimensional（2-D） box-shaped floating structure in focused waves. Computations are performed by an enhanced Constrained Interpolation Profile（CIP）-based Cartesian grid model, in which a more accurate VOF（Volume of Fluid） method, the THINC/SW scheme（THINC： tangent of hyperbola for interface capturing; SW： Slope Weighting）, is used for interface capturing. A focusing wave theory is used for the focused wave generation. The wave component of constant steepness is chosen. Comparisons between predictions and physical measurements show good agreement including body motions and free surface profiles. Although the overall agreement is good, some discrepancies are observed for impact pressure on the superstructure due to water on deck. The effect of grid resolution on the results is checked. With a fine grid, no obvious improvement is seen in the global body motions and impact pressures due to water on deck. It is concluded that highly nonlinear phenomena, such as distorted free surface, large-amplitude body motions, and violent impact flow, have been predicted successfully.
文摘The flow pattern of supercritical flow in bend channels is complicated due to the shock wave phenomenon, which creates difficulties with regard to research and design of bend channels. Using the spillway of an actual project as an example, a three-dimensional numerical investigation was conducted to simulate the flow in a steep-slope bend based on the renormalization group(RNG) k-ε turbulence flow model and the volume of fluid(VOF) method. The validity of the numerical simulation was demonstrated by comparison between the results of numerical simulation and physical model tests. An optimal scheme of setting vertical vanes in the bend channel is presented. The results of numerical simulation and physical model tests are in agreement, which demonstrates the effectiveness of optimization of vertical vanes and the validity of the three-dimensional numerical simulation. Water depths along both bend walls were analyzed numerically and theoretically. The formula for calculating supercritical water depth along either bend wall was derived, and the critical condition of flow separation from the inner wall was determined.
基金Supported by the Science and Technology Program on Transportation Construction in Western China Ministry of Communications under Grant No.2004-328-832-51
文摘A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid （VOF） method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force.
文摘In this paper, the solitary wave deformation along a gentle slope and the impact pressure, on the wall are investigated experimentally and the results are compared with numerical results obtained based on the volume of fluid (VOF) method. The topography used in the experiment consists of three segments. The left segment is a 1:4 slope, the middle segment 1 :SO slope and the right segment a horizontal bed. Both the wave heights and breaking points obtained from numerical simulation and experiments are in good agreement. Numerical results give reasonable pressure distributions of breaking waves on the wall.
