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.展开更多
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.展开更多
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.展开更多
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.展开更多
The approximate analytical solution of velocity is presented for incompressible and viscous fluid driven by the oscillation of the periodic pressure, between two slit parallel plates with corrugated walls by employing...The approximate analytical solution of velocity is presented for incompressible and viscous fluid driven by the oscillation of the periodic pressure, between two slit parallel plates with corrugated walls by employing perturbation method. The corrugations of the two walls are described as periodic sinusoidal waves with small amplitude either in phase or half-period out of phase. Based on the analysis, we discuss the influence of the dimensionless parameters on velocity u±and mean velocity parameter φ±numerically, such as Reynolds number Re, nondimensional amplitude A of pressure gradient and wave number k.展开更多
基金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.
基金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.
文摘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.
基金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.
基金Supported by the National Natural Science Foundation of China under Grant No.11472140the Natural Science Foundation of Inner Mongolia Autonomous Region of China under Grant No.2016MS0106the Inner Mongolia Grassland Talent under Grant No.12000-12102013
文摘The approximate analytical solution of velocity is presented for incompressible and viscous fluid driven by the oscillation of the periodic pressure, between two slit parallel plates with corrugated walls by employing perturbation method. The corrugations of the two walls are described as periodic sinusoidal waves with small amplitude either in phase or half-period out of phase. Based on the analysis, we discuss the influence of the dimensionless parameters on velocity u±and mean velocity parameter φ±numerically, such as Reynolds number Re, nondimensional amplitude A of pressure gradient and wave number k.
