A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient meth...A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.展开更多
A new favorable iterative algorithm named as PBiCGSTAB （preconditioned bi-conjugate gradient stabilized） algorithm is presented for solving large sparse complex systems. Based on the orthogonal list, the special tec...A new favorable iterative algorithm named as PBiCGSTAB （preconditioned bi-conjugate gradient stabilized） algorithm is presented for solving large sparse complex systems. Based on the orthogonal list, the special technique of only storing non-zero elements is carried out. The incomplete LU factorization without fill-ins is adopted to reduce the condition number of the coefficient matrix. The BiCGSTAB algorithm is extended from the real system to the complex system and it is used to solve the preconditioned complex linear equations. The locked-rotor state of a single-sided linear induction machine is simulated by the software programmed with the finite element method and the PBiCGSTAB algorithm. Then the results are compared with those from the commercial software ANSYS, showing the validation of the proposed software. The iterative steps required for the proposed algorithm are reduced to about one-third, when compared to the BiCG method, therefore the algorithm is fast.展开更多
The restrictively preconditioned conjugate gradient （RPCG） method is further developed to solve large sparse system of linear equations of a block two-by-two structure. The basic idea of this new approach is that we...The restrictively preconditioned conjugate gradient （RPCG） method is further developed to solve large sparse system of linear equations of a block two-by-two structure. The basic idea of this new approach is that we apply the RPCG method to the normal-residual equation of the block two-by-two linear system and construct each required approximate matrix by making use of the incomplete orthogonal factorization of the involved matrix blocks. Numerical experiments show that the new method, called the restrictively preconditioned conjugate gradient on normal residual （RPCGNR）, is more robust and effective than either the known RPCG method or the standard conjugate gradient on normal residual （CGNR） method when being used for solving the large sparse saddle point problems.展开更多
文摘A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.
文摘A new favorable iterative algorithm named as PBiCGSTAB （preconditioned bi-conjugate gradient stabilized） algorithm is presented for solving large sparse complex systems. Based on the orthogonal list, the special technique of only storing non-zero elements is carried out. The incomplete LU factorization without fill-ins is adopted to reduce the condition number of the coefficient matrix. The BiCGSTAB algorithm is extended from the real system to the complex system and it is used to solve the preconditioned complex linear equations. The locked-rotor state of a single-sided linear induction machine is simulated by the software programmed with the finite element method and the PBiCGSTAB algorithm. Then the results are compared with those from the commercial software ANSYS, showing the validation of the proposed software. The iterative steps required for the proposed algorithm are reduced to about one-third, when compared to the BiCG method, therefore the algorithm is fast.
基金supported by the National Basic Research Program (No.2005CB321702)the China NNSF Outstanding Young Scientist Foundation (No.10525102)the National Natural Science Foundation (No.10471146),P.R.China
文摘The restrictively preconditioned conjugate gradient （RPCG） method is further developed to solve large sparse system of linear equations of a block two-by-two structure. The basic idea of this new approach is that we apply the RPCG method to the normal-residual equation of the block two-by-two linear system and construct each required approximate matrix by making use of the incomplete orthogonal factorization of the involved matrix blocks. Numerical experiments show that the new method, called the restrictively preconditioned conjugate gradient on normal residual （RPCGNR）, is more robust and effective than either the known RPCG method or the standard conjugate gradient on normal residual （CGNR） method when being used for solving the large sparse saddle point problems.
