This paper presents numerical simulations of viscous flow past a submarine model in steady turn by solving the Reynolds-Averaged Navier-Stokes Equations(RANSE) for incompressible, steady flows. The rotating coordina...This paper presents numerical simulations of viscous flow past a submarine model in steady turn by solving the Reynolds-Averaged Navier-Stokes Equations(RANSE) for incompressible, steady flows. The rotating coordinate system was adopted to deal with the rotation problem. The Coriolis force and centrifugal force due to the computation in a bodyfixed rotating frame of reference were treated explicitly and added to momentum equations as source terms. Furthermore, velocities of entrances were coded to give the correct magnitude and direction needed. Two turbulence closure models(TCMs), the RNG k-ε model with wall functions and curvature correction and the Shear Stress Transport(SST) k-ω model without the use of wall functions, but with curvature correction and low-Re correction were introduced, respectively. Take DARPA SUBOFF model as the test case, a series of drift angle varying between 0° and 16° at a Reynolds number of 6.53×10^6 undergoing rotating arm test simulations were conducted. The computed forces and moment as a function of drift angle during the steady turn are mostly in close agreement with available experimental data. Though the difference between the pressure coefficients around the hull form was observed, they always show the same trend. It was demonstrated that using sufficiently fine grids and advanced turbulence models will lead to accurate prediction of the flow field as well as the forces and moments on the hull.展开更多
This paper is concerned with the numerical solution of turbulent flows on the concave surfaces of spillway dams. Orthogonal curvilinear coordinates are used to deal with the complicated computational region and the ef...This paper is concerned with the numerical solution of turbulent flows on the concave surfaces of spillway dams. Orthogonal curvilinear coordinates are used to deal with the complicated computational region and the effects of streamline curvature on turbulent flows are included. The SIMPLEC procedure has been used for the transformed governing equations in the transformed domain. The comparison between computed results and experimental data shows a satisfactory agreement.展开更多
This paper proposes a vertex-estimation-based, feature-preserving smoothingtechnique for meshes. A robust mesh smoothing operator called mean value coordinates flow isintroduced to modify mean curvature flow and make ...This paper proposes a vertex-estimation-based, feature-preserving smoothingtechnique for meshes. A robust mesh smoothing operator called mean value coordinates flow isintroduced to modify mean curvature flow and make it more stable. Also the paper proposes athree-pass vertex estimation based on bilateral filtering of local neighbors which is transferredfrom image processing settings and a Quasi-Laplacian operation, derived from the standard Laplacianoperator, is performed to increase the smoothness order of the mesh rapidly whilst denoising meshesefficiently, preventing volume shrinkage as well as preserving sharp features of the mesh. Comparedwith previous algorithms, the result shows it is simple, efficient and robust.展开更多
As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the ...As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics(SPH) method is used to model the compressible fluid, the natural coordinate formulation(NCF) and absolute nodal coordinate formulation(ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit(GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant No.51179199)
文摘This paper presents numerical simulations of viscous flow past a submarine model in steady turn by solving the Reynolds-Averaged Navier-Stokes Equations(RANSE) for incompressible, steady flows. The rotating coordinate system was adopted to deal with the rotation problem. The Coriolis force and centrifugal force due to the computation in a bodyfixed rotating frame of reference were treated explicitly and added to momentum equations as source terms. Furthermore, velocities of entrances were coded to give the correct magnitude and direction needed. Two turbulence closure models(TCMs), the RNG k-ε model with wall functions and curvature correction and the Shear Stress Transport(SST) k-ω model without the use of wall functions, but with curvature correction and low-Re correction were introduced, respectively. Take DARPA SUBOFF model as the test case, a series of drift angle varying between 0° and 16° at a Reynolds number of 6.53×10^6 undergoing rotating arm test simulations were conducted. The computed forces and moment as a function of drift angle during the steady turn are mostly in close agreement with available experimental data. Though the difference between the pressure coefficients around the hull form was observed, they always show the same trend. It was demonstrated that using sufficiently fine grids and advanced turbulence models will lead to accurate prediction of the flow field as well as the forces and moments on the hull.
文摘This paper is concerned with the numerical solution of turbulent flows on the concave surfaces of spillway dams. Orthogonal curvilinear coordinates are used to deal with the complicated computational region and the effects of streamline curvature on turbulent flows are included. The SIMPLEC procedure has been used for the transformed governing equations in the transformed domain. The comparison between computed results and experimental data shows a satisfactory agreement.
文摘This paper proposes a vertex-estimation-based, feature-preserving smoothingtechnique for meshes. A robust mesh smoothing operator called mean value coordinates flow isintroduced to modify mean curvature flow and make it more stable. Also the paper proposes athree-pass vertex estimation based on bilateral filtering of local neighbors which is transferredfrom image processing settings and a Quasi-Laplacian operation, derived from the standard Laplacianoperator, is performed to increase the smoothness order of the mesh rapidly whilst denoising meshesefficiently, preventing volume shrinkage as well as preserving sharp features of the mesh. Comparedwith previous algorithms, the result shows it is simple, efficient and robust.
基金supported by the 111 China Project(Grant No.B16003)the National Natural Science Foundation of China(Grant Nos.11290151,11702022,and 11221202)
文摘As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics(SPH) method is used to model the compressible fluid, the natural coordinate formulation(NCF) and absolute nodal coordinate formulation(ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit(GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.