Investigating the wave field near structures in coastal and offshore engineering is of increasing significance. In the present study, simulation is done of the wave profile and flow field for waves propagating over su...Investigating the wave field near structures in coastal and offshore engineering is of increasing significance. In the present study, simulation is done of the wave profile and flow field for waves propagating over submerged bars, using PLIC-VOF （Pieeewise Linear Interface Construction） to trace the free surface of wave and finite difference method to solve vertical 2D Navier-Stokes （N-S） equations. A comparison of the numerical results for two kinds of submerged bars with the experimental ones shows that the PLIC-VOF model used in this study is effective and can compute the wave field precisely.展开更多
The PLIC/SN method that combines the second-order volume tracking method (PLIC-VOF) with the equation of surface normal (SN) vector was recently proposed (M. Sun, “Volume Tracking of Subgrid Particles,” Internationa...The PLIC/SN method that combines the second-order volume tracking method (PLIC-VOF) with the equation of surface normal (SN) vector was recently proposed (M. Sun, “Volume Tracking of Subgrid Particles,” International Journal for Numerical Methods in Fluids, Vol. 66, No. 12, 2011, pp. 1530-1554). The method is able to track the motion of a subgrid particle, but the accuracy is not as good as expected on high resolution grids for vortical flows. In this paper, a simple unsplit multidimensional advection algorithm is coupled with the equation of SN vector. The advection algorithm is formulated as the finite volume method, so that it can be used readily for both structured and unstructured grids while maintaining the exact mass conservation. The new method improves the accuracy significantly for high resolution grids. In the well-known test of the time-resolved vortex problem of T = 2, the circular interface is resolved with an accuracy better than ever using the equation of SN vector.展开更多
该文采用VOF(Volume-of-Fluid)中的PLIC(Piecewise Linear Interface Calculation)界面重构方法模拟研究了不同放置方式气泡之间的相互作用,重点分析了外围流体黏性及气泡间距对其影响。结果表明,外围流体黏性及气泡间距对气泡的融合及...该文采用VOF(Volume-of-Fluid)中的PLIC(Piecewise Linear Interface Calculation)界面重构方法模拟研究了不同放置方式气泡之间的相互作用,重点分析了外围流体黏性及气泡间距对其影响。结果表明,外围流体黏性及气泡间距对气泡的融合及上升速度都具有不同程度的影响,气泡水平放置时,泡间距较小时,气泡上升速度比单个要大。在气泡开始融合到此过程结束的同时,上升速度出现了先降低后反弹的现象,且随着雷诺数的增加速度降低开始时间滞后。展开更多
本文研究的目的在于,使用分段线性PLIC(Piecewise Linear Interface Construction)方式构造界面,推导三维空间中的VOF(Volume of Fluid)算法PLIC3D,给出详尽的计算公式和过程。然后,用球形流体匀速运动作为分析PLIC3D算法的数值算例,比...本文研究的目的在于,使用分段线性PLIC(Piecewise Linear Interface Construction)方式构造界面,推导三维空间中的VOF(Volume of Fluid)算法PLIC3D,给出详尽的计算公式和过程。然后,用球形流体匀速运动作为分析PLIC3D算法的数值算例,比较了PLIC3D和分段常数PCIC(Piecewise Constant Interface Construction)构造法中DA3D(Donor-Acceptor in 3-Dimension)在三维空间计算界面的差别。结果显示PLIC3D优于DA3D,它能够更好地保持流体在平移运动中的形状,证实了PLIC方式的空间构造精度高于PCIC的结论。最后,利用PLIC3D计算自由面,模拟容器中流体泄漏问题,得到了其液面形状变化和速度场。展开更多
In this study, characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical two-dimensional model are Rey...In this study, characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical two-dimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k - ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propagates over breakwater. When wave crest propagates over breakwater, the anticlockwise vortex may generate. On the contrary, when wave hollow propagates over breakwater, the clockwise vortex may generate. Meanwhile, the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy, turbulent dissipation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.展开更多
In the paper, a weak coupling numerical model is developed for the study of the nonlinear dynamic interaction between water waves and permeable sandy seabed. The wave field solveris based on the VOF （Volume of Fluid...In the paper, a weak coupling numerical model is developed for the study of the nonlinear dynamic interaction between water waves and permeable sandy seabed. The wave field solveris based on the VOF （Volume of Fluid） method for continuity equation and the two-dimensional Reynolds Averaged Navier Stokes （RANS） equations with a k-ε closure. The free surface of cnoidal wave is traced through the PLIC-VOF （P/ecewise Linear/nterface Construction）. Blot＇s equations have been applied to solve the sandy seabed, and the u-p fmite dement formulations are derived by the application of the Galerkin weighted-residual procedure. The continuity of the pressure on the interface between fluid and porous medium domains is considered. Laboratory tests were performed to verify the proposed numerical model, and it is shown that the pore-water pressures and the wave heights computed by the VOF-FEM models are in good agreement with the experimental results. It is found that the proposed model is effective in predicting the seabed-nonlinear wave interaction and is able to handle the wave-breakwater-seabed interaction problem.展开更多
Ship bow wave breaking contains complex flow mechanism,which is very important for ship performance.In this study,a practical numerical simulation scheme for bow wave breaking is proposed and the scheme is applied to ...Ship bow wave breaking contains complex flow mechanism,which is very important for ship performance.In this study,a practical numerical simulation scheme for bow wave breaking is proposed and the scheme is applied to the simulation of bow wave breaking of KCS ship model with Fr=0.26,0.30,0.35,0.40,analyzing the impact of speed on the bow wave breaking.The results indicate that an increase in speed leads to a significant rise in viscous pressure resistance and more pronounced bow wave breaking.Moreover,it is found that the traditional wave height function in OpenFOAM is not suitable for detailed studies of bow wave breaking.This study extracts different free surfaces through top and bottom views to further analyze the free surface overturning,droplet splashing,and cavity entrainment in bow wave breaking.Additionally,the spatial and temporal distribution of cavities at Fr=0.40 is analyzed,revealing that cavity distribution is closely related to vortex structures and exhibits a periodic pulsation characteristic of approximately 12 s.展开更多
油水界面不稳定是围油栏失效的主要原因.该文采用油水二维湍流模型,基于PLIC-VOF(Piecewise Linear Interface Construction-Volume of Fluid)方法模拟油水界面运动,以油水密度比值、水流流速、围油栏吃水深度和油的运动黏度为参数,研...油水界面不稳定是围油栏失效的主要原因.该文采用油水二维湍流模型,基于PLIC-VOF(Piecewise Linear Interface Construction-Volume of Fluid)方法模拟油水界面运动,以油水密度比值、水流流速、围油栏吃水深度和油的运动黏度为参数,研究了围油栏前压力梯度的变化及其对油水界面稳定性影响的规律.研究结果表明围油栏前压力垂直结构直接决定了油水界面稳定性;临界压力梯度作为判断油水界面稳定性的标准在多数情况下都是很有效的,但是由于没有考虑运动黏度的影响,当运动黏度小于0.15cm2s-1时,不能再运用临界压力梯度标准判断油水界面稳定性.展开更多
文摘Investigating the wave field near structures in coastal and offshore engineering is of increasing significance. In the present study, simulation is done of the wave profile and flow field for waves propagating over submerged bars, using PLIC-VOF （Pieeewise Linear Interface Construction） to trace the free surface of wave and finite difference method to solve vertical 2D Navier-Stokes （N-S） equations. A comparison of the numerical results for two kinds of submerged bars with the experimental ones shows that the PLIC-VOF model used in this study is effective and can compute the wave field precisely.
文摘The PLIC/SN method that combines the second-order volume tracking method (PLIC-VOF) with the equation of surface normal (SN) vector was recently proposed (M. Sun, “Volume Tracking of Subgrid Particles,” International Journal for Numerical Methods in Fluids, Vol. 66, No. 12, 2011, pp. 1530-1554). The method is able to track the motion of a subgrid particle, but the accuracy is not as good as expected on high resolution grids for vortical flows. In this paper, a simple unsplit multidimensional advection algorithm is coupled with the equation of SN vector. The advection algorithm is formulated as the finite volume method, so that it can be used readily for both structured and unstructured grids while maintaining the exact mass conservation. The new method improves the accuracy significantly for high resolution grids. In the well-known test of the time-resolved vortex problem of T = 2, the circular interface is resolved with an accuracy better than ever using the equation of SN vector.
文摘该文采用VOF(Volume-of-Fluid)中的PLIC(Piecewise Linear Interface Calculation)界面重构方法模拟研究了不同放置方式气泡之间的相互作用,重点分析了外围流体黏性及气泡间距对其影响。结果表明,外围流体黏性及气泡间距对气泡的融合及上升速度都具有不同程度的影响,气泡水平放置时,泡间距较小时,气泡上升速度比单个要大。在气泡开始融合到此过程结束的同时,上升速度出现了先降低后反弹的现象,且随着雷诺数的增加速度降低开始时间滞后。
文摘本文研究的目的在于,使用分段线性PLIC(Piecewise Linear Interface Construction)方式构造界面,推导三维空间中的VOF(Volume of Fluid)算法PLIC3D,给出详尽的计算公式和过程。然后,用球形流体匀速运动作为分析PLIC3D算法的数值算例,比较了PLIC3D和分段常数PCIC(Piecewise Constant Interface Construction)构造法中DA3D(Donor-Acceptor in 3-Dimension)在三维空间计算界面的差别。结果显示PLIC3D优于DA3D,它能够更好地保持流体在平移运动中的形状,证实了PLIC方式的空间构造精度高于PCIC的结论。最后,利用PLIC3D计算自由面,模拟容器中流体泄漏问题,得到了其液面形状变化和速度场。
基金the National Natural Science Foundation of China(Grant No.10472032,50879030,51179077)the Project Founded by Priority Academic Program Development of Jiangsu Higher Education Institutions
基金The National Natural Science Foundation of China under contract Nos 50979008 and 50909009Program for Hunan Province Key Laboratory of WaterSediment Sciences & Flood Hazard Prevention and Open Research Fund Program of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai University No.2008490911
文摘In this study, characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical two-dimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k - ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propagates over breakwater. When wave crest propagates over breakwater, the anticlockwise vortex may generate. On the contrary, when wave hollow propagates over breakwater, the clockwise vortex may generate. Meanwhile, the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy, turbulent dissipation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.
基金The study was financially supported by the National Natural Science Foundation of China(Grant Nos.10202003 and 50479015)Program for New Century Excellent Talents in University(NCET-05-0710)
文摘In the paper, a weak coupling numerical model is developed for the study of the nonlinear dynamic interaction between water waves and permeable sandy seabed. The wave field solveris based on the VOF （Volume of Fluid） method for continuity equation and the two-dimensional Reynolds Averaged Navier Stokes （RANS） equations with a k-ε closure. The free surface of cnoidal wave is traced through the PLIC-VOF （P/ecewise Linear/nterface Construction）. Blot＇s equations have been applied to solve the sandy seabed, and the u-p fmite dement formulations are derived by the application of the Galerkin weighted-residual procedure. The continuity of the pressure on the interface between fluid and porous medium domains is considered. Laboratory tests were performed to verify the proposed numerical model, and it is shown that the pore-water pressures and the wave heights computed by the VOF-FEM models are in good agreement with the experimental results. It is found that the proposed model is effective in predicting the seabed-nonlinear wave interaction and is able to handle the wave-breakwater-seabed interaction problem.
基金Project supported by the National Natural Science Foundation of China(Grant No.52131102).
文摘Ship bow wave breaking contains complex flow mechanism,which is very important for ship performance.In this study,a practical numerical simulation scheme for bow wave breaking is proposed and the scheme is applied to the simulation of bow wave breaking of KCS ship model with Fr=0.26,0.30,0.35,0.40,analyzing the impact of speed on the bow wave breaking.The results indicate that an increase in speed leads to a significant rise in viscous pressure resistance and more pronounced bow wave breaking.Moreover,it is found that the traditional wave height function in OpenFOAM is not suitable for detailed studies of bow wave breaking.This study extracts different free surfaces through top and bottom views to further analyze the free surface overturning,droplet splashing,and cavity entrainment in bow wave breaking.Additionally,the spatial and temporal distribution of cavities at Fr=0.40 is analyzed,revealing that cavity distribution is closely related to vortex structures and exhibits a periodic pulsation characteristic of approximately 12 s.
文摘油水界面不稳定是围油栏失效的主要原因.该文采用油水二维湍流模型,基于PLIC-VOF(Piecewise Linear Interface Construction-Volume of Fluid)方法模拟油水界面运动,以油水密度比值、水流流速、围油栏吃水深度和油的运动黏度为参数,研究了围油栏前压力梯度的变化及其对油水界面稳定性影响的规律.研究结果表明围油栏前压力垂直结构直接决定了油水界面稳定性;临界压力梯度作为判断油水界面稳定性的标准在多数情况下都是很有效的,但是由于没有考虑运动黏度的影响,当运动黏度小于0.15cm2s-1时,不能再运用临界压力梯度标准判断油水界面稳定性.
