Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompress...Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual (GMRES) method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation.展开更多
The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore...The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.展开更多
Added mass and damping coefficients are very important in hydrodynamic analysis of naval structures. In this paper,a double submerged inclined plates with ‘/\’ configuration is firstly considered. By use of the boun...Added mass and damping coefficients are very important in hydrodynamic analysis of naval structures. In this paper,a double submerged inclined plates with ‘/\’ configuration is firstly considered. By use of the boundary element method(BEM) based on Green function with the wave term, the radiation problem of this special type structure is investigated. The added mass and damping coefficients due to different plate lengths and inclined angles are obtained. The results show that: the added mass and damping coefficients for sway are the largest. Heave is the most sensitive mode to inclined angles. The wave frequencies of the maximal added mass and damping coefficients for sway and roll are the same.展开更多
In this paper, an unstructured, collocated finite volume method for solvingthe Navier-Stokes equations was developed by virtue of auxiliary points. The derivatives weredetermined by the Gauss theorem. The proposed met...In this paper, an unstructured, collocated finite volume method for solvingthe Navier-Stokes equations was developed by virtue of auxiliary points. The derivatives weredetermined by the Gauss theorem. The proposed method could provide control volumes with arbitrarygeometry and preserve the second-order accuracy even if highly distorted grids are used. Althougharbitrary number of cell faces can be used, the hybrid quadrilateral/triangular grids are moredesirable for the simplicity of implementation and applications to engineering problems. Thepressure-velocity coupling was treated using a SIMPLE-like algorithm. The Generalized MinimumResidual (GMRES) method with the Incomplete LU (ILU) preconditioner was used to solve linearequations. Four test cases were studied for validating the proposed method. In using this method,grid quality is not important. Thus, engineers can pay mostly attention to physical mechanism ofproblems. Turbulence models can be simply integrated and the method can be straightforwardlyextended to treat three-dimensional problems.展开更多
Recently, Bal proposed a block-counter-diagonal and a block-counter-triangular precon- ditioning matrices to precondition the GMRES method for solving the structured system of linear equations arising from the Galerki...Recently, Bal proposed a block-counter-diagonal and a block-counter-triangular precon- ditioning matrices to precondition the GMRES method for solving the structured system of linear equations arising from the Galerkin finite-element discretizations of the distributed control problems in (Computing 91 (2011) 379-395). He analyzed the spectral properties and derived explicit expressions of the eigenvalues and eigenvectors of the preconditioned matrices. By applying the special structures and properties of the eigenvector matrices of the preconditioned matrices, we derive upper bounds for the 2-norm condition numbers of the eigenvector matrices and give asymptotic convergence factors of the preconditioned GMRES methods with the block-counter-diagonal and the block-counter-triangular pre- conditioners. Experimental results show that the convergence analyses match well with the numerical results.展开更多
An augmented method is proposed for solving stationary incompressible Stokes equations with a Dirichlet boundary condition along parts of the boundary.In this approach,the normal derivative of the pressure along the p...An augmented method is proposed for solving stationary incompressible Stokes equations with a Dirichlet boundary condition along parts of the boundary.In this approach,the normal derivative of the pressure along the parts of the boundary is introduced as an additional variable and it is solved by the GMRES iterative method.The dimension of the augmented variable in discretization is the number of grid points along the boundary which is O(N).Each GMRES iteration(or one matrix-vector multiplication)requires three fast Poisson solvers for the pressure and the velocity.In our numerical experiments,only a few iterations are needed.We have also combined the augmented approach for Stokes equations involving interfaces,discontinuities,and singularities.展开更多
Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme...Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.展开更多
基金supported by the National Natural Science Foundation of China (No. 50839003)the Natural Science Foundation of Yunnan Province (No. 2008GA027)
文摘Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual (GMRES) method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation.
基金supported by the National Natural Science Foundation of China(Nos. 41130418 and 40925013)the National Basic Research Program(973 Program)(No.2009CB219403)
文摘The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.
基金financially supported by the National Key Basic Research Program of China(Grant No.2013CB036101)the National Natural Science Foundation of China(Grant No.51379037)
文摘Added mass and damping coefficients are very important in hydrodynamic analysis of naval structures. In this paper,a double submerged inclined plates with ‘/\’ configuration is firstly considered. By use of the boundary element method(BEM) based on Green function with the wave term, the radiation problem of this special type structure is investigated. The added mass and damping coefficients due to different plate lengths and inclined angles are obtained. The results show that: the added mass and damping coefficients for sway are the largest. Heave is the most sensitive mode to inclined angles. The wave frequencies of the maximal added mass and damping coefficients for sway and roll are the same.
文摘In this paper, an unstructured, collocated finite volume method for solvingthe Navier-Stokes equations was developed by virtue of auxiliary points. The derivatives weredetermined by the Gauss theorem. The proposed method could provide control volumes with arbitrarygeometry and preserve the second-order accuracy even if highly distorted grids are used. Althougharbitrary number of cell faces can be used, the hybrid quadrilateral/triangular grids are moredesirable for the simplicity of implementation and applications to engineering problems. Thepressure-velocity coupling was treated using a SIMPLE-like algorithm. The Generalized MinimumResidual (GMRES) method with the Incomplete LU (ILU) preconditioner was used to solve linearequations. Four test cases were studied for validating the proposed method. In using this method,grid quality is not important. Thus, engineers can pay mostly attention to physical mechanism ofproblems. Turbulence models can be simply integrated and the method can be straightforwardlyextended to treat three-dimensional problems.
文摘Recently, Bal proposed a block-counter-diagonal and a block-counter-triangular precon- ditioning matrices to precondition the GMRES method for solving the structured system of linear equations arising from the Galerkin finite-element discretizations of the distributed control problems in (Computing 91 (2011) 379-395). He analyzed the spectral properties and derived explicit expressions of the eigenvalues and eigenvectors of the preconditioned matrices. By applying the special structures and properties of the eigenvector matrices of the preconditioned matrices, we derive upper bounds for the 2-norm condition numbers of the eigenvector matrices and give asymptotic convergence factors of the preconditioned GMRES methods with the block-counter-diagonal and the block-counter-triangular pre- conditioners. Experimental results show that the convergence analyses match well with the numerical results.
基金USA NSF-NIH grant#0201094,the USA-ARO under grant number 43751-MAan NSF grant DMS-0412654。
文摘An augmented method is proposed for solving stationary incompressible Stokes equations with a Dirichlet boundary condition along parts of the boundary.In this approach,the normal derivative of the pressure along the parts of the boundary is introduced as an additional variable and it is solved by the GMRES iterative method.The dimension of the augmented variable in discretization is the number of grid points along the boundary which is O(N).Each GMRES iteration(or one matrix-vector multiplication)requires three fast Poisson solvers for the pressure and the velocity.In our numerical experiments,only a few iterations are needed.We have also combined the augmented approach for Stokes equations involving interfaces,discontinuities,and singularities.
文摘Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.