Based on the numerical governing formulation and non-linear complementary conditions of contact and impact problems, a reduced projection augmented Lagrange bi- conjugate gradient method is proposed for contact and im...Based on the numerical governing formulation and non-linear complementary conditions of contact and impact problems, a reduced projection augmented Lagrange bi- conjugate gradient method is proposed for contact and impact problems by translating non-linear complementary conditions into equivalent formulation of non-linear program- ming. For contact-impact problems, a larger time-step can be adopted arriving at numer- ical convergence compared with penalty method. By establishment of the impact-contact formulations which are equivalent with original non-linear complementary conditions, a reduced projection augmented Lagrange bi-conjugate gradient method is deduced to im- prove precision and efficiency of numerical solutions. A numerical example shows that the algorithm we suggested is valid and exact.展开更多
The bi-conjugate gradients(Bi-CG)and bi-conjugate residual(Bi-CR)methods are powerful tools for solving nonsymmetric linear systems Ax=b.By using Kronecker product and vectorization operator,this paper develops the Bi...The bi-conjugate gradients(Bi-CG)and bi-conjugate residual(Bi-CR)methods are powerful tools for solving nonsymmetric linear systems Ax=b.By using Kronecker product and vectorization operator,this paper develops the Bi-CG and Bi-CR methods for the solution of the generalized Sylvester-transpose matrix equationp i=1(Ai X Bi+Ci XTDi)=E(including Lyapunov,Sylvester and Sylvester-transpose matrix equations as special cases).Numerical results validate that the proposed algorithms are much more efcient than some existing algorithms.展开更多
We develop a fast stochastic Galerkin method for an optimal control problem governed by a random space-fractional diffusion equation with deterministic constrained control. Optimal control problems governed by a fract...We develop a fast stochastic Galerkin method for an optimal control problem governed by a random space-fractional diffusion equation with deterministic constrained control. Optimal control problems governed by a fractional diffusion equation tends to provide a better description for transport or conduction processes in heterogeneous media. Howev- er, the fractional control problem introduces significant computation complexity due to the nonlocal nature of fractional differential operators, and this is further worsen by the large number of random space dimensions to discretize the probability space. We ap- proximate the optimality system by a gradient algorithm combined with the stochastic Galerkin method through the discretization with respect to both the spatial space and the probability space. The resulting linear system can be decoupled for the random and spatial variable, and thus solved separately. A fast preconditioned Bi-Conjugate Gradient Stabilized method is developed to efficiently solve the decoupled systems derived from the fractional diffusion operators in the spatial space. Numerical experiments show the utility of the method.展开更多
文摘Based on the numerical governing formulation and non-linear complementary conditions of contact and impact problems, a reduced projection augmented Lagrange bi- conjugate gradient method is proposed for contact and impact problems by translating non-linear complementary conditions into equivalent formulation of non-linear program- ming. For contact-impact problems, a larger time-step can be adopted arriving at numer- ical convergence compared with penalty method. By establishment of the impact-contact formulations which are equivalent with original non-linear complementary conditions, a reduced projection augmented Lagrange bi-conjugate gradient method is deduced to im- prove precision and efficiency of numerical solutions. A numerical example shows that the algorithm we suggested is valid and exact.
文摘The bi-conjugate gradients(Bi-CG)and bi-conjugate residual(Bi-CR)methods are powerful tools for solving nonsymmetric linear systems Ax=b.By using Kronecker product and vectorization operator,this paper develops the Bi-CG and Bi-CR methods for the solution of the generalized Sylvester-transpose matrix equationp i=1(Ai X Bi+Ci XTDi)=E(including Lyapunov,Sylvester and Sylvester-transpose matrix equations as special cases).Numerical results validate that the proposed algorithms are much more efcient than some existing algorithms.
基金This work was supported by the National Natural Science Foundation of China under grants 11371229, 11571026 and 11501326, and by the China Scholarship Council (File No. 2013083Y0102).
文摘We develop a fast stochastic Galerkin method for an optimal control problem governed by a random space-fractional diffusion equation with deterministic constrained control. Optimal control problems governed by a fractional diffusion equation tends to provide a better description for transport or conduction processes in heterogeneous media. Howev- er, the fractional control problem introduces significant computation complexity due to the nonlocal nature of fractional differential operators, and this is further worsen by the large number of random space dimensions to discretize the probability space. We ap- proximate the optimality system by a gradient algorithm combined with the stochastic Galerkin method through the discretization with respect to both the spatial space and the probability space. The resulting linear system can be decoupled for the random and spatial variable, and thus solved separately. A fast preconditioned Bi-Conjugate Gradient Stabilized method is developed to efficiently solve the decoupled systems derived from the fractional diffusion operators in the spatial space. Numerical experiments show the utility of the method.