A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevo...A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.展开更多
The idea that the collapse proceeds from the outer boundary of the cavity cloud towards its center for the ultrasonic cavitation proposed by Hasson and Morch in 1980s is further developed for calculating the collapse ...The idea that the collapse proceeds from the outer boundary of the cavity cloud towards its center for the ultrasonic cavitation proposed by Hasson and Morch in 1980s is further developed for calculating the collapse pressure and boundaries of cavity cloud at the collapse stage of bubbles for hydraulic cavitation flow in Venturi in present research. The numerical simulation is carried out based on Gilmore＇s eouations of bubble dynamics, which take account of the compressibility of fluid besides the viscosity and interfacial tension. The collapse of the cavity cloud is considered to proceed layer by layer from the outer cloud towards its inner part. The simulation results indicate that thepredicted boundaries of the cavity cloudat the collapse stage agree.well with the exPerimental ones.It is also found that the maximum collapse pressure of the cavity cloud is several times as high as the collapse pressure of outside boundary, and it is located at a point in the axis, where the cavity cloud disappears completely. This means that a cavity cloud has higher collapse pressure or strength than that of a single bubble due to the interactions of the bubbles. The effects of operation and structural parameters on the collapse pressure are also analyzed in detail.展开更多
Vortex methods have been alternative tools of finite element and finite difference methods for several decades. This paper presents a brief review of vortex method development in the last decades and introduces effici...Vortex methods have been alternative tools of finite element and finite difference methods for several decades. This paper presents a brief review of vortex method development in the last decades and introduces efficient vortex methods developed for high Reynolds number bluff body flows and suitable for running on parallel computer architectures. Included in this study are particle strength exchange methods, core-spreading method, deterministic particle method and hybrid vortex methods. Combined with conservative methods, vortex methods can comprise the most available tools for simulations of three-dimensional complex bluff body flows at high Reynolds numbers.展开更多
In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NAC...In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.展开更多
We investigate the symmetry reduction for the two-dimensional incompressible Navier-Stokes equationin conventional stream function form through Lie symmetry method and construct some similarity reduction solutions.Two...We investigate the symmetry reduction for the two-dimensional incompressible Navier-Stokes equationin conventional stream function form through Lie symmetry method and construct some similarity reduction solutions.Two special cases in [D.K.Ludlow,P.A.Clarkson,and A.P.Bassom,Stud.Appl.Math.103 (1999) 183] and a theoremin [S.Y.Lou,M.Jia,X.Y.Tang,and F.Huang,Phys.Rev.E 75 (2007) 056318] are retrieved.展开更多
The temperature separation was discovered inside the short vortex chamber （H/D = 0.18）. Experiments revealed that the highest temperature of the periphery was 465 ℃, and the lowest temperature of the central zone w...The temperature separation was discovered inside the short vortex chamber （H/D = 0.18）. Experiments revealed that the highest temperature of the periphery was 465 ℃, and the lowest temperature of the central zone was -45 ℃ （the compressed air was pumped into the chamber at room temperature）. The objective of this paper is to proof that this temperature separation effect cannot be explained by conventional heat transfer processes. To explain this phenomenon, the concept of PGEW （Pressure Gradient Elastic Waves） is proposed. PGEW are kind of elastic waves, which operate in compressible fluids with pressure gradients and density fluctuations. The result of PGEW propagation is a heat transfer from area of low pressure to high pressure zone. The physical model of a gas in a strong field of mass forces is proposed to substantiate the PGEW existence. This physical model is intended for the construction of a theory of PGEW. Understanding the processes associated with the PGEW permits the possibility of creating new devices for energy saving and low potential heat utilization, which have unique properties.展开更多
This paper studies the phase effect in mode coupling of Kelvin-Helmholtz instability in two-dimensionalincompressible fluid.It is found that there is an important growth phenomenon of every mode in the mode couplingpr...This paper studies the phase effect in mode coupling of Kelvin-Helmholtz instability in two-dimensionalincompressible fluid.It is found that there is an important growth phenomenon of every mode in the mode couplingprocess.The growth changes periodically with phase difference and in the condition of our simulation the period is about0.7π.The period characteristic is apparent in all stage of the mode coupling process,especially in the relatively laterstage.展开更多
In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann （LB） model for shacked compressible flows. The proposed flux limit...In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann （LB） model for shacked compressible flows. The proposed flux limiter scheme is efficient in decreasing the artificial oscillations and numerical diffusion around the interface. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. Numerical simulations for the Richtmyer-Meshkov instability show that with the new model the computed interfaces are smoother and more consistent with physical analysis. The growth rates of bubble and spike present a satisfying agreement with the theoretical predictions and other numerical simulations.展开更多
The equations of motion of a bubble, expanding adiabatically through an incompressible viscous fluid, are deduced when the centre of the bubble moves in a vertical plane in the presence of gravitational acceleration, ...The equations of motion of a bubble, expanding adiabatically through an incompressible viscous fluid, are deduced when the centre of the bubble moves in a vertical plane in the presence of gravitational acceleration, acting vertically downwards. The non-linear equations of motion obtained are solved numerically for different values of the various parameters of the problem. The path traced by the centre of the bubble and velocity of the centre, the change of radius R with time, and the influence of the buoyancy force, which is experienced by the expanding bubble for different values of the gravitational acceleration on these quantities, are investigated. The radius R（t） of the bubble is found to vary periodically with time when the acceleration due to gravity is small. But when the acceleration due to gravity increases, this periodicity in the value of R（t） with t is lost. The influence of viscosity in determining the periodicity of the bubble motion is also investigated.展开更多
We investigate the nonlinear instability of a smooth steady density profile solution to the threedimensional nonhomogeneous incompressible Navier-Stokes equations in the presence of a uniform gravitational field,inclu...We investigate the nonlinear instability of a smooth steady density profile solution to the threedimensional nonhomogeneous incompressible Navier-Stokes equations in the presence of a uniform gravitational field,including a Rayleigh-Taylor steady-state solution with heavier density with increasing height(referred to the Rayleigh-Taylor instability).We first analyze the equations obtained from linearization around the steady density profile solution.Then we construct solutions to the linearized problem that grow in time in the Sobolev space H k,thus leading to a global instability result for the linearized problem.With the help of the constructed unstable solutions and an existence theorem of classical solutions to the original nonlinear equations,we can then demonstrate the instability of the nonlinear problem in some sense.Our analysis shows that the third component of the velocity already induces the instability,which is different from the previous known results.展开更多
The large eddy simulation (LES) of compressible turbulent channel flows at three different Mach numbers is performed in the present work,by extending the dynamic mixed subgrid-scale (SGS) model to compressible flows.T...The large eddy simulation (LES) of compressible turbulent channel flows at three different Mach numbers is performed in the present work,by extending the dynamic mixed subgrid-scale (SGS) model to compressible flows.The turbulent statistics agree well with those from the existing direct numerical simulation (DNS) results,indicating that the LES method established in the present work is reliable.The analysis of the turbulent fluctuations computed by the present LES reveals that the flows considered in this work follow the Morkovin's hypothesis.Thus,the compressibility effects are dominated by the mean field properties,and the relevant statistical ratios are invariant to the variation of Mach number.The near-wall streamwise streaks are more coherent and the spacing between streaks is wider as the Mach number increases.This can be regarded as a direct feature characterized by the compressibility effects.The restrained influences of compressibility effects on the production and dissipation of the turbulence kinetic energy are also identified based on the present LES results.展开更多
Based on the standard k-ε turbulence model, a new compressible k-ε model considering the pressure expansion influence due to the compressibility of fluid is developed and applied to the simulation of 3D transonic tu...Based on the standard k-ε turbulence model, a new compressible k-ε model considering the pressure expansion influence due to the compressibility of fluid is developed and applied to the simulation of 3D transonic turbulent flows in a nozzle and a cascade. The Reynolds avenged N-S equations in generalized curvilinear coordinates are solved with implementation of the new model. The high resolution TVD scheme is used to discretize the convective terms. The numerical results show that the compressible k-ε model behaves well in the simulation of transonic internal turbulent flows.展开更多
A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) ...A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) stabilization scheme is used for the formulation of the Navier-Stokes equations. For the spatial discretization, the convection term is treated explicitly, while the viscous term is treated implicitly, and for the temporal discretization, a three-step method is employed. The present method is applied to simulate the lid driven cavity problems with different geometries at low and high Reynolds numbers. The results compared with other numerical experiments are found to be feasible and satisfactory.展开更多
Based on the immersed boundary method,a fast simulation for solving unsteady,incompressible,viscous flow associated with the oscillating cascade is established on a quasi-three-dimensional coordinate system.The numeri...Based on the immersed boundary method,a fast simulation for solving unsteady,incompressible,viscous flow associated with the oscillating cascade is established on a quasi-three-dimensional coordinate system.The numerical method is applied to the simulation of the flow passing an oscillating circular cylinder which is forced to move in X direction under prescribed motions in water at rest at low Keulegan-Carpenter numbers.Then vor-tex-induced vibration of a cylinder with two degrees of freedom which oscillates in in-line direction and transverse direction is simulated using this method.The results are in good agreement with the previous research.Then the method is extended to the oscillating cascade simulation of making various comparisons.It is found that the IBPA(inter blade phase angle) will change as the time goes on,because of the non-uniformity of the flow in the circumferential direction,until the oscillating cascade goes to a stable situation.The reduced velocity and the number of blades are chosen to investigate the effects of them on IBPA.The results indicate that both the reduced velocity and the number of blades are the main factors which influence IBPA.It is worth noting that the coupling process is not necessary to generate any body-fitting grids,which makes it much faster in computational process for such a complicated fluid-structure interaction problem.展开更多
In this paper, the flow fields of underexpanded impinging jet issued from rectangular nozzles of aspect ratio 1, 3 and 5 are numerically and experimentally studied. Two dimensional temperature and pressure distributio...In this paper, the flow fields of underexpanded impinging jet issued from rectangular nozzles of aspect ratio 1, 3 and 5 are numerically and experimentally studied. Two dimensional temperature and pressure distributions are measured by using infrared camera and the combination of a pressure scanning device and a stepping motor, respectively. The variation of the stagnation pressure on the impinging plate reveals that a hysteretic phenomenon exists during the increasing and decreasing of the pressure ratio for the aspect ratio of 3.0 and 5.0. It is also found that the nozzle of aspect ratio 1.0 caused the largest total pressure loss Pc /p0= 0.27 at the pressure ratio of Po /p0 = 6.5, where Pc is the stagnation center pressure on the wall, P0 the upstream stagnation pressure, Pb the ambient pressure. The other two nozzles showed that the pressure loss Pc / P0 =0.52 and 0.55 were achieved by the nozzles of the aspect ratio 3,0 and 5.0, respectively. The comparison between the calculations and experiments is fairly good, showing the three dimensional streamlines and structures of the shock waves in the jets. However, the hysteresis of the pressure variations observed in the experiments between the pressure ratio of 3.5 and 4.5 cannot be confirmed in the calculations.展开更多
The crossflow instability of a three-dimensional (3-D) boundary layer is an important factor which affects the transition over a swept-wing.In this report,the primary instability of the incompressible flow over a swep...The crossflow instability of a three-dimensional (3-D) boundary layer is an important factor which affects the transition over a swept-wing.In this report,the primary instability of the incompressible flow over a swept wing is investigated by solving nonlinear parabolized stability equations (NPSE).The Floquet theory is applied to study the dependence of the secondary and high-frequency instabilities on curvature,Reynolds number and angle of swept (AOS).The computational results show that the curvature in the present case has no significant effect on the secondary instabilities.It is generally believed that the secondary instability growth rate increases with the magnitude of the nonlinear mode of crossflow vortex.But,at a certain state,when the Reynolds number is 3.2 million,we find that the secondary instability growth rate becomes smaller even when the magnitude of the nonlinear mode of the crossflow vortex is larger.The effect of the angle of swept at 35,45 and 55 degrees,respectively,is also studied in the framework of the secondary linear stability theory.Larger angles of swept tend to decrease the spanwise spacing of the crossflow vortices,which correspondingly helps the stimulation of 'z' mode.展开更多
We investigate initial-boundary-value problem for three-dimensional magnetohydrodynamic （MHD） system of compressible viscous heat-conductive flows and the three-dimensional full compressible Navier-Stokes equations....We investigate initial-boundary-value problem for three-dimensional magnetohydrodynamic （MHD） system of compressible viscous heat-conductive flows and the three-dimensional full compressible Navier-Stokes equations. We establish a blowup criterion only in terms of the derivative of velocity field, similar to the Beale^Kato-Majda type criterion for compressible viscous barotropic flows by Huang et al. （2011）. The results indicate that the nature of the blowup for compressible MHD models of viscous media is similar to the barotropic compressible Navier-Stokes equations and does not depend on further sophistication of the MHD model, in particular, it is independent of the temperature and magnetic field. It also reveals that the deformation tensor of the velocity field plays a more dominant role than the electromagnetic field and the temperature in regularity theory. Especially, the similar results also hold for compressible viscous heat-conductive Navier-Stokes flows, which extend the results established by Fan et al. （2010）, and I-Iuang and Li （2009）. In addition, the viscous coefficients are only restricted by the physical conditions in this paper.展开更多
The authors study the Rayleigh-Taylor instability for two incompressible immiscible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian c...The authors study the Rayleigh-Taylor instability for two incompressible immiscible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian coordinates. To deal with the free surface, instead of using the transformation to Lagrangian coordinates, the perturbed equations in Eulerian coordinates are transformed to an integral form and the two-fluid flow is formulated as a single-fluid flow in a fixed domain, thus offering an alternative approach to deal with the jump conditions at the free interface. First, the linearized problem around the steady state which describes a denser immiscible fluid lying above a light one with a free interface separating the two fluids, both fluids being in(unstable) equilibrium is analyzed. By a general method of studying a family of modes, the smooth(when restricted to each fluid domain) solutions to the linearized problem that grow exponentially fast in time in Sobolev spaces are constructed, thus leading to a global instability result for the linearized problem.Then, by using these pathological solutions, the global instability for the corresponding nonlinear problem in an appropriate sense is demonstrated.展开更多
The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to ...The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to understand the detailed flow physics associated with the twin impulse wave, compared with those in a single impulse wave. In the current study, the merging phenomena and propagation characteristics of the impulse waves are investigated using a shock tube experiment and by numerical computations. The Harten-Yee''s total variation diminishing (TVD) scheme is used to solve the unsteady two-dimensional compressible Euler equations. The Mach number Ms of incident shock wave is changed below 1.5 and the distance between two-parallel tubes, L/d, is changed from 1.2 to 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. The results obtained show that on the symmetric axis between two-parallel tubes, the peak pressure produced by the twin impulse waves and its location strongly depend upon the distance between two-parallel tubes, L/d and the incident shock Mach number, Ms. The predicted Schlieren images represent the measured twin-impulse wave with a good accuracy.展开更多
Consider the systemwhich can be used to model the adiabatic gas flow through porous media. Here v is specific volume, u denotes velocity, s stands for entropy, p denotes pressure with pv <0 for v >0. It is prove...Consider the systemwhich can be used to model the adiabatic gas flow through porous media. Here v is specific volume, u denotes velocity, s stands for entropy, p denotes pressure with pv <0 for v >0. It is proved that the solutions of (1) tend to those of the following nonlinear parabolic equation time-asymptotically:展开更多
基金Supported by Marie Curie International Incoming Fellowship (No. PIIF-GA-2009-253453)
文摘A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.
基金Supported by the National Natural Science Foundation of China （10472024）.
文摘The idea that the collapse proceeds from the outer boundary of the cavity cloud towards its center for the ultrasonic cavitation proposed by Hasson and Morch in 1980s is further developed for calculating the collapse pressure and boundaries of cavity cloud at the collapse stage of bubbles for hydraulic cavitation flow in Venturi in present research. The numerical simulation is carried out based on Gilmore＇s eouations of bubble dynamics, which take account of the compressibility of fluid besides the viscosity and interfacial tension. The collapse of the cavity cloud is considered to proceed layer by layer from the outer cloud towards its inner part. The simulation results indicate that thepredicted boundaries of the cavity cloudat the collapse stage agree.well with the exPerimental ones.It is also found that the maximum collapse pressure of the cavity cloud is several times as high as the collapse pressure of outside boundary, and it is located at a point in the axis, where the cavity cloud disappears completely. This means that a cavity cloud has higher collapse pressure or strength than that of a single bubble due to the interactions of the bubbles. The effects of operation and structural parameters on the collapse pressure are also analyzed in detail.
基金Project (No. 50236030) supported by the National Natural Science Foundation of China
文摘Vortex methods have been alternative tools of finite element and finite difference methods for several decades. This paper presents a brief review of vortex method development in the last decades and introduces efficient vortex methods developed for high Reynolds number bluff body flows and suitable for running on parallel computer architectures. Included in this study are particle strength exchange methods, core-spreading method, deterministic particle method and hybrid vortex methods. Combined with conservative methods, vortex methods can comprise the most available tools for simulations of three-dimensional complex bluff body flows at high Reynolds numbers.
基金Supported by National 863 Plan Project of Ministry of Science and Technology of China under Grant No. 2006AA09Z354National Natural Science Foundation of China under Grant No. 10672101.
文摘In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.
基金Supported by National Natural Science Foundations of China under Grant Nos.10735030,10475055,10675065,and 90503006National Basic Research Program of China (973 Program) under Grant No.2007CB814800+2 种基金PCSIRT (IRT0734)the Research Fund of Postdoctoral of China under Grant No.20070410727Specialized Research Fund for the Doctoral Program of Higher Education under Grant No.20070248120
文摘We investigate the symmetry reduction for the two-dimensional incompressible Navier-Stokes equationin conventional stream function form through Lie symmetry method and construct some similarity reduction solutions.Two special cases in [D.K.Ludlow,P.A.Clarkson,and A.P.Bassom,Stud.Appl.Math.103 (1999) 183] and a theoremin [S.Y.Lou,M.Jia,X.Y.Tang,and F.Huang,Phys.Rev.E 75 (2007) 056318] are retrieved.
文摘The temperature separation was discovered inside the short vortex chamber （H/D = 0.18）. Experiments revealed that the highest temperature of the periphery was 465 ℃, and the lowest temperature of the central zone was -45 ℃ （the compressed air was pumped into the chamber at room temperature）. The objective of this paper is to proof that this temperature separation effect cannot be explained by conventional heat transfer processes. To explain this phenomenon, the concept of PGEW （Pressure Gradient Elastic Waves） is proposed. PGEW are kind of elastic waves, which operate in compressible fluids with pressure gradients and density fluctuations. The result of PGEW propagation is a heat transfer from area of low pressure to high pressure zone. The physical model of a gas in a strong field of mass forces is proposed to substantiate the PGEW existence. This physical model is intended for the construction of a theory of PGEW. Understanding the processes associated with the PGEW permits the possibility of creating new devices for energy saving and low potential heat utilization, which have unique properties.
基金Supported by the National Basic Research Program of China under Grant No.2007CB815100the Research Fund for the Doctoral Program of Higher Education of China under Grant No.20070290008the National Natural Science Foundation of China under Grant Nos.10775020 and 10874242
文摘This paper studies the phase effect in mode coupling of Kelvin-Helmholtz instability in two-dimensionalincompressible fluid.It is found that there is an important growth phenomenon of every mode in the mode couplingprocess.The growth changes periodically with phase difference and in the condition of our simulation the period is about0.7π.The period characteristic is apparent in all stage of the mode coupling process,especially in the relatively laterstage.
基金Supported by the Science Foundation of Laboratory of Computational Physics, Science Foundation of China Academy of Engineering Physics under Grant Nos. 2009A0102005, 2009B0101012National Basic Research Program of China under Grant No. 2007CB815105+1 种基金National Natural Science Foundation of China under Grant Nos. 11074300, 11075021, and 11074303the Fundamental Research Funds for the Central Universities under Grant No. 2010YS03
文摘In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann （LB） model for shacked compressible flows. The proposed flux limiter scheme is efficient in decreasing the artificial oscillations and numerical diffusion around the interface. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. Numerical simulations for the Richtmyer-Meshkov instability show that with the new model the computed interfaces are smoother and more consistent with physical analysis. The growth rates of bubble and spike present a satisfying agreement with the theoretical predictions and other numerical simulations.
文摘The equations of motion of a bubble, expanding adiabatically through an incompressible viscous fluid, are deduced when the centre of the bubble moves in a vertical plane in the presence of gravitational acceleration, acting vertically downwards. The non-linear equations of motion obtained are solved numerically for different values of the various parameters of the problem. The path traced by the centre of the bubble and velocity of the centre, the change of radius R with time, and the influence of the buoyancy force, which is experienced by the expanding bubble for different values of the gravitational acceleration on these quantities, are investigated. The radius R（t） of the bubble is found to vary periodically with time when the acceleration due to gravity is small. But when the acceleration due to gravity increases, this periodicity in the value of R（t） with t is lost. The influence of viscosity in determining the periodicity of the bubble motion is also investigated.
基金supported by National Natural Science Foundation of China (Grant Nos. 11101044,11271051,11229101 and 91130020)National Basic Research Program of China (Grant No.2011CB309705)
文摘We investigate the nonlinear instability of a smooth steady density profile solution to the threedimensional nonhomogeneous incompressible Navier-Stokes equations in the presence of a uniform gravitational field,including a Rayleigh-Taylor steady-state solution with heavier density with increasing height(referred to the Rayleigh-Taylor instability).We first analyze the equations obtained from linearization around the steady density profile solution.Then we construct solutions to the linearized problem that grow in time in the Sobolev space H k,thus leading to a global instability result for the linearized problem.With the help of the constructed unstable solutions and an existence theorem of classical solutions to the original nonlinear equations,we can then demonstrate the instability of the nonlinear problem in some sense.Our analysis shows that the third component of the velocity already induces the instability,which is different from the previous known results.
文摘The large eddy simulation (LES) of compressible turbulent channel flows at three different Mach numbers is performed in the present work,by extending the dynamic mixed subgrid-scale (SGS) model to compressible flows.The turbulent statistics agree well with those from the existing direct numerical simulation (DNS) results,indicating that the LES method established in the present work is reliable.The analysis of the turbulent fluctuations computed by the present LES reveals that the flows considered in this work follow the Morkovin's hypothesis.Thus,the compressibility effects are dominated by the mean field properties,and the relevant statistical ratios are invariant to the variation of Mach number.The near-wall streamwise streaks are more coherent and the spacing between streaks is wider as the Mach number increases.This can be regarded as a direct feature characterized by the compressibility effects.The restrained influences of compressibility effects on the production and dissipation of the turbulence kinetic energy are also identified based on the present LES results.
文摘Based on the standard k-ε turbulence model, a new compressible k-ε model considering the pressure expansion influence due to the compressibility of fluid is developed and applied to the simulation of 3D transonic turbulent flows in a nozzle and a cascade. The Reynolds avenged N-S equations in generalized curvilinear coordinates are solved with implementation of the new model. The high resolution TVD scheme is used to discretize the convective terms. The numerical results show that the compressible k-ε model behaves well in the simulation of transonic internal turbulent flows.
基金Project supported by the National Natural Science Foundation of China (No.51078230)the Research Fund for the Doctoral Program of Higher Education of China (No.200802480056)the Key Project of Fund of Science and Technology Development of Shanghai (No.10JC1407900),China
文摘A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) stabilization scheme is used for the formulation of the Navier-Stokes equations. For the spatial discretization, the convection term is treated explicitly, while the viscous term is treated implicitly, and for the temporal discretization, a three-step method is employed. The present method is applied to simulate the lid driven cavity problems with different geometries at low and high Reynolds numbers. The results compared with other numerical experiments are found to be feasible and satisfactory.
文摘Based on the immersed boundary method,a fast simulation for solving unsteady,incompressible,viscous flow associated with the oscillating cascade is established on a quasi-three-dimensional coordinate system.The numerical method is applied to the simulation of the flow passing an oscillating circular cylinder which is forced to move in X direction under prescribed motions in water at rest at low Keulegan-Carpenter numbers.Then vor-tex-induced vibration of a cylinder with two degrees of freedom which oscillates in in-line direction and transverse direction is simulated using this method.The results are in good agreement with the previous research.Then the method is extended to the oscillating cascade simulation of making various comparisons.It is found that the IBPA(inter blade phase angle) will change as the time goes on,because of the non-uniformity of the flow in the circumferential direction,until the oscillating cascade goes to a stable situation.The reduced velocity and the number of blades are chosen to investigate the effects of them on IBPA.The results indicate that both the reduced velocity and the number of blades are the main factors which influence IBPA.It is worth noting that the coupling process is not necessary to generate any body-fitting grids,which makes it much faster in computational process for such a complicated fluid-structure interaction problem.
文摘In this paper, the flow fields of underexpanded impinging jet issued from rectangular nozzles of aspect ratio 1, 3 and 5 are numerically and experimentally studied. Two dimensional temperature and pressure distributions are measured by using infrared camera and the combination of a pressure scanning device and a stepping motor, respectively. The variation of the stagnation pressure on the impinging plate reveals that a hysteretic phenomenon exists during the increasing and decreasing of the pressure ratio for the aspect ratio of 3.0 and 5.0. It is also found that the nozzle of aspect ratio 1.0 caused the largest total pressure loss Pc /p0= 0.27 at the pressure ratio of Po /p0 = 6.5, where Pc is the stagnation center pressure on the wall, P0 the upstream stagnation pressure, Pb the ambient pressure. The other two nozzles showed that the pressure loss Pc / P0 =0.52 and 0.55 were achieved by the nozzles of the aspect ratio 3,0 and 5.0, respectively. The comparison between the calculations and experiments is fairly good, showing the three dimensional streamlines and structures of the shock waves in the jets. However, the hysteresis of the pressure variations observed in the experiments between the pressure ratio of 3.5 and 4.5 cannot be confirmed in the calculations.
基金supported by the National Natural Science Foundation of China(Grant Nos. 90505005 and 10932005)
文摘The crossflow instability of a three-dimensional (3-D) boundary layer is an important factor which affects the transition over a swept-wing.In this report,the primary instability of the incompressible flow over a swept wing is investigated by solving nonlinear parabolized stability equations (NPSE).The Floquet theory is applied to study the dependence of the secondary and high-frequency instabilities on curvature,Reynolds number and angle of swept (AOS).The computational results show that the curvature in the present case has no significant effect on the secondary instabilities.It is generally believed that the secondary instability growth rate increases with the magnitude of the nonlinear mode of crossflow vortex.But,at a certain state,when the Reynolds number is 3.2 million,we find that the secondary instability growth rate becomes smaller even when the magnitude of the nonlinear mode of the crossflow vortex is larger.The effect of the angle of swept at 35,45 and 55 degrees,respectively,is also studied in the framework of the secondary linear stability theory.Larger angles of swept tend to decrease the spanwise spacing of the crossflow vortices,which correspondingly helps the stimulation of 'z' mode.
基金supported by National Natural Science Foundation of China(Grant Nos.11171236 and 71372189)Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT1273)+1 种基金Sichuan Youth Science and Technology Foundation(Grant No.2014JQ0003)China Postdoctoral Science Foundation(Grant No.2013M542285)
文摘We investigate initial-boundary-value problem for three-dimensional magnetohydrodynamic （MHD） system of compressible viscous heat-conductive flows and the three-dimensional full compressible Navier-Stokes equations. We establish a blowup criterion only in terms of the derivative of velocity field, similar to the Beale^Kato-Majda type criterion for compressible viscous barotropic flows by Huang et al. （2011）. The results indicate that the nature of the blowup for compressible MHD models of viscous media is similar to the barotropic compressible Navier-Stokes equations and does not depend on further sophistication of the MHD model, in particular, it is independent of the temperature and magnetic field. It also reveals that the deformation tensor of the velocity field plays a more dominant role than the electromagnetic field and the temperature in regularity theory. Especially, the similar results also hold for compressible viscous heat-conductive Navier-Stokes flows, which extend the results established by Fan et al. （2010）, and I-Iuang and Li （2009）. In addition, the viscous coefficients are only restricted by the physical conditions in this paper.
基金supported by the National Natural Science Foundation of China(Nos.11101044,11271051,11229101,11301083,11371065,11471134)the Fujian Provincial Natural Science Foundation of China(No.2014J01011)+1 种基金the National Basic Research Program(No.2011CB309705)the Beijing Center for Mathematics and Information Interdisciplinary Sciences
文摘The authors study the Rayleigh-Taylor instability for two incompressible immiscible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian coordinates. To deal with the free surface, instead of using the transformation to Lagrangian coordinates, the perturbed equations in Eulerian coordinates are transformed to an integral form and the two-fluid flow is formulated as a single-fluid flow in a fixed domain, thus offering an alternative approach to deal with the jump conditions at the free interface. First, the linearized problem around the steady state which describes a denser immiscible fluid lying above a light one with a free interface separating the two fluids, both fluids being in(unstable) equilibrium is analyzed. By a general method of studying a family of modes, the smooth(when restricted to each fluid domain) solutions to the linearized problem that grow exponentially fast in time in Sobolev spaces are constructed, thus leading to a global instability result for the linearized problem.Then, by using these pathological solutions, the global instability for the corresponding nonlinear problem in an appropriate sense is demonstrated.
文摘The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to understand the detailed flow physics associated with the twin impulse wave, compared with those in a single impulse wave. In the current study, the merging phenomena and propagation characteristics of the impulse waves are investigated using a shock tube experiment and by numerical computations. The Harten-Yee''s total variation diminishing (TVD) scheme is used to solve the unsteady two-dimensional compressible Euler equations. The Mach number Ms of incident shock wave is changed below 1.5 and the distance between two-parallel tubes, L/d, is changed from 1.2 to 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. The results obtained show that on the symmetric axis between two-parallel tubes, the peak pressure produced by the twin impulse waves and its location strongly depend upon the distance between two-parallel tubes, L/d and the incident shock Mach number, Ms. The predicted Schlieren images represent the measured twin-impulse wave with a good accuracy.
文摘Consider the systemwhich can be used to model the adiabatic gas flow through porous media. Here v is specific volume, u denotes velocity, s stands for entropy, p denotes pressure with pv <0 for v >0. It is proved that the solutions of (1) tend to those of the following nonlinear parabolic equation time-asymptotically:
