Verification and validation(V&V)is a helpful tool for evaluating simulation errors,but its application in unsteady cavitating flow remains a challenging issue due to the difficulty in meeting the requirement of an...Verification and validation(V&V)is a helpful tool for evaluating simulation errors,but its application in unsteady cavitating flow remains a challenging issue due to the difficulty in meeting the requirement of an asymptotic range.Hence,a new V&V approach for large eddy simulation(LES)is proposed.This approach offers a viable solution for the error estimation of simulation data that are unable to satisfy the asymptotic range.The simulation errors of cavitating flow around a projectile near the free surface are assessed using the new V&V method.The evident error values are primarily dispersed around the cavity region and free surface.The increasingly intense cavitating flow increases the error magnitudes.In addition,the modeling error magnitudes of the Dynamic Smagorinsky-Lilly model are substantially smaller than that of the Smagorinsky-Lilly model.The present V&V method can capture the decrease in the modeling errors due to model enhancements,further exhibiting its applicability in cavitating flow simulations.Moreover,the monitoring points where the simulation data are beyond the asymptotic range are primarily dispersed near the cavity region,and the number of such points grows as the cavitating flow intensifies.The simulation outcomes also suggest that the re-entrant jet and shedding cavity collapse are the chief sources of vorticity motions,which remarkably affect the simulation accuracy.The results of this study provide a valuable reference for V&V research.展开更多
The electric-static body force model is obtained by solving Maxwell＇s electromagnetic equations. Based on the electro-static model, numerical modeling of flow around a cylinder with a dielectric barrier discharge （...The electric-static body force model is obtained by solving Maxwell＇s electromagnetic equations. Based on the electro-static model, numerical modeling of flow around a cylinder with a dielectric barrier discharge （DBD） plasma effect is also presented. The flow streamlines between the numerical simulation and the particle image velocimetry （PIV） experiment are consistent. According to the numerical simulation, DBD plasma can reduce the drag coefficient and change the vortex shedding frequencies of flow around tile cylinder.展开更多
In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteris...In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.展开更多
A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently so...A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently solved by the modified incomplete LU decomposition conjugate gradient scheme ( MILU-CG). The effects of surface suction or blowing' s position and strength on the vortex structures in the cylinder wake, as well as on the drag and lift forces at Reynoldes number Re = 100 were investigated numerically. The results show that the suction on the shoulder of the cylinder or the blowing on the rear of the cylinder can effeciently suppress the asymmetry of the vortex wake in the transverse direction and greatly reduce the lift force; the suction on the shoulder of the cylinder, when its strength is properly chosen, can reduce the drag force significantly, too.展开更多
ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged u...ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.展开更多
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.展开更多
This article investigates the separated flow around the bent-nose biconic, in which the second ordered TVD(Total Variation Diminishing) method is applied. The governing cqua-tion is a thin layer Navier-Stokes equation...This article investigates the separated flow around the bent-nose biconic, in which the second ordered TVD(Total Variation Diminishing) method is applied. The governing cqua-tion is a thin layer Navier-Stokes equation. The surface pressure distribution at the leeward side is a little higher in the vicinity of the symmetrical plane. This seems to be caused by the development of the vortex at the leeward side due to the separated flow. The location of a separation line and the size are in very good agreement with the experiment. Also the secondary separation is captured very clearly.展开更多
A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusio...A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.展开更多
Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of th...Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.展开更多
In the traditional machining process for diffusers, blades are easily deformed, and methods suffer from high tool wear and low efficiency. Electrochemical machining(ECM) possesses unique advantages when applied to the...In the traditional machining process for diffusers, blades are easily deformed, and methods suffer from high tool wear and low efficiency. Electrochemical machining(ECM) possesses unique advantages when applied to these difficult-to-machine materials. In the ECM process, theflow field plays a crucial role. Here, an electrolyte flow mode that supplies uniform flow around the entire blade profile was adopted for electrochemical trepanning of diffusers. Various flow rates were employed to obtain the optimal flow field. Simulations were conducted using ANSYS software, and results indicated that increasing the flow rate substantially afforded a more uniform flowfield. A series of experiments was then performed, and results revealed that increasing the flow rate greatly improved both the machining efficiency and the surface quality of the diffusers. The maximum feeding rate of the cathode reached 4 mm/min, the blade taper of the concave part decreased to 0.02, and the blade roughness was reduced to 1.216 lm. The results of this study demonstrated the high feasibility of this method and its potential for machining other complex components for engineering applications.展开更多
Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the co...Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the continuity and momentum equations are satisfied simultaneously at each time step for oscillating flow. The numerical results show that the motion of vortices around a rolling ship hull is cyclical. It is found that the location of the vortices is very similar to the existing experimental result. Using these simulation results, we can calculate the roll damping of ships including viscous effects.展开更多
With the increase of petroleum and gas production in deep ocean, marine risers of circular cylinder shape are widely used in the offshore oil and gas platform. In order to research the hydrodynamic performance of mari...With the increase of petroleum and gas production in deep ocean, marine risers of circular cylinder shape are widely used in the offshore oil and gas platform. In order to research the hydrodynamic performance of marine risers, the dynamic mesh technique and User-Defined Function(UDF) are used to simulate the circular cylinder motion. The motion of a transversely oscillating circular cylinder in combination of uniform flow and oscillating flow is simulated. The uniform flow and oscillating flow both are in x direction. SIMPLE algorithm is used to solve the Navier-Stokes equations. The User-Defined Function is used to control the cylinder transverse vibration and the inlet flow. The lift and drag coefficient changing with time and the map of vorticity isolines at different phase angle are obtained. Force time histories are shown for uniform flow at Reynolds number(Re) of 200 and for the combination of uniform and oscillating flows. With the increase of amplitude of oscillating flow in combined flow, the change of lift amplitude is not sensitive to the the change of cylinder oscillating frequency. Lift amplitude increases with the increase of oscillating flow amplitude in the combined flow, but there is no definite periodicity of the lift coefficient. The drag and inertia force coefficients change when the maximum velocity of the oscillating flow increases in the combined flow. The vortex shedding near the circular cylinder shows different characteristics.展开更多
In order to predict the effects of cavitation on a hydrofoil, the state equations of the cavitation model were combined with a linear viscous turbulent method for mixed fluids in the computational fluid dynamics (CFD)...In order to predict the effects of cavitation on a hydrofoil, the state equations of the cavitation model were combined with a linear viscous turbulent method for mixed fluids in the computational fluid dynamics (CFD) software FLUENT to simulate steady cavitating flow. At a fixed attack angle, pressure distributions and volume fractions of vapor at different cavitation numbers were simulated, and the results on foil sections agreed well with experimental data. In addition, at the various cavitation numbers, the vapor fractions at different attack angles were also predicted. The vapor region moved towards the front of the airfoil and the length of the cavity grew with increased attack angle. The results show that this method of applying FLUENT to simulate cavitation is reliable.展开更多
In this paper, flow around two circular cylinders in tandem arrangement with unequal diameters has been investigated using the particle image velocimetry technique(PIV) in a water channel. The upstream to downstream d...In this paper, flow around two circular cylinders in tandem arrangement with unequal diameters has been investigated using the particle image velocimetry technique(PIV) in a water channel. The upstream to downstream diameter ratio was kept constant at d/D = 2/3, the centre-to-centre distance was varied from 1.2D to 5D and the Reynolds number was varied from 1200 to 4800. The flow characteristics were analyzed through ensemble-averaged patterns of velocity, vorticity, normalized Reynolds stress contours and streamlines. Based on ensemble-averaged and instantaneous flow fields, different flow patterns, including single-wakeshedding at small spacing ratio, bi-stable flow behavior(alternating behavior of reattachment and vortex shedding) at intermediate spacing ratio and co-shedding pattern at large spacing ratio were observed. The effects of Reynolds number and the centre-to-centre spacing ratio on flow patterns and turbulent characteristics were also investigated. It was found that the diameter ratio appears to have a certain effect on the flow patterns at intermediate spacing ratios, where the reattachment of shear layer depends on the lateral width of the wake flow in the lee of the upstream cylinder. Extensive discussion on the distributions of Reynolds stress and turbulent kinetic energy was presented.展开更多
The flexible transmission shaft and wheel propeller are combined as the kinetic source equipment, which realizes the nmlti-motion modes of the autonomous underwater vehicle （AUV） such as vectored thruster and wheele...The flexible transmission shaft and wheel propeller are combined as the kinetic source equipment, which realizes the nmlti-motion modes of the autonomous underwater vehicle （AUV） such as vectored thruster and wheeled movement. In order to study the interactional principle between the hull and the wheel propellers while the AUV navigating in water, the computational fluid dynamics （CFD） method is used to simulate numerically the unsteady viscous flow around AUV with propellers by using the Reynolds-averaged Navier-Stokes （RANS） equations, shear-stress transport （SST） k-w model and pressure with splitting of operators （PISO） algorithm based on sliding mesh. The hydrodynamic parameters of AUV with propellers such as resistance, pressure and velocity are got, which reflect well the real ambient flow field of AUV with propellers. Then, the semi-implicit method for pressure-linked equations （SIMPLE） algorithm is used to compute the steady viscous flow field of AUV hull and propellers, respectively. The computational results agree well with the experimental data, which shows that the numerical method has good accuracy in the prediction of hydrodynamic performance. The interaction between AUV hull and wheel propellers is predicted qualitatively and quantitatively by comparing the hydrodynamic parameters such as resistance, pressure and velocity with those from integral computation and partial computation of the viscous flow around AUV with propellers, which provides an effective reference to the shady on noise and vibration of AUV hull and propellers in real environment. It also provides technical support for the design of new AUVs.展开更多
In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge mod...In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.展开更多
In order to verify the effectiveness and superiority of the dynamic hybrid RANS/LES(DHRL)model,the flow around a cylinder with sinusoidal fluctuating velocity at the inlet was used as the test case.The latest computat...In order to verify the effectiveness and superiority of the dynamic hybrid RANS/LES(DHRL)model,the flow around a cylinder with sinusoidal fluctuating velocity at the inlet was used as the test case.The latest computational fluid dynamics(CFD)model can flexibly choose any existing large-eddy simulation(LES)method combined with RANS method to calculate the flow field.In addition,the DLES model and DDES model are selected as typical representatives of the turbulence model to compare the capture ability of the flow field mechanism.The internal flow field including the y+value,velocity distribution,turbulent kinetic energy and vortex structures is comprehensively analyzed.Finally,the results show that the new model has enough sensitivity to capture the information of the flow field and has more consistent velocity distribution with the experimental value,which shows its potential in practical engineering applications to some extent.展开更多
In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Backgroun...In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Background Oriented Schlieren technique (BOS) allows the measurement of the light deflection caused by density gradients in a compressible flow. For this purpose the distortion of the image of a background pattern observed through the flow is used. In order to increase the performance of the conventional Background Oriented Schlieren technique, the monochromatic background is replaced by a colored dot pattern. The different colors are treated separately using suitable correlation algorithms. Therefore, the precision and the spatial resolution can be highly increased. Furthermore a special arrangement of the different colored dot patterns in the background allows astigmatism in the region with high density gradients to be overcome. For the first time an algebraic reconstruction technique (ART) is then used to reconstruct the density field of unsteady flows around a spike-tipped model from CBOS measurements. The obtained images reveal the interaction between the free-stream flow and the high-pressure region in front of the model, which leads to large-scale instabilities in the flow.展开更多
This paper presents the results of an experimental study on the influences of floodplain impermeable groynes on flow structure, velocity, and water depth around the groyne（s）. A wooden symmetrical compound channel w...This paper presents the results of an experimental study on the influences of floodplain impermeable groynes on flow structure, velocity, and water depth around the groyne（s）. A wooden symmetrical compound channel was used. Groyne models with three different groyne relative lengths, 0.5, 0.75, and 1.0, were used on one floodplain with single and series arrangements. Analysis of the experimental results using the measured flow velocity and water depth values showed that flow structure, velocity, and water depth mainly depend on groyne relative length and the relative distance between series groynes. The flow velocity at the main channel centerline increased by about 40%, 60%, and 85%, and in other parts on the horizontal plane at the floodplain mid-water, depth by about 75%, 125%, and 175% of its original value in eases of one-side floodplain groyne（s） with relative lengths of 0.5, 0.75, and 1.0, respectively. The effective distance between two groynes in series arrangement ranges from 3 to 4 times the groyne length. Using an impermeable groyne with a large relative length in river floodplains increases the generation of eddy and roller zones downstream of the groyne, leading to more scouring and deposition. To avoid that, the groyne relative length must be kept below half the floodplain width,展开更多
A method to evaluate the properties of turbulent flow in proximity to the vehicle and close to the ground surface has been elaborated.Numerical simulations have been performed on the basis of the Unsteady Reynolds-ave...A method to evaluate the properties of turbulent flow in proximity to the vehicle and close to the ground surface has been elaborated.Numerical simulations have been performed on the basis of the Unsteady Reynolds-averaged Navier-Stokes equations(URANS)written with respect to an arbitrary curvilinear coordinate system.These equations have been solved using the Spalart-Allmaras differential one-parametric turbulence model.The method of artificial compressibility has been used to improve the coupling of pressure and velocity in the framework of a finite volume approach.Time-averaged distributions of pressure fields,velocity components,streamlines in the entire area and near the tractor-trailer,as well as integral and distributed characteristic parameters(such as coefficients of pressure,friction and drag force)are presented.According to our results,the turbulent flow accelerates in the area of the tractor cabin and in the gap between surfaces.Above the driver’s cabin,a pressure drop occurs due to a sharp acceleration of flow in this area.Downstream,pressure is restored and becomes almost constant in proximity to the edge of the trailer.The dimensions of the separation area exceed the length of the transport system several times.Though agreement with experimental results is relatively limited due to the two-dimensional nature of the numerical simulations,the present approach succeeds in identifying the main physical effects involved in the considered dynamics.It might be used in future studies for initial approximate assessments of the influence of the vehicle shape on its aerodynamic characteristics.展开更多
基金Supported by the National Key R&D Program of China(2022YFB3303501)the National Natural Science Foundation of China(Project Nos.52176041 and 12102308)the Fundamental Research Funds for the Central Universities(Project Nos.2042023kf0208 and 2042023kf0159).
文摘Verification and validation(V&V)is a helpful tool for evaluating simulation errors,but its application in unsteady cavitating flow remains a challenging issue due to the difficulty in meeting the requirement of an asymptotic range.Hence,a new V&V approach for large eddy simulation(LES)is proposed.This approach offers a viable solution for the error estimation of simulation data that are unable to satisfy the asymptotic range.The simulation errors of cavitating flow around a projectile near the free surface are assessed using the new V&V method.The evident error values are primarily dispersed around the cavity region and free surface.The increasingly intense cavitating flow increases the error magnitudes.In addition,the modeling error magnitudes of the Dynamic Smagorinsky-Lilly model are substantially smaller than that of the Smagorinsky-Lilly model.The present V&V method can capture the decrease in the modeling errors due to model enhancements,further exhibiting its applicability in cavitating flow simulations.Moreover,the monitoring points where the simulation data are beyond the asymptotic range are primarily dispersed near the cavity region,and the number of such points grows as the cavitating flow intensifies.The simulation outcomes also suggest that the re-entrant jet and shedding cavity collapse are the chief sources of vorticity motions,which remarkably affect the simulation accuracy.The results of this study provide a valuable reference for V&V research.
文摘The electric-static body force model is obtained by solving Maxwell＇s electromagnetic equations. Based on the electro-static model, numerical modeling of flow around a cylinder with a dielectric barrier discharge （DBD） plasma effect is also presented. The flow streamlines between the numerical simulation and the particle image velocimetry （PIV） experiment are consistent. According to the numerical simulation, DBD plasma can reduce the drag coefficient and change the vortex shedding frequencies of flow around tile cylinder.
基金supported by the National High-Tec Research and Development Program of China(2006AA09A104)
文摘In order to investigate the scale effect of turbulent flow around a circular cylinder, two similarity numbers （criteria） based on turbulent kinetic and dissipation rates associ- ated with the fluctuation characteristics of turbulence wake are deduced by analyzing the Reynolds averaged NavierStokes equations （RANS）. The RNG k-s models and finite volume method are used to solve the governing equations and the second-order implicit time and upwind space discretization algorithms are used to discrete the governing equations. A numerical computation of flow parameters around a two-dimensional circular cylinder with Reynolds numbers ranging from 102 to l07 is accomplished and the result indicates that the fluctuation of turbulence flow along the center line in the wake of circular cylinder can never be changed with increasing Reynolds numbers when Re ≥ 3 × 10^6. This conclusion is useful for controlling the scale of numerical calculations and for applying model test data to engineering practice.
基金Foundation item:the Natural Science Foundation of Jiangsu Province(BK97056109)
文摘A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently solved by the modified incomplete LU decomposition conjugate gradient scheme ( MILU-CG). The effects of surface suction or blowing' s position and strength on the vortex structures in the cylinder wake, as well as on the drag and lift forces at Reynoldes number Re = 100 were investigated numerically. The results show that the suction on the shoulder of the cylinder or the blowing on the rear of the cylinder can effeciently suppress the asymmetry of the vortex wake in the transverse direction and greatly reduce the lift force; the suction on the shoulder of the cylinder, when its strength is properly chosen, can reduce the drag force significantly, too.
基金Project(2006BAJ04B04)supported by the National Science and Technology Pillar Program in the Eleventh Five-year Plan PeriodProject(2006AA05Z229)supported by the National High Technology Research and Development Program of China+1 种基金Project supportedby the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education MinistryProject(06wk3023)supported by Hunan Science and Technology Office
文摘ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.
文摘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.
文摘This article investigates the separated flow around the bent-nose biconic, in which the second ordered TVD(Total Variation Diminishing) method is applied. The governing cqua-tion is a thin layer Navier-Stokes equation. The surface pressure distribution at the leeward side is a little higher in the vicinity of the symmetrical plane. This seems to be caused by the development of the vortex at the leeward side due to the separated flow. The location of a separation line and the size are in very good agreement with the experiment. Also the secondary separation is captured very clearly.
文摘A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.
基金China Postdoctoral Science Foundation (20100481368)National Key Laboratory Foundation of China
文摘Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.
基金co-supported by the National Natural Science Foundation of China(51675271)the Natural Science Foundation of Jiangsu Province(BK20171413)+1 种基金the Fundamental Research Funds for the Central Universities(NE 2017003)the Qing Lan Project
文摘In the traditional machining process for diffusers, blades are easily deformed, and methods suffer from high tool wear and low efficiency. Electrochemical machining(ECM) possesses unique advantages when applied to these difficult-to-machine materials. In the ECM process, theflow field plays a crucial role. Here, an electrolyte flow mode that supplies uniform flow around the entire blade profile was adopted for electrochemical trepanning of diffusers. Various flow rates were employed to obtain the optimal flow field. Simulations were conducted using ANSYS software, and results indicated that increasing the flow rate substantially afforded a more uniform flowfield. A series of experiments was then performed, and results revealed that increasing the flow rate greatly improved both the machining efficiency and the surface quality of the diffusers. The maximum feeding rate of the cathode reached 4 mm/min, the blade taper of the concave part decreased to 0.02, and the blade roughness was reduced to 1.216 lm. The results of this study demonstrated the high feasibility of this method and its potential for machining other complex components for engineering applications.
基金This project is supported by Fok Ying Tung Education Foundation and National Natural Science Foundation of China
文摘Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the continuity and momentum equations are satisfied simultaneously at each time step for oscillating flow. The numerical results show that the motion of vortices around a rolling ship hull is cyclical. It is found that the location of the vortices is very similar to the existing experimental result. Using these simulation results, we can calculate the roll damping of ships including viscous effects.
基金supported financially by the Natural Science Foundation of China (No. 51079136/51179179/51239008)
文摘With the increase of petroleum and gas production in deep ocean, marine risers of circular cylinder shape are widely used in the offshore oil and gas platform. In order to research the hydrodynamic performance of marine risers, the dynamic mesh technique and User-Defined Function(UDF) are used to simulate the circular cylinder motion. The motion of a transversely oscillating circular cylinder in combination of uniform flow and oscillating flow is simulated. The uniform flow and oscillating flow both are in x direction. SIMPLE algorithm is used to solve the Navier-Stokes equations. The User-Defined Function is used to control the cylinder transverse vibration and the inlet flow. The lift and drag coefficient changing with time and the map of vorticity isolines at different phase angle are obtained. Force time histories are shown for uniform flow at Reynolds number(Re) of 200 and for the combination of uniform and oscillating flows. With the increase of amplitude of oscillating flow in combined flow, the change of lift amplitude is not sensitive to the the change of cylinder oscillating frequency. Lift amplitude increases with the increase of oscillating flow amplitude in the combined flow, but there is no definite periodicity of the lift coefficient. The drag and inertia force coefficients change when the maximum velocity of the oscillating flow increases in the combined flow. The vortex shedding near the circular cylinder shows different characteristics.
文摘In order to predict the effects of cavitation on a hydrofoil, the state equations of the cavitation model were combined with a linear viscous turbulent method for mixed fluids in the computational fluid dynamics (CFD) software FLUENT to simulate steady cavitating flow. At a fixed attack angle, pressure distributions and volume fractions of vapor at different cavitation numbers were simulated, and the results on foil sections agreed well with experimental data. In addition, at the various cavitation numbers, the vapor fractions at different attack angles were also predicted. The vapor region moved towards the front of the airfoil and the length of the cavity grew with increased attack angle. The results show that this method of applying FLUENT to simulate cavitation is reliable.
基金supported by the Zhejiang Provincial Natural Science Foundation of China under Grant No. LY14E090009State Key Laboratory of Satellite Ocean Environment Dynamics (Second Institute of Oceanography, SOA), State Key Laboratory of Fluid Power Transmission and Control (GZKF-201310)+1 种基金State Key Laboratory of Ocean Engineering, China. The National Research Foundation of Singapore (NRF-CRP5-2009-01)Maritime Research Centre and Nanyang Environment and Water Research Institute, Nanyang Technological University, Singapore, is acknowledged
文摘In this paper, flow around two circular cylinders in tandem arrangement with unequal diameters has been investigated using the particle image velocimetry technique(PIV) in a water channel. The upstream to downstream diameter ratio was kept constant at d/D = 2/3, the centre-to-centre distance was varied from 1.2D to 5D and the Reynolds number was varied from 1200 to 4800. The flow characteristics were analyzed through ensemble-averaged patterns of velocity, vorticity, normalized Reynolds stress contours and streamlines. Based on ensemble-averaged and instantaneous flow fields, different flow patterns, including single-wakeshedding at small spacing ratio, bi-stable flow behavior(alternating behavior of reattachment and vortex shedding) at intermediate spacing ratio and co-shedding pattern at large spacing ratio were observed. The effects of Reynolds number and the centre-to-centre spacing ratio on flow patterns and turbulent characteristics were also investigated. It was found that the diameter ratio appears to have a certain effect on the flow patterns at intermediate spacing ratios, where the reattachment of shear layer depends on the lateral width of the wake flow in the lee of the upstream cylinder. Extensive discussion on the distributions of Reynolds stress and turbulent kinetic energy was presented.
基金Project(2006AA09Z235) supported by National High Technology Research and Development Program of ChinaProject(CX2009B003) supported by Hunan Provincial Innovation Foundation For Postgraduate,China
文摘The flexible transmission shaft and wheel propeller are combined as the kinetic source equipment, which realizes the nmlti-motion modes of the autonomous underwater vehicle （AUV） such as vectored thruster and wheeled movement. In order to study the interactional principle between the hull and the wheel propellers while the AUV navigating in water, the computational fluid dynamics （CFD） method is used to simulate numerically the unsteady viscous flow around AUV with propellers by using the Reynolds-averaged Navier-Stokes （RANS） equations, shear-stress transport （SST） k-w model and pressure with splitting of operators （PISO） algorithm based on sliding mesh. The hydrodynamic parameters of AUV with propellers such as resistance, pressure and velocity are got, which reflect well the real ambient flow field of AUV with propellers. Then, the semi-implicit method for pressure-linked equations （SIMPLE） algorithm is used to compute the steady viscous flow field of AUV hull and propellers, respectively. The computational results agree well with the experimental data, which shows that the numerical method has good accuracy in the prediction of hydrodynamic performance. The interaction between AUV hull and wheel propellers is predicted qualitatively and quantitatively by comparing the hydrodynamic parameters such as resistance, pressure and velocity with those from integral computation and partial computation of the viscous flow around AUV with propellers, which provides an effective reference to the shady on noise and vibration of AUV hull and propellers in real environment. It also provides technical support for the design of new AUVs.
基金Project(U1534210)supported by the National Natural Science Foundation of ChinaProject(14JJ1003)supported by the Natural Science Foundation of Hunan Province,China+2 种基金Project(2015CX003)supported by the Project of Innovation-driven Plan in Central South University,ChinaProject(14JC1003)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2015T002-A)supported by the Technological Research and Development program of China Railways Cooperation
文摘In order to describe an investigation of the flow around high-speed train on a bridge under cross winds using detached-eddy simulation(DES), a 1/8th scale model of a three-car high-speed train and a typical bridge model are employed, Numerical wind tunnel technology based on computational fluid dynamics(CFD) is used, and the CFD models are set as stationary models. The Reynolds number of the flow, based on the inflow velocity and the height of the vehicle, is 1.9×10~6. The computations are conducted under three cases, train on the windward track on the bridge(WWC), train on the leeward track on the bridge(LWC) and train on the flat ground(FGC). Commercial software FLUENT is used and the mesh sensitivity research is carried out by three different grids: coarse, medium and fine. Results show that compared with FGC case, the side force coefficients of the head cars for the WWC and LWC cases increases by 14% and 29%, respectively; the coefficients of middle cars for the WWC and LWC increase by 32% and 10%, respectively; and that of the tail car increases by 45% for the WWC whereas decreases by 2% for the LWC case. The most notable thing is that the side force and the rolling moment of the head car are greater for the LWC, while the side force and the rolling moment of the middle car and the tail car are greater for the WWC. Comparing the velocity profiles at different locations, the flow is significantly influenced by the bridge-train system when the air is close to it. For the three cases(WWC, LWC and FGC), the pressure on the windward side of train is mostly positive while that of the leeward side is negative. The discrepancy of train's aerodynamic force is due to the different surface area of positive pressure and negative pressure zone. Many vortices are born on the leeward edge of the roofs. Theses vortices develop downstream, detach and dissipate into the wake region. The eddies develop irregularly, leading to a noticeably turbulent flow at leeward side of train.
基金Supported by the Open Fund of Key Laboratory of Road Construction Technology and Equipment of Chang’an University,Ministry of Education(310825171104)the Advanced Manufacturing Projects of Government and University Co-construction Program Funded by Jilin Province(SXGJSF2017-2)
文摘In order to verify the effectiveness and superiority of the dynamic hybrid RANS/LES(DHRL)model,the flow around a cylinder with sinusoidal fluctuating velocity at the inlet was used as the test case.The latest computational fluid dynamics(CFD)model can flexibly choose any existing large-eddy simulation(LES)method combined with RANS method to calculate the flow field.In addition,the DLES model and DDES model are selected as typical representatives of the turbulence model to compare the capture ability of the flow field mechanism.The internal flow field including the y+value,velocity distribution,turbulent kinetic energy and vortex structures is comprehensively analyzed.Finally,the results show that the new model has enough sensitivity to capture the information of the flow field and has more consistent velocity distribution with the experimental value,which shows its potential in practical engineering applications to some extent.
文摘In this paper, the improved Background Oriented Schlieren technique called CBOS (Colored Background Oriented Schlieren) is described and used to reconstruct the density fields of three-dimensional flows. The Background Oriented Schlieren technique (BOS) allows the measurement of the light deflection caused by density gradients in a compressible flow. For this purpose the distortion of the image of a background pattern observed through the flow is used. In order to increase the performance of the conventional Background Oriented Schlieren technique, the monochromatic background is replaced by a colored dot pattern. The different colors are treated separately using suitable correlation algorithms. Therefore, the precision and the spatial resolution can be highly increased. Furthermore a special arrangement of the different colored dot patterns in the background allows astigmatism in the region with high density gradients to be overcome. For the first time an algebraic reconstruction technique (ART) is then used to reconstruct the density field of unsteady flows around a spike-tipped model from CBOS measurements. The obtained images reveal the interaction between the free-stream flow and the high-pressure region in front of the model, which leads to large-scale instabilities in the flow.
文摘This paper presents the results of an experimental study on the influences of floodplain impermeable groynes on flow structure, velocity, and water depth around the groyne（s）. A wooden symmetrical compound channel was used. Groyne models with three different groyne relative lengths, 0.5, 0.75, and 1.0, were used on one floodplain with single and series arrangements. Analysis of the experimental results using the measured flow velocity and water depth values showed that flow structure, velocity, and water depth mainly depend on groyne relative length and the relative distance between series groynes. The flow velocity at the main channel centerline increased by about 40%, 60%, and 85%, and in other parts on the horizontal plane at the floodplain mid-water, depth by about 75%, 125%, and 175% of its original value in eases of one-side floodplain groyne（s） with relative lengths of 0.5, 0.75, and 1.0, respectively. The effective distance between two groynes in series arrangement ranges from 3 to 4 times the groyne length. Using an impermeable groyne with a large relative length in river floodplains increases the generation of eddy and roller zones downstream of the groyne, leading to more scouring and deposition. To avoid that, the groyne relative length must be kept below half the floodplain width,
文摘A method to evaluate the properties of turbulent flow in proximity to the vehicle and close to the ground surface has been elaborated.Numerical simulations have been performed on the basis of the Unsteady Reynolds-averaged Navier-Stokes equations(URANS)written with respect to an arbitrary curvilinear coordinate system.These equations have been solved using the Spalart-Allmaras differential one-parametric turbulence model.The method of artificial compressibility has been used to improve the coupling of pressure and velocity in the framework of a finite volume approach.Time-averaged distributions of pressure fields,velocity components,streamlines in the entire area and near the tractor-trailer,as well as integral and distributed characteristic parameters(such as coefficients of pressure,friction and drag force)are presented.According to our results,the turbulent flow accelerates in the area of the tractor cabin and in the gap between surfaces.Above the driver’s cabin,a pressure drop occurs due to a sharp acceleration of flow in this area.Downstream,pressure is restored and becomes almost constant in proximity to the edge of the trailer.The dimensions of the separation area exceed the length of the transport system several times.Though agreement with experimental results is relatively limited due to the two-dimensional nature of the numerical simulations,the present approach succeeds in identifying the main physical effects involved in the considered dynamics.It might be used in future studies for initial approximate assessments of the influence of the vehicle shape on its aerodynamic characteristics.
