This paper investigates the hydrodynamic performance of a cylindrical-dual or rectangular-single pontoon floating breakwater using the numerical method and experimental study. The numerical simulation work is based on...This paper investigates the hydrodynamic performance of a cylindrical-dual or rectangular-single pontoon floating breakwater using the numerical method and experimental study. The numerical simulation work is based on the multi-physics computational fluid dynamics(CFD) code and an innovative full-structured dynamic grid method applied to update the three-degree-of-freedom(3-DOF) rigid structure motions. As a time-marching scheme, the trapezoid analogue integral method is used to update the time integration combined with remeshing at each time step.The application of full-structured mesh elements can prevent grids distortion or deformation caused by large-scale movement and improve the stability of calculation. In movable regions, each moving zone is specified with particular motion modes(sway, heave and roll). A series of experimental studies are carried out to validate the performance of the floating body and verify the accuracy of the proposed numerical model. The results are systematically assessed in terms of wave coefficients, mooring line forces, velocity streamlines and the 3-DOF motions of the floating breakwater. When compared with the wave coefficient solutions, excellent agreements are achieved between the computed and experimental data, except in the vicinity of resonant frequency. The velocity streamlines and wave profile movement in the fluid field can also be reproduced using this numerical model.展开更多
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 me...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 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, an...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.展开更多
A rotating stopper-rod technique was proposed to suppress the formation of free-surface vortex in the tundish.The large eddy simulation model coupled with volume of fluid model was developed to study the steel–slag–...A rotating stopper-rod technique was proposed to suppress the formation of free-surface vortex in the tundish.The large eddy simulation model coupled with volume of fluid model was developed to study the steel–slag–gas three-phase flow behavior.The critical slag entrapment height of the free-surface vortex and mass of residual steel were predicted at different rotating speeds(30,60,90 and 120 r/min)of the rotating stopper-rod.The numerical model was verified by water model experiment.The results showed that by rotating the stopper-rod in the opposite direction of the vortex above the submerged entry nozzle,the formation of vortex can be effectively disturbed and the critical height of the free-surface vortex can be reduced.Particularly for the 2nd strand,when the rotating speeds are 30,60,90 and 120 r/min,the critical height of the free-surface vortex above the 2nd strand is 7.3,4.7,6.3 and 7.4 cm,respectively.A reasonable rotating speed should be 60 r/min,which can reduce about 2 tons of residual steel.Other rotating speeds just can reduce about 1.6 tons of residual steel.展开更多
With a target treated as the incompressible Tresca and Mohr-Coulomb material, by assuming that cavity expansion produces plastic-elastic and plastic-cracked-elastic response region, the decay function for the free-sur...With a target treated as the incompressible Tresca and Mohr-Coulomb material, by assuming that cavity expansion produces plastic-elastic and plastic-cracked-elastic response region, the decay function for the free-surface effect is constructed for metal and geological tar- gets, respectively. The forcing function for oblique penetration and perforation is obtained by multiplying the forcing function derived on the basis of infinite target assumption with the de- cay function. Then the projectile is modeled with an explicit transient dynamic finite element code and the target is represented by the forcing function as the pressure boundary condition. This methodology eliminates discretizing the target as well as the need for a complex contact algorithm and is implemented in ABAQUS explicit solver via the user subroutine VDLOAD. It is found that the free-surface effect must be considered in terms of the projectile deformation, residual velocity, projectile trajectory, ricochet limits and critical reverse velocity. The numerical predictions are in good agreement with the available experimental data if the free-surface effect is taken into account.展开更多
A raised panel method based on NURBS (Non-Uniform Rational B-Splines) forfree-surface flows with forward speed is presented. In this generalized panel method, NURBS areemployed to represent the body geometry, disturbe...A raised panel method based on NURBS (Non-Uniform Rational B-Splines) forfree-surface flows with forward speed is presented. In this generalized panel method, NURBS areemployed to represent the body geometry, disturbed free surface, and to express the unknown sourcestrength distribution, on the body surface and above the free surface. Compared with commonhigher-order panel methods, it has no need of adopting local coordinates. NURBS make the geometryrepresentation of the body shape and the wave pattern more precise. Raised panels above the freesurface produce less numerical dispersion error, need less CPU consumption and are helpful andcombined with collocation-point shifting up-stream, can satisfy the radiation condition numerically.By using continuous and discrete Fourier analysis, numerical errors of this method are discussedand a general expression for the errors of numerical damping and dispersion, including the effectsof the vertical distance of singularities to the free surface, the order of singularity distributionrepresented by B-splines in panels, and collocation-point shifting is derived.展开更多
This paper introduces and discusses numerical methods for flee-surface flow simulations and applies a large eddy simula- tion (LES) based flee-surface-resolved CFD method to a couple of flows of hydraulic engineerin...This paper introduces and discusses numerical methods for flee-surface flow simulations and applies a large eddy simula- tion (LES) based flee-surface-resolved CFD method to a couple of flows of hydraulic engineering interest. The advantages, dis- advantages and limitations of the various methods are discussed. The review prioritises interface capturing methods over interface tracking methods, as these have shown themselves to be more generally applicable to viscous flows of practical engineering interest, particularly when complex and rapidly changing surface topologies are encountered. Then, a LES solver that employs the level set method to capture flee-surface deformation in 3-D flows is presented, as are results from two example calculations that concern com- plex low submergence turbulent flows over idealised roughness elements and bluff bodies. The results show that the method is capable of predicting very complex flows that are characterised by strong interactions between the bulk flow and the free-surface, and permits the identification of turbulent events and structures that would be very difficult to measure experimentally.展开更多
The free-surface waves and the flow field due to a body moving on the surfaceof an incompressible viscous fluid of infinite depth were studied analytically. The floating bodywas modeled as a normal point pressure on t...The free-surface waves and the flow field due to a body moving on the surfaceof an incompressible viscous fluid of infinite depth were studied analytically. The floating bodywas modeled as a normal point pressure on the free surface. Based on the Oseen approximation forgoverning equations and the linearity assumption for boundary conditions, the exact solutions inintegral form for the free-surface elevation, the velocities and the pressure were given. Byemploying Lighthill's two-stage scheme, the asymptotic representations in far field for largeReynolds numbers were derived explicitly. The effect of viscosity on the wave profiles was expressedby an exponential decay factor, which removes the singular behavior predicted by the potentialtheory.展开更多
In this paper, the transient fluid resonance phenomenon inside a narrow gap between two adjacent boxes excited by the incident focused waves with various spectral peak periods and focused wave amplitudes is simulated ...In this paper, the transient fluid resonance phenomenon inside a narrow gap between two adjacent boxes excited by the incident focused waves with various spectral peak periods and focused wave amplitudes is simulated by utilizing the open-sourced computational fluid dynamics software, Open FOAM. The weather-side box is allowed to heave freely under the action of waves, and the lee-side box keeps fixed. This paper mainly focuses on how both the spectral peak period and the focused wave amplitude affect the free-surface amplification inside the gap, the motion of the weather-side box, and the wave loads(including the vertical and the horizontal wave forces) acting on both boxes.For comparison, another two-box system with both boxes fixed is also considered as a control group. It is found that the motion of the weather-side box significantly changes the characteristics of the transient gap resonance, and it would cause that the fluid resonant period becomes 1.4-1.6 times that of the two-box system with both boxes fixed.All the concerned physical quantities(i.e., the free-surface amplification in the gap, the motion of the weather-side box, the wave loads) are found to closely depend on both the spectral peak period and the focused wave amplitude.展开更多
Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting...Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.展开更多
To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as i...To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as it can form a matrix with Tocplitz and Hankel properties. In this paper, a method is proposed to decompose the finite depth Green function into two terms, which can form matrices with the Toeplitz and a Hankel properties respectively. Then, a pFFT method for finite depth problems is developed. Based on the pFFT method, a numerical code pFFT-HOBEM is developed with the discretization of high order elements. The model is validated, and examinations on the computing efficiency and memory requirement of the new method have also been carried out. It shows that the new method has the same advantages as that for infinite depth.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51579122,51609109,and 51622902)the Natural Science Found of Jiangsu Province(Grant No.BK20160556)+1 种基金the University Natural Science Research Project of Jiangsu Province(Grant No.16kjb70003)the Key Lab Foundation for Advanced Manufacturing Technology of Jiangsu Province(Grant No.CJ1506)
文摘This paper investigates the hydrodynamic performance of a cylindrical-dual or rectangular-single pontoon floating breakwater using the numerical method and experimental study. The numerical simulation work is based on the multi-physics computational fluid dynamics(CFD) code and an innovative full-structured dynamic grid method applied to update the three-degree-of-freedom(3-DOF) rigid structure motions. As a time-marching scheme, the trapezoid analogue integral method is used to update the time integration combined with remeshing at each time step.The application of full-structured mesh elements can prevent grids distortion or deformation caused by large-scale movement and improve the stability of calculation. In movable regions, each moving zone is specified with particular motion modes(sway, heave and roll). A series of experimental studies are carried out to validate the performance of the floating body and verify the accuracy of the proposed numerical model. The results are systematically assessed in terms of wave coefficients, mooring line forces, velocity streamlines and the 3-DOF motions of the floating breakwater. When compared with the wave coefficient solutions, excellent agreements are achieved between the computed and experimental data, except in the vicinity of resonant frequency. The velocity streamlines and wave profile movement in the fluid field can also be reproduced using this numerical model.
文摘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 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.
基金the National Natural Science Foundation of China(52004191)the China Postdoctoral Science Foundation(2022M711120)+1 种基金the Science and Technology Research Project of Education Department of Hubei Province(B2022020)Besides,the numerical calculation is supported by High-Performance Computing Center of Wuhan University of Science and Technology.
文摘A rotating stopper-rod technique was proposed to suppress the formation of free-surface vortex in the tundish.The large eddy simulation model coupled with volume of fluid model was developed to study the steel–slag–gas three-phase flow behavior.The critical slag entrapment height of the free-surface vortex and mass of residual steel were predicted at different rotating speeds(30,60,90 and 120 r/min)of the rotating stopper-rod.The numerical model was verified by water model experiment.The results showed that by rotating the stopper-rod in the opposite direction of the vortex above the submerged entry nozzle,the formation of vortex can be effectively disturbed and the critical height of the free-surface vortex can be reduced.Particularly for the 2nd strand,when the rotating speeds are 30,60,90 and 120 r/min,the critical height of the free-surface vortex above the 2nd strand is 7.3,4.7,6.3 and 7.4 cm,respectively.A reasonable rotating speed should be 60 r/min,which can reduce about 2 tons of residual steel.Other rotating speeds just can reduce about 1.6 tons of residual steel.
基金Project supported by the National Natural Science Foundations of China(Nos.51321064,51178461 and 51378015)
文摘With a target treated as the incompressible Tresca and Mohr-Coulomb material, by assuming that cavity expansion produces plastic-elastic and plastic-cracked-elastic response region, the decay function for the free-surface effect is constructed for metal and geological tar- gets, respectively. The forcing function for oblique penetration and perforation is obtained by multiplying the forcing function derived on the basis of infinite target assumption with the de- cay function. Then the projectile is modeled with an explicit transient dynamic finite element code and the target is represented by the forcing function as the pressure boundary condition. This methodology eliminates discretizing the target as well as the need for a complex contact algorithm and is implemented in ABAQUS explicit solver via the user subroutine VDLOAD. It is found that the free-surface effect must be considered in terms of the projectile deformation, residual velocity, projectile trajectory, ricochet limits and critical reverse velocity. The numerical predictions are in good agreement with the available experimental data if the free-surface effect is taken into account.
文摘A raised panel method based on NURBS (Non-Uniform Rational B-Splines) forfree-surface flows with forward speed is presented. In this generalized panel method, NURBS areemployed to represent the body geometry, disturbed free surface, and to express the unknown sourcestrength distribution, on the body surface and above the free surface. Compared with commonhigher-order panel methods, it has no need of adopting local coordinates. NURBS make the geometryrepresentation of the body shape and the wave pattern more precise. Raised panels above the freesurface produce less numerical dispersion error, need less CPU consumption and are helpful andcombined with collocation-point shifting up-stream, can satisfy the radiation condition numerically.By using continuous and discrete Fourier analysis, numerical errors of this method are discussedand a general expression for the errors of numerical damping and dispersion, including the effectsof the vertical distance of singularities to the free surface, the order of singularity distributionrepresented by B-splines in panels, and collocation-point shifting is derived.
基金supported by the UK Engineering and Physical Sciences Research Council(EPSRC)
文摘This paper introduces and discusses numerical methods for flee-surface flow simulations and applies a large eddy simula- tion (LES) based flee-surface-resolved CFD method to a couple of flows of hydraulic engineering interest. The advantages, dis- advantages and limitations of the various methods are discussed. The review prioritises interface capturing methods over interface tracking methods, as these have shown themselves to be more generally applicable to viscous flows of practical engineering interest, particularly when complex and rapidly changing surface topologies are encountered. Then, a LES solver that employs the level set method to capture flee-surface deformation in 3-D flows is presented, as are results from two example calculations that concern com- plex low submergence turbulent flows over idealised roughness elements and bluff bodies. The results show that the method is capable of predicting very complex flows that are characterised by strong interactions between the bulk flow and the free-surface, and permits the identification of turbulent events and structures that would be very difficult to measure experimentally.
文摘The free-surface waves and the flow field due to a body moving on the surfaceof an incompressible viscous fluid of infinite depth were studied analytically. The floating bodywas modeled as a normal point pressure on the free surface. Based on the Oseen approximation forgoverning equations and the linearity assumption for boundary conditions, the exact solutions inintegral form for the free-surface elevation, the velocities and the pressure were given. Byemploying Lighthill's two-stage scheme, the asymptotic representations in far field for largeReynolds numbers were derived explicitly. The effect of viscosity on the wave profiles was expressedby an exponential decay factor, which removes the singular behavior predicted by the potentialtheory.
基金supported by the National Natural Science Foundation of China (Grant No. 51911530205)the Natural Science Foundation of Jiangsu Province (Grant Nos. BK20201455 and BK20210885)+5 种基金the Key Laboratory of PortWaterway and Sedimentation Engineering of MOT (Grant No. YK222001-2)the Shandong Provincial Key Laboratory of Ocean Engineering (Grant No. kloe202010)the Key R&D Projects in Guangdong Province (Grant No. 2020B1111500001)the Qing Lan Project of Jiangsu Universitiesthe Royal Society (Grant No. IECNSFC181321) for providing partial support for this work
文摘In this paper, the transient fluid resonance phenomenon inside a narrow gap between two adjacent boxes excited by the incident focused waves with various spectral peak periods and focused wave amplitudes is simulated by utilizing the open-sourced computational fluid dynamics software, Open FOAM. The weather-side box is allowed to heave freely under the action of waves, and the lee-side box keeps fixed. This paper mainly focuses on how both the spectral peak period and the focused wave amplitude affect the free-surface amplification inside the gap, the motion of the weather-side box, and the wave loads(including the vertical and the horizontal wave forces) acting on both boxes.For comparison, another two-box system with both boxes fixed is also considered as a control group. It is found that the motion of the weather-side box significantly changes the characteristics of the transient gap resonance, and it would cause that the fluid resonant period becomes 1.4-1.6 times that of the two-box system with both boxes fixed.All the concerned physical quantities(i.e., the free-surface amplification in the gap, the motion of the weather-side box, the wave loads) are found to closely depend on both the spectral peak period and the focused wave amplitude.
基金supported by the National Natural Science Foundation of China(51809169,51879159)Chang Jiang Scholars Program(T2014099)+2 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)National Key Research and Development Program of China(2019YFB1704203,2019YFC0312400).
文摘Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.
基金supported by the National Natural Science Foundation of China(Grant Nos.51490672 and 51379032)
文摘To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as it can form a matrix with Tocplitz and Hankel properties. In this paper, a method is proposed to decompose the finite depth Green function into two terms, which can form matrices with the Toeplitz and a Hankel properties respectively. Then, a pFFT method for finite depth problems is developed. Based on the pFFT method, a numerical code pFFT-HOBEM is developed with the discretization of high order elements. The model is validated, and examinations on the computing efficiency and memory requirement of the new method have also been carried out. It shows that the new method has the same advantages as that for infinite depth.