- document
-
Qiu, Y. (author), Van Gijzen, M.B. (author), Van Wingerden, J.W. (author), Verhaegen, M.H.G. (author), Vuik, C. (author)In this manuscript, we study preconditioning techniques for optimal in-domain control of the Navier-Stokes equation, where the control only acts on a few parts of the domain. Optimization and linearization of the optimal in-domain control problem results in a generalized linear saddle-point system. The Schur complement for the generalized saddle...report 2015
- document
-
Qiu, Y. (author), van Gijzen, M.B. (author), van Wingerden, J.W. (author), Verhaegen, M. (author), Vuik, C. (author)Multilevel sequentially semiseparable (MSSS) matrices form a class of structured matrices that have low-rank off-diagonal structure, which allows the matrix-matrix operations to be performed in linear computational complexity. MSSS preconditioners are computed by replacing the Schur complements in the block LU factorization of the global linear...report 2015
- document
-
Astudillo, R. (author), Van Gijzen, M.B. (author)This paper discusses the solution of large-scale linear matrix equations using the Induced Dimension reduction method (IDR(s)). IDR(s) was originally presented to solve system of linear equations, and is based on the IDR(s) theorem. We generalize the IDR(s) theorem to solve linear problems in any finite-dimensional space. This generalization...report 2015
- document
-
Qiu, Y. (author), Van Gijzen, M.B. (author), Van Wingerden, J.W. (author), Verhaegen, M. (author), Vuik, C. (author)In this manuscript, we study preconditioning techniques for optimal in-domain control of the Navier-Stokes equation, where the control only acts on a few parts of the domain. Optimization and linearization of the optimal in-domain control problem results in a generalized linear saddle-point system. The Schur complement for the generalized saddle...report 2015
- document
-
Zhao, J. (author), Vollebregt, E.A.H. (author), Oosterlee, C.W. (author)This paper presents a fast numerical solver for a nonlinear constrained optimization problem, arising from a 3D frictional contact problem. It incorporates an active set strategy with a nonlinear conjugate gradient method. One novelty is to consider the tractions of each slip element in a polar coordinate system, and use azimuth angles as...report 2014
- document
-
Qiu, Y. (author), Van Gijzen, M.B. (author), Van Wingerden, J.W. (author), Verhaegen, M. (author), Vuik, C. (author)In this paper, we consider preconditioning for PDE-constrained optimization problems. The underlying problems yield a linear saddle-point system. We study a class of preconditioners based on multilevel sequentially semiseparable (MSSS) matrix computations. The novel global preconditioner is to make use of the global structure of the saddle-point...report 2014
- document
-
He, X. (author), Vuik, C. (author)In this report we explore the performance of the SIMPLER , augmented Lagrangian, ’grad-div’ preconditioners and their variants for the two-by-two block systems arising in the incompressible Navier-Stokes equations. The lid-driven cavity and flow over a finite flat plate are chosen as the benchmark problems. For each problem Reynolds number...report 2013
- document
-
Qiu, Y. (author), Van Gijzen, M.B. (author), Van Wingerden, J. (author), Verhaegen, M. (author), Vuik, C. (author)This paper studies a new preconditioning technique for sparse systems arising from discretized partial differential equations (PDEs) in computational fluid dynamics (CFD), which exploit the multilevel sequentially semiseparable (MSSS) structure of the system matrix. MSSS matrix computations give a data-sparse way to approximate the LU...report 2013
- document
-
He, X. (author), Neytcheva, M. (author), Vuik, C. (author)This paper deals with fast and reliable numerical solution methods for the incompressible non-Newtonian Navier-Stokes equations. To handle the nonlinearity of the governing equations, the Picard and Newton methods are used to linearize these coupled partial differential equations. For space discretization we use the finite element method and...report 2013
- document
-
Jönsthövel, T.B. (author), Van Gijzen, M.B. (author), MacLachlan, S. (author), Vuik, C. (author), Scarpas, A. (author)The demand for large FE meshes increases as parallel computing becomes the standard in FE simulations. Direct and iterative solution methods are used to solve the resulting linear systems. Many applications concern composite materials, which are characterized by large discontinuities in the material properties. An example of such a material is...report 2011
- document
-
Jönsthövel, T.B. (author), Van Gijzen, M.B. (author), Vuik, C. (author), Scarpas, A. (author)Large discontinuities in material properties, such as encountered in composite materials, lead to ill-conditioned systems of linear equations. These discontinuities give rise to small eigenvalues that may negatively affect the convergence of iterative solution methods such as the Preconditioned Conjugate Gradient (PCG) method. This paper...report 2011
- document
-
van Gijzen, M.B. (author), Erlangga, Y.A. (author), Vuik, C. (author)Shifted Laplace preconditioners have attracted considerable attention as a technique to speed up convergence of iterative solution methods for the Helmholtz equation. In this paper we present a comprehensive spectral analysis of the Helmholtz operator preconditioned with a shifted Laplacian. Our analysis is valid under general conditions. The...report 2006
- document
-
Nabben, R. (author), Vuik, C. (author)In this paper we compare various preconditioners for the numerical solution of partial differential equations. We compare the well-known balancing Neumann Neumann preconditioner used in domain decomposition methods with a so-called deflation preconditioner. We prove that the effective condition number of the deflated preconditioned system is...report 2004
- document
-
Erlangga, Y.A. (author), Oosterlee, C.W. (author), Vuik, C. (author)An iterative solution method, in the form of a preconditioner for a Krylov subspace method, is presented for the Helmholtz equation. The preconditioner is based on a Helmholtz type differential operator with a complex term. A multigrid iteration is used for approximately inverting the preconditioner. The choice of multigrid components for the...report 2004
- document
-
Li, C. (author), Vuik, C. (author)In this paper, some theoretical results on the eigenvalue analysis of the SIMPLER preconditioning for incompressible now is presented. Some formulations have been derived to characterize the spectrum of the preconditioned matrix. These results could be helpful for the practical use of the SIMPLER preconditioning. Some numerical tests are reported.report 2003
- document
-
Nabben, R. (author), Vuik, C. (author)In this paper we compare various preconditioners for the numerical solution of partial dierential equations. We compare a coarse grid correction preconditioner used in domain decomposition methods with a so-called deflation preconditioner. We prove that the effective condition number of the de ated preconditioned system is always, i.e. for all...report 2003
- document
- Erlangga, Y.A. (author), Vuik, C. (author), Oosterlee, C.W. (author) report 2003
- document
-
Li, C. (author), Vuik, C. (author)In this paper, an eigenvalue analysis of the SIMPLE preconditioning for incompressible flow is presented. Some formulations have been set up to characterize the spectrum of the preconditioned matrix. This leads to a generalized eigenvalue problem. The generalized eigenvalue problem is investigated. Some eigenvalue bounds and the estimation for...report 2002
- document
-
Erlangga, Y.A. (author)In this report, several numerical aspects and diculties for solving a linear system derived from the time-harmonic wave equations are overviewed. The presentation begins with the derivation of the governing equation for waves propagating in general inhomogeneous media. Due to the need of numerical solutions, various discretizations based on...report 2002