|
|
Some practical aspects of parallel adaptive BDDC method
Šístek, Jakub ; Mandel, J. ; Sousedík, B.
We describe a parallel implementation of the Balancing Domain Decomposition by Constraints (BDDC) method enhanced by an adaptive construction of coarse problem. The method is designed for numerically difficult problems, where standard choice of continuity of arithmetic averages across faces and edges of subdomains fails to maintain the low condition number of the preconditioned system. Problems of elasticity analysis of bodies consisting of different materials with rapidly changing stiffness may represent one class of such challenging problems. The adaptive selection of constraints is shown to significantly increase the robustness of the method for this class of problems. However, since the cost of the set-up of the preconditioner with adaptive constraints is considerably larger than for the standard choices, computational feasibility of the presented implementation is obtained only for large contrasts of material coefficients.
|
|
|
Incompressible fluids with temperature dependent viscosity - numerical analysis and computational simulations
Ulrych, Oldřich ; Málek, Josef (advisor) ; Dolejší, Vít (referee) ; Šístek, Jakub (referee)
Title: Incompressible fluids with temperature dependent visco- sity, numerical analysis and computational simulations Author: RNDr. Oldřich Ulrych Department: Mathematical Institute of Charles University Supervisor: prof. RNDr. Josef Málek, CSc., DSc. Abstract: Flows of incompressible fluids connected with significant exchange of ther- mal and mechanical energy and with material moduli varying with the temperature and the shear rate, are described by the balance equations for linear momentum and energy, complemented by suitable constitution equations for the Cauchy stress and the heat flux. Assuming sufficient smoothness of quantities involved, the energy balance equation exhibits several equivalent formulations. However, within the context of weak solution, these formulations are, in general, not equivalent. This thesis is based on the existence theory for the generalized Navier-Stokes-Fourier system describing planar flow of fluids with a shear and temperature dependent vis- cosity. We specify parameters of a generalized power-law model under which weak formulations of balance equations are meaningful and both considered formulations of the energy balance equation are equivalent. Supported by the existence theory, we propose and numerically solve several problems pursuing the aim to systematically compare the...
|
|
|
Proceedings of the International Conference Applications of Mathematics 2015 : Prague, November 18-21, 2015
Brandts, J. ; Korotov, S. ; Křížek, Michal ; Segeth, Karel ; Šístek, Jakub ; Vejchodský, Tomáš
Professors Ivo Babuška, Milan Práger, and Emil Vitásek are renowned experts in numerical analysis and computational methods. Their fruitful scientific careers started in Prague, at the Institute of Mathematics of the Czechoslovak Academy of Sciences (now Czech Academy of Sciences). They collaborated there on various projects including the computational analysis of the construction technology for Orlík Dam. In 1966 they published their joint book entitled Numerical Processes in Differential Equations. It is an honor for the Institute of Mathematics to host a conference on the occasion of their birthdays.
|
|
|
Programs and Algorithms of Numerical Mathematics 17 : Dolní Maxov, June 8-13, 2014 : Proceedings of Seminar
Chleboun, J. ; Přikryl, Petr ; Segeth, Karel ; Šístek, Jakub ; Vejchodský, Tomáš
This volume comprises peer-reviewed papers that are based on invited lectures, survey lectures, short communications, and posters presented at the 17th seminar Programs and Algorithms of Numerical Mathematics (PANM) held in Dolní Maxov, Czech Republic, June 8–13, 2014. The seminar was organized by the Institute of Mathematics of the Academy of Sciences of the Czech Republic. It continued the previous seminars on mathematical software and numerical methods held (with only one exception) biannually in\nAlšovice, Bratříkov, Janov nad Nisou, Kořenov, Lázně Libverda, Dolní Maxov, and Prague in the period 1983–2012. The objective of this series of seminars is to provide a forum for presenting and discussing advanced topics in numerical analysis, singleor multi-processor applications of computational methods, and new approaches to mathematical modeling.
|
|
|
An application of the BDDC method to the Navier-Stokes equations in 3-D cavity
Hanek, M. ; Šístek, Jakub ; Burda, P.
We deal with numerical simulation of incompressible flow governed by the Navier-Stokes equations. The problem is discretised using the finite element method, and the arising system of nonlinear equations is solved by Picard iteration. We explore the applicability of the Balancing Domain Decomposition by Constraints (BDDC) method to nonsymmetric problems arising from such linearisation. One step of BDDC is applied as the preconditioner for the stabilized variant of the biconjugate gradient (BiCGstab) method. We present results for a 3-D cavity problem computed on 32 cores of a parallel supercomputer.
|
|
|
Different approaches to interface weights in the BDDC method in 3D
Čertíková, M. ; Šístek, Jakub ; Burda, P.
In this paper, we discuss the choice of weights in averaging of local (subdomain) solutions on the interface for the BDDC method (Balancing Domain Decomposition by Constraints). We try to find relations among different choices of the interface weights and compare them numerically on model problems of the Poisson equation and linear elasticity in 3D. Problems with jumps in coefficients of material properties are considered and both regular and irregular interfaces between subdomains are tested.
|
|
|
A parallel finite element solver for unsteady incompressible Navier-Stokes equations
Šístek, Jakub
A parallel solver for unsteady incompressible Navier-Stokes equations is presented. It is based on the finite element method combined with the pressure-correction approach. Semi-implicit treatment of the convective term is considered, leading to five systems of linear algebraic equations to be solved in each time step. Krylov subspace iterative methods are employed for the solution of these systems with a particular emphasis on efficient parallel preconditioners. A simulation of a benchmark problem of incompressible viscous flow around a sphere at Reynolds number 300 is presented and compared with literature.
|
|
| |
|
|
On selection of interface weights in domain decomposition methods
Čertíková, M. ; Šístek, Jakub ; Burda, P.
Different choices of the averaging operator within the BDDC method are compared on a series of 2D experiments. Subdomains with irregular interface and with jumps in material coefficients are included into the study. Two new approaches are studied along three standard choices. No approach is shown to be universally superior to others, and the resulting recommendation is that an actual method should be chosen based on properties of the problem.
|
|
| |
guest ::
