Large eddy simulation cooperated with the second order full extension ETG(Euler-Taylor-Galerkin) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ra...Large eddy simulation cooperated with the second order full extension ETG(Euler-Taylor-Galerkin) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ratio of (1.5.) The second order full extension ETG finite element method was developed by Wang and He. By means of Taylor expansion of terms containing time derivative, time derivative is replaced by space derivative. The function of it is equal to introducing an artificial viscosity term. The streamlines of the flow at different moments were obtained. The time history of drag coefficient, lift coefficient and the streamwise velocity on the symmetrical points were presented. Furthermore, the symmetrical problem of the frequency spectrum of flow around two square cylinders arranged side by side were studied by using the spectral analysis technology. The data obtained at the initial stage are excluded in order to avoid the influence of initial condition on the results. The power spectrums of drag coefficient, lift coefficient, the streamwise velocity on the symmetrical points were analyzed respectively. The results show that although the time domain process of dynamic parameters is non-symmetrical, the frequency domain process of them is symmetrical under the symmetrical boundary conditions.展开更多
The numerical method is used to calculate the flow around two square cylinders arranged side-by-side and the mean and fluctuating aerodynamic forces, and Strouhal numbers and power spectrum of lift force and drag forc...The numerical method is used to calculate the flow around two square cylinders arranged side-by-side and the mean and fluctuating aerodynamic forces, and Strouhal numbers and power spectrum of lift force and drag force are obtained. An improved MAC method proposed by Chen Suqin et al.,which uses three order upwind scheme to discretize the convection term and uses multigrid method to solve the Poisson equation for pressure is applied to simulate the flow around two square cylinders arranged side-by-side. Results show that the interference characteristic of two square cylinders arranged side-by-side is completely different with the different spacing ratio. When the spacing ratio is smaller than a certain critical value, the gap flow between two cylinders is biased to one side in a stable or unstable manner.展开更多
Helical strake is a widely-used device for passive flow-induced vibration(FIV)control of cylindrical structures.It is omnidirectional and can effectively reduce FIV response amplitude.Studies on the passive FIV contro...Helical strake is a widely-used device for passive flow-induced vibration(FIV)control of cylindrical structures.It is omnidirectional and can effectively reduce FIV response amplitude.Studies on the passive FIV control for cylindrical structures are mainly concerned with a single isolated cylinder,while the influence of wake interference between multiple cylinders on FIV suppression devices is less considered up to now.In engineering applications,multiple flexible cylinders with large aspect ratios can be subjected to complex flow forces,and the effects of wake interference are obvious.The FIV suppression effect of helical strake of a common configuration(17.5D pitch and 0.25D height,where D is the cylinder diameter)in two staggered cylinders system is still unknown.This paper systematically studied the FIV response of multiple cylinders system fitted with the helical strakes by model tests.The relative spatial position of the two cylinders is fixed at S=3.0D and T=8.0D,which ensures the cylindrical structures in the flow interference region.The experimental results show that the helical strakes effectively reduce the FIV response on staggered upstream cylinder,and the suppression efficiency is barely affected by the smooth or straked downstream cylinder.The corresponding FIV suppression efficiency on the downstream cylinder is remarkably reduced by the influence of the upstream wake flow.The wake-induced vibration(WIV)phenomenon is not observed on the staggered downstream cylinder,which normally occurs on the downstream straked cylinder in a tandem arrangement.展开更多
-Based on the extended Morison Equation and model tests, the in-line forces on small square cylinders caused by waves (regular and irregular) and currents are analyzed in detail in this paper. The hydrodynamic coeffic...-Based on the extended Morison Equation and model tests, the in-line forces on small square cylinders caused by waves (regular and irregular) and currents are analyzed in detail in this paper. The hydrodynamic coefficient CD and Cu related to KC number and the effect of direction of wave incidence are also given, which can be used in engineering practice.展开更多
The spatial evolution of vortices and transition to three-dimensionality in the wake of two circular cylinders in tandem arrangement have been numerically studied. An improved virtual body method developed from the vi...The spatial evolution of vortices and transition to three-dimensionality in the wake of two circular cylinders in tandem arrangement have been numerically studied. An improved virtual body method developed from the virtual boundary method is used here. A Reynolds number range between 220 and 270 has been considered, and the spacing between two cylinders is selected as L/D=3 and L/D=3.5. When L/D=3, the secondary vortices of Mode-A are seen to appear at Re=240 and persist over the range of the Reynolds number of 240~270. When L/D=3.5, the similar critical Reynolds number has been found at Re=250. No obvious discontinuity has been found in the Strouhal-Reynolds number relationship, and this is different from three-dimensional flow around a single cylinder at the critical Reynolds number. The spanwise wavelength is about four times the diameter of the cylinder, and it is the characteristic wavelength for Mode-A instability. This paper can give some foremost insight into the three-dimensional instability of flow by complicated geometrical configuration.展开更多
Few studies jointly investigate thermal and turbulent effects. In general, these subjects are treated separately. The purpose of this paper is to use the Immersed Boundary Method (IBM) coupled with the Virtual Physica...Few studies jointly investigate thermal and turbulent effects. In general, these subjects are treated separately. The purpose of this paper is to use the Immersed Boundary Method (IBM) coupled with the Virtual Physical Model (VPM) to investigate incompressible two-dimensional Newtonian flow around a heated square cylinder at constant temperature on its surface with forced convection and turbulence. The VPM model dynamically evaluates the force that the fluid exerts on the immersed surface and the thermal exchange between both in the Reynolds numbers (Re) window 40 ≤ Re ≤ 5×103 . For simulations of turbulence the Smagorinsky and Spalart-Allmaras models are used. The first model uses the Large Eddy Simulation (LES) methodology and is based on the local equilibrium hypothesis for small scales associated with the Boussinesq hypothesis, such that the energy injected into the spectrum of the turbulence balances the energy dissipated by convective effects. The second model uses the concept Unsteady Reynolds Averaged Navier-Stokes Equations (URANS), with only one transport equation for turbulent viscosity, being calibrated in pressure gradient layers. The goal of this work is to analyse the combination of the heat-transfer phenomena with the turbulence for the thermo-fluid-structure interaction in a square cylinder. For this, it was developed a C/C++ code that requires low computational costs in regards to memory and computer facilities. It is observed that, with the increase of the Reynolds number, an increase of the drag coefficient occurs, as well as reinforces the influence of the pressure distribution downstream of the cylinder, which is strongly influenced by the formation and detachment of vortices on the upper and lower sides of the square cylinder.展开更多
Measurements were performed using Particle Image Velocimetry (PIV) to analyze the modification of flow by the combined effects of the rotation and the Reynolds number on the flow past two rotating circular cylinders...Measurements were performed using Particle Image Velocimetry (PIV) to analyze the modification of flow by the combined effects of the rotation and the Reynolds number on the flow past two rotating circular cylinders in a side-by-side-arrangement at a range of 425 〈 Re ≤ 1130,0 ≤α ≤4 ( α is the rotational speed) at one gap spacing of T / d = 1.11 (T and d are the distance between the centers of two cylinders and the cylinder diameter, respectively). A new Immersed-Lattice Boltzmann Method (ILBM) scheme was used to study the effect of the gap spacing on the flow. The results show that the vortex shedding is suppressed as rotational speed increases. The flow reaches a steady state when the vortex shedding for both cylinders is completely suppressed at critical rotational speed. As the rotational speed further increases, the separation phenomenon in the boundary layers disappears at the attachment rotational speed. The critical rotational speed and attachment rotational speed become small as Reynolds number increases. The absolute rotational speed of cylinders should be large at same critical rotational speed and attachment rotational speed in the case of large Reynolds number. The gap spacing has an important role in changing the pattern of vortex shedding. It is very different in the mechanism of vortex shedding suppression for the flows around two rotating cylinders and single rotating cylinder.展开更多
The wake induced vibration(WIV)of a one-and two-degree-of-freedom(1DOF,2DOF)downstream wave-cone cylinder(WCC)behind a stationary equal-size upstream wave-cone cylinder in the staggered arrangement is numerically inve...The wake induced vibration(WIV)of a one-and two-degree-of-freedom(1DOF,2DOF)downstream wave-cone cylinder(WCC)behind a stationary equal-size upstream wave-cone cylinder in the staggered arrangement is numerically investigated at subcritical Reynolds number of 3900 by using shear stressed transfer(SST)k-ωturbulence model.The streamwise pitch ratios(P/Dm)vary from 4 to 6 with a fixed incident angleα=8°.Experimental measurements were also performed for the validation of the present numerical models.It is found that the largest vibration amplitude in crossflow direction occurred at P/Dm=4,Ur=8 with small difference of streamwise vibration at P/Dm=4,6.Different from single wavy-cone cylinder(SWCC),the downstream flexible one of a pair staggered WCCs got larger vibration amplitude during phase switching stage instead of in-phase stage.The upstream wake will suppress the triple frequency of main frequency in the power spectra density(PSD)functions of Cl but stimulate the double one of that.An intriguing vibration mechanism happened in all 2DOF cases where the trajectory of downstream WCC is a significant ellipse rather than a figure of 8.The transformation of phase switching and the variation of the main frequency of drag coefficient(Cd)can be explained by the vortex-shedding modes of downstream WCC.展开更多
文摘Large eddy simulation cooperated with the second order full extension ETG(Euler-Taylor-Galerkin) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ratio of (1.5.) The second order full extension ETG finite element method was developed by Wang and He. By means of Taylor expansion of terms containing time derivative, time derivative is replaced by space derivative. The function of it is equal to introducing an artificial viscosity term. The streamlines of the flow at different moments were obtained. The time history of drag coefficient, lift coefficient and the streamwise velocity on the symmetrical points were presented. Furthermore, the symmetrical problem of the frequency spectrum of flow around two square cylinders arranged side by side were studied by using the spectral analysis technology. The data obtained at the initial stage are excluded in order to avoid the influence of initial condition on the results. The power spectrums of drag coefficient, lift coefficient, the streamwise velocity on the symmetrical points were analyzed respectively. The results show that although the time domain process of dynamic parameters is non-symmetrical, the frequency domain process of them is symmetrical under the symmetrical boundary conditions.
文摘The numerical method is used to calculate the flow around two square cylinders arranged side-by-side and the mean and fluctuating aerodynamic forces, and Strouhal numbers and power spectrum of lift force and drag force are obtained. An improved MAC method proposed by Chen Suqin et al.,which uses three order upwind scheme to discretize the convection term and uses multigrid method to solve the Poisson equation for pressure is applied to simulate the flow around two square cylinders arranged side-by-side. Results show that the interference characteristic of two square cylinders arranged side-by-side is completely different with the different spacing ratio. When the spacing ratio is smaller than a certain critical value, the gap flow between two cylinders is biased to one side in a stable or unstable manner.
基金the National Natural Science Foundation of China(Grant No.51979193)the Natural Science Foundation of Tianjin(Grant No.20JCYBJC00890).
文摘Helical strake is a widely-used device for passive flow-induced vibration(FIV)control of cylindrical structures.It is omnidirectional and can effectively reduce FIV response amplitude.Studies on the passive FIV control for cylindrical structures are mainly concerned with a single isolated cylinder,while the influence of wake interference between multiple cylinders on FIV suppression devices is less considered up to now.In engineering applications,multiple flexible cylinders with large aspect ratios can be subjected to complex flow forces,and the effects of wake interference are obvious.The FIV suppression effect of helical strake of a common configuration(17.5D pitch and 0.25D height,where D is the cylinder diameter)in two staggered cylinders system is still unknown.This paper systematically studied the FIV response of multiple cylinders system fitted with the helical strakes by model tests.The relative spatial position of the two cylinders is fixed at S=3.0D and T=8.0D,which ensures the cylindrical structures in the flow interference region.The experimental results show that the helical strakes effectively reduce the FIV response on staggered upstream cylinder,and the suppression efficiency is barely affected by the smooth or straked downstream cylinder.The corresponding FIV suppression efficiency on the downstream cylinder is remarkably reduced by the influence of the upstream wake flow.The wake-induced vibration(WIV)phenomenon is not observed on the staggered downstream cylinder,which normally occurs on the downstream straked cylinder in a tandem arrangement.
文摘-Based on the extended Morison Equation and model tests, the in-line forces on small square cylinders caused by waves (regular and irregular) and currents are analyzed in detail in this paper. The hydrodynamic coefficient CD and Cu related to KC number and the effect of direction of wave incidence are also given, which can be used in engineering practice.
基金This work was financially supported by the National Natural Science Foundation of China (Grant No. 10272094)
文摘The spatial evolution of vortices and transition to three-dimensionality in the wake of two circular cylinders in tandem arrangement have been numerically studied. An improved virtual body method developed from the virtual boundary method is used here. A Reynolds number range between 220 and 270 has been considered, and the spacing between two cylinders is selected as L/D=3 and L/D=3.5. When L/D=3, the secondary vortices of Mode-A are seen to appear at Re=240 and persist over the range of the Reynolds number of 240~270. When L/D=3.5, the similar critical Reynolds number has been found at Re=250. No obvious discontinuity has been found in the Strouhal-Reynolds number relationship, and this is different from three-dimensional flow around a single cylinder at the critical Reynolds number. The spanwise wavelength is about four times the diameter of the cylinder, and it is the characteristic wavelength for Mode-A instability. This paper can give some foremost insight into the three-dimensional instability of flow by complicated geometrical configuration.
基金the partial support by CMUP(UID/MAT/00144/2013),which is funded by FCT(Portugal)with national(MCTES)and European structural funds(FEDER),under the partnership agreement PT2020-ext.to 2018the financial support by CAPES(Brazil)SG acknowledges the Project STRIDE-NORTE-01-0145-FEDER-000033,funded by ERDF NORTE 2020.
文摘Few studies jointly investigate thermal and turbulent effects. In general, these subjects are treated separately. The purpose of this paper is to use the Immersed Boundary Method (IBM) coupled with the Virtual Physical Model (VPM) to investigate incompressible two-dimensional Newtonian flow around a heated square cylinder at constant temperature on its surface with forced convection and turbulence. The VPM model dynamically evaluates the force that the fluid exerts on the immersed surface and the thermal exchange between both in the Reynolds numbers (Re) window 40 ≤ Re ≤ 5×103 . For simulations of turbulence the Smagorinsky and Spalart-Allmaras models are used. The first model uses the Large Eddy Simulation (LES) methodology and is based on the local equilibrium hypothesis for small scales associated with the Boussinesq hypothesis, such that the energy injected into the spectrum of the turbulence balances the energy dissipated by convective effects. The second model uses the concept Unsteady Reynolds Averaged Navier-Stokes Equations (URANS), with only one transport equation for turbulent viscosity, being calibrated in pressure gradient layers. The goal of this work is to analyse the combination of the heat-transfer phenomena with the turbulence for the thermo-fluid-structure interaction in a square cylinder. For this, it was developed a C/C++ code that requires low computational costs in regards to memory and computer facilities. It is observed that, with the increase of the Reynolds number, an increase of the drag coefficient occurs, as well as reinforces the influence of the pressure distribution downstream of the cylinder, which is strongly influenced by the formation and detachment of vortices on the upper and lower sides of the square cylinder.
基金Project supported by the Major Program of the National Natural Science Foundation of China (Grant No.10632070)
文摘Measurements were performed using Particle Image Velocimetry (PIV) to analyze the modification of flow by the combined effects of the rotation and the Reynolds number on the flow past two rotating circular cylinders in a side-by-side-arrangement at a range of 425 〈 Re ≤ 1130,0 ≤α ≤4 ( α is the rotational speed) at one gap spacing of T / d = 1.11 (T and d are the distance between the centers of two cylinders and the cylinder diameter, respectively). A new Immersed-Lattice Boltzmann Method (ILBM) scheme was used to study the effect of the gap spacing on the flow. The results show that the vortex shedding is suppressed as rotational speed increases. The flow reaches a steady state when the vortex shedding for both cylinders is completely suppressed at critical rotational speed. As the rotational speed further increases, the separation phenomenon in the boundary layers disappears at the attachment rotational speed. The critical rotational speed and attachment rotational speed become small as Reynolds number increases. The absolute rotational speed of cylinders should be large at same critical rotational speed and attachment rotational speed in the case of large Reynolds number. The gap spacing has an important role in changing the pattern of vortex shedding. It is very different in the mechanism of vortex shedding suppression for the flows around two rotating cylinders and single rotating cylinder.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.11972268,12372232).
文摘The wake induced vibration(WIV)of a one-and two-degree-of-freedom(1DOF,2DOF)downstream wave-cone cylinder(WCC)behind a stationary equal-size upstream wave-cone cylinder in the staggered arrangement is numerically investigated at subcritical Reynolds number of 3900 by using shear stressed transfer(SST)k-ωturbulence model.The streamwise pitch ratios(P/Dm)vary from 4 to 6 with a fixed incident angleα=8°.Experimental measurements were also performed for the validation of the present numerical models.It is found that the largest vibration amplitude in crossflow direction occurred at P/Dm=4,Ur=8 with small difference of streamwise vibration at P/Dm=4,6.Different from single wavy-cone cylinder(SWCC),the downstream flexible one of a pair staggered WCCs got larger vibration amplitude during phase switching stage instead of in-phase stage.The upstream wake will suppress the triple frequency of main frequency in the power spectra density(PSD)functions of Cl but stimulate the double one of that.An intriguing vibration mechanism happened in all 2DOF cases where the trajectory of downstream WCC is a significant ellipse rather than a figure of 8.The transformation of phase switching and the variation of the main frequency of drag coefficient(Cd)can be explained by the vortex-shedding modes of downstream WCC.