The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA00...The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA0012 airfoil is used to verify the efficiency of the proposed method. Two cases of the low Mach number flows around the multi-element airfoil and the circular cylinder are also used to test the proposed method. Numerical results show that the methods combined the preconditioning method and compressible Navier-Stokes equations are efficient to solve low Mach number flows.展开更多
In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed...In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.展开更多
In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
Central discontinuous Galerkin（CDG）method is used to solve the Navier-Stokes equations for viscous flow in this paper.The CDG method involves two pieces of approximate solutions defined on overlapping meshes.Taking ...Central discontinuous Galerkin（CDG）method is used to solve the Navier-Stokes equations for viscous flow in this paper.The CDG method involves two pieces of approximate solutions defined on overlapping meshes.Taking advantages of the redundant representation of the solution on the overlapping meshes,the cell interface of one computational mesh is right inside the staggered mesh,hence approximate Riemann solvers are not needed at cell interfaces.Third order total variation diminishing（TVD）Runge-Kutta（RK）methods are applied in time discretization.Numerical examples for 1D and2 D viscous flow simulations are presented to validate the accuracy and robustness of the CDG method.展开更多
This article is concerned with the time periodic solution to the isentropic compressible Navier-Stokes equations in a periodic domain. Using an approach of parabolic regularization, we first obtain the existence of th...This article is concerned with the time periodic solution to the isentropic compressible Navier-Stokes equations in a periodic domain. Using an approach of parabolic regularization, we first obtain the existence of the time periodic solution to a regularized problem under some smallness and symmetry assumptions on the external force. The result for the original compressible Navier-Stokes equations is then obtained by a limiting process. The uniqueness of the periodic solution is also given.展开更多
In this article, we prove the existence and uniqueness of solutions of the NavierStokes equations with Navier slip boundary condition for incompressible fluid in a bounded domain of R^3. The results are established by...In this article, we prove the existence and uniqueness of solutions of the NavierStokes equations with Navier slip boundary condition for incompressible fluid in a bounded domain of R^3. The results are established by the Galerkin approximation method and improved the existing results.展开更多
Although the two-grid finite element decoupled scheme for mixed Navier-Stokes/ Darcy model in literatures has given the numerical results of optimal convergence order, the theoretical analysis only obtain the optimal ...Although the two-grid finite element decoupled scheme for mixed Navier-Stokes/ Darcy model in literatures has given the numerical results of optimal convergence order, the theoretical analysis only obtain the optimal error order for the porous media flow and a non-optimal error order for the fluid flow. In this article, we give a more rigorous of the error analysis for the fluid flow and obtain the optimal error estimates of the velocity and the pressure.展开更多
Hybrid grids are used for the solution of 3D turbulence Navier Stokes equations. The prismatic grids are generated near the wall, and the tetrahedron grids are generated in the other field. A Navier Stokes solver usin...Hybrid grids are used for the solution of 3D turbulence Navier Stokes equations. The prismatic grids are generated near the wall, and the tetrahedron grids are generated in the other field. A Navier Stokes solver using Jamson′s finite volume method is developed. The algebraic Baldwin Lomax turbulence model is adopted. The numerical tests show that the above method is very efficient.展开更多
We are concerned with the Cauchy problem regarding the full compressible Navier-Stokes equations in R^(d)(d=2,3).By exploiting the intrinsic structure of the equations and using harmonic analysis tools(especially the ...We are concerned with the Cauchy problem regarding the full compressible Navier-Stokes equations in R^(d)(d=2,3).By exploiting the intrinsic structure of the equations and using harmonic analysis tools(especially the Littlewood-Paley theory),we prove the global solutions to this system with small initial data restricted in the Sobolev spaces.Moreover,the initial temperature may vanish at infinity.展开更多
This paper considers the Calerkin finite element method for the incompressible Navier-Stokes equations in two dimensions, where the finite-dimensional spaces employed consist of piecewise polynomials enriched with res...This paper considers the Calerkin finite element method for the incompressible Navier-Stokes equations in two dimensions, where the finite-dimensional spaces employed consist of piecewise polynomials enriched with residual-free bubble （RFB） functions. The stability features of the residual-free bubble functions for the linearized Navier-Stokes equations are analyzed in this work. It is shown that the enrichment of the velocity space by bubble functions stabilizes the numerical method for any value of the viscosity parameter for triangular elements and for values of the viscosity parameter in the vanishing limit case for quadrilateral elements.展开更多
Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous f...Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.展开更多
A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as th...A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.展开更多
There are many ways of describing a solid,porous or fluid region of the computational domain when solving the Navier-Stokes equations(NSE)for flow motions.Amongst these the porous cell method is one of the most flexib...There are many ways of describing a solid,porous or fluid region of the computational domain when solving the Navier-Stokes equations(NSE)for flow motions.Amongst these the porous cell method is one of the most flexible approaches.In this method,a parameter is defined as a ratio of the volume open to water and air in a calculation cell to its cell volume.In the calculation,the same numerical procedure is applied to every cell and no explicit boundary conditions are needed at solid boundaries.The method is used to simulate flow through porous media,around solid bodies and over a moving seabed.The results compare well with experimental data and other numerical results.In our future work the porous cell method will be applied to more complex fluid-solid interaction situations.展开更多
Generally the incompressible viscous flow problem is described by the Navier-Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the s...Generally the incompressible viscous flow problem is described by the Navier-Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the step-bystep computation in the field of time, and adopting the least-square estimation of the-same-order shift, this paper has calculated both velocity and pressure from the decoupling independent equations. Each time fraction Newton-Raphson iterative method is applied for the velocity and pressure. Finally, this paper puts the method into practice of the shear-drive cavity flow, verifying the validity, high accuracy and stability.展开更多
The discontinuous Galerkin(DG) method is established and innovatively conducted on accurately simulating the evolution of blade-tip vortex and the aerodynamic characteristics of helicopter rotor. Firstly,the Reynolds-...The discontinuous Galerkin(DG) method is established and innovatively conducted on accurately simulating the evolution of blade-tip vortex and the aerodynamic characteristics of helicopter rotor. Firstly,the Reynolds-Averaged Navier-Stokes(RANS)equations in rotating reference frame are employed,and the embedded grid system is developed with the finite volume method(FVM)and the DG method conducted on the blade grid and background grid respectively. Besides,the Harten-Lax-Van Leer contact(HLLC)scheme with high-resolution and low-dissipation is employed for spatial discretization,and the explicit third-order Runge-Kutta scheme is used to accomplish the temporal discretization. Secondly,the aerodynamic characteristics and the evolution of blade-tip vortex for Caradonna-Tung rotor are simulated by the established CFD method,and the numerical results are in good agreement with experimental data,which well validates the accuracy of the DG method and shows the advantages of DG method on capturing the detailed blade-tip vortex compared with the FVM method. Finally,the evolution of tip vortex at different blade tip Mach numbers and collective pitches is discussed.展开更多
An unsteady load calculation method for the support configuration of a monopile-supported offshore wind turbine is developed based on the Fluent software platform.Firstly,the water wave is generated by imposing the in...An unsteady load calculation method for the support configuration of a monopile-supported offshore wind turbine is developed based on the Fluent software platform.Firstly,the water wave is generated by imposing the inlet boundary conditions according to the exact potential flow solution.Then the wave evolution is simulated by solving the unsteady incompressible Navier-Stokes(N-S)equations coupled with the volume of fluid method.For the small amplitude wave with reasonable wave parameters,the numerical wave result agrees well with that of the given wave model.Finally,a monopile support configuration is introduced and a CFD-based load calculation method is established to accurately calculate the unsteady load under the combined action of wave and wind.The computed unsteady wave load on a small-size monopile support located in the small amplitude wave flow coincides with that of the Morison formula.The load calculations are also performed on a large-size monopile support and a monopile-supported offshore wind turbine under the combined action of small amplitude wave and wind.展开更多
In this study,the performance of a twin-screw propeller under the influence of the wake field of a fully appended ship was investigated using a coupled Reynolds-averaged Navier–Stokes(RANS)/boundary element method(BE...In this study,the performance of a twin-screw propeller under the influence of the wake field of a fully appended ship was investigated using a coupled Reynolds-averaged Navier–Stokes(RANS)/boundary element method(BEM)code.The unsteady BEM is an efficient approach to predicting propeller performance.By applying the time-stepping method in the BEM solver,the trailing vortex sheet pattern of the propeller can be accurately captured at each time step.This is the main innovation of the coupled strategy.Furthermore,to ascertain the effect of the wake field of the ship with acceptable accuracy,a RANS solver was developed.A finite volume method was used to discretize the Navier–Stokes equations on fully unstructured grids.To simulate ship motions,the volume of the fluid method was applied to the RANS solver.The validation of each solver(BEM/RANS)was separately performed,and the results were compared with experimental data.Ultimately,the BEM and RANS solvers were coupled to estimate the performance of a twin-screw propeller,which was affected by the wake field of the fully appended hull.The proposed model was applied to a twin-screw oceanography research vessel.The results demonstrated that the presented model can estimate the thrust coefficient of a propeller with good accuracy as compared to an experimental self-propulsion test.The wake sheet pattern of the propeller in open water(uniform flow)was also compared with the propeller in a real wake field.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
Constructive interference between tidal stream turbines in multi-rotor fence configurations arrayed normally to the flow has been shown analytically, computationally, and experimentally to enhance turbine performance....Constructive interference between tidal stream turbines in multi-rotor fence configurations arrayed normally to the flow has been shown analytically, computationally, and experimentally to enhance turbine performance. The increased resistance to bypass flow due to the presence of neighbouring turbines allows a static pressure difference to develop in the channel and entrains a greater flow rate through the rotor swept area. Exploiting the potential improvement in turbine performance requires that turbines either be operated at higher tip speed ratios or that turbines are redesigned in order to increase thrust. Recent studies have demonstrated that multi-scale flow dynamics, in which a distinction is made between device-scale and fence-scale flow events, have an important role in the physics of flow past tidal turbine fences partially spanning larger channels. Although the reduction in flow rate through the fence as the turbine thrust level increases has been previously demonstrated, the within-fence variation in turbine performance, and the consequences for overall farm performance, is less well understood. The impact of turbine design and operating conditions, on the performance of a multi-rotor tidal fence is investigated using Reynolds-Averaged Navier-Stokes embedded blade element actuator disk simulations. Fences consisting of four, six, and eight turbines are simulated, and it is demonstrated that the combination of device-and fence-scale flow effects gives rise to cross-fence thrust and power variation. These cross-fence variations are also a function of turbine thrust, and hence design conditions,although it is shown simple turbine control strategies can be adopted in order to reduce the cross-fence variations and improve overall fence performance. As the number of turbines in the fence, and hence fence length, increases, it is shown that the turbines may be designed or operated to achieve higher thrust levels than if the turbines were not deployed in a fence configuration.展开更多
A discontinuous Galerkin(DG)-based lattice Boltzmann method is employed to solve the Euler and Navier-Stokes equations.Instead of adopting the widely used local Lax-Friedrichs flux and Roe Flux etc.,a hybrid lattice B...A discontinuous Galerkin(DG)-based lattice Boltzmann method is employed to solve the Euler and Navier-Stokes equations.Instead of adopting the widely used local Lax-Friedrichs flux and Roe Flux etc.,a hybrid lattice Boltzmann flux solver(LBFS)is employed to evaluate the inviscid flux across the cell interfaces.The main advantage of the hybrid LBFS is its flexibility for capturing both strong shocks and thin boundary layers through introducing a function which varies from zero to one to control the artificial viscosity.Numerical results indicate that the hybrid lattice Boltzmann flux solver behaves very well combining with the high-order DG method when simulating both inviscid and viscous flows.展开更多
文摘The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA0012 airfoil is used to verify the efficiency of the proposed method. Two cases of the low Mach number flows around the multi-element airfoil and the circular cylinder are also used to test the proposed method. Numerical results show that the methods combined the preconditioning method and compressible Navier-Stokes equations are efficient to solve low Mach number flows.
文摘In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
基金Supported by the National Natural Science Foundation of China(11602262)
文摘Central discontinuous Galerkin（CDG）method is used to solve the Navier-Stokes equations for viscous flow in this paper.The CDG method involves two pieces of approximate solutions defined on overlapping meshes.Taking advantages of the redundant representation of the solution on the overlapping meshes,the cell interface of one computational mesh is right inside the staggered mesh,hence approximate Riemann solvers are not needed at cell interfaces.Third order total variation diminishing（TVD）Runge-Kutta（RK）methods are applied in time discretization.Numerical examples for 1D and2 D viscous flow simulations are presented to validate the accuracy and robustness of the CDG method.
基金supported by the Program for New Century Excellent Talents in University of the Ministry of Education(NCET-13-0804)NSFC(11471127)+3 种基金Guangdong Natural Science Funds for Distinguished Young Scholar(2015A030306029)The Excellent Young Teachers Program of Guangdong Province(HS2015007)Pearl River S&T Nova Program of Guangzhou(2013J2200064)supported by the General Research Fund of Hong Kong,City U 104511
文摘This article is concerned with the time periodic solution to the isentropic compressible Navier-Stokes equations in a periodic domain. Using an approach of parabolic regularization, we first obtain the existence of the time periodic solution to a regularized problem under some smallness and symmetry assumptions on the external force. The result for the original compressible Navier-Stokes equations is then obtained by a limiting process. The uniqueness of the periodic solution is also given.
基金supported by National Board for Higher Mathematics(02011/9/2019NBHM(R.P.)/R and D Ⅱ/1324)
文摘In this article, we prove the existence and uniqueness of solutions of the NavierStokes equations with Navier slip boundary condition for incompressible fluid in a bounded domain of R^3. The results are established by the Galerkin approximation method and improved the existing results.
基金Subsidized by NSFC(11571274 and 11171269)the Ph.D.Programs Foundation of Ministry of Education of China(20110201110027)
文摘Although the two-grid finite element decoupled scheme for mixed Navier-Stokes/ Darcy model in literatures has given the numerical results of optimal convergence order, the theoretical analysis only obtain the optimal error order for the porous media flow and a non-optimal error order for the fluid flow. In this article, we give a more rigorous of the error analysis for the fluid flow and obtain the optimal error estimates of the velocity and the pressure.
文摘Hybrid grids are used for the solution of 3D turbulence Navier Stokes equations. The prismatic grids are generated near the wall, and the tetrahedron grids are generated in the other field. A Navier Stokes solver using Jamson′s finite volume method is developed. The algebraic Baldwin Lomax turbulence model is adopted. The numerical tests show that the above method is very efficient.
基金supported by the National Natural Science Foundation of China(11801090 and 12161004)Jiangxi Provincial Natural Science Foundation,China(20212BAB211004)+3 种基金supported by the National Natural Science Foundation of China(12171493)supported by the National Natural Science Foundation of China(11601533)Guangdong Provincial Natural Science Foundation,China(2022A1515011977)the Science and Technology Program of Shenzhen under grant 20200806104726001.
文摘We are concerned with the Cauchy problem regarding the full compressible Navier-Stokes equations in R^(d)(d=2,3).By exploiting the intrinsic structure of the equations and using harmonic analysis tools(especially the Littlewood-Paley theory),we prove the global solutions to this system with small initial data restricted in the Sobolev spaces.Moreover,the initial temperature may vanish at infinity.
文摘This paper considers the Calerkin finite element method for the incompressible Navier-Stokes equations in two dimensions, where the finite-dimensional spaces employed consist of piecewise polynomials enriched with residual-free bubble （RFB） functions. The stability features of the residual-free bubble functions for the linearized Navier-Stokes equations are analyzed in this work. It is shown that the enrichment of the velocity space by bubble functions stabilizes the numerical method for any value of the viscosity parameter for triangular elements and for values of the viscosity parameter in the vanishing limit case for quadrilateral elements.
基金Supported by the National Natural Science Foundation of China(10372043,11172134)the Fundingof Jiangsu Innovation Program for Graduate Education(CXZZ11-0192)~~
文摘Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.
文摘A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.
文摘There are many ways of describing a solid,porous or fluid region of the computational domain when solving the Navier-Stokes equations(NSE)for flow motions.Amongst these the porous cell method is one of the most flexible approaches.In this method,a parameter is defined as a ratio of the volume open to water and air in a calculation cell to its cell volume.In the calculation,the same numerical procedure is applied to every cell and no explicit boundary conditions are needed at solid boundaries.The method is used to simulate flow through porous media,around solid bodies and over a moving seabed.The results compare well with experimental data and other numerical results.In our future work the porous cell method will be applied to more complex fluid-solid interaction situations.
文摘Generally the incompressible viscous flow problem is described by the Navier-Stokes equation. Based on the weighted residual method the discrete formulation of element-free Galerkin is inferred in this paper. By the step-bystep computation in the field of time, and adopting the least-square estimation of the-same-order shift, this paper has calculated both velocity and pressure from the decoupling independent equations. Each time fraction Newton-Raphson iterative method is applied for the velocity and pressure. Finally, this paper puts the method into practice of the shear-drive cavity flow, verifying the validity, high accuracy and stability.
基金supported by the National Natural Science Foundation of China(Nos.12072156, 12032012)the Foundation of Rotor Aerodynamic Key Laboratory (No.RAL20190102)the Priority Academic Program Development Project of Jiangsu Higher Education Institutions(PAPD)。
文摘The discontinuous Galerkin(DG) method is established and innovatively conducted on accurately simulating the evolution of blade-tip vortex and the aerodynamic characteristics of helicopter rotor. Firstly,the Reynolds-Averaged Navier-Stokes(RANS)equations in rotating reference frame are employed,and the embedded grid system is developed with the finite volume method(FVM)and the DG method conducted on the blade grid and background grid respectively. Besides,the Harten-Lax-Van Leer contact(HLLC)scheme with high-resolution and low-dissipation is employed for spatial discretization,and the explicit third-order Runge-Kutta scheme is used to accomplish the temporal discretization. Secondly,the aerodynamic characteristics and the evolution of blade-tip vortex for Caradonna-Tung rotor are simulated by the established CFD method,and the numerical results are in good agreement with experimental data,which well validates the accuracy of the DG method and shows the advantages of DG method on capturing the detailed blade-tip vortex compared with the FVM method. Finally,the evolution of tip vortex at different blade tip Mach numbers and collective pitches is discussed.
基金supported partly by the National Basic Research Program of China("973"Program)(No.2014CB046200)the National Natural Science Foundation of China(No.11372135)the NUAA Fundamental Research Funds(No.NS2013005)
文摘An unsteady load calculation method for the support configuration of a monopile-supported offshore wind turbine is developed based on the Fluent software platform.Firstly,the water wave is generated by imposing the inlet boundary conditions according to the exact potential flow solution.Then the wave evolution is simulated by solving the unsteady incompressible Navier-Stokes(N-S)equations coupled with the volume of fluid method.For the small amplitude wave with reasonable wave parameters,the numerical wave result agrees well with that of the given wave model.Finally,a monopile support configuration is introduced and a CFD-based load calculation method is established to accurately calculate the unsteady load under the combined action of wave and wind.The computed unsteady wave load on a small-size monopile support located in the small amplitude wave flow coincides with that of the Morison formula.The load calculations are also performed on a large-size monopile support and a monopile-supported offshore wind turbine under the combined action of small amplitude wave and wind.
文摘In this study,the performance of a twin-screw propeller under the influence of the wake field of a fully appended ship was investigated using a coupled Reynolds-averaged Navier–Stokes(RANS)/boundary element method(BEM)code.The unsteady BEM is an efficient approach to predicting propeller performance.By applying the time-stepping method in the BEM solver,the trailing vortex sheet pattern of the propeller can be accurately captured at each time step.This is the main innovation of the coupled strategy.Furthermore,to ascertain the effect of the wake field of the ship with acceptable accuracy,a RANS solver was developed.A finite volume method was used to discretize the Navier–Stokes equations on fully unstructured grids.To simulate ship motions,the volume of the fluid method was applied to the RANS solver.The validation of each solver(BEM/RANS)was separately performed,and the results were compared with experimental data.Ultimately,the BEM and RANS solvers were coupled to estimate the performance of a twin-screw propeller,which was affected by the wake field of the fully appended hull.The proposed model was applied to a twin-screw oceanography research vessel.The results demonstrated that the presented model can estimate the thrust coefficient of a propeller with good accuracy as compared to an experimental self-propulsion test.The wake sheet pattern of the propeller in open water(uniform flow)was also compared with the propeller in a real wake field.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
文摘Constructive interference between tidal stream turbines in multi-rotor fence configurations arrayed normally to the flow has been shown analytically, computationally, and experimentally to enhance turbine performance. The increased resistance to bypass flow due to the presence of neighbouring turbines allows a static pressure difference to develop in the channel and entrains a greater flow rate through the rotor swept area. Exploiting the potential improvement in turbine performance requires that turbines either be operated at higher tip speed ratios or that turbines are redesigned in order to increase thrust. Recent studies have demonstrated that multi-scale flow dynamics, in which a distinction is made between device-scale and fence-scale flow events, have an important role in the physics of flow past tidal turbine fences partially spanning larger channels. Although the reduction in flow rate through the fence as the turbine thrust level increases has been previously demonstrated, the within-fence variation in turbine performance, and the consequences for overall farm performance, is less well understood. The impact of turbine design and operating conditions, on the performance of a multi-rotor tidal fence is investigated using Reynolds-Averaged Navier-Stokes embedded blade element actuator disk simulations. Fences consisting of four, six, and eight turbines are simulated, and it is demonstrated that the combination of device-and fence-scale flow effects gives rise to cross-fence thrust and power variation. These cross-fence variations are also a function of turbine thrust, and hence design conditions,although it is shown simple turbine control strategies can be adopted in order to reduce the cross-fence variations and improve overall fence performance. As the number of turbines in the fence, and hence fence length, increases, it is shown that the turbines may be designed or operated to achieve higher thrust levels than if the turbines were not deployed in a fence configuration.
文摘A discontinuous Galerkin(DG)-based lattice Boltzmann method is employed to solve the Euler and Navier-Stokes equations.Instead of adopting the widely used local Lax-Friedrichs flux and Roe Flux etc.,a hybrid lattice Boltzmann flux solver(LBFS)is employed to evaluate the inviscid flux across the cell interfaces.The main advantage of the hybrid LBFS is its flexibility for capturing both strong shocks and thin boundary layers through introducing a function which varies from zero to one to control the artificial viscosity.Numerical results indicate that the hybrid lattice Boltzmann flux solver behaves very well combining with the high-order DG method when simulating both inviscid and viscous flows.
