By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improv...By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.展开更多
To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the ch...To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us.We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of(13944.7021,13946.5333)by the method of bisection.Through Fourier analysis,it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval.Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram,Kolmogorov entropy and maximal Lyapunov exponent.The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.展开更多
A streamline upwind/Petrov-Galerkin (SUPG) finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for-...A streamline upwind/Petrov-Galerkin (SUPG) finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.展开更多
To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of...To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of theoretical analyses and numerical simulations. In this paper, we study the characteristics of the second Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme recently developed by us. We numerically identify the critical Reynolds number of the second Hopf bifurcation located in the interval of(11093.75, 11094.3604) by bisection. In addition, we find that there are two dominant frequencies in its spectral diagram when the flow is in the status of the second Hopf bifurcation, while only one dominant frequency is identified if the flow is in the first Hopf bifurcation via the Fourier analysis. More interestingly, the flow phase portrait of velocity components is found to make transition from a regular elliptical closed form for the first Hopf bifurcation to a non-elliptical closed form with self-intersection for the second Hopf bifurcation. Such characteristics disclose flow in a quasi-periodic state when the second Hopf bifurcation occurs.展开更多
Particle based methods can be used for both the simulations of solid and fluid phases in multiphase medium, such as the discrete-element method for solid phase and the smoothed particle hydrodynamics for fluid phase. ...Particle based methods can be used for both the simulations of solid and fluid phases in multiphase medium, such as the discrete-element method for solid phase and the smoothed particle hydrodynamics for fluid phase. This paper presents a computational method combining these two methods for solid-liquid medium. The two phases are coupled by using an improved model from a reported Lagrangian-Eulerian method. The technique is verified by simulating liquid-solid flows in a two-dimensional lid-driven cavity.展开更多
The present study refers to a cavitating Venturi type section geometry characterized by a convergent angle of 18°?and a divergent angle of about 8°?where the sheet cavity presents typical self-oscillation be...The present study refers to a cavitating Venturi type section geometry characterized by a convergent angle of 18°?and a divergent angle of about 8°?where the sheet cavity presents typical self-oscillation behavior with quasi-periodic vapor clouds shedding. This work is an extension of previous works concerning void ratio measurements and velocity fields using double optical probe and constitutes a complete analysis of the two-phase structure of unsteady cavitating flow. This paper provides a new method based on conditional and phase averaging technique with wall pressure signal to treat experimental data in order to evaluate more precisely time-averaged and rms values of the void ratio and instantaneous velocity fields. Conditional analysis shows a different behavior of the two-phase flow dynamics leading to highlight high void ratio events linked to the break-off cycle. Unsteady phase averaging of the optical probe signal gives the evolution of the void ratio at each studied location in the venturi and shows that the fluctuations close to the wall (where the re-entrant jet is predominant) are in phase with the upper part of the cavity instead of the thickness of the cavity which is unchanged.展开更多
An investigation on flow and heat transfer due to mixed convection, in a lid-driven rectangular cavity filled with Cu- water nanofluids and submitted to uniform heat flux along with its vertical short sides, has been ...An investigation on flow and heat transfer due to mixed convection, in a lid-driven rectangular cavity filled with Cu- water nanofluids and submitted to uniform heat flux along with its vertical short sides, has been conducted numerically by solving the full governing equations with the finite volume method and the SIMPLER algorithm. In the case of a slender enclosure, these equations are considerably reduced by using the parallel flow concept. Solutions, for the flow and temperature fields, and the heat transfer rate, have been obtained depending on the governing parameters, which are the Reynolds, the Richardson numbers and the solid volume fraction of nanoparticles. A perfect agreement has been found between the results of the two approaches for a wide range of the abovementioned parameters. It has been shown that at low and high Richardson numbers, the convection is ensured by lid and buoyancy-driven effects, respectively, whereas between these extremes, both mechanisms compete. Moreover, the addition of Cu-nanoparticles, into the pure water, has been seen enhancing and degrading heat transfer by lid and buoyancy-driven effects, respectively.展开更多
As an integral part of the internal air system of aero-engines,the axial throughflow of the cooling air can interact with the cavity flow between the rotating compressor disks,forming a threedimensional,unsteady,and u...As an integral part of the internal air system of aero-engines,the axial throughflow of the cooling air can interact with the cavity flow between the rotating compressor disks,forming a threedimensional,unsteady,and unstable flow field.The flow characteristics in an engine-like rotating multi-stage cavity with throughflow were investigated using particle image velocimetry,flow visualization technology and three-dimensional unsteady Reynolds-Averaged Navier-Stokes (RANS)simulations.The focus of current research was to understand the distribution of the mean swirl ratio and its variation with a wide range of non-dimensional parameters in the co-rotating cavity with high inlet pre-swirl axial throughflow.The maximum axial Reynolds number and rotational Reynolds numbers could reach 4.41×10^(4)and 1.24×10^(6),respectively.The velocity measurement results indicate that the mean swirl ratio is greater than 1 and decreases with an increase in the radial position.The flow structure is dominated by the Rossby number,and two different flow patterns (flow penetration and flow stratification) are identified and confirmed by flow visualization images.In the absence of buoyancy,the flow penetration caused by the precession of the throughflow makes it easier for the throughflow to reach a high radius region.Satisfactory consistency of results between measurements and numerical calculations is obtained.This study provides a theoretical basis and data support for toroidal vortex breakdown,which is of practical significance for the design of high-pressure compressor cavities.展开更多
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.展开更多
A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physica...A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physical domain, and the particle streaming process is carried out in the corresponding computational domain, which is transferred from the physical domain using interpolation method. For the interpolation calculation, a second-order upwind interpolation method is adopted on internal lattice nodes in flow fields while a second-order central interpolation algorithm is employed at neighbor-boundary lattice nodes. Then the above-mentioned model and algorithms are used to numerically simulate the 3D flows in the lid-driven cavity at Reynolds numbers of 100, 400 and 1000 on non-uniform meshes. Various vortices on the x-y, y-z and x-z symmetrical planes are successfully predicted, and their changes in position with the Reynolds number increasing are obtained. The velocity profiles of u component along the vertical centerline and w component along the horizontal centerline are both in good agreement with the data in literature and the calculated results on uniform meshes. Besides, the velocity vector distributions on various cross sections in lid-driven cavity predicted on non-uniform meshes are compared with those simulated on uniform meshes and those in the literature. All the comparisons and validations show that the 3D lattice Boltzmann model and all the numerical algorithms on non-uniform meshes are accurate and reliable to predict effectively flow fields.展开更多
This paper investigates the chaotic lid-driven square cavity flows at extreme Reynolds numbers.Several observations have been made from this study.Firstly,at extreme Reynolds numbers two principles add at the genesis ...This paper investigates the chaotic lid-driven square cavity flows at extreme Reynolds numbers.Several observations have been made from this study.Firstly,at extreme Reynolds numbers two principles add at the genesis of tiny,loose counterclockwise-or clockwise-rotating eddies.One concerns the arousing of them owing to the influence of the clockwise-or counterclockwise currents nearby;the other,the arousing of counterclockwise-rotating eddies near attached to the moving(lid)top wall which moves from left to right.Secondly,unexpectedly,the kinetic energy soon reaches the qualitative temporal limit’s pace,fluctuating briskly,randomly inside the total kinetic energy range,fluctuations which concentrate on two distinct fragments:one on its upper side,the upper fragment,the other on its lower side,the lower fragment,switching briskly,randomly from each other;and further on many small fragments arousing randomly within both,switching briskly,randomly from one another.As the Reynolds number Re→∞,both distance and then close,and the kinetic energy fluctuates shorter and shorter at the upper fragment and longer and longer at the lower fragment,displaying tall high spikes which enlarge and then disappear.As the time t→∞(at the Reynolds number Re fixed)they recur from time to time with roughly the same amplitude.For the most part,at the upper fragment the leading eddy rotates clockwise,and at the lower fragment,in stark contrast,it rotates counterclockwise.At Re=109 the leading eddy-at its qualitative temporal limit’s pace—appears to rotate solely counterclockwise.展开更多
In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split i...In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.展开更多
A two-phase mixture model was established to study unsteady cavitating flows. A local compressible system of equations was derived by introducing a density-pressure function to account for the two-phase flow of water/...A two-phase mixture model was established to study unsteady cavitating flows. A local compressible system of equations was derived by introducing a density-pressure function to account for the two-phase flow of water/vapor and the transition from one phase to the other. An algorithm for solving the variable-density Navier-Stokes equations of cavitating flow problem was put forward. The numerical results for unsteady characteristics of cavitating flows on a 2D NACA hydrofoil coincide well with experimental data.展开更多
In this study,an experiment was performed to clarify the flow field,in which the jets were normally injected into a main supersonic flow surrounded by a porous cavity,and this report figures out interaction between st...In this study,an experiment was performed to clarify the flow field,in which the jets were normally injected into a main supersonic flow surrounded by a porous cavity,and this report figures out interaction between starting shock wave and porous cavity.In the experiment,a porous cavity is attached to a main duct and jets and rods are inserted to the main duct on the porous cavity.To reveal this flow field,the thermal tuft probe was adopted to experimentally investigate the flow in the cavity.In the experiments,the effect of the porous cavity with jets or rods on the flow field is studied by means of visualization of schlieren method with a high speed camera and measurement of cavity flow with thermal tuft probe.As a results,frequency analysis of output of the thermal tuft probe revealed that some clear dominant frequencies were confirmed when the starting shock wave existed around the porous cavity in all cases of jets and rods arrangements.Moreover,visualization of schlieren method with a high speed camera clarified that a starting shock wave had the same dominant frequencies as that of the flow fluctuation in the cavity only around the cavity.展开更多
Supersonic cavity flows are characterized by compression and expansion waves, shear layer, and oscillations inside the cavity. For decades, investigations into cavity flows have been conducted, mostly with flows at ze...Supersonic cavity flows are characterized by compression and expansion waves, shear layer, and oscillations inside the cavity. For decades, investigations into cavity flows have been conducted, mostly with flows at zero pressure gradient entering the cavity in straight walls. Since cavity flows on curved walls exert centrifugal force, the features of these flows are likely to differ from those of straight wall flows. The aim of the present work is to study the flow physics of a cavity that is cut out on a curved wall. Steady and unsteady numerical simulations were carried out for supersonic flow through curved channels over the cavity with L/H = 1. A straight channel flow was also analyzed which serves as the base model. The velocity gradient along the width of the channel was observed to increase with increasing the channel curvature for both concave and convex channels. The pressure on the cavity floor increases with the increase in channel curvature for concave channels and decreases for convex channels. Moreover, unsteady flow characteristics are more dependent on channel curvature under supersonic free stream conditions.展开更多
基金Project supported by the National Natural Science Foundation of China (Grant No 70271069).
文摘By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.
基金Project supported by the National Natural Science Foundation of China(Grant No.12162001)the Natural Science Foundation of Ningxia(Grant No.2019AAC03129)the Construction Project of First-Class Disciplines in Ningxia Higher Education(Grant No.NXYLXK2017B09)。
文摘To date,there are very few studies on the transition beyond second Hopf bifurcation in a lid-driven square cavity,due to the difficulties in theoretical analysis and numerical simulations.In this paper,we study the characteristics of the third Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme rectently developed by us.We numerically identify the critical Reynolds number of the third Hopf bifurcation located in the interval of(13944.7021,13946.5333)by the method of bisection.Through Fourier analysis,it is discovered that the flow becomes chaotic with a characteristic of period-doubling bifurcation when the Reynolds number is beyond the third bifurcation critical interval.Nonlinear time series analysis further ascertains the flow chaotic behaviors via the phase diagram,Kolmogorov entropy and maximal Lyapunov exponent.The phase diagram changes interestingly from a closed curve with self-intersection to an unclosed curve and the attractor eventually becomes strange when the flow becomes chaotic.
基金the National Natural Science Foundation of China (Grants 41372301 and 51349011)the Preeminent Youth Talent Project of Southwest University of Science and Technology (Grant 13zx9109)
文摘A streamline upwind/Petrov-Galerkin (SUPG) finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11601013 and 91530325)。
文摘To date, there are very few studies on the second Hopf bifurcation in a driven square cavity, although there are intensive investigations focused on the first Hopf bifurcation in literature, due to the difficulties of theoretical analyses and numerical simulations. In this paper, we study the characteristics of the second Hopf bifurcation in a driven square cavity by applying a consistent fourth-order compact finite difference scheme recently developed by us. We numerically identify the critical Reynolds number of the second Hopf bifurcation located in the interval of(11093.75, 11094.3604) by bisection. In addition, we find that there are two dominant frequencies in its spectral diagram when the flow is in the status of the second Hopf bifurcation, while only one dominant frequency is identified if the flow is in the first Hopf bifurcation via the Fourier analysis. More interestingly, the flow phase portrait of velocity components is found to make transition from a regular elliptical closed form for the first Hopf bifurcation to a non-elliptical closed form with self-intersection for the second Hopf bifurcation. Such characteristics disclose flow in a quasi-periodic state when the second Hopf bifurcation occurs.
基金supported by Department of Energy and Process Engineering,Norwegian University of Science and TechnologyInstitute for Energy Technology and SINTEF through the FACE(Multiphase Flow Assurance Innovation Center) Project
文摘Particle based methods can be used for both the simulations of solid and fluid phases in multiphase medium, such as the discrete-element method for solid phase and the smoothed particle hydrodynamics for fluid phase. This paper presents a computational method combining these two methods for solid-liquid medium. The two phases are coupled by using an improved model from a reported Lagrangian-Eulerian method. The technique is verified by simulating liquid-solid flows in a two-dimensional lid-driven cavity.
文摘The present study refers to a cavitating Venturi type section geometry characterized by a convergent angle of 18°?and a divergent angle of about 8°?where the sheet cavity presents typical self-oscillation behavior with quasi-periodic vapor clouds shedding. This work is an extension of previous works concerning void ratio measurements and velocity fields using double optical probe and constitutes a complete analysis of the two-phase structure of unsteady cavitating flow. This paper provides a new method based on conditional and phase averaging technique with wall pressure signal to treat experimental data in order to evaluate more precisely time-averaged and rms values of the void ratio and instantaneous velocity fields. Conditional analysis shows a different behavior of the two-phase flow dynamics leading to highlight high void ratio events linked to the break-off cycle. Unsteady phase averaging of the optical probe signal gives the evolution of the void ratio at each studied location in the venturi and shows that the fluctuations close to the wall (where the re-entrant jet is predominant) are in phase with the upper part of the cavity instead of the thickness of the cavity which is unchanged.
文摘An investigation on flow and heat transfer due to mixed convection, in a lid-driven rectangular cavity filled with Cu- water nanofluids and submitted to uniform heat flux along with its vertical short sides, has been conducted numerically by solving the full governing equations with the finite volume method and the SIMPLER algorithm. In the case of a slender enclosure, these equations are considerably reduced by using the parallel flow concept. Solutions, for the flow and temperature fields, and the heat transfer rate, have been obtained depending on the governing parameters, which are the Reynolds, the Richardson numbers and the solid volume fraction of nanoparticles. A perfect agreement has been found between the results of the two approaches for a wide range of the abovementioned parameters. It has been shown that at low and high Richardson numbers, the convection is ensured by lid and buoyancy-driven effects, respectively, whereas between these extremes, both mechanisms compete. Moreover, the addition of Cu-nanoparticles, into the pure water, has been seen enhancing and degrading heat transfer by lid and buoyancy-driven effects, respectively.
文摘As an integral part of the internal air system of aero-engines,the axial throughflow of the cooling air can interact with the cavity flow between the rotating compressor disks,forming a threedimensional,unsteady,and unstable flow field.The flow characteristics in an engine-like rotating multi-stage cavity with throughflow were investigated using particle image velocimetry,flow visualization technology and three-dimensional unsteady Reynolds-Averaged Navier-Stokes (RANS)simulations.The focus of current research was to understand the distribution of the mean swirl ratio and its variation with a wide range of non-dimensional parameters in the co-rotating cavity with high inlet pre-swirl axial throughflow.The maximum axial Reynolds number and rotational Reynolds numbers could reach 4.41×10^(4)and 1.24×10^(6),respectively.The velocity measurement results indicate that the mean swirl ratio is greater than 1 and decreases with an increase in the radial position.The flow structure is dominated by the Rossby number,and two different flow patterns (flow penetration and flow stratification) are identified and confirmed by flow visualization images.In the absence of buoyancy,the flow penetration caused by the precession of the throughflow makes it easier for the throughflow to reach a high radius region.Satisfactory consistency of results between measurements and numerical calculations is obtained.This study provides a theoretical basis and data support for toroidal vortex breakdown,which is of practical significance for the design of high-pressure compressor cavities.
基金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 the National Natural Science Foundation of China (Grant Nos. 51179192, 50779069, 51139007)the Program for New Century Excellent Talents in University (NCET) (Grant No. NETC-10-0784)+1 种基金the National Hi-Tech Research and Development Program of China ("863" Project) (Grant No. 2011AA100505)the Chinese Universities Scientific Fund (Grant No. 2013RC045)
文摘A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physical domain, and the particle streaming process is carried out in the corresponding computational domain, which is transferred from the physical domain using interpolation method. For the interpolation calculation, a second-order upwind interpolation method is adopted on internal lattice nodes in flow fields while a second-order central interpolation algorithm is employed at neighbor-boundary lattice nodes. Then the above-mentioned model and algorithms are used to numerically simulate the 3D flows in the lid-driven cavity at Reynolds numbers of 100, 400 and 1000 on non-uniform meshes. Various vortices on the x-y, y-z and x-z symmetrical planes are successfully predicted, and their changes in position with the Reynolds number increasing are obtained. The velocity profiles of u component along the vertical centerline and w component along the horizontal centerline are both in good agreement with the data in literature and the calculated results on uniform meshes. Besides, the velocity vector distributions on various cross sections in lid-driven cavity predicted on non-uniform meshes are compared with those simulated on uniform meshes and those in the literature. All the comparisons and validations show that the 3D lattice Boltzmann model and all the numerical algorithms on non-uniform meshes are accurate and reliable to predict effectively flow fields.
基金supported in part by the National Science Foundation Grants No.DMS-0906440 and No.DMS-1206438.
文摘This paper investigates the chaotic lid-driven square cavity flows at extreme Reynolds numbers.Several observations have been made from this study.Firstly,at extreme Reynolds numbers two principles add at the genesis of tiny,loose counterclockwise-or clockwise-rotating eddies.One concerns the arousing of them owing to the influence of the clockwise-or counterclockwise currents nearby;the other,the arousing of counterclockwise-rotating eddies near attached to the moving(lid)top wall which moves from left to right.Secondly,unexpectedly,the kinetic energy soon reaches the qualitative temporal limit’s pace,fluctuating briskly,randomly inside the total kinetic energy range,fluctuations which concentrate on two distinct fragments:one on its upper side,the upper fragment,the other on its lower side,the lower fragment,switching briskly,randomly from each other;and further on many small fragments arousing randomly within both,switching briskly,randomly from one another.As the Reynolds number Re→∞,both distance and then close,and the kinetic energy fluctuates shorter and shorter at the upper fragment and longer and longer at the lower fragment,displaying tall high spikes which enlarge and then disappear.As the time t→∞(at the Reynolds number Re fixed)they recur from time to time with roughly the same amplitude.For the most part,at the upper fragment the leading eddy rotates clockwise,and at the lower fragment,in stark contrast,it rotates counterclockwise.At Re=109 the leading eddy-at its qualitative temporal limit’s pace—appears to rotate solely counterclockwise.
基金financially supported by the National Natural Science Foundation of China(Grant No.51349011)the Foundation of Si’chuan Educational Committee(Grant No.17ZB0452)+1 种基金the Innovation Team Project of Si’chuan Educational Committee(Grant No.18TD0019)the Longshan Academic Talent Research Support Program of the Southwest of Science and Technology(Grant Nos.18LZX715 and 18LZX410)
文摘In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.
基金Project supported by the National Natural Science Foundation of China (Grant No: 10372061) and the National Defense Key Laboratory on Hydrodynamics
文摘A two-phase mixture model was established to study unsteady cavitating flows. A local compressible system of equations was derived by introducing a density-pressure function to account for the two-phase flow of water/vapor and the transition from one phase to the other. An algorithm for solving the variable-density Navier-Stokes equations of cavitating flow problem was put forward. The numerical results for unsteady characteristics of cavitating flows on a 2D NACA hydrofoil coincide well with experimental data.
文摘In this study,an experiment was performed to clarify the flow field,in which the jets were normally injected into a main supersonic flow surrounded by a porous cavity,and this report figures out interaction between starting shock wave and porous cavity.In the experiment,a porous cavity is attached to a main duct and jets and rods are inserted to the main duct on the porous cavity.To reveal this flow field,the thermal tuft probe was adopted to experimentally investigate the flow in the cavity.In the experiments,the effect of the porous cavity with jets or rods on the flow field is studied by means of visualization of schlieren method with a high speed camera and measurement of cavity flow with thermal tuft probe.As a results,frequency analysis of output of the thermal tuft probe revealed that some clear dominant frequencies were confirmed when the starting shock wave existed around the porous cavity in all cases of jets and rods arrangements.Moreover,visualization of schlieren method with a high speed camera clarified that a starting shock wave had the same dominant frequencies as that of the flow fluctuation in the cavity only around the cavity.
基金supported by Advanced Research Center Program(NRF-2013R1A5A1073861)through the National Research Foundation of Korea(NRF)
文摘Supersonic cavity flows are characterized by compression and expansion waves, shear layer, and oscillations inside the cavity. For decades, investigations into cavity flows have been conducted, mostly with flows at zero pressure gradient entering the cavity in straight walls. Since cavity flows on curved walls exert centrifugal force, the features of these flows are likely to differ from those of straight wall flows. The aim of the present work is to study the flow physics of a cavity that is cut out on a curved wall. Steady and unsteady numerical simulations were carried out for supersonic flow through curved channels over the cavity with L/H = 1. A straight channel flow was also analyzed which serves as the base model. The velocity gradient along the width of the channel was observed to increase with increasing the channel curvature for both concave and convex channels. The pressure on the cavity floor increases with the increase in channel curvature for concave channels and decreases for convex channels. Moreover, unsteady flow characteristics are more dependent on channel curvature under supersonic free stream conditions.