The free-surface vortex is a rotational flow phenomenon characterized by two-phase coupling,formed by the rupture of surface fluid in the final stage of discharge.It is a significant concept with broad applications in...The free-surface vortex is a rotational flow phenomenon characterized by two-phase coupling,formed by the rupture of surface fluid in the final stage of discharge.It is a significant concept with broad applications in engineering fields like metallurgy and hydraulics.The basic concepts and characteristics of free-surface vortices were introduced,and their hazards in various fields were discussed.The development of theoretical and numerical models over recent decades was reviewed,and the factors affecting vortex formation and existing suppression methods were outlined.Finally,the key challenges and focus areas on the study of free-surface vortex were summarized.With the ongoing advancements in computational fluid dynamics and experimental technology,research on free-surface vortices will become more in depth and precise.Additionally,interdisciplinary cooperation and technological innovation are expected to achieve precise control and optimal design of free-surface vortices,offering more efficient and sustainable solutions for metallurgy and related engineering fields.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Floating zone technique is a crucible-free process for growth of high quality single crystals. Unstable thermocapillary convection is a typical phenomenon during the process under microgravity. Therefore, it is very i...Floating zone technique is a crucible-free process for growth of high quality single crystals. Unstable thermocapillary convection is a typical phenomenon during the process under microgravity. Therefore, it is very important to investigate the instability of thermocapillary convection in liquid bridges with deformable free-surface under microgravity. In this works, the Volume of Fluid(VOF) method is employed to track the free-surface movement. The results are presented as the behavior of flow structure and temperature distribution of the molten zone. The impact of Marangoni number(Ma) is also investigated on free-surface deformation as well as the instability of thermocapillary convection. The free-surface exhibits a noticeable axisymmetric(but it is non-centrosymmetric) and elliptical shape along the circumferential direction. This specific surface shape presents a typical narrow ‘neck-shaped' structure with convex at two ends of the zone and concave at the mid-plane along the axial direction. At both θ = 0° and θ = 90°, the deformation ratio ξ increases rapidly with Ma at first, and then increases slowly. Moreover, the hydrothermal wave number m and the instability of thermocapillary convection increase with Ma.展开更多
The interaction of laminar flows with free surface waves generated by submerged bodies in an incompressible viscous fluid of infinite depth is investigated analytically. The analysis is based on the linearized Navier-...The interaction of laminar flows with free surface waves generated by submerged bodies in an incompressible viscous fluid of infinite depth is investigated analytically. The analysis is based on the linearized Navier-Stokes equations for disturbed flows. The kinematic and dynamic boundary conditions are linearized for the small-amplitude free-surface waves, and the initial values of the flow are taken to be those of the steady state cases. The submerged bodies are mathematically represented by fundamental singularities of viscous flows. The asymptotic representations for unsteady free-surface waves produced by the Stokeslets and Oseenlets are derived analytically. It is found that the unsteady waves generated by a body consist of steady-state and transient responses. As time tends to infinity, the transient waves vanish due to the presence of a viscous decay factor. Thus, an ultimate steady state can be attained.展开更多
Making full use of functional variations with variable domain for handling both the free surface and unknown part of the dam wall,a family of variational principles(VPs)for the hybrid(direct-inverse)problem of free su...Making full use of functional variations with variable domain for handling both the free surface and unknown part of the dam wall,a family of variational principles(VPs)for the hybrid(direct-inverse)problem of free surface spillway flow under gravity is formulated in terms of the potential function Φ.This theory offers new rational and versatile ways for dam de- sign and a sound theoretical basis for numerical solution,e.g.by the finite element method (FEM).展开更多
We develop a lattice Boltzmann method for modeling free-surface temperature dispersion in the shallow water flows.The governing equations are derived from the incompressible Navier-Stokes equations with assumptions of...We develop a lattice Boltzmann method for modeling free-surface temperature dispersion in the shallow water flows.The governing equations are derived from the incompressible Navier-Stokes equations with assumptions of shallow water flows including bed frictions,eddy viscosity,wind shear stresses and Coriolis forces.The thermal effects are incorporated in the momentum equation by using a Boussinesq approximation.The dispersion of free-surface temperature is modelled by an advection-diffusion equation.Two distribution functions are used in the lattice Boltzmann method to recover the flow and temperature variables using the same lattice structure.Neither upwind discretization procedures nor Riemann problem solvers are needed in discretizing the shallow water equations.In addition,the source terms are straightforwardly included in the model without relying on well-balanced techniques to treat flux gradients and source terms.We validate the model for a class of problems with known analytical solutions and we also present numerical results for sea-surface temperature distribution in the Strait of Gibraltar.展开更多
In this work the feasibility of a numerical wave tank using a hybrid particle-mesh method is investigated.Based on the fluid implicit particle method（FLIP）a formulation for the hybrid method is presented for incompr...In this work the feasibility of a numerical wave tank using a hybrid particle-mesh method is investigated.Based on the fluid implicit particle method（FLIP）a formulation for the hybrid method is presented for incompressible multiphase flows involving large density jumps and wave generating boundaries.The performance of the method is assessed for a standing wave and for the generation and propagation of a solitary wave over a flat and a sloping bed.A comparison is made with results obtained with a well-established SPH package.The tests demonstrate that the method is a promising and attractive tool for simulating the nearshore propagation of waves.展开更多
A length scale refinement study is a standard practice to ensure the independence of a numerical model on spatial approximations.For smoothed particle hydrodynamics(SPH),the process of length scale refinement study te...A length scale refinement study is a standard practice to ensure the independence of a numerical model on spatial approximations.For smoothed particle hydrodynamics(SPH),the process of length scale refinement study tends to be conducted based on experience.A challenge of defining a universal length scale refinement strategy is the existence of two length scales–particle spacing and smoothing length.Despite the challenge,further investigations of the impact of different refinement strategies should be continually conducted to improve the reliability of practical SPH applications on 3-D free-surface flows.In this study,a conventional strategy and a novel coupled refinement strategy are used to investigate the convergence of SPH simulations for free-surface flows using a standard SPH scheme available in an open-source framework.The two case studies are a dam break flow and a lesser-known stable regime water flow inside a rotating drum with lifters.Validations are conducted using existing data from literature for the dam break flow and laser Doppler velocimetry(LDV)measurements for the rotating drum flow.The investigation shows that the proposed coupled length scale refinement strategy does not offer a significant improvement for the SPH model of the dam break flow comparing with the conventional strategy.On the other hand,the stable regime rotating drum fluid flow shows that both refinement strategies are not sufficient to tackle SPH’s on-going fundamental challenge of accurately predicting the flow field of complex 3-D turbulent flows with free surfaces.展开更多
This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunzia...This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.展开更多
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.展开更多
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.展开更多
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.展开更多
The numerical simulation of wake and flee-surface flow around ships is a complex topic that involves multiple tasks： the generation of an optimal computational grid and the development of numerical algorithms capable...The numerical simulation of wake and flee-surface flow around ships is a complex topic that involves multiple tasks： the generation of an optimal computational grid and the development of numerical algorithms capable to predict the flow field around a hull. In this paper, a numerical framework is developed aimed at high-resolution CFD simulations of turbulent, free-surface flows around ship hulls. The framework consists in the concatenation of ＂tools＂, partly available in the open-source finite volume library OpenFOAM. A novel, flexible mesh-generation algorithm is presented, capable of producing high-quality computational grids for free-surface ship hydrodynamics. The numerical frame work is used to solve some benchmark problems, providing results that are in excellent agreement with the experimental measures.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.52474339)Anhui Province Outstanding Research and Innovation Team in Higher Education Institutions(Grant No.2022AH010024).
文摘The free-surface vortex is a rotational flow phenomenon characterized by two-phase coupling,formed by the rupture of surface fluid in the final stage of discharge.It is a significant concept with broad applications in engineering fields like metallurgy and hydraulics.The basic concepts and characteristics of free-surface vortices were introduced,and their hazards in various fields were discussed.The development of theoretical and numerical models over recent decades was reviewed,and the factors affecting vortex formation and existing suppression methods were outlined.Finally,the key challenges and focus areas on the study of free-surface vortex were summarized.With the ongoing advancements in computational fluid dynamics and experimental technology,research on free-surface vortices will become more in depth and precise.Additionally,interdisciplinary cooperation and technological innovation are expected to achieve precise control and optimal design of free-surface vortices,offering more efficient and sustainable solutions for metallurgy and related engineering fields.
基金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.
基金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.
基金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.
文摘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.
基金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.
文摘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.
文摘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 National Natural Science Foundation of China(Grant Number 51276089)
文摘Floating zone technique is a crucible-free process for growth of high quality single crystals. Unstable thermocapillary convection is a typical phenomenon during the process under microgravity. Therefore, it is very important to investigate the instability of thermocapillary convection in liquid bridges with deformable free-surface under microgravity. In this works, the Volume of Fluid(VOF) method is employed to track the free-surface movement. The results are presented as the behavior of flow structure and temperature distribution of the molten zone. The impact of Marangoni number(Ma) is also investigated on free-surface deformation as well as the instability of thermocapillary convection. The free-surface exhibits a noticeable axisymmetric(but it is non-centrosymmetric) and elliptical shape along the circumferential direction. This specific surface shape presents a typical narrow ‘neck-shaped' structure with convex at two ends of the zone and concave at the mid-plane along the axial direction. At both θ = 0° and θ = 90°, the deformation ratio ξ increases rapidly with Ma at first, and then increases slowly. Moreover, the hydrothermal wave number m and the instability of thermocapillary convection increase with Ma.
文摘The interaction of laminar flows with free surface waves generated by submerged bodies in an incompressible viscous fluid of infinite depth is investigated analytically. The analysis is based on the linearized Navier-Stokes equations for disturbed flows. The kinematic and dynamic boundary conditions are linearized for the small-amplitude free-surface waves, and the initial values of the flow are taken to be those of the steady state cases. The submerged bodies are mathematically represented by fundamental singularities of viscous flows. The asymptotic representations for unsteady free-surface waves produced by the Stokeslets and Oseenlets are derived analytically. It is found that the unsteady waves generated by a body consist of steady-state and transient responses. As time tends to infinity, the transient waves vanish due to the presence of a viscous decay factor. Thus, an ultimate steady state can be attained.
文摘Making full use of functional variations with variable domain for handling both the free surface and unknown part of the dam wall,a family of variational principles(VPs)for the hybrid(direct-inverse)problem of free surface spillway flow under gravity is formulated in terms of the potential function Φ.This theory offers new rational and versatile ways for dam de- sign and a sound theoretical basis for numerical solution,e.g.by the finite element method (FEM).
基金support by the DeutscheForschungsGemeinschaft(DFG)under grant No.KL 1105/9.
文摘We develop a lattice Boltzmann method for modeling free-surface temperature dispersion in the shallow water flows.The governing equations are derived from the incompressible Navier-Stokes equations with assumptions of shallow water flows including bed frictions,eddy viscosity,wind shear stresses and Coriolis forces.The thermal effects are incorporated in the momentum equation by using a Boussinesq approximation.The dispersion of free-surface temperature is modelled by an advection-diffusion equation.Two distribution functions are used in the lattice Boltzmann method to recover the flow and temperature variables using the same lattice structure.Neither upwind discretization procedures nor Riemann problem solvers are needed in discretizing the shallow water equations.In addition,the source terms are straightforwardly included in the model without relying on well-balanced techniques to treat flux gradients and source terms.We validate the model for a class of problems with known analytical solutions and we also present numerical results for sea-surface temperature distribution in the Strait of Gibraltar.
基金The NWO(Netherlands Organisation for Scientific Research)is acknowledged for their support through the JMBC-EM Graduate Programme research grant
文摘In this work the feasibility of a numerical wave tank using a hybrid particle-mesh method is investigated.Based on the fluid implicit particle method（FLIP）a formulation for the hybrid method is presented for incompressible multiphase flows involving large density jumps and wave generating boundaries.The performance of the method is assessed for a standing wave and for the generation and propagation of a solitary wave over a flat and a sloping bed.A comparison is made with results obtained with a well-established SPH package.The tests demonstrate that the method is a promising and attractive tool for simulating the nearshore propagation of waves.
文摘A length scale refinement study is a standard practice to ensure the independence of a numerical model on spatial approximations.For smoothed particle hydrodynamics(SPH),the process of length scale refinement study tends to be conducted based on experience.A challenge of defining a universal length scale refinement strategy is the existence of two length scales–particle spacing and smoothing length.Despite the challenge,further investigations of the impact of different refinement strategies should be continually conducted to improve the reliability of practical SPH applications on 3-D free-surface flows.In this study,a conventional strategy and a novel coupled refinement strategy are used to investigate the convergence of SPH simulations for free-surface flows using a standard SPH scheme available in an open-source framework.The two case studies are a dam break flow and a lesser-known stable regime water flow inside a rotating drum with lifters.Validations are conducted using existing data from literature for the dam break flow and laser Doppler velocimetry(LDV)measurements for the rotating drum flow.The investigation shows that the proposed coupled length scale refinement strategy does not offer a significant improvement for the SPH model of the dam break flow comparing with the conventional strategy.On the other hand,the stable regime rotating drum fluid flow shows that both refinement strategies are not sufficient to tackle SPH’s on-going fundamental challenge of accurately predicting the flow field of complex 3-D turbulent flows with free surfaces.
基金funded by the Italian Ministry of Education,University and Research(MIUR)in the frame of the Departments of Excellence Initiative 2018-2027 attributed to DICAM of the University of Trento(grant L.232/2016)in the frame of the PRIN 2017 project Innovative numerical methods for evolutionary partial differential equations and applications,the PRIN 2022 project High order structure-preserving semi-implicit schemes for hyperbolic equations.D.is member of INdAM GNCS and was also co-funded by the European Union NextGenerationEU(PNRR,Spoke 7 CN HPC).Views and opinions expressed are however those of the author(s)only and do not necessarily reflect those of the European Union or the European Research Council.Neither the European Union nor the granting authority can be held responsible for them.
文摘This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.
基金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 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(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.
文摘The numerical simulation of wake and flee-surface flow around ships is a complex topic that involves multiple tasks： the generation of an optimal computational grid and the development of numerical algorithms capable to predict the flow field around a hull. In this paper, a numerical framework is developed aimed at high-resolution CFD simulations of turbulent, free-surface flows around ship hulls. The framework consists in the concatenation of ＂tools＂, partly available in the open-source finite volume library OpenFOAM. A novel, flexible mesh-generation algorithm is presented, capable of producing high-quality computational grids for free-surface ship hydrodynamics. The numerical frame work is used to solve some benchmark problems, providing results that are in excellent agreement with the experimental measures.
