Cavitation generation methods have been used in multifarious directions because of their diversity,and numerous studies and discussions have been conducted on cavitation generation methods.This study aims to explore t...Cavitation generation methods have been used in multifarious directions because of their diversity,and numerous studies and discussions have been conducted on cavitation generation methods.This study aims to explore the generating mechanism and evolution law of volume alternate cavitation(VAC).In the VAC,liquid water is placed in an airtight container with a variable volume.As the volume alternately changes,the liquid water inside the container continues to cavitate.Then,the mixture turbulence model and in-cylinder dynamic grid model are adopted to conduct computational fluid dynamics simulation of volume alternate cavitation.In the simulation,the cloud images at seven heights on the central axis are monitored,and the phenomenon and mechanism of height and eccentricity are analyzed in detail.By employing the cavitation flow visualization method,the generating mechanism and evolution law of cavitation are revealed.The synergistic effects of experiments and high-speed camera capturing confirm the correctness of the simulation results.In the experiment,the volume change stroke of the airtight container is set to 20 mm,the volume change frequency is 18 Hz,and the shooting frequency of the high-speed camera is set to 10000 FPS.The experimental results indicate that the position of the cavitation phenomenon has a reasonable law during the whole evolution cycle of the cavitation cloud.Also,the volume alternation cycle corresponds to the generation,development,and collapse stages of cavitation bubbles.展开更多
The cavitating flow around a Delft Twist-11 hydrofoil is simulated using the large eddy simulation approach.The volume-of-fluid method incorporated with the Schnerr-Sauer cavitation model is utilized to track the wate...The cavitating flow around a Delft Twist-11 hydrofoil is simulated using the large eddy simulation approach.The volume-of-fluid method incorporated with the Schnerr-Sauer cavitation model is utilized to track the water-vapor interface.Adaptive mesh refinement(AMR)is also applied to improve the simulation accuracy automatically.Two refinement levels are conducted to verify the dominance of AMR in predicting cavitating flows.Results show that cavitation features,including the U-type structure of shedding clouds,are consistent with experimental observations.Even a coarse mesh can precisely capture the phase field without increasing the total cell number significantly using mesh adaption.The predicted shedding frequency agrees fairly well with the experimental data under refinement level 2.This study illustrates that AMR is a promising approach to achieve accurate simulations for multiscale cavitating flows within limited computational costs.Finally,the force element method is currently adopted to investigate the lift and drag fluctuations during the evolution of cavitation structure.The mechanisms of lift and drag fluctuations due to cavitation and the interaction between vorticity forces and cavitation are explicitly revealed.展开更多
The aim of this study is to suppress a cavitation near the orifice of arelief valve by changing the shape of a poppet. An experimental flow visualization technique and anumerical cavitating flow simulation, using a RN...The aim of this study is to suppress a cavitation near the orifice of arelief valve by changing the shape of a poppet. An experimental flow visualization technique and anumerical cavitating flow simulation, using a RNG kappa-epsilon turbulence model and a cavitationmodel, are employed to achieve the purpose. In the flow visualization, the cavitation phenomenonnear the orifice of a relief valve is observed using a transparent test valve body model and acamera. On the other hand, a three dimensional cavitating flow simulation is conducted to predictthe cavitation near the orifice of a relief valve. Six types of poppets are designed by changing theshape of a traditional poppet shape, which is expected to influence the cavitating flow in anorifice. In addition, the cavitation noise of a relief valve is measured and the noise spectrum isanalyzed. In conclusion, the cavitation intensity and the cavitation noise are reduced for anoptimal poppet geometry obtained in the present study.展开更多
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.展开更多
Computational Fluid Dynamics (CFD) simulations of cavitating flow through water hydraulic poppet valves were performed using advanced RNG k-epsilon turbulence model. The flow was turbulent, incompressible and unsteady...Computational Fluid Dynamics (CFD) simulations of cavitating flow through water hydraulic poppet valves were performed using advanced RNG k-epsilon turbulence model. The flow was turbulent, incompressible and unsteady, for Reynolds numbers greater than 43 000. The working fluid was water, and the structure of the valve was simplified as a two dimensional axisymmetric geometrical model. Flow field visualization was numerically achieved. The effects of inlet velocity, outlet pressure, opening size as well as poppet angle on cavitation intensity in the poppet valve were numerically investigated. Experimental flow visualization was conducted to capture cavitation images near the orifice in the poppet valve with 30° poppet angle using high speed video camera. The binary cavitating flow field distribution obtained from digital processing of the original cavitation image showed a good agreement with the numerical result.展开更多
The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resour...The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.展开更多
The cavitating flow in different regimes has the intricate flow structure with multiple time and space scales.The present work develops a multiscale model by coupling the volume of fluid(VOF)method and a discrete bubb...The cavitating flow in different regimes has the intricate flow structure with multiple time and space scales.The present work develops a multiscale model by coupling the volume of fluid(VOF)method and a discrete bubble model(DBM),to simulate the cavitating flow in a convergent-divergent test section.The Schnerr-Sauer cavitation model is used to calculate the mass transfer rate to obtain the macroscale phase structure,and the simplified Rayleigh-Plesset equation is applied to simulate the growing and collapsing of discrete bubbles.An algorithm for bridging between the macroscale cavities and microscale bubbles is also developed to achieve the multiscale simulation.For the flow field,the very large eddy simulation(VLES)approach is applied.Conditions from inception to sheet/cloud cavitation regimes are taken into account and simulations are conducted.Compared with the experimental observations,it is shown that the cavitation inception,bubble clouds formation and glass cavity generation are all well represented,indicating that the proposed VOF-DBM model is a promising approach to accurately and comprehensively reveal the multiscale phase field induced by cavitation.展开更多
Based on a self-developed hydrodynamic cavitation device with different geometric parameters for circular multi-orifice plates,turbulence characteristics of cavitating flow behind multi-orifice plates,including the ef...Based on a self-developed hydrodynamic cavitation device with different geometric parameters for circular multi-orifice plates,turbulence characteristics of cavitating flow behind multi-orifice plates,including the effects of orifice number and orifice layout on longitudinal velocity,turbulence intensity,and Reynolds stress,were measured with the particle image velocimetry(PIV)technique.Flow regimes of the cavitating flow were also observed with high-speed photography.The experimental results showed the following:(1)high-velocity multiple cavitating jets occurred behind the multi-orifice plates,and the cavitating flow fields were characterized by topological structures;(2)the longitudinal velocity at each cross-section exhibited a sawtooth-like distribution close to the multi-orifice plate,and each sawtooth indicated one jet issuing from one orifice;(3)there were similar magnitudes and forms for the longitudinal and vertical turbulence intensities at the same cross-section;(4)the variation in amplitude of Reynolds stress increased with an increase in orifice number;and(5)the cavitation clouds in the flow fields became denser with the increase in orifice number,and the clouds generated by the staggered layout of orifices were greater in number than those generated by the checkerboard-type one for the same orifice number.The experimental results can be used to analyze the mechanism of killing pathogenic microorganisms through hydrodynamic cavitation.展开更多
Nowadays,most researchers focus on the cavity shedding mechanisms of unsteady cavitating flows over different objects,such as 2D/3D hydrofoils,venturi-type section,axisymmetric bodies with different headforms,and so o...Nowadays,most researchers focus on the cavity shedding mechanisms of unsteady cavitating flows over different objects,such as 2D/3D hydrofoils,venturi-type section,axisymmetric bodies with different headforms,and so on.But few of them pay attention to the differences of cavity shedding modality under different cavitation numbers in unsteady cavitating flows over the same object.In the present study,two kinds of shedding patterns are investigated experimentally.A high speed camera system is used to observe the cavitating flows over an axisymmetric blunt body and the velocity fields are measured by a particle image velocimetry(PIV)technique in a water tunnel for different cavitation conditions.The U-type cavitating vortex shedding is observed in unsteady cavitating flows.When the cavitation number is 0.7,there is a large scale cavity rolling up and shedding,which cause the instability and dramatic fluctuation of the flows,while at cavitation number of 0.6,the detached cavities can be conjunct with the attached part to induce the break-off behavior again at the tail of the attached cavity,as a result,the final shedding is in the form of small scale cavity and keeps a relatively steady flow field.It is also found that the interaction between the re-entrant flow and the attached cavity plays an important role in the unsteady cavity shedding modality.When the attached cavity scale is insufficient to overcome the re-entrant flow,it deserves the large cavity rolling up and shedding just as that at cavitation number of 0.7.Otherwise,the re-entrant flow is defeated by large enough cavity to induce the cavity-combined process and small scale cavity vortexes shedding just as that of the cavitation number of0.6.This research shows the details of two different cavity shedding modalities which is worthful and meaningful for the further study of unsteady cavitation.展开更多
The water-lubricated bearings are usually the state of turbulent cavitating flow under high-speed conditions. And the distribution of cavitation bubbles and the interface effect between the two phases have not been in...The water-lubricated bearings are usually the state of turbulent cavitating flow under high-speed conditions. And the distribution of cavitation bubbles and the interface effect between the two phases have not been included in previous studies on high-speed water-lubricated bearings. In order to study the influence of interface effect and cavitation bubble distribution on the dynamic characteristics of high-speed water-lubricated spiral groove thrust bearings(SGTB).A turbulent cavitating flow lubrication model based on two-phase fluid and population balance equation of bubbles was established in this paper. Stiffness and the damping coefficients of the SGTB were calculated using the perturbation pressure equations. An experimental apparatus was developed to verify the theoretical model. Simulating and experimental results show that the small-sized bubbles tend to generate in the turbulent cavitating flow when at a high rotary speed, and the bubbles mainly locate at the edges of the spiral groove. The simulating results also show that the direct stiffness coefficients are increased due to cavitation effect, and cross stiffness coefficients and damping coefficients are hardly affected by the cavitation effect. Turbulent effect on the dynamic characteristics of SGTB is much stronger than the cavitating effect.展开更多
Quasi-periodical evolutions such as shedding and collapsing of unsteady cloud cavitating flow, induce strong pressure fluctuations, what may deteriorate maneuvering stability and corrode surfaces of underwater vehicle...Quasi-periodical evolutions such as shedding and collapsing of unsteady cloud cavitating flow, induce strong pressure fluctuations, what may deteriorate maneuvering stability and corrode surfaces of underwater vehicles. This paper analyzed effects on cavitation stability of a trip bar arranged on high-speed underwater projectile. Small scale water tank experiment and large eddy simulation using the open source software Open FOAM were used, and the results agree well with each other. Results also indicate that trip bar can obstruct downstream re-entrant jet and pressure wave propagation caused by collapse, resulting in a relatively stable sheet cavity between trip bar and shoulder of projectiles.展开更多
The transient cavitating flow around the Clark-Y hydrofoil is numerically investigated by the dynamic mode decomposition with criterion.Based on the ranking dominant modes,frequencies of the first four modes are in go...The transient cavitating flow around the Clark-Y hydrofoil is numerically investigated by the dynamic mode decomposition with criterion.Based on the ranking dominant modes,frequencies of the first four modes are in good accordance with those obtained by fast Fourier transform.Furthermore,the cavitating flow field is reconstructed by the first four modes,and the dominant flow features are well captured with the reconstructed error below 12%when compared to the simulated flow field.This paper offers a reference for observing and reconstructing the flow fields,and gives a novel insight into the transient cavitating flow features.展开更多
Cavitation within the tip vortex(TV)flow remains a challenging issue in the design of high-speed and low-noise hydraulic machinery.In this paper,the TV cavitating flow around an elliptical hydrofoil is calculated by u...Cavitation within the tip vortex(TV)flow remains a challenging issue in the design of high-speed and low-noise hydraulic machinery.In this paper,the TV cavitating flow around an elliptical hydrofoil is calculated by using large eddy simulation(LES)combined with a modified Schnerr-Sauer(S-S)cavitation model.The original S-S cavitation model is modified by taking into account the typical effect of vortex flow.The partial pressure term which can describe the vortex quantitatively and qualitatively is confirmed asρ_(m)ω^(2) x r _(c)^(2) ,and is considered into the R-P equation of the modified S-S cavitation model.Comparison between the numerical and experimental results shows good agreement in the form and evolution of cavities,including attached cavities(AC)and tip vortex cavities(TVC).The vorticity transport equation is utilized to investigate the dynamic mechanisms of the vortex development around the TVC.Further analyses indicate that cavitation in the TV flow influences the pressure in the core of the cavity and the local flow patterns.Typical vortex structures in the TV cavitating flow include TV,secondary vortex(SV)and wake vortex(WV).The direction and magnitude of the rotation effect can be described by axial vorticity which is drawn on the iso-surface of Q=1×105 s−2.The development of the TV cavitating flow can be divided into two stages:Stage I,the development and fusion of TV,SV,stage II,the dissipation of SV.The stretching term dominates the evolution of TV,and the dilatation term is the main reason in the mergence process of SV.展开更多
In this paper,after the successful applications to open water propeller performance estimations,the influence of transition sensitive and modified mass transfer models tuned to account for the laminar flow in the pred...In this paper,after the successful applications to open water propeller performance estimations,the influence of transition sensitive and modified mass transfer models tuned to account for the laminar flow in the prediction of the cavitation inception of marine propulsors is investigated from the point of view of the unsteady functioning and induced pressure pulses.The VP1304(also known as PPTC)test case,for which dedicated data were collected during several workshops,is considered first.After preliminary analyses using RANS,also Detached Eddy Simulations(DES)are included to better account for the vortex dynamics and its influence on pressure pulses.Similarly to what observed in uniform inflow,results show a better agreement with the available measurements of propeller performances and confirm the reliability of the proposed approaches for unsteady,non-cavitating,model scale propeller predictions.The overall cavitation pattern is improved too by the application of the transition sensitive correction to the mass transfer model,but the complex dynamics of bubble cavitation observed in experiments prevents quantitatively better predictions in terms of thrust/torque breakdown and induced pressure pulses levels regardless the use of RANS or DES methods.展开更多
In the present paper,the unsteady cavitating turbulent flow over the twisted NACA66 hydrofoil is investigated based on an modified shear stress transfer k-ωpartially averaged Navier-Stokes(MSST PANS)model,i.e.,new MS...In the present paper,the unsteady cavitating turbulent flow over the twisted NACA66 hydrofoil is investigated based on an modified shear stress transfer k-ωpartially averaged Navier-Stokes(MSST PANS)model,i.e.,new MSST PANS(NMSST PANS)model,where the production term of kinetic energy in the turbulence model is modified with helicity.Compared with the experimental data,cavitation evolution and its characteristic frequency are satisfactorily predicted by the proposed NMSST PANS model.It is revealed that the interaction among the main flow,the reentrant jets,and sheet cavitation causes the formation of the primary shedding cavity near the mid-span and the secondary shedding cavity at each side of the twisted hydrofoil,and further induces the remarkable pressure gradient around shedding cavities.Along with the development of the primary and the secondary shedding cavities,the great pressure gradient associated with large cavity volume variation promotes the vortical flow generation and the spatial deformation of vortex structure during cavitation evolution,and results in the primary and the secondary U-type vortices.Further,dynamic mode decomposition(DMD)analysis is utilized to confirm the interaction among the main flow,the main reentrant jet and two side reentrant jets,and cavitation.These results indicate that the proposed NMSST PANS model is suitable to simulate the complicated cavitating turbulent flow for various engineering applications.展开更多
In this paper, by introducing the flow velocity item into the classical Rayleigh-Plesset dynamic equation, a newequation, which does not involve the time term and can describe the motion of cavitation bubble in the st...In this paper, by introducing the flow velocity item into the classical Rayleigh-Plesset dynamic equation, a newequation, which does not involve the time term and can describe the motion of cavitation bubble in the steadycavitating flow, has been obtained. By solving the new motion equation using Runge-Kutta fourth order methodwith adaptive step size control, the dynamic behaviors of cavitation bubble driven by the varying pressure fielddownstream of a venturi cavitation reactor are numerically simulated. The effects of liquid temperature (correspondingto the saturated vapor pressure of liquid), cavitation number and inlet pressure of venturi on radial motionof bubble and pressure pulse due to the radial motion are analyzed and discussed in detail. Some dynamicbehaviors of bubble different from those in previous papers are displayed. In addition, the internal relationshipbetween bubble dynamics and process intensification is also discussed. The simulation results reported in thiswork reveal the variation laws of cavitation intensity with the flow conditions of liquid, and will lay a foundationfor the practical application of hydrodynamic cavitation technology.展开更多
Developing a robust computational strategy to address the rich physical characteristic involved in the thermcdynamic effects on the cryogenic cavitation remains a challenge in research. The objective of the present st...Developing a robust computational strategy to address the rich physical characteristic involved in the thermcdynamic effects on the cryogenic cavitation remains a challenge in research. The objective of the present study is to focus on developing mod- elling strategy to simulate cavitating flows in liquid nitrogen. For this purpose, numerical simulation over a 2D quarter caliber hydrofoil is investigated by calibrating cavitation model parameters and implementing the thermodynamic effects to the Zwart cavitation model. Experimental measurements of pressure and temperature are utilized to validate the extensional Zwart cavi- tation model. The results show that the cavitation dynamics characteristic under the cryogenic environment ale different from that under the isothermal conditions: the cryogenic case yields a substantially shorter cavity around the hydrofoil, and the pre- dicted pressure and temperature inside the cavity are steeper under the cryogenic conditions. Compared with the experimental data, the computational predictions with the modified evaporation and condensation parameters display better results than the default parameters from the room temperature liquids. Based on a wide range of computations and comparisons, the extension- al Zwart cavitation model may predict more accurately the quasi-steady cavitation over a hydrofoil in liquid nitrogen by pri- marily altering the evaporation rate near the leading edge and the condensation rate in the cavity closure region.展开更多
Successful modeling and/or design of engineering systems often requires one to address the impact of multiple "design variables" on the prescribed outcome.There are often multiple,competing objectives based on which...Successful modeling and/or design of engineering systems often requires one to address the impact of multiple "design variables" on the prescribed outcome.There are often multiple,competing objectives based on which we assess the outcome of optimization.Since accurate,high fidelity models are typically time consuming and computationally expensive,comprehensive evaluations can be conducted only if an efficient framework is available.Furthermore,informed decisions of the model/hardware's overall performance rely on an adequate understanding of the global,not local,sensitivity of the individual design variables on the objectives.The surrogate-based approach,which involves approximating the objectives as continuous functions of design variables from limited data,offers a rational framework to reduce the number of important input variables,i.e.,the dimension of a design or modeling space.In this paper,we review the fundamental issues that arise in surrogate-based analysis and optimization,highlighting concepts,methods,techniques,as well as modeling implications for mechanics problems.To aid the discussions of the issues involved,we summarize recent efforts in investigating cryogenic cavitating flows,active flow control based on dielectric barrier discharge concepts,and lithium(Li)-ion batteries.It is also stressed that many multi-scale mechanics problems can naturally benefit from the surrogate approach for "scale bridging."展开更多
Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated ac...Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated accuracy has been evaluated using the LES verification and validation method.Various cavitation numbers are simulated to study important flow characteristics in the sphere wake,e.g.periodic cavity growth/contraction,interactions between the cloud and sheet cavitations and the vortex structure evolution.The spectral characteristics of the wake for typical cloud cavitation conditions were classified as the periodic cavitation mode,high Strouhal number mode and low Strouhal number mode.Main frequency distributions in the wake were analyzed and different dominant flow structures were identified for each of the three modes.Further,the cavitation and vortex relationship was also studied,which is an important issue associated with complex cavitating sphere wakes.Three types of cavitating vortex structures alternate,which indicates that three different cavity shedding regimes may exist in the wake.Analysis of vorticity transport equation shows a significant vorticity increase at the cavitation closure region and in the vortex cavitation region.This study provides a physical perspective to further understand the flow mechanisms in cavitating sphere wakes.展开更多
The paper proposes a methodology within the Reynolds averaged Navier Stokes(RANS) solvers for cavitating flows capable of predicting the flow regions of bubble collapse and the potential aggressiveness to material d...The paper proposes a methodology within the Reynolds averaged Navier Stokes(RANS) solvers for cavitating flows capable of predicting the flow regions of bubble collapse and the potential aggressiveness to material damage. An aggressiveness index is introduced, called cavitation aggressiveness index(CAI) based on the total derivative of pressure which identifies surface areas exposed to bubble collapses, the index is tested in two known cases documented in the open literature and seems to identify regions of potential cavitation damage.展开更多
基金Supported by National Nature Science Foundation of China(Grant No.51575245)Jiangsu Provincial Key research and development program(Grant No.BE2015134)Zhenjiang Municipal Key Research and Development Project(Grant No.KZ2020001).
文摘Cavitation generation methods have been used in multifarious directions because of their diversity,and numerous studies and discussions have been conducted on cavitation generation methods.This study aims to explore the generating mechanism and evolution law of volume alternate cavitation(VAC).In the VAC,liquid water is placed in an airtight container with a variable volume.As the volume alternately changes,the liquid water inside the container continues to cavitate.Then,the mixture turbulence model and in-cylinder dynamic grid model are adopted to conduct computational fluid dynamics simulation of volume alternate cavitation.In the simulation,the cloud images at seven heights on the central axis are monitored,and the phenomenon and mechanism of height and eccentricity are analyzed in detail.By employing the cavitation flow visualization method,the generating mechanism and evolution law of cavitation are revealed.The synergistic effects of experiments and high-speed camera capturing confirm the correctness of the simulation results.In the experiment,the volume change stroke of the airtight container is set to 20 mm,the volume change frequency is 18 Hz,and the shooting frequency of the high-speed camera is set to 10000 FPS.The experimental results indicate that the position of the cavitation phenomenon has a reasonable law during the whole evolution cycle of the cavitation cloud.Also,the volume alternation cycle corresponds to the generation,development,and collapse stages of cavitation bubbles.
基金financially supported by the National Natural Science Foundation of China(Nos.U21A20126 and 52006197)the National Science Foundation of Zhejiang Province(Nos.LQ21E060012 and LR20E090001)the Key Research and Development Program of Zhejiang Province(No.2021C05006)。
文摘The cavitating flow around a Delft Twist-11 hydrofoil is simulated using the large eddy simulation approach.The volume-of-fluid method incorporated with the Schnerr-Sauer cavitation model is utilized to track the water-vapor interface.Adaptive mesh refinement(AMR)is also applied to improve the simulation accuracy automatically.Two refinement levels are conducted to verify the dominance of AMR in predicting cavitating flows.Results show that cavitation features,including the U-type structure of shedding clouds,are consistent with experimental observations.Even a coarse mesh can precisely capture the phase field without increasing the total cell number significantly using mesh adaption.The predicted shedding frequency agrees fairly well with the experimental data under refinement level 2.This study illustrates that AMR is a promising approach to achieve accurate simulations for multiscale cavitating flows within limited computational costs.Finally,the force element method is currently adopted to investigate the lift and drag fluctuations during the evolution of cavitation structure.The mechanisms of lift and drag fluctuations due to cavitation and the interaction between vorticity forces and cavitation are explicitly revealed.
基金This project is supported by Teaching and Research Award Program for Outstanding Young Teachers in Higher Education Institutions of MOE, China (No.2002123) National Natural Science Foundation of China (No.50175097).
文摘The aim of this study is to suppress a cavitation near the orifice of arelief valve by changing the shape of a poppet. An experimental flow visualization technique and anumerical cavitating flow simulation, using a RNG kappa-epsilon turbulence model and a cavitationmodel, are employed to achieve the purpose. In the flow visualization, the cavitation phenomenonnear the orifice of a relief valve is observed using a transparent test valve body model and acamera. On the other hand, a three dimensional cavitating flow simulation is conducted to predictthe cavitation near the orifice of a relief valve. Six types of poppets are designed by changing theshape of a traditional poppet shape, which is expected to influence the cavitating flow in anorifice. In addition, the cavitation noise of a relief valve is measured and the noise spectrum isanalyzed. In conclusion, the cavitation intensity and the cavitation noise are reduced for anoptimal poppet geometry obtained in the present study.
文摘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.
文摘Computational Fluid Dynamics (CFD) simulations of cavitating flow through water hydraulic poppet valves were performed using advanced RNG k-epsilon turbulence model. The flow was turbulent, incompressible and unsteady, for Reynolds numbers greater than 43 000. The working fluid was water, and the structure of the valve was simplified as a two dimensional axisymmetric geometrical model. Flow field visualization was numerically achieved. The effects of inlet velocity, outlet pressure, opening size as well as poppet angle on cavitation intensity in the poppet valve were numerically investigated. Experimental flow visualization was conducted to capture cavitation images near the orifice in the poppet valve with 30° poppet angle using high speed video camera. The binary cavitating flow field distribution obtained from digital processing of the original cavitation image showed a good agreement with the numerical result.
基金supported by National Natural Science Foundation of China (Grant No.11172040)
文摘The Reynolds-averaged Navier-Stokes(RANS),such as the original k-ω two-equation closures,have been very popular in providing good prediction for a wide variety of flows with presently available computational resource.But for cavitating flows,the above equations noticeably over-predict turbulent production and hence effective viscosity.In this paper,the detached eddy simulation(DES) method for time-dependent turbulent cavitating flows is investigated.To assess the state-of-the-art of computational capabilities,different turbulence models including the widely used RANS model and DES model are conducted.Firstly,in order to investigate the grid dependency in computations,different grid sizes are adopted in the computation.Furthermore,the credibility of DES model is supported by the unsteady cavitating flows over a 2D hydrofoil.The results show that the DES model can effectively reduce the eddy viscosities.From the experimental validations regarding the force analysis,frequency and the unsteady cavity visualizations,more favorable agreement with experimental visualizations and measurements are obtained by DES model.DES model is better able to capture unsteady phenomena including cavity length and the resulting hydrodynamic characteristics,reproduces the time-averaged velocity quantitatively around the hydrofoil,and yields more acceptable and unsteady dynamics features.The DES model has shown to be effective in improving the overall predictive capability of unsteady cavitating flows.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52006197 and 51676174)the Natural Science Foundation of Jiangsu Province(Grant No.BK20180505)+1 种基金the National Science Foundation of Zhejiang Province(Grant Nos.LQ21E060012 and LR20E090001)the Key Research and Development Program of Zhejiang Province(Grant No.2020C01027)。
文摘The cavitating flow in different regimes has the intricate flow structure with multiple time and space scales.The present work develops a multiscale model by coupling the volume of fluid(VOF)method and a discrete bubble model(DBM),to simulate the cavitating flow in a convergent-divergent test section.The Schnerr-Sauer cavitation model is used to calculate the mass transfer rate to obtain the macroscale phase structure,and the simplified Rayleigh-Plesset equation is applied to simulate the growing and collapsing of discrete bubbles.An algorithm for bridging between the macroscale cavities and microscale bubbles is also developed to achieve the multiscale simulation.For the flow field,the very large eddy simulation(VLES)approach is applied.Conditions from inception to sheet/cloud cavitation regimes are taken into account and simulations are conducted.Compared with the experimental observations,it is shown that the cavitation inception,bubble clouds formation and glass cavity generation are all well represented,indicating that the proposed VOF-DBM model is a promising approach to accurately and comprehensively reveal the multiscale phase field induced by cavitation.
基金supported by the National Natural Science Foundation of China(Grant No.51479177).
文摘Based on a self-developed hydrodynamic cavitation device with different geometric parameters for circular multi-orifice plates,turbulence characteristics of cavitating flow behind multi-orifice plates,including the effects of orifice number and orifice layout on longitudinal velocity,turbulence intensity,and Reynolds stress,were measured with the particle image velocimetry(PIV)technique.Flow regimes of the cavitating flow were also observed with high-speed photography.The experimental results showed the following:(1)high-velocity multiple cavitating jets occurred behind the multi-orifice plates,and the cavitating flow fields were characterized by topological structures;(2)the longitudinal velocity at each cross-section exhibited a sawtooth-like distribution close to the multi-orifice plate,and each sawtooth indicated one jet issuing from one orifice;(3)there were similar magnitudes and forms for the longitudinal and vertical turbulence intensities at the same cross-section;(4)the variation in amplitude of Reynolds stress increased with an increase in orifice number;and(5)the cavitation clouds in the flow fields became denser with the increase in orifice number,and the clouds generated by the staggered layout of orifices were greater in number than those generated by the checkerboard-type one for the same orifice number.The experimental results can be used to analyze the mechanism of killing pathogenic microorganisms through hydrodynamic cavitation.
基金Supported by National Natural Science Foundation of China(Grant Nos.51209004,11172040)Beijing Municipal Natural Science Foundation of China(Grant No.3144043)
文摘Nowadays,most researchers focus on the cavity shedding mechanisms of unsteady cavitating flows over different objects,such as 2D/3D hydrofoils,venturi-type section,axisymmetric bodies with different headforms,and so on.But few of them pay attention to the differences of cavity shedding modality under different cavitation numbers in unsteady cavitating flows over the same object.In the present study,two kinds of shedding patterns are investigated experimentally.A high speed camera system is used to observe the cavitating flows over an axisymmetric blunt body and the velocity fields are measured by a particle image velocimetry(PIV)technique in a water tunnel for different cavitation conditions.The U-type cavitating vortex shedding is observed in unsteady cavitating flows.When the cavitation number is 0.7,there is a large scale cavity rolling up and shedding,which cause the instability and dramatic fluctuation of the flows,while at cavitation number of 0.6,the detached cavities can be conjunct with the attached part to induce the break-off behavior again at the tail of the attached cavity,as a result,the final shedding is in the form of small scale cavity and keeps a relatively steady flow field.It is also found that the interaction between the re-entrant flow and the attached cavity plays an important role in the unsteady cavity shedding modality.When the attached cavity scale is insufficient to overcome the re-entrant flow,it deserves the large cavity rolling up and shedding just as that at cavitation number of 0.7.Otherwise,the re-entrant flow is defeated by large enough cavity to induce the cavity-combined process and small scale cavity vortexes shedding just as that of the cavitation number of0.6.This research shows the details of two different cavity shedding modalities which is worthful and meaningful for the further study of unsteady cavitation.
基金Supported by National Natural Science Foundation of China (Grant Nos. 51635004, 11472078)。
文摘The water-lubricated bearings are usually the state of turbulent cavitating flow under high-speed conditions. And the distribution of cavitation bubbles and the interface effect between the two phases have not been included in previous studies on high-speed water-lubricated bearings. In order to study the influence of interface effect and cavitation bubble distribution on the dynamic characteristics of high-speed water-lubricated spiral groove thrust bearings(SGTB).A turbulent cavitating flow lubrication model based on two-phase fluid and population balance equation of bubbles was established in this paper. Stiffness and the damping coefficients of the SGTB were calculated using the perturbation pressure equations. An experimental apparatus was developed to verify the theoretical model. Simulating and experimental results show that the small-sized bubbles tend to generate in the turbulent cavitating flow when at a high rotary speed, and the bubbles mainly locate at the edges of the spiral groove. The simulating results also show that the direct stiffness coefficients are increased due to cavitation effect, and cross stiffness coefficients and damping coefficients are hardly affected by the cavitation effect. Turbulent effect on the dynamic characteristics of SGTB is much stronger than the cavitating effect.
基金the National Nature Science Foundation of China (11332011 and 11202215)
文摘Quasi-periodical evolutions such as shedding and collapsing of unsteady cloud cavitating flow, induce strong pressure fluctuations, what may deteriorate maneuvering stability and corrode surfaces of underwater vehicles. This paper analyzed effects on cavitation stability of a trip bar arranged on high-speed underwater projectile. Small scale water tank experiment and large eddy simulation using the open source software Open FOAM were used, and the results agree well with each other. Results also indicate that trip bar can obstruct downstream re-entrant jet and pressure wave propagation caused by collapse, resulting in a relatively stable sheet cavity between trip bar and shoulder of projectiles.
基金the National Key R&D Program of China(Grants 2016YFC0300800 and 2016YFC0300802)the National Natural Science Foundation of China(Grants 11772340 and 11672315)the Science and Technology on Water Jet Propulsion Laboratory(Grant 6142223190101).
文摘The transient cavitating flow around the Clark-Y hydrofoil is numerically investigated by the dynamic mode decomposition with criterion.Based on the ranking dominant modes,frequencies of the first four modes are in good accordance with those obtained by fast Fourier transform.Furthermore,the cavitating flow field is reconstructed by the first four modes,and the dominant flow features are well captured with the reconstructed error below 12%when compared to the simulated flow field.This paper offers a reference for observing and reconstructing the flow fields,and gives a novel insight into the transient cavitating flow features.
基金supported by the National Natural Science Foundation of China(Grant No.52009001)supported by the Postdoctoral Research Foundation of China(Granr No.2020M680380)+1 种基金the Natural Science Foundation Projectof Chongqing,Chongqing Science and Technology Commission(Grant No.cstc2021jcyj-msxmX1046)the Beijing Institute of Technology Research Fund Program for Young Scholars(Grant No.XSQD-202003008).
文摘Cavitation within the tip vortex(TV)flow remains a challenging issue in the design of high-speed and low-noise hydraulic machinery.In this paper,the TV cavitating flow around an elliptical hydrofoil is calculated by using large eddy simulation(LES)combined with a modified Schnerr-Sauer(S-S)cavitation model.The original S-S cavitation model is modified by taking into account the typical effect of vortex flow.The partial pressure term which can describe the vortex quantitatively and qualitatively is confirmed asρ_(m)ω^(2) x r _(c)^(2) ,and is considered into the R-P equation of the modified S-S cavitation model.Comparison between the numerical and experimental results shows good agreement in the form and evolution of cavities,including attached cavities(AC)and tip vortex cavities(TVC).The vorticity transport equation is utilized to investigate the dynamic mechanisms of the vortex development around the TVC.Further analyses indicate that cavitation in the TV flow influences the pressure in the core of the cavity and the local flow patterns.Typical vortex structures in the TV cavitating flow include TV,secondary vortex(SV)and wake vortex(WV).The direction and magnitude of the rotation effect can be described by axial vorticity which is drawn on the iso-surface of Q=1×105 s−2.The development of the TV cavitating flow can be divided into two stages:Stage I,the development and fusion of TV,SV,stage II,the dissipation of SV.The stretching term dominates the evolution of TV,and the dilatation term is the main reason in the mergence process of SV.
文摘In this paper,after the successful applications to open water propeller performance estimations,the influence of transition sensitive and modified mass transfer models tuned to account for the laminar flow in the prediction of the cavitation inception of marine propulsors is investigated from the point of view of the unsteady functioning and induced pressure pulses.The VP1304(also known as PPTC)test case,for which dedicated data were collected during several workshops,is considered first.After preliminary analyses using RANS,also Detached Eddy Simulations(DES)are included to better account for the vortex dynamics and its influence on pressure pulses.Similarly to what observed in uniform inflow,results show a better agreement with the available measurements of propeller performances and confirm the reliability of the proposed approaches for unsteady,non-cavitating,model scale propeller predictions.The overall cavitation pattern is improved too by the application of the transition sensitive correction to the mass transfer model,but the complex dynamics of bubble cavitation observed in experiments prevents quantitatively better predictions in terms of thrust/torque breakdown and induced pressure pulses levels regardless the use of RANS or DES methods.
基金Project supported by the National Natural Science Foundation of China (Grant No.52336001).
文摘In the present paper,the unsteady cavitating turbulent flow over the twisted NACA66 hydrofoil is investigated based on an modified shear stress transfer k-ωpartially averaged Navier-Stokes(MSST PANS)model,i.e.,new MSST PANS(NMSST PANS)model,where the production term of kinetic energy in the turbulence model is modified with helicity.Compared with the experimental data,cavitation evolution and its characteristic frequency are satisfactorily predicted by the proposed NMSST PANS model.It is revealed that the interaction among the main flow,the reentrant jets,and sheet cavitation causes the formation of the primary shedding cavity near the mid-span and the secondary shedding cavity at each side of the twisted hydrofoil,and further induces the remarkable pressure gradient around shedding cavities.Along with the development of the primary and the secondary shedding cavities,the great pressure gradient associated with large cavity volume variation promotes the vortical flow generation and the spatial deformation of vortex structure during cavitation evolution,and results in the primary and the secondary U-type vortices.Further,dynamic mode decomposition(DMD)analysis is utilized to confirm the interaction among the main flow,the main reentrant jet and two side reentrant jets,and cavitation.These results indicate that the proposed NMSST PANS model is suitable to simulate the complicated cavitating turbulent flow for various engineering applications.
基金support of the National Natural Science Foundation of China (Grant No. 50806078)the National High-Tech Research and Development Program of China (863 Program, Grant No. 2006AA05Z203).
文摘In this paper, by introducing the flow velocity item into the classical Rayleigh-Plesset dynamic equation, a newequation, which does not involve the time term and can describe the motion of cavitation bubble in the steadycavitating flow, has been obtained. By solving the new motion equation using Runge-Kutta fourth order methodwith adaptive step size control, the dynamic behaviors of cavitation bubble driven by the varying pressure fielddownstream of a venturi cavitation reactor are numerically simulated. The effects of liquid temperature (correspondingto the saturated vapor pressure of liquid), cavitation number and inlet pressure of venturi on radial motionof bubble and pressure pulse due to the radial motion are analyzed and discussed in detail. Some dynamicbehaviors of bubble different from those in previous papers are displayed. In addition, the internal relationshipbetween bubble dynamics and process intensification is also discussed. The simulation results reported in thiswork reveal the variation laws of cavitation intensity with the flow conditions of liquid, and will lay a foundationfor the practical application of hydrodynamic cavitation technology.
基金supported by the Natural Science Foundation of Heilongjiang Province of China(Grant No.A201409)the Special Fund Project for Technology Innovation Talent of Harbin(Grant No.2013RFLXJ007)the Fundamental Research Funds for the Central Universities(Grant No.HIT.NSRIF.201159)
文摘Developing a robust computational strategy to address the rich physical characteristic involved in the thermcdynamic effects on the cryogenic cavitation remains a challenge in research. The objective of the present study is to focus on developing mod- elling strategy to simulate cavitating flows in liquid nitrogen. For this purpose, numerical simulation over a 2D quarter caliber hydrofoil is investigated by calibrating cavitation model parameters and implementing the thermodynamic effects to the Zwart cavitation model. Experimental measurements of pressure and temperature are utilized to validate the extensional Zwart cavi- tation model. The results show that the cavitation dynamics characteristic under the cryogenic environment ale different from that under the isothermal conditions: the cryogenic case yields a substantially shorter cavity around the hydrofoil, and the pre- dicted pressure and temperature inside the cavity are steeper under the cryogenic conditions. Compared with the experimental data, the computational predictions with the modified evaporation and condensation parameters display better results than the default parameters from the room temperature liquids. Based on a wide range of computations and comparisons, the extension- al Zwart cavitation model may predict more accurately the quasi-steady cavitation over a hydrofoil in liquid nitrogen by pri- marily altering the evaporation rate near the leading edge and the condensation rate in the cavity closure region.
文摘Successful modeling and/or design of engineering systems often requires one to address the impact of multiple "design variables" on the prescribed outcome.There are often multiple,competing objectives based on which we assess the outcome of optimization.Since accurate,high fidelity models are typically time consuming and computationally expensive,comprehensive evaluations can be conducted only if an efficient framework is available.Furthermore,informed decisions of the model/hardware's overall performance rely on an adequate understanding of the global,not local,sensitivity of the individual design variables on the objectives.The surrogate-based approach,which involves approximating the objectives as continuous functions of design variables from limited data,offers a rational framework to reduce the number of important input variables,i.e.,the dimension of a design or modeling space.In this paper,we review the fundamental issues that arise in surrogate-based analysis and optimization,highlighting concepts,methods,techniques,as well as modeling implications for mechanics problems.To aid the discussions of the issues involved,we summarize recent efforts in investigating cryogenic cavitating flows,active flow control based on dielectric barrier discharge concepts,and lithium(Li)-ion batteries.It is also stressed that many multi-scale mechanics problems can naturally benefit from the surrogate approach for "scale bridging."
基金This work was financially supported by the National Natural Science Foundation of China(Grants 51822903 and 11772239)the Natural Science Foundation of Hubei Province(Grant 2018CFA010).The numerical calculations were done on the supercomputing system in the Supercomputing Center of Wuhan University.
文摘Large eddy simulation(LES)was coupled with a homogeneous cavitation model to study turbulent cavitating flows around a sphere.The simulations are in good agreement with available experimental data and the simulated accuracy has been evaluated using the LES verification and validation method.Various cavitation numbers are simulated to study important flow characteristics in the sphere wake,e.g.periodic cavity growth/contraction,interactions between the cloud and sheet cavitations and the vortex structure evolution.The spectral characteristics of the wake for typical cloud cavitation conditions were classified as the periodic cavitation mode,high Strouhal number mode and low Strouhal number mode.Main frequency distributions in the wake were analyzed and different dominant flow structures were identified for each of the three modes.Further,the cavitation and vortex relationship was also studied,which is an important issue associated with complex cavitating sphere wakes.Three types of cavitating vortex structures alternate,which indicates that three different cavity shedding regimes may exist in the wake.Analysis of vorticity transport equation shows a significant vorticity increase at the cavitation closure region and in the vortex cavitation region.This study provides a physical perspective to further understand the flow mechanisms in cavitating sphere wakes.
基金funding from the People Programme(Marie Curie Actions)of the European Union's Seventh Framework Programme FP7/2007-2013/under REA grant agreement No.324313
文摘The paper proposes a methodology within the Reynolds averaged Navier Stokes(RANS) solvers for cavitating flows capable of predicting the flow regions of bubble collapse and the potential aggressiveness to material damage. An aggressiveness index is introduced, called cavitation aggressiveness index(CAI) based on the total derivative of pressure which identifies surface areas exposed to bubble collapses, the index is tested in two known cases documented in the open literature and seems to identify regions of potential cavitation damage.